#include <G4VMultipleScattering.hh>

Inheritance diagram for G4VMultipleScattering:

Public Member Functions
	G4VMultipleScattering (const G4String &name="msc", G4ProcessType type=fElectromagnetic)

virtual	~G4VMultipleScattering ()

virtual G4bool	IsApplicable (const G4ParticleDefinition &p)=0

virtual void	PrintInfo ()=0

void	PreparePhysicsTable (const G4ParticleDefinition &)

void	BuildPhysicsTable (const G4ParticleDefinition &)

void	PrintInfoDefinition ()

G4bool	StorePhysicsTable (const G4ParticleDefinition *, const G4String &directory, G4bool ascii=false)

G4bool	RetrievePhysicsTable (const G4ParticleDefinition *, const G4String &directory, G4bool ascii)

void	StartTracking (G4Track *)

G4double	AlongStepGetPhysicalInteractionLength (const G4Track &, G4double previousStepSize, G4double currentMinimalStep, G4double &currentSafety, G4GPILSelection *selection)

G4double	PostStepGetPhysicalInteractionLength (const G4Track &, G4double previousStepSize, G4ForceCondition *condition)

G4VParticleChange *	AlongStepDoIt (const G4Track &, const G4Step &)

G4VParticleChange *	PostStepDoIt (const G4Track &, const G4Step &)

G4double	ContinuousStepLimit (const G4Track &track, G4double previousStepSize, G4double currentMinimalStep, G4double &currentSafety)

G4VEmModel *	SelectModel (G4double kinEnergy, size_t idx)

void	AddEmModel (G4int order, G4VEmModel , const G4Region region=0)

void	SetModel (G4VMscModel *, G4int index=1)

G4VMscModel *	Model (G4int index=1)

void	SetEmModel (G4VMscModel *, G4int index=1)

G4VMscModel *	EmModel (G4int index=1)

G4VEmModel *	GetModelByIndex (G4int idx=0, G4bool ver=false) const

void	SetIonisation (G4VEnergyLossProcess *)

G4bool	LateralDisplasmentFlag () const

void	SetLateralDisplasmentFlag (G4bool val)

G4double	Skin () const

void	SetSkin (G4double val)

G4double	RangeFactor () const

void	SetRangeFactor (G4double val)

G4double	GeomFactor () const

void	SetGeomFactor (G4double val)

G4double	PolarAngleLimit () const

void	SetPolarAngleLimit (G4double val)

G4MscStepLimitType	StepLimitType () const

void	SetStepLimitType (G4MscStepLimitType val)

const G4ParticleDefinition *	FirstParticle () const

Public Member Functions inherited from G4VContinuousDiscreteProcess
	G4VContinuousDiscreteProcess (const G4String &, G4ProcessType aType=fNotDefined)

	G4VContinuousDiscreteProcess (G4VContinuousDiscreteProcess &)

virtual	~G4VContinuousDiscreteProcess ()

virtual G4double	PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)

virtual G4VParticleChange *	PostStepDoIt (const G4Track &, const G4Step &)

virtual G4double	AlongStepGetPhysicalInteractionLength (const G4Track &, G4double previousStepSize, G4double currentMinimumStep, G4double &currentSafety, G4GPILSelection *selection)

virtual G4VParticleChange *	AlongStepDoIt (const G4Track &, const G4Step &)

virtual G4double	AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)

virtual G4VParticleChange *	AtRestDoIt (const G4Track &, const G4Step &)

Public Member Functions inherited from G4VProcess
	G4VProcess (const G4String &aName="NoName", G4ProcessType aType=fNotDefined)

	G4VProcess (const G4VProcess &right)

virtual	~G4VProcess ()

G4int	operator== (const G4VProcess &right) const

G4int	operator!= (const G4VProcess &right) const

virtual G4VParticleChange *	PostStepDoIt (const G4Track &track, const G4Step &stepData)=0

virtual G4VParticleChange *	AlongStepDoIt (const G4Track &track, const G4Step &stepData)=0

virtual G4VParticleChange *	AtRestDoIt (const G4Track &track, const G4Step &stepData)=0

virtual G4double	AlongStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)=0

virtual G4double	AtRestGetPhysicalInteractionLength (const G4Track &track, G4ForceCondition *condition)=0

virtual G4double	PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)=0

G4double	GetCurrentInteractionLength () const

void	SetPILfactor (G4double value)

G4double	GetPILfactor () const

G4double	AlongStepGPIL (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)

G4double	AtRestGPIL (const G4Track &track, G4ForceCondition *condition)

