G4ParticleHPChannel Class Reference

#include <G4ParticleHPChannel.hh>

Public Member Functions
G4HadFinalState *	ApplyYourself (const G4HadProjectile &theTrack, G4int isoNumber=-1)
void	DumpInfo ()
	G4ParticleHPChannel ()
	G4ParticleHPChannel (G4ParticleDefinition *projectile)
G4ParticleHPFinalState **	GetFinalStates () const
G4double	GetFSCrossSection (G4double energy, G4int isoNumber)
G4String	GetFSType () const
G4double	GetM (G4int i)
G4double	GetN (G4int i)
G4int	GetNiso ()
G4double	GetWeightedXsec (G4double energy, G4int isoNumber)
G4double	GetXsec (G4double energy)
G4double	GetZ (G4int i)
void	Harmonise (G4ParticleHPVector *&theStore, G4ParticleHPVector *theNew)
G4bool	HasAnyData (G4int isoNumber)
G4bool	HasDataInAnyFinalState ()
G4bool	HasFSData (G4int isoNumber)
void	Init (G4Element *theElement, const G4String dirName)
void	Init (G4Element *theElement, const G4String dirName, const G4String fsType)
G4bool	IsActive (G4int isoNumber)
G4bool	Register (G4ParticleHPFinalState *theFS)
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)
	~G4ParticleHPChannel ()

Private Attributes
G4bool *	active
G4int	niso
G4int	registerCount
G4ParticleHPVector *	theBuffer
G4ParticleHPVector *	theChannelData
G4String	theDir
G4Element *	theElement
G4ParticleHPFinalState **	theFinalStates
G4String	theFSType
G4ParticleHPIsoData *	theIsotopeWiseData
G4ParticleDefinition *	theProjectile
G4StableIsotopes	theStableOnes
G4WendtFissionFragmentGenerator *const	wendtFissionGenerator

Detailed Description

Definition at line 55 of file G4ParticleHPChannel.hh.

Constructor & Destructor Documentation

G4ParticleHPChannel() [1/2]

G4ParticleHPChannel::G4ParticleHPChannel ( G4ParticleDefinition *  projectile )

inline

Definition at line 59 of file G4ParticleHPChannel.hh.

                                                        :
    wendtFissionGenerator( G4ParticleHPManager::GetInstance()->GetUseWendtFissionModel() ?
               G4WendtFissionFragmentGenerator::GetInstance() : nullptr )
  {
    if ( G4ParticleHPManager::GetInstance()->GetUseWendtFissionModel() ) {
      // Make sure both fission fragment models are not active at same time
      G4ParticleHPManager::GetInstance()->SetProduceFissionFragments( false );
    }
    theProjectile = const_cast<G4ParticleDefinition*>(projectile);
    theChannelData = new G4ParticleHPVector; 
    theBuffer = 0;
    theIsotopeWiseData = 0;
    theFinalStates = 0;
    active = 0;
    registerCount = -1;
    niso = -1;
    theElement = NULL;
  }

References active, G4ParticleHPManager::GetInstance(), niso, registerCount, G4ParticleHPManager::SetProduceFissionFragments(), theBuffer, theChannelData, theElement, theFinalStates, theIsotopeWiseData, and theProjectile.

G4ParticleHPChannel() [2/2]

G4ParticleHPChannel::G4ParticleHPChannel ( )

inline

Definition at line 78 of file G4ParticleHPChannel.hh.

                        : wendtFissionGenerator( G4ParticleHPManager::GetInstance()->GetUseWendtFissionModel() ?
                         G4WendtFissionFragmentGenerator::GetInstance() : nullptr )
  {
    if ( G4ParticleHPManager::GetInstance()->GetUseWendtFissionModel() ) {
      // Make sure both fission fragment models are not active at same time
      G4ParticleHPManager::GetInstance()->SetProduceFissionFragments( false );
    }
    theProjectile = G4Neutron::Neutron();
    theChannelData = new G4ParticleHPVector; 
    theBuffer = 0;
    theIsotopeWiseData = 0;
    theFinalStates = 0;
    active = 0;
    registerCount = -1;
    niso = -1;
    theElement = NULL;
  }

