64 #include "HistoManager.hh"
76 G4cout <<
" >>> G4HIJING_Model default constructor" <<
G4endl;
100 InitialiseDataTables();
189 if (AP1>1 ||definitionP==anti_deu ||definitionP==anti_he3
190 ||definitionP==anti_tri ||definitionP==anti_alp)
197 }
else if (
id==2212) {
203 }
else if(
id==-2212){
213 }
else if(
id==-2112){
224 }
else if(
id==-211) {
235 }
else if(
id==-321) {
242 G4cout <<
" Sorry, No definition for PROJECTLE for HIJING::"
246 #if G4VERSION_NUMBER>=950
250 "Sorry, no definition for PROJECTILE for HIJING");
318 Psum.boost(0.0,0.0,-1.0*beta_rel);
374 G4cout <<
"===============Warning====================================="<<
G4endl;
375 G4cout <<
"-----------------------------------------------------------"<<
G4endl;
377 G4cout <<
"------------------------------------------------------------"<<
G4endl;
378 G4cout <<
"============================================================"<<
G4endl;
425 lorenzCM.boost(0.0,0.0,-1.0*betann);
427 -1.0*lorenzCM.pz(),lorenzCM.e());
437 #ifdef G4ANALYSIS_USE //khaled new
455 for (
G4int j=0; j<nSecondaries; j++) {
457 GetParticle()->GetTotalEnergy()/
GeV;
460 GetParticle()->GetMomentum()/
GeV;
463 GetParticle()->GetDefinition();
470 G4cout <<
"----------------------------------------"
471 <<
"----------------------------------------"
473 G4cout <<
"Total energy before collision = " <<TotalEbefore
475 G4cout <<
"Total energy after collision = " <<TotalEafter/nSecondaries
479 G4cout <<
"----------------------------------------"<<
G4endl;
481 G4cout <<
"Total momentum before collision = " <<TotalPbefore
484 G4cout <<
"Total momentum after collision = " <<TotalPafter.z()/nSecondaries
487 G4cout <<
"----------------------------------------"<<
G4endl;
490 G4cout <<
"Total charge before collision = " <<(ZP+ZT)
492 G4cout <<
"Total charge after collision = " <<charge
495 G4cout <<
"----------------------------------------"<<
G4endl;
497 G4cout <<
"Total baryon number before collision = "<<AP+AT
499 G4cout <<
"Total baryon number after collision = "<<baryon
501 G4cout <<
"----------------------------------------"<<
G4endl;
505 G4cout <<
"----------------------------------------"
506 <<
"----------------------------------------"
526 G4cout <<
" *****************************************************************"
528 G4cout <<
" Interface to G4HIJING_Model activated"
530 G4cout <<
" Version number : 01.00.0B File date : 10/09/2013" <<
G4endl;
531 G4cout <<
" Interface written by Khaled Abdel-Waged "
533 G4cout <<
" Umm Al-Qura University "
538 G4cout <<
" *****************************************************************"
static G4AntiHe3 * AntiHe3()
virtual G4HadFinalState * ApplyYourself(const G4HadProjectile &, G4Nucleus &)
G4ParticleDefinition * FindParticle(G4int PDGEncoding)
G4HadSecondary * GetSecondary(size_t i)
CLHEP::Hep3Vector G4ThreeVector
static G4AntiDeuteron * AntiDeuteron()
void WelcomeMessage() const
G4int GetPDGEncoding() const
static G4AntiAlpha * AntiAlpha()
virtual ~G4HIJING_Model()
struct cchijinghimain1 himain1_
void SetStatusChange(G4HadFinalStateStatus aS)
G4HIJING_Model(const G4String &name="HIJING")
std::vector< G4ReactionProduct * > G4ReactionProductVector
G4GLOB_DLL std::ostream G4cout
const G4ParticleDefinition * GetDefinition() const
struct cchijinghimain2 himain2_
G4double Eplab(G4double, G4double)
struct cchijinghiparnt hiparnt_
G4ErrorTarget * theTarget
const G4LorentzVector & Get4Momentum() const
G4ReactionProductVector * Propagate(G4KineticTrackVector *theSecondaries, G4V3DNucleus *theTarget)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
G4HadFinalState theResult
static G4ParticleTable * GetParticleTable()
T max(const T t1, const T t2)
brief Return the largest of the two arguments
void hijing_(float *, float *)
void AddSecondary(G4DynamicParticle *aP, G4int mod=-1)
void SetMinEnergy(G4double aM)
static HistoManager * GetPointer()
static G4AntiTriton * AntiTriton()
G4double GetPDGCharge() const
HistoManager * fHistoManager
void InitialiseDataTables()
G4int GetNumberOfSecondaries() const
G4int GetBaryonNumber() const
CLHEP::HepLorentzVector G4LorentzVector
void g4hijingblockdata_()