G4double	PostStepGPIL (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)

virtual G4bool	IsApplicable (const G4ParticleDefinition &)

virtual void	BuildPhysicsTable (const G4ParticleDefinition &)

virtual void	PreparePhysicsTable (const G4ParticleDefinition &)

virtual G4bool	StorePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)

virtual G4bool	RetrievePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)

const G4String &	GetPhysicsTableFileName (const G4ParticleDefinition *, const G4String &directory, const G4String &tableName, G4bool ascii=false)

const G4String &	GetProcessName () const

G4ProcessType	GetProcessType () const

void	SetProcessType (G4ProcessType)

G4int	GetProcessSubType () const

void	SetProcessSubType (G4int)

virtual void	StartTracking (G4Track *)

virtual void	EndTracking ()

virtual void	SetProcessManager (const G4ProcessManager *)

virtual const G4ProcessManager *	GetProcessManager ()

virtual void	ResetNumberOfInteractionLengthLeft ()

G4double	GetNumberOfInteractionLengthLeft () const

G4double	GetTotalNumberOfInteractionLengthTraversed () const

G4bool	isAtRestDoItIsEnabled () const

G4bool	isAlongStepDoItIsEnabled () const

G4bool	isPostStepDoItIsEnabled () const

virtual void	DumpInfo () const

void	SetVerboseLevel (G4int value)

G4int	GetVerboseLevel () const

Protected Member Functions
virtual void	InitialiseProcess (const G4ParticleDefinition *)=0

G4double	GetMeanFreePath (const G4Track &track, G4double, G4ForceCondition *condition)

G4double	GetContinuousStepLimit (const G4Track &track, G4double previousStepSize, G4double currentMinimalStep, G4double &currentSafety)

Protected Member Functions inherited from G4VContinuousDiscreteProcess
virtual G4double	GetMeanFreePath (const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *condition)=0

virtual G4double	GetContinuousStepLimit (const G4Track &aTrack, G4double previousStepSize, G4double currentMinimumStep, G4double &currentSafety)=0

void	SetGPILSelection (G4GPILSelection selection)

G4GPILSelection	GetGPILSelection () const

Protected Member Functions inherited from G4VProcess
void	SubtractNumberOfInteractionLengthLeft (G4double previousStepSize)

void	ClearNumberOfInteractionLengthLeft ()

Additional Inherited Members
Static Public Member Functions inherited from G4VProcess
static const G4String &	GetProcessTypeName (G4ProcessType)

Detailed Description

Definition at line 91 of file G4VMultipleScattering.hh.

Constructor & Destructor Documentation

◆ G4VMultipleScattering()

G4VMultipleScattering::G4VMultipleScattering	(	const G4String &	name = `"msc"`,
		G4ProcessType	type = `fElectromagnetic`
	)

Definition at line 86 of file G4VMultipleScattering.cc.

                                                                               :
  G4VContinuousDiscreteProcess("msc", fElectromagnetic),
  numberOfModels(0),
  firstParticle(0),
  currParticle(0),
  stepLimit(fUseSafety),
  skin(1.0),
  facrange(0.04),
  facgeom(2.5),
  latDisplasment(true),
  isIon(false)
{
  SetVerboseLevel(1);
  SetProcessSubType(fMultipleScattering);
  if("ionmsc" == name) { firstParticle = G4GenericIon::GenericIon(); }
 
  geomMin = 1.e-6*CLHEP::mm;
  lowestKinEnergy = 1*eV;
 
  // default limit on polar angle
  polarAngleLimit = 0.0;
 
  physStepLimit = gPathLength = tPathLength = 0.0;
  fIonisation = 0;
 
  pParticleChange = &fParticleChange;
  safetyHelper = 0;
  fPositionChanged = false;
  isActive = false;
 
  modelManager = new G4EmModelManager();
  emManager = G4LossTableManager::Instance();
  emManager->Register(this);
 
  warn = 0;
}

◆ ~G4VMultipleScattering()

G4VMultipleScattering::~G4VMultipleScattering ( )

virtual

Definition at line 125 of file G4VMultipleScattering.cc.

{
  if(1 < verboseLevel) {
    G4cout << "G4VMultipleScattering destruct " << GetProcessName() 
       << G4endl;
  }
  delete modelManager;
  emManager->DeRegister(this);
}

Member Function Documentation

◆ AddEmModel()

void G4VMultipleScattering::AddEmModel	(	G4int	order,
		G4VEmModel *	p,
		const G4Region *	region = `0`
	)

Definition at line 137 of file G4VMultipleScattering.cc.

