G4BigBanger Class Reference

#include <G4BigBanger.hh>

Inheritance diagram for G4BigBanger:

Public Member Functions
	G4BigBanger ()
virtual	~G4BigBanger ()
void	collide (G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &output)
Public Member Functions inherited from G4CascadeColliderBase
	G4CascadeColliderBase (const char *name, G4int verbose=0)
virtual	~G4CascadeColliderBase ()
virtual void	rescatter (G4InuclParticle *, G4KineticTrackVector *, G4V3DNucleus *, G4CollisionOutput &)
virtual void	setVerboseLevel (G4int verbose=0)
virtual void	setConservationChecks (G4bool doBalance=true)
Public Member Functions inherited from G4VCascadeCollider
	G4VCascadeCollider (const char *name, G4int verbose=0)
virtual	~G4VCascadeCollider ()
virtual void	collide (G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &output)=0
virtual void	setVerboseLevel (G4int verbose=0)

Additional Inherited Members
Protected Member Functions inherited from G4CascadeColliderBase
virtual G4bool	useEPCollider (G4InuclParticle *bullet, G4InuclParticle *target) const
virtual G4bool	explosion (G4InuclNuclei *target) const
virtual G4bool	explosion (G4Fragment *target) const
virtual G4bool	explosion (G4int A, G4int Z, G4double excitation) const
virtual G4bool	inelasticInteractionPossible (G4InuclParticle *bullet, G4InuclParticle *target, G4double ekin) const
virtual G4bool	validateOutput (G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &output)
virtual G4bool	validateOutput (G4InuclParticle *bullet, G4InuclParticle *target, const std::vector< G4InuclElementaryParticle > &particles)
virtual G4bool	validateOutput (G4InuclParticle *bullet, G4InuclParticle *target, const std::vector< G4InuclNuclei > &fragments)
Protected Member Functions inherited from G4VCascadeCollider
virtual void	setName (const char *name)
Protected Attributes inherited from G4CascadeColliderBase
G4InteractionCase	interCase
G4bool	doConservationChecks
G4CascadeCheckBalance *	balance
Protected Attributes inherited from G4VCascadeCollider
const char *	theName
G4int	verboseLevel

Detailed Description

Definition at line 47 of file G4BigBanger.hh.

Constructor & Destructor Documentation

G4BigBanger()

G4BigBanger::G4BigBanger

(

)

Definition at line 62 of file G4BigBanger.cc.

: G4CascadeColliderBase("G4BigBanger") {}

~G4BigBanger()

virtual G4BigBanger::~G4BigBanger ( )

inlinevirtual

Definition at line 50 of file G4BigBanger.hh.

{};

Member Function Documentation

collide()

void G4BigBanger::collide ( G4InuclParticle *  bullet, G4InuclParticle *  target, G4CollisionOutput &  output  )

virtual

Implements G4VCascadeCollider.

Definition at line 65 of file G4BigBanger.cc.

                                        {
 
  if (verboseLevel) G4cout << " >>> G4BigBanger::collide" << G4endl;
 
  // primitive explosion model A -> nucleons to prevent too exotic evaporation
 
  G4InuclNuclei* nuclei_target = dynamic_cast<G4InuclNuclei*>(target);
  if (!nuclei_target) {
    G4cerr << " BigBanger -> try to bang not nuclei " << G4endl;
    return;
  }
 
  G4int A = nuclei_target->getA();
  G4int Z = nuclei_target->getZ();
 
  G4LorentzVector PEX = nuclei_target->getMomentum();
  G4double EEXS = nuclei_target->getExitationEnergy();
 
  G4ThreeVector toTheLabFrame = PEX.boostVector();  // From rest to lab
 
  // This "should" be difference between E-target and sum of m(nucleons)
  G4double etot = (EEXS - bindingEnergy(A,Z)) * MeV/GeV;  // To Bertini units
  if (etot < 0.0) etot = 0.0;
  
  if (verboseLevel > 2) {
    G4cout << " BigBanger: target\n" << *nuclei_target
       << "\n etot " << etot << G4endl;
  }
 
  if (verboseLevel > 3) {
    G4LorentzVector PEXrest = PEX;
    PEXrest.boost(-toTheLabFrame);
    G4cout << " target rest frame: px " << PEXrest.px() << " py "
       << PEXrest.py() << " pz " << PEXrest.pz() << " E " << PEXrest.e()
       << G4endl;
  }
 
  generateBangInSCM(etot, A, Z);
  
  if (verboseLevel > 2) {
    G4cout << " particles " << particles.size() << G4endl;
    for(G4int i = 0; i < G4int(particles.size()); i++) 
      G4cout << particles[i] << G4endl;
  }
 
  if (particles.empty()) {  // No bang!  Don't know why...
    G4cerr << " >>> G4BigBanger unable to process fragment "
       << nuclei_target->getDefinition()->GetParticleName() << G4endl;
 
    // FIXME:  This will violate baryon number, momentum, energy, etc.
    return;
  }
 
  // convert back to Lab
  G4LorentzVector totscm;
  G4LorentzVector totlab;
 
  if (verboseLevel > 2) G4cout << " BigBanger: boosting to lab" << G4endl;
 
  particleIterator ipart;
  for(ipart = particles.begin(); ipart != particles.end(); ipart++) {
    G4LorentzVector mom = ipart->getMomentum();
    if (verboseLevel > 2) totscm += mom;
 
    mom.boost(toTheLabFrame);
    if (verboseLevel > 2) totlab += mom;
 
    ipart->setMomentum(mom); 
    if (verboseLevel > 2) G4cout << *ipart << G4endl;
  }
  
  std::sort(particles.begin(), particles.end(), G4ParticleLargerEkin());
 
  validateOutput(0, target, particles);     // Checks <vector> directly
  
  if (verboseLevel > 2) {
    G4cout << " In SCM: total outgoing momentum " << G4endl 
       << " E " << totscm.e() << " px " << totscm.x()
       << " py " << totscm.y() << " pz " << totscm.z() << G4endl; 
    G4cout << " In Lab: mom cons " << G4endl 
       << " E " << PEX.e() - totlab.e() // PEX now includes EEXS
       << " px " << PEX.x() - totlab.x()
       << " py " << PEX.y() - totlab.y() 
       << " pz " << PEX.z() - totlab.z() << G4endl; 
  }
 
  output.addOutgoingParticles(particles);
}            

Referenced by G4CascadeDeexcitation::collide(), and G4EquilibriumEvaporator::collide().

---

The documentation for this class was generated from the following files:

• G4BigBanger.hh
• G4BigBanger.cc