#include <G4ParticleHPPhotonDist.hh>

Public Member Functions
	G4ParticleHPPhotonDist ()

G4double	GetLevelEnergy ()

G4ReactionProductVector *	GetPhotons (G4double anEnergy)

G4double	GetTargetMass ()

void	InitAngular (std::istream &aDataFile)

void	InitEnergies (std::istream &aDataFile)

G4bool	InitMean (std::istream &aDataFile)

void	InitPartials (std::istream &aDataFile, G4ParticleHPVector *theXsec=0)

G4bool	NeedsCascade ()

	~G4ParticleHPPhotonDist ()

Private Attributes
G4Cache< std::vector< G4int > * >	actualMult

G4int *	distribution

G4int *	disType

G4double *	energy

G4int	isoFlag

G4int *	isPrimary

G4int	nDiscrete

G4int	nDiscrete2

G4int	nGammaEnergies

G4int	nIso

G4int *	nNeu

G4int	nPartials

G4ParticleHPPartial **	partials

G4ParticleHPVector *	probs

G4int	repFlag

G4int	tabulationType

G4double	targetMass

G4ParticleHPAngularP **	theAngular

G4double	theBaseEnergy

G4double *	theGammas

G4ParticleHPInterpolator	theInt

G4int	theInternalConversionFlag

G4ParticleHPFastLegendre	theLegend

G4ParticleHPLegendreTable **	theLegendre

G4InterpolationManager	theLegendreManager

G4double *	theLevelEnergies

G4ParticleHPVector *	thePartialXsec

G4double *	thePhotonTransitionFraction

G4ParticleHPVector *	theReactionXsec

G4double *	theShells

G4ParticleHPVector	theTotalXsec

G4double *	theTransitionProbabilities

G4ParticleHPVector *	theYield

Detailed Description

Definition at line 57 of file G4ParticleHPPhotonDist.hh.

Constructor & Destructor Documentation

◆ G4ParticleHPPhotonDist()

G4ParticleHPPhotonDist::G4ParticleHPPhotonDist ( )

inline

Definition at line 61 of file G4ParticleHPPhotonDist.hh.

      : repFlag( 0 ) 
      , targetMass( 0.0 ) 
      , nDiscrete( 0 ) 
      , isoFlag( 0 )
      , tabulationType( 0 )
      , nDiscrete2( 0 ) 
      , nIso( 0 ) 
      , nPartials( 0 )
      , theInternalConversionFlag( 0 )
      , nGammaEnergies( 0 )
      , theBaseEnergy( 0.0 )
  {
 
     disType = 0;
     energy = 0;
     theYield = 0;
     thePartialXsec = 0;
     theReactionXsec = 0;
     isPrimary = 0;
     theShells = 0;
     theGammas = 0;
     nNeu = 0;
     theLegendre = 0;
     theAngular = 0;
     distribution = 0;
     probs = 0;
     partials = 0;
     actualMult.Put( 0 );
 
     theLevelEnergies = 0;
     theTransitionProbabilities = 0;
     thePhotonTransitionFraction = 0;
 
  }

References actualMult, distribution, disType, energy, isPrimary, nNeu, partials, probs, G4Cache< VALTYPE >::Put(), theAngular, theGammas, theLegendre, theLevelEnergies, thePartialXsec, thePhotonTransitionFraction, theReactionXsec, theShells, theTransitionProbabilities, and theYield.

◆ ~G4ParticleHPPhotonDist()

G4ParticleHPPhotonDist::~G4ParticleHPPhotonDist ( )

inline

Definition at line 97 of file G4ParticleHPPhotonDist.hh.

  {
     delete [] disType;
     delete [] energy;
     delete [] theYield;
     delete [] thePartialXsec;
//     delete [] theReactionXsec;   
//     DHW: not created in this class 
     delete [] isPrimary;
     delete [] theShells;
     delete [] theGammas;
     delete [] nNeu;
     delete [] theAngular;
     delete [] distribution;
     delete [] probs;
 
     if ( theLegendre != 0 )
     {
        for ( G4int i = 0 ; i < (nDiscrete2-nIso) ; i++ )
           if ( theLegendre[i] != 0 ) delete[] theLegendre[i]; 
 
        delete [] theLegendre;
     }
 
     if ( partials != 0 ) 
     {
        for ( G4int i = 0 ; i < nPartials ; i++ )
           { delete partials[i]; }
 
        delete [] partials;
     }
 
     delete [] theLevelEnergies;
     delete [] theTransitionProbabilities;
     delete [] thePhotonTransitionFraction;
     if (actualMult.Get() != 0) delete actualMult.Get();
  }

References actualMult, distribution, disType, energy, G4Cache< VALTYPE >::Get(), isPrimary, nDiscrete2, nIso, nNeu, nPartials, partials, probs, theAngular, theGammas, theLegendre, theLevelEnergies, thePartialXsec, thePhotonTransitionFraction, theShells, theTransitionProbabilities, and theYield.

Member Function Documentation

◆ GetLevelEnergy()

G4double G4ParticleHPPhotonDist::GetLevelEnergy ( )

inline

Definition at line 149 of file G4ParticleHPPhotonDist.hh.

149{return theBaseEnergy;}

References theBaseEnergy.

Referenced by G4ParticleHPInelasticCompFS::CompositeApply().

◆ GetPhotons()

G4ReactionProductVector * G4ParticleHPPhotonDist::GetPhotons ( G4double anEnergy )

Definition at line 285 of file G4ParticleHPPhotonDist.cc.

