G4CascadeRecoilMaker Class Reference

#include <G4CascadeRecoilMaker.hh>

Inheritance diagram for G4CascadeRecoilMaker:

Public Member Functions
	G4CascadeRecoilMaker (G4double tolerance=0.001 *CLHEP::MeV)
virtual	~G4CascadeRecoilMaker ()
void	collide (G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &output)
void	collide (G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &output, const std::vector< G4CascadParticle > &cparticles)
void	setTolerance (G4double tolerance)
void	setRecoilExcitation (G4double Eexc)
G4InuclNuclei *	makeRecoilNuclei (G4InuclParticle::Model model=G4InuclParticle::DefaultModel)
G4Fragment *	makeRecoilFragment ()
void	addExcitonConfiguration (const G4ExitonConfiguration exciton)
G4int	getRecoilA () const
G4int	getRecoilZ () const
G4double	getRecoilExcitation () const
const G4LorentzVector &	getRecoilMomentum () const
G4bool	goodFragment () const
G4bool	goodRecoil () const
G4bool	wholeEvent () const
G4bool	unphysicalRecoil () const
G4bool	goodNucleus () const
Public Member Functions inherited from G4VCascadeCollider
	G4VCascadeCollider (const G4String &name, G4int verbose=0)
virtual	~G4VCascadeCollider ()
virtual void	collide (G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &output)=0
virtual void	setVerboseLevel (G4int verbose=0)

Protected Member Functions
void	fillRecoil ()
G4double	deltaM () const
Protected Member Functions inherited from G4VCascadeCollider
virtual void	setName (const G4String &name)

Additional Inherited Members
Protected Attributes inherited from G4VCascadeCollider
G4String	theName
G4int	verboseLevel

Detailed Description

Definition at line 64 of file G4CascadeRecoilMaker.hh.

Constructor & Destructor Documentation

G4CascadeRecoilMaker()

G4CascadeRecoilMaker::G4CascadeRecoilMaker ( G4double  tolerance = 0.001*CLHEP::MeV )

explicit

Definition at line 66 of file G4CascadeRecoilMaker.cc.

  : G4VCascadeCollider("G4CascadeRecoilMaker"),
    excTolerance(tolerance), inputEkin(0.),
    recoilA(0), recoilZ(0), excitationEnergy(0.) {
  balance = new G4CascadeCheckBalance(tolerance, tolerance, theName);
}

~G4CascadeRecoilMaker()

G4CascadeRecoilMaker::~G4CascadeRecoilMaker ( )

virtual

Definition at line 73 of file G4CascadeRecoilMaker.cc.

                                            {
  delete balance;
}

Member Function Documentation

addExcitonConfiguration()

void G4CascadeRecoilMaker::addExcitonConfiguration ( const G4ExitonConfiguration  exciton )

inline

Definition at line 88 of file G4CascadeRecoilMaker.hh.

                                                                    {
    theExcitons = exciton;
  }

Referenced by G4IntraNucleiCascader::finishCascade().

collide() [1/2]

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

virtual

Implements G4VCascadeCollider.

Definition at line 80 of file G4CascadeRecoilMaker.cc.

                                               {
  if (verboseLevel > 1)
    G4cout << " >>> G4CascadeRecoilMaker::collide" << G4endl;
 
  // Available energy needed for "goodNucleus()" test at end
  inputEkin = bullet ? bullet->getKineticEnergy() : 0.;
 
  balance->setVerboseLevel(verboseLevel);
  balance->collide(bullet, target, output);
  fillRecoil();
}

Referenced by G4IntraNucleiCascader::finishCascade(), and G4IntraNucleiCascader::generateCascade().

collide() [2/2]

void G4CascadeRecoilMaker::collide	(	G4InuclParticle *	bullet,
		G4InuclParticle *	target,
		G4CollisionOutput &	output,
		const std::vector< G4CascadParticle > &	cparticles
	)

Definition at line 96 of file G4CascadeRecoilMaker.cc.

                                                        {
  if (verboseLevel > 1)
    G4cout << " >>> G4CascadeRecoilMaker::collide(<EP>,<CP>)" << G4endl;
 
  // Available energy needed for "goodNucleus()" test at end
  inputEkin = bullet ? bullet->getKineticEnergy() : 0.;
 
  balance->setVerboseLevel(verboseLevel);
  balance->collide(bullet, target, output, cparticles);
  fillRecoil();
}

deltaM()

G4double G4CascadeRecoilMaker::deltaM ( ) const

protected

Definition at line 200 of file G4CascadeRecoilMaker.cc.

                                            {
  G4double nucMass = G4InuclNuclei::getNucleiMass(recoilA,recoilZ);
  return (recoilMomentum.m() - nucMass);
}

Referenced by fillRecoil().

