G4CascadeRecoilMaker Class Reference

#include <G4CascadeRecoilMaker.hh>

Inheritance diagram for G4CascadeRecoilMaker:

Collaboration 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	setVerboseLevel (G4int verbose=0)

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

Private Member Functions
	G4CascadeRecoilMaker (const G4CascadeRecoilMaker &)
G4CascadeRecoilMaker &	operator= (const G4CascadeRecoilMaker &)

Private Attributes
G4CascadeCheckBalance *	balance
G4double	excTolerance
G4double	inputEkin
G4int	recoilA
G4int	recoilZ
G4LorentzVector	recoilMomentum
G4double	excitationEnergy
G4ExitonConfiguration	theExcitons
G4InuclNuclei	theRecoilNuclei
G4Fragment	theRecoilFragment

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

Detailed Description

Definition at line 65 of file G4CascadeRecoilMaker.hh.

Constructor & Destructor Documentation

G4CascadeRecoilMaker() [1/2]

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

explicit

Definition at line 67 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 74 of file G4CascadeRecoilMaker.cc.

                                            {
  delete balance;
}

G4CascadeRecoilMaker() [2/2]

G4CascadeRecoilMaker::G4CascadeRecoilMaker ( const G4CascadeRecoilMaker &  )

private

Member Function Documentation

addExcitonConfiguration()

void G4CascadeRecoilMaker::addExcitonConfiguration ( const G4ExitonConfiguration  exciton )

inline

Definition at line 89 of file G4CascadeRecoilMaker.hh.

                                                                    {
    theExcitons = exciton;
  }

Here is the caller graph for this function:

collide() [1/2]

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

virtual

Implements G4VCascadeCollider.

Definition at line 81 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();
}

Here is the call graph for this function:

Here is the caller graph for this function:

collide() [2/2]

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

Definition at line 97 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();
}

Here is the call graph for this function:

deltaM()

G4double G4CascadeRecoilMaker::deltaM ( ) const

protected

Definition at line 201 of file G4CascadeRecoilMaker.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

fillRecoil()

void G4CascadeRecoilMaker::fillRecoil ( )

protected

Definition at line 116 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;
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

getRecoilA()

G4int G4CascadeRecoilMaker::getRecoilA ( ) const

inline

Definition at line 94 of file G4CascadeRecoilMaker.hh.

{ return recoilA; }

Here is the caller graph for this function:

getRecoilExcitation()

G4double G4CascadeRecoilMaker::getRecoilExcitation ( ) const

inline

Definition at line 96 of file G4CascadeRecoilMaker.hh.

{ return excitationEnergy; }

Here is the caller graph for this function:

getRecoilMomentum()

const G4LorentzVector& G4CascadeRecoilMaker::getRecoilMomentum ( ) const

inline

Definition at line 97 of file G4CascadeRecoilMaker.hh.

{ return recoilMomentum; }

Here is the call graph for this function:

Here is the caller graph for this function:

getRecoilZ()

G4int G4CascadeRecoilMaker::getRecoilZ ( ) const

inline

Definition at line 95 of file G4CascadeRecoilMaker.hh.

{ return recoilZ; }

Here is the caller graph for this function:

goodFragment()

G4bool G4CascadeRecoilMaker::goodFragment

(

)

const

Definition at line 209 of file G4CascadeRecoilMaker.cc.

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

Here is the caller graph for this function:

goodNucleus()

G4bool G4CascadeRecoilMaker::goodNucleus

(

)

const

Definition at line 235 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
}

Here is the call graph for this function:

Here is the caller graph for this function:

goodRecoil()

G4bool G4CascadeRecoilMaker::goodRecoil

(

)

const

Definition at line 213 of file G4CascadeRecoilMaker.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

makeRecoilFragment()

G4Fragment * G4CascadeRecoilMaker::makeRecoilFragment

(

)

Definition at line 165 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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

makeRecoilNuclei()

G4InuclNuclei * G4CascadeRecoilMaker::makeRecoilNuclei

(

G4InuclParticle::Model

model = G4InuclParticle::DefaultModel

)

Definition at line 144 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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

operator=()

G4CascadeRecoilMaker& G4CascadeRecoilMaker::operator= ( const G4CascadeRecoilMaker &  )

private

setRecoilExcitation()

void G4CascadeRecoilMaker::setRecoilExcitation ( G4double  Eexc )

inline

Definition at line 82 of file G4CascadeRecoilMaker.hh.

{ excitationEnergy = Eexc; }

Here is the call graph for this function:

Here is the caller graph for this function:

setTolerance()

void G4CascadeRecoilMaker::setTolerance ( G4double  tolerance )

inline

Definition at line 80 of file G4CascadeRecoilMaker.hh.

{ excTolerance = tolerance; }

Here is the caller graph for this function:

unphysicalRecoil()

G4bool G4CascadeRecoilMaker::unphysicalRecoil ( ) const

inline

Definition at line 103 of file G4CascadeRecoilMaker.hh.

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

Here is the call graph for this function:

wholeEvent()

G4bool G4CascadeRecoilMaker::wholeEvent

(

)

const

Definition at line 217 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);
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

balance

G4CascadeCheckBalance* G4CascadeRecoilMaker::balance

private

Definition at line 112 of file G4CascadeRecoilMaker.hh.

excitationEnergy

G4double G4CascadeRecoilMaker::excitationEnergy

private

Definition at line 121 of file G4CascadeRecoilMaker.hh.

excTolerance

G4double G4CascadeRecoilMaker::excTolerance

private

Definition at line 114 of file G4CascadeRecoilMaker.hh.

inputEkin

G4double G4CascadeRecoilMaker::inputEkin

private

Definition at line 116 of file G4CascadeRecoilMaker.hh.

recoilA

G4int G4CascadeRecoilMaker::recoilA

private

Definition at line 118 of file G4CascadeRecoilMaker.hh.

recoilMomentum

G4LorentzVector G4CascadeRecoilMaker::recoilMomentum

private

Definition at line 120 of file G4CascadeRecoilMaker.hh.

recoilZ

G4int G4CascadeRecoilMaker::recoilZ

private

Definition at line 119 of file G4CascadeRecoilMaker.hh.

theExcitons

G4ExitonConfiguration G4CascadeRecoilMaker::theExcitons

private

Definition at line 123 of file G4CascadeRecoilMaker.hh.

theRecoilFragment

G4Fragment G4CascadeRecoilMaker::theRecoilFragment

private

Definition at line 126 of file G4CascadeRecoilMaker.hh.

theRecoilNuclei

G4InuclNuclei G4CascadeRecoilMaker::theRecoilNuclei

private

Definition at line 125 of file G4CascadeRecoilMaker.hh.

---

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

• G4CascadeRecoilMaker.hh
• G4CascadeRecoilMaker.cc