{
  G4VEmFluctuationModel* fm = 0;
  modelManager->AddEmModel(order, p, fm, region);
  if(p) { p->SetParticleChange(pParticleChange); }
}

◆ AlongStepDoIt()

G4VParticleChange * G4VMultipleScattering::AlongStepDoIt	(	const G4Track &	track,
		const G4Step &	step
	)

virtual

Reimplemented from G4VContinuousDiscreteProcess.

Definition at line 403 of file G4VMultipleScattering.cc.

{
  fParticleChange.ProposeMomentumDirection(step.GetPostStepPoint()->GetMomentumDirection());
  fNewPosition = step.GetPostStepPoint()->GetPosition();
  fParticleChange.ProposePosition(fNewPosition);
  fPositionChanged = false;
 
  G4double geomLength = step.GetStepLength();
 
  // very small step - no msc
  if(!isActive) {
    tPathLength = geomLength;
 
    // sample msc
  } else {
    G4double range = 
      currentModel->GetRange(currParticle,track.GetKineticEnergy(),
                 track.GetMaterialCutsCouple());
 
    G4double trueLength = currentModel->ComputeTrueStepLength(geomLength);
 
    
    // protection against wrong t->g->t conversion
    //    if(trueLength > tPathLength) 
    /*
    if(currParticle->GetPDGMass() > GeV)    
    G4cout << "G4VMsc::AlongStepDoIt: GeomLength= " 
       << geomLength 
       << " trueLenght= " << trueLength 
       << " tPathLength= " << tPathLength
       << " dr= " << range - trueLength 
       << " ekin= " << track.GetKineticEnergy() << G4endl;
    */
    if (trueLength <= physStepLimit) {
      tPathLength = trueLength; 
    } else {
      tPathLength = physStepLimit - 0.5*geomMin; 
    }
 
    // do not sample scattering at the last or at a small step
    if(tPathLength + geomMin < range && tPathLength > geomMin) {
 
      G4double preSafety = step.GetPreStepPoint()->GetSafety();
      G4double postSafety= preSafety - geomLength; 
      G4bool safetyRecomputed = false;
      if( postSafety < geomMin ) {
    safetyRecomputed = true;
    postSafety = currentModel->ComputeSafety(fNewPosition,0.0); 
      } 
      G4ThreeVector displacement = 
    currentModel->SampleScattering(track.GetDynamicParticle(),postSafety);
 
      G4double r2 = displacement.mag2();
 
      //G4cout << "R= " << sqrt(r2) << " postSafety= " << postSafety << G4endl;
 
      // make correction for displacement
      if(r2 > 0.0) {
 
    fPositionChanged = true;
        G4double fac = 1.0;
 
    // displaced point is definitely within the volume
    if(r2 > postSafety*postSafety) {
          if(!safetyRecomputed) {
        postSafety = currentModel->ComputeSafety(fNewPosition, 0.0);
      } 
      // add a factor which ensure numerical stability
      if(r2 > postSafety*postSafety) { fac = 0.99*postSafety/std::sqrt(r2); }
    }
    // compute new endpoint of the Step
    fNewPosition += fac*displacement;
    //safetyHelper->ReLocateWithinVolume(fNewPosition);
      }
    }
  }
  fParticleChange.ProposeTrueStepLength(tPathLength);
  //fParticleChange.ProposePosition(fNewPosition);
  return &fParticleChange;
}

◆ AlongStepGetPhysicalInteractionLength()

G4double G4VMultipleScattering::AlongStepGetPhysicalInteractionLength	(	const G4Track &	track,
		G4double	previousStepSize,
		G4double	currentMinimalStep,
		G4double &	currentSafety,
		G4GPILSelection *	selection
	)

virtual

Reimplemented from G4VContinuousDiscreteProcess.

Definition at line 346 of file G4VMultipleScattering.cc.

{
  // get Step limit proposed by the process
  *selection = NotCandidateForSelection;
  physStepLimit = gPathLength = tPathLength = currentMinimalStep;
 
  G4double ekin = track.GetKineticEnergy();
  // isIon flag is used only to select a model
  if(isIon) { 
    ekin *= proton_mass_c2/track.GetParticleDefinition()->GetPDGMass(); 
  }
  
  // select new model
  if(1 < numberOfModels) {
    currentModel = static_cast<G4VMscModel*>(
      SelectModel(ekin,track.GetMaterialCutsCouple()->GetIndex()));
  }
 
  // step limit
  if(currentModel->IsActive(ekin) && gPathLength >= geomMin 
     && ekin >= lowestKinEnergy) {
    isActive = true;
    tPathLength = currentModel->ComputeTruePathLengthLimit(track, gPathLength);
    if (tPathLength < physStepLimit) { 
      *selection = CandidateForSelection; 
    }
  } else { isActive = false; }
  /*  
  if(currParticle->GetPDGMass() > GeV)    
  G4cout << "MSC::AlongStepGPIL: Ekin= " << ekin
     << " gPathLength= " << gPathLength
     << " tPathLength= " << tPathLength
     << " currentMinimalStep= " << currentMinimalStep
     << " isActive " << isActive << G4endl;
  */
  return gPathLength;
}

