Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
G4ParticleHPFinalState Class Referenceabstract

#include <G4ParticleHPFinalState.hh>

Inheritance diagram for G4ParticleHPFinalState:
Inheritance graph
[legend]
Collaboration diagram for G4ParticleHPFinalState:
Collaboration graph
[legend]

Public Member Functions

 G4ParticleHPFinalState ()
 
virtual ~G4ParticleHPFinalState ()
 
void Init (G4double A, G4double Z, G4String &dirName, G4String &aFSType, G4ParticleDefinition *projectile)
 
virtual void Init (G4double A, G4double Z, G4int M, G4String &dirName, G4String &aFSType, G4ParticleDefinition *)=0
 
virtual G4HadFinalState * ApplyYourself (const G4HadProjectile &)
 
virtual G4ParticleHPFinalState * New ()=0
 
G4bool HasXsec ()
 
G4bool HasFSData ()
 
G4bool HasAnyData ()
 
virtual G4double GetXsec (G4double)
 
virtual G4ParticleHPVector * GetXsec ()
 
void SetA_Z (G4double anA, G4double aZ, G4int aM=0)
 
G4double GetZ ()
 
G4double GetN ()
 
G4double GetA ()
 
G4int GetM ()
 
void SetAZMs (G4double anA, G4double aZ, G4int aM, G4ParticleHPDataUsed used)
 
void SetProjectile (G4ParticleDefinition *projectile)
 

Protected Member Functions

void adjust_final_state (G4LorentzVector)
 
G4bool DoNotAdjustFinalState ()
 

Protected Attributes

G4bool hasXsec
 
G4bool hasFSData
 
G4bool hasAnyData
 
G4ParticleHPNames theNames
 
G4Cache< G4HadFinalState * > theResult
 
G4double theBaseA
 
G4double theBaseZ
 
G4int theBaseM
 
G4int theNDLDataZ
 
G4int theNDLDataA
 
G4int theNDLDataM
 
G4ParticleDefinition * theProjectile
 

Detailed Description

Definition at line 47 of file G4ParticleHPFinalState.hh.

Constructor & Destructor Documentation

G4ParticleHPFinalState::G4ParticleHPFinalState ( )
inline

Definition at line 51 of file G4ParticleHPFinalState.hh.

52  {
53  hasFSData = true;
54  hasXsec = true;
55  hasAnyData = true;
56  theBaseZ = 0;
57  theBaseA = 0;
58  theBaseM = 0;
59 
60  theNDLDataZ = 0;
61  theNDLDataA = 0;
62  theNDLDataM = 0;
63 
64  adjustResult = true;
65  if ( getenv( "G4PHP_DO_NOT_ADJUST_FINAL_STATE" ) ) adjustResult = false;
66 
67  theProjectile = G4Neutron::Neutron();
68 
69  theResult.Put( NULL );
70 
71  };
G4ParticleHPFinalState::theResult
G4Cache< G4HadFinalState * > theResult
Definition: G4ParticleHPFinalState.hh:117
G4ParticleHPFinalState::theProjectile
G4ParticleDefinition * theProjectile
Definition: G4ParticleHPFinalState.hh:136
G4Neutron::Neutron
static G4Neutron * Neutron()
Definition: G4Neutron.cc:104
G4Cache::Put
void Put(const value_type &val) const
Definition: G4Cache.hh:286

Here is the call graph for this function:

virtual G4ParticleHPFinalState::~G4ParticleHPFinalState ( )
inlinevirtual

Definition at line 73 of file G4ParticleHPFinalState.hh.

73 {};

Member Function Documentation

void G4ParticleHPFinalState::adjust_final_state ( G4LorentzVector  init_4p_lab)
protected

Definition at line 46 of file G4ParticleHPFinalState.cc.