References active, G4ParticleHPManager::GetInstance(), G4Neutron::Neutron(), niso, registerCount, G4ParticleHPManager::SetProduceFissionFragments(), theBuffer, theChannelData, theElement, theFinalStates, theIsotopeWiseData, and theProjectile.

~G4ParticleHPChannel()

G4ParticleHPChannel::~G4ParticleHPChannel ( )

inline

Definition at line 96 of file G4ParticleHPChannel.hh.

  {
    delete theChannelData; 
    // Following statement disabled to avoid SEGV
    // theBuffer is also deleted as "theChannelData" in
    // ~G4ParticleHPIsoData.  FWJ 06-Jul-1999
    //if(theBuffer != 0) delete theBuffer; 
    if(theIsotopeWiseData != 0) delete [] theIsotopeWiseData;
    // Deletion of FinalStates disabled to avoid endless looping
    // in the destructor heirarchy.  FWJ 06-Jul-1999
    //if(theFinalStates != 0)
    //{
    //  for(i=0; i<niso; i++)
    //  {
    //    delete theFinalStates[i];
    //  }
    //  delete [] theFinalStates;
    //}
    // FWJ experiment
    //if(active!=0) delete [] active;
    // T.K. 
    if ( theFinalStates != 0 )
    {
      for ( G4int i = 0 ; i < niso ; i++ )
      {
         delete theFinalStates[i];
      }
      delete [] theFinalStates;
    }
    if ( active != 0 ) delete [] active;
  }

References active, niso, theChannelData, theFinalStates, and theIsotopeWiseData.

Member Function Documentation

ApplyYourself()

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

Definition at line 228 of file G4ParticleHPChannel.cc.

  {
//    G4cout << "G4ParticleHPChannel::ApplyYourself+"<<niso<<G4endl;
    if ( anIsotope != -1 && anIsotope != -2 ) 
    {
       //Inelastic Case
       //G4cout << "G4ParticleHPChannel Inelastic Case" 
       //<< " Z= " << this->GetZ(it) << " A = " << this->GetN(it) << G4endl;
       G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->SetTargA( (G4int)this->GetN(anIsotope) ); 
       G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->SetTargZ( (G4int)this->GetZ(anIsotope) ); 
       return theFinalStates[anIsotope]->ApplyYourself(theTrack);
    }
    G4double sum=0;
    G4int it=0;
    G4double * xsec = new G4double[niso];
    G4ParticleHPThermalBoost aThermalE;
    for (G4int i=0; i<niso; i++)
    {
      if(theFinalStates[i]->HasAnyData())
      {
        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) 
    {
//      G4cout << "G4ParticleHPChannel::ApplyYourself theFinalState->Initialize+"<<G4endl;
//      G4cout << "G4ParticleHPChannel::ApplyYourself theFinalState->Initialize-"<<G4endl;
      it = static_cast<G4int>(niso*G4UniformRand());
    }
    else
    {
//      G4cout << "Are we still here? "<<sum<<G4endl;
//      G4cout << "TESTHP 23 NISO="<<niso<<G4endl;
      G4double random = G4UniformRand();
      G4double running=0;
//      G4cout << "G4ParticleHPChannel::ApplyYourself Done the sum"<<niso<<G4endl;
//      G4cout << "TESTHP 24 NISO="<<niso<<G4endl;
      for (G4int ix=0; ix<niso; ix++)
      {
        running += xsec[ix];
        //if(random<=running/sum) 
        if( sum == 0 || random <= running/sum ) 
        { 
          it = ix;
      break;
        }
      }
      if(it==niso) it--;
    }
    delete [] xsec;
    G4HadFinalState * theFinalState=0;
    const G4int A = (G4int)this->GetN(it);
    const G4int Z = (G4int)this->GetZ(it);
    const G4int M = (G4int)this->GetM(it);
 
                                       //-2:Marker for Fission
    if(wendtFissionGenerator&&anIsotope==-2)
    {
        theFinalState = wendtFissionGenerator->ApplyYourself(theTrack, Z, A);
    }
 
    // Use the standard procedure if the G4FissionFragmentGenerator model fails
    if (!theFinalState)
    {
 
       G4int icounter=0;
       G4int icounter_max=1024;
       while(theFinalState==0) // 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;
          }
//        G4cout << "TESTHP 24 it="<<it<<G4endl;
          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;
    G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->SetTargA( A );
    G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->SetTargZ( Z );
    G4ParticleHPManager::GetInstance()->GetReactionWhiteBoard()->SetTargM( M );
 
    return theFinalState;
  }

References A, G4ParticleHPFinalState::ApplyYourself(), G4WendtFissionFragmentGenerator::ApplyYourself(), G4cout, G4endl, G4UniformRand, G4ParticleHPManager::GetInstance(), GetM(), G4HadProjectile::GetMaterial(), G4ParticleHPFinalState::GetN(), GetN(), G4ParticleHPManager::GetReactionWhiteBoard(), G4Material::GetTemperature(), G4ParticleHPThermalBoost::GetThermalEnergy(), G4ParticleHPIsoData::GetXsec(), G4ParticleHPFinalState::GetZ(), GetZ(), HasAnyData(), M, niso, G4ParticleHPReactionWhiteBoard::SetTargA(), G4ParticleHPReactionWhiteBoard::SetTargM(), G4ParticleHPReactionWhiteBoard::SetTargZ(), theFinalStates, theIsotopeWiseData, wendtFissionGenerator, and Z.

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

DumpInfo()

void G4ParticleHPChannel::DumpInfo

(

)

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){
    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;
    }
  }
 
}

