Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer 3 // * License and Disclaimer * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. 9 // * include a list of copyright holders. * 10 // * 10 // * * 11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. * 17 // * 17 // * * 18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************* 24 // ******************************************************************** 25 // 25 // 26 // 26 // 27 // 27 // 28 28 29 // ------------------------------------------- 29 // ------------------------------------------------------------ 30 // GEANT 4 class implementation file 30 // GEANT 4 class implementation file 31 // 31 // 32 // ---------------- G4ExcitedString ----- 32 // ---------------- G4ExcitedString ---------------- 33 // by Gunter Folger, June 1998. 33 // by Gunter Folger, June 1998. 34 // class for an excited string used by P 34 // class for an excited string used by Parton String Models 35 // ------------------------------------------- 35 // ------------------------------------------------------------ 36 36 37 37 38 // G4ExcitedString 38 // G4ExcitedString 39 #include "G4ExcitedString.hh" 39 #include "G4ExcitedString.hh" 40 #include <algorithm> 40 #include <algorithm> 41 41 42 42 43 G4ExcitedString::G4ExcitedString(G4Parton* Col 43 G4ExcitedString::G4ExcitedString(G4Parton* Color, G4Parton* AntiColor, G4int Direction) 44 { 44 { 45 thePartons.push_back(Color); 45 thePartons.push_back(Color); 46 thePartons.push_back(AntiColor); 46 thePartons.push_back(AntiColor); 47 theTimeOfCreation = 0.; 47 theTimeOfCreation = 0.; // Uzhi 15.05.08 48 thePosition = Color->GetPosition(); 48 thePosition = Color->GetPosition(); 49 theDirection = Direction; 49 theDirection = Direction; 50 theTrack=0; 50 theTrack=0; 51 } 51 } 52 52 53 G4ExcitedString::G4ExcitedString(G4Parton* Col 53 G4ExcitedString::G4ExcitedString(G4Parton* Color, G4Parton* Gluon, G4Parton* AntiColor, G4int Direction) 54 { 54 { 55 thePartons.push_back(Color); 55 thePartons.push_back(Color); 56 thePartons.push_back(Gluon); 56 thePartons.push_back(Gluon); 57 thePartons.push_back(AntiColor); 57 thePartons.push_back(AntiColor); 58 theTimeOfCreation = 0.; 58 theTimeOfCreation = 0.; // Uzhi 15.05.08 59 thePosition = Color->GetPosition(); 59 thePosition = Color->GetPosition(); 60 theDirection = Direction; 60 theDirection = Direction; 61 theTrack=0; 61 theTrack=0; 62 } 62 } 63 63 64 G4ExcitedString::G4ExcitedString(G4KineticTrac 64 G4ExcitedString::G4ExcitedString(G4KineticTrack * track) 65 { 65 { 66 theTimeOfCreation = track->GetFormatio 66 theTimeOfCreation = track->GetFormationTime(); // Uzhi 15.05.08 67 thePosition = track->GetPosition(); 67 thePosition = track->GetPosition(); 68 theTrack= track; 68 theTrack= track; 69 theDirection=0; 69 theDirection=0; 70 } 70 } 71 71 72 G4ExcitedString::~G4ExcitedString() 72 G4ExcitedString::~G4ExcitedString() 73 { 73 { 74 std::for_each(thePartons.begin(), thePartons 74 std::for_each(thePartons.begin(), thePartons.end(), DeleteParton()); 75 if ( theTrack ) { 75 if ( theTrack ) { 76 delete theTrack; 76 delete theTrack; 77 theTrack=0; 77 theTrack=0; 78 } 78 } 79 } 79 } 80 80 81 81 82 //const G4ExcitedString & G4ExcitedString::ope 82 //const G4ExcitedString & G4ExcitedString::operator=(const G4ExcitedString &right) 83 //{} 83 //{} 84 84 85 85 86 //G4bool G4ExcitedString::operator==(const G4E 86 //G4bool G4ExcitedString::operator==(const G4ExcitedString &right) const 87 //{} 87 //{} 88 88 89 //G4bool G4ExcitedString::operator!=(const G4E 89 //G4bool G4ExcitedString::operator!=(const G4ExcitedString &right) const 90 //{} 90 //{} 91 91 92 92 93 93 94 // Additional Declarations 94 // Additional Declarations 95 95 96 96 97 void G4ExcitedString::Boost(G4ThreeVector& Vel 97 void G4ExcitedString::Boost(G4ThreeVector& Velocity) 98 { 98 { 99 for(unsigned int cParton = 0; cParton < th 99 for(unsigned int cParton = 0; cParton < thePartons.size() ; cParton++ ) 100 { 100 { 101 G4LorentzVector Mom = thePartons[cPart 101 G4LorentzVector Mom = thePartons[cParton]->Get4Momentum(); 102 Mom.boost(Velocity); 102 Mom.boost(Velocity); 103 thePartons[cParton]->Set4Momentum(Mom) 103 thePartons[cParton]->Set4Momentum(Mom); 104 } 104 } 105 } 105 } 106 106 107 //-------------------------------------------- 107 //--------------------------------------------------------------------------------- 108 108 109 G4Parton* G4ExcitedString::GetColorParton(void 109 G4Parton* G4ExcitedString::GetColorParton(void) const 110 { 110 { 111 G4Parton * start = *(thePartons.begin()); 111 G4Parton * start = *(thePartons.begin()); 112 G4Parton * end = *(thePartons.end()-1); 112 G4Parton * end = *(thePartons.end()-1); 113 G4int Encoding = start->GetPDGcode(); 113 G4int Encoding = start->GetPDGcode(); 114 if (Encoding < -1000 || ((Encoding < 1000 114 if (Encoding < -1000 || ((Encoding < 1000) && (Encoding > 0))) 115 return start; 115 return start; 116 return end; 116 return end; 117 } 117 } 118 118 119 //-------------------------------------------- 119 //--------------------------------------------------------------------------------- 120 120 121 G4Parton* G4ExcitedString::GetGluon(void) cons 121 G4Parton* G4ExcitedString::GetGluon(void) const 122 { 122 { 123 return thePartons[1]; 123 return thePartons[1]; 124 } 124 } 125 125 126 //-------------------------------------------- 126 //--------------------------------------------------------------------------------- 127 127 128 G4Parton* G4ExcitedString::GetGluon(G4int Gluo 128 G4Parton* G4ExcitedString::GetGluon(G4int GluonPos) const 129 { 129 { 130 return thePartons[1 + GluonPos]; 130 return thePartons[1 + GluonPos]; 131 } 131 } 132 132 133 //-------------------------------------------- 133 //--------------------------------------------------------------------------------- 134 134 135 G4Parton* G4ExcitedString::GetAntiColorParton( 135 G4Parton* G4ExcitedString::GetAntiColorParton(void) const 136 { 136 { 137 G4Parton * start = *(thePartons.begin()); 137 G4Parton * start = *(thePartons.begin()); 138 G4Parton * end = *(thePartons.end()-1); 138 G4Parton * end = *(thePartons.end()-1); 139 G4int Encoding = start->GetPDGcode(); 139 G4int Encoding = start->GetPDGcode(); 140 if (Encoding < -1000 || ((Encoding < 1000 140 if (Encoding < -1000 || ((Encoding < 1000) && (Encoding > 0))) 141 return end; 141 return end; 142 return start; 142 return start; 143 } 143 } 144 144 145 //-------------------------------------------- 145 //--------------------------------------------------------------------------------- 146 146 147 G4bool G4ExcitedString::IsItKinkyString(void) 147 G4bool G4ExcitedString::IsItKinkyString(void) const 148 { 148 { 149 return (thePartons.size() > 2); 149 return (thePartons.size() > 2); 150 } 150 } 151 151 152 //-------------------------------------------- 152 //--------------------------------------------------------------------------------- 153 153 154 G4int G4ExcitedString::GetDirection(void) cons 154 G4int G4ExcitedString::GetDirection(void) const 155 { 155 { 156 return theDirection; 156 return theDirection; 157 } 157 } 158 158 159 //******************************************** 159 //********************************************************************************* 160 160 161 G4Parton* G4ExcitedString::GetLeftParton(void) 161 G4Parton* G4ExcitedString::GetLeftParton(void) const 162 { 162 { 163 return *thePartons.begin(); 163 return *thePartons.begin(); 164 } 164 } 165 165 166 //-------------------------------------------- 166 //--------------------------------------------------------------------------------- 167 167 168 G4Parton* G4ExcitedString::GetRightParton(void 168 G4Parton* G4ExcitedString::GetRightParton(void) const 169 { 169 { 170 return *(thePartons.end()-1); 170 return *(thePartons.end()-1); 171 } 171 } 172 172 173 //******************************************** 173 //********************************************************************************* 174 174