#include <G4ParticleHPProduct.hh>

Data Structures
struct	toBeCached

Public Member Functions
	G4ParticleHPProduct ()

G4double	GetMass ()

G4double	GetMassCode ()

G4double	GetMeanYield (G4double anEnergy)

G4int	GetMultiplicity (G4double anEnergy)

G4ReactionProduct *	GetProjectileRP ()

G4double	GetQValue ()

G4ReactionProduct *	GetTarget ()

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

G4double	MeanEnergyOfThisInteraction ()

G4ReactionProductVector *	Sample (G4double anEnergy, G4int nParticles)

void	SetProjectileRP (G4ReactionProduct *aIncidentPart)

void	SetTarget (G4ReactionProduct *aTarget)

	~G4ParticleHPProduct ()

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

References fCache, G4HPMultiBetweenInts, G4HPMultiPoisson, theActualStateQValue, theDist, theDistLaw, theGroundStateQValue, theIsomerFlag, theMass, theMassCode, and theMultiplicityMethod.

◆ ~G4ParticleHPProduct()

G4ParticleHPProduct::~G4ParticleHPProduct ( )

inline

Definition at line 94 of file G4ParticleHPProduct.hh.

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

References theDist.

Member Function Documentation

◆ GetMass()

G4double G4ParticleHPProduct::GetMass ( )

inline

Definition at line 222 of file G4ParticleHPProduct.hh.

222{return theMass;}

References theMass.

◆ GetMassCode()

G4double G4ParticleHPProduct::GetMassCode ( )

inline

Definition at line 221 of file G4ParticleHPProduct.hh.

221{return theMassCode;}

References theMassCode.

◆ GetMeanYield()

G4double G4ParticleHPProduct::GetMeanYield ( G4double anEnergy )

inline

Definition at line 174 of file G4ParticleHPProduct.hh.

  {
    return theYield.GetY(anEnergy);
  }

References G4ParticleHPVector::GetY(), and theYield.

◆ GetMultiplicity()

G4int G4ParticleHPProduct::GetMultiplicity ( G4double anEnergy )

Definition at line 45 of file G4ParticleHPProduct.cc.

{
  //if(theDist == 0) { return 0; }
  //151120 TK Modified for solving reproducibility problem 
  if ( theDist == 0 ) { 
     fCache.Get().theCurrentMultiplicity = 0;
     return 0; 
  }
 
  G4double mean = theYield.GetY(anEnergy);
  //g  G4cout << "G4ParticleHPProduct MEAN NUMBER OF PARTICLES " << mean << " for " << theMass << G4endl;
  //if( mean <= 0. ) return 0;
  //151120 TK Modified for solving reproducibility problem 
  //This is also a real fix
  if ( mean <= 0. )  {
     fCache.Get().theCurrentMultiplicity = 0;
     return 0; 
  }
 
  G4int multi;
  multi = G4int(mean+0.0001);
  //if(theMassCode==0) multi = G4Poisson(mean); // @@@@gammas. please X-check this
  //080718
#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 = G4Poisson ( mean );
     #ifdef G4VERBOSE
     if( std::getenv("G4PHPTEST") && G4HadronicParameters::Instance()->GetVerboseLevel() > 0 )
       G4cout << " MULTIPLICITY MULTIPLIED " << multi << " " << theMassCode << G4endl;
     #endif
       } else { // if( theMultiplicityMethod == G4HPMultiBetweenInts ) {
     G4double radnf = CLHEP::RandFlat::shoot();
     G4int imulti = G4int(mean);
     multi = imulti + G4int(radnf < mean-imulti);
     //  G4cout << theMass << " multi " << multi << " mean " << mean 
     //     << " radnf " << radnf << " mean-imulti " << mean-imulti << G4endl;
       }
#endif
    //       multi = int(mean);
    //       if( CLHEP::RandFlat::shoot() > mean-multi ) multi++;
     }
#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;
}

References fCache, G4cout, G4endl, G4HPMultiPoisson, G4Poisson(), G4ParticleHPVector::GetY(), G4HadronicParameters::Instance(), CLHEP::RandFlat::shoot(), theDist, theMass, theMassCode, theMultiplicityMethod, and theYield.

Referenced by G4ParticleHPEnAngCorrelation::Sample().