Referenced by GetContinuousStepLimit().

◆ BuildPhysicsTable()

void G4VMultipleScattering::BuildPhysicsTable ( const G4ParticleDefinition & part )

virtual

Reimplemented from G4VProcess.

Definition at line 267 of file G4VMultipleScattering.cc.

{
  G4String num = part.GetParticleName();
  if(1 < verboseLevel) {
    G4cout << "### G4VMultipleScattering::BuildPhysicsTable() for "
           << GetProcessName()
           << " and particle " << num
           << G4endl;
  }
 
  emManager->BuildPhysicsTable(firstParticle);
 
  // explicitly defined printout by particle name
  if(1 < verboseLevel || 
     (0 < verboseLevel && (num == "e-" || 
               num == "e+"    || num == "mu+" || 
               num == "mu-"   || num == "proton"|| 
               num == "pi+"   || num == "pi-" || 
               num == "kaon+" || num == "kaon-" || 
               num == "alpha" || num == "anti_proton" || 
               num == "GenericIon")))
    { 
      G4cout << G4endl << GetProcessName() 
         << ":   for " << num
         << "    SubType= " << GetProcessSubType() 
         << G4endl;
      PrintInfo();
      modelManager->DumpModelList(verboseLevel);
    }
 
  if(1 < verboseLevel) {
    G4cout << "### G4VMultipleScattering::BuildPhysicsTable() done for "
           << GetProcessName()
           << " and particle " << num
           << G4endl;
  }
}

◆ ContinuousStepLimit()

G4double G4VMultipleScattering::ContinuousStepLimit	(	const G4Track &	track,
		G4double	previousStepSize,
		G4double	currentMinimalStep,
		G4double &	currentSafety
	)

Definition at line 514 of file G4VMultipleScattering.cc.

{
  return GetContinuousStepLimit(track,previousStepSize,currentMinimalStep,
                currentSafety);
}

◆ EmModel()

G4VMscModel * G4VMultipleScattering::EmModel ( G4int index = 1 )

Definition at line 186 of file G4VMultipleScattering.cc.

{
  G4VMscModel* p = 0;
  if(index >= 0 && index <  G4int(mscModels.size())) { p = mscModels[index]; }
  return p;
}

Referenced by G4MuMultipleScattering::InitialiseProcess(), G4eMultipleScattering::InitialiseProcess(), and G4hMultipleScattering::InitialiseProcess().

◆ FirstParticle()

const G4ParticleDefinition * G4VMultipleScattering::FirstParticle ( ) const

inline

Definition at line 402 of file G4VMultipleScattering.hh.

{
  return firstParticle;
}

◆ GeomFactor()

G4double G4VMultipleScattering::GeomFactor ( ) const

inline

Definition at line 357 of file G4VMultipleScattering.hh.

{
  return facgeom;
}

Referenced by G4AdjointhMultipleScattering::InitialiseProcess(), PreparePhysicsTable(), G4eMultipleScattering::PrintInfo(), and G4hMultipleScattering::PrintInfo().

◆ GetContinuousStepLimit()

G4double G4VMultipleScattering::GetContinuousStepLimit	(	const G4Track &	track,
		G4double	previousStepSize,
		G4double	currentMinimalStep,
		G4double &	currentSafety
	)

protectedvirtual

Implements G4VContinuousDiscreteProcess.

Definition at line 499 of file G4VMultipleScattering.cc.

{
  G4GPILSelection selection = NotCandidateForSelection;
  G4double x = AlongStepGetPhysicalInteractionLength(track,previousStepSize,
                             currentMinimalStep,
                             currentSafety, &selection);
  return x;
}

Referenced by ContinuousStepLimit().

◆ GetMeanFreePath()

G4double G4VMultipleScattering::GetMeanFreePath	(	const G4Track &	track,
		G4double	,
		G4ForceCondition *	condition
	)

protectedvirtual

Implements G4VContinuousDiscreteProcess.

Definition at line 526 of file G4VMultipleScattering.cc.

