G4ParticleHPChannel Class Reference

#include <G4ParticleHPChannel.hh>

Public Member Functions
	G4ParticleHPChannel (G4ParticleDefinition *projectile=nullptr)
	~G4ParticleHPChannel ()
G4double	GetXsec (G4double energy)
G4double	GetWeightedXsec (G4double energy, G4int isoNumber)
G4double	GetFSCrossSection (G4double energy, G4int isoNumber)
G4bool	IsActive (G4int isoNumber) const
G4bool	HasFSData (G4int isoNumber)
G4bool	HasAnyData (G4int isoNumber)
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 ()
void	DumpInfo ()
G4String &	GetFSType ()
G4ParticleHPFinalState **	GetFinalStates () 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;
}

~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 227 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 G4FissLib::ApplyYourself(), G4NeutronHPElasticVI::ApplyYourself(), and G4ParticleHPChannelList::ApplyYourself().

DumpInfo()

void G4ParticleHPChannel::DumpInfo

(

)

Definition at line 320 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;
    }
  }
}

Referenced by G4ParticleHPChannelList::DumpInfo().

GetFinalStates()

G4ParticleHPFinalState ** G4ParticleHPChannel::GetFinalStates ( ) const

inline

Definition at line 128 of file G4ParticleHPChannel.hh.

{ return theFinalStates; }

GetFSCrossSection()

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

Definition at line 90 of file G4ParticleHPChannel.cc.

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

Referenced by G4ParticleHPChannelList::ApplyYourself().

GetFSType()

G4String & G4ParticleHPChannel::GetFSType ( )

inline

Definition at line 126 of file G4ParticleHPChannel.hh.

{ return theFSType; }

GetM()

G4double G4ParticleHPChannel::GetM ( G4int i ) const

inline

Definition at line 110 of file G4ParticleHPChannel.hh.

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

Referenced by ApplyYourself().

GetN()

G4double G4ParticleHPChannel::GetN ( G4int i ) const

inline

Definition at line 108 of file G4ParticleHPChannel.hh.

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

Referenced by ApplyYourself(), and G4ParticleHPChannelList::ApplyYourself().

GetNiso()

G4int G4ParticleHPChannel::GetNiso ( ) const

inline

Definition at line 106 of file G4ParticleHPChannel.hh.

{ return niso; }

Referenced by G4ParticleHPChannelList::ApplyYourself().

GetWeightedXsec()

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

Definition at line 84 of file G4ParticleHPChannel.cc.

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

Referenced by G4ParticleHPChannelList::ApplyYourself().

GetXsec()

G4double G4ParticleHPChannel::GetXsec

(

G4double

energy

)

Definition at line 79 of file G4ParticleHPChannel.cc.

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

Referenced by G4FissLib::ApplyYourself().

GetZ()

G4double G4ParticleHPChannel::GetZ ( G4int i ) const

inline

Definition at line 109 of file G4ParticleHPChannel.hh.

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

Referenced by ApplyYourself(), and G4ParticleHPChannelList::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);
      if (std::abs(std::abs(xp - xa) / xa) < 0.001) {
        ++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 )

inline

Definition at line 79 of file G4ParticleHPChannel.hh.

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

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

HasDataInAnyFinalState()

G4bool G4ParticleHPChannel::HasDataInAnyFinalState ( )

inline

Definition at line 112 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 )

inline

Definition at line 74 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 G4FissLib::G4FissLib(), Init(), and G4ParticleHPChannelList::Register().

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 69 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;
}

Referenced by G4FissLib::G4FissLib(), and G4ParticleHPChannelList::Register().

UpdateData() [1/2]

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

inline

Definition at line 91 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
  {
    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 135 of file G4ParticleHPChannel.hh.

Referenced by ApplyYourself(), and G4ParticleHPChannel().

---

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

• G4ParticleHPChannel.hh
• G4ParticleHPChannel.cc