Geant4-11
Public Member Functions | Protected Member Functions | Private Attributes
G4ParticleHPContEnergyAngular Class Reference

#include <G4ParticleHPContEnergyAngular.hh>

Inheritance diagram for G4ParticleHPContEnergyAngular:
G4VParticleHPEnergyAngular

Public Member Functions

void ClearHistories ()
 
 G4ParticleHPContEnergyAngular (G4ParticleDefinition *proj)
 
G4ReactionProduct * GetCMS ()
 
G4ReactionProduct * GetProjectileRP ()
 
G4ReactionProduct * GetTarget ()
 
void Init (std::istream &aDataFile)
 
G4double MeanEnergyOfThisInteraction ()
 
G4ReactionProduct * Sample (G4double anEnergy, G4double massCode, G4double mass)
 
void SetProjectileRP (G4ReactionProduct *aIncidentParticleRP)
 
void SetQValue (G4double aValue)
 
void SetTarget (G4ReactionProduct *aTarget)
 
 ~G4ParticleHPContEnergyAngular ()
 

Protected Member Functions

G4double GetQValue ()
 

Private Attributes

G4Cache< G4double > currentMeanEnergy
 
G4Cache< toBeCached > fCache
 
G4Cache< G4ParticleHPContAngularPar * > fCacheAngular
 
G4int nEnergy
 
G4ParticleHPContAngularPar * theAngular
 
G4int theAngularRep
 
G4int theInterpolation
 
G4InterpolationManager theManager
 
G4ParticleDefinition * theProjectile
 
G4double theQValue
 
G4double theTargetCode
 

Detailed Description

Definition at line 48 of file G4ParticleHPContEnergyAngular.hh.

Constructor & Destructor Documentation

◆ G4ParticleHPContEnergyAngular()

G4ParticleHPContEnergyAngular::G4ParticleHPContEnergyAngular ( G4ParticleDefinition *  proj)
inline

Definition at line 52 of file G4ParticleHPContEnergyAngular.hh.

53 : theProjectile(proj)
54 {
55 theAngular = 0;
56 currentMeanEnergy.Put( -2 );
57 theTargetCode = -1.0;
58 theAngularRep = -1;
59 nEnergy = -1;
60 theInterpolation = -1;
61 fCacheAngular.Put(0); //fix
62 }
G4Cache::Put
void Put(const value_type &val) const
Definition: G4Cache.hh:321
G4ParticleHPContEnergyAngular::fCacheAngular
G4Cache< G4ParticleHPContAngularPar * > fCacheAngular
Definition: G4ParticleHPContEnergyAngular.hh:105

References currentMeanEnergy, fCacheAngular, nEnergy, G4Cache< VALTYPE >::Put(), theAngular, theAngularRep, theInterpolation, and theTargetCode.

◆ ~G4ParticleHPContEnergyAngular()

G4ParticleHPContEnergyAngular::~G4ParticleHPContEnergyAngular ( )
inline

Definition at line 64 of file G4ParticleHPContEnergyAngular.hh.

65 {
66 if(theAngular!=0) delete [] theAngular;
67 if (fCacheAngular.Get() != 0) delete fCacheAngular.Get(); //fix
68 }
G4Cache::Get
value_type & Get() const
Definition: G4Cache.hh:315

References fCacheAngular, G4Cache< VALTYPE >::Get(), and theAngular.

Member Function Documentation

◆ ClearHistories()

void G4ParticleHPContEnergyAngular::ClearHistories ( )
virtual

Reimplemented from G4VParticleHPEnergyAngular.

Definition at line 128 of file G4ParticleHPContEnergyAngular.cc.

129{
130 if (theAngular!= NULL) {
131 for (G4int i = 0; i < nEnergy; i++) theAngular[i].ClearHistories();
132 }
133
134 // Added fCacheAngular ClearHistories() - this is the one actually used!
135 // Maybe theAngular does not even need ClearHistories()?
136 if (fCacheAngular.Get() != 0) fCacheAngular.Get()->ClearHistories();
137}
G4int
int G4int
Definition: G4Types.hh:85

