ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Run.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file Run.hh
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 // This example is provided by the Geant4-DNA collaboration
27 // Any report or published results obtained using the Geant4-DNA software
28 // shall cite the following Geant4-DNA collaboration publication:
29 // Med. Phys. 37 (2010) 4692-4708
30 // and papers
31 // M. Batmunkh et al. J Radiat Res Appl Sci 8 (2015) 498-507
32 // O. Belov et al. Physica Medica 32 (2016) 1510-1520
33 // The Geant4-DNA web site is available at http://geant4-dna.org
34 //
35 // -------------------------------------------------------------------
36 // November 2016
37 // -------------------------------------------------------------------
38 //
41 
42 #ifndef Run_h
43 #define Run_h 1
44 
45 #include "G4Run.hh"
46 #include "G4VProcess.hh"
47 #include "globals.hh"
48 #include <map>
49 #include "G4Molecule.hh"
50 
51 class DetectorConstruction;
52 class G4ParticleDefinition;
53 //class NeuronLoadDataFile;
54 class G4Molecule;
55 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
56 
57 class Run : public G4Run
58 {
59  public:
60  Run(DetectorConstruction*);
61  ~Run();
62 
63  public:
64  void SetPrimary(G4ParticleDefinition* particle, G4double energy);
65  void AddPrimaryLET(G4double let);
66  void SetTrackLength(G4double tracklen);
67  void CountProcesses(const G4VProcess* process);
68  void ParticleCount(G4String, G4double);
69  void ParticleCountNeuron(G4String, G4double);
70  void MoleculeCount(G4String, G4double);
71  void MoleculeCountNeuron(G4Molecule* molecule);
72  void AddEdep (G4double edep);
73  void AddEflow (G4double eflow);
74  void ParticleFlux(G4String, G4double);
75 
76  virtual void Merge(const G4Run*);
77  void EndOfRun();
78 
79 // Edep in all volume
80  G4double GetEdepALL(){return fEdepAll;}
81  void SetEdepALL(G4double vall){ fEdepAll = vall;}
82  void AddEdepALL (G4double vall)
83  {
84  fEdepAll += vall;
85  fEdepAll_err += vall*vall;
86  }
87 // 0. Edep in homogeneous Medium
88  G4double GetEdepMedium(){return fEdepMedium;}
89  void SetEdepMedium(G4double vall){ fEdepMedium = vall;}
90  void AddEdepMedium (G4double vall)
91  {
92  fEdepMedium += vall;
93  fEdepMedium_err += vall*vall;
94  }
95 // 1. Edep in Bounding Slice Volume
96  G4double GetEdepSlice(){return fEdepSlice;}
97  void SetEdepSlice(G4double vall){ fEdepSlice = vall;}
98  void AddEdepSlice (G4double vall)
99  {
100  fEdepSlice += vall;
101  fEdepSlice_err += vall*vall;
102  }
103 // 2. Edep in Soma volume
104  G4double GetEdepSoma(){return fEdepSoma;}
105  void SetEdepSoma(G4double vall){ fEdepSoma = vall;}
106  void AddEdepSoma (G4double vall)
107  {
108  fEdepSoma += vall;
109  fEdepSoma_err += vall*vall;
110  }
111 
112  void AddSomaCompart(G4int i, G4double x){ fSoma3DEdep[i] +=x;}
113  G4double GetSomaCompart(G4int i){ return fSoma3DEdep[i];}
114  //G4double GetSomalength(G4int i){ return fSoma3DEdep[i];}
115  //G4ThreeVector GetVectSoma(G4int i) {return fSomaCompartments[i];}
116 
117 // 3. Edep in Dendrites volume
118  G4double GetEdepDend(){return fEdepDend;}
119  void SetEdepDend(G4double vall){ fEdepDend = vall;}
120  void AddEdepDend (G4double vall)
121  {
122  fEdepDend += vall;
123  fEdepDend_err += vall*vall;
124  }
125 
126  void AddDendCompart(G4int i, G4double x){ fDend3DEdep[i] +=x;}
127  G4double GetDendCompart(G4int i){ return fDend3DEdep[i];}
128 
129 // 4. Edep in Axon volume
130  G4double GetEdepAxon(){return fEdepAxon;}
131  void SetEdepAxon(G4double vall){ fEdepAxon = vall;}
132  void AddEdepAxon (G4double vall)
133  {
134  fEdepAxon += vall;
135  fEdepAxon_err += vall*vall;
136  }
137 
138  void AddAxonCompart(G4int i, G4double x){ fAxon3DEdep[i] +=x;}
139  G4double GetAxonCompart(G4int i){ return fAxon3DEdep[i];}
140 
141 // 5. Edep in whole Neuron volume
142  G4double GetEdepNeuron(){return fEdepNeuron;}
143  void SetEdepNeuron(G4double vall){ fEdepNeuron = vall;}
144  void AddEdepNeuron (G4double vall)
145  {
146  fEdepNeuron += vall;
147  fEdepNeuron_err += vall*vall;
148  }
149 
150  private:
151  struct ParticleData {
152  ParticleData()
153  : fCount(0), fEmean(0.), fEmin(0.), fEmax(0.) {}
154  ParticleData(G4int count, G4double ekin, G4double emin, G4double emax)
155  : fCount(count), fEmean(ekin), fEmin(emin), fEmax(emax) {}
156  G4int fCount;
157  G4double fEmean;
158  G4double fEmin;
159  G4double fEmax;
160  };
161 
162  //NeuronLoadDataFile * fNeuronLoadParamz;
163  DetectorConstruction* fDetector;
164  G4ParticleDefinition* fParticle;
165  G4double fEkin;
166  G4double fLET, fLET2;
167  G4double ftrackLength;
168  G4double fTrackLen, fTrackLen2;
169  G4int fNParticle;
170 
171  G4double * fSoma3DEdep;
172  G4double * fDend3DEdep ;
173  G4double * fAxon3DEdep ;
174  G4double * fNeuron3DEdep;
175  G4double fEdepAll, fEdepAll_err,fEdepMedium, fEdepMedium_err, fEdepSlice,
176  fEdepSlice_err,fEdepSoma,fEdepSoma_err, fEdepDend, fEdepDend_err,fEdepAxon,
177  fEdepAxon_err,fEdepNeuron, fEdepNeuron_err;
178 
179  G4int fmolNum, fmolNum2;
180  G4double fEnergyFlow, fEnergyFlow2;
181  std::map<G4String,G4int> fMoleculeNumber;
182  std::map<G4String,G4int> fProcCounter;
183  std::map<G4String,ParticleData> fParticleDataMap1;
184  std::map<G4String,ParticleData> fParticleDataMap2;
185 };
186 
187 #endif