69 excTolerance(tolerance), inputEkin(0.),
70 recoilA(0), recoilZ(0), excitationEnergy(0.) {
85 G4cout <<
" >>> G4CascadeRecoilMaker::collide" <<
G4endl;
100 const std::vector<G4CascadParticle>& cparticles) {
102 G4cout <<
" >>> G4CascadeRecoilMaker::collide(<EP>,<CP>)" <<
G4endl;
146 G4cout <<
" >>> G4CascadeRecoilMaker::makeRecoilNuclei" <<
G4endl;
167 G4cout <<
" >>> G4CascadeRecoilMaker::makeRecoilFragment" <<
G4endl;
219 G4cout <<
" >>> G4CascadeRecoilMaker::wholeEvent:"
222 <<
"\n wholeEvent returns "
237 G4cout <<
" >>> G4CascadeRecoilMaker::goodNucleus" <<
G4endl;
241 const G4double reasonableExcitation = 7.0;
242 const G4double fractionalExcitation = 0.2;
248 G4cerr <<
" goodNucleus: negative excitation" <<
G4endl;
258 G4double exc_dm = reasonableExcitation * dm;
259 G4double exc_max = (exc_max0z > exc_dm) ? exc_max0z : exc_dm;
263 <<
" dm " << dm <<
G4endl;
267 G4cerr <<
" goodNucleus: too much excitation" <<
G4endl;
void fill(G4int a, G4int z, G4double exc=0., Model model=DefaultModel)
G4bool goodNucleus() const
G4LorentzVector deltaLV() const
void collide(G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &output)
G4InuclNuclei theRecoilNuclei
void SetNumberOfHoles(G4int valueTot, G4int valueP=0)
G4double excitationEnergy
G4bool wholeEvent() const
static const G4double tolerance
virtual void setVerboseLevel(G4int verbose=0)
G4double getKineticEnergy() const
void collide(G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &output)
void SetNumberOfExcitedParticle(G4int valueTot, G4int valueP)
virtual ~G4CascadeRecoilMaker()
G4bool goodRecoil() const
G4GLOB_DLL std::ostream G4cout
G4bool goodFragment() const
G4Fragment theRecoilFragment
void SetMomentum(const G4LorentzVector &value)
G4double getNucleiMass() const
G4InuclNuclei * makeRecoilNuclei(G4InuclParticle::Model model=G4InuclParticle::DefaultModel)
void setExitonConfiguration(const G4ExitonConfiguration &config)
G4ExitonConfiguration theExcitons
void SetZandA_asInt(G4int Znew, G4int Anew)
G4CascadeCheckBalance * balance
G4LorentzVector recoilMomentum
G4CascadeRecoilMaker(G4double tolerance=0.001 *CLHEP::MeV)
G4double bindingEnergy(G4int A, G4int Z)
G4Fragment * makeRecoilFragment()
G4GLOB_DLL std::ostream G4cerr
G4int neutronQuasiParticles
CLHEP::HepLorentzVector G4LorentzVector
G4int protonQuasiParticles