ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
B1RunAction.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file B1RunAction.cc
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
29 
30 #include "B1RunAction.hh"
31 #include "B1PrimaryGeneratorAction.hh"
32 #include "B1DetectorConstruction.hh"
33 // #include "B1Run.hh"
34 
35 #include "G4RunManager.hh"
36 #include "G4Run.hh"
37 #include "G4AccumulableManager.hh"
38 #include "G4LogicalVolumeStore.hh"
39 #include "G4LogicalVolume.hh"
40 #include "G4UnitsTable.hh"
41 #include "G4SystemOfUnits.hh"
42 
43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
44 
46 : G4UserRunAction(),
47  fEdep(0.),
48  fEdep2(0.)
49 {
50  // add new units for dose
51  //
52  const G4double milligray = 1.e-3*gray;
53  const G4double microgray = 1.e-6*gray;
54  const G4double nanogray = 1.e-9*gray;
55  const G4double picogray = 1.e-12*gray;
56 
57  new G4UnitDefinition("milligray", "milliGy" , "Dose", milligray);
58  new G4UnitDefinition("microgray", "microGy" , "Dose", microgray);
59  new G4UnitDefinition("nanogray" , "nanoGy" , "Dose", nanogray);
60  new G4UnitDefinition("picogray" , "picoGy" , "Dose", picogray);
61 
62  // Register accumulable to the accumulable manager
64  accumulableManager->RegisterAccumulable(fEdep);
65  accumulableManager->RegisterAccumulable(fEdep2);
66 }
67 
68 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
69 
71 {}
72 
73 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
74 
76 {
77  // inform the runManager to save random number seed
79 
80  // reset accumulables to their initial values
82  accumulableManager->Reset();
83 
84 }
85 
86 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
87 
89 {
90  G4int nofEvents = run->GetNumberOfEvent();
91  if (nofEvents == 0) return;
92 
93  // Merge accumulables
95  accumulableManager->Merge();
96 
97  // Compute dose = total energy deposit in a run and its variance
98  //
100  G4double edep2 = fEdep2.GetValue();
101 
102  G4double rms = edep2 - edep*edep/nofEvents;
103  if (rms > 0.) rms = std::sqrt(rms); else rms = 0.;
104 
105  const B1DetectorConstruction* detectorConstruction
106  = static_cast<const B1DetectorConstruction*>
108  G4double mass = detectorConstruction->GetScoringVolume()->GetMass();
109  G4double dose = edep/mass;
110  G4double rmsDose = rms/mass;
111 
112  // Run conditions
113  // note: There is no primary generator action object for "master"
114  // run manager for multi-threaded mode.
115  const B1PrimaryGeneratorAction* generatorAction
116  = static_cast<const B1PrimaryGeneratorAction*>
118  G4String runCondition;
119  if (generatorAction)
120  {
121  const G4ParticleGun* particleGun = generatorAction->GetParticleGun();
122  runCondition += particleGun->GetParticleDefinition()->GetParticleName();
123  runCondition += " of ";
124  G4double particleEnergy = particleGun->GetParticleEnergy();
125  runCondition += G4BestUnit(particleEnergy,"Energy");
126  }
127 
128  // Print
129  //
130  if (IsMaster()) {
131  G4cout
132  << G4endl
133  << "--------------------End of Global Run-----------------------";
134  }
135  else {
136  G4cout
137  << G4endl
138  << "--------------------End of Local Run------------------------";
139  }
140 
141  G4cout
142  << G4endl
143  << " The run consists of " << nofEvents << " "<< runCondition
144  << G4endl
145  << " Cumulated dose per run, in scoring volume : "
146  << G4BestUnit(dose,"Dose") << " rms = " << G4BestUnit(rmsDose,"Dose")
147  << G4endl
148  << "------------------------------------------------------------"
149  << G4endl
150  << G4endl;
151 }
152 
153 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
154 
156 {
157  fEdep += edep;
158  fEdep2 += edep*edep;
159 }
160 
161 
162 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
163