#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	SetEmModel (G4VMscModel *, G4int index=1)

G4VMscModel *	EmModel (G4int index=1) const

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)

G4double	LowestKinEnergy () const

void	SetLowestKinEnergy (G4double val)

const G4ParticleDefinition *	FirstParticle () const

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

	G4VContinuousDiscreteProcess (G4VContinuousDiscreteProcess &)

virtual	~G4VContinuousDiscreteProcess ()

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

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)

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	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

virtual void	SetMasterProcess (G4VProcess *masterP)

const G4VProcess *	GetMasterProcess () const

virtual void	BuildWorkerPhysicsTable (const G4ParticleDefinition &part)

virtual void	PrepareWorkerPhysicsTable (const G4ParticleDefinition &)

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
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	(	const G4String &	name = `"msc"`,
		G4ProcessType	type = `fElectromagnetic`
	)

Definition at line 90 of file G4VMultipleScattering.cc.

References python.hepunit::eV, fMultipleScattering, fParticleChange, G4GenericIon::GenericIon(), G4LossTableManager::Instance(), G4VProcess::pParticleChange, G4LossTableManager::Register(), G4VProcess::SetProcessSubType(), and G4VProcess::SetVerboseLevel().

                                        :
   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 = 0.01*CLHEP::nm;
   lowestKinEnergy = 10*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 ( )

virtual

Definition at line 130 of file G4VMultipleScattering.cc.

References G4LossTableManager::DeRegister().

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

Member Function Documentation

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

Definition at line 144 of file G4VMultipleScattering.cc.

References G4EmModelManager::AddEmModel(), fm, G4VProcess::pParticleChange, and G4VEmModel::SetParticleChange().

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

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

virtual

Reimplemented from G4VContinuousDiscreteProcess.

Definition at line 457 of file G4VMultipleScattering.cc.

References G4SafetyHelper::ComputeSafety(), G4VMscModel::ComputeTrueStepLength(), fParticleChange, G4Track::GetKineticEnergy(), G4Track::GetMaterialCutsCouple(), G4StepPoint::GetMomentumDirection(), G4StepPoint::GetPosition(), G4Step::GetPostStepPoint(), G4Step::GetPreStepPoint(), G4VMscModel::GetRange(), G4StepPoint::GetSafety(), G4Step::GetStepLength(), CLHEP::Hep3Vector::mag2(), G4ParticleChangeForMSC::ProposeMomentumDirection(), G4ParticleChangeForMSC::ProposePosition(), G4VParticleChange::ProposeTrueStepLength(), and G4VMscModel::SampleScattering().

 {
   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());
 
     tPathLength = currentModel->ComputeTrueStepLength(geomLength);
   
     // protection against wrong t->g->t conversion
     /*
     if(currParticle->GetPDGMass() > GeV)    
     G4cout << "G4VMsc::AlongStepDoIt: GeomLength= " 
        << geomLength 
        << " tPathLength= " << tPathLength
        << " physStepLimit= " << physStepLimit
        << " dr= " << range - trueLength 
        << " ekin= " << track.GetKineticEnergy() << G4endl;
     */
     if (tPathLength > physStepLimit) {
       tPathLength = physStepLimit; 
     }
 
     // do not sample scattering at the last or at a small step
     if(tPathLength + geomMin < range && tPathLength > geomMin) {
 
       G4double preSafety = step.GetPreStepPoint()->GetSafety();
       G4double maxDisp   = (tPathLength + geomLength)*0.5; 
       G4double postSafety= preSafety - maxDisp; 
       G4bool safetyRecomputed = false;
       if(postSafety < maxDisp) {
     safetyRecomputed = true;
     postSafety = safetyHelper->ComputeSafety(fNewPosition,maxDisp); 
       } 
       G4ThreeVector displacement = currentModel->SampleScattering(
         step.GetPostStepPoint()->GetMomentumDirection(), 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) {
           G4double dispR = std::sqrt(r2);
           if(!safetyRecomputed) {
         postSafety = safetyHelper->ComputeSafety(fNewPosition, dispR);
       } 
 
       if(dispR > postSafety) { 
         fac = 0.99*postSafety/dispR; 
       }
     }
     // compute new endpoint of the Step
     fNewPosition += fac*displacement;
     //safetyHelper->ReLocateWithinVolume(fNewPosition);
       }
     }
   }
   fParticleChange.ProposeTrueStepLength(tPathLength);
   //fParticleChange.ProposePosition(fNewPosition);
   return &fParticleChange;
 }

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