47 {
48 
49  G4double minimum_energy = 1*keV;
50 
51  if ( G4ParticleHPManager::GetInstance()->GetDoNotAdjustFinalState() ) return;
52 
53  G4int nSecondaries = theResult.Get()->GetNumberOfSecondaries();
54 
55  G4int sum_Z = 0;
56  G4int sum_A = 0;
57  G4int max_SecZ = 0;
58  G4int max_SecA = 0;
59  G4int imaxA = -1;
60  for ( int i = 0 ; i < nSecondaries ; i++ )
61  {
62  //G4cout << "G4ParticleHPFinalState::adjust_final_state theResult.Get()->GetSecondary( i )->GetParticle()->GetDefinition()->GetParticleName() = " << theResult.Get()->GetSecondary( i )->GetParticle()->GetDefinition()->GetParticleName() << G4endl;
63  sum_Z += theResult.Get()->GetSecondary( i )->GetParticle()->GetDefinition()->GetAtomicNumber();
64  max_SecZ = std::max ( max_SecZ , theResult.Get()->GetSecondary( i )->GetParticle()->GetDefinition()->GetAtomicNumber() );
65  sum_A += theResult.Get()->GetSecondary( i )->GetParticle()->GetDefinition()->GetAtomicMass();
66  max_SecA = std::max ( max_SecA , theResult.Get()->GetSecondary( i )->GetParticle()->GetDefinition()->GetAtomicMass() );
67  if ( theResult.Get()->GetSecondary( i )->GetParticle()->GetDefinition()->GetAtomicMass() == max_SecA ) imaxA = i;
68 #ifdef G4PHPDEBUG
69  if( getenv("G4ParticleHPDebug")) G4cout << "G4ParticleHPFinalState::adjust_final_stat SECO " << i << " " <<theResult.Get()->GetSecondary( i )->GetParticle()->GetDefinition()->GetParticleName() << G4endl;
70 #endif
71 
72  }
73 
74  G4ParticleDefinition* resi_pd = NULL;
75 
76  G4double baseZNew = theBaseZ;
77  G4double baseANew = theBaseA;
78  if( theProjectile == G4Neutron::Neutron() ) {
79  baseANew ++;
80  } else if( theProjectile == G4Proton::Proton() ) {
81  baseZNew ++;
82  baseANew ++;
83  } else if( theProjectile == G4Deuteron::Deuteron() ) {
84  baseZNew ++;
85  baseANew += 2;
86  } else if( theProjectile == G4Triton::Triton() ) {
87  baseZNew ++;
88  baseANew += 3;
89  } else if( theProjectile == G4He3::He3() ) {
90  baseZNew += 2;
91  baseANew += 3;
92  } else if( theProjectile == G4Alpha::Alpha() ) {
93  baseZNew += 2;
94  baseANew += 4;
95  }
96 
97 #ifdef G4PHPDEBUG
98  if( getenv("G4ParticleHPDebug")) G4cout << "G4ParticleHPFinalState::adjust_final_stat BaseZ " << baseZNew << " BaseA " << baseANew << " sum_Z " << sum_Z << " sum_A " << sum_A << G4endl;
99 #endif
100 
101  G4bool needOneMoreSec = false;
102  G4ParticleDefinition* oneMoreSec_pd = NULL;
103  if ( (int)(baseZNew - sum_Z) == 0 && (int)(baseANew - sum_A) == 0 )
104  {
105  //All secondaries are already created;
106  resi_pd = theResult.Get()->GetSecondary( imaxA )->GetParticle()->GetDefinition();
107  }
108  else
109  {
110  if ( max_SecA >= int(baseANew - sum_A) )
111  {
112  //Most heavy secondary is interpreted as residual
113  resi_pd = theResult.Get()->GetSecondary( imaxA )->GetParticle()->GetDefinition();
114  needOneMoreSec = true;
115  }
116  else
117  {
118  //creation of residual is requierd
119  resi_pd = G4IonTable::GetIonTable()->GetIon ( int(baseZNew - sum_Z) , (int)(baseANew - sum_A) , 0.0 );
120  }
121 
122  if ( needOneMoreSec )
123  {
124  if ( int(baseZNew - sum_Z) == 0 && (int)(baseANew - sum_A) > 0 )
125  {
126  //In this case, one neutron is added to secondaries
127  if ( int(baseANew - sum_A) > 1 ) G4cout << "More than one neutron is required for the balance of baryon number!" << G4endl;
128  oneMoreSec_pd = G4Neutron::Neutron();
129  }
130  else
131  {
132 #ifdef G4PHPDEBUG
133  if( getenv("G4ParticleHPDebug")) G4cout << this << "G4ParticleHPFinalState oneMoreSec_pd Z " << baseZNew << " - " << sum_Z << " A " << baseANew << " - " << sum_A << " projectile " << theProjectile->GetParticleName() << G4endl;
134 #endif
135  oneMoreSec_pd = G4IonTable::GetIonTable()->GetIon ( int(baseZNew - sum_Z) , (int)(baseANew - sum_A) , 0.0 );
136  if( !oneMoreSec_pd ) {
137  G4cerr << this << "G4ParticleHPFinalState oneMoreSec_pd Z " << baseZNew << " - " << sum_Z << " A " << baseANew << " - " << sum_A << " projectile " << theProjectile->GetParticleName() << G4endl;
138  G4Exception("G4ParticleHPFinalState:adjust_final_state",
139  "Warning",
140  JustWarning,
141  "No adjustment will be done!");
142  return;
143  }
144  }
145  }
146 
147  if ( resi_pd == NULL )
148  {
149  // theNDLDataZ,A has the Z and A of used NDL file
150  G4double ndlZNew = theNDLDataZ;
151  G4double ndlANew = theNDLDataA;
152  if( theProjectile == G4Neutron::Neutron() ) {
153  ndlANew ++;
154  } else if( theProjectile == G4Proton::Proton() ) {
155  ndlZNew ++;
156  ndlANew ++;
157  } else if( theProjectile == G4Deuteron::Deuteron() ) {
158  ndlZNew ++;
159  ndlANew += 2;
160  } else if( theProjectile == G4Triton::Triton() ) {
161  ndlZNew ++;
162  ndlANew += 3;
163  } else if( theProjectile == G4He3::He3() ) {
164  ndlZNew += 2;
165  ndlANew += 3;
166  } else if( theProjectile == G4Alpha::Alpha() ) {
167  ndlZNew += 2;
168  ndlANew += 4;
169  }
170  // theNDLDataZ,A has the Z and A of used NDL file
171  if ( (int)(ndlZNew - sum_Z) == 0 && (int)(ndlANew - sum_A) == 0 )
172  {
173  G4int dif_Z = ( int ) ( theNDLDataZ - theBaseZ );
174  G4int dif_A = ( int ) ( theNDLDataA - theBaseA );
175  resi_pd = G4IonTable::GetIonTable()->GetIon ( max_SecZ - dif_Z , max_SecA - dif_A , 0.0 );
176  if( !resi_pd ) {
177  G4cerr << "G4ParticleHPFinalState resi_pd Z " << max_SecZ << " - " << dif_Z << " A " << max_SecA << " - " << dif_A << " projectile " << theProjectile->GetParticleName() << G4endl;
178  G4Exception("G4ParticleHPFinalState:adjust_final_state",
179  "Warning",
180  JustWarning,
181  "No adjustment will be done!");
182  return;
183  }
184 
185  for ( int i = 0 ; i < nSecondaries ; i++ )
186  {
187  if ( theResult.Get()->GetSecondary( i )->GetParticle()->GetDefinition()->GetAtomicNumber() == max_SecZ
188  && theResult.Get()->GetSecondary( i )->GetParticle()->GetDefinition()->GetAtomicMass() == max_SecA )
189  {
190  G4ThreeVector p = theResult.Get()->GetSecondary( i )->GetParticle()->GetMomentum();
191  p = p * resi_pd->GetPDGMass()/ G4IonTable::GetIonTable()->GetIon ( max_SecZ , max_SecA , 0.0 )->GetPDGMass();
192  theResult.Get()->GetSecondary( i )->GetParticle()->SetDefinition( resi_pd );
193  theResult.Get()->GetSecondary( i )->GetParticle()->SetMomentum( p );
194  }
195  }
196  }
197  }
198  }
199 
200 
201  G4LorentzVector secs_4p_lab( 0.0 );
202 
203  G4int n_sec = theResult.Get()->GetNumberOfSecondaries();
204  G4double fast = 0;
205  G4double slow = 1;
206  G4int ifast = 0;
207  G4int islow = 0;
208  G4int ires = -1;
209 
210  for ( G4int i = 0 ; i < n_sec ; i++ )
211  {
212 
213  //G4cout << "HP_DB " << i
214  // << " " << theResult.GetSecondary( i )->GetParticle()->GetDefinition()->GetParticleName()
215  // << " 4p " << theResult.GetSecondary( i )->GetParticle()->Get4Momentum()
216  // << " ke " << theResult.GetSecondary( i )->GetParticle()->Get4Momentum().e() - theResult.GetSecondary( i )->GetParticle()->GetDefinition()->GetPDGMass()
217  // << G4endl;
218 
219  secs_4p_lab += theResult.Get()->GetSecondary( i )->GetParticle()->Get4Momentum();
220 
221  G4double beta = 0;
222  if ( theResult.Get()->GetSecondary( i )->GetParticle()->GetDefinition() != G4Gamma::Gamma() )
223  {
224  beta = theResult.Get()->GetSecondary( i )->GetParticle()->Get4Momentum().beta();
225  }
226  else
227  {
228  beta = 1;
229  }
230 
231  if ( theResult.Get()->GetSecondary( i )->GetParticle()->GetDefinition() == resi_pd ) ires = i;
232 
233  if ( slow > beta && beta != 0 )
234  {
235  slow = beta;
236  islow = i;
237  }
238 
239  if ( fast <= beta )
240  {
241  if ( fast != 1 )
242  {
243  fast = beta;
244  ifast = i;
245  }
246  else
247  {
248 // fast is already photon then check E
249  G4double e = theResult.Get()->GetSecondary( i )->GetParticle()->Get4Momentum().e();
250  if ( e > theResult.Get()->GetSecondary( ifast )->GetParticle()->Get4Momentum().e() )
251  {
252 // among photons, the highest E becomes the fastest
253  ifast = i;
254  }
255  }
256  }
257  }
258 
259 
260  G4LorentzVector dif_4p = init_4p_lab - secs_4p_lab;
261 
262  //G4cout << "HP_DB dif_4p " << init_4p_lab - secs_4p_lab << G4endl;
263  //G4cout << "HP_DB dif_3p mag " << ( dif_4p.v() ).mag() << G4endl;
264  //G4cout << "HP_DB dif_e " << dif_4p.e() - ( dif_4p.v() ).mag()<< G4endl;
265 
266  G4LorentzVector p4(0);
267  if ( ires == -1 )
268  {
269 // Create and Add Residual Nucleus
270  ires = nSecondaries;
271  nSecondaries += 1;
272 
273  G4DynamicParticle* res = new G4DynamicParticle ( resi_pd , dif_4p.v() );
274  theResult.Get()->AddSecondary ( res );
275 
276  p4 = res->Get4Momentum();
277  if ( slow > p4.beta() )
278  {
279  slow = p4.beta();
280  islow = ires;
281  }
282  dif_4p = init_4p_lab - ( secs_4p_lab + p4 );
283  }
284 
285  if ( needOneMoreSec && oneMoreSec_pd)
286  //
287  // fca: this is not a fix, this is a crash avoidance...
288  // fca: the baryon number is still wrong, most probably because it
289  // fca: should have been decreased, but since we could not create a particle
290  // fca: we just do not add it
291  //
292  {
293  nSecondaries += 1;
294  G4DynamicParticle* one = new G4DynamicParticle ( oneMoreSec_pd , dif_4p.v() );
295  theResult.Get()->AddSecondary ( one );
296  p4 = one->Get4Momentum();
297  if ( slow > p4.beta() )
298  {
299  slow = p4.beta();
300  islow = nSecondaries-1; //Because the first is 0th, so the last becomes "nSecondaries-1"
301  }
302  dif_4p = init_4p_lab - ( secs_4p_lab + p4 );
303  }
304 
305  //Which is bigger dif_p or dif_e
306 
307  if ( dif_4p.v().mag() < std::abs( dif_4p.e() ) )
308  {
309 
310  // Adjust p
311  //if ( dif_4p.v().mag() < 1*MeV )
312  if ( minimum_energy < dif_4p.v().mag() && dif_4p.v().mag() < 1*MeV )
313  {
314 
315  nSecondaries += 1;
316  theResult.Get()->AddSecondary ( new G4DynamicParticle ( G4Gamma::Gamma() , dif_4p.v() ) );
317 
318  }
319  else
320  {
321  //G4cout << "HP_DB Difference in dif_p is too large (>1MeV) or too small(<1keV) to adjust, so that give up tuning" << G4endl;
322  }
323 
324  }
325  else
326  {
327 
328  // dif_p > dif_e
329  // at first momentum
330  // Move residual momentum
331 
332  p4 = theResult.Get()->GetSecondary( ires )->GetParticle()->Get4Momentum();
333  theResult.Get()->GetSecondary( ires )->GetParticle()->SetMomentum( p4.v() + dif_4p.v() );
334  dif_4p = init_4p_lab - ( secs_4p_lab - p4 + theResult.Get()->GetSecondary( ires )->GetParticle()->Get4Momentum() );
335 
336  //G4cout << "HP_DB new residual kinetic energy " << theResult.GetSecondary( ires )->GetParticle()->GetKineticEnergy() << G4endl;
337 
338  }
339 
340  G4double dif_e = dif_4p.e() - ( dif_4p.v() ).mag();
341  //G4cout << "HP_DB dif_e " << dif_e << G4endl;
342 
343  if ( dif_e > 0 )
344  {
345 
346 // create 2 gamma
347 
348  nSecondaries += 2;
349  G4double e1 = ( dif_4p.e() -dif_4p.v().mag() ) / 2;
350 
351  if ( minimum_energy < e1 )
352  {
353  G4double costh = 2.*G4UniformRand()-1.;
354  G4double phi = twopi*G4UniformRand();
355  G4ThreeVector dir( std::sin(std::acos(costh))*std::cos(phi),
356  std::sin(std::acos(costh))*std::sin(phi),
357  costh);
358  theResult.Get()->AddSecondary ( new G4DynamicParticle ( G4Gamma::Gamma() , e1*dir ) );
359  theResult.Get()->AddSecondary ( new G4DynamicParticle ( G4Gamma::Gamma() , -e1*dir ) );
360  }
361  else
362  {
363  //G4cout << "HP_DB Difference is too small(<1keV) to adjust, so that neglect it" << G4endl;
364  }
365 
366  }
367  else //dif_e < 0
368  {
369 
370 // At first reduce KE of the fastest secondary;
371  G4double ke0 = theResult.Get()->GetSecondary( ifast )->GetParticle()->GetKineticEnergy();
372  G4ThreeVector p0 = theResult.Get()->GetSecondary( ifast )->GetParticle()->GetMomentum();
373  G4ThreeVector dir = ( theResult.Get()->GetSecondary( ifast )->GetParticle()->GetMomentum() ).unit();
374 
375  //G4cout << "HP_DB ifast " << ifast << " ke0 " << ke0 << G4endl;
376 
377  if ( ke0 + dif_e > 0 )
378  {
379  theResult.Get()->GetSecondary( ifast )->GetParticle()->SetKineticEnergy( ke0 + dif_e );
380  G4ThreeVector dp = p0 - theResult.Get()->GetSecondary( ifast )->GetParticle()->GetMomentum();
381 
382  G4ThreeVector p = theResult.Get()->GetSecondary( islow )->GetParticle()->GetMomentum();
383  //theResult.GetSecondary( islow )->GetParticle()->SetMomentum( p - dif_e*dir );
384  theResult.Get()->GetSecondary( islow )->GetParticle()->SetMomentum( p + dp );
385  }
386  else
387  {
388  //G4cout << "HP_DB Difference in dif_e too large ( <0MeV ) to adjust, so that give up tuning" << G4endl;
389  }
390 
391  }
392 
393 }
G4ParticleHPManager::GetInstance
static G4ParticleHPManager * GetInstance()
Definition: G4ParticleHPManager.hh:56
G4ParticleHPFinalState::theResult
G4Cache< G4HadFinalState * > theResult
Definition: G4ParticleHPFinalState.hh:117
G4DynamicParticle::SetMomentum
void SetMomentum(const G4ThreeVector &momentum)
Definition: G4DynamicParticle.cc:371
G4HadFinalState::GetSecondary
G4HadSecondary * GetSecondary(size_t i)
Definition: G4HadFinalState.cc:79
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
p
const char * p
Definition: xmltok.h:285
G4IonTable::GetIon
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
Definition: G4IonTable.cc:503
G4Cache::Get
value_type & Get() const
Definition: G4Cache.hh:282
CLHEP::HepLorentzVector::v
Hep3Vector v() const
G4DynamicParticle::GetDefinition
G4ParticleDefinition * GetDefinition() const
JustWarning
Definition: G4ExceptionSeverity.hh:64
G4int
int G4int
Definition: G4Types.hh:78
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition: G4ParticleDefinition.hh:120
G4ParticleDefinition::GetAtomicNumber
G4int GetAtomicNumber() const
G4ParticleHPFinalState::theProjectile
G4ParticleDefinition * theProjectile
Definition: G4ParticleHPFinalState.hh:136
twopi
static constexpr double twopi
Definition: G4SIunits.hh:76
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:97
G4cout
G4GLOB_DLL std::ostream G4cout
G4bool
bool G4bool
Definition: G4Types.hh:79
int
typedef int(XMLCALL *XML_NotStandaloneHandler)(void *userData)
G4Triton::Triton
static G4Triton * Triton()
Definition: G4Triton.cc:95
G4Proton::Proton
static G4Proton * Proton()
Definition: G4Proton.cc:93
G4Neutron::Neutron
static G4Neutron * Neutron()
Definition: G4Neutron.cc:104
G4Gamma::Gamma
static G4Gamma * Gamma()
Definition: G4Gamma.cc:86
G4Deuteron::Deuteron
static G4Deuteron * Deuteron()
Definition: G4Deuteron.cc:94
G4ParticleDefinition::GetAtomicMass
G4int GetAtomicMass() const
G4DynamicParticle::SetKineticEnergy
void SetKineticEnergy(G4double aEnergy)
G4DynamicParticle::Get4Momentum
G4LorentzVector Get4Momentum() const
G4IonTable::GetIonTable
static G4IonTable * GetIonTable()
Definition: G4IonTable.hh:78
G4Exception
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
Definition: G4Exception.cc:41
G4ParticleDefinition::GetPDGMass
G4double GetPDGMass() const
Definition: G4ParticleDefinition.hh:122
G4INCL::Math::max
T max(const T t1, const T t2)
brief Return the largest of the two arguments
Definition: G4INCLGlobals.hh:112
G4HadSecondary::GetParticle
G4DynamicParticle * GetParticle()
Definition: G4HadSecondary.hh:42
G4endl
#define G4endl
Definition: G4ios.hh:61
MeV
static constexpr double MeV
Definition: G4SIunits.hh:214
G4Alpha::Alpha
static G4Alpha * Alpha()
Definition: G4Alpha.cc:89
G4HadFinalState::AddSecondary
void AddSecondary(G4DynamicParticle *aP, G4int mod=-1)
Definition: G4HadFinalState.hh:61
G4double
double G4double
Definition: G4Types.hh:76
G4DynamicParticle::SetDefinition
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
Definition: G4DynamicParticle.cc:305
G4He3::He3
static G4He3 * He3()
Definition: G4He3.cc:94
CLHEP::Hep3Vector::mag
double mag() const
keV
static constexpr double keV
Definition: G4SIunits.hh:216
G4HadFinalState::GetNumberOfSecondaries
G4int GetNumberOfSecondaries() const
Definition: G4HadFinalState.hh:51
G4DynamicParticle::GetMomentum
G4ThreeVector GetMomentum() const
G4cerr
G4GLOB_DLL std::ostream G4cerr