{
 
  //G4cout << "G4ParticleHPPhotonDist::GetPhotons repFlag " << repFlag << G4endl;
  // the partial cross-section case is not in this yet. @@@@  << 070601 TK add partial 
  if ( actualMult.Get() == NULL ) {
     actualMult.Get() = new std::vector<G4int>( nDiscrete );
  }
  G4int i, ii, iii;
  G4int nSecondaries = 0;
  G4ReactionProductVector* thePhotons = new G4ReactionProductVector;
 
  if (repFlag==1) {
    G4double current=0;
    for (i = 0; i < nDiscrete; i++) {
      current = theYield[i].GetY(anEnergy);
      actualMult.Get()->at(i) = G4Poisson(current); // max cut-off still missing @@@
      if (nDiscrete == 1 && current < 1.0001) {
        actualMult.Get()->at(i) = static_cast<G4int>(current);
        if(current<1) 
        {
          actualMult.Get()->at(i) = 0;
          if(G4UniformRand()<current) actualMult.Get()->at(i) = 1;
        }
      }
      nSecondaries += actualMult.Get()->at(i);
    }
    //G4cout << "nSecondaries " << nSecondaries  << " anEnergy " << anEnergy/eV << G4endl;
    for (i = 0; i < nSecondaries; i++) {
      G4ReactionProduct * theOne = new G4ReactionProduct;
      theOne->SetDefinition(G4Gamma::Gamma());
      thePhotons->push_back(theOne);
    }
 
    G4int count = 0;
 
    if (nDiscrete == 1 && nPartials == 1) {
      if (actualMult.Get()->at(0) > 0) {
        if (disType[0] == 1) {
          // continuum
          G4ParticleHPVector* temp;
          temp = partials[ 0 ]->GetY(anEnergy); //@@@ look at, seems fishy
          G4double maximumE = temp->GetX( temp->GetVectorLength()-1 ); // This is an assumption.
 
          //G4cout << "start " << actualMult[ 0 ] << " maximumE " << maximumE/eV << G4endl;
 
          std::vector< G4double > photons_e_best( actualMult.Get()->at(0) , 0.0 );
          G4double best = DBL_MAX;
          G4int maxTry = 1000; 
          for (G4int j = 0; j < maxTry; j++) {
            std::vector<G4double> photons_e(actualMult.Get()->at(0), 0.0);
            for (std::vector<G4double>::iterator it = photons_e.begin(); it < photons_e.end(); it++) {
              *it = temp->Sample();
            }
 
            if (std::accumulate(photons_e.begin(), photons_e.end(), 0.0) > maximumE) {
              if (std::accumulate(photons_e.begin(), photons_e.end(), 0.0) < best)
                photons_e_best = photons_e;
              continue;
 
            } else {
              G4int iphot = 0;
              for (std::vector<G4double>::iterator it = photons_e.begin(); it < photons_e.end(); it++) {
                thePhotons->operator[](iphot)->SetKineticEnergy(*it);   // Replace index count, which was not incremented, 
                                                                        // with iphot, which is, as per Artem Zontikov,
                                                                        // bug report 2167
                iphot++;
              }  
              // G4cout << "OK " << actualMult[0] << " j " << j << " total photons E  " 
              //        << std::accumulate(photons_e.begin(), photons_e.end(), 0.0)/eV << " ratio " 
              //        << std::accumulate(photons_e.begin(), photons_e.end(), 0.0)/maximumE 
              //        << G4endl;
                    
              break;
            }
            //G4cout << "Not Good " << actualMult[0] << " j " << j << " total photons E  " 
            //       << best/eV << " ratio " << best / maximumE 
            //       << G4endl;
          }
               // TKDB
          delete temp;
 
        } else {
          // discrete
          thePhotons->operator[](count)->SetKineticEnergy(energy[i]);
        }
        count++;
    if (count > nSecondaries) throw G4HadronicException(__FILE__, __LINE__, "G4ParticleHPPhotonDist::GetPhotons inconsistency");
      }
         
    } else {  // nDiscrete != 1 or nPartials != 1
      for (i=0; i<nDiscrete; i++) { 
        for (ii=0; ii< actualMult.Get()->at(i); ii++) {   
          if (disType[i] == 1) {
            // continuum
            G4double  sum=0, run=0;
            for (iii = 0; iii < nPartials; iii++) sum+=probs[iii].GetY(anEnergy);
            G4double random = G4UniformRand();
            G4int theP = 0;
            for (iii = 0; iii < nPartials; iii++) {
              run+=probs[iii].GetY(anEnergy);
              theP = iii;
              if(random<run/sum) break;
            }
 
            if (theP == nPartials) theP=nPartials-1; // das sortiert J aus.
            sum = 0; 
            G4ParticleHPVector * temp;
            temp = partials[theP]->GetY(anEnergy); //@@@ look at, seems fishy
            G4double eGamm = temp->Sample();
            thePhotons->operator[](count)->SetKineticEnergy(eGamm);
            delete temp;
 
          } else { 
            // discrete
            thePhotons->operator[](count)->SetKineticEnergy(energy[i]);
      }
          count++;
          if (count > nSecondaries) throw G4HadronicException(__FILE__, __LINE__, "G4ParticleHPPhotonDist::GetPhotons inconsistency");
        }
      }
    }
 
    // now do the angular distributions...
    if (isoFlag == 1) {
      for (i=0; i< nSecondaries; i++)
      {
    G4double costheta = 2.*G4UniformRand()-1;
    G4double theta = std::acos(costheta);
    G4double phi = twopi*G4UniformRand();
    G4double sinth = std::sin(theta);
    G4double en = thePhotons->operator[](i)->GetTotalEnergy();
    G4ThreeVector temp(en*sinth*std::cos(phi), en*sinth*std::sin(phi), en*std::cos(theta) );
    thePhotons->operator[](i)->SetMomentum( temp ) ;
  //      G4cout << "Isotropic distribution in PhotonDist"<<temp<<G4endl;
      }
    }
    else
    {
      for(i=0; i<nSecondaries; i++)
      { 
    G4double currentEnergy = thePhotons->operator[](i)->GetTotalEnergy();
    for(ii=0; ii<nDiscrete2; ii++) 
    {
          if (std::abs(currentEnergy-theGammas[ii])<0.1*keV) break;
    }
    if(ii==nDiscrete2) ii--; // fix for what seems an (file12 vs file 14) inconsistency found in the ENDF 7N14 data. @@
    if(ii<nIso)
    {
          // isotropic distribution
          //
          //Fix Bugzilla report #1745
          //G4double theta = pi*G4UniformRand();
      G4double costheta = 2.*G4UniformRand()-1;
      G4double theta = std::acos(costheta);
          G4double phi = twopi*G4UniformRand();
          G4double sinth = std::sin(theta);
          G4double en = thePhotons->operator[](i)->GetTotalEnergy();
          // DHW G4ThreeVector tempVector(en*sinth*std::cos(phi), en*sinth*std::sin(phi), en*std::cos(theta) );
          G4ThreeVector tempVector(en*sinth*std::cos(phi), en*sinth*std::sin(phi), en*costheta );
          thePhotons->operator[](i)->SetMomentum( tempVector ) ;
    }
    else if(tabulationType==1)
    {
          // legendre polynomials
          G4int it(0);
          for (iii=0; iii<nNeu[ii-nIso]; iii++) // find the neutron energy
          {
            it = iii;
        if(theLegendre[ii-nIso][iii].GetEnergy()>anEnergy)
              break;
          }
          G4ParticleHPLegendreStore aStore(2);
          aStore.SetCoeff(1, &(theLegendre[ii-nIso][it]));  
          //aStore.SetCoeff(0, &(theLegendre[ii-nIso][it-1])); 
          //TKDB 110512
          if ( it > 0 ) 
          {
             aStore.SetCoeff(0, &(theLegendre[ii-nIso][it-1])); 
          }
          else
          {
             aStore.SetCoeff(0, &(theLegendre[ii-nIso][it])); 
          }
          G4double cosTh = aStore.SampleMax(anEnergy);
          G4double theta = std::acos(cosTh);
          G4double phi = twopi*G4UniformRand();
          G4double sinth = std::sin(theta);
          G4double en = thePhotons->operator[](i)->GetTotalEnergy();
          G4ThreeVector tempVector(en*sinth*std::cos(phi), en*sinth*std::sin(phi), en*std::cos(theta) );
          thePhotons->operator[](i)->SetMomentum( tempVector ) ;
    }
    else
    {
          // tabulation of probabilities.
          G4int it(0);
          for (iii=0; iii<nNeu[ii-nIso]; iii++) // find the neutron energy
          {
            it = iii;
        if(theAngular[ii-nIso][iii].GetEnergy()>anEnergy)
              break;
          }
          G4double costh = theAngular[ii-nIso][it].GetCosTh(); // no interpolation yet @@
          G4double theta = std::acos(costh);
          G4double phi = twopi*G4UniformRand();
          G4double sinth = std::sin(theta);
          G4double en = thePhotons->operator[](i)->GetTotalEnergy();
          G4ThreeVector tmpVector(en*sinth*std::cos(phi), en*sinth*std::sin(phi), en*costh );
          thePhotons->operator[](i)->SetMomentum( tmpVector ) ;
    }
      }  
    }
 
  } else if (repFlag == 2) {
    G4double * running = new G4double[nGammaEnergies];
    running[0]=theTransitionProbabilities[0];
    //G4int i; //declaration at 284th
    for(i=1; i<nGammaEnergies; i++)
    {
      running[i]=running[i-1]+theTransitionProbabilities[i];
    }
    G4double random = G4UniformRand();
    G4int it=0;
    for(i=0; i<nGammaEnergies; i++)
    {
      it = i;
      if(random < running[i]/running[nGammaEnergies-1]) break;
    }
    delete [] running;
    G4double totalEnergy = theBaseEnergy - theLevelEnergies[it];
    G4ReactionProduct * theOne = new G4ReactionProduct;
    theOne->SetDefinition(G4Gamma::Gamma());
    random = G4UniformRand();
    if(theInternalConversionFlag==2 && random>thePhotonTransitionFraction[it])
    {
      theOne->SetDefinition(G4Electron::Electron());
      //Bug reported Chao Zhang (Chao.Zhang@usd.edu), Dongming Mei(Dongming.Mei@usd.edu) Feb. 25, 2009 
      //But never enter at least with G4NDL3.13
      totalEnergy += G4Electron::Electron()->GetPDGMass(); //proposed correction: add this line for electron
    }
    theOne->SetTotalEnergy(totalEnergy);
    if( isoFlag == 1 )
    {
      G4double costheta = 2.*G4UniformRand()-1;
      G4double theta = std::acos(costheta);
      G4double phi = twopi*G4UniformRand();
      G4double sinth = std::sin(theta);
      //Bug reported Chao Zhang (Chao.Zhang@usd.edu), Dongming Mei(Dongming.Mei@usd.edu) Feb. 25, 2009 
      //G4double en = theOne->GetTotalEnergy();
      G4double en = theOne->GetTotalMomentum();
      //But never cause real effect at least with G4NDL3.13 TK
      G4ThreeVector temp(en*sinth*std::cos(phi), en*sinth*std::sin(phi), en*std::cos(theta) );
      theOne->SetMomentum( temp ) ;
    }
    else
    {
      G4double currentEnergy = theOne->GetTotalEnergy();
      for(ii=0; ii<nDiscrete2; ii++) 
      {
        if (std::abs(currentEnergy-theGammas[ii])<0.1*keV) break;
      }
      if(ii==nDiscrete2) ii--; // fix for what seems an (file12 vs file 14) inconsistency found in the ENDF 7N14 data. @@
      if(ii<nIso)
      {
        //Bug reported Chao Zhang (Chao.Zhang@usd.edu), Dongming Mei(Dongming.Mei@usd.edu) Feb. 25, 2009 
        // isotropic distribution
        //G4double theta = pi*G4UniformRand();
        G4double theta = std::acos(2.*G4UniformRand()-1.);
        //But this is alos never cause real effect at least with G4NDL3.13 TK  not repFlag == 2 AND isoFlag != 1
        G4double phi = twopi*G4UniformRand();
        G4double sinth = std::sin(theta);
        //Bug reported Chao Zhang (Chao.Zhang@usd.edu), Dongming Mei(Dongming.Mei@usd.edu) Feb. 25, 2009 
        //G4double en = theOne->GetTotalEnergy();
        G4double en = theOne->GetTotalMomentum();
        //But never cause real effect at least with G4NDL3.13 TK
        G4ThreeVector tempVector(en*sinth*std::cos(phi), en*sinth*std::sin(phi), en*std::cos(theta) );
        theOne->SetMomentum( tempVector ) ;
      }
      else if(tabulationType==1)
      {
        // legendre polynomials
        G4int itt(0);
        for (iii=0; iii<nNeu[ii-nIso]; iii++) // find the neutron energy
        {
          itt = iii;
      if(theLegendre[ii-nIso][iii].GetEnergy()>anEnergy)
            break;
        }
        G4ParticleHPLegendreStore aStore(2);
        aStore.SetCoeff(1, &(theLegendre[ii-nIso][itt]));  
        //aStore.SetCoeff(0, &(theLegendre[ii-nIso][it-1])); 
        //TKDB 110512
        if ( itt > 0 ) 
        {
           aStore.SetCoeff(0, &(theLegendre[ii-nIso][itt-1])); 
        }
        else
        {
           aStore.SetCoeff(0, &(theLegendre[ii-nIso][itt])); 
        }
        G4double cosTh = aStore.SampleMax(anEnergy);
        G4double theta = std::acos(cosTh);
        G4double phi = twopi*G4UniformRand();
        G4double sinth = std::sin(theta);
        //Bug reported Chao Zhang (Chao.Zhang@usd.edu), Dongming Mei(Dongming.Mei@usd.edu) Feb. 25, 2009 
        //G4double en = theOne->GetTotalEnergy();
        G4double en = theOne->GetTotalMomentum();
        //But never cause real effect at least with G4NDL3.13 TK
        G4ThreeVector tempVector(en*sinth*std::cos(phi), en*sinth*std::sin(phi), en*std::cos(theta) );
        theOne->SetMomentum( tempVector ) ;
      }
      else
      {
        // tabulation of probabilities.
        G4int itt(0);
        for (iii=0; iii<nNeu[ii-nIso]; iii++) // find the neutron energy
        {
          itt = iii;
      if(theAngular[ii-nIso][iii].GetEnergy()>anEnergy)
            break;
        }
        G4double costh = theAngular[ii-nIso][itt].GetCosTh(); // no interpolation yet @@
        G4double theta = std::acos(costh);
        G4double phi = twopi*G4UniformRand();
        G4double sinth = std::sin(theta);
        //Bug reported Chao Zhang (Chao.Zhang@usd.edu), Dongming Mei(Dongming.Mei@usd.edu) Feb. 25, 2009 
        //G4double en = theOne->GetTotalEnergy();
        G4double en = theOne->GetTotalMomentum();
        //But never cause real effect at least with G4NDL3.13 TK
        G4ThreeVector tmpVector(en*sinth*std::cos(phi), en*sinth*std::sin(phi), en*costh );
        theOne->SetMomentum( tmpVector ) ;
      }
    }
    thePhotons->push_back(theOne);
  }
  else if( repFlag==0 )
  {
 
// TK add 
      if ( thePartialXsec == 0 ) 
      {
         //G4cout << "repFlag is 0, but no PartialXsec data" << G4endl;
         //G4cout << "This is not support yet." << G4endl;
         return thePhotons;
      }
 
// Partial Case 
 
      G4ReactionProduct * theOne = new G4ReactionProduct;
      theOne->SetDefinition( G4Gamma::Gamma() );
      thePhotons->push_back( theOne );
 
// Energy 
 
      //G4cout << "Partial Case nDiscrete " << nDiscrete << G4endl;
      G4double sum = 0.0; 
      std::vector < G4double > dif( nDiscrete , 0.0 ); 
      for ( G4int j = 0 ; j < nDiscrete ; j++ ) 
      {
         G4double x = thePartialXsec[ j ].GetXsec( anEnergy );  // x in barn 
         if ( x > 0 ) 
         {
            sum += x;   
         } 
         dif [ j ] = sum; 
         //G4cout << "j " << j << ", x " << x << ", dif " << dif [ j ] << G4endl;
      }
      
      G4double rand = G4UniformRand();
 
      G4int iphoton = 0; 
      for ( G4int j = 0 ; j < nDiscrete ; j++ ) 
      {
         G4double y = rand*sum; 
         if ( dif [ j ] > y ) 
         {
            iphoton = j; 
            break;  
         } 
      }
      //G4cout << "iphoton " << iphoton << G4endl;
      //G4cout << "photon energy " << theGammas[ iphoton ] /eV  << G4endl;
 
      // Statistically suppress the photon according to reaction cross section  
      // Fix proposed by Artem Zontikov, Bug report #1824
      if (theReactionXsec) {
        if (thePartialXsec[iphoton].GetXsec(anEnergy)/theReactionXsec->GetXsec(anEnergy) < G4UniformRand() ) {
          delete thePhotons;
          thePhotons = 0;
          return thePhotons;
        }
      }
 
// Angle 
      G4double cosTheta = 0.0; // mu
 
      if ( isoFlag == 1 )
      {
 
//       Isotropic Case
 
         cosTheta = 2.*G4UniformRand()-1;
 
      }
      else
      {
 
         if ( iphoton < nIso )
         {
 
//          still Isotropic 
 
            cosTheta = 2.*G4UniformRand()-1;
 
         }
         else
         {
 
            //G4cout << "Not Isotropic and isoFlag " << isoFlag  << G4endl;
            //G4cout << "tabulationType " << tabulationType << G4endl;
            //G4cout << "nDiscrete2  " << nDiscrete2 << G4endl;
            //G4cout << "nIso " << nIso << G4endl;
            //G4cout << "size of nNeu " << nDiscrete2-nIso << G4endl;
            //G4cout << "nNeu[iphoton-nIso] " << nNeu[iphoton-nIso] << G4endl;
 
            if ( tabulationType == 1 )
            {
//             legendre polynomials
 
               G4int iangle = 0; 
               for ( G4int j = 0 ; j < nNeu [ iphoton - nIso ] ; j++ )
               {
                  iangle = j;
                  if ( theLegendre[ iphoton - nIso ][ j ].GetEnergy() > anEnergy ) break;
               }
 
               G4ParticleHPLegendreStore aStore( 2 );
               aStore.SetCoeff( 1 , &( theLegendre[ iphoton - nIso ][ iangle ] ) );  
               aStore.SetCoeff( 0 , &( theLegendre[ iphoton - nIso ][ iangle - 1 ] ) ); 
 
               cosTheta = aStore.SampleMax( anEnergy );
 
            }
            else if ( tabulationType == 2 )
            {
        
//             tabulation of probabilities.
 
               G4int iangle = 0; 
               for ( G4int j = 0 ; j < nNeu [ iphoton - nIso ] ; j++ )
               {
                  iangle = j;
                  if ( theAngular[ iphoton - nIso ][ j ].GetEnergy() > anEnergy ) break;
               }
 
               cosTheta = theAngular[iphoton-nIso][ iangle ].GetCosTh(); // no interpolation yet @@
 
            }
         }
      }
      
// Set 
      G4double phi = twopi*G4UniformRand();
      G4double theta = std::acos( cosTheta );
      G4double sinTheta = std::sin( theta );
 
      G4double photonE = theGammas[ iphoton ];
      G4ThreeVector direction ( sinTheta*std::cos( phi ) , sinTheta * std::sin( phi ) , cosTheta );
      G4ThreeVector photonP = photonE * direction;
      thePhotons->operator[]( 0 )->SetMomentum( photonP ) ;
    
  }
  else
  {
    delete thePhotons;
    thePhotons = 0; // no gamma data available; some work needed @@@@@@@
  }    
  return thePhotons;
}

