| 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: F02RunAction.hh,v 1.4 2006/06/29 18:27:30 gunter Exp $
|
|---|
| 28 | // GEANT4 tag $Name: HEAD $
|
|---|
| 29 | //
|
|---|
| 30 | //
|
|---|
| 31 |
|
|---|
| 32 | //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
|
|---|
| 33 | //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
|
|---|
| 34 |
|
|---|
| 35 | #ifndef F02RunAction_h
|
|---|
| 36 | #define F02RunAction_h 1
|
|---|
| 37 |
|
|---|
| 38 | #include "G4UserRunAction.hh"
|
|---|
| 39 | #include "globals.hh"
|
|---|
| 40 | #include <iostream>
|
|---|
| 41 |
|
|---|
| 42 | //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
|
|---|
| 43 |
|
|---|
| 44 | class F02RunMessenger;
|
|---|
| 45 | class G4Run;
|
|---|
| 46 |
|
|---|
| 47 | #ifndef G4NOHIST
|
|---|
| 48 | class HepTupleManager;
|
|---|
| 49 | class HepHistogram;
|
|---|
| 50 | #endif
|
|---|
| 51 |
|
|---|
| 52 | //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
|
|---|
| 53 |
|
|---|
| 54 | class F02RunAction : public G4UserRunAction
|
|---|
| 55 | {
|
|---|
| 56 | public:
|
|---|
| 57 | F02RunAction();
|
|---|
| 58 | ~F02RunAction();
|
|---|
| 59 |
|
|---|
| 60 | public:
|
|---|
| 61 | void BeginOfRunAction(const G4Run*);
|
|---|
| 62 | void EndOfRunAction(const G4Run*);
|
|---|
| 63 |
|
|---|
| 64 | void CountEvent();
|
|---|
| 65 | void CountParticles(G4double,G4double);
|
|---|
| 66 | void AddEP(G4double,G4double);
|
|---|
| 67 | void AddEdeps(G4double Eabs);
|
|---|
| 68 | void AddTrackLength(G4double tlabs);
|
|---|
| 69 | void AddnStepsCharged(G4double ns);
|
|---|
| 70 | void AddnStepsNeutral(G4double ns);
|
|---|
| 71 | void AddnTimeGlob(G4double timeGlob);
|
|---|
| 72 |
|
|---|
| 73 | void AddTrRef(G4double tr,G4double ref);
|
|---|
| 74 |
|
|---|
| 75 | void FillEn(G4double En);
|
|---|
| 76 | void FillTh(G4double Th);
|
|---|
| 77 | void FillThBack(G4double Th);
|
|---|
| 78 | void FillR(G4double R);
|
|---|
| 79 | void FillTt(G4double Tt);
|
|---|
| 80 | void FillTb(G4double Tt);
|
|---|
| 81 | void FillTsec(G4double T);
|
|---|
| 82 | void FillGammaSpectrum(G4double E);
|
|---|
| 83 | void FillNbOfSteps(G4double nstep);
|
|---|
| 84 | void Fillvertexz(G4double z);
|
|---|
| 85 | void FillTglob(G4double timeGlob);
|
|---|
| 86 |
|
|---|
| 87 | void SethistName(G4String name) ;
|
|---|
| 88 |
|
|---|
| 89 | void SetnbinStep(G4int nbin);
|
|---|
| 90 | void SetSteplow(G4double Slow);
|
|---|
| 91 | void SetStephigh(G4double Shigh);
|
|---|
| 92 |
|
|---|
| 93 | void SetnbinEn(G4int nbin);
|
|---|
| 94 | void SetEnlow(G4double Elow);
|
|---|
| 95 | void SetEnhigh(G4double Enhigh);
|
|---|
| 96 |
|
|---|
| 97 | void SetnbinTt(G4int nbin);
|
|---|
| 98 | void SetTtlow(G4double Ttlow);
|
|---|
| 99 | void SetTthigh(G4double Tthigh);
|
|---|
| 100 |
|
|---|
| 101 | void SetnbinTb(G4int nbin);
|
|---|
| 102 | void SetTblow(G4double Tblow);
|
|---|
| 103 | void SetTbhigh(G4double Tbhigh);
|
|---|
| 104 |
|
|---|
| 105 | void SetnbinTsec(G4int nbin);
|
|---|
| 106 | void SetTseclow(G4double Tlow);
|
|---|
| 107 | void SetTsechigh(G4double Thigh);
|
|---|
| 108 |
|
|---|
| 109 | void SetnbinTglob(G4int nbin);
|
|---|
| 110 | void SetTgloblow(G4double Tgloblow);
|
|---|
| 111 | void SetTglobhigh(G4double Tglobhigh);
|
|---|
| 112 |
|
|---|
| 113 | void SetnbinTh(G4int nbin);
|
|---|
| 114 | void SetThlow(G4double Thlow);
|
|---|
| 115 | void SetThhigh(G4double Thhigh);
|
|---|
| 116 |
|
|---|
| 117 | void SetnbinThBack(G4int nbin);
|
|---|
| 118 | void SetThlowBack(G4double Thlow);
|
|---|
| 119 | void SetThhighBack(G4double Thhigh);
|
|---|
| 120 |
|
|---|
| 121 | void SetnbinR(G4int nbin);
|
|---|
| 122 | void SetRlow(G4double Rlow);
|
|---|
| 123 | void SetRhigh(G4double Rhigh);
|
|---|
| 124 |
|
|---|
| 125 | void SetnbinGamma(G4int nbin);
|
|---|
| 126 | void SetElowGamma(G4double Elow);
|
|---|
| 127 | void SetEhighGamma(G4double Ehigh);
|
|---|
| 128 |
|
|---|
| 129 | void Setnbinzvertex(G4int nbin);
|
|---|
| 130 | void Setzlow(G4double z);
|
|---|
| 131 | void Setzhigh(G4double z);
|
|---|
| 132 |
|
|---|
| 133 | void SetRndmFreq(G4int val) {saveRndm = val;}
|
|---|
| 134 | G4int GetRndmFreq() {return saveRndm;}
|
|---|
| 135 |
|
|---|
| 136 | private:
|
|---|
| 137 |
|
|---|
| 138 | void bookHisto();
|
|---|
| 139 |
|
|---|
| 140 | private:
|
|---|
| 141 |
|
|---|
| 142 | G4String histName ;
|
|---|
| 143 |
|
|---|
| 144 | #ifndef G4NOHIST
|
|---|
| 145 | HepTupleManager* hbookManager;
|
|---|
| 146 | HepHistogram *histo1, *histo2, *histo3, *histo4, *histo5 ;
|
|---|
| 147 | HepHistogram *histo6, *histo7, *histo8, *histo9, *histo10, *histo11;
|
|---|
| 148 | #endif
|
|---|
| 149 |
|
|---|
| 150 | G4double EnergySumAbs,EnergySquareSumAbs;
|
|---|
| 151 | G4double tlSumAbs,tlsquareSumAbs;
|
|---|
| 152 | G4double nStepSumCharged,nStepSum2Charged ;
|
|---|
| 153 | G4double nStepSumNeutral,nStepSum2Neutral ;
|
|---|
| 154 | G4double TotNbofEvents;
|
|---|
| 155 | G4double SumCharged,Sum2Charged,SumNeutral,Sum2Neutral;
|
|---|
| 156 | G4double Selectron,Spositron;
|
|---|
| 157 |
|
|---|
| 158 | G4double Transmitted,Reflected ;
|
|---|
| 159 |
|
|---|
| 160 | G4double entryStep,underStep,overStep,distStep[200];
|
|---|
| 161 | G4double Steplow,Stephigh,dStep;
|
|---|
| 162 | G4int nbinStep;
|
|---|
| 163 |
|
|---|
| 164 | G4double entryEn,underEn,overEn,distEn[200];
|
|---|
| 165 | G4double Enlow,Enhigh,dEn;
|
|---|
| 166 | G4int nbinEn;
|
|---|
| 167 |
|
|---|
| 168 | G4double entryTt,underTt,overTt,distTt[200];
|
|---|
| 169 | G4double Ttlow,Tthigh,dTt;
|
|---|
| 170 | G4int nbinTt;
|
|---|
| 171 | G4double Ttmean,Tt2mean;
|
|---|
| 172 |
|
|---|
| 173 | G4double entryTb,underTb,overTb,distTb[200];
|
|---|
| 174 | G4double Tblow,Tbhigh,dTb;
|
|---|
| 175 | G4int nbinTb;
|
|---|
| 176 | G4double Tbmean,Tb2mean;
|
|---|
| 177 |
|
|---|
| 178 | G4double entryTsec,underTsec,overTsec,distTsec[200];
|
|---|
| 179 | G4double Tseclow,Tsechigh,dTsec;
|
|---|
| 180 | G4int nbinTsec;
|
|---|
| 181 |
|
|---|
| 182 | G4double entryTglob,underTglob,overTglob,distTglob[200];
|
|---|
| 183 | G4double Tgloblow,Tglobhigh,dTglob;
|
|---|
| 184 | G4int nbinTglob;
|
|---|
| 185 |
|
|---|
| 186 | G4double entryTh,underTh,overTh,distTh[200];
|
|---|
| 187 | G4double Thlow,Thhigh,dTh;
|
|---|
| 188 | G4int nbinTh;
|
|---|
| 189 |
|
|---|
| 190 | G4double entryThback,underThback,overThback,distThback[200];
|
|---|
| 191 | G4double Thlowback,Thhighback,dThback;
|
|---|
| 192 | G4int nbinThback;
|
|---|
| 193 |
|
|---|
| 194 | G4double entryR ,underR ,overR ,distR[200];
|
|---|
| 195 | G4double Rlow,Rhigh,dR;
|
|---|
| 196 | G4int nbinR;
|
|---|
| 197 |
|
|---|
| 198 | G4double Rmean,R2mean;
|
|---|
| 199 | G4double entryGamma,underGamma,overGamma,distGamma[200];
|
|---|
| 200 | G4double ElowGamma,EhighGamma,dEGamma;
|
|---|
| 201 | G4int nbinGamma ;
|
|---|
| 202 |
|
|---|
| 203 | G4double entryvertexz,undervertexz,oververtexz,distvertexz[200];
|
|---|
| 204 | G4double zlow,zhigh,dz;
|
|---|
| 205 | G4int nbinvertexz;
|
|---|
| 206 |
|
|---|
| 207 | F02RunMessenger* runMessenger;
|
|---|
| 208 | G4int saveRndm;
|
|---|
| 209 | };
|
|---|
| 210 |
|
|---|
| 211 | #endif
|
|---|
| 212 |
|
|---|