Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // 27 28 #include "globals.hh" 29 #include "G4ios.hh" 30 #include "G4VCollision.hh" 31 #include "G4VCrossSectionSource.hh" 32 #include "G4KineticTrack.hh" 33 34 G4VCollision::G4VCollision() 35 { } 36 37 G4VCollision::~G4VCollision() 38 { 39 } 40 41 42 G4bool G4VCollision::operator==(const G4VColli 43 { 44 return (this == (G4VCollision *) &right); 45 } 46 47 48 G4bool G4VCollision::operator!=(const G4VColli 49 { 50 return (this != (G4VCollision *) &right); 51 } 52 53 54 G4double G4VCollision::CrossSection(const G4Ki 55 const G4Ki 56 { 57 G4double sigma = 0.; 58 59 const G4VCrossSectionSource* xSource = GetCr 60 61 if (xSource != nullptr) 62 { 63 // There is a cross section for this Colli 64 sigma = xSource->CrossSection(aTrk1,aTrk2) 65 } 66 return sigma; 67 } 68 69 70 void G4VCollision::Print() const 71 { 72 G4String name = GetName(); 73 74 G4cout << "---- " << name << "---- Cross sec 75 76 const G4VCrossSectionSource* xSource = GetCr 77 if (xSource) xSource->Print(); 78 79 G4int nComponents = 0; 80 const G4CollisionVector* components = GetCom 81 if (components) 82 { 83 nComponents = (G4int)components->size(); 84 } 85 G4cout << "---- " << name << "---- has " << 86 G4int i = 0; 87 if (components) 88 { 89 for (auto iter = components->cbegin(); ite 90 { 91 G4cout << "---- " << name << " ---- Comp 92 ((*iter))->Print(); 93 ++i; 94 } 95 } 96 } 97 98 99 void G4VCollision::Print(const G4KineticTrack& 100 const G4KineticTrack& 101 { 102 G4String name = GetName(); 103 104 if (IsInCharge(trk1,trk2)) 105 { 106 G4cout << "---- " << name << "is in charge 107 } 108 else 109 { 110 G4cout << "---- " << name << "is not in ch 111 } 112 113 G4cout << "---- " << name << "---- Cross sec 114 115 const G4VCrossSectionSource* xSource = GetCr 116 if (xSource) xSource->Print(); 117 G4cout << "Cross section = " << CrossSection 118 119 G4int nComponents = 0; 120 const G4CollisionVector* components = GetCom 121 if (components) 122 { 123 nComponents = (G4int)components->size(); 124 } 125 G4cout << "---- " << name << "has " << nComp 126 127 G4int i = 0; 128 if (components) 129 { 130 for (auto iter = components->cbegin(); ite 131 { 132 G4cout << "Component " << i << G4endl; 133 ((*iter))->Print(); 134 ++i; 135 } 136 } 137 } 138 139 G4VCollision::G4VCollision(void*, void*, void 140 { 141 } 142 143 void G4VCollision::establish_G4MT_TLS_G4VColli 144 { 145 } 146