References eV, G4cout, G4endl, G4ParticleHPVector::GetEnergy(), G4Element::GetName(), G4ParticleHPVector::GetVectorLength(), G4ParticleHPVector::GetXsec(), G4ParticleHPFinalState::HasXsec(), niso, theChannelData, theDir, theElement, theFinalStates, and theFSType.

Referenced by G4ParticleHPChannelList::DumpInfo().

GetFinalStates()

G4ParticleHPFinalState ** G4ParticleHPChannel::GetFinalStates ( ) const

inline

Definition at line 177 of file G4ParticleHPChannel.hh.

                                                   {
    return theFinalStates; 
  }

References theFinalStates.

GetFSCrossSection()

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

Definition at line 61 of file G4ParticleHPChannel.cc.

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

References G4INCL::KinematicsUtils::energy(), G4ParticleHPFinalState::GetXsec(), and theFinalStates.

Referenced by G4ParticleHPChannelList::ApplyYourself().

GetFSType()

G4String G4ParticleHPChannel::GetFSType ( ) const

inline

Definition at line 173 of file G4ParticleHPChannel.hh.

                             {
    return theFSType;
  }

References theFSType.

GetM()

G4double G4ParticleHPChannel::GetM ( G4int  i )

inline

Definition at line 158 of file G4ParticleHPChannel.hh.

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

References G4ParticleHPFinalState::GetM(), and theFinalStates.

Referenced by ApplyYourself().

GetN()

G4double G4ParticleHPChannel::GetN ( G4int  i )

inline

Definition at line 156 of file G4ParticleHPChannel.hh.

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

References G4ParticleHPFinalState::GetN(), and theFinalStates.

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

GetNiso()

G4int G4ParticleHPChannel::GetNiso ( )

inline

Definition at line 154 of file G4ParticleHPChannel.hh.

{return niso;}

References niso.

Referenced by G4ParticleHPChannelList::ApplyYourself().

