G4ParticleHPEnAngCorrelation Class Reference

#include <G4ParticleHPEnAngCorrelation.hh>

Data Structures
struct	toBeCached

Public Member Functions
	G4ParticleHPEnAngCorrelation ()
	G4ParticleHPEnAngCorrelation (G4ParticleDefinition *proj)
G4double	GetTargetMass ()
G4double	GetTotalMeanEnergy ()
G4bool	InCharge ()
void	Init (std::istream &aDataFile)
G4ReactionProductVector *	Sample (G4double anEnergy)
G4ReactionProduct *	SampleOne (G4double anEnergy)
void	SetProjectileRP (G4ReactionProduct &aIncidentPart)
void	SetTarget (G4ReactionProduct &aTarget)
	~G4ParticleHPEnAngCorrelation ()

Private Attributes
G4Cache< toBeCached >	fCache
G4int	frameFlag
G4bool	inCharge
G4int	nProducts
G4double	targetMass
G4ParticleHPProduct *	theProducts
G4ParticleDefinition *	theProjectile

Detailed Description

Definition at line 46 of file G4ParticleHPEnAngCorrelation.hh.

Constructor & Destructor Documentation

G4ParticleHPEnAngCorrelation() [1/2]

G4ParticleHPEnAngCorrelation::G4ParticleHPEnAngCorrelation ( )

inline

Definition at line 59 of file G4ParticleHPEnAngCorrelation.hh.

  {
    theProjectile = G4Neutron::Neutron();
    theProducts = 0;
    inCharge = false;
    toBeCached val;
    fCache.Put( val );
    //theTotalMeanEnergy = -1.;
    fCache.Get().theTotalMeanEnergy = -1.;
    targetMass = 0.0;
    frameFlag = 0;
    nProducts = 0;
  }

References fCache, frameFlag, inCharge, G4Neutron::Neutron(), nProducts, targetMass, theProducts, and theProjectile.

G4ParticleHPEnAngCorrelation() [2/2]

G4ParticleHPEnAngCorrelation::G4ParticleHPEnAngCorrelation ( G4ParticleDefinition *  proj )

inline

Definition at line 73 of file G4ParticleHPEnAngCorrelation.hh.

    : theProjectile(proj)
  {
    theProducts = 0;
    inCharge = false;
    toBeCached val;
    fCache.Put( val );
    //theTotalMeanEnergy = -1.;
    fCache.Get().theTotalMeanEnergy = -1.;
    targetMass = 0.0;
    frameFlag = 0;
    nProducts = 0;
  }

References fCache, frameFlag, inCharge, nProducts, targetMass, and theProducts.

~G4ParticleHPEnAngCorrelation()

G4ParticleHPEnAngCorrelation::~G4ParticleHPEnAngCorrelation ( )

inline

Definition at line 87 of file G4ParticleHPEnAngCorrelation.hh.

  {
    if(theProducts!=0) delete [] theProducts;
  }

References theProducts.

Member Function Documentation

GetTargetMass()

G4double G4ParticleHPEnAngCorrelation::GetTargetMass ( )

inline

Definition at line 126 of file G4ParticleHPEnAngCorrelation.hh.

  {
    return targetMass;
  }

References targetMass.

Referenced by G4ParticleHPInelasticBaseFS::BaseApply().

GetTotalMeanEnergy()

G4double G4ParticleHPEnAngCorrelation::GetTotalMeanEnergy ( )

inline

Definition at line 131 of file G4ParticleHPEnAngCorrelation.hh.

  {
    return fCache.Get().theTotalMeanEnergy;  // cached in 'sample' call
  }

References fCache.

Referenced by G4ParticleHPInelasticBaseFS::BaseApply().

InCharge()

G4bool G4ParticleHPEnAngCorrelation::InCharge ( )

inline

Definition at line 121 of file G4ParticleHPEnAngCorrelation.hh.

  {
    return inCharge;
  }

References inCharge.

Init()

void G4ParticleHPEnAngCorrelation::Init ( std::istream &  aDataFile )

inline

Definition at line 92 of file G4ParticleHPEnAngCorrelation.hh.

  {
    inCharge = true;
    aDataFile>>targetMass>>frameFlag>>nProducts;
    theProducts = new G4ParticleHPProduct[nProducts];
    for(G4int i=0; i<nProducts; i++)
    {
      theProducts[i].Init(aDataFile,theProjectile);
    }
  }

References frameFlag, inCharge, G4ParticleHPProduct::Init(), nProducts, targetMass, theProducts, and theProjectile.

