Possibility to switch off heavy flavor production added.

diff --git a/THijing/AliGenHijing.cxx b/THijing/AliGenHijing.cxx
index 2cd1c1500e151f84638230ea20b49a0c855aa127..f456abceef2779c3ac6039b02c726c53f5a3ba3d 100644 (file)
--- a/THijing/AliGenHijing.cxx
+++ b/THijing/AliGenHijing.cxx
@@ -76,6 +76,8 @@ AliGenHijing::AliGenHijing(Int_t npart)
     //
     SetSimpleJets();
     SetNoGammas();
     //
     SetSimpleJets();
     SetNoGammas();

+    SetRandomPz();
+    SwitchOffHeavyQuarks(kFALSE);

 //
     fParticles = new TClonesArray("TParticle",10000);    
 //
 //
     fParticles = new TClonesArray("TParticle",10000);    
 //


@@ -151,6 +153,15 @@ void AliGenHijing::Init()
        fHijing->SetHIPR1(11, 2.5);
     }
     
        fHijing->SetHIPR1(11, 2.5);
     }
     

+//
+// Heavy quarks
+//    
+    if (fNoHeavyQuarks) {
+       fHijing->SetIHPR2(49, 1);
+    } else {
+       fHijing->SetIHPR2(49, 0);
+    }
+    

     
     AliGenMC::Init();
     
     
     AliGenMC::Init();
     

diff --git a/THijing/AliGenHijing.h b/THijing/AliGenHijing.h
index d930efe486da4d42e3ff98b0386f555b72b98f97..6d341b0c4288060b9ff1ef3ed149accf21726e2a 100644 (file)
--- a/THijing/AliGenHijing.h
+++ b/THijing/AliGenHijing.h
@@ -55,6 +55,9 @@ class AliGenHijing : public AliGenMC
        {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;}
        {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;}
 // Getters
     virtual Float_t GetEnergyCMS()       const {return fEnergyCMS;}
     virtual TString GetReferenceFrame()  const {return fFrame;}


@@ -125,7 +128,9 @@ class AliGenHijing : public AliGenMC
     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      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
  private:
     void Copy(TObject &rhs) const;
     // adjust the weight from kinematic cuts