Here is the call graph for this function:

Here is the caller graph for this function:

virtual G4HadFinalState* G4ParticleHPFinalState::ApplyYourself ( const G4HadProjectile &  )
inlinevirtual

Reimplemented in G4FissionLibrary, G4ParticleHPInelasticCompFS, G4ParticleHPInelasticBaseFS, G4ParticleHPElasticFS, G4ParticleHPCaptureFS, G4ParticleHPFissionFS, G4ParticleHP2AInelasticFS, G4ParticleHP3NAInelasticFS, G4ParticleHP3NInelasticFS, G4ParticleHP3NPInelasticFS, G4ParticleHP4NInelasticFS, G4ParticleHPAInelasticFS, G4ParticleHP2N2AInelasticFS, G4ParticleHP2NAInelasticFS, G4ParticleHP2NDInelasticFS, G4ParticleHP2NInelasticFS, G4ParticleHP2NPInelasticFS, G4ParticleHP2PInelasticFS, G4ParticleHP3AInelasticFS, G4ParticleHPDAInelasticFS, G4ParticleHPDInelasticFS, G4ParticleHPHe3InelasticFS, G4ParticleHPN2AInelasticFS, G4ParticleHPN2PInelasticFS, G4ParticleHPN3AInelasticFS, G4ParticleHPNAInelasticFS, G4ParticleHPND2AInelasticFS, G4ParticleHPNDInelasticFS, G4ParticleHPNHe3InelasticFS, G4ParticleHPNInelasticFS, G4ParticleHPNPAInelasticFS, G4ParticleHPNPInelasticFS, G4ParticleHPNT2AInelasticFS, G4ParticleHPNTInelasticFS, G4ParticleHPNXInelasticFS, G4ParticleHPPAInelasticFS, G4ParticleHPPDInelasticFS, G4ParticleHPPInelasticFS, G4ParticleHPPTInelasticFS, G4ParticleHPTInelasticFS, G4ParticleHPD2AInelasticFS, and G4ParticleHPT2AInelasticFS.