{
  *condition = Forced;
  return DBL_MAX;
}

◆ GetModelByIndex()

G4VEmModel * G4VMultipleScattering::GetModelByIndex	(	G4int	idx = `0`,
		G4bool	ver = `false`
	)		const

Definition at line 196 of file G4VMultipleScattering.cc.

{
  return modelManager->GetModel(idx, ver);
}

◆ InitialiseProcess()

virtual void G4VMultipleScattering::InitialiseProcess ( const G4ParticleDefinition * )

protectedpure virtual

Implemented in G4AdjointhMultipleScattering, G4MuMultipleScattering, G4eMultipleScattering, and G4hMultipleScattering.

Referenced by PreparePhysicsTable().

◆ IsApplicable()

virtual G4bool G4VMultipleScattering::IsApplicable ( const G4ParticleDefinition & p )

pure virtual

Reimplemented from G4VProcess.

Implemented in G4AdjointhMultipleScattering, G4MuMultipleScattering, G4eMultipleScattering, and G4hMultipleScattering.

◆ LateralDisplasmentFlag()

G4bool G4VMultipleScattering::LateralDisplasmentFlag ( ) const

inline

Definition at line 314 of file G4VMultipleScattering.hh.

{
  return latDisplasment;
}

Referenced by G4AdjointhMultipleScattering::InitialiseProcess(), PreparePhysicsTable(), G4AdjointhMultipleScattering::PrintInfo(), G4MuMultipleScattering::PrintInfo(), G4eMultipleScattering::PrintInfo(), and G4hMultipleScattering::PrintInfo().

◆ Model()

G4VMscModel * G4VMultipleScattering::Model ( G4int index = 1 )

Definition at line 162 of file G4VMultipleScattering.cc.

{
  ++warn;
  if(warn < 10) { 
    G4cout << "### G4VMultipleScattering::Model is obsolete method "
       << "and will be removed for the next release." << G4endl;
    G4cout << "    Please, use EmModel" << G4endl;
  } 
  G4VMscModel* p = 0;
  if(index >= 0 && index <  G4int(mscModels.size())) { p = mscModels[index]; }
  return p;
}

◆ PolarAngleLimit()

G4double G4VMultipleScattering::PolarAngleLimit ( ) const

inline

Definition at line 371 of file G4VMultipleScattering.hh.

{
  return polarAngleLimit;
}

Referenced by G4MuMultipleScattering::PrintInfo().

◆ PostStepDoIt()

G4VParticleChange * G4VMultipleScattering::PostStepDoIt	(	const G4Track &	track,
		const G4Step &
	)

virtual

Reimplemented from G4VContinuousDiscreteProcess.

Definition at line 487 of file G4VMultipleScattering.cc.

{
  fParticleChange.Initialize(track);  
  if(fPositionChanged) { 
    safetyHelper->ReLocateWithinVolume(fNewPosition);
    fParticleChange.ProposePosition(fNewPosition); 
  }
  return &fParticleChange;
}

◆ PostStepGetPhysicalInteractionLength()

G4double G4VMultipleScattering::PostStepGetPhysicalInteractionLength	(	const G4Track &	,
		G4double	previousStepSize,
		G4ForceCondition *	condition
	)

virtual

Reimplemented from G4VContinuousDiscreteProcess.

Definition at line 392 of file G4VMultipleScattering.cc.

{
  *condition = Forced;
  //*condition = NotForced;
  return DBL_MAX;
}

◆ PreparePhysicsTable()

void G4VMultipleScattering::PreparePhysicsTable ( const G4ParticleDefinition & part )

virtual

Reimplemented from G4VProcess.

Definition at line 204 of file G4VMultipleScattering.cc.