References G4ParticleHPContAngularPar::ClearHistories(), ClearHistories(), fCacheAngular, G4Cache< VALTYPE >::Get(), nEnergy, and theAngular.

Referenced by ClearHistories().

◆ GetCMS()

G4ReactionProduct * G4VParticleHPEnergyAngular::GetCMS ( )
inlineinherited

Definition at line 94 of file G4VParticleHPEnergyAngular.hh.

95 {
96 *fCache.Get().theCMS = *fCache.Get().theProjectileRP
97 + *fCache.Get().theTarget;
98 return fCache.Get().theCMS;
99 }

References G4VParticleHPEnergyAngular::fCache.

◆ GetProjectileRP()

G4ReactionProduct * G4VParticleHPEnergyAngular::GetProjectileRP ( )
inlineinherited

Definition at line 89 of file G4VParticleHPEnergyAngular.hh.

90 {
91 return fCache.Get().theProjectileRP;
92 }

References G4VParticleHPEnergyAngular::fCache.

Referenced by G4ParticleHPNBodyPhaseSpace::GetEmax(), Sample(), G4ParticleHPDiscreteTwoBody::Sample(), and G4ParticleHPNBodyPhaseSpace::Sample().

◆ GetQValue()

G4double G4VParticleHPEnergyAngular::GetQValue ( )
inlineprotectedinherited

Definition at line 112 of file G4VParticleHPEnergyAngular.hh.

112{ return theQValue; }

References G4VParticleHPEnergyAngular::theQValue.

Referenced by G4ParticleHPNBodyPhaseSpace::GetEmax(), G4ParticleHPDiscreteTwoBody::Sample(), and G4ParticleHPIsotropic::Sample().

◆ GetTarget()

G4ReactionProduct * G4VParticleHPEnergyAngular::GetTarget ( void  )
inlineinherited

Definition at line 84 of file G4VParticleHPEnergyAngular.hh.

85 {
86 return fCache.Get().theTarget;
87 }

References G4VParticleHPEnergyAngular::fCache.

Referenced by G4ParticleHPNBodyPhaseSpace::GetEmax(), Sample(), G4ParticleHPDiscreteTwoBody::Sample(), and G4ParticleHPNBodyPhaseSpace::Sample().

◆ Init()

void G4ParticleHPContEnergyAngular::Init ( std::istream &  aDataFile)
inlinevirtual

Implements G4VParticleHPEnergyAngular.

Definition at line 70 of file G4ParticleHPContEnergyAngular.hh.

71 {
72 aDataFile >> theTargetCode >> theAngularRep >> theInterpolation >> nEnergy;
73 theAngular = new G4ParticleHPContAngularPar[nEnergy];
74 theManager.Init(aDataFile);
75 for(G4int i=0; i<nEnergy; i++)
76 {
77 theAngular[i].Init(aDataFile, theProjectile);
78 theAngular[i].SetInterpolation(theInterpolation);
79#ifndef PHP_AS_HP
80 if( i != 0 ) {
81 theAngular[i].PrepareTableInterpolation();
82 } else {
83 theAngular[i].PrepareTableInterpolation();
84 }
85#endif
86 }
87 }
G4InterpolationManager::Init
void Init(G4int aScheme, G4int aRange)
Definition: G4InterpolationManager.hh:84
G4ParticleHPContAngularPar::Init
void Init(std::istream &aDataFile, G4ParticleDefinition *projectile)
Definition: G4ParticleHPContAngularPar.cc:86
G4ParticleHPContAngularPar::SetInterpolation
void SetInterpolation(G4int theInterpolation)
Definition: G4ParticleHPContAngularPar.hh:143

References G4InterpolationManager::Init(), G4ParticleHPContAngularPar::Init(), nEnergy, G4ParticleHPContAngularPar::PrepareTableInterpolation(), G4ParticleHPContAngularPar::SetInterpolation(), theAngular, theAngularRep, theInterpolation, theManager, theProjectile, and theTargetCode.

◆ MeanEnergyOfThisInteraction()

G4double G4ParticleHPContEnergyAngular::MeanEnergyOfThisInteraction ( )
virtual