Definition at line 78 of file G4ParticleHPFinalState.hh.

79  {
80  throw G4HadronicException(__FILE__, __LINE__, "G4HadFinalState * ApplyYourself(const G4HadProjectile & theTrack) needs implementation");
81  return 0;
82  };

Here is the caller graph for this function:

G4bool G4ParticleHPFinalState::DoNotAdjustFinalState ( )
inlineprotected

Definition at line 127 of file G4ParticleHPFinalState.hh.

127 { return !G4ParticleHPManager::GetInstance()->GetDoNotAdjustFinalState(); };
G4ParticleHPManager::GetInstance
static G4ParticleHPManager * GetInstance()
Definition: G4ParticleHPManager.hh:56

Here is the caller graph for this function:

G4double G4ParticleHPFinalState::GetA ( void  )
inline

Definition at line 99 of file G4ParticleHPFinalState.hh.

99 { return theBaseA; };
G4int G4ParticleHPFinalState::GetM ( )
inline

Definition at line 100 of file G4ParticleHPFinalState.hh.

100 { return theBaseM; };

Here is the caller graph for this function:

G4double G4ParticleHPFinalState::GetN ( )
inline

Definition at line 98 of file G4ParticleHPFinalState.hh.

98 { return theBaseA; };

Here is the caller graph for this function:

virtual G4double G4ParticleHPFinalState::GetXsec ( G4double  )
inlinevirtual

Reimplemented in G4ParticleHPInelasticCompFS, G4ParticleHPInelasticBaseFS, and G4ParticleHPFissionBaseFS.

Definition at line 93 of file G4ParticleHPFinalState.hh.

93 { return 0; };

Here is the caller graph for this function:

virtual G4ParticleHPVector* G4ParticleHPFinalState::GetXsec ( )
inlinevirtual

Reimplemented in G4ParticleHPInelasticCompFS, G4ParticleHPInelasticBaseFS, and G4ParticleHPFissionBaseFS.

Definition at line 94 of file G4ParticleHPFinalState.hh.

94 { return 0; };
G4double G4ParticleHPFinalState::GetZ ( void  )
inline

Definition at line 97 of file G4ParticleHPFinalState.hh.

97 { return theBaseZ; };

Here is the caller graph for this function:

G4bool G4ParticleHPFinalState::HasAnyData ( )
inline

Definition at line 91 of file G4ParticleHPFinalState.hh.

91 {return hasAnyData;};

Here is the caller graph for this function:

G4bool G4ParticleHPFinalState::HasFSData ( )
inline

Definition at line 90 of file G4ParticleHPFinalState.hh.