Referenced by G4ParticleHPCaptureFS::Init(), G4ParticleHPInelasticCompFS::Init(), and G4ParticleHPInelasticBaseFS::Init().

Sample()

G4ReactionProductVector * G4ParticleHPEnAngCorrelation::Sample

(

G4double

anEnergy

)

Definition at line 77 of file G4ParticleHPEnAngCorrelation.cc.

{
  G4ReactionProductVector * result = new G4ReactionProductVector;
  G4int i;
  G4ReactionProductVector * it;
  G4ReactionProduct theCMS;
  G4LorentzRotation toZ;
 
  if(frameFlag==2 
     || frameFlag==3) // Added for particle HP
  {
    // simplify and double check @
    G4ThreeVector the3IncidentPart = fCache.Get().theProjectileRP->GetMomentum(); //theProjectileRP has value in LAB
    G4double nEnergy = fCache.Get().theProjectileRP->GetTotalEnergy();
    G4ThreeVector the3Target = fCache.Get().theTarget->GetMomentum();  //theTarget has value in LAB
    G4double tEnergy = fCache.Get().theTarget->GetTotalEnergy();
    G4double totE = nEnergy+tEnergy;
    G4ThreeVector the3CMS = the3Target+the3IncidentPart;
    theCMS.SetMomentum(the3CMS);
    G4double cmsMom = std::sqrt(the3CMS*the3CMS);
    G4double sqrts = std::sqrt((totE-cmsMom)*(totE+cmsMom));
    theCMS.SetMass(sqrts);
    theCMS.SetTotalEnergy(totE);
    G4ReactionProduct aIncidentPart;
    aIncidentPart.Lorentz(*fCache.Get().theProjectileRP, theCMS);
    //TKDB 100413 
    //ENDF-6 Formats Manual ENDF-102
    //CHAPTER 6. FILE 6: PRODUCT ENERGY-ANGLE DISTRIBUTIONS
    //LCT Reference system for secondary energy and angle (incident energy is always given in the LAB system)
    //anEnergy = aIncidentPart.GetKineticEnergy();
    anEnergy = fCache.Get().theProjectileRP->GetKineticEnergy(); //should be same argumment of "anEnergy"
 
    G4LorentzVector Ptmp (aIncidentPart.GetMomentum(), aIncidentPart.GetTotalEnergy());
 
    toZ.rotateZ(-1*Ptmp.phi());
    toZ.rotateY(-1*Ptmp.theta());
  }
  fCache.Get().theTotalMeanEnergy=0;
  G4LorentzRotation toLab(toZ.inverse()); //toLab only change axis NOT to LAB system
 
  for(i=0; i<nProducts; i++)
  {
    G4int nPart = theProducts[i].GetMultiplicity(anEnergy);
    //-    if( nParticles[i] == 0 ) continue;
    it = theProducts[i].Sample(anEnergy,nPart); 
    G4double aMeanEnergy = theProducts[i].MeanEnergyOfThisInteraction();
    //    if( getenv("G4PHPTEST") ) G4cout << " EnAnG energy sampled " << it->operator[](0)->GetKineticEnergy() << " aMeanEnergy " << aMeanEnergy << G4endl; // GDEB
    //if(aMeanEnergy>0)
    //151120 TK Modified for solving reproducibility problem 
    //This change may have side effect. 
    if(aMeanEnergy>=0)
    {
      fCache.Get().theTotalMeanEnergy += aMeanEnergy;
    }
    else
    {
      fCache.Get().theTotalMeanEnergy = anEnergy/nProducts+theProducts[i].GetQValue();
    }
    if(it!=0)
    {
      for(unsigned int ii=0; ii<it->size(); ii++)
      {
    //if(!std::getenv("G4PHP_NO_LORENTZ_BOOST")) {
    G4LorentzVector pTmp1 (it->operator[](ii)->GetMomentum(),
                   it->operator[](ii)->GetTotalEnergy());
    pTmp1 = toLab*pTmp1;
    if( std::getenv("G4PHPTEST") )  G4cout << " G4particleHPEnAngCorrelation COS THETA " <<  std::cos(it->operator[](ii)->GetMomentum().theta()) << G4endl;
    it->operator[](ii)->SetMomentum(pTmp1.vect());
    it->operator[](ii)->SetTotalEnergy(pTmp1.e());
    if( std::getenv("G4PHPTEST") ) G4cout << " G4particleHPEnAngCorrelation COS THETA after toLab " <<  std::cos(it->operator[](ii)->GetMomentum().theta()) << G4endl;
 
    if(frameFlag==1) // target rest //TK 100413 should be LAB?
    {
      it->operator[](ii)->Lorentz(*(it->operator[](ii)), -1.*(*fCache.Get().theTarget)); //TK 100413 Is this really need?
    }
    else if(frameFlag==2 ) // CMS
        {
#ifdef G4PHPDEBUG
      if( std::getenv("G4ParticleHPDebug") ) 
        G4cout <<"G4ParticleHPEnAngCorrelation: before Lorentz boost "<<
          it->at(ii)->GetKineticEnergy()<<" "<<
          it->at(ii)->GetMomentum()<<G4endl;
#endif
      it->operator[](ii)->Lorentz(*(it->operator[](ii)), -1.*theCMS);
#ifdef G4PHPDEBUG
      if( std::getenv("G4ParticleHPDebug") ) 
        G4cout <<"G4ParticleHPEnAngCorrelation: after Lorentz boost "<<
        it->at(ii)->GetKineticEnergy()<<" "<<
          it->at(ii)->GetMomentum()<<G4endl;
#endif
    }
        //TK120515 migrate frameFlag (MF6 LCT) = 3 
    else if(frameFlag==3) // CMS A<=4 other LAB
        {
           if ( theProducts[i].GetMassCode() > 2004.5 ) //Alpha AWP 3.96713
           {
              //LAB
              it->operator[](ii)->Lorentz(*(it->operator[](ii)), -1.*(*fCache.Get().theTarget)); //TK 100413 Is this really need?
#ifdef G4PHPDEBUG
      if( std::getenv("G4ParticleHPDebug") ) 
        G4cout <<"G4ParticleHPEnAngCorrelation: after Lorentz boost "<<
        it->at(ii)->GetKineticEnergy()<<" "<<
          it->at(ii)->GetMomentum()<<G4endl;
#endif
           }
           else
           {
              //CMS
              it->operator[](ii)->Lorentz(*(it->operator[](ii)), -1.*theCMS);
#ifdef G4PHPDEBUG
      if( std::getenv("G4ParticleHPDebug") ) 
        G4cout <<"G4ParticleHPEnAngCorrelation: after Lorentz boost "<<
        it->at(ii)->GetKineticEnergy()<<" "<<
          it->at(ii)->GetMomentum()<<G4endl;
#endif
           }
        }
    else
    {
      throw G4HadronicException(__FILE__, __LINE__, "G4ParticleHPEnAngCorrelation::Sample: The frame of the finalstate is not specified");
    }
          if( std::getenv("G4PHPTEST") ) G4cout << frameFlag << " G4particleHPEnAngCorrelation COS THETA after Lorentz " <<  std::cos(it->operator[](ii)->GetMomentum().theta()) << G4endl;
 
    // }//getenv("G4PHP_NO_LORENTZ_BOOST")) 
    //  G4cout <<  ii << " EnAnG energy after boost " << it->operator[](ii)->GetKineticEnergy() << G4endl; //GDEB
    result->push_back(it->operator[](ii));
      }
      delete it;
    }
  }   
 
  return result;
}

