G4ParticleHPProduct Class Reference

#include <G4ParticleHPProduct.hh>

Public Member Functions
	G4ParticleHPProduct ()
	~G4ParticleHPProduct ()
void	Init (std::istream &aDataFile, const G4ParticleDefinition *projectile)
G4int	GetMultiplicity (G4double anEnergy)
G4ReactionProductVector *	Sample (G4double anEnergy, G4int nParticles)
G4double	GetMeanYield (G4double anEnergy)
void	SetProjectileRP (G4ReactionProduct *aIncidentPart)
void	SetTarget (G4ReactionProduct *aTarget)
G4ReactionProduct *	GetTarget ()
G4ReactionProduct *	GetProjectileRP ()
G4double	MeanEnergyOfThisInteraction ()
G4double	GetQValue ()
G4double	GetMassCode ()
G4double	GetMass ()

Detailed Description

Definition at line 50 of file G4ParticleHPProduct.hh.

Constructor & Destructor Documentation

G4ParticleHPProduct()

G4ParticleHPProduct::G4ParticleHPProduct

(

)

Definition at line 54 of file G4ParticleHPProduct.cc.

{
  toBeCached val;
  fCache.Put(val);
 
  if (G4ParticleHPManager::GetInstance()->GetPHCUsePoisson()) {
    theMultiplicityMethod = G4HPMultiPoisson;
  } else {
    theMultiplicityMethod = G4HPMultiBetweenInts;
  }
}

~G4ParticleHPProduct()

G4ParticleHPProduct::~G4ParticleHPProduct

(

)

Definition at line 66 of file G4ParticleHPProduct.cc.

{
  delete theDist;
}

Member Function Documentation

GetMass()

G4double G4ParticleHPProduct::GetMass ( )

inline

Definition at line 97 of file G4ParticleHPProduct.hh.

{ return theMass; }

GetMassCode()

G4double G4ParticleHPProduct::GetMassCode ( )

inline

Definition at line 96 of file G4ParticleHPProduct.hh.

{ return theMassCode; }

GetMeanYield()

G4double G4ParticleHPProduct::GetMeanYield ( G4double anEnergy )

inline

Definition at line 69 of file G4ParticleHPProduct.hh.

{ return theYield.GetY(anEnergy); }

GetMultiplicity()

G4int G4ParticleHPProduct::GetMultiplicity

(

G4double

anEnergy

)

Definition at line 124 of file G4ParticleHPProduct.cc.

{
  if (theDist == nullptr) {
    fCache.Get().theCurrentMultiplicity = 0;
    return 0;
  }
 
  G4double mean = theYield.GetY(anEnergy);
  if (mean <= 0.) {
    fCache.Get().theCurrentMultiplicity = 0;
    return 0;
  }
  G4int multi = (theMultiplicityMethod == G4HPMultiPoisson) ?
    (G4int)G4Poisson(mean) : G4lrint(mean);
 
#ifdef G4VERBOSE
  if (G4ParticleHPManager::GetInstance()->GetDEBUG())
    G4cout << "G4ParticleHPProduct::GetMultiplicity code=" << theMassCode << " M=" << theMass
       << " multi=" << multi << " mean=" << mean << G4endl;
#endif
  fCache.Get().theCurrentMultiplicity = multi;
 
  return multi;
}

Referenced by G4ParticleHPEnAngCorrelation::Sample().

GetProjectileRP()

G4ReactionProduct * G4ParticleHPProduct::GetProjectileRP ( )

inline

Definition at line 80 of file G4ParticleHPProduct.hh.

{ return fCache.Get().theProjectileRP; }

GetQValue()

G4double G4ParticleHPProduct::GetQValue ( )

inline

Definition at line 92 of file G4ParticleHPProduct.hh.

{ return theActualStateQValue; }

Referenced by G4ParticleHPEnAngCorrelation::Sample().

GetTarget()

G4ReactionProduct * G4ParticleHPProduct::GetTarget ( )

inline

Definition at line 78 of file G4ParticleHPProduct.hh.

{ return fCache.Get().theTarget; }

Init()

void G4ParticleHPProduct::Init	(	std::istream &	aDataFile,
		const G4ParticleDefinition *	projectile )

Definition at line 71 of file G4ParticleHPProduct.cc.

