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 G4FragmentVector *	BreakItUp (const G4Fragment &theNucleus)=0
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::G4InuclEvaporation

(

)

Definition at line 72 of file G4InuclEvaporation.cc.

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

~G4InuclEvaporation()

G4InuclEvaporation::~G4InuclEvaporation

(

)

Definition at line 79 of file G4InuclEvaporation.cc.

                                        {
  delete evaporator;
}

Member Function Documentation

BreakItUp()

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

virtual

Implements G4VEvaporation.

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

setVerboseLevel()

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