References CLHEP::HepLorentzVector::e(), fCache, frameFlag, G4cout, G4endl, G4ReactionProduct::GetMomentum(), G4ParticleHPProduct::GetMultiplicity(), G4ParticleHPProduct::GetQValue(), G4ReactionProduct::GetTotalEnergy(), CLHEP::HepLorentzRotation::inverse(), G4ReactionProduct::Lorentz(), G4ParticleHPProduct::MeanEnergyOfThisInteraction(), nProducts, CLHEP::HepLorentzVector::phi(), CLHEP::HepLorentzRotation::rotateY(), CLHEP::HepLorentzRotation::rotateZ(), G4ParticleHPProduct::Sample(), G4ReactionProduct::SetMass(), G4ReactionProduct::SetMomentum(), G4ReactionProduct::SetTotalEnergy(), theProducts, CLHEP::HepLorentzVector::theta(), and CLHEP::HepLorentzVector::vect().

Referenced by G4ParticleHPCaptureFS::ApplyYourself(), G4ParticleHPInelasticBaseFS::BaseApply(), and G4ParticleHPInelasticCompFS::CompositeApply().

SampleOne()

G4ReactionProduct * G4ParticleHPEnAngCorrelation::SampleOne

(

G4double

anEnergy

)

Definition at line 42 of file G4ParticleHPEnAngCorrelation.cc.

