G4ParticleHPChannel Class Reference

#include <G4ParticleHPChannel.hh>

Public Member Functions
	G4ParticleHPChannel (G4ParticleDefinition *projectile=nullptr)
	~G4ParticleHPChannel ()
G4double	GetXsec (G4double energy) const
G4double	GetWeightedXsec (G4double energy, G4int isoNumber) const
G4double	GetFSCrossSection (G4double energy, G4int isoNumber) const
G4bool	IsActive (G4int isoNumber) const
G4bool	HasFSData (G4int isoNumber) const
G4bool	HasAnyData (G4int isoNumber) const
G4bool	Register (G4ParticleHPFinalState *theFS)
void	Init (G4Element *theElement, const G4String &dirName)
void	Init (G4Element *theElement, const G4String &dirName, const G4String &fsType)
void	UpdateData (G4int A, G4int Z, G4int index, G4double abundance, G4ParticleDefinition *projectile)
void	UpdateData (G4int A, G4int Z, G4int M, G4int index, G4double abundance, G4ParticleDefinition *projectile)
void	Harmonise (G4ParticleHPVector *&theStore, G4ParticleHPVector *theNew)
G4HadFinalState *	ApplyYourself (const G4HadProjectile &theTrack, G4int isoNumber=-1, G4bool isElastic=false)
G4int	GetNiso () const
G4double	GetN (G4int i) const
G4double	GetZ (G4int i) const
G4double	GetM (G4int i) const
G4bool	HasDataInAnyFinalState () const
void	DumpInfo () const
const G4String &	GetFSType ()
G4ParticleHPFinalState **	GetFinalStates () const
G4WendtFissionFragmentGenerator *	GetWendtFissionGenerator () const
	G4ParticleHPChannel (G4ParticleHPChannel &)=delete
G4ParticleHPChannel &	operator= (const G4ParticleHPChannel &right)=delete

Protected Attributes
G4ParticleHPManager *	fManager

Detailed Description

Definition at line 56 of file G4ParticleHPChannel.hh.

Constructor & Destructor Documentation

G4ParticleHPChannel() [1/2]

G4ParticleHPChannel::G4ParticleHPChannel

(

G4ParticleDefinition *

projectile = nullptr

)

Definition at line 52 of file G4ParticleHPChannel.cc.

{
  fManager = G4ParticleHPManager::GetInstance();
  if (fManager->GetUseWendtFissionModel()) {
    wendtFissionGenerator = G4WendtFissionFragmentGenerator::GetInstance();
    // Make sure both fission fragment models are not active at same time
    fManager->SetProduceFissionFragments(false);
  }
  theProjectile = (nullptr == p) ? G4Neutron::Neutron() : p;
  theChannelData = new G4ParticleHPVector;
}

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

~G4ParticleHPChannel()

G4ParticleHPChannel::~G4ParticleHPChannel

(

)

Definition at line 64 of file G4ParticleHPChannel.cc.

{
  delete theChannelData;
  // Following statement disabled to avoid SEGV
  // theBuffer is also deleted as "theChannelData" in
  delete[] theIsotopeWiseData;
  if (theFinalStates != nullptr) {
    for (G4int i = 0; i < niso; i++) {
      delete theFinalStates[i];
    }
    delete[] theFinalStates;
  }
  delete[] active;
}

G4ParticleHPChannel() [2/2]

G4ParticleHPChannel::G4ParticleHPChannel ( G4ParticleHPChannel & )

delete

Member Function Documentation

ApplyYourself()

G4HadFinalState * G4ParticleHPChannel::ApplyYourself	(	const G4HadProjectile &	theTrack,
		G4int	isoNumber = -1,
		G4bool	isElastic = false )

Definition at line 234 of file G4ParticleHPChannel.cc.