{
  if(!firstParticle) { firstParticle = &part; }
  if(part.GetParticleType() == "nucleus") {
    SetStepLimitType(fMinimal);
    SetLateralDisplasmentFlag(false);
    SetRangeFactor(0.2);
    if(&part == G4GenericIon::GenericIon()) { firstParticle = &part; }
    isIon = true; 
  }
 
  emManager->PreparePhysicsTable(&part, this);
  currParticle = 0;
 
  if(1 < verboseLevel) {
    G4cout << "### G4VMultipleScattering::PrepearPhysicsTable() for "
           << GetProcessName()
           << " and particle " << part.GetParticleName()
           << " local particle " << firstParticle->GetParticleName()
       << " isIon= " << isIon 
           << G4endl;
  }
 
  if(firstParticle == &part) {
 
    InitialiseProcess(firstParticle);
 
    // initialisation of models
    numberOfModels = modelManager->NumberOfModels();
    for(G4int i=0; i<numberOfModels; ++i) {
      G4VMscModel* msc = static_cast<G4VMscModel*>(modelManager->GetModel(i));
      msc->SetIonisation(0, firstParticle);
      if(0 == i) { currentModel = msc; }
      if(isIon) {
    msc->SetStepLimitType(fMinimal);
    msc->SetLateralDisplasmentFlag(false);
    msc->SetRangeFactor(0.2);
      } else {
    msc->SetStepLimitType(StepLimitType());
    msc->SetLateralDisplasmentFlag(LateralDisplasmentFlag());
    msc->SetSkin(Skin());
    msc->SetRangeFactor(RangeFactor());
    msc->SetGeomFactor(GeomFactor());
      }
      msc->SetPolarAngleLimit(polarAngleLimit);
      G4double emax = 
    std::min(msc->HighEnergyLimit(),emManager->MaxKinEnergy());
      msc->SetHighEnergyLimit(emax);
    }
 
    modelManager->Initialise(firstParticle, G4Electron::Electron(), 
                 10.0, verboseLevel);
 
    if(!safetyHelper) {
      safetyHelper = G4TransportationManager::GetTransportationManager()
    ->GetSafetyHelper();
      safetyHelper->InitialiseHelper();
    }
  }
}

◆ PrintInfo()

virtual void G4VMultipleScattering::PrintInfo ( )

pure virtual

Implemented in G4AdjointhMultipleScattering, G4MuMultipleScattering, G4eMultipleScattering, and G4hMultipleScattering.

Referenced by BuildPhysicsTable(), and PrintInfoDefinition().

◆ PrintInfoDefinition()

void G4VMultipleScattering::PrintInfoDefinition ( )

Definition at line 307 of file G4VMultipleScattering.cc.

{
  if (0 < verboseLevel) {
    G4cout << G4endl << GetProcessName() 
       << ":   for " << firstParticle->GetParticleName()
       << "    SubType= " << GetProcessSubType() 
       << G4endl;
    PrintInfo();
    modelManager->DumpModelList(verboseLevel);
  }
}

◆ RangeFactor()

G4double G4VMultipleScattering::RangeFactor ( ) const

inline

Definition at line 343 of file G4VMultipleScattering.hh.

{
  return facrange;
}

Referenced by G4AdjointhMultipleScattering::InitialiseProcess(), PreparePhysicsTable(), G4AdjointhMultipleScattering::PrintInfo(), G4MuMultipleScattering::PrintInfo(), G4eMultipleScattering::PrintInfo(), and G4hMultipleScattering::PrintInfo().

◆ RetrievePhysicsTable()

G4bool G4VMultipleScattering::RetrievePhysicsTable	(	const G4ParticleDefinition *	,
		const G4String &	directory,
		G4bool	ascii
	)

virtual

Reimplemented from G4VProcess.

Definition at line 574 of file G4VMultipleScattering.cc.

{
  return true;
}

◆ SelectModel()

G4VEmModel * G4VMultipleScattering::SelectModel	(	G4double	kinEnergy,
		size_t	idx
	)

inline

Definition at line 307 of file G4VMultipleScattering.hh.

{
  return modelManager->SelectModel(kinEnergy, coupleIndex);
}

Referenced by AlongStepGetPhysicalInteractionLength().

◆ SetEmModel()

void G4VMultipleScattering::SetEmModel	(	G4VMscModel *	p,
		G4int	index = `1`
	)

Definition at line 177 of file G4VMultipleScattering.cc.

{
  G4int n = mscModels.size();
  if(index >= n) { for(G4int i=n; i<=index; ++i) { mscModels.push_back(0); } }
  mscModels[index] = p;
}

Referenced by G4MuMultipleScattering::InitialiseProcess(), G4eMultipleScattering::InitialiseProcess(), and G4hMultipleScattering::InitialiseProcess().

◆ SetGeomFactor()

void G4VMultipleScattering::SetGeomFactor ( G4double val )

inline

Definition at line 364 of file G4VMultipleScattering.hh.

{
  if(val > 0.0) facgeom = val;
}

◆ SetIonisation()

void G4VMultipleScattering::SetIonisation ( G4VEnergyLossProcess * p )

Definition at line 583 of file G4VMultipleScattering.cc.

{
  for(G4int i=0; i<numberOfModels; ++i) {
    G4VMscModel* msc = static_cast<G4VMscModel*>(modelManager->GetModel(i));
    msc->SetIonisation(p, firstParticle);
  }
}