GetWeightedXsec()

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

Definition at line 56 of file G4ParticleHPChannel.cc.

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

References G4INCL::KinematicsUtils::energy(), G4ParticleHPIsoData::GetXsec(), and theIsotopeWiseData.

Referenced by G4ParticleHPChannelList::ApplyYourself().

GetXsec()

G4double G4ParticleHPChannel::GetXsec

(

G4double

energy

)

Definition at line 51 of file G4ParticleHPChannel.cc.

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

References G4INCL::KinematicsUtils::energy(), G4ParticleHPVector::GetXsec(), G4INCL::Math::max(), and theChannelData.

Referenced by G4FissLib::ApplyYourself().

GetZ()

G4double G4ParticleHPChannel::GetZ ( G4int  i )

inline

Definition at line 157 of file G4ParticleHPChannel.hh.

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

References G4ParticleHPFinalState::GetZ(), and theFinalStates.

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

Harmonise()

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

Definition at line 184 of file G4ParticleHPChannel.cc.

  {
    G4int s_tmp = 0, n=0, m_tmp=0;
    G4ParticleHPVector * 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;
  }

References G4ParticleHPVector::GetEnergy(), G4ParticleHPVector::GetVectorLength(), G4ParticleHPVector::GetXsec(), G4INCL::Math::max(), CLHEP::detail::n, and G4ParticleHPVector::SetData().

Referenced by UpdateData().

HasAnyData()

G4bool G4ParticleHPChannel::HasAnyData ( G4int  isoNumber )

inline

Definition at line 138 of file G4ParticleHPChannel.hh.

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

References G4ParticleHPFinalState::HasAnyData(), and theFinalStates.

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

HasDataInAnyFinalState()

G4bool G4ParticleHPChannel::HasDataInAnyFinalState ( )

inline

Definition at line 160 of file G4ParticleHPChannel.hh.

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

References HasAnyData(), niso, and theFinalStates.

Referenced by Register().

HasFSData()

G4bool G4ParticleHPChannel::HasFSData ( G4int  isoNumber )

inline

Definition at line 136 of file G4ParticleHPChannel.hh.

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

References G4ParticleHPFinalState::HasFSData(), and theFinalStates.

Init() [1/2]

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

Definition at line 73 of file G4ParticleHPChannel.cc.

  {
    theDir = dirName;
    theElement = anElement;
  }

References theDir, and theElement.

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 66 of file G4ParticleHPChannel.cc.

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

References Init(), and theFSType.

IsActive()

G4bool G4ParticleHPChannel::IsActive ( G4int  isoNumber )

inline

Definition at line 134 of file G4ParticleHPChannel.hh.

{ return active[isoNumber]; }

References active.

Register()

G4bool G4ParticleHPChannel::Register

(

G4ParticleHPFinalState *

theFS

)

