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