virtual

Reimplemented from G4VContinuousDiscreteProcess.

Definition at line 391 of file G4VMultipleScattering.cc.

References CandidateForSelection, G4VMscModel::ComputeTruePathLengthLimit(), G4MaterialCutsCouple::GetIndex(), G4Track::GetKineticEnergy(), G4Track::GetMaterialCutsCouple(), G4Track::GetParticleDefinition(), G4ParticleDefinition::GetPDGMass(), G4VEmModel::IsActive(), NotCandidateForSelection, python.hepunit::proton_mass_c2, and SelectModel().

Referenced by GetContinuousStepLimit().

 {
   // get Step limit proposed by the process
   *selection = NotCandidateForSelection;
   physStepLimit = gPathLength = tPathLength = currentMinimalStep;
 
   G4double ekin = track.GetKineticEnergy();
   /*
   G4cout << "MSC::AlongStepGPIL: Ekin= " << ekin
          << "  " << currParticle->GetParticleName() 
      << " currMod " << currentModel 
      << G4endl;
   */
   // 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
          << "  " << currParticle->GetParticleName()
      << " gPathLength= " << gPathLength
      << " tPathLength= " << tPathLength
      << " currentMinimalStep= " << currentMinimalStep
      << " isActive " << isActive << G4endl;
   */
   return gPathLength;
 }

void G4VMultipleScattering::BuildPhysicsTable ( const G4ParticleDefinition & part )

virtual

Reimplemented from G4VProcess.

Definition at line 268 of file G4VMultipleScattering.cc.

References G4LossTableManager::BuildPhysicsTable(), G4EmModelManager::DumpModelList(), G4cout, G4endl, G4VEmModel::GetCrossSectionTable(), G4VProcess::GetMasterProcess(), GetModelByIndex(), G4ParticleDefinition::GetParticleName(), G4VProcess::GetProcessName(), G4VProcess::GetProcessSubType(), G4LossTableManager::Instance(), G4LossTableManager::IsMaster(), G4EmModelManager::NumberOfModels(), PrintInfo(), G4VEmModel::SetCrossSectionTable(), and G4VProcess::verboseLevel.

 {
   G4String num = part.GetParticleName();
   if(1 < verboseLevel) {
     G4cout << "### G4VMultipleScattering::BuildPhysicsTable() for "
            << GetProcessName()
            << " and particle " << num
        << " IsMaster= " << G4LossTableManager::Instance()->IsMaster()
            << G4endl;
   }
   G4bool master = true;
   const G4VMultipleScattering* masterProcess = 
     static_cast<const G4VMultipleScattering*>(GetMasterProcess()); 
   if(masterProcess != this) { master = false; }
 
   if(firstParticle == &part) { 
     /*    
     G4cout << "### G4VMultipleScattering::BuildPhysicsTable() for "
            << GetProcessName()
            << " and particle " << num
        << " IsMaster= " << G4LossTableManager::Instance()->IsMaster()
        << "  " << this
            << G4endl;
     */
     emManager->BuildPhysicsTable(firstParticle);
 
     if(!master) {
       // initialisation of models
       G4bool printing = true;
       numberOfModels = modelManager->NumberOfModels();
       /*
     G4cout << "### G4VMultipleScattering::SlaveBuildPhysicsTable() for "
     << GetProcessName()
     << " and particle " << num
     << " Nmod= " << numberOfModels << "  " << this
     << G4endl;
       */
       for(G4int i=0; i<numberOfModels; ++i) {
     G4VMscModel* msc = 
       static_cast<G4VMscModel*>(GetModelByIndex(i, printing));
     G4VMscModel* msc0= 
       static_cast<G4VMscModel*>(masterProcess->GetModelByIndex(i,printing));
     msc->SetCrossSectionTable(msc0->GetCrossSectionTable(), false);
       }
     }
 
   }
 
   // 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;
   }
 }

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

