G4ParticleHPEnAngCorrelation Class Reference

#include <G4ParticleHPEnAngCorrelation.hh>

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

Detailed Description

Definition at line 46 of file G4ParticleHPEnAngCorrelation.hh.

Constructor & Destructor Documentation

G4ParticleHPEnAngCorrelation() [1/2]

G4ParticleHPEnAngCorrelation::G4ParticleHPEnAngCorrelation ( const G4ParticleDefinition * proj = nullptr )

explicit

Definition at line 50 of file G4ParticleHPEnAngCorrelation.cc.

{
  theProjectile = (nullptr == proj) ? G4Neutron::Neutron() : proj;
  toBeCached val;
  fCache.Put( val );
}

Referenced by G4ParticleHPEnAngCorrelation(), and operator=().

~G4ParticleHPEnAngCorrelation()

G4ParticleHPEnAngCorrelation::~G4ParticleHPEnAngCorrelation

(

)

Definition at line 57 of file G4ParticleHPEnAngCorrelation.cc.

{
  delete[] theProducts;
}

G4ParticleHPEnAngCorrelation() [2/2]

G4ParticleHPEnAngCorrelation::G4ParticleHPEnAngCorrelation ( G4ParticleHPEnAngCorrelation & )

delete

Member Function Documentation

GetNumberOfProducts()

G4double G4ParticleHPEnAngCorrelation::GetNumberOfProducts ( )

inline

Definition at line 84 of file G4ParticleHPEnAngCorrelation.hh.

{ return nProducts; }

GetTargetMass()

G4double G4ParticleHPEnAngCorrelation::GetTargetMass ( )

inline

Definition at line 82 of file G4ParticleHPEnAngCorrelation.hh.

{ return targetMass; }

GetTotalMeanEnergy()

G4double G4ParticleHPEnAngCorrelation::GetTotalMeanEnergy ( )

inline

Definition at line 86 of file G4ParticleHPEnAngCorrelation.hh.

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

InCharge()

G4bool G4ParticleHPEnAngCorrelation::InCharge ( )

inline

Definition at line 80 of file G4ParticleHPEnAngCorrelation.hh.

{ return inCharge; }

Init()

void G4ParticleHPEnAngCorrelation::Init

(

std::istream &

aDataFile

)

Definition at line 62 of file G4ParticleHPEnAngCorrelation.cc.

{
  inCharge = true;
  aDataFile >> targetMass >> frameFlag >> nProducts;
  const std::size_t psize = nProducts > 0 ? nProducts : 1;
  theProducts = new G4ParticleHPProduct[psize];
  for (std::size_t i = 0; i < psize; ++i) {
    theProducts[i].Init(aDataFile, theProjectile);
  }
}

operator=()

G4ParticleHPEnAngCorrelation & G4ParticleHPEnAngCorrelation::operator= ( const G4ParticleHPEnAngCorrelation & right )

delete

Sample()

G4ReactionProductVector * G4ParticleHPEnAngCorrelation::Sample

(

G4double

anEnergy

)

Definition at line 116 of file G4ParticleHPEnAngCorrelation.cc.

{
  auto result = new G4ReactionProductVector;
  G4int i;
  G4ReactionProductVector* it;
  G4ReactionProduct theCMS;
  G4LorentzRotation toZ;
  /*
  G4cout << "G4ParticleHPEnAngCorrelation::Sample for E=" << anEnergy 
     << " flag=" << frameFlag << " nProducts=" << nProducts <<G4endl;
  */
  if (frameFlag == 2 || frameFlag == 3)  // Added for particle HP
  {
    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 = 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);
    it = theProducts[i].Sample(anEnergy, nPart);
    G4double aMeanEnergy = theProducts[i].MeanEnergyOfThisInteraction();
 
    // 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 != nullptr) {
      for (unsigned int ii = 0; ii < it->size(); ++ii) {
        G4LorentzVector pTmp1(it->operator[](ii)->GetMomentum(),
                              it->operator[](ii)->GetTotalEnergy());
        pTmp1 = toLab * pTmp1;
        it->operator[](ii)->SetMomentum(pTmp1.vect());
        it->operator[](ii)->SetTotalEnergy(pTmp1.e());
 
        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 (G4ParticleHPManager::GetInstance()->GetDEBUG())
            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 (G4ParticleHPManager::GetInstance()->GetDEBUG())
            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 (G4ParticleHPManager::GetInstance()->GetDEBUG())
              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 (G4ParticleHPManager::GetInstance()->GetDEBUG())
              G4cout << "G4ParticleHPEnAngCorrelation: after Lorentz boost "
                     << it->at(ii)->GetKineticEnergy() << " "
                     << it->at(ii)->GetMomentum() << G4endl;
#endif
          }
        }
        else
        {
          throw G4HadronicException(__FILE__, __LINE__,
            "Sample: The frame of the final state is not specified");
        }
        result->push_back(it->operator[](ii));
      }
      delete it;
    }
  }
  return result;
}

SampleOne()

G4ReactionProduct * G4ParticleHPEnAngCorrelation::SampleOne

(

G4double

anEnergy

)

Definition at line 73 of file G4ParticleHPEnAngCorrelation.cc.

{
  auto 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 = nullptr;
  G4int i = 0;
 
  G4int icounter = 0;
  G4int icounter_max = 1024;
  while (temp == nullptr)
  {
    ++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);
  }
  if (nullptr != temp)
  {
    if (temp->size() == 1)
    {
      result = (*temp)[0];
    }
    else
    {
      for (auto const& ptr : *temp) { delete ptr; }
      throw G4HadronicException(__FILE__, __LINE__,
                "SampleOne: Yield not correct");
    }
  }
  delete temp;
  return result;
}

SetProjectileRP()

void G4ParticleHPEnAngCorrelation::SetProjectileRP ( G4ReactionProduct & aIncidentPart )

inline

Definition at line 73 of file G4ParticleHPEnAngCorrelation.hh.

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

Referenced by SetProjectileRP().

SetTarget()

void G4ParticleHPEnAngCorrelation::SetTarget ( G4ReactionProduct & aTarget )

inline

Definition at line 66 of file G4ParticleHPEnAngCorrelation.hh.

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

Referenced by SetTarget().

---

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

• G4ParticleHPEnAngCorrelation.hh
• G4ParticleHPEnAngCorrelation.cc