G4NeutronHPEnAngCorrelation Class Reference

#include <G4NeutronHPEnAngCorrelation.hh>

Public Member Functions
	G4NeutronHPEnAngCorrelation ()
	~G4NeutronHPEnAngCorrelation ()
void	Init (std::istream &aDataFile)
G4ReactionProduct *	SampleOne (G4double anEnergy)
G4ReactionProductVector *	Sample (G4double anEnergy)
void	SetTarget (G4ReactionProduct &aTarget)
void	SetNeutron (G4ReactionProduct &aNeutron)
G4bool	InCharge ()
G4double	GetTargetMass ()
G4double	GetTotalMeanEnergy ()

Detailed Description

Definition at line 40 of file G4NeutronHPEnAngCorrelation.hh.

Constructor & Destructor Documentation

G4NeutronHPEnAngCorrelation::G4NeutronHPEnAngCorrelation ( )

inline

Definition at line 43 of file G4NeutronHPEnAngCorrelation.hh.

  {
    theProducts = 0;
    inCharge = false;
    theTotalMeanEnergy = -1.;
  }

G4NeutronHPEnAngCorrelation::~G4NeutronHPEnAngCorrelation ( )

inline

Definition at line 49 of file G4NeutronHPEnAngCorrelation.hh.

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

Member Function Documentation

G4double G4NeutronHPEnAngCorrelation::GetTargetMass ( )

inline

Definition at line 86 of file G4NeutronHPEnAngCorrelation.hh.

Referenced by G4NeutronHPInelasticBaseFS::BaseApply().

{ return targetMass; }

G4double G4NeutronHPEnAngCorrelation::GetTotalMeanEnergy ( )

inline

Definition at line 88 of file G4NeutronHPEnAngCorrelation.hh.

Referenced by G4NeutronHPInelasticBaseFS::BaseApply().

  {
     // cashed in 'sample' call
    return theTotalMeanEnergy; 
  }

G4bool G4NeutronHPEnAngCorrelation::InCharge ( )

inline

Definition at line 81 of file G4NeutronHPEnAngCorrelation.hh.

  {
    return inCharge;
  }

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

inline

Definition at line 54 of file G4NeutronHPEnAngCorrelation.hh.

References G4NeutronHPProduct::Init().

Referenced by G4NeutronHPCaptureFS::Init(), G4NeutronHPInelasticBaseFS::Init(), and G4NeutronHPInelasticCompFS::Init().

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

G4ReactionProductVector * G4NeutronHPEnAngCorrelation::Sample

(

G4double

anEnergy

)

Definition at line 62 of file G4NeutronHPEnAngCorrelation.cc.

References G4endl, G4ReactionProduct::GetKineticEnergy(), G4ReactionProduct::GetMomentum(), G4NeutronHPProduct::GetQValue(), G4ReactionProduct::GetTotalEnergy(), CLHEP::HepLorentzRotation::inverse(), G4ReactionProduct::Lorentz(), G4NeutronHPProduct::MeanEnergyOfThisInteraction(), CLHEP::HepLorentzRotation::rotateY(), CLHEP::HepLorentzRotation::rotateZ(), G4NeutronHPProduct::Sample(), G4ReactionProduct::SetMass(), G4ReactionProduct::SetMomentum(), and G4ReactionProduct::SetTotalEnergy().

Referenced by G4NeutronHPCaptureFS::ApplyYourself(), G4NeutronHPInelasticBaseFS::BaseApply(), and G4NeutronHPInelasticCompFS::CompositeApply().