fillRecoil()

void G4CascadeRecoilMaker::fillRecoil ( )

protected

Definition at line 115 of file G4CascadeRecoilMaker.cc.

                                      {
  recoilZ = -(balance->deltaQ());       // Charge "non-conservation"
  recoilA = -(balance->deltaB());       // Baryon "non-conservation"
  recoilMomentum = -(balance->deltaLV());
 
  theExcitons.clear();      // Discard previous exciton configuraiton
 
  // Bertini uses MeV for excitation energy
  if (!goodFragment()) excitationEnergy = 0.;
  else excitationEnergy = deltaM() * GeV;
 
  // Allow for very small negative mass difference, and round to zero
  if (std::abs(excitationEnergy) < excTolerance) excitationEnergy = 0.;
 
  if (verboseLevel > 2) {
    G4cout << "  recoil px " << recoilMomentum.px()
       << " py " << recoilMomentum.py() << " pz " << recoilMomentum.pz()
       << " E " << recoilMomentum.e() << " baryon " << recoilA
       << " charge " << recoilZ
       << "\n  recoil mass " << recoilMomentum.m()
       << " 'excitation' energy " << excitationEnergy << G4endl;
  }
}

Referenced by collide().

getRecoilA()

G4int G4CascadeRecoilMaker::getRecoilA ( ) const

inline

Definition at line 93 of file G4CascadeRecoilMaker.hh.

{ return recoilA; }

Referenced by G4IntraNucleiCascader::finishCascade(), and G4IntraNucleiCascader::generateCascade().

getRecoilExcitation()

G4double G4CascadeRecoilMaker::getRecoilExcitation ( ) const

inline

Definition at line 95 of file G4CascadeRecoilMaker.hh.

{ return excitationEnergy; }

Referenced by G4IntraNucleiCascader::finishCascade().

getRecoilMomentum()

const G4LorentzVector & G4CascadeRecoilMaker::getRecoilMomentum ( ) const

inline

Definition at line 96 of file G4CascadeRecoilMaker.hh.

{ return recoilMomentum; }

Referenced by G4IntraNucleiCascader::finishCascade().

getRecoilZ()

G4int G4CascadeRecoilMaker::getRecoilZ ( ) const

inline

Definition at line 94 of file G4CascadeRecoilMaker.hh.

{ return recoilZ; }

Referenced by G4IntraNucleiCascader::finishCascade().

goodFragment()

G4bool G4CascadeRecoilMaker::goodFragment

(

)

const

Definition at line 208 of file G4CascadeRecoilMaker.cc.

                                                {
  return (recoilA>0 && recoilZ>=0 && recoilA >= recoilZ);
}

Referenced by fillRecoil(), G4IntraNucleiCascader::finishCascade(), goodNucleus(), and goodRecoil().

goodNucleus()

G4bool G4CascadeRecoilMaker::goodNucleus

(

)

const

Definition at line 234 of file G4CascadeRecoilMaker.cc.

                                               {
  if (verboseLevel > 2) {
    G4cout << " >>> G4CascadeRecoilMaker::goodNucleus" << G4endl;
  }
 
  const G4double minExcitation = 0.1*keV;
  const G4double reasonableExcitation = 7.0;    // Multiple of binding energy
  const G4double fractionalExcitation = 0.2;    // Fraction of input to excite
 
  if (!goodRecoil()) {
    if (verboseLevel>2) {
      if (!goodFragment()) G4cerr << " goodNucleus: invalid A/Z" << G4endl;
      else if (excitationEnergy < -excTolerance) 
    G4cerr << " goodNucleus: negative excitation" << G4endl;
    }
    return false;               // Not a sensible nucleus
  }
 
  if (excitationEnergy <= minExcitation) return true;   // Effectively zero
 
  // Maximum possible excitation energy determined by initial energy
  G4double dm = bindingEnergy(recoilA,recoilZ);
  G4double exc_max0z = fractionalExcitation * inputEkin*GeV;
  G4double exc_dm    = reasonableExcitation * dm;
  G4double exc_max = (exc_max0z > exc_dm) ? exc_max0z : exc_dm;
  
  if (verboseLevel > 3) {
    G4cout << " eexs " << excitationEnergy << " max " << exc_max
       << " dm " << dm << G4endl;
  }
 
  if (verboseLevel > 2 && excitationEnergy >= exc_max)
    G4cerr << " goodNucleus: too much excitation" << G4endl;
 
  return (excitationEnergy < exc_max);      // Below maximum possible
}

Referenced by G4IntraNucleiCascader::finishCascade().

goodRecoil()

G4bool G4CascadeRecoilMaker::goodRecoil

(

)

const

Definition at line 212 of file G4CascadeRecoilMaker.cc.

                                              {
  return (goodFragment() && excitationEnergy > -excTolerance);
}