References actualMult, DBL_MAX, disType, G4Electron::Electron(), energy, G4Poisson(), G4UniformRand, G4Gamma::Gamma(), G4Cache< VALTYPE >::Get(), G4ParticleHPAngularP::GetCosTh(), G4ParticleDefinition::GetPDGMass(), G4ReactionProduct::GetTotalEnergy(), G4ReactionProduct::GetTotalMomentum(), G4ParticleHPVector::GetVectorLength(), G4ParticleHPVector::GetX(), G4ParticleHPVector::GetXsec(), G4ParticleHPVector::GetY(), G4ParticleHPPartial::GetY(), isoFlag, keV, nDiscrete, nDiscrete2, nGammaEnergies, nIso, nNeu, nPartials, partials, probs, repFlag, G4ParticleHPVector::Sample(), G4ParticleHPLegendreStore::SampleMax(), G4ParticleHPLegendreStore::SetCoeff(), G4ReactionProduct::SetDefinition(), G4ReactionProduct::SetMomentum(), G4ReactionProduct::SetTotalEnergy(), tabulationType, theAngular, theBaseEnergy, theGammas, theInternalConversionFlag, theLegendre, theLevelEnergies, thePartialXsec, thePhotonTransitionFraction, theReactionXsec, theTransitionProbabilities, theYield, and twopi.

