Geant4  v4-10.4-release
 모두 클래스 네임스페이스들 파일들 함수 변수 타입정의 열거형 타입 열거형 멤버 Friends 매크로 그룹들 페이지들
G4INCLInteractionAvatar.hh
이 파일의 문서화 페이지로 가기
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // INCL++ intra-nuclear cascade model
27 // Alain Boudard, CEA-Saclay, France
28 // Joseph Cugnon, University of Liege, Belgium
29 // Jean-Christophe David, CEA-Saclay, France
30 // Pekka Kaitaniemi, CEA-Saclay, France, and Helsinki Institute of Physics, Finland
31 // Sylvie Leray, CEA-Saclay, France
32 // Davide Mancusi, CEA-Saclay, France
33 //
34 #define INCLXX_IN_GEANT4_MODE 1
35 
36 #include "globals.hh"
37 
38 /* \file G4INCLInteractionAvatar.hh
39  * \brief Virtual class for interaction avatars.
40  *
41  * This class is inherited by decay and collision avatars. The goal is to
42  * provide a uniform treatment of common physics, such as Pauli blocking,
43  * enforcement of energy conservation, etc.
44  *
45  * \date Mar 1st, 2011
46  * \author Davide Mancusi
47  */
48 
49 #ifndef G4INCLINTERACTIONAVATAR_HH_
50 #define G4INCLINTERACTIONAVATAR_HH_
51 
52 #include "G4INCLIAvatar.hh"
53 #include "G4INCLNucleus.hh"
54 #include "G4INCLFinalState.hh"
55 #include "G4INCLRootFinder.hh"
56 #include "G4INCLKinematicsUtils.hh"
57 #include "G4INCLAllocationPool.hh"
58 
59 namespace G4INCL {
60 
62  public:
65  virtual ~InteractionAvatar();
66 
68  static const G4double locEAccuracy;
70  static const G4int maxIterLocE;
71 
73  static void deleteBackupParticles();
74 
75  protected:
76  virtual G4INCL::IChannel* getChannel() = 0;
77 
79 
84  void preInteractionLocalEnergy(Particle * const p);
85 
93 
94  void preInteraction();
96 
101  void restoreParticles() const;
102 
105 
113 
114  private:
117  public:
122  ViolationEMomentumFunctor(Nucleus * const nucleus, ParticleList const &modAndCre, const G4double totalEnergyBeforeInteraction, ThreeVector const &boost, const G4bool localE);
123  virtual ~ViolationEMomentumFunctor();
124 
130  G4double operator()(const G4double x) const;
131 
133  void cleanUp(const G4bool success) const;
134 
135  private:
139  std::vector<ThreeVector> particleMomenta;
146 
149 
158  void scaleParticleMomenta(const G4double alpha) const;
159 
160  };
161 
164  public:
169  ViolationEEnergyFunctor(Nucleus * const nucleus, Particle * const aParticle, const G4double totalEnergyBeforeInteraction, const G4bool localE);
171 
177  G4double operator()(const G4double x) const;
178 
180  void cleanUp(const G4bool success) const;
181 
186  void setParticleEnergy(const G4double energy) const;
187 
188  private:
206  };
207 
209 
210  protected:
222 
224 
226  };
227 
228 }
229 
230 #endif /* G4INCLINTERACTIONAVATAR_HH_ */
static const G4int maxIterLocE
Max number of iterations for the determination of the local-energy Q-value.
Float_t x
Definition: compare.C:6
const G4bool shouldUseLocalEnergy
True if we should use local energy.
ThreeVector const & boostVector
Pointer to the boost vector.
Static root-finder algorithm.
void preInteractionLocalEnergy(Particle *const p)
Apply local-energy transformation, if appropriate.
#define INCL_DECLARE_ALLOCATION_POOL(T)
static G4ThreadLocal Particle * backupParticle1
G4bool enforceEnergyConservation(FinalState *const fs)
Enforce energy conservation.
G4bool bringParticleInside(Particle *const p)
G4double initialEnergy
Total energy before the interaction.
const char * p
Definition: xmltok.h:285
static void deleteBackupParticles()
Release the memory allocated for the backup particles.
G4bool shouldUseLocalEnergy() const
true if the given avatar should use local energy
G4double energyThreshold
Threshold for the energy of the particle.
#define G4ThreadLocal
Definition: tls.hh:69
InteractionAvatar(G4double, G4INCL::Nucleus *, G4INCL::Particle *)
G4double operator()(const G4double x) const
Compute the energy-conservation violation.
double G4double
Definition: G4Types.hh:76
bool G4bool
Definition: G4Types.hh:79
RootFunctor-derived object for enforcing energy conservation in delta production. ...
ParticleList finalParticles
List of final-state particles.
ViolationEEnergyFunctor(Nucleus *const nucleus, Particle *const aParticle, const G4double totalEnergyBeforeInteraction, const G4bool localE)
Prepare for calling the () operator and setParticleEnergy.
virtual G4INCL::IChannel * getChannel()=0
void setParticleEnergy(const G4double energy) const
Set the energy of the particle.
void preInteractionBlocking()
Store the state of the particles before the interaction.
const G4bool shouldUseLocalEnergy
Whether we should use local energy.
static const G4double alpha
double energy
Definition: plottest35.C:25
Singleton for recycling allocation of instances of a given class.
void cleanUp(const G4bool success) const
Clean up after root finding.
ViolationEMomentumFunctor(Nucleus *const nucleus, ParticleList const &modAndCre, const G4double totalEnergyBeforeInteraction, ThreeVector const &boost, const G4bool localE)
Prepare for calling the () operator and scaleParticleMomenta.
G4double operator()(const G4double x) const
Compute the energy-conservation violation.
int G4int
Definition: G4Types.hh:78
ThreeVector theMomentum
The initial momentum of the particle.
std::vector< ThreeVector > particleMomenta
CM particle momenta, as determined by the channel.
void scaleParticleMomenta(const G4double alpha) const
Scale the momenta of the modified and created particles.
static const G4double locEAccuracy
Target accuracy in the determination of the local-energy Q-value.
void cleanUp(const G4bool success) const
Clean up after root finding.
static G4ThreadLocal Particle * backupParticle2
RootFunctor-derived object for enforcing energy conservation in N-N.
void restoreParticles() const
Restore the state of both particles.
G4double theEnergy
The initial energy of the particle.
G4double initialEnergy
Total energy before the interaction.