G4BigBanger Class Reference

#include <G4BigBanger.hh>

Inheritance diagram for G4BigBanger:

Public Member Functions
	G4BigBanger ()
virtual	~G4BigBanger ()
void	collide (G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &output)

---

Detailed Description

Definition at line 47 of file G4BigBanger.hh.

---

Constructor & Destructor Documentation

G4BigBanger::G4BigBanger

(

)

Definition at line 62 of file G4BigBanger.cc.

: G4CascadeColliderBase("G4BigBanger") {}

virtual G4BigBanger::~G4BigBanger

(

)

[inline, virtual]

Definition at line 50 of file G4BigBanger.hh.

{};

---

Member Function Documentation

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

Implements G4VCascadeCollider.

Definition at line 65 of file G4BigBanger.cc.

References G4CollisionOutput::addOutgoingParticles(), G4InuclSpecialFunctions::bindingEnergy(), G4cerr, G4cout, G4endl, G4InuclNuclei::getA(), G4InuclParticle::getDefinition(), G4InuclNuclei::getExitationEnergy(), G4InuclParticle::getMomentum(), G4ParticleDefinition::GetParticleName(), G4InuclNuclei::getZ(), G4CascadeColliderBase::validateOutput(), and G4VCascadeCollider::verboseLevel.

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

                                                {

  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);
}                    

---

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

• G4BigBanger.hh
• G4BigBanger.cc

---