{
  //G4cout << "G4ParticleHPChannel::ApplyYourself niso=" << niso
  //     << " ni=" << anIsotope << " isElastic=" << isElastic <<G4endl;
  if (anIsotope != -1 && anIsotope != -2) {
    // Inelastic Case
    //G4cout << "G4ParticleHPChannel Inelastic Case"
    //<< " Z= " << GetZ(anIsotope) << " A = " << GetN(anIsotope) << G4endl;
    fManager->GetReactionWhiteBoard()->SetTargA((G4int)GetN(anIsotope));
    fManager->GetReactionWhiteBoard()->SetTargZ((G4int)GetZ(anIsotope));
    return theFinalStates[anIsotope]->ApplyYourself(theTrack);
  }
  G4double sum = 0;
  G4int it = 0;
  auto xsec = new G4double[niso];
  G4ParticleHPThermalBoost aThermalE;
  for (G4int i = 0; i < niso; i++) {
    if (theFinalStates[i]->HasAnyData()) {
      /*
      G4cout << "FS: " << i << theTrack.GetDefinition()->GetParticleName()
         << " Z=" << theFinalStates[i]->GetZ() 
         << " A=" << theFinalStates[i]->GetN() 
         << G4endl;
      */
      xsec[i] = theIsotopeWiseData[i].GetXsec(
        aThermalE.GetThermalEnergy(theTrack, theFinalStates[i]->GetN(),
                                   theFinalStates[i]->GetZ(),
                                   theTrack.GetMaterial()->GetTemperature()));
      sum += xsec[i];
    }
    else {
      xsec[i] = 0;
    }
  }
  if (sum == 0) {
    it = G4lrint(niso * G4UniformRand());
  }
  else {
    G4double random = G4UniformRand();
    G4double running = 0;
    for (G4int ix = 0; ix < niso; ix++) {
      running += xsec[ix];
      if (sum == 0 || random <= running / sum) {
        it = ix;
        break;
      }
    }
    if (it == niso) it--;
  }
  delete[] xsec;
  G4HadFinalState* theFinalState = nullptr;
  const auto A = (G4int)this->GetN(it);
  const auto Z = (G4int)this->GetZ(it);
  const auto M = (G4int)this->GetM(it);
 
  //-2:Marker for Fission
  if ((wendtFissionGenerator != nullptr) && anIsotope == -2) {
    theFinalState = wendtFissionGenerator->ApplyYourself(theTrack, Z, A);
  }
 
  // Use the standard procedure if the G4FissionFragmentGenerator model fails
  if (theFinalState == nullptr) {
    G4int icounter = 0;
    G4int icounter_max = 1024;
    while (theFinalState == nullptr)  // Loop checking, 11.05.2015, T. Koi
    {
      icounter++;
      if (icounter > icounter_max) {
        G4cout << "Loop-counter exceeded the threshold value at " 
               << __LINE__ << "th line of " << __FILE__ << "." << G4endl;
        break;
      }
      if (isElastic) {
        // Register 0 K cross-section for DBRC for Doppler broadened elastic scattering kernel
        G4ParticleHPVector* xsforFS = theIsotopeWiseData[it].MakeChannelData();
        // Only G4ParticleHPElasticFS has the RegisterCrossSection method
        static_cast<G4ParticleHPElasticFS*>(theFinalStates[it])->RegisterCrossSection(xsforFS);
      }
      theFinalState = theFinalStates[it]->ApplyYourself(theTrack);
    }
  }
 
  // G4cout <<"THE IMPORTANT RETURN"<<G4endl;
  // G4cout << "TK G4ParticleHPChannel Elastic, Capture and Fission Cases "
  //<< " Z= " << this->GetZ(it) << " A = " << this->GetN(it) << G4endl;
  fManager->GetReactionWhiteBoard()->SetTargA(A);
  fManager->GetReactionWhiteBoard()->SetTargZ(Z);
  fManager->GetReactionWhiteBoard()->SetTargM(M);
 
  return theFinalState;
}

Referenced by G4NeutronHPElasticVI::ApplyYourself().

DumpInfo()

void G4ParticleHPChannel::DumpInfo

(

)

const

Definition at line 327 of file G4ParticleHPChannel.cc.

{
  G4cout << " Element: " << theElement->GetName() << G4endl;
  G4cout << " Directory name: " << theDir << G4endl;
  G4cout << " FS name: " << theFSType << G4endl;
  G4cout << " Number of Isotopes: " << niso << G4endl;
  G4cout << " Have cross sections: " << G4endl;
  for (int i = 0; i < niso; i++) {
    G4cout << theFinalStates[i]->HasXsec() << "  ";
  }
  G4cout << G4endl;
  if (theChannelData != nullptr) {
    G4cout << " Cross Section (total for this channel):" << G4endl;
    int np = theChannelData->GetVectorLength();
    G4cout << np << G4endl;
    for (int i = 0; i < np; i++) {
      G4cout << theChannelData->GetEnergy(i) / eV << "  " << theChannelData->GetXsec(i) << G4endl;
    }
  }
}

GetFinalStates()

G4ParticleHPFinalState ** G4ParticleHPChannel::GetFinalStates ( ) const

inline

Definition at line 129 of file G4ParticleHPChannel.hh.

{ return theFinalStates; }

GetFSCrossSection()

G4double G4ParticleHPChannel::GetFSCrossSection	(	G4double	energy,
		G4int	isoNumber ) const

Definition at line 90 of file G4ParticleHPChannel.cc.

{
  return theFinalStates[isoNumber]->GetXsec(energy);
}

GetFSType()

const G4String & G4ParticleHPChannel::GetFSType ( )

inline

Definition at line 127 of file G4ParticleHPChannel.hh.

{ return theFSType; }

GetM()

