29 #ifndef G4DNAChampionElasticModel_h
30 #define G4DNAChampionElasticModel_h 1
47 const G4String& nam =
"DNAChampionElasticModel");
84 typedef std::map<G4String,G4String,std::less<G4String> >
MapFile;
87 typedef std::map<G4String,G4DNACrossSectionDataSet*,std::less<G4String> >
MapData;
120 typedef std::map<double, std::vector<double> >
VecMap;
146 G4Exception (
"*** WARNING : G4DNAChampionElasticModel::SetKillBelowThreshold INACTIVE for now",
"",
JustWarning,
"") ;
G4double RandomizeCosTheta(G4double k)
std::map< G4String, G4String, std::less< G4String > > MapFile
G4DNAChampionElasticModel(const G4ParticleDefinition *p=0, const G4String &nam="DNAChampionElasticModel")
void SetKillBelowThreshold(G4double threshold)
virtual ~G4DNAChampionElasticModel()
std::vector< double > eTdummyVec
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)
std::map< G4String, G4DNACrossSectionDataSet *, std::less< G4String > > MapData
G4double LinLogInterpolate(G4double e1, G4double e2, G4double e, G4double xs1, G4double xs2)
virtual G4double CrossSectionPerVolume(const G4Material *material, const G4ParticleDefinition *p, G4double ekin, G4double emin, G4double emax)
std::map< double, std::map< double, double > > TriDimensionMap
G4double GetKillBelowThreshold()
G4double Theta(G4ParticleDefinition *aParticleDefinition, G4double k, G4double integrDiff)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
G4double QuadInterpolator(G4double e11, G4double e12, G4double e21, G4double e22, G4double x11, G4double x12, G4double x21, G4double x22, G4double t1, G4double t2, G4double t, G4double e)
G4DNAChampionElasticModel & operator=(const G4DNAChampionElasticModel &right)
G4ParticleChangeForGamma * fParticleChangeForGamma
TriDimensionMap eDiffCrossSectionData
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
std::map< double, std::vector< double > > VecMap
G4double LinLinInterpolate(G4double e1, G4double e2, G4double e, G4double xs1, G4double xs2)
G4double LogLogInterpolate(G4double e1, G4double e2, G4double e, G4double xs1, G4double xs2)
const std::vector< G4double > * fpMolWaterDensity