G4ParticleHPChannelList Class Reference

#include <G4ParticleHPChannelList.hh>

Public Member Functions
	G4ParticleHPChannelList (G4int n=0, G4ParticleDefinition *p=nullptr)
void	Init (G4int n)
	~G4ParticleHPChannelList ()
G4HadFinalState *	ApplyYourself (const G4Element *theElement, const G4HadProjectile &aTrack)
void	Init (G4Element *anElement, const G4String &dirName, G4ParticleDefinition *projectile)
void	Register (G4ParticleHPFinalState *theFS, const G4String &aName)
G4double	GetXsec (G4double anEnergy)
G4int	GetNumberOfChannels ()
G4bool	HasDataInAnyFinalState ()
void	DumpInfo ()
	G4ParticleHPChannelList (G4ParticleHPChannelList &)=delete
G4ParticleHPChannelList &	operator= (const G4ParticleHPChannelList &right)=delete

Detailed Description

Definition at line 46 of file G4ParticleHPChannelList.hh.

Constructor & Destructor Documentation

G4ParticleHPChannelList() [1/2]

G4ParticleHPChannelList::G4ParticleHPChannelList	(	G4int	n = 0,
		G4ParticleDefinition *	p = nullptr )

Definition at line 46 of file G4ParticleHPChannelList.cc.

{
  nChannels = n;
  theChannels = new G4ParticleHPChannel*[n];
  theProjectile = (nullptr == p) ? G4Neutron::Neutron() : p;
}

~G4ParticleHPChannelList()

G4ParticleHPChannelList::~G4ParticleHPChannelList

(

)

Definition at line 53 of file G4ParticleHPChannelList.cc.

{
  if (theChannels != nullptr) {
    for (G4int i = 0; i < nChannels; ++i) {
      delete theChannels[i];
    }
    delete[] theChannels;
  }
}

G4ParticleHPChannelList() [2/2]

G4ParticleHPChannelList::G4ParticleHPChannelList ( G4ParticleHPChannelList & )

delete

Member Function Documentation

ApplyYourself()

G4HadFinalState * G4ParticleHPChannelList::ApplyYourself	(	const G4Element *	theElement,
		const G4HadProjectile &	aTrack )

Definition at line 63 of file G4ParticleHPChannelList.cc.

{
  G4ParticleHPThermalBoost aThermalE;
  G4int i, ii;
 
  // decide on the isotope
  G4int numberOfIsos(0);
  for (ii = 0; ii < nChannels; ++ii) {
    numberOfIsos = theChannels[ii]->GetNiso();
    if (numberOfIsos != 0) break;
  }
  auto running = new G4double[numberOfIsos];
  running[0] = 0;
  for (i = 0; i < numberOfIsos; i++) {
    if (i != 0) running[i] = running[i - 1];
    for (ii = 0; ii < nChannels; ii++) {
      if (theChannels[ii]->HasAnyData(i)) {
        running[i] += theChannels[ii]->GetWeightedXsec(
          aThermalE.GetThermalEnergy(aTrack, theChannels[ii]->GetN(i), theChannels[ii]->GetZ(i),
                                     aTrack.GetMaterial()->GetTemperature()), i);
      }
    }
  }
  G4int isotope = nChannels - 1;
  G4double random = G4UniformRand();
  for (i = 0; i < numberOfIsos; ++i) {
    isotope = i;
    if (running[numberOfIsos - 1] != 0)
      if (random < running[i] / running[numberOfIsos - 1]) break;
  }
  delete[] running;
 
  // decide on the channel
  running = new G4double[nChannels];
  running[0] = 0;
  G4int targA = -1;  // For production of unChanged
  G4int targZ = -1;
  for (i = 0; i < nChannels; i++) {
    if (i != 0) running[i] = running[i - 1];
    if (theChannels[i]->HasAnyData(isotope)) {
      targA = (G4int)theChannels[i]->GetN(isotope);  // Will be simply used the last valid value
      targZ = (G4int)theChannels[i]->GetZ(isotope);
      running[i] += theChannels[i]->GetFSCrossSection(
        aThermalE.GetThermalEnergy(aTrack, targA, targZ, aTrack.GetMaterial()->GetTemperature()),
        isotope);
      targA = (G4int)theChannels[i]->GetN(isotope);  // Will be simply used the last valid value
      targZ = (G4int)theChannels[i]->GetZ(isotope);
      //  G4cout << " G4ParticleHPChannelList " << i << " targA " << targA << " targZ " << targZ <<
      //" running " << running[i] << G4endl;
    }
  }
 
  // TK120607
  if (running[nChannels - 1] == 0) {
    // This happened usually by the miss match between the cross section data and model
    if (targA == -1 && targZ == -1) {
      throw G4HadronicException(
        __FILE__, __LINE__,
        "ParticleHP model encounter lethal discrepancy with cross section data");
    }
 
    // TK121106
    G4cout << "Warning from NeutronHP: could not find proper reaction channel. This may cause by "
              "inconsistency between cross section and model.  Unchanged final states are returned."
           << G4endl;
    unChanged.Clear();
 
    // For Ep Check create unchanged final state including rest target
    G4ParticleDefinition* targ_pd = G4IonTable::GetIonTable()->GetIon(targZ, targA, 0.0);
    auto targ_dp = new G4DynamicParticle(targ_pd, G4ThreeVector(1, 0, 0), 0.0);
    unChanged.SetEnergyChange(aTrack.GetKineticEnergy());
    unChanged.SetMomentumChange(aTrack.Get4Momentum().vect().unit());
    unChanged.AddSecondary(targ_dp);
    // TK121106
    G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->SetTargA(targA);
    G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->SetTargZ(targZ);
    delete[] running;
    return &unChanged;
  }
  // TK120607
 
  G4int lChan = 0;
  random = G4UniformRand();
  for (i = 0; i < nChannels; i++) {
    lChan = i;
    if (running[nChannels - 1] != 0)
      if (random < running[i] / running[nChannels - 1]) break;
  }
  delete[] running;
#ifdef G4PHPDEBUG
  if (G4ParticleHPManager::GetInstance()->GetDEBUG())
    G4cout << " G4ParticleHPChannelList SELECTED ISOTOPE " << isotope << " SELECTED CHANNEL "
           << lChan << G4endl;
#endif
  return theChannels[lChan]->ApplyYourself(aTrack, isotope);
}