{
  G4ReactionProductVector * result = new G4ReactionProductVector;
  G4int i;
  G4ReactionProductVector * it;
  G4ReactionProduct theCMS;
  G4LorentzRotation toZ;
  //TK120515 migrate frameFlag (MF6 LCT) = 3 
  //if(frameFlag==2)
  if(frameFlag==2||frameFlag==3)
  {
    // simplify and double check @
    G4ThreeVector the3Neutron = theNeutron.GetMomentum(); //theNeutron has value in LAB
    G4double nEnergy = theNeutron.GetTotalEnergy();
    G4ThreeVector the3Target = theTarget.GetMomentum();  //theTarget has value in LAB
    G4double tEnergy = theTarget.GetTotalEnergy();
    G4double totE = nEnergy+tEnergy;
    G4ThreeVector the3CMS = the3Target+the3Neutron;
    theCMS.SetMomentum(the3CMS);
    G4double cmsMom = std::sqrt(the3CMS*the3CMS);
    G4double sqrts = std::sqrt((totE-cmsMom)*(totE+cmsMom));
    theCMS.SetMass(sqrts);
    theCMS.SetTotalEnergy(totE);
    G4ReactionProduct aNeutron;
    aNeutron.Lorentz(theNeutron, 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 = aNeutron.GetKineticEnergy();
    anEnergy = theNeutron.GetKineticEnergy(); //should be same argumment of "anEnergy"

    G4LorentzVector Ptmp (aNeutron.GetMomentum(), aNeutron.GetTotalEnergy());

    toZ.rotateZ(-1*Ptmp.phi());
    toZ.rotateY(-1*Ptmp.theta());
  }
  theTotalMeanEnergy=0;
  G4LorentzRotation toLab(toZ.inverse()); //toLab only change axis NOT to LAB system
  for(i=0; i<nProducts; i++)
  {
    it = theProducts[i].Sample(anEnergy);
    G4double aMeanEnergy = theProducts[i].MeanEnergyOfThisInteraction();
    if(aMeanEnergy>0)
    {
      theTotalMeanEnergy += aMeanEnergy;
    }
    else
    {
      theTotalMeanEnergy = anEnergy/nProducts+theProducts[i].GetQValue();
    }
    if(it!=0)
    {
      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.*theTarget); //TK 100413 Is this really need?
    }
    else if(frameFlag==2) // CMS
    {
#ifdef G4_NHP_DEBUG
          cout <<"G4NeutronHPEnAngCorrelation: "<<
             it->at(ii)->GetTotalEnergy()<<" "<<
             it->at(ii)->GetMomentum()<<G4endl;
#endif
          it->operator[](ii)->Lorentz(*(it->operator[](ii)), -1.*theCMS);
    }
        //TK120515 migrate frameFlag (MF6 LCT) = 3 
    else if(frameFlag==3) // CMS A<=4 other LAB
        {
           if ( theProducts[i].GetMassCode() > 4 ) //Alpha AWP 3.96713
           {
              //LAB
              it->operator[](ii)->Lorentz(*(it->operator[](ii)), -1.*theTarget); //TK 100413 Is this really need?
           }
           else
           {
              //CMS
              it->operator[](ii)->Lorentz(*(it->operator[](ii)), -1.*theCMS);
           }
        }
    else
    {
          throw G4HadronicException(__FILE__, __LINE__, "G4NeutronHPEnAngCorrelation::Sample: The frame of the finalstate is not specified");
    }
    result->push_back(it->operator[](ii));
      }
    delete it;
    }
  }   
  return result;
}

G4ReactionProduct * G4NeutronHPEnAngCorrelation::SampleOne

(

G4double

anEnergy

)

Definition at line 37 of file G4NeutronHPEnAngCorrelation.cc.

References G4NeutronHPProduct::Sample().

{  
  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;
  while(temp == 0) temp = theProducts[i++].Sample(anEnergy);
  
  // 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;
}

void G4NeutronHPEnAngCorrelation::SetNeutron ( G4ReactionProduct &  aNeutron )

inline

Definition at line 75 of file G4NeutronHPEnAngCorrelation.hh.

Referenced by G4NeutronHPCaptureFS::ApplyYourself(), G4NeutronHPInelasticBaseFS::BaseApply(), and G4NeutronHPInelasticCompFS::InitDistributionInitialState().

  {
    theNeutron = aNeutron;
    for(G4int i=0;i<nProducts;i++)theProducts[i].SetNeutron(&theNeutron);
  }

void G4NeutronHPEnAngCorrelation::SetTarget ( G4ReactionProduct &  aTarget )

inline

Definition at line 69 of file G4NeutronHPEnAngCorrelation.hh.

Referenced by G4NeutronHPCaptureFS::ApplyYourself(), G4NeutronHPInelasticBaseFS::BaseApply(), and G4NeutronHPInelasticCompFS::InitDistributionInitialState().

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

---

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

• G4NeutronHPEnAngCorrelation.hh
• G4NeutronHPEnAngCorrelation.cc