Implements G4VParticleHPEnergyAngular.

Definition at line 114 of file G4ParticleHPContEnergyAngular.cc.

115{
116 G4double result(0);
117 if (currentMeanEnergy.Get() < -1) {
118 throw G4HadronicException(__FILE__, __LINE__,
119 "G4ParticleHPContEnergyAngular: Logical error in Product class");
120 } else {
121 result = currentMeanEnergy.Get();
122 }
123 currentMeanEnergy.Put(-2);
124 return result;
125}
G4double
double G4double
Definition: G4Types.hh:83

References currentMeanEnergy, G4Cache< VALTYPE >::Get(), and G4Cache< VALTYPE >::Put().

Referenced by Sample().

◆ Sample()

G4ReactionProduct * G4ParticleHPContEnergyAngular::Sample ( G4double  anEnergy,
G4double  massCode,
G4double  mass 
)
virtual

Implements G4VParticleHPEnergyAngular.

Definition at line 37 of file G4ParticleHPContEnergyAngular.cc.

39{
40 G4ReactionProduct* result;
41 G4int i(0);
42 G4int it(0);
43 for (i = 0; i < nEnergy; i++) {
44 it = i;
45#ifdef PHP_AS_HP
46 if(theAngular[i].GetEnergy() > anEnergy) break;
47#else
48 if(theAngular[i].GetEnergy() >= anEnergy) break;
49#endif
50 }
51
52 if (std::getenv("G4PHPTEST") )
53 G4cout << i << " G4ParticleHPContEnergyAngular dataE "
54 << theAngular[i].GetEnergy()
55 << " > " << anEnergy << " it_theAngular " << it << " interpolation "
56 << theInterpolation << G4endl; //GDEB
57
58 G4double targetMass = GetTarget()->GetMass();
59 if (it == 0) {
60 theAngular[0].SetTarget(GetTarget());
61 theAngular[0].SetTargetCode(theTargetCode);
62 theAngular[0].SetPrimary(GetProjectileRP());
63 result = theAngular[0].Sample(anEnergy, massCode, targetMass,
64 theAngularRep, theInterpolation);
65 currentMeanEnergy.Put(theAngular[0].MeanEnergyOfThisInteraction() );
66
67 } else {
68 // interpolation through alternating sampling. This needs improvement @@@
69 // This is the cause of the He3 problem !!!!!!!!
70 // See to it, if you can improve this.
71 //G4double random = G4UniformRand();
72 //G4double deltaE = theAngular[it].GetEnergy()-theAngular[it-1].GetEnergy();
73 //G4double offset = theAngular[it].GetEnergy()-anEnergy;
74 //if(random<offset/deltaE) it--;
75 //--- create new
76 // if (theManager.GetScheme(0) != LINLIN) {
77 // // asserted in G4ParticleHPContEnergyAngular::init there is only one range
78#ifdef PHP_AS_HP
79 theAngular[it].SetTarget(GetTarget());
80 theAngular[it].SetTargetCode(theTargetCode);
81 theAngular[it].SetPrimary(GetProjectileRP());
82 result = theAngular[it].Sample(anEnergy, massCode, targetMass,
83 theAngularRep, theInterpolation);
84 currentMeanEnergy.Put(theAngular[it].MeanEnergyOfThisInteraction() );
85#else
86 if (std::getenv("G4PHPTEST") )
87 G4cout << i << " G4ParticleHPContEnergyAngular To BUILDBYINTERPOLATION "
88 << it << " : " << theAngular[it].GetEnergy()<< " , "
89 << theAngular[it].GetNEnergies() << " " << it-1 << " : "
90 << theAngular[it-1].GetEnergy()<< " : "
91 << theAngular[it-1].GetNEnergies() << G4endl; //GDEB
92
93 if (fCacheAngular.Get() == NULL) {
94 G4ParticleHPContAngularPar* angpar = new G4ParticleHPContAngularPar(theProjectile);
95 fCacheAngular.Put(angpar);
96 }
97 fCacheAngular.Get()->SetInterpolation(theInterpolation);
98 fCacheAngular.Get()->BuildByInterpolation(anEnergy, theManager.GetScheme(0),
99 (theAngular[it-1]), (theAngular[it]) );
100 fCacheAngular.Get()->SetTarget(GetTarget());
101 fCacheAngular.Get()->SetTargetCode(theTargetCode);
102 fCacheAngular.Get()->SetPrimary(GetProjectileRP());
103
104 result = fCacheAngular.Get()->Sample(anEnergy, massCode, targetMass,
105 theAngularRep, theInterpolation);
106 currentMeanEnergy.Put(fCacheAngular.Get()->MeanEnergyOfThisInteraction() );
107#endif
108 } // end (it != 0) branch
109
110 return result;
111}
G4endl
#define G4endl
Definition: G4ios.hh:57
G4cout
G4GLOB_DLL std::ostream G4cout
G4InterpolationManager::GetScheme
G4InterpolationScheme GetScheme(G4int index) const
Definition: G4InterpolationManager.hh:131
G4ParticleHPContAngularPar::Sample
G4ReactionProduct * Sample(G4double anEnergy, G4double massCode, G4double mass, G4int angularRep, G4int interpol)
Definition: G4ParticleHPContAngularPar.cc:110
G4ParticleHPContAngularPar::SetPrimary
void SetPrimary(G4ReactionProduct *aPrimary)
Definition: G4ParticleHPContAngularPar.hh:128
G4ParticleHPContAngularPar::SetTarget
void SetTarget(G4ReactionProduct *aTarget)
Definition: G4ParticleHPContAngularPar.hh:133
G4ParticleHPContAngularPar::SetTargetCode
void SetTargetCode(G4double aTargetCode)
Definition: G4ParticleHPContAngularPar.hh:138
G4ParticleHPContAngularPar::BuildByInterpolation
void BuildByInterpolation(G4double anEnergy, G4InterpolationScheme aScheme, G4ParticleHPContAngularPar &store1, G4ParticleHPContAngularPar &store2)
Definition: G4ParticleHPContAngularPar.cc:714
G4ReactionProduct::GetMass
G4double GetMass() const
Definition: G4ReactionProduct.hh:150