90 {return hasFSData;};

Here is the caller graph for this function:

G4bool G4ParticleHPFinalState::HasXsec ( )
inline

Definition at line 89 of file G4ParticleHPFinalState.hh.

89 {return hasXsec;};

Here is the caller graph for this function:

void G4ParticleHPFinalState::Init ( G4double  A,
G4double  Z,
G4String &  dirName,
G4String &  aFSType,
G4ParticleDefinition *  projectile 
)
inline

Definition at line 76 of file G4ParticleHPFinalState.hh.

76 { G4int M = 0; Init ( A, Z, M, dirName, aFSType,const_cast<G4ParticleDefinition*>(projectile)); };
G4ParticleHPFinalState::Init
void Init(G4double A, G4double Z, G4String &dirName, G4String &aFSType, G4ParticleDefinition *projectile)
Definition: G4ParticleHPFinalState.hh:76
G4int
int G4int
Definition: G4Types.hh:78
A
double A(double temperature)
Definition: G4DNAElectronHoleRecombination.cc:60
Z
G4int Z
Definition: G4ParticleHPJENDLHEData.cc:262

Here is the call graph for this function:

Here is the caller graph for this function:

virtual void G4ParticleHPFinalState::Init ( G4double  A,
G4double  Z,
G4int  M,
G4String &  dirName,
G4String &  aFSType,
G4ParticleDefinition *   
)
pure virtual