◆ SetLateralDisplasmentFlag()

void G4VMultipleScattering::SetLateralDisplasmentFlag ( G4bool val )

inline

Definition at line 321 of file G4VMultipleScattering.hh.

{
  latDisplasment = val;
}

Referenced by G4AdjointhMultipleScattering::InitialiseProcess(), and PreparePhysicsTable().

◆ SetModel()

void G4VMultipleScattering::SetModel	(	G4VMscModel *	p,
		G4int	index = `1`
	)

Definition at line 147 of file G4VMultipleScattering.cc.

{
  ++warn;
  if(warn < 10) { 
    G4cout << "### G4VMultipleScattering::SetModel is obsolete method "
       << "and will be removed for the next release." << G4endl;
    G4cout << "    Please, use SetEmModel" << G4endl;
  } 
  G4int n = mscModels.size();
  if(index >= n) { for(G4int i=n; i<=index; ++i) {mscModels.push_back(0);} }
  mscModels[index] = p;
}

◆ SetPolarAngleLimit()

void G4VMultipleScattering::SetPolarAngleLimit ( G4double val )

inline

Definition at line 378 of file G4VMultipleScattering.hh.

{
  if(val < 0.0)            { polarAngleLimit = 0.0; }
  else if(val > CLHEP::pi) { polarAngleLimit = CLHEP::pi; }
  else                     { polarAngleLimit = val; }
}

◆ SetRangeFactor()

void G4VMultipleScattering::SetRangeFactor ( G4double val )

inline

Definition at line 350 of file G4VMultipleScattering.hh.

{
  if(val > 0.0) facrange = val;
}

Referenced by PreparePhysicsTable().

◆ SetSkin()

void G4VMultipleScattering::SetSkin ( G4double val )

inline

Definition at line 335 of file G4VMultipleScattering.hh.

{
  if(val < 1.0) { skin = 0.0; }
  else          { skin = val; }
}

◆ SetStepLimitType()

void G4VMultipleScattering::SetStepLimitType ( G4MscStepLimitType val )

inline

Definition at line 394 of file G4VMultipleScattering.hh.

{
  stepLimit = val;
  if(val == fMinimal) { facrange = 0.2; }
}

Referenced by G4EmDNAPhysics::ConstructProcess(), G4EmDNAPhysicsChemistry::ConstructProcess(), G4EmLivermorePhysics::ConstructProcess(), G4EmLivermorePolarizedPhysics::ConstructProcess(), G4EmLowEPPhysics::ConstructProcess(), G4EmPenelopePhysics::ConstructProcess(), G4EmStandardPhysics_option1::ConstructProcess(), G4EmStandardPhysics_option2::ConstructProcess(), G4EmStandardPhysics_option3::ConstructProcess(), G4EmStandardPhysics_option4::ConstructProcess(), G4AdjointhMultipleScattering::G4AdjointhMultipleScattering(), G4hMultipleScattering::G4hMultipleScattering(), G4MuMultipleScattering::G4MuMultipleScattering(), G4AdjointhMultipleScattering::InitialiseProcess(), and PreparePhysicsTable().

◆ Skin()

G4double G4VMultipleScattering::Skin ( ) const

inline

Definition at line 328 of file G4VMultipleScattering.hh.

{
  return skin;
}

Referenced by G4AdjointhMultipleScattering::InitialiseProcess(), PreparePhysicsTable(), G4AdjointhMultipleScattering::PrintInfo(), G4eMultipleScattering::PrintInfo(), and G4hMultipleScattering::PrintInfo().

◆ StartTracking()

void G4VMultipleScattering::StartTracking ( G4Track * track )

virtual

Reimplemented from G4VProcess.

Definition at line 321 of file G4VMultipleScattering.cc.

{
  G4VEnergyLossProcess* eloss = 0;
  if(track->GetParticleDefinition() != currParticle) {
    currParticle = track->GetParticleDefinition();
    fIonisation = emManager->GetEnergyLossProcess(currParticle);
    eloss = fIonisation;
  }
  // one model
  if(1 == numberOfModels) {
    currentModel->StartTracking(track);
    if(eloss) { currentModel->SetIonisation(fIonisation, currParticle); }
 
    // many models
  } else { 
    for(G4int i=0; i<numberOfModels; ++i) {
      G4VMscModel* msc = static_cast<G4VMscModel*>(modelManager->GetModel(i));
      msc->StartTracking(track);
      if(eloss) { msc->SetIonisation(fIonisation, currParticle); }
    }
  }
}

◆ StepLimitType()