Referenced by G4ParticleHPCaptureFS::ApplyYourself(), G4ParticleHPInelasticBaseFS::BaseApply(), G4ParticleHPInelasticCompFS::CompositeApply(), and G4ParticleHPFSFissionFS::GetPhotons().

◆ GetTargetMass()

G4double G4ParticleHPPhotonDist::GetTargetMass ( )

inline

Definition at line 145 of file G4ParticleHPPhotonDist.hh.

145{return targetMass;}

References targetMass.

Referenced by G4ParticleHPCaptureFS::Init().

◆ InitAngular()

void G4ParticleHPPhotonDist::InitAngular ( std::istream & aDataFile )

Definition at line 122 of file G4ParticleHPPhotonDist.cc.

{
  G4int i, ii;
  //angular distributions
  aDataFile >> isoFlag;
  if (isoFlag != 1)
  {
    if (repFlag == 2) G4cout << "G4ParticleHPPhotonDist: repFlag == 2 && isoFlag != 1 is unexpected! If you use G4ND3.x, then please report to Geant4 HyperNews. " << G4endl;
    aDataFile >> tabulationType >> nDiscrete2 >> nIso;
//080731
      if (theGammas != NULL && nDiscrete2 != nDiscrete) 
        G4cout << "080731c G4ParticleHPPhotonDist nDiscrete2 != nDiscrete, It looks like something wrong in your NDL files. Please update the latest. If you still have this messages after the update, then please report to Geant4 Hyper News." << G4endl;
 
      // The order of cross section (InitPartials) and distribution
      // (InitAngular here) data are different, we have to re-coordinate
      // consistent data order.
      std::vector < G4double > vct_gammas_par; 
      std::vector < G4double > vct_shells_par; 
      std::vector < G4int > vct_primary_par; 
      std::vector < G4int > vct_distype_par; 
      std::vector < G4ParticleHPVector* > vct_pXS_par;
      if ( theGammas != NULL && theShells != NULL ) 
      {
         //copy the cross section data 
         for ( i = 0 ; i < nDiscrete ; i++ )
         {
            vct_gammas_par.push_back( theGammas[ i ] );
            vct_shells_par.push_back( theShells[ i ] );
            vct_primary_par.push_back( isPrimary[ i ] );
            vct_distype_par.push_back( disType[ i ] );
            G4ParticleHPVector* hpv = new G4ParticleHPVector;
            *hpv = thePartialXsec[ i ];
            vct_pXS_par.push_back( hpv );
         }
      }
     if ( theGammas == NULL ) theGammas = new G4double[nDiscrete2];
     if ( theShells == NULL ) theShells = new G4double[nDiscrete2];
//080731
 
    for (i=0; i< nIso; i++) // isotropic photons
    {
       aDataFile >> theGammas[i] >> theShells[i];
       theGammas[i]*=eV;
       theShells[i]*=eV;
    }
    nNeu = new G4int [nDiscrete2-nIso];
    if(tabulationType==1)theLegendre=new G4ParticleHPLegendreTable *[nDiscrete2-nIso];
    if(tabulationType==2)theAngular =new G4ParticleHPAngularP *[nDiscrete2-nIso];
    for(i=nIso; i< nDiscrete2; i++)
    {
      if(tabulationType==1) 
      {
        aDataFile >> theGammas[i] >> theShells[i] >> nNeu[i-nIso];
        theGammas[i]*=eV;
        theShells[i]*=eV;
        theLegendre[i-nIso]=new G4ParticleHPLegendreTable[nNeu[i-nIso]];
        theLegendreManager.Init(aDataFile); 
        for (ii=0; ii<nNeu[i-nIso]; ii++)
        {
          theLegendre[i-nIso][ii].Init(aDataFile);
        }
      }
      else if(tabulationType==2)
      {
        aDataFile >> theGammas[i] >> theShells[i] >> nNeu[i-nIso];
        theGammas[i]*=eV;
        theShells[i]*=eV;
        theAngular[i-nIso]=new G4ParticleHPAngularP[nNeu[i-nIso]];
        for (ii=0; ii<nNeu[i-nIso]; ii++)
        {
          theAngular[i-nIso][ii].Init(aDataFile);
        }
      }
      else
      {
        G4cout << "tabulation type: tabulationType"<<G4endl;
        throw G4HadronicException(__FILE__, __LINE__, "cannot deal with this tabulation type for angular distributions.");
      }
    }
//080731
     if ( vct_gammas_par.size() > 0 ) 
     {
        //Reordering cross section data to corrsponding distribution data 
        for ( i = 0 ; i < nDiscrete ; i++ ) 
        {
           for ( G4int j = 0 ; j < nDiscrete ; j++ )  
           {
              // Checking gamma and shell to identification 
              if ( theGammas[ i ] == vct_gammas_par [ j ] && theShells [ i ] == vct_shells_par[ j ] )
              {
                 isPrimary [ i ] = vct_primary_par [ j ];
                 disType [ i ] = vct_distype_par [ j ];
                 thePartialXsec[ i ] = ( *( vct_pXS_par[ j ] ) );
              }
           }
        }
        //Garbage collection 
        for ( std::vector < G4ParticleHPVector* >::iterator 
           it = vct_pXS_par.begin() ; it != vct_pXS_par.end() ; it++ )
        {
           delete *it;
        }
     }
//080731
  }
}