G4double G4ParticleHPChannel::GetM ( G4int i ) const

inline

Definition at line 111 of file G4ParticleHPChannel.hh.

{ return theFinalStates[i]->GetM(); }

Referenced by ApplyYourself().

GetN()

G4double G4ParticleHPChannel::GetN ( G4int i ) const

inline

Definition at line 109 of file G4ParticleHPChannel.hh.

{ return theFinalStates[i]->GetN(); }

Referenced by ApplyYourself().

GetNiso()

G4int G4ParticleHPChannel::GetNiso ( ) const

inline

Definition at line 107 of file G4ParticleHPChannel.hh.

{ return niso; }

GetWeightedXsec()

G4double G4ParticleHPChannel::GetWeightedXsec	(	G4double	energy,
		G4int	isoNumber ) const

Definition at line 84 of file G4ParticleHPChannel.cc.

{
  return theIsotopeWiseData[isoNumber].GetXsec(energy);
}

GetWendtFissionGenerator()

G4WendtFissionFragmentGenerator * G4ParticleHPChannel::GetWendtFissionGenerator

(

)

const

Definition at line 228 of file G4ParticleHPChannel.cc.

                                                                                     {
  if ( wendtFissionGenerator ) return wendtFissionGenerator;
  else                         return nullptr;
}

GetXsec()

G4double G4ParticleHPChannel::GetXsec

(

G4double

energy

)

const

Definition at line 79 of file G4ParticleHPChannel.cc.

{
  return std::max(0., theChannelData->GetXsec(energy));
}

GetZ()

G4double G4ParticleHPChannel::GetZ ( G4int i ) const

inline

Definition at line 110 of file G4ParticleHPChannel.hh.

{ return theFinalStates[i]->GetZ(); }

Referenced by ApplyYourself().

Harmonise()

void G4ParticleHPChannel::Harmonise	(	G4ParticleHPVector *&	theStore,
		G4ParticleHPVector *	theNew )

Definition at line 179 of file G4ParticleHPChannel.cc.

{
  G4int s_tmp = 0, n = 0, m_tmp = 0;
  auto theMerge = new G4ParticleHPVector;
  G4ParticleHPVector* anActive = theStore;
  G4ParticleHPVector* aPassive = theNew;
  G4ParticleHPVector* tmp;
  G4int a = s_tmp, p = n, t;
  while (a < anActive->GetVectorLength() && p < aPassive->GetVectorLength())
    // Loop checking, 11.05.2015, T. Koi
  {
    if (anActive->GetEnergy(a) <= aPassive->GetEnergy(p)) {
      G4double xa = anActive->GetEnergy(a);
      theMerge->SetData(m_tmp, xa, anActive->GetXsec(a) + std::max(0., aPassive->GetXsec(xa)));
      m_tmp++;
      a++;
      G4double xp = aPassive->GetEnergy(p);
      xa = std::max(xa, 0.0);
      xp = std::max(xp, 0.0);
      if (std::abs(xp - xa) < 0.001*xa) {
        ++p;
      }
    }
    else {
      tmp = anActive;
      t = a;
      anActive = aPassive;
      a = p;
      aPassive = tmp;
      p = t;
    }
  }
  while (a != anActive->GetVectorLength())  // Loop checking, 11.05.2015, T. Koi
  {
    theMerge->SetData(m_tmp++, anActive->GetEnergy(a), anActive->GetXsec(a));
    ++a;
  }
  while (p != aPassive->GetVectorLength())  // Loop checking, 11.05.2015, T. Koi
  {
    if (std::abs(theMerge->GetEnergy(std::max(0, m_tmp - 1)) -
         aPassive->GetEnergy(p)) / aPassive->GetEnergy(p) > 0.001)
      theMerge->SetData(m_tmp++, aPassive->GetEnergy(p), aPassive->GetXsec(p));
    ++p;
  }
  delete theStore;
  theStore = theMerge;
}

Referenced by UpdateData().

HasAnyData()

G4bool G4ParticleHPChannel::HasAnyData ( G4int isoNumber ) const

inline

Definition at line 80 of file G4ParticleHPChannel.hh.

  {
    return theFinalStates[isoNumber]->HasAnyData();
  }

Referenced by ApplyYourself(), HasDataInAnyFinalState(), and UpdateData().

HasDataInAnyFinalState()

G4bool G4ParticleHPChannel::HasDataInAnyFinalState ( ) const

inline

Definition at line 113 of file G4ParticleHPChannel.hh.

  {
    G4bool result = false;
    for (G4int i = 0; i < niso; ++i) {
      if (theFinalStates[i]->HasAnyData()) {
    result = true;
    break;
      }
    }
    return result;
  }

Referenced by Register().

HasFSData()

G4bool G4ParticleHPChannel::HasFSData ( G4int isoNumber ) const

inline

