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 //....oooOO0OOooo........oooOO0OOooo........oo 29 //....oooOO0OOooo........oooOO0OOooo........oo 30 31 #include "Run.hh" 32 33 #include "PrimaryGeneratorAction.hh" 34 35 //....oooOO0OOooo........oooOO0OOooo........oo 36 Run::Run(G4bool isP, G4int nbProjs, G4int nbS, 37 : G4Run(), 38 isPIXE(isP), 39 fProjectionIndex(0), 40 fSliceIndex(0), 41 fPixelIndex(0), 42 fNbProjections(nbProjs), 43 fNbSlices(nbS), 44 fNbPixels(nbP), 45 fResumeProjIndex(resumeProjIndex) 46 {} 47 48 //....oooOO0OOooo........oooOO0OOooo........oo 49 50 Run::~Run() = default; 51 52 //....oooOO0OOooo........oooOO0OOooo........oo 53 54 void Run::Merge(const G4Run* run) 55 { 56 G4cout << "--------------Merge is called---- 57 const Run* localRun = static_cast<const Run* 58 59 for (size_t i = 0; i != localRun->gammaAtExi 60 gammaAtExit.push_back(localRun->gammaAtExi 61 } 62 for (size_t i = 0; i != localRun->gammaAtCre 63 gammaAtCreation.push_back(localRun->gammaA 64 } 65 66 for (size_t i = 0; i != localRun->protonAtEx 67 protonAtExit.push_back(localRun->protonAtE 68 } 69 70 G4Run::Merge(run); 71 } 72 73 //....oooOO0OOooo........oooOO0OOooo........oo 74 75 void Run::EndOfRun() 76 { 77 // Write X-ray data or proton data into file 78 // G4int runID = GetRunID(); 79 fProjectionIndex = GetCurrentProjection(); 80 fSliceIndex = GetCurrentSlice(); 81 fPixelIndex = GetCurrentPixel(); 82 83 RunInfo runInfo((uint8_t)fProjectionIndex, ( 84 85 if (isPIXE) { 86 runInfo.nbParticle = (uint32_t)gammaAtCrea 87 WriteFile("GammaAtCreation.dat", runInfo, 88 89 runInfo.nbParticle = (uint32_t)gammaAtExit 90 WriteFile("GammaAtExit.dat", runInfo, gamm 91 } 92 else { 93 runInfo.nbParticle = (uint32_t)protonAtExi 94 WriteFile("ProtonAtExit.dat", runInfo, pro 95 } 96 97 ClearVecs(); 98 } 99 100 //....oooOO0OOooo........oooOO0OOooo........oo 101 102 G4bool Run::GetIsPIXE() 103 { 104 return isPIXE; 105 } 106 107 //....oooOO0OOooo........oooOO0OOooo........oo 108 109 void Run::FillGammaAtExit(ParticleInfo gammaIn 110 { 111 gammaAtExit.push_back(gammaInfo); 112 } 113 114 void Run::FillGammaAtCreation(ParticleInfo gam 115 { 116 gammaAtCreation.push_back(gammaInfo); 117 } 118 119 void Run::FillProtonAtExit(ParticleInfo proton 120 { 121 protonAtExit.push_back(protonInfo); 122 } 123 124 void Run::ClearVecs() 125 { 126 if (gammaAtExit.size()) { 127 gammaAtExit.clear(); 128 } 129 if (gammaAtCreation.size()) { 130 gammaAtCreation.clear(); 131 } 132 if (protonAtExit.size()) { 133 protonAtExit.clear(); 134 } 135 } 136 137 void Run::WriteFile(const std::string fName, R 138 { 139 std::ofstream ofs; 140 if (runID) { 141 ofs.open(fName.c_str(), 142 std::ios::out | std::ios::app | s 143 } 144 else { 145 // if runId ==0, delete the existing file 146 ofs.open(fName.c_str(), std::ios::out | st 147 } 148 ofs.write((const char*)&runInfo, sizeof(RunI 149 if (vec.size()) { 150 ofs.write((const char*)vec.data(), sizeof( 151 } 152 153 ofs.close(); 154 } 155 G4int Run::GetCurrentProjection() 156 { 157 if (fResumeProjIndex > fNbProjections - 1) { 158 G4Exception("fResumeProjIndex", "Run::GetC 159 "To resume a simulation, the s 160 "maximal index"); 161 } 162 G4int projIndex = fResumeProjIndex + runID / 163 return projIndex; 164 } 165 G4int Run::GetCurrentSlice() 166 { 167 G4int remain = runID % (fNbSlices * fNbPixel 168 return remain / fNbPixels; 169 } 170 171 G4int Run::GetCurrentPixel() 172 { 173 G4int remain = runID % (fNbSlices * fNbPixel 174 return remain % fNbPixels; 175 } 176 177 //....oooOO0OOooo........oooOO0OOooo........oo 178