#include <G4UnstableFragmentBreakUp.hh>

Inheritance diagram for G4UnstableFragmentBreakUp:

Public Member Functions
	G4UnstableFragmentBreakUp ()

virtual	~G4UnstableFragmentBreakUp ()

virtual G4FragmentVector *	BreakUpFragment (G4Fragment *fragment)

virtual G4Fragment *	EmittedFragment (G4Fragment *fragment)

virtual G4FragmentVector *	BreakUp (const G4Fragment &fragment)

virtual G4double	GetEmissionProbability (G4Fragment *fragment)

void	SetVerboseLevel (G4int val)

Public Member Functions inherited from G4VEvaporationChannel
	G4VEvaporationChannel (const G4String &aName="Anonymous", G4EvaporationChannelType timeType=fDelayedEmission)

virtual	~G4VEvaporationChannel ()

virtual G4double	GetLifeTime (G4Fragment *theNucleus)

G4String	GetName () const

void	SetName (const G4String &aName)

void	SetOPTxs (G4int opt)

void	UseSICB (G4bool use)

Additional Inherited Members
Protected Attributes inherited from G4VEvaporationChannel
G4EvaporationChannelType	sampleDecayTime

G4int	OPTxs

G4bool	useSICB

Detailed Description

Definition at line 57 of file G4UnstableFragmentBreakUp.hh.

Constructor & Destructor Documentation

G4UnstableFragmentBreakUp::G4UnstableFragmentBreakUp ( )

Definition at line 53 of file G4UnstableFragmentBreakUp.cc.

References test::a, G4NucleiProperties::GetNuclearMass(), G4NistManager::Instance(), and z.

   :verbose(0)
 { 
   fNistManager = G4NistManager::Instance();
   const G4int z[6] = {0, 1, 1, 1, 2, 2};
   const G4int a[6] = {1, 1, 2, 3, 3, 4};
   for(G4int i=0; i<6; ++i) {
     Zfr[i] = z[i];
     Afr[i] = a[i];
     masses[i] = G4NucleiProperties::GetNuclearMass(a[i], z[i]);
   }
 }

G4UnstableFragmentBreakUp::~G4UnstableFragmentBreakUp ( )

virtual

Definition at line 66 of file G4UnstableFragmentBreakUp.cc.

67 {}

Member Function Documentation

G4FragmentVector * G4UnstableFragmentBreakUp::BreakUp ( const G4Fragment & fragment )

virtual

Implements G4VEvaporationChannel.

Definition at line 148 of file G4UnstableFragmentBreakUp.cc.

 {
   return 0;
 }

G4FragmentVector * G4UnstableFragmentBreakUp::BreakUpFragment ( G4Fragment * fragment )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 74 of file G4UnstableFragmentBreakUp.cc.

References CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), G4UniformRand, G4Fragment::GetA_asInt(), G4Fragment::GetCreationTime(), G4Fragment::GetMomentum(), G4NucleiProperties::GetNuclearMass(), G4Fragment::GetZ_asInt(), CLHEP::HepLorentzVector::mag(), G4Fragment::SetCreationTime(), G4Fragment::SetMomentum(), G4Fragment::SetZandA_asInt(), and python.hepunit::twopi.

Referenced by G4Evaporation::BreakItUp().

 {
   //G4cout << "G4UnstableFragmentBreakUp::BreakUpFragment" << G4endl;
   G4FragmentVector * theResult = new G4FragmentVector();
 
   G4int Z = nucleus->GetZ_asInt();
   G4int A = nucleus->GetA_asInt();
   G4int Amax = A;
   G4LorentzVector lv = nucleus->GetMomentum();
   G4double time = nucleus->GetCreationTime();
 
   G4double deltaE, mass, mass1(0.0), mass2(0.0);
   G4int i, index;
 
   // Starts loop over evaporated particles, loop is limited by number
   // of nucleons
   for(G4int ia=0; ia<Amax; ++ia) {
  
     mass = lv.mag();
     deltaE = 0.0;
     index  = -1;
     for(i=0; i<6; ++i) {
       G4int Zres = Z - Zfr[i];
       G4int Ares = A - Afr[i];
       if(Zres >= 0 && Ares >= Zres && Ares > 0) {
     G4double m1 = G4NucleiProperties::GetNuclearMass(Ares, Zres);
     G4double de = mass - m1 - masses[i];
         if(de > deltaE) {
       mass1 = m1;
           mass2 = masses[i];
           deltaE= de;
           index = i;
     }
       }
     }
 
     // no decay channels
     if(index < 0) { break; }
 
     // compute energy of light fragment
     G4double e2 = 0.5*((mass - mass1)*(mass + mass1) + mass2*mass2)/mass;
     if(e2 < mass2) { break; }
 
     // sample decay
     G4ThreeVector bst  = lv.boostVector();
 
     G4double cosTheta = 1. - 2. * G4UniformRand(); 
     G4double sinTheta = std::sqrt(1. - cosTheta * cosTheta);
     G4double phi = twopi * G4UniformRand();
     G4double mom = std::sqrt((e2 - mass2)*(e2 + mass2));
     G4LorentzVector mom2(mom * sinTheta * std::cos(phi),
              mom * sinTheta * std::sin(phi),
              mom * cosTheta,
              e2);
     mom2.boost(bst);  
     G4Fragment* fr = new G4Fragment(Afr[index], Zfr[index], mom2);
     fr->SetCreationTime(time);
     theResult->push_back(fr);
 
     // residual
     lv -= mom2;
     Z  -= Zfr[index];
     A  -= Afr[index];
   }
 
   // updated fragment
   if( theResult->size() > 0) {
     nucleus->SetZandA_asInt(Z, A);
     nucleus->SetMomentum(lv);
   }
 
   return theResult;
 }