Definition at line 75 of file G4ParticleHPChannel.hh.

  { 
    return theFinalStates[isoNumber]->HasFSData();
  }

Init() [1/2]

void G4ParticleHPChannel::Init	(	G4Element *	theElement,
		const G4String &	dirName )

Definition at line 103 of file G4ParticleHPChannel.cc.

{
  theDir = dirName;
  theElement = anElement;
}

Referenced by Init().

Init() [2/2]

void G4ParticleHPChannel::Init	(	G4Element *	theElement,
		const G4String &	dirName,
		const G4String &	fsType )

Definition at line 96 of file G4ParticleHPChannel.cc.

{
  theFSType = aFSType;
  Init(anElement, dirName);
}

IsActive()

G4bool G4ParticleHPChannel::IsActive ( G4int isoNumber ) const

inline

Definition at line 70 of file G4ParticleHPChannel.hh.

  {
    return active[isoNumber];
  }

operator=()

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

delete

Register()

G4bool G4ParticleHPChannel::Register

(

G4ParticleHPFinalState *

theFS

)

Definition at line 109 of file G4ParticleHPChannel.cc.

{
  ++registerCount;
  G4int Z = theElement->GetZasInt();
 
  niso = (G4int)theElement->GetNumberOfIsotopes();
  const std::size_t nsize = niso > 0 ? niso : 1;
 
  delete[] theIsotopeWiseData;
  theIsotopeWiseData = new G4ParticleHPIsoData[nsize];
  delete[] active;
  active = new G4bool[nsize];
 
  delete[] theFinalStates;
  theFinalStates = new G4ParticleHPFinalState*[nsize];
  delete theChannelData;
  theChannelData = new G4ParticleHPVector;
  for (G4int i = 0; i < niso; ++i) {
    theFinalStates[i] = theFS->New();
    theFinalStates[i]->SetProjectile(theProjectile);
  }
  if (niso != 0 && registerCount == 0) {
    for (G4int i1 = 0; i1 < niso; ++i1) {
      G4int A = theElement->GetIsotope(i1)->GetN();
      G4int M = theElement->GetIsotope(i1)->Getm();
      //G4cout <<" Init: normal case i=" << i1 
      //     << " Z=" << Z << " A=" << A << G4endl;
      G4double frac = theElement->GetRelativeAbundanceVector()[i1] / perCent;
      theFinalStates[i1]->SetA_Z(A, Z, M);
      UpdateData(A, Z, M, i1, frac, theProjectile);
    }
  }
  G4bool result = HasDataInAnyFinalState();
 
  // To avoid issuing hash by worker threads
  if (result) theChannelData->Hash();
 
  return result;
}

UpdateData() [1/2]

void G4ParticleHPChannel::UpdateData ( G4int A, G4int Z, G4int index, G4double abundance, G4ParticleDefinition * projectile )

inline

Definition at line 92 of file G4ParticleHPChannel.hh.

  {
    UpdateData(A, Z, 0, index, abundance, projectile);
  }

Referenced by Register(), and UpdateData().

UpdateData() [2/2]

void G4ParticleHPChannel::UpdateData	(	G4int	A,
		G4int	Z,
		G4int	M,
		G4int	index,
		G4double	abundance,
		G4ParticleDefinition *	projectile )

Definition at line 149 of file G4ParticleHPChannel.cc.

{
  // Initialze the G4FissionFragment generator for this isomer if needed
  if (wendtFissionGenerator != nullptr) {
    wendtFissionGenerator->InitializeANucleus(A, Z, M, theDir);
  }
 
  theFinalStates[index]->Init(A, Z, M, theDir, theFSType, projectile);
  if (!theFinalStates[index]->HasAnyData()) return;
  // nothing there for exactly this isotope.
 
  // the above has put the X-sec into the FS
  theBuffer = nullptr;
  if (theFinalStates[index]->HasXsec()) {
    theBuffer = theFinalStates[index]->GetXsec();
    theBuffer->Times(abundance / 100.);
    theIsotopeWiseData[index].FillChannelData(theBuffer);
  }
  else  // get data from CrossSection directory
  {
    const G4String& tString = "/CrossSection";
    active[index] = theIsotopeWiseData[index].Init(A, Z, M, abundance,
                                                   theDir, tString);
    if (active[index]) theBuffer = theIsotopeWiseData[index].MakeChannelData();
  }
  if (theBuffer != nullptr) Harmonise(theChannelData, theBuffer);
}

Member Data Documentation

fManager

G4ParticleHPManager* G4ParticleHPChannel::fManager

protected

Definition at line 139 of file G4ParticleHPChannel.hh.

Referenced by ApplyYourself(), and G4ParticleHPChannel().

---

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

• G4ParticleHPChannel.hh
• G4ParticleHPChannel.cc