Implemented in G4FissionLibrary, G4ParticleHPInelasticCompFS, G4ParticleHPInelasticBaseFS, G4ParticleHPElasticFS, G4ParticleHPFissionBaseFS, G4ParticleHPFSFissionFS, G4ParticleHPCaptureFS, G4ParticleHPFissionFS, G4ParticleHP2AInelasticFS, G4ParticleHP3NAInelasticFS, G4ParticleHP3NInelasticFS, G4ParticleHP3NPInelasticFS, G4ParticleHP4NInelasticFS, G4ParticleHPAInelasticFS, G4ParticleHP2N2AInelasticFS, G4ParticleHP2NAInelasticFS, G4ParticleHP2NDInelasticFS, G4ParticleHP2NInelasticFS, G4ParticleHP2NPInelasticFS, G4ParticleHP2PInelasticFS, G4ParticleHP3AInelasticFS, G4ParticleHPDAInelasticFS, G4ParticleHPDInelasticFS, G4ParticleHPHe3InelasticFS, G4ParticleHPN2AInelasticFS, G4ParticleHPN2PInelasticFS, G4ParticleHPN3AInelasticFS, G4ParticleHPNAInelasticFS, G4ParticleHPND2AInelasticFS, G4ParticleHPNDInelasticFS, G4ParticleHPNHe3InelasticFS, G4ParticleHPNInelasticFS, G4ParticleHPNPAInelasticFS, G4ParticleHPNPInelasticFS, G4ParticleHPNT2AInelasticFS, G4ParticleHPNTInelasticFS, G4ParticleHPNXInelasticFS, G4ParticleHPPAInelasticFS, G4ParticleHPPDInelasticFS, G4ParticleHPPInelasticFS, G4ParticleHPPTInelasticFS, G4ParticleHPTInelasticFS, G4ParticleHPFFFissionFS, G4ParticleHPD2AInelasticFS, G4ParticleHPFCFissionFS, G4ParticleHPLCFissionFS, G4ParticleHPSCFissionFS, G4ParticleHPT2AInelasticFS, and G4ParticleHPTCFissionFS.

virtual G4ParticleHPFinalState* G4ParticleHPFinalState::New ( )
pure virtual

