57 G4cout <<
"Born excitation model is constructed " <<
G4endl;
68 std::map<G4String, G4DNACrossSectionDataSet*, std::less<G4String> >::iterator
pos;
84 G4cout <<
"Calling G4DNABornExcitationModel::Initialise()" <<
G4endl;
86 G4String fileElectron(
"dna/sigma_excitation_e_born");
87 G4String fileProton(
"dna/sigma_excitation_p_born");
131 if (particle==electronDef)
137 if (particle==protonDef)
145 G4cout <<
"Born excitation model is initialized " << G4endl
171 G4cout <<
"Calling CrossSectionPerVolume() of G4DNABornExcitationModel"
190 if(waterDensity!= 0.0)
195 std::map< G4String,G4double,std::less<G4String> >::iterator pos1;
199 lowLim = pos1->second;
202 std::map< G4String,G4double,std::less<G4String> >::iterator pos2;
206 highLim = pos2->second;
209 if (ekin >= lowLim && ekin < highLim)
211 std::map< G4String,G4DNACrossSectionDataSet*,std::less<G4String> >::iterator
pos;
224 G4Exception(
"G4DNABornExcitationModel::CrossSectionPerVolume",
"em0002",
231 G4cout <<
"__________________________________" <<
G4endl;
232 G4cout <<
"=== G4DNABornExcitationModel - XS INFO START" <<
G4endl;
233 G4cout <<
"=== Kinetic energy(eV)=" << ekin/
eV <<
" particle : " << particleName <<
G4endl;
234 G4cout <<
"=== Cross section per water molecule (cm^2)=" << sigma/
cm/
cm <<
G4endl;
235 G4cout <<
"=== Cross section per water molecule (cm^-1)=" << sigma*waterDensity/(1./
cm) << G4endl;
237 G4cout <<
"=== G4DNABornExcitationModel - XS INFO END" <<
G4endl;
242 return sigma*waterDensity;
256 G4cout <<
"Calling SampleSecondaries() of G4DNABornExcitationModel" <<
G4endl;
264 G4double newEnergy = k - excitationEnergy;
286 std::map<G4String, G4DNACrossSectionDataSet*, std::less<G4String> >::iterator
pos;
304 value += valuesBuffer[i];
315 if (valuesBuffer[i] > value)
317 delete[] valuesBuffer;
320 value -= valuesBuffer[i];
323 if (valuesBuffer)
delete[] valuesBuffer;
329 G4Exception(
"G4DNABornExcitationModel::RandomSelect",
"em0002",
static G4Electron * ElectronDefinition()
G4ParticleChangeForGamma * fParticleChangeForGamma
G4double LowEnergyLimit() const
virtual G4double FindValue(G4double x, G4int componentId=0) const =0
std::map< G4String, G4double, std::less< G4String > > lowEnergyLimit
G4double ExcitationEnergy(G4int level)
G4double GetKineticEnergy() const
G4double HighEnergyLimit() const
virtual const G4VEMDataSet * GetComponent(G4int componentId) const
static G4Material * GetMaterial(const G4String &name, G4bool warning=true)
virtual ~G4DNABornExcitationModel()
static G4Proton * ProtonDefinition()
virtual G4double CrossSectionPerVolume(const G4Material *material, const G4ParticleDefinition *p, G4double ekin, G4double emin, G4double emax)
virtual G4bool LoadData(const G4String &argFileName)
const std::vector< G4double > * fpMolWaterDensity
G4ParticleDefinition * GetDefinition() const
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
const G4String & GetParticleName() const
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
void SetHighEnergyLimit(G4double)
G4DNAWaterExcitationStructure waterStructure
G4DNABornExcitationModel(const G4ParticleDefinition *p=0, const G4String &nam="DNABornExcitationModel")
G4GLOB_DLL std::ostream G4cout
const std::vector< double > * GetNumMolPerVolTableFor(const G4Material *) const
G4int RandomSelect(G4double energy, const G4String &particle)
const G4ThreeVector & GetMomentumDirection() const
virtual G4double FindValue(G4double e, G4int componentId=0) const
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &=*(new G4DataVector()))
virtual size_t NumberOfComponents(void) const
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
static G4DNAChemistryManager * Instance()
static G4DNAMolecularMaterial * Instance()
void CreateWaterMolecule(ElectronicModification, G4int, const G4Track *)
Method used by DNA physics model to create a water molecule.
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)
std::map< G4String, G4double, std::less< G4String > > highEnergyLimit
const G4Track * GetCurrentTrack() const
void SetProposedKineticEnergy(G4double proposedKinEnergy)
void SetLowEnergyLimit(G4double)
static const G4double pos
G4ParticleChangeForGamma * GetParticleChangeForGamma()