Definition at line 573 of file G4VMultipleScattering.cc.

References GetContinuousStepLimit().

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

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

Definition at line 163 of file G4VMultipleScattering.cc.

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

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

const G4ParticleDefinition * G4VMultipleScattering::FirstParticle ( ) const

inline

Definition at line 413 of file G4VMultipleScattering.hh.

 {
   return firstParticle;
 }

G4double G4VMultipleScattering::GeomFactor ( ) const

inline

Definition at line 354 of file G4VMultipleScattering.hh.

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

 {
   return facgeom;
 }

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

protectedvirtual

Implements G4VContinuousDiscreteProcess.

Definition at line 557 of file G4VMultipleScattering.cc.

References AlongStepGetPhysicalInteractionLength(), NotCandidateForSelection, and test::x.

Referenced by ContinuousStepLimit().

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

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

protectedvirtual

Implements G4VContinuousDiscreteProcess.

Definition at line 585 of file G4VMultipleScattering.cc.

References DBL_MAX, and Forced.

 {
   *condition = Forced;
   return DBL_MAX;
 }

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

Definition at line 173 of file G4VMultipleScattering.cc.

References G4EmModelManager::GetModel().

Referenced by BuildPhysicsTable(), SetIonisation(), and StartTracking().

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

virtual void G4VMultipleScattering::InitialiseProcess ( const G4ParticleDefinition * )

protectedpure virtual

Implemented in G4AdjointhMultipleScattering, G4hMultipleScattering, G4MuMultipleScattering, and G4eMultipleScattering.

Referenced by PreparePhysicsTable().

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

pure virtual

Reimplemented from G4VProcess.

Implemented in G4AdjointhMultipleScattering, G4hMultipleScattering, G4MuMultipleScattering, and G4eMultipleScattering.

G4bool G4VMultipleScattering::LateralDisplasmentFlag ( ) const

inline

Definition at line 311 of file G4VMultipleScattering.hh.

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

 {
   return latDisplasment;
 }

G4double G4VMultipleScattering::LowestKinEnergy ( ) const

inline

Definition at line 399 of file G4VMultipleScattering.hh.

 {
   return lowestKinEnergy;
 }

G4double G4VMultipleScattering::PolarAngleLimit ( ) const

inline

Definition at line 368 of file G4VMultipleScattering.hh.

Referenced by G4MuMultipleScattering::PrintInfo().

 {
   return polarAngleLimit;
 }

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

virtual

Reimplemented from G4VContinuousDiscreteProcess.

Definition at line 543 of file G4VMultipleScattering.cc.

References fParticleChange, G4ParticleChangeForMSC::Initialize(), G4ParticleChangeForMSC::ProposePosition(), and G4SafetyHelper::ReLocateWithinVolume().

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

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

virtual

Reimplemented from G4VContinuousDiscreteProcess.

Definition at line 446 of file G4VMultipleScattering.cc.

References DBL_MAX, and Forced.

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

void G4VMultipleScattering::PreparePhysicsTable ( const G4ParticleDefinition & part )

virtual

Reimplemented from G4VProcess.

Definition at line 181 of file G4VMultipleScattering.cc.