References G4ParticleHPContAngularPar::BuildByInterpolation(), currentMeanEnergy, fCacheAngular, G4cout, G4endl, G4Cache< VALTYPE >::Get(), G4ParticleHPContAngularPar::GetEnergy(), G4ReactionProduct::GetMass(), G4ParticleHPContAngularPar::GetNEnergies(), G4VParticleHPEnergyAngular::GetProjectileRP(), G4InterpolationManager::GetScheme(), G4VParticleHPEnergyAngular::GetTarget(), G4ParticleHPContAngularPar::MeanEnergyOfThisInteraction(), MeanEnergyOfThisInteraction(), nEnergy, G4Cache< VALTYPE >::Put(), G4ParticleHPContAngularPar::Sample(), G4ParticleHPContAngularPar::SetInterpolation(), G4ParticleHPContAngularPar::SetPrimary(), G4ParticleHPContAngularPar::SetTarget(), G4ParticleHPContAngularPar::SetTargetCode(), theAngular, theAngularRep, theInterpolation, theManager, theProjectile, and theTargetCode.

◆ SetProjectileRP()

void G4VParticleHPEnergyAngular::SetProjectileRP ( G4ReactionProduct *  aIncidentParticleRP)
inlineinherited

Definition at line 74 of file G4VParticleHPEnergyAngular.hh.

75 {
76 fCache.Get().theProjectileRP = aIncidentParticleRP;
77 }

References G4VParticleHPEnergyAngular::fCache.

Referenced by G4ParticleHPInelasticBaseFS::BaseApply(), and G4ParticleHPProduct::Sample().

◆ SetQValue()

void G4VParticleHPEnergyAngular::SetQValue ( G4double  aValue)
inlineinherited

Definition at line 101 of file G4VParticleHPEnergyAngular.hh.

102 {
103 theQValue = aValue;
104 }

References G4VParticleHPEnergyAngular::theQValue.

Referenced by G4ParticleHPInelasticBaseFS::BaseApply(), and G4ParticleHPProduct::Init().

