//
SetSimpleJets();
SetNoGammas();
+ SetRandomPz();
+ SwitchOffHeavyQuarks(kFALSE);
//
fParticles = new TClonesArray("TParticle",10000);
//
fHijing->SetHIPR1(11, 2.5);
}
+//
+// Heavy quarks
+//
+ if (fNoHeavyQuarks) {
+ fHijing->SetIHPR2(49, 1);
+ } else {
+ fHijing->SetIHPR2(49, 0);
+ }
+
AliGenMC::Init();
{fPhiMinJet = TMath::Pi()*phimin/180.; fPhiMaxJet = TMath::Pi()*phimax/180.;}
virtual void SetBoostLHC(Int_t flag = 0) {fLHC = flag;}
virtual void SetRandomPz(Bool_t flag = 0) {fRandomPz = flag;}
+ virtual void SwitchOffHeavyQuarks(Bool_t flag = kTRUE) {fNoHeavyQuarks = flag;}
+
+
// Getters
virtual Float_t GetEnergyCMS() const {return fEnergyCMS;}
virtual TString GetReferenceFrame() const {return fFrame;}
Int_t fTargetSpecn; // Num. of spectator neutrons from target nucleus
Int_t fTargetSpecp; // Num. of spectator protons from target nucleus
Int_t fLHC; // Assume LHC as lab frame
- Int_t fRandomPz; // Randomise sign of pz event by event
+ Bool_t fRandomPz; // Randomise sign of pz event by event
+ Bool_t fNoHeavyQuarks; // If true no heavy quarks are produced
+
private:
void Copy(TObject &rhs) const;
// adjust the weight from kinematic cuts