Definition at line 79 of file G4ParticleHPChannel.cc.

  {
    registerCount++;
    G4int Z = G4lrint(theElement->GetZ());
 
    Z = Z-registerCount;
    if ( registerCount > 5 ) throw G4HadronicException(__FILE__, __LINE__, "Channel: Do not know what to do with this material"); // for Elastic, Capture, Fission case 
    if ( Z < 1 ) return false; 
/*
    if(registerCount<5)
    {
      Z = Z-registerCount;
    }
*/
    //if(Z=theElement->GetZ()-5) throw G4HadronicException(__FILE__, __LINE__, "Channel: Do not know what to do with this material");
    // Bug fix by TK on behalf of AH
    //if ( Z <=theElement->GetZ()-5 ) throw G4HadronicException(__FILE__, __LINE__, "Channel: Do not know what to do with this material");
    G4int count = 0;
    if(registerCount==0) count = theElement->GetNumberOfIsotopes();
    if(count == 0||registerCount!=0) count +=
         theStableOnes.GetNumberOfIsotopes(Z);
    niso = count;
    delete [] theIsotopeWiseData;
    theIsotopeWiseData = new G4ParticleHPIsoData [niso];
    delete [] active;
    active = new G4bool[niso];
 
    delete [] theFinalStates;
    theFinalStates = new G4ParticleHPFinalState * [niso];
    delete theChannelData;
    theChannelData = new G4ParticleHPVector; 
    for(G4int i=0; i<niso; i++)
    {
      theFinalStates[i] = theFS->New();
      theFinalStates[i]->SetProjectile(theProjectile);
    }
    count = 0;
    G4int nIsos = niso;
    if(theElement->GetNumberOfIsotopes()!=0&&registerCount==0)
    {
      for (G4int i1=0; i1<nIsos; i1++)
      {
        // G4cout <<" Init: normal case"<<G4endl;
        G4int A = theElement->GetIsotope(i1)->GetN();
        G4int M = theElement->GetIsotope(i1)->Getm();
        G4double frac = theElement->GetRelativeAbundanceVector()[i1]/perCent;
        //theFinalStates[i1]->SetA_Z(A, Z);
    //UpdateData(A, Z, count++, frac);
        theFinalStates[i1]->SetA_Z(A, Z, M);
    UpdateData(A, Z, M, count++, frac, theProjectile);
      }
    } else {
      //G4cout <<" Init: mean case: "
      //       <<theStableOnes.GetNumberOfIsotopes(Z)<<" "
    //     <<Z<<" "<<theElement
    //     << G4endl;
      G4int first = theStableOnes.GetFirstIsotope(Z);
      for(G4int i1=0; 
        i1<theStableOnes.GetNumberOfIsotopes(Z);
        i1++)
      {
        G4int A = theStableOnes.GetIsotopeNucleonCount(first+i1);
        G4double frac = theStableOnes.GetAbundance(first+i1);
        theFinalStates[i1]->SetA_Z(A, Z);
        UpdateData(A, Z, count++, frac, theProjectile);
      }
    }
    G4bool result = HasDataInAnyFinalState();
 
    //To avoid issuing hash by worker threads
    if ( result ) theChannelData->Hash();
 
    return result;
  }

References A, active, G4lrint(), G4StableIsotopes::GetAbundance(), G4StableIsotopes::GetFirstIsotope(), G4Element::GetIsotope(), G4StableIsotopes::GetIsotopeNucleonCount(), G4Isotope::Getm(), G4Isotope::GetN(), G4Element::GetNumberOfIsotopes(), G4StableIsotopes::GetNumberOfIsotopes(), G4Element::GetRelativeAbundanceVector(), G4Element::GetZ(), HasDataInAnyFinalState(), G4ParticleHPVector::Hash(), M, G4ParticleHPFinalState::New(), niso, perCent, registerCount, G4ParticleHPFinalState::SetA_Z(), G4ParticleHPFinalState::SetProjectile(), theChannelData, theElement, theFinalStates, theIsotopeWiseData, theProjectile, theStableOnes, UpdateData(), and Z.

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 147 of file G4ParticleHPChannel.hh.

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

References A, M, UpdateData(), and Z.

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 155 of file G4ParticleHPChannel.cc.

  {
    // Initialze the G4FissionFragment generator for this isomer if needed
    if(wendtFissionGenerator)
    {
      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 = 0;
    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, abundance, theDir, tString);
      active[index] = theIsotopeWiseData[index].Init(A, Z, M, abundance, theDir, tString);
      if(active[index]) theBuffer = theIsotopeWiseData[index].MakeChannelData();
    }
    if(theBuffer != 0) Harmonise(theChannelData, theBuffer);
  }

