G4NeutronHPorLCapture Class Reference

#include <G4NeutronHPorLCapture.hh>

Inheritance diagram for G4NeutronHPorLCapture:

Public Member Functions
	G4NeutronHPorLCapture ()
	~G4NeutronHPorLCapture ()
G4HadFinalState *	ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus)
virtual const std::pair< G4double, G4double >	GetFatalEnergyCheckLevels () const
G4int	GetNiso ()
G4bool	IsThisElementOK (G4String)
G4VCrossSectionDataSet *	GiveXSectionDataSet ()

---

Detailed Description

Definition at line 57 of file G4NeutronHPorLCapture.hh.

---

Constructor & Destructor Documentation

G4NeutronHPorLCapture::G4NeutronHPorLCapture

(

)

Definition at line 43 of file G4NeutronHPorLCapture.cc.

References G4cout, G4endl, G4Element::GetElementTable(), G4Element::GetNumberOfElements(), G4HadronicInteraction::SetMaxEnergy(), and G4HadronicInteraction::SetMinEnergy().

  :G4HadronicInteraction("NeutronHPorLCapture")
{
   G4NeutronHPCaptureFS * theFS = new G4NeutronHPCaptureFS;
   if(!getenv("G4NEUTRONHPDATA")) 
       throw G4HadronicException(__FILE__, __LINE__, "Please setenv G4NEUTRONHPDATA to point to the neutron cross-section files.");
   dirName = getenv("G4NEUTRONHPDATA");
   G4String tString = "/Capture/";
   dirName = dirName + tString;
//    G4cout <<"G4NeutronHPorLCapture::G4NeutronHPorLCapture testit "<<dirName<<G4endl;
   numEle = G4Element::GetNumberOfElements();
   theCapture = new G4NeutronHPChannel[numEle];
   unavailable_elements.clear();
   for (G4int i=0; i<numEle; i++)
   {
      theCapture[i].Init((*(G4Element::GetElementTable()))[i], dirName);
      //G4cout << (*(G4Element::GetElementTable()))[i] -> GetName()  << G4endl;
      //while(!theCapture[i].Register(theFS));
      try { while(!theCapture[i].Register(theFS)) ; }
      catch ( G4HadronicException )
      {
          unavailable_elements.insert ( (*(G4Element::GetElementTable()))[i]->GetName() ); 
      }
   }
   delete theFS;
   SetMinEnergy(0.*eV);
   SetMaxEnergy(20.*MeV);
   if ( unavailable_elements.size() > 0 )
   {
      std::set< G4String>::iterator it;
      G4cout << "HP Capture data are not available for thess  elements "<< G4endl;
      for ( it = unavailable_elements.begin() ; it != unavailable_elements.end() ; it++ )
         G4cout << *it << G4endl;
      G4cout << "Low Energy Parameterization Models will be used."<< G4endl;
   }

   createXSectionDataSet();
}

G4NeutronHPorLCapture::~G4NeutronHPorLCapture

(

)

Definition at line 82 of file G4NeutronHPorLCapture.cc.

{
   delete [] theCapture;
   delete theDataSet; 
}

---

Member Function Documentation

G4HadFinalState * G4NeutronHPorLCapture::ApplyYourself	(	const G4HadProjectile &	aTrack,
		G4Nucleus &	aTargetNucleus
	)			[virtual]

Implements G4HadronicInteraction.

Definition at line 90 of file G4NeutronHPorLCapture.cc.

References G4NeutronHPChannel::ApplyYourself(), G4UniformRand, G4Material::GetElement(), G4Element::GetIndex(), G4HadProjectile::GetMaterial(), G4Material::GetNumberOfElements(), G4Material::GetTemperature(), G4NeutronHPThermalBoost::GetThermalEnergy(), G4Material::GetVecNbOfAtomsPerVolume(), G4NeutronHPChannel::GetXsec(), and CLHEP::detail::n.

Referenced by G4NeutronHPorLCaptureModel::ApplyYourself().

{
   const G4Material * theMaterial = aTrack.GetMaterial();
   G4int n = theMaterial->GetNumberOfElements();
   G4int index = theMaterial->GetElement(0)->GetIndex();
   if(n!=1)
   {
      G4int i;
      xSec = new G4double[n];
      G4double sum=0;
      const G4double * NumAtomsPerVolume = theMaterial->GetVecNbOfAtomsPerVolume();
      G4double rWeight;    
      G4NeutronHPThermalBoost aThermalE;
      for (i=0; i<n; i++)
      {
        index = theMaterial->GetElement(i)->GetIndex();
        rWeight = NumAtomsPerVolume[i];
        G4double x = aThermalE.GetThermalEnergy(aTrack, theMaterial->GetElement(i), theMaterial->GetTemperature());

        //xSec[i] = theCapture[index].GetXsec(aThermalE.GetThermalEnergy(aTrack,
        //                                                                   theMaterial->GetElement(i),
        //                                                                   theMaterial->GetTemperature()));
        xSec[i] = theCapture[index].GetXsec(x);

        xSec[i] *= rWeight;
        sum+=xSec[i];
      }
      G4double random = G4UniformRand();
      G4double running = 0;
      for (i=0; i<n; i++)
      {
        running += xSec[i];
        index = theMaterial->GetElement(i)->GetIndex();
        if(random<=running/sum) break;
      }
      delete [] xSec;
      // it is element-wise initialised.
    }
    return theCapture[index].ApplyYourself(aTrack); 
}

const std::pair< G4double, G4double > G4NeutronHPorLCapture::GetFatalEnergyCheckLevels

(

)

const [virtual]

Reimplemented from G4HadronicInteraction.

Definition at line 147 of file G4NeutronHPorLCapture.cc.

References DBL_MAX.

{
   //return std::pair<G4double, G4double>(10*perCent,10*GeV);
   return std::pair<G4double, G4double>(10*perCent,DBL_MAX);
}

G4int G4NeutronHPorLCapture::GetNiso

(

)

[inline]

Definition at line 69 of file G4NeutronHPorLCapture.hh.

References G4NeutronHPChannel::GetNiso().

{return theCapture[0].GetNiso();}

G4VCrossSectionDataSet* G4NeutronHPorLCapture::GiveXSectionDataSet

(

)

[inline]

Definition at line 83 of file G4NeutronHPorLCapture.hh.

Referenced by G4NeutronHPorLCaptureModel::GiveHPXSectionDataSet().

{ return theDataSet; }; 

G4bool G4NeutronHPorLCapture::IsThisElementOK

(

G4String

)

Definition at line 133 of file G4NeutronHPorLCapture.cc.

References FALSE, and TRUE.

Referenced by G4NeutronHPorLCaptureModel::ApplyYourself().

{
   if ( unavailable_elements.find( name ) == unavailable_elements.end() ) 
      return TRUE;
   else 
      return FALSE; 
}

---

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

• G4NeutronHPorLCapture.hh
• G4NeutronHPorLCapture.cc

---