References G4Electron::Electron(), G4ParticleTable::FindParticle(), fMinimal, G4cout, G4endl, GeomFactor(), G4ProcessManager::GetAlongStepProcessVector(), G4VProcess::GetMasterProcess(), G4EmModelManager::GetModel(), G4ParticleDefinition::GetParticleName(), G4ParticleTable::GetParticleTable(), G4ParticleDefinition::GetParticleType(), G4ParticleDefinition::GetProcessManager(), G4VProcess::GetProcessName(), G4TransportationManager::GetSafetyHelper(), G4TransportationManager::GetTransportationManager(), G4VEmModel::HighEnergyLimit(), G4EmModelManager::Initialise(), G4SafetyHelper::InitialiseHelper(), InitialiseProcess(), LateralDisplasmentFlag(), G4LossTableManager::MaxKinEnergy(), G4INCL::Math::min(), n, G4EmModelManager::NumberOfModels(), eplot::pname, G4LossTableManager::PreparePhysicsTable(), RangeFactor(), G4VMscModel::SetGeomFactor(), G4VEmModel::SetHighEnergyLimit(), G4VMscModel::SetIonisation(), G4VMscModel::SetLateralDisplasmentFlag(), SetLateralDisplasmentFlag(), G4VEmModel::SetMasterThread(), G4VEmModel::SetPolarAngleLimit(), G4VMscModel::SetRangeFactor(), SetRangeFactor(), G4VMscModel::SetSkin(), G4VMscModel::SetStepLimitType(), SetStepLimitType(), G4ProcessVector::size(), Skin(), StepLimitType(), test::v, and G4VProcess::verboseLevel.

 {
   G4bool master = true;
   if(GetMasterProcess() != this) { master = false; }
 
   if(!firstParticle) { firstParticle = &part; }
   if(part.GetParticleType() == "nucleus") {
     SetStepLimitType(fMinimal);
     SetLateralDisplasmentFlag(false);
     SetRangeFactor(0.2);
     G4String pname = part.GetParticleName();
     if(pname != "deuteron" && pname != "triton" &&
        pname != "alpha+"   && pname != "helium" &&
        pname != "alpha"    && pname != "He3" &&
        pname != "hydrogen") {
 
       const G4ParticleDefinition* theGenericIon = 
     G4ParticleTable::GetParticleTable()->FindParticle("GenericIon");
      
       if(theGenericIon && firstParticle != theGenericIon) {
     G4ProcessManager* pm =  theGenericIon->GetProcessManager();
     G4ProcessVector* v = pm->GetAlongStepProcessVector();
     size_t n = v->size();
     for(size_t j=0; j<n; ++j) {
       if((*v)[j] == this) {
         firstParticle = theGenericIon;
         isIon = true; 
         break; 
       }
     }
       }
     }
   }
 
   emManager->PreparePhysicsTable(&part, this, master);
   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);
       msc->SetMasterThread(master);
       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();
     }
   }
 }

virtual void G4VMultipleScattering::PrintInfo ( )

pure virtual

Implemented in G4AdjointhMultipleScattering, G4hMultipleScattering, G4MuMultipleScattering, and G4eMultipleScattering.

Referenced by BuildPhysicsTable(), and PrintInfoDefinition().

void G4VMultipleScattering::PrintInfoDefinition ( )

Definition at line 344 of file G4VMultipleScattering.cc.

References G4EmModelManager::DumpModelList(), G4cout, G4endl, G4ParticleDefinition::GetParticleName(), G4VProcess::GetProcessName(), G4VProcess::GetProcessSubType(), PrintInfo(), and G4VProcess::verboseLevel.

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

G4double G4VMultipleScattering::RangeFactor ( ) const

inline

Definition at line 340 of file G4VMultipleScattering.hh.

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

 {
   return facrange;
 }

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

virtual

Reimplemented from G4VProcess.

Definition at line 635 of file G4VMultipleScattering.cc.

 {
   return true;
 }

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

inline

Definition at line 304 of file G4VMultipleScattering.hh.

References G4EmModelManager::SelectModel().

Referenced by AlongStepGetPhysicalInteractionLength().

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

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

Definition at line 154 of file G4VMultipleScattering.cc.

References n.

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

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

void G4VMultipleScattering::SetGeomFactor ( G4double val )

inline

Definition at line 361 of file G4VMultipleScattering.hh.

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

void G4VMultipleScattering::SetIonisation ( G4VEnergyLossProcess * p )