References disType, eV, G4cout, G4endl, G4InterpolationManager::Init(), G4ParticleHPAngularP::Init(), G4ParticleHPLegendreTable::Init(), isoFlag, isPrimary, nDiscrete, nDiscrete2, nIso, nNeu, repFlag, tabulationType, theAngular, theGammas, theLegendre, theLegendreManager, thePartialXsec, and theShells.

Referenced by G4FissionLibrary::Init(), G4ParticleHPCaptureFS::Init(), G4ParticleHPFSFissionFS::Init(), G4ParticleHPInelasticCompFS::Init(), and G4ParticleHPInelasticBaseFS::Init().

◆ InitEnergies()

void G4ParticleHPPhotonDist::InitEnergies ( std::istream & aDataFile )

Definition at line 230 of file G4ParticleHPPhotonDist.cc.

{
  G4int i, energyDistributionsNeeded = 0;
  for (i=0; i<nDiscrete; i++)
  {
    if( disType[i]==1) energyDistributionsNeeded =1;
  }
  if(!energyDistributionsNeeded) return;
  aDataFile >>  nPartials;
  distribution = new G4int[nPartials];
  probs = new G4ParticleHPVector[nPartials];
  partials = new G4ParticleHPPartial * [nPartials];
  G4int nen;
  G4int dummy;
  for (i=0; i<nPartials; i++)
  {
    aDataFile >> dummy;
    probs[i].Init(aDataFile, eV);
    aDataFile >> nen;
    partials[i] = new G4ParticleHPPartial(nen);
    partials[i]->InitInterpolation(aDataFile);
    partials[i]->Init(aDataFile);
  }
}

