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 char *name, G4int verbose=0)
virtual	~G4VCascadeCollider ()
virtual void	setVerboseLevel (G4int verbose=0)

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

Additional Inherited Members
Protected Attributes inherited from G4VCascadeCollider
const char *	theName
G4int	verboseLevel

Detailed Description

Definition at line 65 of file G4CascadeRecoilMaker.hh.

Constructor & Destructor Documentation

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

explicit

Definition at line 67 of file G4CascadeRecoilMaker.cc.

References G4VCascadeCollider::theName.

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

G4CascadeRecoilMaker::~G4CascadeRecoilMaker ( )

virtual

Definition at line 74 of file G4CascadeRecoilMaker.cc.

                                            {
  delete balance;
}

Member Function Documentation

void G4CascadeRecoilMaker::addExcitonConfiguration ( const G4ExitonConfiguration  exciton )

inline

Definition at line 89 of file G4CascadeRecoilMaker.hh.

Referenced by G4IntraNucleiCascader::finishCascade().

                                                                    {
    theExcitons = exciton;
  }

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

virtual

Implements G4VCascadeCollider.

Definition at line 81 of file G4CascadeRecoilMaker.cc.

References G4CascadeCheckBalance::collide(), fillRecoil(), G4cout, G4endl, G4InuclParticle::getKineticEnergy(), G4VCascadeCollider::setVerboseLevel(), and G4VCascadeCollider::verboseLevel.

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

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

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

Definition at line 97 of file G4CascadeRecoilMaker.cc.

References G4CascadeCheckBalance::collide(), fillRecoil(), G4cout, G4endl, G4InuclParticle::getKineticEnergy(), G4VCascadeCollider::setVerboseLevel(), and G4VCascadeCollider::verboseLevel.

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

G4double G4CascadeRecoilMaker::deltaM ( ) const

protected

Definition at line 201 of file G4CascadeRecoilMaker.cc.

References G4InuclNuclei::getNucleiMass(), and CLHEP::HepLorentzVector::m().

Referenced by fillRecoil().

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

void G4CascadeRecoilMaker::fillRecoil ( )

protected

Definition at line 116 of file G4CascadeRecoilMaker.cc.

References G4ExitonConfiguration::clear(), G4CascadeCheckBalance::deltaB(), G4CascadeCheckBalance::deltaLV(), deltaM(), G4CascadeCheckBalance::deltaQ(), CLHEP::HepLorentzVector::e(), G4cout, G4endl, python.hepunit::GeV, goodFragment(), CLHEP::HepLorentzVector::m(), CLHEP::HepLorentzVector::px(), CLHEP::HepLorentzVector::py(), CLHEP::HepLorentzVector::pz(), and G4VCascadeCollider::verboseLevel.

Referenced by collide().

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

G4int G4CascadeRecoilMaker::getRecoilA ( ) const

inline

Definition at line 94 of file G4CascadeRecoilMaker.hh.

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

{ return recoilA; }

G4double G4CascadeRecoilMaker::getRecoilExcitation ( ) const

inline

Definition at line 96 of file G4CascadeRecoilMaker.hh.

Referenced by G4IntraNucleiCascader::finishCascade().

{ return excitationEnergy; }

const G4LorentzVector& G4CascadeRecoilMaker::getRecoilMomentum ( ) const

inline

Definition at line 97 of file G4CascadeRecoilMaker.hh.

Referenced by G4IntraNucleiCascader::finishCascade().

{ return recoilMomentum; }

G4int G4CascadeRecoilMaker::getRecoilZ ( ) const

inline

Definition at line 95 of file G4CascadeRecoilMaker.hh.

Referenced by G4IntraNucleiCascader::finishCascade().

{ return recoilZ; }

G4bool G4CascadeRecoilMaker::goodFragment

(

)

const

Definition at line 209 of file G4CascadeRecoilMaker.cc.

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

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

G4bool G4CascadeRecoilMaker::goodNucleus

(

)

const

Definition at line 235 of file G4CascadeRecoilMaker.cc.

References G4InuclSpecialFunctions::bindingEnergy(), G4cerr, G4cout, G4endl, python.hepunit::GeV, goodFragment(), goodRecoil(), python.hepunit::keV, and G4VCascadeCollider::verboseLevel.

Referenced by G4IntraNucleiCascader::finishCascade().

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

G4bool G4CascadeRecoilMaker::goodRecoil

(

)

const

Definition at line 213 of file G4CascadeRecoilMaker.cc.

References goodFragment().

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

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

G4Fragment * G4CascadeRecoilMaker::makeRecoilFragment

(

)

Definition at line 165 of file G4CascadeRecoilMaker.cc.

References G4cout, G4endl, G4InuclNuclei::getNucleiMass(), python.hepunit::GeV, goodRecoil(), G4ExitonConfiguration::neutronHoles, G4ExitonConfiguration::neutronQuasiParticles, G4ExitonConfiguration::protonHoles, G4ExitonConfiguration::protonQuasiParticles, G4Fragment::SetMomentum(), G4Fragment::SetNumberOfExcitedParticle(), G4Fragment::SetNumberOfHoles(), CLHEP::HepLorentzVector::setVectM(), G4Fragment::SetZandA_asInt(), G4VCascadeCollider::theName, CLHEP::HepLorentzVector::vect(), G4VCascadeCollider::verboseLevel, and wholeEvent().

Referenced by G4IntraNucleiCascader::finishCascade().

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

G4InuclNuclei * G4CascadeRecoilMaker::makeRecoilNuclei

(

G4InuclParticle::Model

model = G4InuclParticle::DefaultModel

)

Definition at line 144 of file G4CascadeRecoilMaker.cc.

References G4InuclNuclei::fill(), G4cout, G4endl, goodRecoil(), G4InuclNuclei::setExitonConfiguration(), G4VCascadeCollider::theName, G4VCascadeCollider::verboseLevel, and wholeEvent().

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

void G4CascadeRecoilMaker::setRecoilExcitation ( G4double  Eexc )

inline

Definition at line 82 of file G4CascadeRecoilMaker.hh.

Referenced by G4IntraNucleiCascader::finishCascade().

{ excitationEnergy = Eexc; }

void G4CascadeRecoilMaker::setTolerance ( G4double  tolerance )

inline

Definition at line 80 of file G4CascadeRecoilMaker.hh.

Referenced by G4IntraNucleiCascader::initialize().

{ excTolerance = tolerance; }

G4bool G4CascadeRecoilMaker::unphysicalRecoil ( ) const

inline

Definition at line 103 of file G4CascadeRecoilMaker.hh.

References goodRecoil(), and wholeEvent().

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

G4bool G4CascadeRecoilMaker::wholeEvent

(

)

const

Definition at line 217 of file G4CascadeRecoilMaker.cc.

References CLHEP::HepLorentzVector::e(), G4cout, G4endl, python.hepunit::GeV, CLHEP::HepLorentzVector::rho(), and G4VCascadeCollider::verboseLevel.

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

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

---

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

• G4CascadeRecoilMaker.hh
• G4CascadeRecoilMaker.cc