G4ParticleHPChannel Class Reference

#include <G4ParticleHPChannel.hh>

Collaboration diagram for G4ParticleHPChannel:

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

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

Detailed Description

Definition at line 57 of file G4ParticleHPChannel.hh.

Constructor & Destructor Documentation

G4ParticleHPChannel() [1/2]

G4ParticleHPChannel::G4ParticleHPChannel ( G4ParticleDefinition *  projectile )

inline

Definition at line 61 of file G4ParticleHPChannel.hh.

  : wendtFissionGenerator(getenv("G4NEUTRON_HP_USE_WENDT_FISSION_MODEL") == NULL ? NULL : G4WendtFissionFragmentGenerator::GetInstance())
  {
    theProjectile = const_cast<G4ParticleDefinition*>(projectile);
    theChannelData = new G4ParticleHPVector; 
    theBuffer = 0;
    theIsotopeWiseData = 0;
    theFinalStates = 0;
    active = 0;
    registerCount = -1;
  }

G4ParticleHPChannel() [2/2]

G4ParticleHPChannel::G4ParticleHPChannel ( )

inline

Definition at line 73 of file G4ParticleHPChannel.hh.

  : wendtFissionGenerator(getenv("G4NEUTRON_HP_USE_WENDT_FISSION_MODEL") == NULL ? NULL : G4WendtFissionFragmentGenerator::GetInstance())
  {
    theProjectile = G4Neutron::Neutron();
    theChannelData = new G4ParticleHPVector; 
    theBuffer = 0;
    theIsotopeWiseData = 0;
    theFinalStates = 0;
    active = 0;
    registerCount = -1;
  }

Here is the call graph for this function:

~G4ParticleHPChannel()

G4ParticleHPChannel::~G4ParticleHPChannel ( )

inline

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

Here is the call graph for this function:

Member Function Documentation

ApplyYourself()

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

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

Here is the call graph for this function:

Here is the caller graph for this function:

DumpInfo()

void G4ParticleHPChannel::DumpInfo

(

)

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

}

Here is the call graph for this function:

Here is the caller graph for this function:

GetFinalStates()

G4ParticleHPFinalState** G4ParticleHPChannel::GetFinalStates ( ) const

inline

Definition at line 167 of file G4ParticleHPChannel.hh.

                                                   {
    return theFinalStates; 
  }

GetFSCrossSection()

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

Definition at line 57 of file G4ParticleHPChannel.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

GetFSType()

G4String G4ParticleHPChannel::GetFSType ( ) const

inline

Definition at line 163 of file G4ParticleHPChannel.hh.

                             {
    return theFSType;
  }

GetM()

G4double G4ParticleHPChannel::GetM ( G4int  i )

inline

Definition at line 148 of file G4ParticleHPChannel.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

GetN()

G4double G4ParticleHPChannel::GetN ( G4int  i )

inline

Definition at line 146 of file G4ParticleHPChannel.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

GetNiso()

G4int G4ParticleHPChannel::GetNiso ( )

inline

Definition at line 144 of file G4ParticleHPChannel.hh.

{return niso;}

Here is the caller graph for this function:

GetWeightedXsec()

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

Definition at line 52 of file G4ParticleHPChannel.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

GetXsec()

G4double G4ParticleHPChannel::GetXsec

(

G4double

energy

)

Definition at line 47 of file G4ParticleHPChannel.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

GetZ()

G4double G4ParticleHPChannel::GetZ ( G4int  i )

inline

Definition at line 147 of file G4ParticleHPChannel.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

Harmonise()

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

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

Here is the call graph for this function:

Here is the caller graph for this function:

HasAnyData()

G4bool G4ParticleHPChannel::HasAnyData ( G4int  isoNumber )

inline

Definition at line 128 of file G4ParticleHPChannel.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

HasDataInAnyFinalState()

G4bool G4ParticleHPChannel::HasDataInAnyFinalState ( )

inline

Definition at line 150 of file G4ParticleHPChannel.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

HasFSData()

G4bool G4ParticleHPChannel::HasFSData ( G4int  isoNumber )

inline

Definition at line 126 of file G4ParticleHPChannel.hh.

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

Here is the call graph for this function:

Init() [1/2]

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

Definition at line 69 of file G4ParticleHPChannel.cc.

  {
    theDir = dirName;
    theElement = anElement;
  }

Here is the caller graph for this function:

Init() [2/2]

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

Definition at line 63 of file G4ParticleHPChannel.cc.

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

Here is the call graph for this function:

IsActive()

G4bool G4ParticleHPChannel::IsActive ( G4int  isoNumber )

inline

Definition at line 124 of file G4ParticleHPChannel.hh.

{ return active[isoNumber]; }

Register()

G4bool G4ParticleHPChannel::Register

(

G4ParticleHPFinalState *

theFS

)

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

Here is the call graph for this function:

Here is the caller graph for this function:

UpdateData() [1/2]

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

inline

Definition at line 137 of file G4ParticleHPChannel.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

UpdateData() [2/2]

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

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

Here is the call graph for this function:

Member Data Documentation

active

G4bool* G4ParticleHPChannel::active

private

Definition at line 179 of file G4ParticleHPChannel.hh.

niso

G4int G4ParticleHPChannel::niso

private

Definition at line 180 of file G4ParticleHPChannel.hh.

registerCount

G4int G4ParticleHPChannel::registerCount

private

Definition at line 188 of file G4ParticleHPChannel.hh.

theBuffer

G4ParticleHPVector* G4ParticleHPChannel::theBuffer

private

Definition at line 175 of file G4ParticleHPChannel.hh.

theChannelData

G4ParticleHPVector* G4ParticleHPChannel::theChannelData

private

Definition at line 174 of file G4ParticleHPChannel.hh.

theDir

G4String G4ParticleHPChannel::theDir

private

Definition at line 184 of file G4ParticleHPChannel.hh.

theElement

G4Element* G4ParticleHPChannel::theElement

private

Definition at line 186 of file G4ParticleHPChannel.hh.

theFinalStates

G4ParticleHPFinalState** G4ParticleHPChannel::theFinalStates

private

Definition at line 178 of file G4ParticleHPChannel.hh.

theFSType

G4String G4ParticleHPChannel::theFSType

private

Definition at line 185 of file G4ParticleHPChannel.hh.

theIsotopeWiseData

G4ParticleHPIsoData* G4ParticleHPChannel::theIsotopeWiseData

private

Definition at line 177 of file G4ParticleHPChannel.hh.

theProjectile

G4ParticleDefinition* G4ParticleHPChannel::theProjectile

private

Definition at line 172 of file G4ParticleHPChannel.hh.

theStableOnes

G4StableIsotopes G4ParticleHPChannel::theStableOnes

private

Definition at line 182 of file G4ParticleHPChannel.hh.

wendtFissionGenerator

G4WendtFissionFragmentGenerator* const G4ParticleHPChannel::wendtFissionGenerator

private

Definition at line 190 of file G4ParticleHPChannel.hh.

---

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

• G4ParticleHPChannel.hh
• G4ParticleHPChannel.cc