source: trunk/examples/advanced/raredecay_calorimetry/src/PhotInStackingAction.cc @ 1288

Last change on this file since 1288 was 807, checked in by garnier, 16 years ago

update
File size: 12.4 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
27//#define debug
28//#define errors
29
30#include "PhotInStackingAction.hh"
31
32// Can be moved to the PhotInConstants.hh
33G4int PhotInStackingAction::PhotInNGamma     [PhotInDiNSections] = {0,0,0,0,0,0};
34G4int PhotInStackingAction::PhotInNElectron  [PhotInDiNSections] = {0,0,0,0,0,0};
35G4int PhotInStackingAction::PhotInNPositron  [PhotInDiNSections] = {0,0,0,0,0,0};
36G4int PhotInStackingAction::PhotInNNeutrons  [PhotInDiNSections] = {0,0,0,0,0,0};
37G4int PhotInStackingAction::PhotInNProtons   [PhotInDiNSections] = {0,0,0,0,0,0};
38G4int PhotInStackingAction::PhotInNDeuterons [PhotInDiNSections] = {0,0,0,0,0,0};
39G4int PhotInStackingAction::PhotInNTritons   [PhotInDiNSections] = {0,0,0,0,0,0};
40G4int PhotInStackingAction::PhotInNHe3s      [PhotInDiNSections] = {0,0,0,0,0,0};
41G4int PhotInStackingAction::PhotInNAlphas    [PhotInDiNSections] = {0,0,0,0,0,0};
42G4int PhotInStackingAction::PhotInNLambdas   [PhotInDiNSections] = {0,0,0,0,0,0};
43G4int PhotInStackingAction::PhotInNHeavyFrags[PhotInDiNSections] = {0,0,0,0,0,0};
44G4int PhotInStackingAction::PhotInNMesons    [PhotInDiNSections] = {0,0,0,0,0,0};
45G4double PhotInStackingAction::PhotInEMinGamma     [PhotInDiNSections] = {DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX};
46G4double PhotInStackingAction::PhotInEMinElectron  [PhotInDiNSections] = {DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX};
47G4double PhotInStackingAction::PhotInEMinPositron  [PhotInDiNSections] = {DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX};
48G4double PhotInStackingAction::PhotInEMinNeutrons  [PhotInDiNSections] = {DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX};
49G4double PhotInStackingAction::PhotInEMinProtons   [PhotInDiNSections] = {DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX};
50G4double PhotInStackingAction::PhotInEMinDeuterons [PhotInDiNSections] = {DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX};
51G4double PhotInStackingAction::PhotInEMinTritons   [PhotInDiNSections] = {DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX};
52G4double PhotInStackingAction::PhotInEMinHe3s      [PhotInDiNSections] = {DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX};
53G4double PhotInStackingAction::PhotInEMinAlphas    [PhotInDiNSections] = {DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX};
54G4double PhotInStackingAction::PhotInEMinLambdas   [PhotInDiNSections] = {DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX};
55G4double PhotInStackingAction::PhotInEMinHeavyFrags[PhotInDiNSections] = {DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX};
56G4double PhotInStackingAction::PhotInEMinMesons    [PhotInDiNSections] = {DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX,DBL_MAX};
57G4double PhotInStackingAction::PhotInEMaxGamma     [PhotInDiNSections]={0.,0.,0.,0.,0.,0.};
58G4double PhotInStackingAction::PhotInEMaxElectron  [PhotInDiNSections]={0.,0.,0.,0.,0.,0.};
59G4double PhotInStackingAction::PhotInEMaxPositron  [PhotInDiNSections]={0.,0.,0.,0.,0.,0.};
60G4double PhotInStackingAction::PhotInEMaxNeutrons  [PhotInDiNSections]={0.,0.,0.,0.,0.,0.};
61G4double PhotInStackingAction::PhotInEMaxProtons   [PhotInDiNSections]={0.,0.,0.,0.,0.,0.};
62G4double PhotInStackingAction::PhotInEMaxDeuterons [PhotInDiNSections]={0.,0.,0.,0.,0.,0.};
63G4double PhotInStackingAction::PhotInEMaxTritons   [PhotInDiNSections]={0.,0.,0.,0.,0.,0.};
64G4double PhotInStackingAction::PhotInEMaxHe3s      [PhotInDiNSections]={0.,0.,0.,0.,0.,0.};
65G4double PhotInStackingAction::PhotInEMaxAlphas    [PhotInDiNSections]={0.,0.,0.,0.,0.,0.};
66G4double PhotInStackingAction::PhotInEMaxLambdas   [PhotInDiNSections]={0.,0.,0.,0.,0.,0.};
67G4double PhotInStackingAction::PhotInEMaxHeavyFrags[PhotInDiNSections]={0.,0.,0.,0.,0.,0.};
68G4double PhotInStackingAction::PhotInEMaxMesons    [PhotInDiNSections]={0.,0.,0.,0.,0.,0.};
69
70PhotInStackingAction::PhotInStackingAction() {}
71
72PhotInStackingAction::~PhotInStackingAction(){}
73
74G4ClassificationOfNewTrack PhotInStackingAction::ClassifyNewTrack(const G4Track* aTrack)
75{
76  G4TouchableHistory* theTouchable = (G4TouchableHistory*)(aTrack->GetTouchable());
77  G4int depth=0;                  // Depth of geometrical hierarchy of volumes
78  if(theTouchable)depth=theTouchable->GetHistoryDepth();
79#ifdef debug
80  G4cout<<"PhotInStackingAction::ClassifyNewTrack: --------->>>>>>> depth="<<depth<<G4endl;
81#endif
82#ifdef debug
83                G4cout<<"PhotInStackingAction::ClassifyNewTrack: XYZ="<<aTrack->GetPosition()<<G4endl;
84#endif
85  if(depth>1)                     // It is in one of the Sections, not in the Mother volume
86  {
87    G4VPhysicalVolume* calPhys = theTouchable->GetVolume(depth-1);// PhysVolume of Sections
88    if(calPhys)
89    {
90      G4String sectName=calPhys->GetName();
91#ifdef debug
92      G4cout<<"PhotInStackingAction::ClassifyNewTrack: SectionName="<<sectName<<G4endl;
93#endif
94      G4int i = -1;                // 2 * Section_number
95      if     (sectName==PhotInCalName[0]) { i = 0; }
96      else if(sectName==PhotInCalName[1]) { i = 2; }
97      else if(sectName==PhotInCalName[2]) { i = 4; }
98#ifdef errors
99      if(i<0) G4cerr<<"*PhotInStackingAction::ClassifyNewTrack: sensitive i="<<i<<G4endl;
100#endif
101#ifdef debug
102      G4int layerNumber=theTouchable->GetReplicaNumber(depth-2);//WhatIfLayers notFill Sect
103      G4cout<<"PhotInStackingAction::ClassifyNewTrack: Layer#"<<layerNumber<<G4endl;
104#endif
105      G4int slabNumber=-1;        // Prototype of the slab number (-1 means absorber)
106      if(depth>2) slabNumber=theTouchable->GetReplicaNumber(depth-3);
107      else i++;                   // Collect information in the absorber
108#ifdef debug
109      G4cout<<"PhotInStackingAction::ClassifyNewTrack:Slab#"<<slabNumber<<",i="<<i<<G4endl;
110#endif
111      G4ParticleDefinition* particleType = aTrack->GetDefinition();
112      if(particleType==G4Gamma::GammaDefinition())
113      {
114        PhotInNGamma[i]++;
115        if(PhotInEMinGamma[i]>aTrack->GetKineticEnergy())
116          PhotInEMinGamma[i]=aTrack->GetKineticEnergy();
117        if(PhotInEMaxGamma[i]<aTrack->GetKineticEnergy())
118          PhotInEMaxGamma[i]=aTrack->GetKineticEnergy();
119      }
120      else if(particleType==G4Electron::ElectronDefinition())
121      {
122        PhotInNElectron[i]++;
123        if(PhotInEMinElectron[i]>aTrack->GetKineticEnergy())
124          PhotInEMinElectron[i]=aTrack->GetKineticEnergy();
125        if(PhotInEMaxElectron[i]<aTrack->GetKineticEnergy())
126          PhotInEMaxElectron[i]=aTrack->GetKineticEnergy();
127      }
128      else if(particleType==G4Positron::PositronDefinition())
129      {
130        PhotInNPositron[i]++;
131        if(PhotInEMinPositron[i]>aTrack->GetKineticEnergy())
132          PhotInEMinPositron[i]=aTrack->GetKineticEnergy();
133        if(PhotInEMaxPositron[i]<aTrack->GetKineticEnergy())
134          PhotInEMaxPositron[i]=aTrack->GetKineticEnergy();
135      }
136      else if(particleType==G4Neutron::NeutronDefinition())
137      {
138        PhotInNNeutrons[i]++;
139        if(PhotInEMinNeutrons[i]>aTrack->GetKineticEnergy())
140          PhotInEMinNeutrons[i]=aTrack->GetKineticEnergy();
141        if(PhotInEMaxNeutrons[i]<aTrack->GetKineticEnergy())
142          PhotInEMaxNeutrons[i]=aTrack->GetKineticEnergy();
143      }
144      else if(particleType==G4Proton::ProtonDefinition())
145      {
146        PhotInNProtons[i]++;
147        if(PhotInEMinProtons[i]>aTrack->GetKineticEnergy())
148          PhotInEMinProtons[i]=aTrack->GetKineticEnergy();
149        if(PhotInEMaxProtons[i]<aTrack->GetKineticEnergy())
150          PhotInEMaxProtons[i]=aTrack->GetKineticEnergy();
151      }
152      else if(particleType==G4Deuteron::DeuteronDefinition())
153      {
154        PhotInNDeuterons[i]++;
155        if(PhotInEMinDeuterons[i]>aTrack->GetKineticEnergy())
156          PhotInEMinDeuterons[i]=aTrack->GetKineticEnergy();
157        if(PhotInEMaxDeuterons[i]<aTrack->GetKineticEnergy())
158          PhotInEMaxDeuterons[i]=aTrack->GetKineticEnergy();
159      }
160      else if(particleType==G4Triton::TritonDefinition())
161      {
162        PhotInNTritons[i]++;
163        if(PhotInEMinTritons[i]>aTrack->GetKineticEnergy())
164          PhotInEMinTritons[i]=aTrack->GetKineticEnergy();
165        if(PhotInEMaxTritons[i]<aTrack->GetKineticEnergy())
166          PhotInEMaxTritons[i]=aTrack->GetKineticEnergy();
167      }
168      else if(particleType==G4He3::He3Definition())
169      {
170        PhotInNHe3s[i]++;
171        if(PhotInEMinHe3s[i]>aTrack->GetKineticEnergy())
172          PhotInEMinHe3s[i]=aTrack->GetKineticEnergy();
173        if(PhotInEMaxHe3s[i]<aTrack->GetKineticEnergy())
174          PhotInEMaxHe3s[i]=aTrack->GetKineticEnergy();
175      }
176      else if(particleType==G4Alpha::AlphaDefinition())
177      {
178        PhotInNAlphas[i]++;
179        if(PhotInEMinAlphas[i]>aTrack->GetKineticEnergy())
180          PhotInEMinAlphas[i]=aTrack->GetKineticEnergy();
181        if(PhotInEMaxAlphas[i]<aTrack->GetKineticEnergy())
182          PhotInEMaxAlphas[i]=aTrack->GetKineticEnergy();
183      }
184      else if(particleType==G4Lambda::LambdaDefinition())
185      {
186        PhotInNLambdas[i]++;
187        if(PhotInEMinLambdas[i]>aTrack->GetKineticEnergy())
188          PhotInEMinLambdas[i]=aTrack->GetKineticEnergy();
189        if(PhotInEMaxLambdas[i]<aTrack->GetKineticEnergy())
190          PhotInEMaxLambdas[i]=aTrack->GetKineticEnergy();
191      }
192      else if(particleType->GetBaryonNumber()>4)
193      {
194        PhotInNHeavyFrags[i]++;
195        if(PhotInEMinHeavyFrags[i]>aTrack->GetKineticEnergy())
196          PhotInEMinHeavyFrags[i]=aTrack->GetKineticEnergy();
197        if(PhotInEMaxHeavyFrags[i]<aTrack->GetKineticEnergy())
198          PhotInEMaxHeavyFrags[i]=aTrack->GetKineticEnergy();
199      }
200      else if(!particleType->GetBaryonNumber())
201      {
202        PhotInNMesons[i]++;
203        if(PhotInEMinMesons[i]>aTrack->GetKineticEnergy())
204          PhotInEMinMesons[i]=aTrack->GetKineticEnergy();
205        if(PhotInEMaxMesons[i]<aTrack->GetKineticEnergy())
206          PhotInEMaxMesons[i]=aTrack->GetKineticEnergy();
207      }
208    }
209  }
210  return fUrgent;
211}
212
213void PhotInStackingAction::PrepareNewEvent()
214{ 
215  for(int i=0; i<PhotInDiNSections ;i++)
216  {
217    PhotInNGamma     [i] = 0; 
218    PhotInNElectron  [i] = 0;
219    PhotInNPositron  [i] = 0;
220    PhotInNNeutrons  [i] = 0; 
221    PhotInNProtons   [i] = 0;
222    PhotInNDeuterons [i] = 0;
223    PhotInNTritons   [i] = 0; 
224    PhotInNHe3s      [i] = 0;
225    PhotInNAlphas    [i] = 0;
226    PhotInNLambdas   [i] = 0; 
227    PhotInNHeavyFrags[i] = 0;
228    PhotInNMesons    [i] = 0;
229
230    PhotInEMinGamma     [i] = DBL_MAX;
231    PhotInEMinElectron  [i] = DBL_MAX;
232    PhotInEMinPositron  [i] = DBL_MAX;
233    PhotInEMinNeutrons  [i] = DBL_MAX;
234    PhotInEMinProtons   [i] = DBL_MAX;
235    PhotInEMinDeuterons [i] = DBL_MAX;
236    PhotInEMinTritons   [i] = DBL_MAX;
237    PhotInEMinHe3s      [i] = DBL_MAX;
238    PhotInEMinAlphas    [i] = DBL_MAX;
239    PhotInEMinLambdas   [i] = DBL_MAX;
240    PhotInEMinHeavyFrags[i] = DBL_MAX;
241    PhotInEMinMesons    [i] = DBL_MAX;
242
243    PhotInEMaxGamma     [i] = 0.;
244    PhotInEMaxElectron  [i] = 0.;
245    PhotInEMaxPositron  [i] = 0.;
246    PhotInEMaxNeutrons  [i] = 0.;
247    PhotInEMaxProtons   [i] = 0.;
248    PhotInEMaxDeuterons [i] = 0.;
249    PhotInEMaxTritons   [i] = 0.;
250    PhotInEMaxHe3s      [i] = 0.;
251    PhotInEMaxAlphas    [i] = 0.;
252    PhotInEMaxLambdas   [i] = 0.;
253    PhotInEMaxHeavyFrags[i] = 0.;
254    PhotInEMaxMesons    [i] = 0.;
255  }
256}
257
258
Note: See TracBrowser for help on using the repository browser.