G4QDiscProcessMixer Class Reference

#include <G4QDiscProcessMixer.hh>

Inheritance diagram for G4QDiscProcessMixer:

Public Member Functions
	G4QDiscProcessMixer (const G4String &processName="Mixed Discrete Process", const G4ParticleDefinition *proj=G4Gamma::Gamma(), G4ProcessType pType=fHadronic)
	~G4QDiscProcessMixer ()
G4bool	IsApplicable (const G4ParticleDefinition &particle)
G4double	PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
G4double	GetMeanFreePath (const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *condition)
G4VParticleChange *	PostStepDoIt (const G4Track &aTrack, const G4Step &aStep)
void	AddDiscreteProcess (G4VDiscreteProcess *DP, G4double MaxE)

---

Detailed Description

Definition at line 53 of file G4QDiscProcessMixer.hh.

---

Constructor & Destructor Documentation

G4QDiscProcessMixer::G4QDiscProcessMixer	(	const G4String &	processName = "Mixed Discrete Process",
		const G4ParticleDefinition *	proj = G4Gamma::Gamma(),
		G4ProcessType	pType = fHadronic
	)

Definition at line 45 of file G4QDiscProcessMixer.cc.

References G4cout, G4endl, G4HadronicDeprecate, G4VProcess::GetProcessName(), and G4VProcess::verboseLevel.

                                                             :
  G4VDiscreteProcess(name, pType), DPParticle(particle)
{
  G4HadronicDeprecate("G4QDiscProcessMixer");

#ifdef debug
  G4cout<<"G4QDiscProcessMixer::Constructor is called processName="<<name<<G4endl;
#endif
  if (verboseLevel>0) G4cout<<GetProcessName()<<" process is created "<<G4endl;
}

G4QDiscProcessMixer::~G4QDiscProcessMixer

(

)

Definition at line 59 of file G4QDiscProcessMixer.cc.

{
  // Now the responsibility of deleting is deligated to the user, who created them
  //for_each(theDPVector.begin(), theDPVector.end(), DeleteDiscreteProcess());
}

---

Member Function Documentation

void G4QDiscProcessMixer::AddDiscreteProcess	(	G4VDiscreteProcess *	DP,
		G4double	MaxE
	)

Definition at line 65 of file G4QDiscProcessMixer.cc.

References FatalException, G4endl, and G4Exception().

Referenced by G4QProtonBuilder::Build(), G4QNeutronHPBuilder::Build(), G4QNeutronBuilder::Build(), and G4QInelasticCHIPS_HPBuilder::Build().

{
  static const G4double maxEn = 1.E8*megaelectronvolt; // Conditional INF
  if(!theDPVector.size()) // The first process in the DiscreteProcessVector (MaxE=INF)
  {
    std::pair<G4VDiscreteProcess*, G4double>* QDiscProc =
                                   new std::pair<G4VDiscreteProcess*, G4double>(DP, maxEn);
    theDPVector.push_back( QDiscProc );
  }
  else
  {
    if(ME < theDPVector[theDPVector.size()-1]->second)
    {
      std::pair<G4VDiscreteProcess*, G4double>* QDiscProc =
                                      new std::pair<G4VDiscreteProcess*, G4double>(DP, ME);
      theDPVector.push_back( QDiscProc );
    }
    else // Wrong Max Energy Order for the new process in the sequence of processes
    {
      // G4cerr<<"G4QDiscProcessMixer::AddDiscreteProcess:LastMaxE("<<theDPVector.size()-1
      //       <<")="<<theDPVector[theDPVector.size()-1]->second<<" <= MaxE="<<ME<<G4endl;
      // G4Exception("G4QDiscProcessMixer::AddDiscreteProcess: Wrong Max Energy Order");
      G4ExceptionDescription ed;
      ed << " LastMaxE(" << theDPVector.size()-1 << ")="
         << theDPVector[theDPVector.size()-1]->second << " <= MaxE=" << ME << G4endl;
      G4Exception("G4QDiscProcessMixer::AddDiscreteProcess()", "HAD_CHPS_0000",
                  FatalException, ed);
    }
  }
}

G4double G4QDiscProcessMixer::GetMeanFreePath	(	const G4Track &	aTrack,
		G4double	previousStepSize,
		G4ForceCondition *	condition
	)			[virtual]

Implements G4VDiscreteProcess.

Definition at line 125 of file G4QDiscProcessMixer.cc.

References DBL_MAX.

{
  return DBL_MAX;
}

G4bool G4QDiscProcessMixer::IsApplicable

(

const G4ParticleDefinition &

particle

)

[virtual]

Reimplemented from G4VProcess.

Definition at line 96 of file G4QDiscProcessMixer.cc.

References G4cout, G4endl, and G4ParticleDefinition::GetParticleName().

{
#ifdef debug
  G4cout<<"G4QDiscProcessMixer::IsApplicable: part="<<particle.GetParticleName()<<" = "
        <<DPParticle->GetParticleName()<<G4endl;
#endif
  if(particle == *DPParticle) return true;
  return false;
}

G4VParticleChange * G4QDiscProcessMixer::PostStepDoIt	(	const G4Track &	aTrack,
		const G4Step &	aStep
	)			[virtual]

Reimplemented from G4VDiscreteProcess.

Definition at line 130 of file G4QDiscProcessMixer.cc.

References G4cout, G4endl, G4Track::GetDynamicParticle(), G4DynamicParticle::GetKineticEnergy(), and G4VDiscreteProcess::PostStepDoIt().

{
  G4double kEn=Track.GetDynamicParticle()->GetKineticEnergy(); // Projectile kinetic energy
  G4int maxDP=theDPVector.size();
  if(maxDP) for(G4int i=maxDP-1; i>-1; i--) if(kEn < theDPVector[i]->second)
  {
#ifdef debug
    G4cout<<"G4QDPMix::PSDoIt: i="<<i<<",kEn="<<kEn<<" < "<<theDPVector[i]->second<<G4endl;
#endif
    //EnMomConservation= theDPVector[i]->first->GetEnegryMomentumConservation();
    //nOfNeutrons      = theDPVector[i]->first->GetNumberOfNeutronsInTarget();
    return theDPVector[i]->first->PostStepDoIt(Track, Step);
  }
  return G4VDiscreteProcess::PostStepDoIt(Track, Step);
}

G4double G4QDiscProcessMixer::PostStepGetPhysicalInteractionLength	(	const G4Track &	track,
		G4double	previousStepSize,
		G4ForceCondition *	condition
	)			[virtual]

Reimplemented from G4VDiscreteProcess.

Definition at line 106 of file G4QDiscProcessMixer.cc.

References DBL_MAX, G4cout, G4endl, G4Track::GetDynamicParticle(), and G4DynamicParticle::GetKineticEnergy().

{
  G4double kEn=Track.GetDynamicParticle()->GetKineticEnergy(); // Projectile kinetic energy
  G4int maxDP=theDPVector.size();
  if(maxDP) for(G4int i=maxDP-1; i>-1; i--) if(kEn < theDPVector[i]->second)
  {
#ifdef debug
    G4cout<<"G4QDPMix::PSGetPIL:"<<i<<",kEn="<<kEn<<" < "<<theDPVector[i]->second<<", PIL="
        <<theDPVector[i]->first->PostStepGetPhysicalInteractionLength(Track,PrevStSize,F)
        <<G4endl;
#endif
    return theDPVector[i]->first->PostStepGetPhysicalInteractionLength(Track,PrevStSize,F);
  }
  return DBL_MAX;
}

---

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

• G4QDiscProcessMixer.hh
• G4QDiscProcessMixer.cc

---