Implemented in G4FissionLibrary, G4ParticleHPInelasticCompFS, G4ParticleHPInelasticBaseFS, G4ParticleHPElasticFS, G4ParticleHPFSFissionFS, G4ParticleHPCaptureFS, G4ParticleHPFissionFS, G4ParticleHP2AInelasticFS, G4ParticleHP3NAInelasticFS, G4ParticleHP3NInelasticFS, G4ParticleHP3NPInelasticFS, G4ParticleHP4NInelasticFS, G4ParticleHPAInelasticFS, G4ParticleHP2N2AInelasticFS, G4ParticleHP2NAInelasticFS, G4ParticleHP2NDInelasticFS, G4ParticleHP2NInelasticFS, G4ParticleHP2NPInelasticFS, G4ParticleHP2PInelasticFS, G4ParticleHP3AInelasticFS, G4ParticleHPDAInelasticFS, G4ParticleHPDInelasticFS, G4ParticleHPHe3InelasticFS, G4ParticleHPN2AInelasticFS, G4ParticleHPN2PInelasticFS, G4ParticleHPN3AInelasticFS, G4ParticleHPNAInelasticFS, G4ParticleHPND2AInelasticFS, G4ParticleHPNDInelasticFS, G4ParticleHPNHe3InelasticFS, G4ParticleHPNInelasticFS, G4ParticleHPNPAInelasticFS, G4ParticleHPNPInelasticFS, G4ParticleHPNT2AInelasticFS, G4ParticleHPNTInelasticFS, G4ParticleHPNXInelasticFS, G4ParticleHPPAInelasticFS, G4ParticleHPPDInelasticFS, G4ParticleHPPInelasticFS, G4ParticleHPPTInelasticFS, G4ParticleHPTInelasticFS, G4ParticleHPFFFissionFS, G4ParticleHPD2AInelasticFS, G4ParticleHPFCFissionFS, G4ParticleHPLCFissionFS, G4ParticleHPSCFissionFS, G4ParticleHPT2AInelasticFS, and G4ParticleHPTCFissionFS.

Here is the caller graph for this function:

void G4ParticleHPFinalState::SetA_Z ( G4double  anA,
G4double  aZ,
G4int  aM = 0 
)
inline

Definition at line 96 of file G4ParticleHPFinalState.hh.

96 {theBaseA = anA; theBaseZ = aZ; theBaseM=aM; };

Here is the caller graph for this function:

void G4ParticleHPFinalState::SetAZMs ( G4double  anA,
G4double  aZ,
G4int  aM,
G4ParticleHPDataUsed  used 
)
inline

Definition at line 102 of file G4ParticleHPFinalState.hh.

103  { theBaseA = anA; theBaseZ = aZ; theBaseM=aM;
104  theNDLDataA=(G4int)used.GetA(); theNDLDataZ=(G4int)used.GetZ(); theNDLDataM=used.GetM(); };
G4int
int G4int
Definition: G4Types.hh:78

Here is the call graph for this function:

Here is the caller graph for this function:

void G4ParticleHPFinalState::SetProjectile ( G4ParticleDefinition *  projectile)
inline

Definition at line 106 of file G4ParticleHPFinalState.hh.

106  {
107  theProjectile = projectile; }
G4ParticleHPFinalState::theProjectile
G4ParticleDefinition * theProjectile
Definition: G4ParticleHPFinalState.hh:136

Here is the caller graph for this function:

Member Data Documentation

G4bool G4ParticleHPFinalState::hasAnyData
protected

Definition at line 113 of file G4ParticleHPFinalState.hh.

G4bool G4ParticleHPFinalState::hasFSData
protected

Definition at line 112 of file G4ParticleHPFinalState.hh.

G4bool G4ParticleHPFinalState::hasXsec
protected

Definition at line 111 of file G4ParticleHPFinalState.hh.

G4double G4ParticleHPFinalState::theBaseA
protected

Definition at line 120 of file G4ParticleHPFinalState.hh.

G4int G4ParticleHPFinalState::theBaseM
protected

Definition at line 122 of file G4ParticleHPFinalState.hh.

G4double G4ParticleHPFinalState::theBaseZ
protected

Definition at line 121 of file G4ParticleHPFinalState.hh.

G4ParticleHPNames G4ParticleHPFinalState::theNames
protected

Definition at line 114 of file G4ParticleHPFinalState.hh.

G4int G4ParticleHPFinalState::theNDLDataA
protected

Definition at line 129 of file G4ParticleHPFinalState.hh.

G4int G4ParticleHPFinalState::theNDLDataM
protected

Definition at line 130 of file G4ParticleHPFinalState.hh.

G4int G4ParticleHPFinalState::theNDLDataZ
protected

Definition at line 127 of file G4ParticleHPFinalState.hh.

G4ParticleDefinition* G4ParticleHPFinalState::theProjectile
protected

Definition at line 136 of file G4ParticleHPFinalState.hh.

G4Cache< G4HadFinalState* > G4ParticleHPFinalState::theResult
protected

Definition at line 117 of file G4ParticleHPFinalState.hh.

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