G4InuclEvaporation Class Reference

#include <G4InuclEvaporation.hh>

Inheritance diagram for G4InuclEvaporation:

Collaboration 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	BreakFragment (G4FragmentVector *, G4Fragment *theNucleus)
virtual void	InitialiseChannels ()
virtual void	SetPhotonEvaporation (G4VEvaporationChannel *ptr)
void	SetFermiBreakUp (G4VFermiBreakUp *ptr)
G4VFermiBreakUp *	GetFermiBreakUp () const
G4VEvaporationChannel *	GetPhotonEvaporation ()
G4VEvaporationChannel *	GetFissionChannel ()
void	SetOPTxs (G4int opt)
void	UseSICB (G4bool use)
size_t	GetNumberOfChannels () const

Additional Inherited Members
Protected Member Functions inherited from G4VEvaporation
void	CleanChannels ()
Protected Attributes inherited from G4VEvaporation
G4VEvaporationChannel *	thePhotonEvaporation
G4VFermiBreakUp *	theFBU
G4int	OPTxs
G4bool	useSICB
std::vector < G4VEvaporationChannel * > *	theChannels
G4VEvaporationFactory *	theChannelFactory

Detailed Description

Definition at line 44 of file G4InuclEvaporation.hh.

Constructor & Destructor Documentation

G4InuclEvaporation::G4InuclEvaporation

(

)

Definition at line 73 of file G4InuclEvaporation.cc.

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

G4InuclEvaporation::~G4InuclEvaporation

(

)

Definition at line 80 of file G4InuclEvaporation.cc.

                                        {
  delete evaporator;
}

Member Function Documentation

G4FragmentVector * G4InuclEvaporation::BreakItUp

(

const G4Fragment &

theNucleus

)

Definition at line 101 of file G4InuclEvaporation.cc.

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

      // TEMPORARY:  Remove constness on PD until G4Fragment is fixed
      theResult->push_back( new G4Fragment(vlab, ipart->getDefinition()) );
    }  
  }

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

  return theResult;
}

Here is the call graph for this function:

void G4InuclEvaporation::setVerboseLevel

(

const G4int

verbose

)

Definition at line 97 of file G4InuclEvaporation.cc.

                                                              {
  verboseLevel = verbose;
}

---

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

• geant4.10.03.p01/source/processes/hadronic/models/cascade/cascade/include/G4InuclEvaporation.hh
• geant4.10.03.p01/source/processes/hadronic/models/cascade/cascade/src/G4InuclEvaporation.cc