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

Constructor & Destructor Documentation

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;
    niso = -1;
    theElement = NULL;
  }

G4ParticleHPChannel::G4ParticleHPChannel ( )

inline

Definition at line 75 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;
    niso = -1;
    theElement = NULL;
  }

Here is the call graph for this function:

G4ParticleHPChannel::~G4ParticleHPChannel ( )

inline

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

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:

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:

G4ParticleHPFinalState** G4ParticleHPChannel::GetFinalStates ( ) const

inline

Definition at line 171 of file G4ParticleHPChannel.hh.

                                                   {
    return theFinalStates; 
  }

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:

G4String G4ParticleHPChannel::GetFSType ( ) const

inline

Definition at line 167 of file G4ParticleHPChannel.hh.

                             {
    return theFSType;
  }

G4double G4ParticleHPChannel::GetM ( G4int  i )

inline

Definition at line 152 of file G4ParticleHPChannel.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4ParticleHPChannel::GetN ( G4int  i )

inline

Definition at line 150 of file G4ParticleHPChannel.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4int G4ParticleHPChannel::GetNiso ( )

inline

Definition at line 148 of file G4ParticleHPChannel.hh.

{return niso;}

Here is the caller graph for this function:

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:

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:

G4double G4ParticleHPChannel::GetZ ( G4int  i )

inline

Definition at line 151 of file G4ParticleHPChannel.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

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:

G4bool G4ParticleHPChannel::HasAnyData ( G4int  isoNumber )

inline

Definition at line 132 of file G4ParticleHPChannel.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool G4ParticleHPChannel::HasDataInAnyFinalState ( )

inline

Definition at line 154 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:

G4bool G4ParticleHPChannel::HasFSData ( G4int  isoNumber )

inline

Definition at line 130 of file G4ParticleHPChannel.hh.

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

Here is the call graph for this function:

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:

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:

G4bool G4ParticleHPChannel::IsActive ( G4int  isoNumber )

inline

Definition at line 128 of file G4ParticleHPChannel.hh.

{ return active[isoNumber]; }

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:

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

inline

Definition at line 141 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:

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:

---

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

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