G4InuclEvaporation Class Reference

#include <G4InuclEvaporation.hh>

Inheritance diagram for G4InuclEvaporation:

Public Member Functions
	G4InuclEvaporation ()
	~G4InuclEvaporation ()
G4FragmentVector *	BreakItUp (const G4Fragment &theNucleus)
void	setVerboseLevel (const G4int verbose)
Public Member Functions inherited from G4VEvaporation
	G4VEvaporation ()
virtual	~G4VEvaporation ()
virtual void	Initialise ()
virtual void	SetPhotonEvaporation (G4VEvaporationChannel *ptr)
G4VEvaporationChannel *	GetPhotonEvaporation ()
void	SetOPTxs (G4int opt)
void	UseSICB (G4bool use)

Additional Inherited Members
Protected Attributes inherited from G4VEvaporation
G4VEvaporationChannel *	thePhotonEvaporation
G4int	OPTxs
G4bool	useSICB

Detailed Description

Definition at line 44 of file G4InuclEvaporation.hh.

Constructor & Destructor Documentation

G4InuclEvaporation::G4InuclEvaporation

(

)

Definition at line 72 of file G4InuclEvaporation.cc.

  : verboseLevel(0),   evaporator(new G4EvaporationInuclCollider) {}

G4InuclEvaporation::~G4InuclEvaporation

(

)

Definition at line 79 of file G4InuclEvaporation.cc.

                                        {
  delete evaporator;
}

Member Function Documentation

G4FragmentVector * G4InuclEvaporation::BreakItUp ( const G4Fragment &  theNucleus )

virtual

Implements G4VEvaporation.

Definition at line 100 of file G4InuclEvaporation.cc.

References CLHEP::HepLorentzVector::boost(), G4VCascadeDeexcitation::collide(), G4cout, G4endl, G4Fragment::GetA_asInt(), G4Fragment::GetExcitationEnergy(), G4Fragment::GetMomentum(), G4NucleiProperties::GetNuclearMass(), G4CollisionOutput::getOutgoingNuclei(), G4CollisionOutput::getOutgoingParticles(), CLHEP::HepLorentzVector::getZ(), G4Fragment::GetZ_asInt(), python.hepunit::GeV, G4InuclNuclei::setExitationEnergy(), and CLHEP::HepLorentzVector::vect().

                                                                            {
  G4FragmentVector* theResult = new G4FragmentVector;

  if (theNucleus.GetExcitationEnergy() <= 0.0) { // Check that Excitation Energy > 0
    theResult->push_back(new G4Fragment(theNucleus));
    return theResult;
  }

  G4int A = theNucleus.GetA_asInt();
  G4int Z = theNucleus.GetZ_asInt();
  G4double mTar  = G4NucleiProperties::GetNuclearMass(A, Z); // Mass of the target nucleus

  G4ThreeVector momentum =  theNucleus.GetMomentum().vect();
  G4double exitationE = theNucleus.GetExcitationEnergy();

  G4double mass = mTar;
  G4ThreeVector boostToLab( momentum/mass ); 

  if ( verboseLevel > 2 )
    G4cout << " G4InuclEvaporation : initial kinematics : boostToLab vector = " << boostToLab << G4endl
       << "                     excitation energy  : " << exitationE << G4endl;

  if (verboseLevel > 2) {
    G4cout << "G4InuclEvaporation::BreakItUp >>> A: " << A << " Z: " << Z
       << " exitation E: " << exitationE << " mass: " << mTar/GeV << " GeV"
       << G4endl;
  };

  G4InuclNuclei* nucleus = new G4InuclNuclei(A, Z);
  nucleus->setExitationEnergy(exitationE);

  G4CollisionOutput output;
  evaporator->collide(0, nucleus, output);

  const std::vector<G4InuclNuclei>& outgoingNuclei = output.getOutgoingNuclei();
  const std::vector<G4InuclElementaryParticle>& particles = output.getOutgoingParticles();

  G4double eTot=0.0;
  G4int  i=1;

  if (!particles.empty()) { 
    G4int outgoingType;
    particleIterator ipart = particles.begin();
    for (; ipart != particles.end(); ipart++) {
      outgoingType = ipart->type();

      if (verboseLevel > 2) {
    G4cout << "Evaporated particle:  " << i << " of type: "
           << outgoingType << G4endl;
        i++;
      }

      eTot += ipart->getEnergy();
      
      G4LorentzVector vlab = ipart->getMomentum().boost(boostToLab);

      theResult->push_back( new G4Fragment(vlab, ipart->getDefinition()) );
    }  
  }

  //  G4cout << "# fragments " << output.getOutgoingNuclei().size() << G4endl;
  if (!outgoingNuclei.empty()) { 
    nucleiIterator ifrag = outgoingNuclei.begin();
    for (i=1; ifrag != outgoingNuclei.end(); ifrag++) {
      if (verboseLevel > 2) {
    G4cout << " Nuclei fragment: " << i << G4endl; i++;
      }

      eTot += ifrag->getEnergy();

      G4LorentzVector vlab = ifrag->getMomentum().boost(boostToLab);
 
      G4int fragA = ifrag->getA();
      G4int fragZ = ifrag->getZ();
      if (verboseLevel > 2) {
    G4cout << "boosted v" << vlab << G4endl;
      }
      theResult->push_back( new G4Fragment(fragA, fragZ, vlab) ); 
    }
  }

  //G4cout << ">>>> G4InuclEvaporation::BreakItUp end " << G4endl;
  return theResult;
}

void G4InuclEvaporation::setVerboseLevel

(

const G4int

verbose

)

Definition at line 96 of file G4InuclEvaporation.cc.

                                                              {
  verboseLevel = verbose;
}

---

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

• G4InuclEvaporation.hh
• G4InuclEvaporation.cc