00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033 #define INCLXX_IN_GEANT4_MODE 1
00034
00035 #include "globals.hh"
00036
00037
00038
00039
00040
00041
00042
00043
00044 #ifndef G4INCLStandardPropagationModel_hh
00045 #define G4INCLStandardPropagationModel_hh 1
00046
00047 #include "G4INCLNucleus.hh"
00048 #include "G4INCLIPropagationModel.hh"
00049 #include "G4INCLIAvatar.hh"
00050 #include "G4INCLConfigEnums.hh"
00051
00052 #include <iterator>
00053
00054 namespace G4INCL {
00055
00068 class StandardPropagationModel: public G4INCL::IPropagationModel {
00069 public:
00070 StandardPropagationModel(LocalEnergyType localEnergyType, LocalEnergyType localEnergyDeltaType);
00071 virtual ~StandardPropagationModel();
00072
00073 G4double getCurrentTime();
00077 void setNucleus(G4INCL::Nucleus *nucleus);
00078
00082 G4INCL::Nucleus* getNucleus();
00083
00084 G4double shoot(ParticleSpecies const projectileSpecies, const G4double kineticEnergy, const G4double impactParameter, const G4double phi);
00085 G4double shootParticle(ParticleType const t, const G4double kineticEnergy, const G4double impactParameter, const G4double phi);
00086 G4double shootComposite(ParticleSpecies const s, const G4double kineticEnergy, const G4double impactParameter, const G4double phi);
00087
00091 void setStoppingTime(G4double);
00092
00096 G4double getStoppingTime();
00097
00101 void registerAvatar(G4INCL::IAvatar *anAvatar);
00102
00108 IAvatar *generateBinaryCollisionAvatar(Particle * const p1, Particle * const p2) const;
00109
00116 G4double getReflectionTime(G4INCL::Particle const * const aParticle);
00117
00121 G4double getTime(G4INCL::Particle const * const particleA,
00122 G4INCL::Particle const * const particleB, G4double *minDistOfApproach) const;
00123
00137 void generateUpdatedCollisions(const ParticleList &updatedParticles, const ParticleList &particles);
00138
00151 void generateCollisions(const ParticleList &particles, const ParticleList &except);
00152
00160 void generateDecays(const ParticleList &particles);
00161
00165 void updateAvatars(const ParticleList &particles);
00166
00171 void generateAllAvatars(G4bool excludeUpdated=false);
00172
00176 G4INCL::IAvatar* propagate();
00177
00178 private:
00179 G4INCL::Nucleus *theNucleus;
00180 G4double maximumTime;
00181 G4double currentTime;
00182 G4bool firstAvatar;
00183 LocalEnergyType theLocalEnergyType, theLocalEnergyDeltaType;
00184
00186 void putSpectatorsOnShell(IAvatarList const &entryAvatars, ParticleList const &spectators);
00187 };
00188
00189 }
00190
00191
00192 #endif