Changeset 1315 for trunk/source/processes/hadronic/models/cascade/cascade/include/G4ElementaryParticleCollider.hh

Timestamp:

Jun 18, 2010, 11:42:07 AM (14 years ago)

Author:

garnier

Message:

update geant4-09-04-beta-cand-01 interfaces-V09-03-09 vis-V09-03-08

File:

• trunk/source/processes/hadronic/models/cascade/cascade/include/G4ElementaryParticleCollider.hh (modified) (1 diff)

trunk/source/processes/hadronic/models/cascade/cascade/include/G4ElementaryParticleCollider.hh

@@ -23 +23 @@
 // * acceptance of all terms of the Geant4 Software license.          *
 // ********************************************************************
+// $Id: G4ElementaryParticleCollider.hh,v 1.31 2010/05/21 17:56:34 mkelsey Exp $
+// Geant4 tag: $Name: geant4-09-04-beta-cand-01 $
 //
+// 20100114  M. Kelsey -- Remove G4CascadeMomentum, use G4LorentzVector directly
+// 20100315  M. Kelsey -- Remove "using" directive and unnecessary #includes.
+// 20100407  M. Kelsey -- Eliminate return-by-value std::vector<> by creating
+//              three data buffers for particles, momenta, and particle types.
+//              The following functions now return void and are non-const:
+//                ::generateSCMfinalState()
+//                ::generateMomModules() (also remove input vector)
+//                ::generateStrangeChannelPartTypes()
+//                ::generateSCMpionAbsorption()
+// 20100413  M. Kelsey -- Pass G4CollisionOutput by ref to ::collide(); merge
+//              public vs. private ::collide() functions.
+// 20100511  M. Kelsey -- Remove G4PionSampler and G4NucleonSampler.  Expand
+//              particle-types selector to all modes, not just strangeness.
+// 20100517  M. Kelsey -- Inherit from common base class, make arrays static
+
 #ifndef G4ELEMENTARY_PARTICLE_COLLIDER_HH
 #define G4ELEMENTARY_PARTICLE_COLLIDER_HH
 
-#include "G4Collider.hh"
+#include "G4VCascadeCollider.hh"
 #include "G4InuclElementaryParticle.hh"
-#include "G4InuclSpecialFunctions.hh"
-#include "G4CascadSpecialFunctions.hh"
-#include "G4LorentzConvertor.hh"
-#include "G4NucleonSampler.hh"
-#include "G4PionSampler.hh"
+#include "G4LorentzVector.hh"
+#include <vector>
 
- 
-using namespace G4InuclSpecialFunctions;
-using namespace G4CascadSpecialFunctions;
+class G4LorentzConvertor;
+class G4CollisionOutput;
 
 
-class G4ElementaryParticleCollider {
-
+class G4ElementaryParticleCollider : public G4VCascadeCollider {
 public:
-
   G4ElementaryParticleCollider();
-
-  G4CollisionOutput collide(G4InuclParticle* bullet,
-                            G4InuclParticle* target);
+  virtual ~G4ElementaryParticleCollider() {};
+  
+  void collide(G4InuclParticle* bullet, G4InuclParticle* target,
+               G4CollisionOutput& output);
 
 private:
-
-  G4int verboseLevel;
-
   void initializeArrays();
 
   G4int generateMultiplicity(G4int is, G4double ekin) const;
 
-  void collide(G4InuclElementaryParticle* bullet,
-               G4InuclElementaryParticle* target,
-               G4CollisionOutput& output);
+  void generateOutgoingPartTypes(G4int is, G4int mult, G4double ekin);
 
-      
-  std::vector<G4InuclElementaryParticle> 
-  generateSCMfinalState(G4double ekin, G4double etot_scm, G4double pscm,             
-                        G4InuclElementaryParticle* particle1,
-                        G4InuclElementaryParticle* particle2, 
-                        G4LorentzConvertor* toSCM) const; 
+  void generateSCMfinalState(G4double ekin, G4double etot_scm, G4double pscm,
+                             G4InuclElementaryParticle* particle1,
+                             G4InuclElementaryParticle* particle2, 
+                             G4LorentzConvertor* toSCM); 
+
+  void generateSCMpionAbsorption(G4double etot_scm,
+                                 G4InuclElementaryParticle* particle1,
+                                 G4InuclElementaryParticle* particle2); 
+
+  void generateMomModules(G4int mult, G4int is, G4double ekin,
+                          G4double etot_cm); 
+
+  // Samples the CM momentum for elastic and charge exchange scattering
+  //
+  G4LorentzVector
+  sampleCMmomentumFor2to2(G4int is, G4int kw, G4double ekin,
+                          G4double pscm) const; 
 
 
-  std::vector<G4double>
-  generateMomModules(const std::vector<G4int>& kinds, G4int mult,
-                     G4int is, G4double ekin, G4double etot_cm) const; 
+  // Samples cos(theta) in the CM for elastic and charge exchange scattering
+  //
+  G4double sampleCMcosFor2to2(G4double pscm, G4double pFrac,
+                              G4double pA, G4double pC, G4double pCos) const;
 
 
-  G4CascadeMomentum
-  particleSCMmomentumFor2to2(G4int is, G4int kw, G4double ekin,
-                             G4double pscm) const; 
-
-
-  G4int getElasticCase(G4int is, G4int kw, G4double ekin) const;
-
-
-  std::vector<G4int>
-  generateStrangeChannelPartTypes(G4int is, G4int mult, 
-                                  G4double ekin) const;
-
-
-  G4double
-  getMomModuleFor2toMany(G4int is, G4int mult, G4int knd, 
-                         G4double ekin) const; 
+  G4double getMomModuleFor2toMany(G4int is, G4int mult, G4int knd, 
+                                  G4double ekin) const; 
 
 
   G4bool satisfyTriangle(const std::vector<G4double>& modules) const; 
         
-  G4CascadeMomentum
+  G4LorentzVector
   particleSCMmomentumFor2to3(G4int is, G4int knd, G4double ekin, 
                              G4double pmod) const; 
 
-
-  std::pair<G4double, G4double> 
-  adjustIntervalForElastic(G4double ekin, G4double ak, G4double ae,
-                           G4int k, G4int l, const std::vector<G4double>& ssv, 
-                           G4double st) const;
- 
-  std::vector<G4InuclElementaryParticle> 
-  generateSCMpionAbsorption(G4double etot_scm,
-                            G4InuclElementaryParticle* particle1,
-                            G4InuclElementaryParticle* particle2) const; 
-
-  G4NucleonSampler nucSampler;
-  G4PionSampler piSampler;
+  // Internal buffers for lists of secondaries
+  std::vector<G4InuclElementaryParticle> particles;
+  std::vector<G4double> modules;
+  std::vector<G4int> particle_kinds;
 
   // Parameter arrays
-
-  G4double rmn[14][10][2];    
-  G4double ang[4][4][13];
-  G4double abn[4][4][4];
-
+  static const G4double rmn[14][10][2];    
+  static const G4double abn[4][4][4];
 };
 
-#endif // G4ELEMENTARY_PARTICLE_COLLIDER_HH
+#endif  /* G4ELEMENTARY_PARTICLE_COLLIDER_HH */