References distribution, disType, eV, G4ParticleHPPartial::Init(), G4ParticleHPVector::Init(), G4ParticleHPPartial::InitInterpolation(), nDiscrete, nPartials, partials, and probs.

Referenced by G4FissionLibrary::Init(), G4ParticleHPCaptureFS::Init(), G4ParticleHPFSFissionFS::Init(), G4ParticleHPInelasticCompFS::Init(), and G4ParticleHPInelasticBaseFS::Init().

◆ InitMean()

G4bool G4ParticleHPPhotonDist::InitMean ( std::istream & aDataFile )

Definition at line 55 of file G4ParticleHPPhotonDist.cc.

{
 
  G4bool result = true;
  if(aDataFile >> repFlag)
  {
 
    aDataFile >> targetMass;
    if(repFlag==1)
    {
    // multiplicities
      aDataFile >> nDiscrete;
      disType = new G4int[nDiscrete];
      energy = new G4double[nDiscrete];
      //actualMult = new G4int[nDiscrete];
      theYield = new G4ParticleHPVector[nDiscrete];
      for (G4int i=0; i<nDiscrete; i++)
      {
    aDataFile >> disType[i]>>energy[i];
    energy[i]*=eV;
    theYield[i].Init(aDataFile, eV);
      }
    }
    else if(repFlag == 2)
    {
       aDataFile >> theInternalConversionFlag;
       aDataFile >> theBaseEnergy;
       theBaseEnergy*=eV;
       aDataFile >> theInternalConversionFlag;
       // theInternalConversionFlag == 1 No IC, theInternalConversionFlag == 2 with IC 
       aDataFile >> nGammaEnergies;
       theLevelEnergies = new G4double[nGammaEnergies];
       theTransitionProbabilities = new G4double[nGammaEnergies];
       if(theInternalConversionFlag == 2) thePhotonTransitionFraction = new G4double[nGammaEnergies];
       for(G4int  ii=0; ii<nGammaEnergies; ii++)
       {
     if(theInternalConversionFlag == 1)
     {
           aDataFile >> theLevelEnergies[ii] >> theTransitionProbabilities[ii];
       theLevelEnergies[ii]*=eV;
     }
     else if(theInternalConversionFlag == 2)
     {
           aDataFile >> theLevelEnergies[ii] >> theTransitionProbabilities[ii] >> thePhotonTransitionFraction[ii];
       theLevelEnergies[ii]*=eV;
     }
     else
     {
           throw G4HadronicException(__FILE__, __LINE__, "G4ParticleHPPhotonDist: Unknown conversion flag");
     }
      }
       // Note, that this is equivalent to using the 'Gamma' classes.
      // throw G4HadronicException(__FILE__, __LINE__, "G4ParticleHPPhotonDist: Transition probability array not sampled for the moment.");
    }
    else
    {
      G4cout << "Data representation in G4ParticleHPPhotonDist: "<<repFlag<<G4endl;
      throw G4HadronicException(__FILE__, __LINE__, "G4ParticleHPPhotonDist: This data representation is not implemented.");
    }
  }
  else
  {
    result = false;
  }
  return result;
}

References disType, energy, eV, G4cout, G4endl, G4ParticleHPVector::Init(), nDiscrete, nGammaEnergies, repFlag, targetMass, theBaseEnergy, theInternalConversionFlag, theLevelEnergies, thePhotonTransitionFraction, theTransitionProbabilities, and theYield.

