100 latDisplasment(true),
158 if(index >= n) {
for(
G4int i=n; i<=index; ++i) {
mscModels.push_back(0); } }
187 if(masterProc && masterProc !=
this) { master =
false; }
195 if(pname !=
"deuteron" && pname !=
"triton" &&
196 pname !=
"alpha+" && pname !=
"helium" &&
197 pname !=
"alpha" && pname !=
"He3" &&
198 pname !=
"hydrogen") {
206 size_t n = v->
size();
207 for(
size_t j=0; j<
n; ++j) {
208 if((*v)[j] ==
this) {
222 G4cout <<
"### G4VMultipleScattering::PrepearPhysicsTable() for "
226 <<
" isIon= " <<
isIon
275 G4cout <<
"### G4VMultipleScattering::BuildPhysicsTable() for "
277 <<
" and particle " << num
284 if(masterProcess && masterProcess !=
this) { master =
false; }
322 num ==
"e+" || num ==
"mu+" ||
323 num ==
"mu-" || num ==
"proton"||
324 num ==
"pi+" || num ==
"pi-" ||
325 num ==
"kaon+" || num ==
"kaon-" ||
326 num ==
"alpha" || num ==
"anti_proton" ||
327 num ==
"GenericIon")))
338 G4cout <<
"### G4VMultipleScattering::BuildPhysicsTable() done for "
340 <<
" and particle " << num
501 G4double postSafety= preSafety - maxDisp;
502 G4bool safetyRecomputed =
false;
503 if(postSafety < maxDisp) {
504 safetyRecomputed =
true;
522 if(r2 > postSafety*postSafety) {
524 if(!safetyRecomputed) {
528 if(dispR > postSafety) {
529 fac = 0.99*postSafety/dispR;
605 static const G4String ss[4] = {
"1",
"2",
"3",
"4"};
606 for(
G4int i=0; i<nmod; ++i) {
618 G4cout <<
"Physics table are stored for "
621 <<
" with a name <" << name <<
"> " <<
G4endl;
624 G4cout <<
"Fail to store Physics Table for "
627 <<
" in the directory <" << directory
G4double condition(const G4ErrorSymMatrix &m)
G4bool RetrievePhysicsTable(const G4ParticleDefinition *, const G4String &directory, G4bool ascii)
const G4VProcess * GetMasterProcess() const
G4VEmModel * SelectModel(G4double kinEnergy, size_t idx)
virtual void Initialize(const G4Track &)
G4SafetyHelper * GetSafetyHelper() const
G4ParticleDefinition * FindParticle(G4int PDGEncoding)
void ReLocateWithinVolume(const G4ThreeVector &pGlobalPoint)
virtual void StartTracking(G4Track *)
static G4LossTableManager * Instance()
G4EmModelManager * modelManager
static const G4double fac
CLHEP::Hep3Vector G4ThreeVector
G4double GetStepLength() const
G4double HighEnergyLimit() const
void SetLateralDisplasmentFlag(G4bool val)
G4ThreeVector fNewPosition
G4double GeomFactor() const
G4VEmModel * GetModel(G4int, G4bool ver=false)
void DeRegister(G4VEnergyLossProcess *p)
const G4MaterialCutsCouple * GetMaterialCutsCouple() const
G4double GetMeanFreePath(const G4Track &track, G4double, G4ForceCondition *condition)
void SetLateralDisplasmentFlag(G4bool val)
G4double PostStepGetPhysicalInteractionLength(const G4Track &, G4double previousStepSize, G4ForceCondition *condition)
void SetEmModel(G4VMscModel *, G4int index=1)
void AddEmModel(G4int, G4VEmModel *, G4VEmFluctuationModel *, const G4Region *)
G4PhysicsTable * GetCrossSectionTable()
virtual ~G4VMultipleScattering()
void SetStepLimitType(G4MscStepLimitType)
G4ProcessManager * GetProcessManager() const
G4VEnergyLossProcess * fIonisation
const G4String & GetPhysicsTableFileName(const G4ParticleDefinition *, const G4String &directory, const G4String &tableName, G4bool ascii=false)
const G4String & GetParticleName() const
void SetHighEnergyLimit(G4double)
G4VMultipleScattering(const G4String &name="msc", G4ProcessType type=fElectromagnetic)
G4StepPoint * GetPreStepPoint() const
const G4ParticleDefinition * currParticle
G4MscStepLimitType StepLimitType() const
G4ParticleChangeForMSC fParticleChange
virtual G4double ComputeTrueStepLength(G4double geomPathLength)
const G4ThreeVector & GetMomentumDirection() const
G4double RangeFactor() const
G4double GetKineticEnergy() const
G4GLOB_DLL std::ostream G4cout
virtual G4double ComputeTruePathLengthLimit(const G4Track &track, G4double &stepLimit)
G4VParticleChange * PostStepDoIt(const G4Track &, const G4Step &)
void SetParticleChange(G4VParticleChange *, G4VEmFluctuationModel *f=0)
const G4ThreeVector & GetPosition() const
G4double GetRange(const G4ParticleDefinition *part, G4double kineticEnergy, const G4MaterialCutsCouple *couple)
G4bool StorePhysicsTable(const G4ParticleDefinition *, const G4String &directory, G4bool ascii=false)
void SetRangeFactor(G4double)
void ProposeTrueStepLength(G4double truePathLength)
const G4ParticleDefinition * GetParticleDefinition() const
void SetProcessSubType(G4int)
const G4String & GetParticleType() const
void SetCrossSectionTable(G4PhysicsTable *, G4bool isLocal)
void Register(G4VEnergyLossProcess *p)
G4LossTableManager * emManager
const G4String & GetProcessName() const
G4bool IsActive(G4double kinEnergy)
void DumpModelList(G4int verb)
G4bool LateralDisplasmentFlag() const
const G4DataVector * Initialise(const G4ParticleDefinition *part, const G4ParticleDefinition *secPart, G4double minSubRange, G4int verb)
G4double ComputeSafety(const G4ThreeVector &pGlobalPoint, G4double maxRadius=DBL_MAX)
static G4TransportationManager * GetTransportationManager()
void SetMasterThread(G4bool val)
G4VEnergyLossProcess * GetEnergyLossProcess(const G4ParticleDefinition *)
void ProposeMomentumDirection(const G4ThreeVector &Pfinal)
void SetIonisation(G4VEnergyLossProcess *, const G4ParticleDefinition *part)
static G4GenericIon * GenericIon()
G4int NumberOfModels() const
void ProposePosition(const G4ThreeVector &finalPosition)
void BuildPhysicsTable(const G4ParticleDefinition *aParticle)
G4double GetPDGMass() const
static G4ParticleTable * GetParticleTable()
G4VParticleChange * AlongStepDoIt(const G4Track &, const G4Step &)
G4double ContinuousStepLimit(const G4Track &track, G4double previousStepSize, G4double currentMinimalStep, G4double ¤tSafety)
const G4ParticleDefinition * firstParticle
G4StepPoint * GetPostStepPoint() const
void AddEmModel(G4int order, G4VEmModel *, const G4Region *region=0)
G4VParticleChange * pParticleChange
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4double GetSafety() const
void PreparePhysicsTable(const G4ParticleDefinition *aParticle, G4VEnergyLossProcess *p, G4bool theMaster)
virtual void InitialiseProcess(const G4ParticleDefinition *)=0
void SetGeomFactor(G4double)
static G4Electron * Electron()
G4bool StorePhysicsTable(const G4String &filename, G4bool ascii=false)
virtual void PrintInfo()=0
void PrintInfoDefinition()
void BuildPhysicsTable(const G4ParticleDefinition &)
G4ProcessVector * GetAlongStepProcessVector(G4ProcessVectorTypeIndex typ=typeGPIL) const
void SetRangeFactor(G4double val)
G4double AlongStepGetPhysicalInteractionLength(const G4Track &, G4double previousStepSize, G4double currentMinimalStep, G4double ¤tSafety, G4GPILSelection *selection)
void StartTracking(G4Track *)
G4double MaxKinEnergy() const
G4VMscModel * EmModel(G4int index=1) const
G4int GetProcessSubType() const
std::vector< G4VMscModel * > mscModels
void SetPolarAngleLimit(G4double)
void SetVerboseLevel(G4int value)
G4double GetContinuousStepLimit(const G4Track &track, G4double previousStepSize, G4double currentMinimalStep, G4double ¤tSafety)
G4VMscModel * currentModel
void SetStepLimitType(G4MscStepLimitType val)
G4VEmModel * GetModelByIndex(G4int idx=0, G4bool ver=false) const
void PreparePhysicsTable(const G4ParticleDefinition &)
G4SafetyHelper * safetyHelper
virtual G4ThreeVector & SampleScattering(const G4ThreeVector &, G4double safety)
void SetIonisation(G4VEnergyLossProcess *)