Referenced by G4NeutronHPInelasticVI::ApplyYourself().

DumpInfo()

void G4ParticleHPChannelList::DumpInfo

(

)

Definition at line 179 of file G4ParticleHPChannelList.cc.

{
  G4cout << "================================================================" << G4endl;
  G4cout << " Element: " << theElement->GetName() << G4endl;
  G4cout << " Number of channels: " << nChannels << G4endl;
  G4cout << " Projectile: " << theProjectile->GetParticleName() << G4endl;
  G4cout << " Directory name: " << theDir << G4endl;
  for (G4int i = 0; i < nChannels; ++i) {
    if (theChannels[i]->HasDataInAnyFinalState()) {
      G4cout << "----------------------------------------------------------------" << G4endl;
      theChannels[i]->DumpInfo();
      G4cout << "----------------------------------------------------------------" << G4endl;
    }
  }
  G4cout << "================================================================" << G4endl;
}

GetNumberOfChannels()

G4int G4ParticleHPChannelList::GetNumberOfChannels ( )

inline

Definition at line 72 of file G4ParticleHPChannelList.hh.

{ return nChannels; }

GetXsec()

G4double G4ParticleHPChannelList::GetXsec ( G4double anEnergy )

inline

Definition at line 63 of file G4ParticleHPChannelList.hh.

    {
      G4double result = 0.0;
      for (G4int i = 0; i < nChannels; ++i) {
        result += std::max(0., theChannels[i]->GetXsec(anEnergy));
      }
      return result;
    }

Referenced by GetXsec().

HasDataInAnyFinalState()

G4bool G4ParticleHPChannelList::HasDataInAnyFinalState ( )

inline

Definition at line 74 of file G4ParticleHPChannelList.hh.

    {
      G4bool result = false;
      for (G4int i = 0; i < nChannels; ++i) {
        if (theChannels[i]->HasDataInAnyFinalState())
    {
          result = true;
      break;
    }
      }
      return result;
    }

Referenced by DumpInfo(), and HasDataInAnyFinalState().

Init() [1/2]

void G4ParticleHPChannelList::Init	(	G4Element *	anElement,
		const G4String &	dirName,
		G4ParticleDefinition *	projectile )

Definition at line 161 of file G4ParticleHPChannelList.cc.

{
  theDir = dirName;
  // G4cout << theDir << G4endl;
  theElement = anElement;
  // G4cout << theElement->GetName() << G4endl;
  theProjectile = projectile;
}

Init() [2/2]

void G4ParticleHPChannelList::Init

(

G4int

n

)

operator=()

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

delete

Register()

void G4ParticleHPChannelList::Register	(	G4ParticleHPFinalState *	theFS,
		const G4String &	aName )

Definition at line 171 of file G4ParticleHPChannelList.cc.

{
  theChannels[idx] = new G4ParticleHPChannel(theProjectile);
  theChannels[idx]->Init(theElement, theDir, aName);
  theChannels[idx]->Register(theFS);
  ++idx;
}

---

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

• G4ParticleHPChannelList.hh
• G4ParticleHPChannelList.cc