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 G4String &name)
virtual	~G4VCascadeDeexcitation ()
virtual void	collide (G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &globalOutput)
Public Member Functions inherited from G4VCascadeCollider
	G4VCascadeCollider (const G4String &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 G4String &name)
Protected Attributes inherited from G4CascadeDeexciteBase
G4CascadeCheckBalance *	balance
G4int	A
G4int	Z
G4LorentzVector	PEX
G4double	EEXS
G4Fragment	aFragment
Protected Attributes inherited from G4VCascadeCollider
G4String	theName
G4int	verboseLevel

Detailed Description

Definition at line 47 of file G4BigBanger.hh.

Constructor & Destructor Documentation

G4BigBanger()

G4BigBanger::G4BigBanger

(

)

Definition at line 66 of file G4BigBanger.cc.

: G4CascadeDeexciteBase("G4BigBanger") {}

~G4BigBanger()

virtual G4BigBanger::~G4BigBanger ( )

inlinevirtual

Definition at line 50 of file G4BigBanger.hh.

{};

Member Function Documentation

deExcite()

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

virtual

Implements G4VCascadeDeexcitation.

Definition at line 68 of file G4BigBanger.cc.

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

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

---

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

• G4BigBanger.hh
• G4BigBanger.cc