source: trunk/examples/extended/electromagnetic/TestEm3/src/SteppingAction.cc@ 1036

Last change on this file since 1036 was 807, checked in by garnier, 17 years ago

update
File size: 5.8 KB
Line 
1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26// $Id: SteppingAction.cc,v 1.25 2006/06/29 16:53:23 gunter Exp $
27// GEANT4 tag $Name: geant4-09-01-patch-02 $
28//
29//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
30//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
31
32#include "SteppingAction.hh"
33
34#include "DetectorConstruction.hh"
35#include "RunAction.hh"
36#include "EventAction.hh"
37#include "HistoManager.hh"
38
39#include "G4Step.hh"
40#include "G4Positron.hh"
41#include "G4RunManager.hh"
42
43//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
44
45SteppingAction::SteppingAction(DetectorConstruction* det, RunAction* run,
46 EventAction* evt, HistoManager* hist)
47:G4UserSteppingAction(),detector(det),runAct(run),eventAct(evt),
48 histoManager(hist)
49{}
50
51//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
52
53SteppingAction::~SteppingAction()
54{}
55
56//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
57
58void SteppingAction::UserSteppingAction(const G4Step* aStep)
59{
60 //track informations
61 const G4StepPoint* prePoint = aStep->GetPreStepPoint();
62 const G4StepPoint* endPoint = aStep->GetPostStepPoint();
63 const G4ParticleDefinition* particle = aStep->GetTrack()->GetDefinition();
64
65 //if World, return
66 //
67 G4VPhysicalVolume* volume = prePoint->GetTouchableHandle()->GetVolume();
68 //if sum of absorbers do not fill exactly a layer: check material, not volume.
69 G4Material* mat = volume->GetLogicalVolume()->GetMaterial();
70 if (mat == detector->GetWorldMaterial()) return;
71
72 //here we are in an absorber. Locate it
73 //
74 G4int absorNum = prePoint->GetTouchableHandle()->GetCopyNumber(0);
75 G4int layerNum = prePoint->GetTouchableHandle()->GetCopyNumber(1);
76
77 // collect energy deposit
78 G4double edep = aStep->GetTotalEnergyDeposit();
79
80 // collect step length of charged particles
81 G4double stepl = 0.;
82 if (particle->GetPDGCharge() != 0.) stepl = aStep->GetStepLength();
83
84 // sum up per event
85 eventAct->SumEnergy(absorNum,edep,stepl);
86
87 //longitudinal profile of edep per absorber
88 if (edep>0.) histoManager->FillHisto(MaxAbsor+absorNum, layerNum+1., edep);
89
90 //energy flow
91 //
92 // unique identificator of layer+absorber
93 G4int Idnow = (detector->GetNbOfAbsor())*layerNum + absorNum;
94 G4int plane;
95 //
96 //leaving the absorber ?
97 if (endPoint->GetStepStatus() == fGeomBoundary) {
98 G4ThreeVector position = endPoint->GetPosition();
99 G4ThreeVector direction = endPoint->GetMomentumDirection();
100 G4double sizeYZ = 0.5*detector->GetCalorSizeYZ();
101 G4double Eflow = endPoint->GetKineticEnergy();
102 if (particle == G4Positron::Positron()) Eflow += 2*electron_mass_c2;
103 if ((std::abs(position.y()) >= sizeYZ) || (std::abs(position.z()) >= sizeYZ))
104 runAct->sumLateralEleak(Idnow, Eflow);
105 else if (direction.x() >= 0.) runAct->sumEnergyFlow(plane=Idnow+1, Eflow);
106 else runAct->sumEnergyFlow(plane=Idnow, -Eflow);
107 }
108
109//// example of Birk attenuation
110//// G4double destep = aStep->GetTotalEnergyDeposit();
111//// G4double response = BirkAttenuation(aStep);
112//// G4cout << " Destep: " << destep/keV << " keV"
113//// << " response after Birk: " << response/keV << " keV" << G4endl;
114}
115
116//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
117
118G4double SteppingAction::BirkAttenuation(const G4Step* aStep)
119{
120 //Example of Birk attenuation law in organic scintillators.
121 //adapted from Geant3 PHYS337. See MIN 80 (1970) 239-244
122 //
123 const G4String myMaterial = "Scintillator";
124 const G4double birk1 = 0.013*g/(MeV*cm2);
125 //
126 G4double destep = aStep->GetTotalEnergyDeposit();
127 G4Material* material = aStep->GetTrack()->GetMaterial();
128 G4double charge = aStep->GetTrack()->GetDefinition()->GetPDGCharge();
129 //
130 G4double response = destep;
131 if ((material->GetName()==myMaterial)&&(charge!=0.))
132 {
133 G4double correction =
134 birk1*destep/((material->GetDensity())*(aStep->GetStepLength()));
135 response = destep/(1. + correction);
136 }
137 return response;
138}
139
140//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
141
Note: See TracBrowser for help on using the repository browser.