{ 
  G4ReactionProduct * result = new G4ReactionProduct;
  
  // do we have an appropriate distribution
   if(nProducts!=1) throw G4HadronicException(__FILE__, __LINE__, "More than one product in SampleOne");
  
  // get the result
  G4ReactionProductVector * temp=0;
  G4int i=0;
 
  G4int icounter=0;
  G4int icounter_max=1024;
  while(temp == 0) {
      icounter++;
      if ( icounter > icounter_max ) {
     G4cout << "Loop-counter exceeded the threshold value at " << __LINE__ << "th line of " << __FILE__ << "." << G4endl;
         break;
      }
     temp = theProducts[i++].Sample(anEnergy,1); 
  }
  
  // is the multiplicity correct
  if(temp->size()!=1) throw G4HadronicException(__FILE__, __LINE__, "SampleOne: Yield not correct");
  
  // fill result
  result = temp->operator[](0);
  
  // some garbage collection
  delete temp;
  
  // return result
  return result;
}

References G4cout, G4endl, nProducts, G4ParticleHPProduct::Sample(), and theProducts.

SetProjectileRP()

void G4ParticleHPEnAngCorrelation::SetProjectileRP ( G4ReactionProduct &  aIncidentPart )

inline

Definition at line 114 of file G4ParticleHPEnAngCorrelation.hh.

  {
    fCache.Get().theProjectileRP = &aIncidentPart;
    for (G4int i=0;i<nProducts;i++)
      theProducts[i].SetProjectileRP(fCache.Get().theProjectileRP);
  }

References fCache, nProducts, SetProjectileRP(), and theProducts.

Referenced by G4ParticleHPCaptureFS::ApplyYourself(), G4ParticleHPInelasticBaseFS::BaseApply(), G4ParticleHPInelasticCompFS::InitDistributionInitialState(), and SetProjectileRP().

SetTarget()

void G4ParticleHPEnAngCorrelation::SetTarget ( G4ReactionProduct &  aTarget )

inline

Definition at line 107 of file G4ParticleHPEnAngCorrelation.hh.

  {
    fCache.Get().theTarget = &aTarget;
    for (G4int i=0;i<nProducts;i++)
      theProducts[i].SetTarget(fCache.Get().theTarget);
  }

References fCache, nProducts, SetTarget(), and theProducts.

Referenced by G4ParticleHPCaptureFS::ApplyYourself(), G4ParticleHPInelasticBaseFS::BaseApply(), G4ParticleHPInelasticCompFS::InitDistributionInitialState(), and SetTarget().

Field Documentation

fCache

G4Cache<toBeCached> G4ParticleHPEnAngCorrelation::fCache

private

Definition at line 148 of file G4ParticleHPEnAngCorrelation.hh.

Referenced by G4ParticleHPEnAngCorrelation(), GetTotalMeanEnergy(), Sample(), SetProjectileRP(), and SetTarget().

frameFlag

G4int G4ParticleHPEnAngCorrelation::frameFlag

private

Definition at line 141 of file G4ParticleHPEnAngCorrelation.hh.

Referenced by G4ParticleHPEnAngCorrelation(), Init(), and Sample().

inCharge

G4bool G4ParticleHPEnAngCorrelation::inCharge

private

Definition at line 144 of file G4ParticleHPEnAngCorrelation.hh.

Referenced by G4ParticleHPEnAngCorrelation(), InCharge(), and Init().

nProducts

G4int G4ParticleHPEnAngCorrelation::nProducts

private

Definition at line 142 of file G4ParticleHPEnAngCorrelation.hh.

Referenced by G4ParticleHPEnAngCorrelation(), Init(), Sample(), SampleOne(), SetProjectileRP(), and SetTarget().

targetMass

G4double G4ParticleHPEnAngCorrelation::targetMass

private

Definition at line 140 of file G4ParticleHPEnAngCorrelation.hh.

Referenced by G4ParticleHPEnAngCorrelation(), GetTargetMass(), and Init().

theProducts

G4ParticleHPProduct* G4ParticleHPEnAngCorrelation::theProducts

private

Definition at line 143 of file G4ParticleHPEnAngCorrelation.hh.

Referenced by G4ParticleHPEnAngCorrelation(), Init(), Sample(), SampleOne(), SetProjectileRP(), SetTarget(), and ~G4ParticleHPEnAngCorrelation().

theProjectile

G4ParticleDefinition* G4ParticleHPEnAngCorrelation::theProjectile

private

Definition at line 150 of file G4ParticleHPEnAngCorrelation.hh.

Referenced by G4ParticleHPEnAngCorrelation(), and Init().

---

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

• source/processes/hadronic/models/particle_hp/include/G4ParticleHPEnAngCorrelation.hh
• source/processes/hadronic/models/particle_hp/src/G4ParticleHPEnAngCorrelation.cc