G4ParticleHPProduct Class Reference

#include <G4ParticleHPProduct.hh>

Public Member Functions
	G4ParticleHPProduct ()
	~G4ParticleHPProduct ()
void	Init (std::istream &aDataFile, 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 52 of file G4ParticleHPProduct.hh.

Constructor & Destructor Documentation

G4ParticleHPProduct()

G4ParticleHPProduct::G4ParticleHPProduct ( )

inline

Definition at line 65 of file G4ParticleHPProduct.hh.

  {
    theDist = 0;
    toBeCached val;
    fCache.Put( val );
 
    char * method = std::getenv( "G4PHP_MULTIPLICITY_METHOD" );
    if( method )
    {
      if( G4String(method) == "Poisson" ) {
    theMultiplicityMethod = G4HPMultiPoisson;
      } else if( G4String(method) == "BetweenInts" ) {
    theMultiplicityMethod = G4HPMultiBetweenInts;
      } else {
    throw G4HadronicException(__FILE__, __LINE__, ("multiplicity method unknown to G4ParticleHPProduct" + G4String(method)).c_str());
      }
    }
    else
    {
      theMultiplicityMethod = G4HPMultiPoisson;
    }
    theMassCode = 0.0;
    theMass = 0.0;
    theIsomerFlag = 0;
    theGroundStateQValue = 0.0;
    theActualStateQValue = 0.0;
    theDistLaw = -1;
  }

~G4ParticleHPProduct()

G4ParticleHPProduct::~G4ParticleHPProduct ( )

inline

Definition at line 94 of file G4ParticleHPProduct.hh.

  { 
    if(theDist != 0) delete theDist;
  }

Member Function Documentation

GetMass()

G4double G4ParticleHPProduct::GetMass ( )

inline

Definition at line 215 of file G4ParticleHPProduct.hh.

{return theMass;}

GetMassCode()

G4double G4ParticleHPProduct::GetMassCode ( )

inline

Definition at line 214 of file G4ParticleHPProduct.hh.

{return theMassCode;}

GetMeanYield()

G4double G4ParticleHPProduct::GetMeanYield ( G4double  anEnergy )

inline

Definition at line 167 of file G4ParticleHPProduct.hh.

  {
    return theYield.GetY(anEnergy);
  }

GetMultiplicity()

G4int G4ParticleHPProduct::GetMultiplicity

(

G4double

anEnergy

)

Definition at line 45 of file G4ParticleHPProduct.cc.

{
  if ( theDist == 0 ) { 
     fCache.Get().theCurrentMultiplicity = 0;
     return 0; 
  }
 
  G4double mean = theYield.GetY(anEnergy);
  if ( mean <= 0. )  {
     fCache.Get().theCurrentMultiplicity = 0;
     return 0; 
  }
 
  G4int multi;
  multi = G4int(mean+0.0001);
#ifdef PHP_AS_HP
  if ( theMassCode == 0 ) // DELETE THIS: IT MUST BE DONE FOR ALL PARTICLES
#endif
  { 
     if ( G4int ( mean ) == mean )
     {
        multi = (G4int) mean;
     }
     else
     {
#ifdef PHP_AS_HP
     multi = G4Poisson ( mean ); 
#else 
       if( theMultiplicityMethod == G4HPMultiPoisson )
       {
     multi = (G4int)G4Poisson ( mean );
       }
       else
       {
     G4double radnf = CLHEP::RandFlat::shoot();
     G4int imulti = G4int(mean);
     multi = imulti + G4int(radnf < mean-imulti);
       }
#endif
     }
#ifdef G4PHPDEBUG
     #ifdef G4VERBOSE
     if( std::getenv("G4ParticleHPDebug") && G4HadronicParameters::Instance()->GetVerboseLevel() > 0 )
       G4cout << "G4ParticleHPProduct::GetMultiplicity " << theMassCode
              << " " << theMass << " multi " << multi << " mean " << mean
              << G4endl;
     #endif
#endif
  }
 
  fCache.Get().theCurrentMultiplicity = static_cast<G4int>(mean);
 
  return multi;
}

Referenced by G4ParticleHPEnAngCorrelation::Sample().

GetProjectileRP()

G4ReactionProduct * G4ParticleHPProduct::GetProjectileRP ( )

inline

Definition at line 187 of file G4ParticleHPProduct.hh.

  {
    return fCache.Get().theProjectileRP;
  }

GetQValue()

G4double G4ParticleHPProduct::GetQValue ( )

inline

Definition at line 207 of file G4ParticleHPProduct.hh.

  {
    return theActualStateQValue;
  }

Referenced by G4ParticleHPEnAngCorrelation::Sample().

GetTarget()

G4ReactionProduct * G4ParticleHPProduct::GetTarget ( )

inline

Definition at line 182 of file G4ParticleHPProduct.hh.

  {
    return fCache.Get().theTarget;
  }

Init()

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

inline

Definition at line 99 of file G4ParticleHPProduct.hh.

  {
    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
      delete theDist;
      theDist = 0;
    }
    //    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!=0)
    {
      theDist->SetQValue(theActualStateQValue);      
      theDist->Init(aDataFile);
    }
  }

Referenced by G4ParticleHPEnAngCorrelation::Init().

MeanEnergyOfThisInteraction()

G4double G4ParticleHPProduct::MeanEnergyOfThisInteraction ( )

inline

Definition at line 192 of file G4ParticleHPProduct.hh.

  { 
    G4double result;
    if(theDist == 0)
    {
      result = 0;
    }
    else
    {
      result=theDist->MeanEnergyOfThisInteraction();
      result *= fCache.Get().theCurrentMultiplicity;
    }
    return result;
  }

Referenced by G4ParticleHPEnAngCorrelation::Sample().

Sample()

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

Definition at line 101 of file G4ParticleHPProduct.cc.

{
  if(theDist == 0) { return nullptr; }
  G4ReactionProductVector * result = new G4ReactionProductVector;
 
  theDist->SetTarget(fCache.Get().theTarget);
  theDist->SetProjectileRP(fCache.Get().theProjectileRP);
  G4int i;
  G4ReactionProduct * tmp;
  theDist->ClearHistories();
 
  for(i=0; i<multi; ++i)
  {
#ifdef G4PHPDEBUG
 
#endif
    tmp = theDist->Sample(anEnergy, theMassCode, theMass);
    if(tmp != 0) { result->push_back(tmp); }
#ifndef G4PHPDEBUG //GDEB
    #ifdef G4VERBOSE
    if( std::getenv("G4ParticleHPDebug")
     && tmp != 0 && G4HadronicParameters::Instance()->GetVerboseLevel() > 0 )
      G4cout << multi << " " << i << " @@@ G4ParticleHPProduct::Sample "
             << tmp->GetDefinition()->GetParticleName() << " E= "
             << tmp->GetKineticEnergy() << G4endl;
    #endif
#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 172 of file G4ParticleHPProduct.hh.

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

SetTarget()

void G4ParticleHPProduct::SetTarget ( G4ReactionProduct *  aTarget )

inline

Definition at line 177 of file G4ParticleHPProduct.hh.

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

---

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

• G4ParticleHPProduct.hh
• G4ParticleHPProduct.cc