Referenced by G4FissionLibrary::Init(), G4ParticleHPCaptureFS::Init(), G4ParticleHPFSFissionFS::Init(), G4ParticleHPInelasticCompFS::Init(), and G4ParticleHPInelasticBaseFS::Init().

◆ InitPartials()

void G4ParticleHPPhotonDist::InitPartials	(	std::istream &	aDataFile,
		G4ParticleHPVector *	theXsec = `0`
	)

Definition at line 255 of file G4ParticleHPPhotonDist.cc.

{
  if (theXsec) theReactionXsec = theXsec;
 
  aDataFile >> nDiscrete >> targetMass;
  if(nDiscrete != 1)
  {
    theTotalXsec.Init(aDataFile, eV);
  }
  G4int i;
  theGammas = new G4double[nDiscrete];
  theShells = new G4double[nDiscrete];
  isPrimary = new G4int[nDiscrete];
  disType = new G4int[nDiscrete];
  thePartialXsec = new G4ParticleHPVector[nDiscrete];
  for(i=0; i<nDiscrete; i++)
  {
    aDataFile>>theGammas[i]>>theShells[i]>>isPrimary[i]>>disType[i];
    theGammas[i]*=eV;
    theShells[i]*=eV;
    thePartialXsec[i].Init(aDataFile, eV);
  }  
 
  //G4cout << "G4ParticleHPPhotonDist::InitPartials Test " << G4endl;
  //G4cout << "G4ParticleHPPhotonDist::InitPartials nDiscrete " << nDiscrete << G4endl;
  //G4ParticleHPVector* aHP = new G4ParticleHPVector;
  //aHP->Check(1);
}

References disType, eV, G4ParticleHPVector::Init(), isPrimary, nDiscrete, targetMass, theGammas, thePartialXsec, theReactionXsec, theShells, and theTotalXsec.

Referenced by G4ParticleHPInelasticCompFS::Init(), and G4ParticleHPInelasticBaseFS::Init().

◆ NeedsCascade()

G4bool G4ParticleHPPhotonDist::NeedsCascade ( )

inline

Definition at line 147 of file G4ParticleHPPhotonDist.hh.

147{return repFlag==2;}

References repFlag.

Field Documentation

◆ actualMult

G4Cache< std::vector<G4int>* > G4ParticleHPPhotonDist::actualMult

private

Definition at line 181 of file G4ParticleHPPhotonDist.hh.

Referenced by G4ParticleHPPhotonDist(), GetPhotons(), and ~G4ParticleHPPhotonDist().

◆ distribution

G4int* G4ParticleHPPhotonDist::distribution

private

Definition at line 176 of file G4ParticleHPPhotonDist.hh.

Referenced by G4ParticleHPPhotonDist(), InitEnergies(), and ~G4ParticleHPPhotonDist().

◆ disType

G4int* G4ParticleHPPhotonDist::disType

private

Definition at line 157 of file G4ParticleHPPhotonDist.hh.

Referenced by G4ParticleHPPhotonDist(), GetPhotons(), InitAngular(), InitEnergies(), InitMean(), InitPartials(), and ~G4ParticleHPPhotonDist().

◆ energy

G4double* G4ParticleHPPhotonDist::energy

private

Definition at line 158 of file G4ParticleHPPhotonDist.hh.

Referenced by G4ParticleHPPhotonDist(), GetPhotons(), InitMean(), and ~G4ParticleHPPhotonDist().

◆ isoFlag

G4int G4ParticleHPPhotonDist::isoFlag

private

Definition at line 165 of file G4ParticleHPPhotonDist.hh.

Referenced by GetPhotons(), and InitAngular().

◆ isPrimary

G4int* G4ParticleHPPhotonDist::isPrimary

private

Definition at line 163 of file G4ParticleHPPhotonDist.hh.

Referenced by G4ParticleHPPhotonDist(), InitAngular(), InitPartials(), and ~G4ParticleHPPhotonDist().

◆ nDiscrete

G4int G4ParticleHPPhotonDist::nDiscrete

private

Definition at line 156 of file G4ParticleHPPhotonDist.hh.

Referenced by GetPhotons(), InitAngular(), InitEnergies(), InitMean(), and InitPartials().

◆ nDiscrete2

G4int G4ParticleHPPhotonDist::nDiscrete2

private

Definition at line 167 of file G4ParticleHPPhotonDist.hh.

Referenced by GetPhotons(), InitAngular(), and ~G4ParticleHPPhotonDist().

◆ nGammaEnergies

G4int G4ParticleHPPhotonDist::nGammaEnergies

private

Definition at line 185 of file G4ParticleHPPhotonDist.hh.

Referenced by GetPhotons(), and InitMean().

◆ nIso

G4int G4ParticleHPPhotonDist::nIso

private

Definition at line 168 of file G4ParticleHPPhotonDist.hh.

Referenced by GetPhotons(), InitAngular(), and ~G4ParticleHPPhotonDist().

◆ nNeu

G4int* G4ParticleHPPhotonDist::nNeu

private

Definition at line 171 of file G4ParticleHPPhotonDist.hh.

Referenced by G4ParticleHPPhotonDist(), GetPhotons(), InitAngular(), and ~G4ParticleHPPhotonDist().

◆ nPartials

G4int G4ParticleHPPhotonDist::nPartials

private

Definition at line 177 of file G4ParticleHPPhotonDist.hh.

Referenced by GetPhotons(), InitEnergies(), and ~G4ParticleHPPhotonDist().

◆ partials

G4ParticleHPPartial** G4ParticleHPPhotonDist::partials

private

Definition at line 179 of file G4ParticleHPPhotonDist.hh.

Referenced by G4ParticleHPPhotonDist(), GetPhotons(), InitEnergies(), and ~G4ParticleHPPhotonDist().

◆ probs

G4ParticleHPVector* G4ParticleHPPhotonDist::probs

private

Definition at line 178 of file G4ParticleHPPhotonDist.hh.

Referenced by G4ParticleHPPhotonDist(), GetPhotons(), InitEnergies(), and ~G4ParticleHPPhotonDist().

◆ repFlag

G4int G4ParticleHPPhotonDist::repFlag

private

Definition at line 153 of file G4ParticleHPPhotonDist.hh.