◆ GetProjectileRP()

G4ReactionProduct * G4ParticleHPProduct::GetProjectileRP ( )

inline

Definition at line 194 of file G4ParticleHPProduct.hh.

  {
    return fCache.Get().theProjectileRP;
  }

References fCache.

◆ GetQValue()

G4double G4ParticleHPProduct::GetQValue ( )

inline

Definition at line 214 of file G4ParticleHPProduct.hh.

  {
    return theActualStateQValue;
  }

References theActualStateQValue.

Referenced by G4ParticleHPEnAngCorrelation::Sample().

◆ GetTarget()

G4ReactionProduct * G4ParticleHPProduct::GetTarget ( void )

inline

Definition at line 189 of file G4ParticleHPProduct.hh.

  {
    return fCache.Get().theTarget;
  }

References fCache.

◆ 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;
    if( std::getenv("G4PHPTEST") )
      G4cout << " G4ParticleHPProduct :: Init MassCode "
             << theMassCode << " " << theMass << " theActualStateQValue "
             << theActualStateQValue << G4endl;// GDEB
    if( std::getenv("G4PHPTEST") )
      G4cout << " G4ParticleHPProduct :: Init theActualStateQValue "
             << theActualStateQValue << G4endl;// GDEB
    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);
    }
  }

References CLHEP::eV, G4cout, G4endl, G4ParticleHPVector::Hash(), G4VParticleHPEnergyAngular::Init(), G4ParticleHPVector::Init(), G4VParticleHPEnergyAngular::SetQValue(), theActualStateQValue, theDist, theDistLaw, theGroundStateQValue, theIsomerFlag, theMass, theMassCode, and theYield.

Referenced by G4ParticleHPEnAngCorrelation::Init().

◆ MeanEnergyOfThisInteraction()

G4double G4ParticleHPProduct::MeanEnergyOfThisInteraction ( )

inline

Definition at line 199 of file G4ParticleHPProduct.hh.

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

References fCache, G4VParticleHPEnergyAngular::MeanEnergyOfThisInteraction(), and theDist.

Referenced by G4ParticleHPEnAngCorrelation::Sample().

◆ Sample()

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

Definition at line 112 of file G4ParticleHPProduct.cc.

{
  if(theDist == 0) { return 0; }
  G4ReactionProductVector * result = new G4ReactionProductVector;
 
  theDist->SetTarget(fCache.Get().theTarget);
  theDist->SetProjectileRP(fCache.Get().theProjectileRP);
  G4int i;
//  G4double eMax = GetTarget()->GetMass()+GetNeutron()->GetMass()
//                  - theActualStateQValue;
  G4ReactionProduct * tmp;
  theDist->ClearHistories();
 
  for(i=0;i<multi;i++)
  {
#ifdef G4PHPDEBUG
 if( std::getenv("G4PHPTEST") )
    #ifdef G4VERBOSE
    if( std::getenv("G4ParticleHPDebug") && tmp != 0 && G4HadronicParameters::Instance()->GetVerboseLevel() > 0 )
      G4cout << multi << " " << i << " @@@ G4ParticleHPProduct::Sample " << anEnergy << " Mass " << theMassCode << " " << theMass << G4endl;
    #endif
#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;
    }
  /*
  //080901 TK Comment out, too many secondaries are produced in deuteron reactions
  if(theTarget->GetMass()<2*GeV) // @@@ take care of residuals in all cases
  {
  tmp = theDist->Sample(anEnergy, theMassCode, theMass);
  tmp->SetDefinition(G4Proton::Proton());
  if(tmp != 0) { result->push_back(tmp); }
  }
  */
 
  return result;
}

References G4VParticleHPEnergyAngular::ClearHistories(), fCache, G4cout, G4endl, G4ReactionProduct::GetDefinition(), G4ReactionProduct::GetKineticEnergy(), G4ParticleDefinition::GetParticleName(), G4HadronicParameters::Instance(), G4VParticleHPEnergyAngular::Sample(), G4VParticleHPEnergyAngular::SetProjectileRP(), G4VParticleHPEnergyAngular::SetTarget(), theDist, theMass, and theMassCode.

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

◆ SetProjectileRP()

void G4ParticleHPProduct::SetProjectileRP ( G4ReactionProduct * aIncidentPart )

