G4BigBanger Class Reference

#include <G4BigBanger.hh>

Inheritance diagram for G4BigBanger:

Public Member Functions
	G4BigBanger ()
virtual	~G4BigBanger ()
virtual void	deExcite (const G4Fragment &target, G4CollisionOutput &output)
Public Member Functions inherited from G4CascadeDeexciteBase
	G4CascadeDeexciteBase (const char *name)
virtual	~G4CascadeDeexciteBase ()
virtual void	setVerboseLevel (G4int verbose=0)
Public Member Functions inherited from G4VCascadeDeexcitation
	G4VCascadeDeexcitation (const char *name)
virtual	~G4VCascadeDeexcitation ()
virtual void	collide (G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &globalOutput)
Public Member Functions inherited from G4VCascadeCollider
	G4VCascadeCollider (const char *name, G4int verbose=0)
virtual	~G4VCascadeCollider ()

Additional Inherited Members
Protected Member Functions inherited from G4CascadeDeexciteBase
virtual G4bool	explosion (const G4Fragment &target) const
virtual G4bool	explosion (G4int A, G4int Z, G4double excitation) const
virtual G4bool	validateOutput (const G4Fragment &target, G4CollisionOutput &output)
virtual G4bool	validateOutput (const G4Fragment &target, const std::vector< G4InuclElementaryParticle > &particles)
virtual G4bool	validateOutput (const G4Fragment &target, const std::vector< G4InuclNuclei > &fragments)
void	getTargetData (const G4Fragment &target)
const G4Fragment &	makeFragment (G4LorentzVector mom, G4int A, G4int Z, G4double EX=0.)
const G4Fragment &	makeFragment (G4int A, G4int Z, G4double EX=0.)
Protected Member Functions inherited from G4VCascadeCollider
virtual void	setName (const char *name)
Protected Attributes inherited from G4CascadeDeexciteBase
G4CascadeCheckBalance *	balance
G4int	A
G4int	Z
G4LorentzVector	PEX
G4double	EEXS
G4Fragment	aFragment
Protected Attributes inherited from G4VCascadeCollider
const char *	theName
G4int	verboseLevel

Detailed Description

Definition at line 48 of file G4BigBanger.hh.

Constructor & Destructor Documentation

G4BigBanger::G4BigBanger

(

)

Definition at line 66 of file G4BigBanger.cc.

: G4CascadeDeexciteBase("G4BigBanger") {}

virtual G4BigBanger::~G4BigBanger ( )

inlinevirtual

Definition at line 51 of file G4BigBanger.hh.

{};

Member Function Documentation

void G4BigBanger::deExcite ( const G4Fragment &  target, G4CollisionOutput &  output  )

virtual

Implements G4VCascadeDeexcitation.

Definition at line 68 of file G4BigBanger.cc.

References G4CascadeDeexciteBase::A, G4CollisionOutput::addOutgoingParticles(), G4InuclSpecialFunctions::bindingEnergy(), CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), CLHEP::HepLorentzVector::e(), G4CascadeDeexciteBase::EEXS, G4cerr, G4cout, G4endl, G4CascadeDeexciteBase::getTargetData(), python.hepunit::GeV, python.hepunit::MeV, G4CascadeDeexciteBase::PEX, CLHEP::HepLorentzVector::px(), CLHEP::HepLorentzVector::py(), CLHEP::HepLorentzVector::pz(), sort(), G4CascadeDeexciteBase::validateOutput(), G4VCascadeCollider::verboseLevel, CLHEP::HepLorentzVector::x(), CLHEP::HepLorentzVector::y(), G4CascadeDeexciteBase::Z, and CLHEP::HepLorentzVector::z().

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

                                                {
  if (verboseLevel) G4cout << " >>> G4BigBanger::deExcite" << G4endl;

  getTargetData(target);
  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" << 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 "
       << target << 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(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; 
  }

  globalOutput.addOutgoingParticles(particles);
}            

---

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

• G4BigBanger.hh
• G4BigBanger.cc