| 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 // $Id: SteppingAction.cc,v 1.9 2009-09-18 17:34:54 maire Exp $
>> 27 // GEANT4 tag $Name: geant4-09-04-patch-02 $
26 // 28 //
27 //....oooOO0OOooo........oooOO0OOooo........oo 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
28 //....oooOO0OOooo........oooOO0OOooo........oo 30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
29 31
30 #include "SteppingAction.hh" 32 #include "SteppingAction.hh"
31 33
32 #include "DetectorConstruction.hh" 34 #include "DetectorConstruction.hh"
>> 35 #include "RunAction.hh"
>> 36 #include "PrimaryGeneratorAction.hh"
33 #include "EventAction.hh" 37 #include "EventAction.hh"
>> 38 #include "HistoManager.hh"
34 39
35 #include "G4Step.hh" 40 #include "G4Step.hh"
36 41
>> 42 #include "G4Geantino.hh"
>> 43
37 //....oooOO0OOooo........oooOO0OOooo........oo 44 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
38 45
39 SteppingAction::SteppingAction(DetectorConstru << 46 SteppingAction::SteppingAction(DetectorConstruction* det, RunAction* run,
40 :G4UserSteppingAction(),detector(det),eventAct << 47 PrimaryGeneratorAction* prim, EventAction* evt,
>> 48 HistoManager* hist)
>> 49 :G4UserSteppingAction(),detector(det),runAct(run),primary(prim),eventAct(evt),
>> 50 histoManager(hist)
41 { 51 {
42 first = true; 52 first = true;
43 lvol_world = lvol_module = lvol_layer = lvol 53 lvol_world = lvol_module = lvol_layer = lvol_fiber = 0;
44 } 54 }
45 55
46 //....oooOO0OOooo........oooOO0OOooo........oo 56 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
47 57
48 SteppingAction::~SteppingAction() 58 SteppingAction::~SteppingAction()
49 {} 59 {}
50 60
51 //....oooOO0OOooo........oooOO0OOooo........oo 61 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
52 62
53 void SteppingAction::UserSteppingAction(const 63 void SteppingAction::UserSteppingAction(const G4Step* step )
54 { 64 {
55 //some initialisation 65 //some initialisation
56 // 66 //
57 if (first) { 67 if (first) {
58 lvol_world = detector->GetLvolWorld(); 68 lvol_world = detector->GetLvolWorld();
59 lvol_module = detector->GetLvolModule(); 69 lvol_module = detector->GetLvolModule();
60 lvol_layer = detector->GetLvolLayer(); 70 lvol_layer = detector->GetLvolLayer();
61 lvol_fiber = detector->GetLvolFiber(); << 71 lvol_fiber = detector->GetLvolFiber();
>> 72
>> 73 calorThickness = detector->GetCalorThickness();
>> 74 calorSizeYZ = detector->GetCalorSizeYZ();
>> 75 moduleThickness = detector->GetModuleThickness();
>> 76 d1Pixel = detector->GetD1Pixel();
>> 77 d2Pixel = detector->GetD2Pixel();
>> 78 n1pxl = detector->GetN1Pixels();
>> 79 n2pxl = detector->GetN2Pixels();
>> 80 n1shift = detector->GetN1Shift();
>> 81
62 first = false; 82 first = false;
63 } 83 }
>> 84
>> 85 //locate point in geometry
>> 86 //
>> 87 G4TouchableHandle touch1 = step->GetPreStepPoint()->GetTouchableHandle();
>> 88 G4LogicalVolume* lvol = touch1->GetVolume()->GetLogicalVolume();
>> 89
>> 90 //if world, return
>> 91 //
>> 92 if (lvol == lvol_world) return;
64 93
>> 94 //sum nb of radiation length of calorimeter with geantino
>> 95 //
>> 96 G4ParticleDefinition* particle = step->GetTrack()->GetDefinition();
>> 97 if (particle == G4Geantino::Geantino()) {
>> 98 G4double radl = lvol->GetMaterial()->GetRadlen();
>> 99 G4double stepl = step->GetStepLength();
>> 100 eventAct->SumNbRadLength(stepl/radl);
>> 101 }
>> 102
65 //if no edep, return 103 //if no edep, return
66 // 104 //
67 G4double edep = step->GetTotalEnergyDeposit() 105 G4double edep = step->GetTotalEnergyDeposit();
68 ///if (edep == 0.) return; << 106 if (edep == 0.) return;
69 107
70 //locate point in geometry << 108 //locate position and compute pixel number
71 // 109 //
72 G4int iModule = 0; << 110 G4ThreeVector point1 = step->GetPreStepPoint()->GetPosition();
73 G4int iLayer = 0; << 111 G4int i1Module = (int) ((point1.x() + 0.5*calorThickness)/moduleThickness);
74 G4int iFiber = 0; << 112 G4int i1Pixel = (int) ((point1.x() + 0.5*calorThickness)/d1Pixel);
75 << 113 if (i1Pixel < 0 ) i1Pixel = 0;
76 G4TouchableHandle touch1 = step->GetPreStepPo << 114 if (i1Pixel >= n1pxl) i1Pixel = n1pxl - 1;
77 G4LogicalVolume* lvol = touch1->GetVolume()-> << 115 G4double point1yz = point1.y();
>> 116 if (i1Module%2 != 0) point1yz = point1.z();
>> 117 G4int i2Pixel = (int) ((point1yz + 0.5*calorSizeYZ)/d2Pixel);
>> 118 if (i2Pixel < 0 ) i2Pixel = 0;
>> 119 if (i2Pixel >= n2pxl) i2Pixel = n2pxl - 1;
78 120
79 if (lvol == lvol_world) return; << 121 G4int iPixel = i1Pixel*n1shift + i2Pixel;
80 else if (lvol == lvol_module) { iModule = tou << 122
81 else if (lvol == lvol_layer) { iLayer = tou << 123 // sum total energy deposit
82 iModule = tou <<
83 else if (lvol == lvol_fiber) { iFiber = tou <<
84 iLayer = tou <<
85 iModule = tou <<
86 <<
87 // sum edep <<
88 // 124 //
89 eventAct->SumDeStep(iModule, iLayer, iFiber, << 125 eventAct->SumTotalEnergy(iPixel, edep);
>> 126
>> 127 //in fiber ?
>> 128 //
>> 129 if (lvol == lvol_fiber) {
>> 130 eventAct->SumVisibleEnergy(iPixel, edep);
>> 131 }
90 } 132 }
91 133
92 //....oooOO0OOooo........oooOO0OOooo........oo 134 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
93 135
94 G4double SteppingAction::BirksAttenuation(cons 136 G4double SteppingAction::BirksAttenuation(const G4Step* aStep)
95 { 137 {
96 //Example of Birk attenuation law in organic 138 //Example of Birk attenuation law in organic scintillators.
97 //adapted from Geant3 PHYS337. See MIN 80 (19 139 //adapted from Geant3 PHYS337. See MIN 80 (1970) 239-244
98 // 140 //
99 const G4Material* material = aStep->GetTrack( << 141 G4Material* material = aStep->GetTrack()->GetMaterial();
100 G4double birk1 = material->GetIonisatio 142 G4double birk1 = material->GetIonisation()->GetBirksConstant();
101 G4double destep = aStep->GetTotalEnergyD 143 G4double destep = aStep->GetTotalEnergyDeposit();
102 G4double stepl = aStep->GetStepLength() 144 G4double stepl = aStep->GetStepLength();
103 G4double charge = aStep->GetTrack()->Get 145 G4double charge = aStep->GetTrack()->GetDefinition()->GetPDGCharge();
104 // 146 //
105 G4double response = destep; 147 G4double response = destep;
106 if (birk1*destep*stepl*charge != 0.) 148 if (birk1*destep*stepl*charge != 0.)
107 { 149 {
108 response = destep/(1. + birk1*destep/step 150 response = destep/(1. + birk1*destep/stepl);
109 } 151 }
110 return response; 152 return response;
111 } 153 }
112 154
113 //....oooOO0OOooo........oooOO0OOooo........oo 155 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
114 156
115 157