Referenced by goodNucleus(), makeRecoilFragment(), makeRecoilNuclei(), and unphysicalRecoil().

makeRecoilFragment()

G4Fragment * G4CascadeRecoilMaker::makeRecoilFragment

(

)

Definition at line 164 of file G4CascadeRecoilMaker.cc.

                                                     {
  if (verboseLevel > 1) 
    G4cout << " >>> G4CascadeRecoilMaker::makeRecoilFragment" << G4endl;
 
  if (!goodRecoil()) {
    if (verboseLevel > 2 && !wholeEvent())
      G4cout << theName << ": event recoil is not a physical nucleus" << G4endl;
 
    return 0;       // Null pointer means no fragment
  }
 
  theRecoilFragment.SetZandA_asInt(recoilZ, recoilA);   // Note convention!
 
  // User may have overridden excitation energy; force four-momentum to match
  G4double fragMass = 
    G4InuclNuclei::getNucleiMass(recoilA,recoilZ) + excitationEnergy/GeV;
 
  G4LorentzVector fragMom; fragMom.setVectM(recoilMomentum.vect(), fragMass);
  theRecoilFragment.SetMomentum(fragMom*GeV);   // Bertini uses GeV!
 
  // Note:  exciton configuration has to be set piece by piece
  //        (arguments are Ntotal,Nproton in both cases)
  G4int nholes = theExcitons.protonHoles+theExcitons.neutronHoles;
  theRecoilFragment.SetNumberOfHoles(nholes, theExcitons.protonHoles);
 
  G4int nexcit = (theExcitons.protonQuasiParticles
          + theExcitons.neutronQuasiParticles);
  theRecoilFragment.SetNumberOfExcitedParticle(nexcit,
                    theExcitons.protonQuasiParticles);
 
  return &theRecoilFragment;
}

Referenced by G4IntraNucleiCascader::finishCascade().

makeRecoilNuclei()

G4InuclNuclei * G4CascadeRecoilMaker::makeRecoilNuclei

(

G4InuclParticle::Model

model = G4InuclParticle::DefaultModel

)

Definition at line 143 of file G4CascadeRecoilMaker.cc.

                                                                 {
  if (verboseLevel > 1) 
    G4cout << " >>> G4CascadeRecoilMaker::makeRecoilNuclei" << G4endl;
 
  if (!goodRecoil()) {
    if (verboseLevel > 2 && !wholeEvent())
      G4cout << theName << ": event recoil is not a physical nucleus" << G4endl;
 
    return 0;       // Null pointer means no fragment
  }
 
  theRecoilNuclei.fill(recoilMomentum, recoilA, recoilZ,
               excitationEnergy, model);
  theRecoilNuclei.setExitonConfiguration(theExcitons);
 
  return &theRecoilNuclei;
}

setRecoilExcitation()

void G4CascadeRecoilMaker::setRecoilExcitation ( G4double  Eexc )

inline

Definition at line 81 of file G4CascadeRecoilMaker.hh.

{ excitationEnergy = Eexc; }

Referenced by G4IntraNucleiCascader::finishCascade().

setTolerance()

void G4CascadeRecoilMaker::setTolerance ( G4double  tolerance )

inline

Definition at line 79 of file G4CascadeRecoilMaker.hh.

{ excTolerance = tolerance; }

Referenced by G4IntraNucleiCascader::initialize().

unphysicalRecoil()

G4bool G4CascadeRecoilMaker::unphysicalRecoil ( ) const

inline

Definition at line 102 of file G4CascadeRecoilMaker.hh.

{ return !wholeEvent() && !goodRecoil(); }

wholeEvent()

G4bool G4CascadeRecoilMaker::wholeEvent

(

)

const

Definition at line 216 of file G4CascadeRecoilMaker.cc.

                                              {
  if (verboseLevel > 2) {
    G4cout << " >>> G4CascadeRecoilMaker::wholeEvent:"
       << " A " << recoilA << " Z " << recoilZ
       << " P " << recoilMomentum.rho() << " E " << recoilMomentum.e()
       << "\n wholeEvent returns "
       << (recoilA==0 && recoilZ==0 && 
           recoilMomentum.rho() < excTolerance/GeV &&
           std::abs(recoilMomentum.e()) < excTolerance/GeV) << G4endl;
  }
 
  return (recoilA==0 && recoilZ==0 && 
      recoilMomentum.rho() < excTolerance/GeV &&
      std::abs(recoilMomentum.e()) < excTolerance/GeV);
}

Referenced by G4IntraNucleiCascader::finishCascade(), makeRecoilFragment(), makeRecoilNuclei(), and unphysicalRecoil().

---

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

• G4CascadeRecoilMaker.hh
• G4CascadeRecoilMaker.cc