G4NeutronHPDeExGammas Class Reference

#include <G4NeutronHPDeExGammas.hh>

Public Member Functions
	G4NeutronHPDeExGammas ()
	~G4NeutronHPDeExGammas ()
void	Init (std::istream &aDataFile)
G4ReactionProductVector *	GetDecayGammas (G4int aLevel)
G4NeutronHPLevel *	GetLevel (G4int i)
G4int	GetNumberOfLevels ()
G4double	GetLevelEnergy (G4int aLevel)

Detailed Description

Definition at line 42 of file G4NeutronHPDeExGammas.hh.

Constructor & Destructor Documentation

G4NeutronHPDeExGammas::G4NeutronHPDeExGammas ( )

inline

Definition at line 46 of file G4NeutronHPDeExGammas.hh.

  {
    levelStart = 0;
    levelSize = 0;
    nLevels = 0;
    theLevels = 0;
  }

G4NeutronHPDeExGammas::~G4NeutronHPDeExGammas ( )

inline

Definition at line 53 of file G4NeutronHPDeExGammas.hh.

  {
    if(levelStart!=0) delete [] levelStart;
    if(levelSize!=0) delete [] levelSize;
    if(theLevels!=0) delete [] theLevels;
  }

Member Function Documentation

G4ReactionProductVector* G4NeutronHPDeExGammas::GetDecayGammas ( G4int  aLevel )

inline

Definition at line 62 of file G4NeutronHPDeExGammas.hh.

References G4UniformRand, G4ReactionProduct::GetTotalMomentum(), and G4ReactionProduct::SetMomentum().

Referenced by G4NeutronHPInelasticBaseFS::BaseApply(), and G4NeutronHPInelasticCompFS::CompositeApply().

  {
    if(aLevel>nLevels-1 || aLevel<0) return 0;
    if(nLevels==0) return new G4ReactionProductVector();
    G4ReactionProductVector * result = new G4ReactionProductVector;
    G4DynamicParticleVector * theResult;

    theResult = theLevels[aLevel]. GetDecayGammas();
    G4ReactionProduct * theCurrent;
    unsigned int i;
    for(i=0; i<theResult->size(); i++)
    {
      theCurrent = new G4ReactionProduct;
      *theCurrent = *(theResult->operator[](i));
      delete theResult->operator[](i);
      G4double costheta = 2.*G4UniformRand()-1;
      G4double theta = std::acos(costheta);
      G4double phi = CLHEP::twopi*G4UniformRand();
      G4double sinth = std::sin(theta);
      G4double en = theCurrent->GetTotalMomentum();
      G4ThreeVector temp(en*sinth*std::cos(phi), en*sinth*std::sin(phi), en*costheta );
      theCurrent->SetMomentum( temp ) ;
      result->push_back(theCurrent);
    }
    delete theResult;
    return result;
  }

G4NeutronHPLevel* G4NeutronHPDeExGammas::GetLevel ( G4int  i )

inline

Definition at line 90 of file G4NeutronHPDeExGammas.hh.

Referenced by G4NeutronHPInelasticCompFS::CompositeApply().

  {
    if(i>nLevels-1) return 0;
    return theLevels+i;
  }

G4double G4NeutronHPDeExGammas::GetLevelEnergy ( G4int  aLevel )

inline

Definition at line 98 of file G4NeutronHPDeExGammas.hh.

References G4NeutronHPLevel::GetLevelEnergy().

Referenced by G4NeutronHPInelasticBaseFS::BaseApply(), G4NeutronHPInelasticCompFS::CompositeApply(), and Init().

  {
    if(aLevel>nLevels-1 || aLevel<0) return 0;
    G4double result = theLevels[aLevel].GetLevelEnergy();
    return result;
  }

G4int G4NeutronHPDeExGammas::GetNumberOfLevels ( )

inline

Definition at line 96 of file G4NeutronHPDeExGammas.hh.

Referenced by G4NeutronHPInelasticBaseFS::BaseApply(), and G4NeutronHPInelasticCompFS::CompositeApply().

{ return nLevels; }

void G4NeutronHPDeExGammas::Init

(

std::istream &

aDataFile

)

Definition at line 35 of file G4NeutronHPDeExGammas.cc.

