G4ParticleHPChannel Class Reference

#include <G4ParticleHPChannel.hh>

Public Member Functions
	G4ParticleHPChannel (G4ParticleDefinition *projectile)
	G4ParticleHPChannel ()
	~G4ParticleHPChannel ()
G4double	GetXsec (G4double energy)
G4double	GetWeightedXsec (G4double energy, G4int isoNumber)
G4double	GetFSCrossSection (G4double energy, G4int isoNumber)
G4bool	IsActive (G4int isoNumber)
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)
G4int	GetNiso ()
G4double	GetN (G4int i)
G4double	GetZ (G4int i)
G4double	GetM (G4int i)
G4bool	HasDataInAnyFinalState ()
void	DumpInfo ()
G4String	GetFSType () const
G4ParticleHPFinalState **	GetFinalStates () const

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

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

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

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

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

Referenced by G4ParticleHPChannelList::DumpInfo().

GetFinalStates()

G4ParticleHPFinalState ** G4ParticleHPChannel::GetFinalStates ( ) const

inline

Definition at line 177 of file G4ParticleHPChannel.hh.

                                                   {
    return theFinalStates; 
  }

GetFSCrossSection()

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

Definition at line 61 of file G4ParticleHPChannel.cc.

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

Referenced by G4ParticleHPChannelList::ApplyYourself().

GetFSType()

G4String G4ParticleHPChannel::GetFSType ( ) const

inline

Definition at line 173 of file G4ParticleHPChannel.hh.

                             {
    return theFSType;
  }

GetM()

G4double G4ParticleHPChannel::GetM ( G4int  i )

inline

Definition at line 158 of file G4ParticleHPChannel.hh.

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

Referenced by ApplyYourself().

GetN()

G4double G4ParticleHPChannel::GetN ( G4int  i )

inline

Definition at line 156 of file G4ParticleHPChannel.hh.

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

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

GetNiso()

G4int G4ParticleHPChannel::GetNiso ( )

inline

Definition at line 154 of file G4ParticleHPChannel.hh.

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

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

Referenced by G4FissLib::ApplyYourself().

GetZ()

G4double G4ParticleHPChannel::GetZ ( G4int  i )

inline

Definition at line 157 of file G4ParticleHPChannel.hh.

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

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

Referenced by UpdateData().

HasAnyData()

G4bool G4ParticleHPChannel::HasAnyData ( G4int  isoNumber )

inline

Definition at line 138 of file G4ParticleHPChannel.hh.

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

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

Referenced by Register().

HasFSData()

G4bool G4ParticleHPChannel::HasFSData ( G4int  isoNumber )

inline

Definition at line 136 of file G4ParticleHPChannel.hh.

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

Init() [1/2]

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

Definition at line 73 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 66 of file G4ParticleHPChannel.cc.

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

IsActive()

G4bool G4ParticleHPChannel::IsActive ( G4int  isoNumber )

inline

Definition at line 134 of file G4ParticleHPChannel.hh.

{ return active[isoNumber]; }

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

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

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

---

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

• G4ParticleHPChannel.hh
• G4ParticleHPChannel.cc