ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
B1ConRunAction.cc
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file B1ConRunAction.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 "B1ConRunAction.hh"
31 #include "B1PrimaryGeneratorAction.hh"
32 #include "B1DetectorConstruction.hh"
33 #include "B1ConRun.hh"
34 
35 #include "G4RunManager.hh"
36 #include "G4LogicalVolumeStore.hh"
37 #include "G4LogicalVolume.hh"
38 #include "G4UnitsTable.hh"
39 #include "G4SystemOfUnits.hh"
40 
41 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
42 
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" , "Dose", milligray);
54  new G4UnitDefinition("microgray", "microGy" , "Dose", microgray);
55  new G4UnitDefinition("nanogray" , "nanoGy" , "Dose", nanogray);
56  new G4UnitDefinition("picogray" , "picoGy" , "Dose", picogray);
57 
58 }
59 
60 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
61 
63 {}
64 
65 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
66 
68 {
69  return new B1ConRun;
70 }
71 
72 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
73 
75 {
76  G4cout << "### Run " << aRun->GetRunID() << " start." << G4endl;
77 
78  //inform the runManager to save random number seed
80 
81  if (IsMaster()) {
82  fdose_tally = new G4ConvergenceTester("DOSE_TALLY");
83  //fdose_tally = new G4ConvergenceTester();
84  }
85 
86 }
87 
88 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
89 
91 {
92  G4int nofEvents = aRun->GetNumberOfEvent();
93  if (nofEvents == 0) return;
94 
95  const B1ConRun* b1ConRun = static_cast<const B1ConRun*>(aRun);
96 
97  // Compute dose
98  //
99  G4double edep = b1ConRun->GetEdep();
100  G4double edep2 = b1ConRun->GetEdep2();
101  G4double rms = edep2 - edep*edep/nofEvents;
102  if (rms > 0.) rms = std::sqrt(rms); else rms = 0.;
103 
104  const B1DetectorConstruction* detectorConstruction
105  = static_cast<const B1DetectorConstruction*>
107  G4double mass = detectorConstruction->GetScoringVolume()->GetMass();
108  G4double dose = edep/mass;
109  G4double rmsDose = rms/mass;
110 
111  // Run conditions
112  // note: There is no primary generator action object for "master"
113  // run manager for multi-threaded mode.
114  const B1PrimaryGeneratorAction* generatorAction
115  = static_cast<const B1PrimaryGeneratorAction*>
117  G4String runCondition;
118  if (generatorAction)
119  {
120  const G4ParticleGun* particleGun = generatorAction->GetParticleGun();
121  runCondition += particleGun->GetParticleDefinition()->GetParticleName();
122  runCondition += " of ";
123  G4double particleEnergy = particleGun->GetParticleEnergy();
124  runCondition += G4BestUnit(particleEnergy,"Energy");
125  }
126 
127  // Print
128  //
129  if (IsMaster())
130  {
131 
132  for ( G4int i = 0 ; i != b1ConRun->GetNumberOfEvent(); i++ ) {
133  G4double aDose = b1ConRun->GetEdepPerEvent(i)/mass/gray;
134  fdose_tally->AddScore( aDose );
135  }
136 
139  delete fdose_tally;
140 
141  G4cout
142  << "\n--------------------End of Global Run-----------------------";
143  }
144  else
145  {
146  G4cout
147  << "\n--------------------End of Local Run------------------------";
148  }
149  G4cout
150  << "\n The run consists of " << nofEvents << " "<< runCondition
151  << "\n Dose in scoring volume : "
152  << G4BestUnit(dose,"Dose") << " +- " << G4BestUnit(rmsDose,"Dose")
153  << "\n------------------------------------------------------------\n"
154  << G4endl;
155 
156 }
157 
158 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......