Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/xray_fluorescence/src/XrayFluoEventAction.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/advanced/xray_fluorescence/src/XrayFluoEventAction.cc (Version 11.3.0) and /examples/advanced/xray_fluorescence/src/XrayFluoEventAction.cc (Version 10.6.p1)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 //                                                 26 //
 27 //                                                 27 //
 28 // Author: Elena Guardincerri (Elena.Guardince     28 // Author: Elena Guardincerri (Elena.Guardincerri@ge.infn.it)
 29 //                                                 29 //
 30 // History:                                        30 // History:
 31 // -----------                                     31 // -----------
 32 // 28 Nov 2001 Elena Guardincerri     Created      32 // 28 Nov 2001 Elena Guardincerri     Created
 33 // 15 Jul 2003 Alfonso Mantero        "Detecto     33 // 15 Jul 2003 Alfonso Mantero        "DetectorType" use integration
 34 // 23 Sep 2003 Alfonso Mantero        differnt     34 // 23 Sep 2003 Alfonso Mantero        differnt geometries integration
 35 //                                                 35 //
 36 // -------------------------------------------     36 // -------------------------------------------------------------------
 37                                                    37 
 38 #include "XrayFluoEventAction.hh"                  38 #include "XrayFluoEventAction.hh"
 39 #include "XrayFluoSensorHit.hh"                    39 #include "XrayFluoSensorHit.hh"
 40 #include "XrayFluoEventActionMessenger.hh"         40 #include "XrayFluoEventActionMessenger.hh"
 41 #include "XrayFluoDataSet.hh"                      41 #include "XrayFluoDataSet.hh"
 42 #include "XrayFluoAnalysisManager.hh"              42 #include "XrayFluoAnalysisManager.hh"
 43                                                    43 
 44 #include "G4SystemOfUnits.hh"                      44 #include "G4SystemOfUnits.hh"
 45 #include "G4Event.hh"                              45 #include "G4Event.hh"
 46 #include "G4EventManager.hh"                       46 #include "G4EventManager.hh"
 47 #include "G4HCofThisEvent.hh"                      47 #include "G4HCofThisEvent.hh"
 48 #include "G4VHitsCollection.hh"                    48 #include "G4VHitsCollection.hh"
 49 #include "G4SDManager.hh"                          49 #include "G4SDManager.hh"
 50 #include "G4UImanager.hh"                          50 #include "G4UImanager.hh"
 51 #include "G4ios.hh"                                51 #include "G4ios.hh"
 52 #include "G4UnitsTable.hh"                         52 #include "G4UnitsTable.hh"
 53 #include "Randomize.hh"                            53 #include "Randomize.hh"
 54 #include <fstream>                                 54 #include <fstream>
 55                                                    55 
 56 //....oooOO0OOooo........oooOO0OOooo........oo     56 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 57                                                    57 
 58 XrayFluoEventAction::XrayFluoEventAction(const     58 XrayFluoEventAction::XrayFluoEventAction(const XrayFluoDetectorConstruction* det)
 59   :drawFlag("all"),                                59   :drawFlag("all"),
 60    HPGeCollID(0),                                  60    HPGeCollID(0),
 61    eventMessenger(0),                              61    eventMessenger(0),
 62    printModulo(1),                                 62    printModulo(1),
 63    detectorType(0)                                 63    detectorType(0)
 64 {                                                  64 {
 65   eventMessenger = new XrayFluoEventActionMess     65   eventMessenger = new XrayFluoEventActionMessenger(this);
 66                                                    66   
 67   if (!(det->GetPhaseSpaceFlag()) ){               67   if (!(det->GetPhaseSpaceFlag()) ){
 68     detectorType = det->GetDetectorType();         68     detectorType = det->GetDetectorType();
 69     HPGeCollID=-1;                                 69     HPGeCollID=-1;
 70   }                                                70   }
 71 }                                                  71 }
 72                                                    72 
 73 //....oooOO0OOooo........oooOO0OOooo........oo     73 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 74                                                    74 
 75 XrayFluoEventAction::XrayFluoEventAction(const     75 XrayFluoEventAction::XrayFluoEventAction(const XrayFluoPlaneDetectorConstruction* det)
 76   :drawFlag("all"),                                76   :drawFlag("all"),
 77    HPGeCollID(-1),                                 77    HPGeCollID(-1),
 78    eventMessenger(0),                              78    eventMessenger(0),
 79    printModulo(1),                                 79    printModulo(1),
 80    detectorType(0)                                 80    detectorType(0)
 81 {                                                  81 {
 82   eventMessenger = new XrayFluoEventActionMess     82   eventMessenger = new XrayFluoEventActionMessenger(this);
 83   detectorType = det->GetDetectorType();           83   detectorType = det->GetDetectorType();
 84 }                                                  84 }
 85                                                    85 
 86 //....oooOO0OOooo........oooOO0OOooo........oo     86 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 87                                                    87 
 88 XrayFluoEventAction::XrayFluoEventAction(const     88 XrayFluoEventAction::XrayFluoEventAction(const XrayFluoMercuryDetectorConstruction* det)
 89   :drawFlag("all"),                                89   :drawFlag("all"),
 90    HPGeCollID(-1),                                 90    HPGeCollID(-1),
 91    eventMessenger(0),                              91    eventMessenger(0),
 92    printModulo(1),                                 92    printModulo(1),
 93    detectorType(0)                                 93    detectorType(0)
 94 {                                                  94 {
 95   eventMessenger = new XrayFluoEventActionMess     95   eventMessenger = new XrayFluoEventActionMessenger(this);
 96   detectorType = det->GetDetectorType();           96   detectorType = det->GetDetectorType();
 97 }                                                  97 }
 98                                                    98 
 99 //....oooOO0OOooo........oooOO0OOooo........oo     99 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
100                                                   100 
101 XrayFluoEventAction::~XrayFluoEventAction()       101 XrayFluoEventAction::~XrayFluoEventAction()
102 {                                                 102 {
103    delete eventMessenger;                         103    delete eventMessenger;
104    eventMessenger = 0;                            104    eventMessenger = 0;
105 }                                                 105 }
106                                                   106 
107 //....oooOO0OOooo........oooOO0OOooo........oo    107 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
108                                                   108 
109 void XrayFluoEventAction::BeginOfEventAction(c    109 void XrayFluoEventAction::BeginOfEventAction(const G4Event* evt)
110 {                                                 110 {
111                                                   111 
112   G4int eventNumber = (evt->GetEventID())+1;      112   G4int eventNumber = (evt->GetEventID())+1;
113   if ( eventNumber == 1){                         113   if ( eventNumber == 1){
114                                                   114 
115   G4cout << "# = 100000 events" << G4endl;        115   G4cout << "# = 100000 events" << G4endl;
116   G4cout << "1--------+---------+---------+---    116   G4cout << "1--------+---------+---------+---------+---------5e6"<<G4endl;
117   }                                               117   }
118                                                   118 
119   if ( ((eventNumber) % 100000) == 0 )  {         119   if ( ((eventNumber) % 100000) == 0 )  {
120                                                   120 
121     if ( eventNumber % (G4int)5e6 != 0 ) G4cou    121     if ( eventNumber % (G4int)5e6 != 0 ) G4cout << "#" << std::flush;
122     else G4cout << "#"<< G4endl;                  122     else G4cout << "#"<< G4endl;   
123   }                                               123   }
124                                                   124 
125   if (HPGeCollID==-1)                             125   if (HPGeCollID==-1)    
126     {                                             126     {
127       G4SDManager * SDman = G4SDManager::GetSD    127       G4SDManager * SDman = G4SDManager::GetSDMpointer();
128       HPGeCollID = SDman->GetCollectionID("HPG    128       HPGeCollID = SDman->GetCollectionID("HPGeCollection");
129       //the pointer points to the ID number of    129       //the pointer points to the ID number of the sensitive detector
130     }                                             130     }
131 }                                                 131 }
132                                                   132   
133 //....oooOO0OOooo........oooOO0OOooo........oo    133 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
134                                                   134 
135 void XrayFluoEventAction::EndOfEventAction(con    135 void XrayFluoEventAction::EndOfEventAction(const G4Event* evt)
136 {                                                 136 {
                                                   >> 137   
137   if (detectorType) {                             138   if (detectorType) {
138                                                   139     
139     // extracted from hits, compute the total     140     // extracted from hits, compute the total energy deposit (and total charged
140     // track length)                              141     // track length) 
141     G4HCofThisEvent* HCE = evt->GetHCofThisEve    142     G4HCofThisEvent* HCE = evt->GetHCofThisEvent();
142                                                   143     
143     XrayFluoSensorHitsCollection* HPGeHC = 0;     144     XrayFluoSensorHitsCollection* HPGeHC = 0;
144     G4int n_hit = 0;                              145     G4int n_hit = 0;
145     G4double totEnergy=0., energyD=0.;         << 146     G4double totEnergyDetect=0., totEnergy=0., energyD=0.;
146                                                   147     
147     if (HCE) HPGeHC =                             148     if (HCE) HPGeHC = 
148       (XrayFluoSensorHitsCollection*)(HCE->Get    149       (XrayFluoSensorHitsCollection*)(HCE->GetHC(HPGeCollID));
149                                                   150 
150     if(HPGeHC)                                    151     if(HPGeHC)            
151       {                                           152       {
152   n_hit = HPGeHC->entries();                      153   n_hit = HPGeHC->entries();
153                                                   154     
154   for (G4int i=0;i<n_hit;i++)                     155   for (G4int i=0;i<n_hit;i++)
155     {                                             156     {
156       totEnergy += (*HPGeHC)[i]->GetEdepTot();    157       totEnergy += (*HPGeHC)[i]->GetEdepTot(); 
157       energyD = detectorType->ResponseFunction    158       energyD = detectorType->ResponseFunction(totEnergy);      
158       XrayFluoAnalysisManager* analysis = Xray    159       XrayFluoAnalysisManager* analysis = XrayFluoAnalysisManager::getInstance();
159       analysis->analyseEnergyDep(energyD);        160       analysis->analyseEnergyDep(energyD);
                                                   >> 161       totEnergyDetect += energyD;           
160     }                                             162     }
161       }                                           163       }
162   }                                               164   }  
163 }                                                 165 }
164                                                   166 
165 //....oooOO0OOooo........oooOO0OOooo........oo    167 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
166                                                   168 
167                                                   169 
168 G4double XrayFluoEventAction::RandomCut(G4doub    170 G4double XrayFluoEventAction::RandomCut(G4double energy)
169                                                   171   
170 {                                                 172 {
171   G4double efficiency = 1.;                       173   G4double efficiency = 1.;
172   G4double F = 0.15;                              174   G4double F = 0.15;
173   G4double epsilon = 2.96 * eV;                   175   G4double epsilon = 2.96 * eV;
174   G4double deltaE = 220 * eV;                     176   G4double deltaE = 220 * eV;
175   G4double EdepDetect = 0.;                       177   G4double EdepDetect = 0.;
176                                                   178 
177   G4double  Random = G4UniformRand();             179   G4double  Random = G4UniformRand(); 
178                                                   180 
179   if ( Random<efficiency )                        181   if ( Random<efficiency )
180     {                                             182     {
181       G4double sigma = std::sqrt(F*epsilon*ene    183       G4double sigma = std::sqrt(F*epsilon*energy+std::pow(deltaE/2355,2));      
182       EdepDetect = G4RandGauss::shoot(energy,     184       EdepDetect = G4RandGauss::shoot(energy, sigma );
183                                                   185       
184     }                                             186     }
185   else {EdepDetect = 0.;}                         187   else {EdepDetect = 0.;}
186                                                   188 
187   return   EdepDetect;                            189   return   EdepDetect;  
188 }                                                 190 }
189                                                   191 
190                                                   192 
191                                                   193