G4MscStepLimitType G4VMultipleScattering::StepLimitType ( ) const

inline

Definition at line 387 of file G4VMultipleScattering.hh.

{
  return stepLimit;
}

Referenced by G4AdjointhMultipleScattering::InitialiseProcess(), PreparePhysicsTable(), G4AdjointhMultipleScattering::PrintInfo(), G4MuMultipleScattering::PrintInfo(), G4eMultipleScattering::PrintInfo(), and G4hMultipleScattering::PrintInfo().

◆ StorePhysicsTable()

G4bool G4VMultipleScattering::StorePhysicsTable	(	const G4ParticleDefinition *	part,
		const G4String &	directory,
		G4bool	ascii = `false`
	)

virtual

Reimplemented from G4VProcess.

Definition at line 536 of file G4VMultipleScattering.cc.

{
  G4bool yes = true;
  if(part != firstParticle) { return yes; }
  G4int nmod = modelManager->NumberOfModels();
  const G4String ss[4] = {"1","2","3","4"};
  for(G4int i=0; i<nmod; ++i) {
    G4VEmModel* msc = modelManager->GetModel(i);
    yes = true;
    G4PhysicsTable* table = msc->GetCrossSectionTable();
    if (table) {
      G4int j = std::min(i,3); 
      G4String name = 
    GetPhysicsTableFileName(part,directory,"LambdaMod"+ss[j],ascii);
      yes = table->StorePhysicsTable(name,ascii);
 
      if ( yes ) {
    if ( verboseLevel>0 ) {
      G4cout << "Physics table are stored for " << part->GetParticleName()
         << " and process " << GetProcessName()
         << " with a name <" << name << "> " << G4endl;
    }
      } else {
    G4cout << "Fail to store Physics Table for " << part->GetParticleName()
           << " and process " << GetProcessName()
           << " in the directory <" << directory
           << "> " << G4endl;
      }
    }
  }
  return yes;
}

Member Data Documentation

◆ fParticleChange

G4ParticleChangeForMSC G4VMultipleScattering::fParticleChange

protected

Definition at line 283 of file G4VMultipleScattering.hh.

Referenced by AlongStepDoIt(), G4VMultipleScattering(), and PostStepDoIt().

◆ valueGPILSelectionMSC

G4GPILSelection G4VMultipleScattering::valueGPILSelectionMSC

protected

Definition at line 282 of file G4VMultipleScattering.hh.

The documentation for this class was generated from the following files:

Protected Attributes
G4GPILSelection	valueGPILSelectionMSC

G4ParticleChangeForMSC	fParticleChange

Protected Attributes inherited from G4VProcess
const G4ProcessManager *	aProcessManager

G4VParticleChange *	pParticleChange

G4ParticleChange	aParticleChange

G4double	theNumberOfInteractionLengthLeft

G4double	currentInteractionLength

G4double	theInitialNumberOfInteractionLength

G4String	theProcessName

G4String	thePhysicsTableFileName

G4ProcessType	theProcessType

G4int	theProcessSubType

G4double	thePILfactor

G4bool	enableAtRestDoIt

G4bool	enableAlongStepDoIt

G4bool	enablePostStepDoIt

G4int	verboseLevel

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4VMultipleScattering()

◆ ~G4VMultipleScattering()

Member Function Documentation

◆ AddEmModel()

◆ AlongStepDoIt()

◆ AlongStepGetPhysicalInteractionLength()

◆ BuildPhysicsTable()

◆ ContinuousStepLimit()

◆ EmModel()

◆ FirstParticle()

◆ GeomFactor()

◆ GetContinuousStepLimit()

◆ GetMeanFreePath()

◆ GetModelByIndex()

◆ InitialiseProcess()

◆ IsApplicable()

◆ LateralDisplasmentFlag()

◆ Model()

◆ PolarAngleLimit()

◆ PostStepDoIt()

◆ PostStepGetPhysicalInteractionLength()

◆ PreparePhysicsTable()

◆ PrintInfo()

◆ PrintInfoDefinition()

◆ RangeFactor()

◆ RetrievePhysicsTable()

◆ SelectModel()

◆ SetEmModel()

◆ SetGeomFactor()

◆ SetIonisation()

◆ SetLateralDisplasmentFlag()

◆ SetModel()

◆ SetPolarAngleLimit()

◆ SetRangeFactor()

◆ SetSkin()

◆ SetStepLimitType()

◆ Skin()

◆ StartTracking()

◆ StepLimitType()

◆ StorePhysicsTable()

Member Data Documentation

◆ fParticleChange

◆ valueGPILSelectionMSC