| Geant4 Cross Reference |
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 // $Id: XrayFluoEventAction.cc
>> 28 // GEANT4 tag $Name: xray_fluo-V03-02-00
27 // 29 //
28 // Author: Elena Guardincerri (Elena.Guardince 30 // Author: Elena Guardincerri (Elena.Guardincerri@ge.infn.it)
29 // 31 //
30 // History: 32 // History:
31 // ----------- 33 // -----------
32 // 28 Nov 2001 Elena Guardincerri Created 34 // 28 Nov 2001 Elena Guardincerri Created
33 // 15 Jul 2003 Alfonso Mantero "Detecto 35 // 15 Jul 2003 Alfonso Mantero "DetectorType" use integration
34 // 23 Sep 2003 Alfonso Mantero differnt 36 // 23 Sep 2003 Alfonso Mantero differnt geometries integration
35 // 37 //
36 // ------------------------------------------- 38 // -------------------------------------------------------------------
37 39
38 #include "XrayFluoEventAction.hh" 40 #include "XrayFluoEventAction.hh"
39 #include "XrayFluoSensorHit.hh" 41 #include "XrayFluoSensorHit.hh"
40 #include "XrayFluoEventActionMessenger.hh" 42 #include "XrayFluoEventActionMessenger.hh"
>> 43
>> 44 //#include "XrayFluoRunAction.hh"
41 #include "XrayFluoDataSet.hh" 45 #include "XrayFluoDataSet.hh"
>> 46
42 #include "XrayFluoAnalysisManager.hh" 47 #include "XrayFluoAnalysisManager.hh"
43 48
44 #include "G4SystemOfUnits.hh" <<
45 #include "G4Event.hh" 49 #include "G4Event.hh"
46 #include "G4EventManager.hh" 50 #include "G4EventManager.hh"
47 #include "G4HCofThisEvent.hh" 51 #include "G4HCofThisEvent.hh"
48 #include "G4VHitsCollection.hh" 52 #include "G4VHitsCollection.hh"
>> 53 #include "G4TrajectoryContainer.hh"
>> 54 #include "G4Trajectory.hh"
>> 55 #include "G4VVisManager.hh"
49 #include "G4SDManager.hh" 56 #include "G4SDManager.hh"
50 #include "G4UImanager.hh" 57 #include "G4UImanager.hh"
51 #include "G4ios.hh" 58 #include "G4ios.hh"
52 #include "G4UnitsTable.hh" 59 #include "G4UnitsTable.hh"
53 #include "Randomize.hh" 60 #include "Randomize.hh"
54 #include <fstream> 61 #include <fstream>
55 <<
56 //....oooOO0OOooo........oooOO0OOooo........oo 62 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
57 63
58 XrayFluoEventAction::XrayFluoEventAction(const << 64 XrayFluoEventAction::XrayFluoEventAction(XrayFluoDetectorConstruction* det)
59 :drawFlag("all"), 65 :drawFlag("all"),
60 HPGeCollID(0), 66 HPGeCollID(0),
61 eventMessenger(0), 67 eventMessenger(0),
62 printModulo(1), 68 printModulo(1),
63 detectorType(0) 69 detectorType(0)
64 { 70 {
65 eventMessenger = new XrayFluoEventActionMess 71 eventMessenger = new XrayFluoEventActionMessenger(this);
66 72
67 if (!(det->GetPhaseSpaceFlag()) ){ 73 if (!(det->GetPhaseSpaceFlag()) ){
68 detectorType = det->GetDetectorType(); 74 detectorType = det->GetDetectorType();
69 HPGeCollID=-1; 75 HPGeCollID=-1;
70 } 76 }
>> 77
>> 78 //runManager = new XrayFluoRunAction();
>> 79 G4cout << "XrayFluoEventAction created" << G4endl;
71 } 80 }
72 81
73 //....oooOO0OOooo........oooOO0OOooo........oo 82 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
74 83
75 XrayFluoEventAction::XrayFluoEventAction(const << 84 XrayFluoEventAction::XrayFluoEventAction(XrayFluoPlaneDetectorConstruction* det)
76 :drawFlag("all"), 85 :drawFlag("all"),
77 HPGeCollID(-1), 86 HPGeCollID(-1),
78 eventMessenger(0), 87 eventMessenger(0),
79 printModulo(1), 88 printModulo(1),
80 detectorType(0) 89 detectorType(0)
81 { 90 {
82 eventMessenger = new XrayFluoEventActionMess 91 eventMessenger = new XrayFluoEventActionMessenger(this);
83 detectorType = det->GetDetectorType(); 92 detectorType = det->GetDetectorType();
>> 93
>> 94 //runManager = new XrayFluoRunAction();
>> 95 G4cout << "XrayFluoEventAction created" << G4endl;
84 } 96 }
85 97
86 //....oooOO0OOooo........oooOO0OOooo........oo 98 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
87 99
88 XrayFluoEventAction::XrayFluoEventAction(const << 100 XrayFluoEventAction::XrayFluoEventAction(XrayFluoMercuryDetectorConstruction* det)
89 :drawFlag("all"), 101 :drawFlag("all"),
90 HPGeCollID(-1), 102 HPGeCollID(-1),
91 eventMessenger(0), 103 eventMessenger(0),
92 printModulo(1), 104 printModulo(1),
93 detectorType(0) 105 detectorType(0)
94 { 106 {
95 eventMessenger = new XrayFluoEventActionMess 107 eventMessenger = new XrayFluoEventActionMessenger(this);
96 detectorType = det->GetDetectorType(); 108 detectorType = det->GetDetectorType();
>> 109
>> 110 //runManager = new XrayFluoRunAction();
>> 111 G4cout << "XrayFluoEventAction created" << G4endl;
97 } 112 }
98 113
99 //....oooOO0OOooo........oooOO0OOooo........oo 114 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
100 115
>> 116
>> 117
>> 118 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
>> 119
101 XrayFluoEventAction::~XrayFluoEventAction() 120 XrayFluoEventAction::~XrayFluoEventAction()
102 { 121 {
103 delete eventMessenger; 122 delete eventMessenger;
104 eventMessenger = 0; 123 eventMessenger = 0;
>> 124 //delete runManager;
>> 125 //runManager = 0;
>> 126 G4cout << "XrayFluoEventAction deleted" << G4endl;
105 } 127 }
106 128
107 //....oooOO0OOooo........oooOO0OOooo........oo 129 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
108 130
109 void XrayFluoEventAction::BeginOfEventAction(c 131 void XrayFluoEventAction::BeginOfEventAction(const G4Event* evt)
110 { 132 {
111 133
112 G4int eventNumber = (evt->GetEventID())+1; 134 G4int eventNumber = (evt->GetEventID())+1;
113 if ( eventNumber == 1){ 135 if ( eventNumber == 1){
114 136
115 G4cout << "# = 100000 events" << G4endl; 137 G4cout << "# = 100000 events" << G4endl;
116 G4cout << "1--------+---------+---------+--- 138 G4cout << "1--------+---------+---------+---------+---------5e6"<<G4endl;
117 } 139 }
118 140
119 if ( ((eventNumber) % 100000) == 0 ) { 141 if ( ((eventNumber) % 100000) == 0 ) {
120 142
121 if ( eventNumber % (G4int)5e6 != 0 ) G4cou 143 if ( eventNumber % (G4int)5e6 != 0 ) G4cout << "#" << std::flush;
122 else G4cout << "#"<< G4endl; << 144 else G4cout << "#"<< G4endl;
>> 145 // if ( eventNumber % 5e6 == 0 ) G4cout << "#"<< G4endl;
>> 146 #ifdef G4ANALYSIS_USE
>> 147 XrayFluoAnalysisManager* analysis = XrayFluoAnalysisManager::getInstance();
>> 148 analysis->PlotCurrentResults();
>> 149 #endif
123 } 150 }
124 151
125 if (HPGeCollID==-1) << 152 if (HPGeCollID==-1)
>> 153
126 { 154 {
127 G4SDManager * SDman = G4SDManager::GetSD 155 G4SDManager * SDman = G4SDManager::GetSDMpointer();
128 HPGeCollID = SDman->GetCollectionID("HPG 156 HPGeCollID = SDman->GetCollectionID("HPGeCollection");
129 //the pointer points to the ID number of 157 //the pointer points to the ID number of the sensitive detector
130 } 158 }
131 } 159 }
132 160
133 //....oooOO0OOooo........oooOO0OOooo........oo 161 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
134 162
135 void XrayFluoEventAction::EndOfEventAction(con 163 void XrayFluoEventAction::EndOfEventAction(const G4Event* evt)
136 { 164 {
>> 165
137 if (detectorType) { 166 if (detectorType) {
138 167
139 // extracted from hits, compute the total 168 // extracted from hits, compute the total energy deposit (and total charged
140 // track length) 169 // track length)
141 G4HCofThisEvent* HCE = evt->GetHCofThisEve 170 G4HCofThisEvent* HCE = evt->GetHCofThisEvent();
142 171
143 XrayFluoSensorHitsCollection* HPGeHC = 0; 172 XrayFluoSensorHitsCollection* HPGeHC = 0;
144 G4int n_hit = 0; 173 G4int n_hit = 0;
145 G4double totEnergy=0., energyD=0.; << 174 G4double totEnergyDetect=0., totEnergy=0., energyD=0.;
146 175
147 if (HCE) HPGeHC = << 176 if (HCE) HPGeHC = (XrayFluoSensorHitsCollection*)(HCE->GetHC(HPGeCollID));
148 (XrayFluoSensorHitsCollection*)(HCE->Get << 177 if(HPGeHC)
149 << 178
150 if(HPGeHC) << 179
151 { 180 {
152 n_hit = HPGeHC->entries(); 181 n_hit = HPGeHC->entries();
153 << 182
>> 183 // if (n_hit) {G4cout << "Ecco quante hit ho nel detector "<< n_hit << G4endl;}
>> 184
154 for (G4int i=0;i<n_hit;i++) 185 for (G4int i=0;i<n_hit;i++)
155 { 186 {
>> 187
156 totEnergy += (*HPGeHC)[i]->GetEdepTot(); 188 totEnergy += (*HPGeHC)[i]->GetEdepTot();
157 energyD = detectorType->ResponseFunction << 189
>> 190
>> 191
>> 192 energyD = detectorType->ResponseFunction(totEnergy);
>> 193 // energyD = totEnergy;
>> 194 // G4cout << "energy deposit: "<< totEnergy << G4endl;
>> 195 #ifdef G4ANALYSIS_USE
158 XrayFluoAnalysisManager* analysis = Xray 196 XrayFluoAnalysisManager* analysis = XrayFluoAnalysisManager::getInstance();
159 analysis->analyseEnergyDep(energyD); 197 analysis->analyseEnergyDep(energyD);
>> 198 #endif
>> 199 totEnergyDetect += energyD;
>> 200
>> 201
160 } 202 }
161 } 203 }
162 } 204 }
>> 205
>> 206 // extract the trajectories and draw them
>> 207
>> 208 if (G4VVisManager::GetConcreteInstance())
>> 209 {
>> 210
>> 211 G4TrajectoryContainer * trajectoryContainer = evt->GetTrajectoryContainer();
>> 212 G4int n_trajectories = 0;
>> 213 if (trajectoryContainer) n_trajectories = trajectoryContainer->size();
>> 214
>> 215 for (G4int i=0; i<n_trajectories; i++)
>> 216 { G4Trajectory* trj = (G4Trajectory*)((*(evt->GetTrajectoryContainer()))[i]);
>> 217 if (drawFlag == "all") trj->DrawTrajectory(50);
>> 218 else if ((drawFlag == "charged")&&(trj->GetCharge() != 0.))
>> 219 trj->DrawTrajectory(50);
>> 220 else if ((drawFlag == "neutral")&&(trj->GetCharge() == 0.))
>> 221 trj->DrawTrajectory(50);
>> 222 }
>> 223 }
163 } 224 }
164 225
165 //....oooOO0OOooo........oooOO0OOooo........oo 226 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
166 227
167 228
>> 229
168 G4double XrayFluoEventAction::RandomCut(G4doub 230 G4double XrayFluoEventAction::RandomCut(G4double energy)
169 231
170 { 232 {
171 G4double efficiency = 1.; 233 G4double efficiency = 1.;
172 G4double F = 0.15; 234 G4double F = 0.15;
173 G4double epsilon = 2.96 * eV; 235 G4double epsilon = 2.96 * eV;
174 G4double deltaE = 220 * eV; 236 G4double deltaE = 220 * eV;
175 G4double EdepDetect = 0.; 237 G4double EdepDetect = 0.;
176 << 238 //const XrayFluoDataSet* dataSet = runManager->GetEfficiencySet();
>> 239
>> 240 //G4double id = 0;
>> 241
>> 242 //efficiency = dataSet->FindValue(energy,id);
>> 243
177 G4double Random = G4UniformRand(); 244 G4double Random = G4UniformRand();
178 245
179 if ( Random<efficiency ) << 246 if ( Random<efficiency )
180 { << 247 {
181 G4double sigma = std::sqrt(F*epsilon*ene << 248 G4double sigma = std::sqrt(F*epsilon*energy+std::pow(deltaE/2355,2));
182 EdepDetect = G4RandGauss::shoot(energy, <<
183 <<
184 } <<
185 else {EdepDetect = 0.;} <<
186 249
187 return EdepDetect; << 250 EdepDetect = G4RandGauss::shoot(energy, sigma );
188 } << 251
>> 252 }
>> 253 else {EdepDetect = 0.;}
>> 254 return EdepDetect;
>> 255
>> 256 };
189 257
190 258
191 259