Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/medical/dna/slowing/src/Run.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /examples/extended/medical/dna/slowing/src/Run.cc (Version 11.3.0) and /examples/extended/medical/dna/slowing/src/Run.cc (Version 9.4.p1)


  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 // This example is provided by the Geant4-DNA     
 27 // Any report or published results obtained us    
 28 // shall cite the following Geant4-DNA collabo    
 29 // Med. Phys. 45 (2018) e722-e739                 
 30 // Phys. Med. 31 (2015) 861-874                   
 31 // Med. Phys. 37 (2010) 4692-4708                 
 32 // Int. J. Model. Simul. Sci. Comput. 1 (2010)    
 33 //                                                
 34 // The Geant4-DNA web site is available at htt    
 35 //                                                
 36 /// \file medical/dna/slowing/src/Run.cc          
 37 /// \brief Implementation of the Run class        
 38                                                   
 39 #include "Run.hh"                                 
 40                                                   
 41 #include "HistoManager.hh"                        
 42 #include "PrimaryGeneratorAction.hh"              
 43                                                   
 44 #include "G4Material.hh"                          
 45 #include "G4SystemOfUnits.hh"                     
 46 #include "G4UnitsTable.hh"                        
 47                                                   
 48 //....oooOO0OOooo........oooOO0OOooo........oo    
 49                                                   
 50 Run::Run(const DetectorConstruction* detector)    
 51   : G4Run(), fDetector(detector), fParticle(0)    
 52 {}                                                
 53                                                   
 54 //....oooOO0OOooo........oooOO0OOooo........oo    
 55                                                   
 56 Run::~Run() {}                                    
 57                                                   
 58 //....oooOO0OOooo........oooOO0OOooo........oo    
 59                                                   
 60 void Run::SetPrimary(G4ParticleDefinition* par    
 61 {                                                 
 62   fParticle = particle;                           
 63   fEkin = energy;                                 
 64 }                                                 
 65                                                   
 66 //....oooOO0OOooo........oooOO0OOooo........oo    
 67                                                   
 68 void Run::AddEdep(G4double e)                     
 69 {                                                 
 70   fEdeposit += e;                                 
 71   fEdeposit2 += e * e;                            
 72 }                                                 
 73                                                   
 74 //....oooOO0OOooo........oooOO0OOooo........oo    
 75                                                   
 76 void Run::Merge(const G4Run* run)                 
 77 {                                                 
 78   const Run* localRun = static_cast<const Run*    
 79                                                   
 80   // Pass information about primary particle      
 81   fParticle = localRun->fParticle;                
 82   fEkin = localRun->fEkin;                        
 83                                                   
 84   // Accumulate sums                              
 85   fEdeposit += localRun->fEdeposit;               
 86   fEdeposit2 += localRun->fEdeposit2;             
 87                                                   
 88   G4Run::Merge(run);                              
 89 }                                                 
 90                                                   
 91 //....oooOO0OOooo........oooOO0OOooo........oo    
 92                                                   
 93 void Run::EndOfRun()                              
 94 {                                                 
 95   std::ios::fmtflags mode = G4cout.flags();       
 96   G4cout.setf(std::ios::fixed, std::ios::float    
 97   G4int prec = G4cout.precision(2);               
 98                                                   
 99   // Run conditions                               
100   G4String partName = fParticle->GetParticleNa    
101   G4Material* material = fDetector->GetAbsorMa    
102   G4double density = material->GetDensity();      
103                                                   
104   G4cout << "\n ======================== run s    
105   G4cout << "\n The run is " << numberOfEvent     
106          << G4BestUnit(fEkin, "Energy") << " t    
107          << " (density: " << G4BestUnit(densit    
108          << G4BestUnit(fDetector->GetAbsorMass    
109                                                   
110   if (numberOfEvent == 0) {                       
111     G4cout.setf(mode, std::ios::floatfield);      
112     G4cout.precision(prec);                       
113     return;                                       
114   }                                               
115                                                   
116   G4cout.precision(3);                            
117   G4cout << "\n Total Energy deposited            
118                                                   
119   /*                                              
120   G4double dose=fEdeposit/fDetector->GetAbsorM    
121   G4double rmsDose=rms/fDetector->GetAbsorMass    
122                                                   
123   G4cout.precision(3);                            
124   G4cout                                          
125     << "\n Dose                         = " <<    
126     << G4endl;                                    
127   */                                              
128                                                   
129   G4cout << G4endl;                               
130                                                   
131   // Normalize histograms                         
132   G4AnalysisManager* analysisManager = G4Analy    
133   G4double fac = 1. / (numberOfEvent * (fEkin     
134   analysisManager->ScaleH1(1, fac);               
135   analysisManager->ScaleH1(2, fac);               
136   analysisManager->ScaleH1(3, fac);               
137 }                                                 
138