Default is pp.
Zaida
fPythia->SetCKIN(7,fYHardMin);
fPythia->SetCKIN(8,fYHardMax);
- if (fAProjectile > 0 && fATarget > 0) fPythia->SetNuclei(fAProjectile, fATarget, fNucPdf);
+ if (fAProjectile > 0 && fATarget > 0) fPythia->SetNuclei(fAProjectile, fATarget, fNucPdf);
+ if (fProjectile != "p" || fTarget != "p") fPythia->SetCollisionSystem(fProjectile,fTarget);
+
// Fragmentation?
if (fFragmentation) {
fPythia->SetMSTP(111,1);
virtual void SetEnergyCMS(Float_t energy = 5500) {fEnergyCMS = energy;}
// Treat protons as inside nuclei with mass numbers a1 and a2
virtual void SetNuclei(Int_t a1, Int_t a2, Int_t pdfset = 0);
+ // Set colliding nuclei ("p","n",...)
+ virtual void SetCollisionSystem(TString projectile, TString target) { fProjectile = projectile; fTarget = target; }
virtual void SetNuclearPDF(Int_t pdf) {fNucPdf = pdf;}
//
// Trigger options
fProcess(kPyMb),
fEcms(0.),
fStrucFunc(kCTEQ5L),
+ fProjectile("p"),
+ fTarget("p"),
fXJet(0.),
fYJet(0.),
fNGmax(30),
fProcess(kPyMb),
fEcms(0.),
fStrucFunc(kCTEQ5L),
+ fProjectile("p"),
+ fTarget("p"),
fXJet(0.),
fYJet(0.),
fNGmax(30),
//
SetMSTP(41,1); // all resonance decays switched on
- Initialize("CMS","p","p",fEcms);
+ Initialize("CMS",fProjectile,fTarget,fEcms);
fOmegaDalitz.Init();
}
(Process_t process, Float_t energy, StrucFunc_t strucfunc, Int_t itune = -1);
// treat protons as inside nuclei
virtual void SetNuclei(Int_t a1, Int_t a2, Int_t pdf);
+ // Set colliding nuclei ("p","n",...)
+ virtual void SetCollisionSystem(TString projectile, TString target) { fProjectile = projectile; fTarget = target; }
// Print particle properties
virtual void PrintParticles();
virtual void ResetDecayTable();
Process_t fProcess; // Process type
Float_t fEcms; // Centre of mass energy
StrucFunc_t fStrucFunc; // Structure function
+ TString fProjectile; // Projectile
+ TString fTarget; // Target
Int_t fDefMDCY[501]; // ! Default decay switches per particle
Int_t fDefMDME[2001]; // ! Default decay switches per mode
Double_t fZQuench[4]; // ! Quenching fractions for this even
fProcess(kPyMb),
fEcms(0.),
fStrucFunc(kCTEQ5L),
+ fProjectile("p"),
+ fTarget("p"),
fXJet(0.),
fYJet(0.),
fNGmax(30),
fProcess(kPyMb),
fEcms(0.),
fStrucFunc(kCTEQ5L),
+ fProjectile("p"),
+ fTarget("p"),
fXJet(0.),
fYJet(0.),
fNGmax(30),
//
// Initialize PYTHIA
SetMSTP(41,1); // all resonance decays switched on
- Initialize("CMS","p","p",fEcms);
+ Initialize("CMS",fProjectile,fTarget,fEcms);
}
virtual Int_t GetParticles(TClonesArray *particles) {return ImportParticles(particles, "All");}
// Treat protons as inside nuclei
virtual void SetNuclei(Int_t a1, Int_t a2);
+ // Set colliding nuclei ("p","n",...)
+ virtual void SetCollisionSystem(TString projectile, TString target) { fProjectile = projectile; fTarget = target; }
// Print particle properties
virtual void PrintParticles();
// Reset the decay table
Process_t fProcess; // Process type
Float_t fEcms; // Centre of mass energy
StrucFunc_t fStrucFunc; // Structure function
+ TString fProjectile; // Projectile
+ TString fTarget; // Target
Int_t fDefMDCY[501]; // ! Default decay switches per particle
Int_t fDefMDME[2001]; // ! Default decay switches per mode
Double_t fZQuench[4]; // ! Quenching fractions for this even