Changeset 1340 for trunk/source/processes/hadronic/models/cascade/cascade/src/G4VCascadeCollider.cc

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/src/G4VCascadeCollider.cc (modified) (2 diffs)

trunk/source/processes/hadronic/models/cascade/cascade/src/G4VCascadeCollider.cc

@@ -23 +23 @@
 // * acceptance of all terms of the Geant4 Software license.          *
 // ********************************************************************
+// $Id: G4VCascadeCollider.cc,v 1.4 2010/07/14 15:41:13 mkelsey Exp $
+// Geant4 tag: $Name: hadr-casc-V09-03-85 $
 //
-// $Id: G4VCascadeCollider.cc,v 1.1.2.1 2010/06/25 09:45:08 gunter Exp $
-// Geant4 tag: $Name: geant4-09-04-beta-01 $
+// 20100615  M. Kelsey -- Split constructor to have verbose separately
+// 20100623  M. Kelsey -- Use old bindingEnergy() wrapper (now returns
+//              G4NucleiProperties::GetBindingEnergy()).
+// 20100714  M. Kelsey -- Move concrete functions to G4CascadeColliderBase
 
 #include "G4VCascadeCollider.hh"
-#include "G4HadTmpUtil.hh"
-#include "G4InteractionCase.hh"
-#include "G4InuclElementaryParticle.hh"
-#include "G4InuclNuclei.hh"
-#include "G4InuclSpecialFunctions.hh"
-#include "G4NucleiProperties.hh"
 #include "G4ios.hh"
-
-using namespace G4InuclSpecialFunctions;
 
 
@@ -43 +39 @@
 G4VCascadeCollider::G4VCascadeCollider(const char* name, G4int verbose)
   : theName(name), verboseLevel(verbose) {
-  if (verbose > 3) G4cout << " >>> " << theName << " ctor " << G4endl;
+  if (verboseLevel) G4cout << " >>> " << theName << " ctor " << G4endl;
 }
-
-
-// Both bullet and target must be hadrons or leptons for this to work
-
-G4bool G4VCascadeCollider::useEPCollider(G4InuclParticle* bullet, 
-                                         G4InuclParticle* target) const {
-  return (dynamic_cast<G4InuclElementaryParticle*>(bullet) &&
-          dynamic_cast<G4InuclElementaryParticle*>(target));
-}
-
-
-// Decide wether nuclear fragment is candidate for G4BigBanger
-
-G4bool G4VCascadeCollider::explosion(G4InuclNuclei* target) const {
-  if (verboseLevel > 3) {
-    G4cout << " >>> " << theName << "::explosion" << G4endl;
-  }
-
-  const G4double a_cut = 20.0;
-  const G4double be_cut = 3.0;
-
-  G4double a = target->getA();
-  G4double z = target->getZ();
-  G4double eexs = target->getExitationEnergy();
-
-  // Only small fragments with high excitations can explode
-  G4bool explo = ((a <= a_cut) && 
-                  (eexs >= be_cut * G4NucleiProperties::GetBindingEnergy(G4lrint(a), G4lrint(z)))
-                  );
-
-  return explo;
-}
-
-
-// Decide whether bullet-target interaction is candidate for cascade
-
-G4bool 
-G4VCascadeCollider::inelasticInteractionPossible(G4InuclParticle* bullet,
-                                                 G4InuclParticle* target, 
-                                                 G4double ekin) const {
-  if (verboseLevel > 3) {
-    G4cout << " >>> " << theName << "::inelasticInteractionPossible" << G4endl;
-  }
-
-  // If hadron-hadron collision, defer to ElementaryParticleCollider
-  if (useEPCollider(bullet, target)) return true;
-
-  // See which one of the two (or both) is a nucleus, get properties
-  // FIXME:  Should set a = baryon() for both, but that's not in base
-  G4InuclNuclei* nuclei_bullet = dynamic_cast<G4InuclNuclei*>(bullet);
-  G4double ab = nuclei_bullet ? nuclei_bullet->getA() : 1;      // FIXME
-  G4double zb = nuclei_bullet ? nuclei_bullet->getZ() : bullet->getCharge();
-  
-  G4InuclNuclei* nuclei_target = dynamic_cast<G4InuclNuclei*>(target);
-  G4double at = nuclei_target ? nuclei_target->getA() : 1;      // FIXME
-  G4double zt = nuclei_target ? nuclei_target->getZ() : target->getCharge();
-  
-  // VCOL (Coulomb barrier) used for testing if elastic collision necessary
-  const G4double coeff = 0.001 * 1.2;
-
-  G4double VCOL = coeff * zt * zb / (G4cbrt(at) + G4cbrt(ab)); 
-  
-  G4bool possible = true;       // Force inelastic; should be (ekin >= VCOL)
-
-  if (verboseLevel > 3) {
-    G4cout << " VCOL: " << VCOL << " ekin: " << ekin << " inelastic possible: "
-           << possible << G4endl;
-  }
-
-  return possible;
-}