◆ SetTarget()

void G4VParticleHPEnergyAngular::SetTarget ( G4ReactionProduct *  aTarget)
inlineinherited

Definition at line 79 of file G4VParticleHPEnergyAngular.hh.

80 {
81 fCache.Get().theTarget = aTarget;
82 }

References G4VParticleHPEnergyAngular::fCache.

Referenced by G4ParticleHPInelasticBaseFS::BaseApply(), and G4ParticleHPProduct::Sample().

Field Documentation

◆ currentMeanEnergy

G4Cache<G4double> G4ParticleHPContEnergyAngular::currentMeanEnergy
private

Definition at line 104 of file G4ParticleHPContEnergyAngular.hh.

Referenced by G4ParticleHPContEnergyAngular(), MeanEnergyOfThisInteraction(), and Sample().

◆ fCache

G4Cache<toBeCached> G4VParticleHPEnergyAngular::fCache
privateinherited

Definition at line 118 of file G4VParticleHPEnergyAngular.hh.

Referenced by G4VParticleHPEnergyAngular::G4VParticleHPEnergyAngular(), G4VParticleHPEnergyAngular::GetCMS(), G4VParticleHPEnergyAngular::GetProjectileRP(), G4VParticleHPEnergyAngular::GetTarget(), G4VParticleHPEnergyAngular::SetProjectileRP(), and G4VParticleHPEnergyAngular::SetTarget().

◆ fCacheAngular

G4Cache<G4ParticleHPContAngularPar*> G4ParticleHPContEnergyAngular::fCacheAngular
private

Definition at line 105 of file G4ParticleHPContEnergyAngular.hh.

Referenced by ClearHistories(), G4ParticleHPContEnergyAngular(), Sample(), and ~G4ParticleHPContEnergyAngular().

◆ nEnergy

G4int G4ParticleHPContEnergyAngular::nEnergy
private

Definition at line 97 of file G4ParticleHPContEnergyAngular.hh.

Referenced by ClearHistories(), G4ParticleHPContEnergyAngular(), Init(), and Sample().

◆ theAngular

G4ParticleHPContAngularPar* G4ParticleHPContEnergyAngular::theAngular
private

Definition at line 102 of file G4ParticleHPContEnergyAngular.hh.

Referenced by ClearHistories(), G4ParticleHPContEnergyAngular(), Init(), Sample(), and ~G4ParticleHPContEnergyAngular().

◆ theAngularRep

G4int G4ParticleHPContEnergyAngular::theAngularRep
private

Definition at line 96 of file G4ParticleHPContEnergyAngular.hh.

Referenced by G4ParticleHPContEnergyAngular(), Init(), and Sample().

◆ theInterpolation

G4int G4ParticleHPContEnergyAngular::theInterpolation
private

Definition at line 99 of file G4ParticleHPContEnergyAngular.hh.

Referenced by G4ParticleHPContEnergyAngular(), Init(), and Sample().

◆ theManager

G4InterpolationManager G4ParticleHPContEnergyAngular::theManager
private

Definition at line 101 of file G4ParticleHPContEnergyAngular.hh.

Referenced by Init(), and Sample().

◆ theProjectile

G4ParticleDefinition* G4ParticleHPContEnergyAngular::theProjectile
private

Definition at line 106 of file G4ParticleHPContEnergyAngular.hh.

Referenced by Init(), and Sample().

◆ theQValue

G4double G4VParticleHPEnergyAngular::theQValue
privateinherited

Definition at line 116 of file G4VParticleHPEnergyAngular.hh.

Referenced by G4VParticleHPEnergyAngular::G4VParticleHPEnergyAngular(), G4VParticleHPEnergyAngular::GetQValue(), and G4VParticleHPEnergyAngular::SetQValue().

◆ theTargetCode

G4double G4ParticleHPContEnergyAngular::theTargetCode
private

Definition at line 95 of file G4ParticleHPContEnergyAngular.hh.

Referenced by G4ParticleHPContEnergyAngular(), Init(), and Sample().

The documentation for this class was generated from the following files:
Generated by doxygen 1.9.3