References A, active, G4ParticleHPIsoData::FillChannelData(), G4ParticleHPFinalState::GetXsec(), Harmonise(), HasAnyData(), G4ParticleHPFinalState::Init(), G4ParticleHPIsoData::Init(), G4WendtFissionFragmentGenerator::InitializeANucleus(), M, G4ParticleHPIsoData::MakeChannelData(), theBuffer, theChannelData, theDir, theFinalStates, theFSType, theIsotopeWiseData, G4ParticleHPVector::Times(), wendtFissionGenerator, and Z.

Field Documentation

active

G4bool* G4ParticleHPChannel::active

private

Definition at line 189 of file G4ParticleHPChannel.hh.

Referenced by G4ParticleHPChannel(), IsActive(), Register(), UpdateData(), and ~G4ParticleHPChannel().

niso

G4int G4ParticleHPChannel::niso

private

Definition at line 190 of file G4ParticleHPChannel.hh.

Referenced by ApplyYourself(), DumpInfo(), G4ParticleHPChannel(), GetNiso(), HasDataInAnyFinalState(), Register(), and ~G4ParticleHPChannel().

registerCount

G4int G4ParticleHPChannel::registerCount

private

Definition at line 198 of file G4ParticleHPChannel.hh.

Referenced by G4ParticleHPChannel(), and Register().

theBuffer

G4ParticleHPVector* G4ParticleHPChannel::theBuffer

private

Definition at line 185 of file G4ParticleHPChannel.hh.

Referenced by G4ParticleHPChannel(), and UpdateData().

theChannelData

G4ParticleHPVector* G4ParticleHPChannel::theChannelData

private

Definition at line 184 of file G4ParticleHPChannel.hh.

Referenced by DumpInfo(), G4ParticleHPChannel(), GetXsec(), Register(), UpdateData(), and ~G4ParticleHPChannel().

theDir

G4String G4ParticleHPChannel::theDir

private

Definition at line 194 of file G4ParticleHPChannel.hh.

Referenced by DumpInfo(), Init(), and UpdateData().

theElement

G4Element* G4ParticleHPChannel::theElement

private

Definition at line 196 of file G4ParticleHPChannel.hh.

Referenced by DumpInfo(), G4ParticleHPChannel(), Init(), and Register().

theFinalStates

G4ParticleHPFinalState** G4ParticleHPChannel::theFinalStates

private

Definition at line 188 of file G4ParticleHPChannel.hh.

Referenced by ApplyYourself(), DumpInfo(), G4ParticleHPChannel(), GetFinalStates(), GetFSCrossSection(), GetM(), GetN(), GetZ(), HasAnyData(), HasDataInAnyFinalState(), HasFSData(), Register(), UpdateData(), and ~G4ParticleHPChannel().

theFSType

G4String G4ParticleHPChannel::theFSType

private

Definition at line 195 of file G4ParticleHPChannel.hh.

Referenced by DumpInfo(), GetFSType(), Init(), and UpdateData().

theIsotopeWiseData

G4ParticleHPIsoData* G4ParticleHPChannel::theIsotopeWiseData

private

Definition at line 187 of file G4ParticleHPChannel.hh.

Referenced by ApplyYourself(), G4ParticleHPChannel(), GetWeightedXsec(), Register(), UpdateData(), and ~G4ParticleHPChannel().

theProjectile

G4ParticleDefinition* G4ParticleHPChannel::theProjectile

private

Definition at line 182 of file G4ParticleHPChannel.hh.

Referenced by G4ParticleHPChannel(), and Register().

theStableOnes

G4StableIsotopes G4ParticleHPChannel::theStableOnes

private

Definition at line 192 of file G4ParticleHPChannel.hh.

Referenced by Register().

wendtFissionGenerator

G4WendtFissionFragmentGenerator* const G4ParticleHPChannel::wendtFissionGenerator

private

Definition at line 200 of file G4ParticleHPChannel.hh.

Referenced by ApplyYourself(), and UpdateData().

---

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

• source/processes/hadronic/models/particle_hp/include/G4ParticleHPChannel.hh
• source/processes/hadronic/models/particle_hp/src/G4ParticleHPChannel.cc