{
  aDataFile >> theMassCode >> theMass >> theIsomerFlag >> theDistLaw >> theGroundStateQValue
        >> theActualStateQValue;
  theGroundStateQValue *= CLHEP::eV;
  theActualStateQValue *= CLHEP::eV;
  theYield.Init(aDataFile, CLHEP::eV);
  theYield.Hash();
  if (theDistLaw == 0) {
    // distribution not known, use E-independent, isotropic
    // angular distribution
    theDist = new G4ParticleHPIsotropic;
  }
  else if (theDistLaw == 1) {
    // Continuum energy-angular distribution
    theDist = new G4ParticleHPContEnergyAngular(projectile);
  }
  else if (theDistLaw == 2) {
    // Discrete 2-body scattering
    theDist = new G4ParticleHPDiscreteTwoBody;
  }
  else if (theDistLaw == 3) {
    // Isotropic emission
    theDist = new G4ParticleHPIsotropic;
  }
  else if (theDistLaw == 4) {
    // Discrete 2-body recoil modification not used for now.
    // theDist = new G4ParticleHPDiscreteTwoBody;
    // the above is only temporary;
    // recoils need to be addressed properly
  }
  //    else if(theDistLaw == 5)
  //    {
  // charged particles only, to be used in a later stage. @@@@
  //    }
  else if (theDistLaw == 6) {
    // N-Body phase space
    theDist = new G4ParticleHPNBodyPhaseSpace;
  }
  else if (theDistLaw == 7) {
    // Laboratory angular energy paraetrisation
    theDist = new G4ParticleHPLabAngularEnergy;
  }
  else {
        throw G4HadronicException(__FILE__, __LINE__,
                                  "distribution law unknown to G4ParticleHPProduct");
  }
  if (theDist != nullptr) {
    theDist->SetQValue(theActualStateQValue);
    theDist->Init(aDataFile);
  }
}

Referenced by G4ParticleHPEnAngCorrelation::Init().

MeanEnergyOfThisInteraction()

G4double G4ParticleHPProduct::MeanEnergyOfThisInteraction ( )

inline

Definition at line 82 of file G4ParticleHPProduct.hh.

    {
      G4double result = 0.0;
      if (theDist != nullptr) {
        result = theDist->MeanEnergyOfThisInteraction();
        result *= fCache.Get().theCurrentMultiplicity;
      }
      return result;
    }

Referenced by G4ParticleHPEnAngCorrelation::Sample().

Sample()

G4ReactionProductVector * G4ParticleHPProduct::Sample	(	G4double	anEnergy,
		G4int	nParticles )

Definition at line 149 of file G4ParticleHPProduct.cc.

{
  if (theDist == nullptr) {
    return nullptr;
  }
  auto result = new G4ReactionProductVector;
 
  theDist->SetTarget(fCache.Get().theTarget);
  theDist->SetProjectileRP(fCache.Get().theProjectileRP);
  G4ReactionProduct* tmp;
  theDist->ClearHistories();
 
  for (G4int i = 0; i < multi; ++i) {
    tmp = theDist->Sample(anEnergy, theMassCode, theMass);
    if (tmp != nullptr) {
      result->push_back(tmp);
#ifdef G4VERBOSE
      if (G4ParticleHPManager::GetInstance()->GetDEBUG())
    G4cout << "multi=" << multi << " i=" << i << " G4ParticleHPProduct::Sample "
           << tmp->GetDefinition()->GetParticleName() << " E=" << tmp->GetKineticEnergy()
           << G4endl;
#endif
    }
  }
  if (multi == 0) {
    tmp = theDist->Sample(anEnergy, theMassCode, theMass);
    delete tmp;
  }
  return result;
}

Referenced by G4ParticleHPEnAngCorrelation::Sample(), and G4ParticleHPEnAngCorrelation::SampleOne().

SetProjectileRP()

void G4ParticleHPProduct::SetProjectileRP ( G4ReactionProduct * aIncidentPart )

inline

Definition at line 71 of file G4ParticleHPProduct.hh.

    {
      fCache.Get().theProjectileRP = aIncidentPart;
    }

SetTarget()

void G4ParticleHPProduct::SetTarget ( G4ReactionProduct * aTarget )

inline

Definition at line 76 of file G4ParticleHPProduct.hh.

{ fCache.Get().theTarget = aTarget; }

---

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

• G4ParticleHPProduct.hh
• G4ParticleHPProduct.cc