#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	BreakFragment (G4FragmentVector , G4Fragment theNucleus)

virtual void	InitialiseChannels ()

virtual void	SetPhotonEvaporation (G4VEvaporationChannel *ptr)

void	SetFermiBreakUp (G4VFermiBreakUp *ptr)

G4VFermiBreakUp *	GetFermiBreakUp () const

G4VEvaporationChannel *	GetPhotonEvaporation () const

G4VEvaporationChannel *	GetFissionChannel () const

G4VEvaporationChannel *	GetChannel (std::size_t idx) const

void	SetOPTxs (G4int opt)

void	UseSICB (G4bool use)

std::size_t	GetNumberOfChannels () const

	G4VEvaporation (const G4VEvaporation &right)=delete

const G4VEvaporation &	operator= (const G4VEvaporation &right)=delete

G4bool	operator== (const G4VEvaporation &right) const =delete

G4bool	operator!= (const G4VEvaporation &right) const =delete

Additional Inherited Members
Protected Member Functions inherited from G4VEvaporation
void	CleanChannels ()

Protected Attributes inherited from G4VEvaporation
G4VEvaporationChannel *	thePhotonEvaporation {nullptr}

G4VFermiBreakUp *	theFBU {nullptr}

G4int	OPTxs {3}

G4bool	useSICB {true}

std::vector< G4VEvaporationChannel * > *	theChannels {nullptr}

G4VEvaporationFactory *	theChannelFactory {nullptr}

Detailed Description

Definition at line 43 of file G4InuclEvaporation.hh.

Constructor & Destructor Documentation

◆ G4InuclEvaporation()

G4InuclEvaporation::G4InuclEvaporation ( )

Definition at line 72 of file G4InuclEvaporation.cc.

73 : 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 )

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();
 
  G4int  i=1;
 
  if (!particles.empty()) { 
    G4int outgoingType;
    particleIterator ipart = particles.cbegin();
    for (; ipart != particles.cend(); ++ipart) {
      outgoingType = ipart->type();
 
      if (verboseLevel > 2) {
    G4cout << "Evaporated particle:  " << i << " of type: "
           << outgoingType << G4endl;
        ++i;
      }
 
      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.cbegin();
    for (i=1; ifrag != outgoingNuclei.cend(); ++ifrag) {
      if (verboseLevel > 2) {
    G4cout << " Nuclei fragment: " << i << G4endl; ++i;
      }
 
      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;
}

◆ 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:

Public Member Functions