Changeset 1340 for trunk/source/processes/hadronic/models/cascade/cascade/include/G4CascadeInterface.hh

Timestamp:

Nov 5, 2010, 3:45:55 PM (14 years ago)

Author:

garnier

Message:

update ti head

File:

• trunk/source/processes/hadronic/models/cascade/cascade/include/G4CascadeInterface.hh (modified) (6 diffs)

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

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4CascadeInterface.hh,v 1.16 2010/05/21 18:07:30 mkelsey Exp $
+// $Id: G4CascadeInterface.hh,v 1.23 2010/09/23 18:13:32 mkelsey Exp $
 // Defines an interface to Bertini (BERT) cascade
 // based on INUCL  intra-nuclear transport.models 
@@ -31 +31 @@
 // 20100405  M. Kelsey -- Fix constness of op== and op!=
 // 20100519  M. Kelsey -- Remove Collider data members
+// 20100617  M. Kelsey -- Make G4InuclCollider a local data member
+// 20100723  M. Kelsey -- Move G4CollisionOutput here for reuse
+// 20100916  M. Kelsey -- Add functions to encapsulate ApplyYourself() actions,
+//              make colliders pointers (don't expose dependencies)
+// 20100922  M. Kelsey -- Add functions to select de-excitation method
 
 #ifndef G4CASCADEINTERFACE_H
@@ -38 +43 @@
 #include "G4FragmentVector.hh"
 #include "G4KineticTrackVector.hh"
+#include "G4LorentzRotation.hh"
 #include "G4Nucleon.hh"
 #include "G4Nucleus.hh"
@@ -43 +49 @@
 #include "G4ReactionProduct.hh"
 #include "G4ReactionProductVector.hh"
+
+class G4InuclCollider;
+class G4InuclParticle;
+class G4CollisionOutput;
+class G4CascadeCheckBalance;
+class G4V3DNucleus;
 
 
@@ -58 +70 @@
                                  G4Nucleus& theNucleus); 
 
+  void setVerboseLevel(G4int verbose) { verboseLevel = verbose; }
+
+  // Select betweeen different post-cascade de-excitation models
+  void useCascadeDeexcitation();
+  void usePreCompoundDeexcitation();
+
+protected:
+  // Convert input projectile and target to Bertini internal types
+  void createBullet(const G4HadProjectile& aTrack);
+  void createTarget(G4Nucleus& theNucleus);
+
+  // Evaluate whether any outgoing particles penetrated Coulomb barrier
+  G4bool coulombBarrierViolation() const;
+
+  // Conditions for rejecting cascade attempt
+  G4bool retryInelasticProton() const;
+  G4bool retryInelasticNucleus() const;
+
+  // Fill sparse array with minimum momenta for inelastic on hydrogen
+  void initializeElasticCuts();
+
+  // Transfer Bertini internal final state to hadronics interface
+  void copyOutputToHadronicResult();
+
+  // Terminate job because of energy/momentum/etc. violations
+  void throwNonConservationFailure();
+
 private:
   G4int operator==(const G4CascadeInterface& right) const {
@@ -67 +106 @@
   }
 
+  static const G4int maximumTries;      // Number of iterations for inelastic
+
+  G4double cutElastic[32];              // Bullet momenta for hydrogen target
+
   G4int verboseLevel;
+  G4int numberOfTries;
 
-private:
-  G4HadFinalState theResult;  
+  G4HadFinalState theResult;
+  G4InuclCollider* collider;
+  G4CascadeCheckBalance* balance;
+
+  G4InuclParticle* bullet;
+  G4InuclParticle* target;
+  G4CollisionOutput* output;
+
+  G4LorentzRotation bulletInLabFrame;
 };