inline

Definition at line 179 of file G4ParticleHPProduct.hh.

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

References fCache.

◆ SetTarget()

void G4ParticleHPProduct::SetTarget ( G4ReactionProduct * aTarget )

inline

Definition at line 184 of file G4ParticleHPProduct.hh.

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

References fCache.

Field Documentation

◆ fCache

G4Cache<toBeCached> G4ParticleHPProduct::fCache

private

Definition at line 239 of file G4ParticleHPProduct.hh.

Referenced by G4ParticleHPProduct(), GetMultiplicity(), GetProjectileRP(), GetTarget(), MeanEnergyOfThisInteraction(), Sample(), SetProjectileRP(), and SetTarget().

◆ theActualStateQValue

G4double G4ParticleHPProduct::theActualStateQValue

private

Definition at line 232 of file G4ParticleHPProduct.hh.

Referenced by G4ParticleHPProduct(), GetQValue(), and Init().

◆ theDist

G4VParticleHPEnergyAngular* G4ParticleHPProduct::theDist

private

Definition at line 235 of file G4ParticleHPProduct.hh.

Referenced by G4ParticleHPProduct(), GetMultiplicity(), Init(), MeanEnergyOfThisInteraction(), Sample(), and ~G4ParticleHPProduct().

◆ theDistLaw

G4int G4ParticleHPProduct::theDistLaw

private

Definition at line 233 of file G4ParticleHPProduct.hh.

Referenced by G4ParticleHPProduct(), and Init().

◆ theGroundStateQValue

G4double G4ParticleHPProduct::theGroundStateQValue

private

Definition at line 231 of file G4ParticleHPProduct.hh.

Referenced by G4ParticleHPProduct(), and Init().

◆ theIsomerFlag

G4int G4ParticleHPProduct::theIsomerFlag

private

Definition at line 230 of file G4ParticleHPProduct.hh.

Referenced by G4ParticleHPProduct(), and Init().

◆ theMass

G4double G4ParticleHPProduct::theMass

private

Definition at line 229 of file G4ParticleHPProduct.hh.

Referenced by G4ParticleHPProduct(), GetMass(), GetMultiplicity(), Init(), and Sample().

◆ theMassCode

G4double G4ParticleHPProduct::theMassCode

private

Definition at line 228 of file G4ParticleHPProduct.hh.

Referenced by G4ParticleHPProduct(), GetMassCode(), GetMultiplicity(), Init(), and Sample().

◆ theMultiplicityMethod

G4HPMultiMethod G4ParticleHPProduct::theMultiplicityMethod

private

Definition at line 241 of file G4ParticleHPProduct.hh.

Referenced by G4ParticleHPProduct(), and GetMultiplicity().

◆ theYield

G4ParticleHPVector G4ParticleHPProduct::theYield

private

Definition at line 234 of file G4ParticleHPProduct.hh.

Referenced by GetMeanYield(), GetMultiplicity(), and Init().

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

source/processes/hadronic/models/particle_hp/include/G4ParticleHPProduct.hh
source/processes/hadronic/models/particle_hp/src/G4ParticleHPProduct.cc

Private Attributes
G4Cache< toBeCached >	fCache

G4double	theActualStateQValue

G4VParticleHPEnergyAngular *	theDist

G4int	theDistLaw

G4double	theGroundStateQValue

G4int	theIsomerFlag

G4double	theMass

G4double	theMassCode

G4HPMultiMethod	theMultiplicityMethod

G4ParticleHPVector	theYield

Data Structures

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4ParticleHPProduct()

◆ ~G4ParticleHPProduct()

Member Function Documentation

◆ GetMass()

◆ GetMassCode()

◆ GetMeanYield()

◆ GetMultiplicity()

◆ GetProjectileRP()

◆ GetQValue()

◆ GetTarget()

◆ Init()

◆ MeanEnergyOfThisInteraction()

◆ Sample()

◆ SetProjectileRP()

◆ SetTarget()

Field Documentation

◆ fCache

◆ theActualStateQValue

◆ theDist

◆ theDistLaw

◆ theGroundStateQValue

◆ theIsomerFlag

◆ theMass

◆ theMassCode

◆ theMultiplicityMethod

◆ theYield