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)
G4HadFinalState *	ApplyYourself (G4int, G4int, G4int, const G4HadProjectile &aTrack)
void	Init (G4Element *anElement, const G4String &dirName, G4ParticleDefinition *projectile)
void	Register (G4ParticleHPFinalState *theFS, const G4String &aName)
G4double	GetXsec (G4double anEnergy)
G4double	GetWeightedXsec (G4double anEnergy, G4int isotopeJ)
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;
}

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

~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() [1/2]

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().

ApplyYourself() [2/2]

G4HadFinalState * G4ParticleHPChannelList::ApplyYourself	(	G4int	isotope,
		G4int	anZ,
		G4int	anA,
		const G4HadProjectile &	aTrack )

Definition at line 161 of file G4ParticleHPChannelList.cc.

                                                                                                                             {
  // do not choose Z and A again and go to selection of channel
  G4ParticleHPThermalBoost aThermalE;
  G4int i;
  G4double random;
  // decide on the channel
  G4double* running = new G4double[nChannels];
  running[0] = 0.0;
  // targA and targZ does not set to -1
  G4int targA = anA;
  G4int targZ = anZ;
  G4double energy = aThermalE.GetThermalEnergy( aTrack, targA, targZ, aTrack.GetMaterial()->GetTemperature() );
  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( energy, isotope );
    }
  }
  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" );
    }
    G4cout << "Warning from NeutronHP: could not find proper reaction channel. "
           << "This may be caused 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 );
    G4DynamicParticle* targ_dp = new G4DynamicParticle( targ_pd, G4ThreeVector(1.0, 0.0, 0.0), 0.0 );
    unChanged.SetEnergyChange( aTrack.GetKineticEnergy() );
    unChanged.SetMomentumChange( aTrack.Get4Momentum().vect().unit() );
    unChanged.AddSecondary( targ_dp );
    G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->SetTargA( targA );
    G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->SetTargZ( targZ );
    delete [] running;
    return &unChanged;
  }
  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 ( G4FindDataDir( "G4ParticleHPDebug" ) != nullptr ) G4cout << " G4ParticleHPChannelList SELECTED ISOTOPE " << isotope 
                                << " SELECTED CHANNEL " << lChan << G4endl;
  #endif
  return theChannels[lChan]->ApplyYourself( aTrack, isotope );
}

DumpInfo()

void G4ParticleHPChannelList::DumpInfo

(

)

Definition at line 233 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 87 of file G4ParticleHPChannelList.hh.

{ return nChannels; }

GetWeightedXsec()

G4double G4ParticleHPChannelList::GetWeightedXsec ( G4double anEnergy, G4int isotopeJ )

inline

Definition at line 76 of file G4ParticleHPChannelList.hh.

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

Referenced by GetWeightedXsec().

GetXsec()

G4double G4ParticleHPChannelList::GetXsec ( G4double anEnergy )

inline

Definition at line 66 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 89 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 215 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 225 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