Geant4_10
G4DNAChampionElasticModel.hh
Go to the documentation of this file.
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: G4DNAChampionElasticModel.hh 70171 2013-05-24 13:34:18Z gcosmo $
27 //
28 
29 #ifndef G4DNAChampionElasticModel_h
30 #define G4DNAChampionElasticModel_h 1
31 
32 #include <map>
34 #include "G4VEmModel.hh"
35 #include "G4Electron.hh"
37 #include "G4LogLogInterpolation.hh"
38 #include "G4ProductionCutsTable.hh"
39 #include "G4NistManager.hh"
40 
42 {
43 
44 public:
45 
47  const G4String& nam = "DNAChampionElasticModel");
48 
50 
51  virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
52 
54  const G4ParticleDefinition* p,
55  G4double ekin,
56  G4double emin,
57  G4double emax);
58 
59  virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
60  const G4MaterialCutsCouple*,
61  const G4DynamicParticle*,
62  G4double tmin,
63  G4double maxEnergy);
64 
65  inline void SetKillBelowThreshold (G4double threshold);
66  G4double GetKillBelowThreshold () { return killBelowEnergy; }
67 
68 protected:
69 
71 
72 private:
73  // Water density table
74  const std::vector<G4double>* fpMolWaterDensity;
75 
76  G4double killBelowEnergy;
77  G4double lowEnergyLimit;
78  G4double highEnergyLimit;
79  G4bool isInitialised;
80  G4int verboseLevel;
81 
82  // Cross section
83 
84  typedef std::map<G4String,G4String,std::less<G4String> > MapFile;
85  MapFile tableFile;
86 
87  typedef std::map<G4String,G4DNACrossSectionDataSet*,std::less<G4String> > MapData;
88  MapData tableData;
89 
90  // Final state
91 
92  //G4double DifferentialCrossSection(G4ParticleDefinition * aParticleDefinition, G4double k, G4double theta);
93 
94  G4double Theta(G4ParticleDefinition * aParticleDefinition, G4double k, G4double integrDiff);
95 
96  G4double LinLinInterpolate(G4double e1, G4double e2, G4double e, G4double xs1, G4double xs2);
97 
98  G4double LinLogInterpolate(G4double e1, G4double e2, G4double e, G4double xs1, G4double xs2);
99 
100  G4double LogLogInterpolate(G4double e1, G4double e2, G4double e, G4double xs1, G4double xs2);
101 
102  G4double QuadInterpolator(G4double e11,
103  G4double e12,
104  G4double e21,
105  G4double e22,
106  G4double x11,
107  G4double x12,
108  G4double x21,
109  G4double x22,
110  G4double t1,
111  G4double t2,
112  G4double t,
113  G4double e);
114 
115  typedef std::map<double, std::map<double, double> > TriDimensionMap;
116 
117  TriDimensionMap eDiffCrossSectionData;
118  std::vector<double> eTdummyVec;
119 
120  typedef std::map<double, std::vector<double> > VecMap;
121  VecMap eVecm;
122 
123  G4double RandomizeCosTheta(G4double k);
124 
125  //
126 
129 
130 };
131 
132 
134 {
135 
136 // SI - commented on 19/06/2013
137 /*
138  killBelowEnergy = threshold;
139 
140  if (threshold < 1*eV)
141  G4Exception ("*** WARNING : the G4DNAChampionElasticModel class is not validated below 1 eV !","",JustWarning,"") ;
142 
143  if (threshold < 0.025*eV) threshold = 0.025*eV;
144 */
145 
146  G4Exception ("*** WARNING : G4DNAChampionElasticModel::SetKillBelowThreshold INACTIVE for now","",JustWarning,"") ;
147 
148 }
149 
150 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
151 
152 #endif
TTree * t1
Definition: plottest35.C:26
G4DNAChampionElasticModel(const G4ParticleDefinition *p=0, const G4String &nam="DNAChampionElasticModel")
void SetKillBelowThreshold(G4double threshold)
const char * p
Definition: xmltok.h:285
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)
int G4int
Definition: G4Types.hh:78
virtual G4double CrossSectionPerVolume(const G4Material *material, const G4ParticleDefinition *p, G4double ekin, G4double emin, G4double emax)
string material
Definition: eplot.py:19
bool G4bool
Definition: G4Types.hh:79
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
Definition: G4Exception.cc:41
G4ParticleChangeForGamma * fParticleChangeForGamma
TTree * t2
Definition: plottest35.C:36
double G4double
Definition: G4Types.hh:76
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)