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 // G4PrimaryVertex 27 // 28 // Authors: G.Cosmo, 2 December 1995 - Design, 29 // M.Asai, 29 January 1996 - First im 30 // ------------------------------------------- 31 32 #include "G4PrimaryVertex.hh" 33 34 #include "G4SystemOfUnits.hh" 35 #include "G4VUserPrimaryVertexInformation.hh" 36 #include "G4ios.hh" 37 38 G4Allocator<G4PrimaryVertex>*& aPrimaryVertexA 39 { 40 G4ThreadLocalStatic G4Allocator<G4PrimaryVer 41 return _instance; 42 } 43 44 G4PrimaryVertex::G4PrimaryVertex(G4double x0, 45 : X0(x0), Y0(y0), Z0(z0), T0(t0) 46 {} 47 48 G4PrimaryVertex::G4PrimaryVertex(G4ThreeVector 49 : X0(xyz0.x()), Y0(xyz0.y()), Z0(xyz0.z()), 50 {} 51 52 G4PrimaryVertex::G4PrimaryVertex(const G4Prima 53 { 54 *this = right; 55 } 56 57 G4PrimaryVertex::~G4PrimaryVertex() 58 { 59 if (theParticle != nullptr) { 60 G4PrimaryParticle* theNext = theParticle; 61 while (theNext != nullptr) { 62 G4PrimaryParticle* thisPrimary = theNext 63 theNext = thisPrimary->GetNext(); 64 thisPrimary->ClearNext(); 65 delete thisPrimary; 66 } 67 theParticle = nullptr; 68 } 69 delete nextVertex; 70 nextVertex = nullptr; 71 72 theTail = nullptr; 73 tailVertex = nullptr; 74 75 delete userInfo; 76 userInfo = nullptr; 77 } 78 79 G4PrimaryVertex& G4PrimaryVertex::operator=(co 80 { 81 if (this != &right) { 82 X0 = right.X0; 83 Y0 = right.Y0; 84 Z0 = right.Z0; 85 T0 = right.T0; 86 Weight0 = right.Weight0; 87 numberOfParticle = right.numberOfParticle; 88 89 delete theParticle; 90 theParticle = nullptr; 91 theTail = nullptr; 92 if (right.theParticle != nullptr) { 93 theParticle = new G4PrimaryParticle(*(ri 94 theTail = theParticle; 95 G4PrimaryParticle* np = theParticle->Get 96 while (np != nullptr) // Loop checking, 97 { 98 theTail = np; 99 np = np->GetNext(); 100 } 101 } 102 103 delete nextVertex; 104 nextVertex = nullptr; 105 tailVertex = nullptr; 106 if (right.nextVertex != nullptr) { 107 nextVertex = new G4PrimaryVertex(*(right 108 tailVertex = nextVertex; 109 G4PrimaryVertex* nv = nextVertex->GetNex 110 while (nv != nullptr) // Loop checking, 111 { 112 tailVertex = nv; 113 nv = nv->GetNext(); 114 } 115 } 116 117 // userInfo cannot be copied 118 userInfo = nullptr; 119 } 120 return *this; 121 } 122 123 G4bool G4PrimaryVertex::operator==(const G4Pri 124 { 125 return (this == &right); 126 } 127 128 G4bool G4PrimaryVertex::operator!=(const G4Pri 129 { 130 return (this != &right); 131 } 132 133 G4PrimaryParticle* G4PrimaryVertex::GetPrimary 134 { 135 if (i >= 0 && i < numberOfParticle) { 136 G4PrimaryParticle* particle = theParticle; 137 for (G4int j = 0; j < i; ++j) { 138 if (particle == nullptr) return nullptr; 139 particle = particle->GetNext(); 140 } 141 return particle; 142 } 143 144 return nullptr; 145 } 146 147 void G4PrimaryVertex::Print() const 148 { 149 G4cout << "Vertex ( " << X0 / mm << "[mm], 150 << T0 / ns << "[ns] )" 151 << " Weight " << Weight0 << G4endl; 152 if (userInfo != nullptr) userInfo->Print(); 153 G4cout << " -- Primary particles :: " 154 << " # of primaries =" << numberOfP 155 if (theParticle != nullptr) theParticle->Pri 156 if (nextVertex != nullptr) { 157 G4cout << "Next Vertex " << G4endl; 158 nextVertex->Print(); 159 } 160 } 161