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 /// \file B1RunAction.cc 28 /// \brief Implementation of the B1RunAction c 29 30 #include "B1RunAction.hh" 31 32 #include "B1DetectorConstruction.hh" 33 #include "B1PrimaryGeneratorAction.hh" 34 #include "B1Run.hh" 35 36 #include "G4LogicalVolume.hh" 37 #include "G4LogicalVolumeStore.hh" 38 #include "G4RunManager.hh" 39 #include "G4SystemOfUnits.hh" 40 #include "G4UnitsTable.hh" 41 42 //....oooOO0OOooo........oooOO0OOooo........oo 43 44 B1RunAction::B1RunAction() : G4UserRunAction() 45 { 46 // add new units for dose 47 // 48 const G4double milligray = 1.e-3 * gray; 49 const G4double microgray = 1.e-6 * gray; 50 const G4double nanogray = 1.e-9 * gray; 51 const G4double picogray = 1.e-12 * gray; 52 53 new G4UnitDefinition("milligray", "milliGy", 54 new G4UnitDefinition("microgray", "microGy", 55 new G4UnitDefinition("nanogray", "nanoGy", " 56 new G4UnitDefinition("picogray", "picoGy", " 57 } 58 59 //....oooOO0OOooo........oooOO0OOooo........oo 60 61 B1RunAction::~B1RunAction() {} 62 63 //....oooOO0OOooo........oooOO0OOooo........oo 64 65 G4Run* B1RunAction::GenerateRun() 66 { 67 return new B1Run; 68 } 69 70 //....oooOO0OOooo........oooOO0OOooo........oo 71 72 void B1RunAction::BeginOfRunAction(const G4Run 73 { 74 // inform the runManager to save random numb 75 G4RunManager::GetRunManager()->SetRandomNumb 76 } 77 78 //....oooOO0OOooo........oooOO0OOooo........oo 79 80 void B1RunAction::EndOfRunAction(const G4Run* 81 { 82 G4int nofEvents = run->GetNumberOfEvent(); 83 if (nofEvents == 0) return; 84 85 const B1Run* b1Run = static_cast<const B1Run 86 87 // Compute dose 88 // 89 G4double edep = b1Run->GetEdep(); 90 G4double edep2 = b1Run->GetEdep2(); 91 G4double rms = edep2 - edep * edep / nofEven 92 if (rms > 0.) 93 rms = std::sqrt(rms); 94 else 95 rms = 0.; 96 97 const B1DetectorConstruction* detectorConstr 98 G4RunManager::GetRunManager()->GetUserDete 99 G4double mass = detectorConstruction->GetSco 100 G4double dose = edep / mass; 101 G4double rmsDose = rms / mass; 102 103 // Run conditions 104 // note: There is no primary generator acti 105 // run manager for multi-threaded mod 106 const B1PrimaryGeneratorAction* generatorAct 107 G4RunManager::GetRunManager()->GetUserPrim 108 G4String runCondition; 109 if (generatorAction) { 110 const G4ParticleGun* particleGun = generat 111 runCondition += particleGun->GetParticleDe 112 runCondition += " of "; 113 G4double particleEnergy = particleGun->Get 114 runCondition += G4BestUnit(particleEnergy, 115 } 116 117 // Print 118 // 119 if (IsMaster()) { 120 G4cout << "\n--------------------End of Gl 121 } 122 else { 123 G4cout << "\n--------------------End of Lo 124 } 125 126 G4cout << "\n The run consists of " << nofEv 127 << "\n Dose in scoring volume : " << 128 << G4BestUnit(rmsDose, "Dose") 129 << "\n------------------------------- 130 << G4endl; 131 } 132 133 //....oooOO0OOooo........oooOO0OOooo........oo 134