source: trunk/source/processes/electromagnetic/lowenergy/test/fluoTest/src/XrayFluoAnalysisManager.cc @ 1350

Last change on this file since 1350 was 1350, checked in by garnier, 13 years ago

update to last version 4.9.4
File size: 8.2 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// $Id: XrayFluoAnalysisManager.cc
28// GEANT4 tag $Name: xray_fluo-V03-02-00
29//
30// Author: Elena Guardincerri (Elena.Guardincerri@ge.infn.it)
31//
32// History:
33// -----------
34// 28 Nov 2001 Elena Guardincerri     Created
35//
36// -------------------------------------------------------------------
37#ifdef  G4ANALYSIS_USE
38
39#include "XrayFluoAnalysisManager.hh"
40#include "XrayFluoAnalysisMessenger.hh"
41#include "G4Step.hh"
42
43XrayFluoAnalysisManager* XrayFluoAnalysisManager::instance = 0;
44
45//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
46
47XrayFluoAnalysisManager::XrayFluoAnalysisManager()
48{
49  analysisMessenger = new XrayFluoAnalysisMessenger(this);
50  histoManager = createIHistoManager(); 
51  factory = Lizard::createNTupleFactory();
52}
53
54//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
55
56XrayFluoAnalysisManager::~XrayFluoAnalysisManager() 
57{
58
59  delete analysisMessenger; 
60  analysisMessenger = 0;
61  delete  histoManager;
62  histoManager = 0;
63  delete factory;
64  factory = 0;
65}
66
67//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
68
69XrayFluoAnalysisManager* XrayFluoAnalysisManager::getInstance()
70{
71  if (instance == 0) instance = new XrayFluoAnalysisManager;
72  return instance;
73}
74//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
75
76
77void XrayFluoAnalysisManager::book()
78
79{
80  histoManager->selectStore("XrayFluo.his");
81
82  // Book histograms
83  histoManager->create1D("1","Energy Deposit", 100,0.,10.);
84 histoManager->create1D("2","Gamma born in the sample", 100,0.,10.);
85 histoManager->create1D("3","Electrons  born in the sample", 100,0.,10.);
86 histoManager->create1D("4","Gammas leaving the sample", 100,0.,10.);
87 histoManager->create1D("5","Electrons leaving the sample ", 100,0.,10.);
88 histoManager->create1D("6","Gammas reaching the detector", 100,0.,10.);
89 histoManager->create1D("7","Spectrum of the incident particles", 100,0.,10.);
90 histoManager->create1D("8","Protons reaching the detector", 100,0.,10.);
91 histoManager->create1D("9","Protons leaving the sample", 100,0.,10.);
92
93 // Book ntuples
94  ntuple = factory->createC("XrayFluo.his::1");
95
96 //  Add and bind the attributes to the ntuple
97  if ( !( ntuple->addAndBind( "energy", eDep) &&
98          ntuple->addAndBind( "counts"     , counts   ) ) )
99
100    {
101      delete ntuple;
102      G4Exception("XrayFluoAnalysisManager::book - Could not addAndBind ntuple");
103    } 
104}
105
106//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
107
108void XrayFluoAnalysisManager::finish()
109{
110
111  histoManager->store("1");
112  histoManager->store("3");
113  histoManager->store("2");
114  histoManager->store("4");
115  histoManager->store("5");
116  histoManager->store("6");
117  histoManager->store("7");
118  histoManager->store("8");
119  histoManager->store("9");
120
121
122  delete ntuple;
123  ntuple = 0;
124  G4cout << "Deleted ntuple" << G4endl;
125}
126
127//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
128
129void XrayFluoAnalysisManager::analyseStepping(const G4Step* aStep)
130{
131  G4double gammaAtTheDetPre=0;
132  G4double protonsAtTheDetPre=0;
133  G4double gammaLeavingSample=0;
134  G4double eleLeavingSample=0;
135  G4double protonsLeavSam=0;
136  G4double gammaBornInSample=0;
137  G4double eleBornInSample=0;
138  if(aStep->GetPreStepPoint()->GetPhysicalVolume()->GetName()=="Sample"){
139   
140    if(aStep->GetTrack()->GetNextVolume()->GetName() == "World" ) 
141      { 
142        if ((aStep->GetTrack()->GetDynamicParticle()
143             ->GetDefinition()-> GetParticleName()) == "gamma" )
144          {
145            gammaLeavingSample = (aStep->GetPreStepPoint()->GetKineticEnergy());
146            IHistogram1D* h1 = histoManager->retrieveHisto1D("4");
147            h1->fill(gammaLeavingSample/keV);
148          }
149      }
150  }
151  if(aStep->GetPreStepPoint()->GetPhysicalVolume()->GetName()=="Sample"){
152   
153    if(aStep->GetTrack()->GetNextVolume()->GetName() == "World" ) 
154      { 
155        if ((aStep->GetTrack()->GetDynamicParticle()
156             ->GetDefinition()-> GetParticleName()) == "e-" ) 
157          {
158            eleLeavingSample = (aStep->GetPreStepPoint()->GetKineticEnergy());
159           
160            IHistogram1D* h2 = histoManager->retrieveHisto1D("5");
161            h2->fill(eleLeavingSample/keV);
162          }
163        else if ((aStep->GetTrack()->GetDynamicParticle()
164                  ->GetDefinition()-> GetParticleName()) == "proton" )
165          {
166            protonsLeavSam = (aStep->GetPreStepPoint()->GetKineticEnergy());
167             IHistogram1D* h3 = histoManager->retrieveHisto1D("9");
168             h3->fill(protonsLeavSam/keV);
169          }
170       
171      }
172  }
173 
174 
175  if((aStep->GetTrack()->GetDynamicParticle()
176      ->GetDefinition()-> GetParticleName()) == "gamma" )
177   
178    {if(1== (aStep->GetTrack()->GetCurrentStepNumber()))
179     
180      {if(0 != aStep->GetTrack()->GetParentID())
181       
182        {if(aStep->GetTrack()->GetVolume()->GetName() == "Sample")
183          {
184            gammaBornInSample = (aStep->GetPreStepPoint()->GetKineticEnergy());
185            IHistogram1D* h4 = histoManager->retrieveHisto1D("2");
186            h4->fill(gammaBornInSample);
187           
188          }
189        }
190      }
191    }
192  if((aStep->GetTrack()->GetDynamicParticle()
193      ->GetDefinition()-> GetParticleName()) == "e-" )
194   
195    {if(1== (aStep->GetTrack()->GetCurrentStepNumber()))
196     
197      {if(0 != aStep->GetTrack()->GetParentID())
198       
199        {if(aStep->GetTrack()->GetVolume()->GetName() == "Sample")
200          {
201            eleBornInSample = (aStep->GetPreStepPoint()->GetKineticEnergy());
202            IHistogram1D* h5 = histoManager->retrieveHisto1D("3");
203            h5->fill(eleBornInSample);
204          }
205        }
206      }
207    }
208 
209  if(aStep->GetTrack()->GetNextVolume()){
210   
211    if(aStep->GetTrack()->GetVolume()->GetName() == "World"){
212     
213      if(aStep->GetTrack()->GetNextVolume()->GetName() == "HPGeDetector")
214       
215        { 
216          if ((aStep->GetTrack()->GetDynamicParticle()
217               ->GetDefinition()-> GetParticleName()) == "gamma" ) 
218            {
219              gammaAtTheDetPre = (aStep->GetPreStepPoint()->GetKineticEnergy());
220              IHistogram1D* h6 = histoManager->retrieveHisto1D("6");
221              h6->fill( gammaAtTheDetPre);
222            }
223          else if ((aStep->GetTrack()->GetDynamicParticle()
224                    ->GetDefinition()-> GetParticleName()) == "proton" ) 
225            {
226              protonsAtTheDetPre = (aStep->GetPreStepPoint()->GetKineticEnergy());
227              IHistogram1D* h7 = histoManager->retrieveHisto1D("8");
228              h7->fill( protonsAtTheDetPre);
229            }
230        }
231    }
232  }
233}
234//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
235
236void XrayFluoAnalysisManager::analyseEnergyDep(G4double energyDep)
237{
238  IHistogram1D* h8 = histoManager->retrieveHisto1D("1");
239  h8->fill(energyDep/keV);
240  counts = 1.;
241  ntuple->addRow();
242
243}
244
245//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
246
247void XrayFluoAnalysisManager::analysePrimaryGenerator(G4double energy)
248{
249 IHistogram1D* h9 = histoManager->retrieveHisto1D("7");
250              h9->fill(energy/keV);
251
252}
253#endif
254
255
256
257
258
259
260
Note: See TracBrowser for help on using the repository browser.