| 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 //
28 // Author: Elena Guardincerri (Elena.Guardince
29 //
30 // History:
31 // -----------
32 // 28 Nov 2001 Elena Guardincerri Created
33 // 15 Jul 2003 Alfonso Mantero "Detecto
34 // 23 Sep 2003 Alfonso Mantero differnt
35 //
36 // -------------------------------------------
37
38 #include "XrayFluoEventAction.hh"
39 #include "XrayFluoSensorHit.hh"
40 #include "XrayFluoEventActionMessenger.hh"
41 #include "XrayFluoDataSet.hh"
42 #include "XrayFluoAnalysisManager.hh"
43
44 #include "G4SystemOfUnits.hh"
45 #include "G4Event.hh"
46 #include "G4EventManager.hh"
47 #include "G4HCofThisEvent.hh"
48 #include "G4VHitsCollection.hh"
49 #include "G4SDManager.hh"
50 #include "G4UImanager.hh"
51 #include "G4ios.hh"
52 #include "G4UnitsTable.hh"
53 #include "Randomize.hh"
54 #include <fstream>
55
56 //....oooOO0OOooo........oooOO0OOooo........oo
57
58 XrayFluoEventAction::XrayFluoEventAction(const
59 :drawFlag("all"),
60 HPGeCollID(0),
61 eventMessenger(0),
62 printModulo(1),
63 detectorType(0)
64 {
65 eventMessenger = new XrayFluoEventActionMess
66
67 if (!(det->GetPhaseSpaceFlag()) ){
68 detectorType = det->GetDetectorType();
69 HPGeCollID=-1;
70 }
71 }
72
73 //....oooOO0OOooo........oooOO0OOooo........oo
74
75 XrayFluoEventAction::XrayFluoEventAction(const
76 :drawFlag("all"),
77 HPGeCollID(-1),
78 eventMessenger(0),
79 printModulo(1),
80 detectorType(0)
81 {
82 eventMessenger = new XrayFluoEventActionMess
83 detectorType = det->GetDetectorType();
84 }
85
86 //....oooOO0OOooo........oooOO0OOooo........oo
87
88 XrayFluoEventAction::XrayFluoEventAction(const
89 :drawFlag("all"),
90 HPGeCollID(-1),
91 eventMessenger(0),
92 printModulo(1),
93 detectorType(0)
94 {
95 eventMessenger = new XrayFluoEventActionMess
96 detectorType = det->GetDetectorType();
97 }
98
99 //....oooOO0OOooo........oooOO0OOooo........oo
100
101 XrayFluoEventAction::~XrayFluoEventAction()
102 {
103 delete eventMessenger;
104 eventMessenger = 0;
105 }
106
107 //....oooOO0OOooo........oooOO0OOooo........oo
108
109 void XrayFluoEventAction::BeginOfEventAction(c
110 {
111
112 G4int eventNumber = (evt->GetEventID())+1;
113 if ( eventNumber == 1){
114
115 G4cout << "# = 100000 events" << G4endl;
116 G4cout << "1--------+---------+---------+---
117 }
118
119 if ( ((eventNumber) % 100000) == 0 ) {
120
121 if ( eventNumber % (G4int)5e6 != 0 ) G4cou
122 else G4cout << "#"<< G4endl;
123 }
124
125 if (HPGeCollID==-1)
126 {
127 G4SDManager * SDman = G4SDManager::GetSD
128 HPGeCollID = SDman->GetCollectionID("HPG
129 //the pointer points to the ID number of
130 }
131 }
132
133 //....oooOO0OOooo........oooOO0OOooo........oo
134
135 void XrayFluoEventAction::EndOfEventAction(con
136 {
137 if (detectorType) {
138
139 // extracted from hits, compute the total
140 // track length)
141 G4HCofThisEvent* HCE = evt->GetHCofThisEve
142
143 XrayFluoSensorHitsCollection* HPGeHC = 0;
144 G4int n_hit = 0;
145 G4double totEnergy=0., energyD=0.;
146
147 if (HCE) HPGeHC =
148 (XrayFluoSensorHitsCollection*)(HCE->Get
149
150 if(HPGeHC)
151 {
152 n_hit = HPGeHC->entries();
153
154 for (G4int i=0;i<n_hit;i++)
155 {
156 totEnergy += (*HPGeHC)[i]->GetEdepTot();
157 energyD = detectorType->ResponseFunction
158 XrayFluoAnalysisManager* analysis = Xray
159 analysis->analyseEnergyDep(energyD);
160 }
161 }
162 }
163 }
164
165 //....oooOO0OOooo........oooOO0OOooo........oo
166
167
168 G4double XrayFluoEventAction::RandomCut(G4doub
169
170 {
171 G4double efficiency = 1.;
172 G4double F = 0.15;
173 G4double epsilon = 2.96 * eV;
174 G4double deltaE = 220 * eV;
175 G4double EdepDetect = 0.;
176
177 G4double Random = G4UniformRand();
178
179 if ( Random<efficiency )
180 {
181 G4double sigma = std::sqrt(F*epsilon*ene
182 EdepDetect = G4RandGauss::shoot(energy,
183
184 }
185 else {EdepDetect = 0.;}
186
187 return EdepDetect;
188 }
189
190
191