64 G4cout <<
"Sanche Excitation model is constructed " <<
G4endl<<
"Energy range: "
86 G4cout <<
"Calling G4DNASancheExcitationModel::Initialise()" <<
G4endl;
92 G4cout <<
"G4DNASancheExcitationModel: low energy limit increased from " <<
99 G4cout <<
"G4DNASancheExcitationModel: high energy limit decreased from " <<
108 G4cout <<
"Sanche Excitation model is initialized " << G4endl
124 char *path = getenv(
"G4LEDATA");
125 std::ostringstream eFullFileName;
126 eFullFileName << path <<
"/dna/sigma_excitationvib_e_sanche.dat";
127 std::ifstream input(eFullFileName.str().c_str());
131 G4Exception(
"G4DNASancheExcitationModel::Initialise",
"em0003",
132 FatalException,
"Missing data file:/dna/sigma_excitationvib_e_sanche.dat");
159 G4cout <<
"Calling CrossSectionPerVolume() of G4DNASancheExcitationModel"
168 if(waterDensity!= 0.0)
182 G4cout <<
"__________________________________" <<
G4endl;
183 G4cout <<
"°°° G4DNASancheExcitationModel - XS INFO START" <<
G4endl;
185 G4cout <<
"°°° Cross section per water molecule (cm^2)=" << sigma/
cm/
cm <<
G4endl;
186 G4cout <<
"°°° Cross section per water molecule (cm^-1)=" << sigma*waterDensity/(1./
cm) << G4endl;
188 G4cout <<
"°°° G4DNASancheExcitationModel - XS INFO END" <<
G4endl;
194 return sigma*2*waterDensity;
210 G4cout <<
"Calling SampleSecondaries() of G4DNASancheExcitationModel"
216 G4double newEnergy = electronEnergy0 - excitationEnergy;
236 if (electronEnergy0 < highEnergyLimit && newEnergy>0.)
251 std::vector<double>::iterator t2 = std::upper_bound(
tdummyVec.begin(),
253 std::vector<double>::iterator t1 = t2 - 1;
257 sigma *= 1e-16 *
cm *
cm;
258 if (sigma == 0.) sigma = 1e-30;
266 G4double energies[9] = { 0.01, 0.024, 0.061, 0.092, 0.204, 0.417, 0.460,
268 return (energies[level] *
eV);
280 std::deque<double> values;
286 values.push_front(partial);
297 if (values[i] > value)
320 return totalCrossSection;
static G4Electron * ElectronDefinition()
G4DNASancheExcitationModel(const G4ParticleDefinition *p=0, const G4String &nam="DNASancheExcitationModel")
G4double LowEnergyLimit() const
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
G4double GetKineticEnergy() const
G4double HighEnergyLimit() const
static G4Material * GetMaterial(const G4String &name, G4bool warning=true)
G4double PartialCrossSection(G4double energy, G4int level)
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
const G4String & GetParticleName() const
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
void SetHighEnergyLimit(G4double)
G4double LinInterpolate(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)
G4GLOB_DLL std::ostream G4cout
const std::vector< G4double > * fpWaterDensity
const G4ThreeVector & GetMomentumDirection() const
G4ParticleChangeForGamma * fParticleChangeForGamma
virtual ~G4DNASancheExcitationModel()
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
G4double VibrationEnergy(G4int level)
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)
static G4DNAMolecularMaterial * Instance()
G4int RandomSelect(G4double energy)
void SetProposedKineticEnergy(G4double proposedKinEnergy)
std::vector< double > tdummyVec
void SetLowEnergyLimit(G4double)
G4ParticleChangeForGamma * GetParticleChangeForGamma()