Referenced by GetPhotons(), InitAngular(), InitMean(), and NeedsCascade().

◆ tabulationType

G4int G4ParticleHPPhotonDist::tabulationType

private

Definition at line 166 of file G4ParticleHPPhotonDist.hh.

Referenced by GetPhotons(), and InitAngular().

◆ targetMass

G4double G4ParticleHPPhotonDist::targetMass

private

Definition at line 154 of file G4ParticleHPPhotonDist.hh.

Referenced by GetTargetMass(), InitMean(), and InitPartials().

◆ theAngular

G4ParticleHPAngularP** G4ParticleHPPhotonDist::theAngular

private

Definition at line 174 of file G4ParticleHPPhotonDist.hh.

Referenced by G4ParticleHPPhotonDist(), GetPhotons(), InitAngular(), and ~G4ParticleHPPhotonDist().

◆ theBaseEnergy

G4double G4ParticleHPPhotonDist::theBaseEnergy

private

Definition at line 186 of file G4ParticleHPPhotonDist.hh.

Referenced by GetLevelEnergy(), GetPhotons(), and InitMean().

◆ theGammas

G4double* G4ParticleHPPhotonDist::theGammas

private

Definition at line 170 of file G4ParticleHPPhotonDist.hh.

Referenced by G4ParticleHPPhotonDist(), GetPhotons(), InitAngular(), InitPartials(), and ~G4ParticleHPPhotonDist().

◆ theInt

G4ParticleHPInterpolator G4ParticleHPPhotonDist::theInt

private

Definition at line 193 of file G4ParticleHPPhotonDist.hh.

◆ theInternalConversionFlag

G4int G4ParticleHPPhotonDist::theInternalConversionFlag

private

Definition at line 184 of file G4ParticleHPPhotonDist.hh.

Referenced by GetPhotons(), and InitMean().

◆ theLegend

G4ParticleHPFastLegendre G4ParticleHPPhotonDist::theLegend

private

Definition at line 192 of file G4ParticleHPPhotonDist.hh.

◆ theLegendre

G4ParticleHPLegendreTable** G4ParticleHPPhotonDist::theLegendre

private

Definition at line 173 of file G4ParticleHPPhotonDist.hh.

Referenced by G4ParticleHPPhotonDist(), GetPhotons(), InitAngular(), and ~G4ParticleHPPhotonDist().

◆ theLegendreManager

G4InterpolationManager G4ParticleHPPhotonDist::theLegendreManager

private

Definition at line 172 of file G4ParticleHPPhotonDist.hh.

Referenced by InitAngular().

◆ theLevelEnergies

G4double* G4ParticleHPPhotonDist::theLevelEnergies

private

Definition at line 187 of file G4ParticleHPPhotonDist.hh.

Referenced by G4ParticleHPPhotonDist(), GetPhotons(), InitMean(), and ~G4ParticleHPPhotonDist().

◆ thePartialXsec

G4ParticleHPVector* G4ParticleHPPhotonDist::thePartialXsec

private

Definition at line 161 of file G4ParticleHPPhotonDist.hh.

Referenced by G4ParticleHPPhotonDist(), GetPhotons(), InitAngular(), InitPartials(), and ~G4ParticleHPPhotonDist().

◆ thePhotonTransitionFraction

G4double* G4ParticleHPPhotonDist::thePhotonTransitionFraction

private

Definition at line 189 of file G4ParticleHPPhotonDist.hh.

Referenced by G4ParticleHPPhotonDist(), GetPhotons(), InitMean(), and ~G4ParticleHPPhotonDist().

◆ theReactionXsec

G4ParticleHPVector* G4ParticleHPPhotonDist::theReactionXsec

private

Definition at line 162 of file G4ParticleHPPhotonDist.hh.

Referenced by G4ParticleHPPhotonDist(), GetPhotons(), and InitPartials().

◆ theShells

G4double* G4ParticleHPPhotonDist::theShells

private

Definition at line 169 of file G4ParticleHPPhotonDist.hh.

Referenced by G4ParticleHPPhotonDist(), InitAngular(), InitPartials(), and ~G4ParticleHPPhotonDist().

◆ theTotalXsec

G4ParticleHPVector G4ParticleHPPhotonDist::theTotalXsec

private

Definition at line 160 of file G4ParticleHPPhotonDist.hh.

Referenced by InitPartials().

◆ theTransitionProbabilities

G4double* G4ParticleHPPhotonDist::theTransitionProbabilities

private

Definition at line 188 of file G4ParticleHPPhotonDist.hh.

Referenced by G4ParticleHPPhotonDist(), GetPhotons(), InitMean(), and ~G4ParticleHPPhotonDist().

◆ theYield

G4ParticleHPVector* G4ParticleHPPhotonDist::theYield

private

Definition at line 159 of file G4ParticleHPPhotonDist.hh.

Referenced by G4ParticleHPPhotonDist(), GetPhotons(), InitMean(), and ~G4ParticleHPPhotonDist().

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

source/processes/hadronic/models/particle_hp/include/G4ParticleHPPhotonDist.hh
source/processes/hadronic/models/particle_hp/src/G4ParticleHPPhotonDist.cc

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4ParticleHPPhotonDist()

◆ ~G4ParticleHPPhotonDist()

Member Function Documentation

◆ GetLevelEnergy()

◆ GetPhotons()

◆ GetTargetMass()

◆ InitAngular()

◆ InitEnergies()

◆ InitMean()

◆ InitPartials()

◆ NeedsCascade()

Field Documentation

◆ actualMult

◆ distribution

◆ disType

◆ energy

◆ isoFlag

◆ isPrimary

◆ nDiscrete

◆ nDiscrete2

◆ nGammaEnergies

◆ nIso

◆ nNeu

◆ nPartials

◆ partials

◆ probs

◆ repFlag

◆ tabulationType

◆ targetMass

◆ theAngular

◆ theBaseEnergy

◆ theGammas

◆ theInt

◆ theInternalConversionFlag

◆ theLegend

◆ theLegendre

◆ theLegendreManager

◆ theLevelEnergies

◆ thePartialXsec

◆ thePhotonTransitionFraction

◆ theReactionXsec

◆ theShells

◆ theTotalXsec

◆ theTransitionProbabilities

◆ theYield