Definition at line 644 of file G4VMultipleScattering.cc.

References GetModelByIndex(), and G4VMscModel::SetIonisation().

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

void G4VMultipleScattering::SetLateralDisplasmentFlag ( G4bool val )

inline

Definition at line 318 of file G4VMultipleScattering.hh.

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

 {
   latDisplasment = val;
 }

void G4VMultipleScattering::SetLowestKinEnergy ( G4double val )

inline

Definition at line 406 of file G4VMultipleScattering.hh.

 {
   lowestKinEnergy = val;
 }

void G4VMultipleScattering::SetPolarAngleLimit ( G4double val )

inline

Definition at line 375 of file G4VMultipleScattering.hh.

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

void G4VMultipleScattering::SetRangeFactor ( G4double val )

inline

Definition at line 347 of file G4VMultipleScattering.hh.

Referenced by G4EmLivermorePhysics::ConstructProcess(), G4EmPenelopePhysics::ConstructProcess(), G4EmStandardPhysics_option4::ConstructProcess(), and PreparePhysicsTable().

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

void G4VMultipleScattering::SetSkin ( G4double val )

inline

Definition at line 332 of file G4VMultipleScattering.hh.

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

void G4VMultipleScattering::SetStepLimitType ( G4MscStepLimitType val )

inline

Definition at line 391 of file G4VMultipleScattering.hh.

References fMinimal.

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

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

G4double G4VMultipleScattering::Skin ( ) const

inline

Definition at line 325 of file G4VMultipleScattering.hh.

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

 {
   return skin;
 }

void G4VMultipleScattering::StartTracking ( G4Track * track )

virtual

Reimplemented from G4VProcess.

Definition at line 358 of file G4VMultipleScattering.cc.

References G4LossTableManager::GetEnergyLossProcess(), GetModelByIndex(), G4Track::GetParticleDefinition(), G4VMscModel::SetIonisation(), and G4VEmModel::StartTracking().

 {
   G4VEnergyLossProcess* eloss = 0;
   if(track->GetParticleDefinition() != currParticle) {
     currParticle = track->GetParticleDefinition();
     fIonisation = emManager->GetEnergyLossProcess(currParticle);
     eloss = fIonisation;
   }
   /*  
   G4cout << "G4VMultipleScattering::StartTracking Nmod= " << numberOfModels
      << "  " << currParticle->GetParticleName() 
      << " E(MeV)= " << track->GetKineticEnergy()
      << "  Ion= " << eloss << "  " << fIonisation << " IsMaster= " 
      << G4LossTableManager::Instance()->IsMaster() 
      << G4endl;
   */
   // 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*>(GetModelByIndex(i,true));
       msc->StartTracking(track);
       if(eloss) { msc->SetIonisation(fIonisation, currParticle); }
     }
   }
 } 

G4MscStepLimitType G4VMultipleScattering::StepLimitType ( ) const

inline

Definition at line 384 of file G4VMultipleScattering.hh.

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

 {
   return stepLimit;
 }

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

virtual

Reimplemented from G4VProcess.

Definition at line 595 of file G4VMultipleScattering.cc.

References G4cout, G4endl, G4VEmModel::GetCrossSectionTable(), G4EmModelManager::GetModel(), G4ParticleDefinition::GetParticleName(), G4VProcess::GetPhysicsTableFileName(), G4VProcess::GetProcessName(), G4INCL::Math::min(), G4EmModelManager::NumberOfModels(), G4PhysicsTable::StorePhysicsTable(), and G4VProcess::verboseLevel.

 {
   G4bool yes = true;
   if(part != firstParticle) { return yes; }
   G4int nmod = modelManager->NumberOfModels();
   static 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;
 }

Field Documentation

G4ParticleChangeForMSC G4VMultipleScattering::fParticleChange

protected

Definition at line 280 of file G4VMultipleScattering.hh.

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

G4GPILSelection G4VMultipleScattering::valueGPILSelectionMSC

protected

Definition at line 279 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

Member Function Documentation

Field Documentation