References buffer, G4NeutronHPGamma::GetGammaEnergy(), G4NeutronHPLevel::GetLevelEnergy(), G4NeutronHPGamma::GetLevelEnergy(), GetLevelEnergy(), G4NeutronHPGamma::Init(), python.hepunit::keV, G4INCL::Math::min(), G4NeutronHPLevel::SetGamma(), G4NeutronHPGamma::SetNext(), and G4NeutronHPLevel::SetNumberOfGammas().

Referenced by G4NeutronHPInelasticBaseFS::InitGammas(), and G4NeutronHPInelasticCompFS::InitGammas().

{
  G4NeutronHPGamma ** theGammas = new G4NeutronHPGamma * [50];
  G4int nGammas = 0;
  G4int nBuff = 50;
  for(;;)
  {
    G4NeutronHPGamma * theNew = new G4NeutronHPGamma;
    if(!theNew->Init(aDataFile))
    {
      delete theNew;
      break;
    }
    else
    {
      if(nGammas==nBuff)
      {
        nBuff+=50;
        G4NeutronHPGamma ** buffer = new G4NeutronHPGamma * [nBuff];
        for(G4int i=0;i<nGammas;i++) buffer[i] = theGammas[i];
        delete [] theGammas;
        theGammas = buffer;
      }
      theGammas[nGammas] = theNew;
      nGammas++;
    }
  }
  // all gammas are in. Now sort them into levels.

  // count the levels

  G4double currentE = 0;
  G4double nextE = 0;
  G4int i;
  G4double epsilon = 0.01*keV;
  for(i=0; i<nGammas; i++)
  {
    nextE = theGammas[i]->GetLevelEnergy();
    if(std::abs(currentE-nextE)>epsilon) nLevels++;
    currentE = nextE;
  }

  // Build the levels

  theLevels = new G4NeutronHPLevel[nLevels];
  levelStart = new G4int [nLevels];
  levelSize = new G4int [nLevels];

  // fill the levels

  currentE = 0;
  nextE = 0;
  G4int levelCounter=-1;
  for(i=0; i<nGammas; i++)
  {
    nextE = theGammas[i]->GetLevelEnergy();
    if(std::abs(currentE-nextE)>epsilon) 
    {
      levelCounter++;
      levelStart[levelCounter] = i;
      levelSize[levelCounter] = 0;
    }
    levelSize[levelCounter]++;
    currentE = nextE;
  }

  for(i=0; i<nLevels; i++)
  {
    theLevels[i].SetNumberOfGammas(levelSize[i]);
    for(G4int ii=levelStart[i]; ii<levelStart[i]+levelSize[i]; ii++)
    {
      theLevels[i].SetGamma(ii-levelStart[i], theGammas[ii]);
    }
  }

// set the next relation in the gammas.
  G4double levelE, gammaE, currentLevelE;
  G4double min; 
  for(i=0; i<nGammas; i++)
  {
    G4int it=-1;
    gammaE = theGammas[i]->GetGammaEnergy();
    currentLevelE = theGammas[i]->GetLevelEnergy();
    min = currentLevelE-gammaE-epsilon;
    for(G4int ii=0; ii<nLevels; ii++)
    {
      levelE = theLevels[ii].GetLevelEnergy();
      if(std::abs(currentLevelE-(levelE+gammaE))<min)
      {
        min = std::abs(currentLevelE-(levelE+gammaE));
        it = ii;
      }
    }
//080728
    if ( it != -1 && currentLevelE == theLevels[it].GetLevelEnergy() )
    {
       //TK Comment; Some data file in /Inelastic/Gammas has inconsistent level data (no level to transit)
       //G4cout << "DeExGammas Transition level error: it " << it << " " << currentLevelE << " " << gammaE << " " << theLevels[it-1].GetLevelEnergy() << " " << currentLevelE - theLevels[it-1].GetLevelEnergy() << G4endl;
       // Forced to connect the next(previous) level 
       it +=-1;
    }
//080728
    if(it!=-1) theGammas[i]->SetNext(&theLevels[it]);
  }
  // some garbage collection

  delete [] theGammas;

  // and we are Done.
}

---

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

• G4NeutronHPDeExGammas.hh
• G4NeutronHPDeExGammas.cc