Setter for Pythia Tune.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 17 Aug 2009 15:49:15 +0000 (15:49 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 17 Aug 2009 15:49:15 +0000 (15:49 +0000)
PYTHIA6/AliGenPythia.cxx
PYTHIA6/AliGenPythia.h
PYTHIA6/AliPythia.cxx
PYTHIA6/AliPythia.h

index a5e7727..947bf64 100644 (file)
@@ -52,6 +52,7 @@ ClassImp(AliGenPythia)
 AliGenPythia::AliGenPythia():
     AliGenMC(),
     fProcess(kPyCharm),          
+    fItune(-1),
     fStrucFunc(kCTEQ5L), 
     fKineBias(0.),
     fTrials(0),
@@ -147,6 +148,7 @@ AliGenPythia::AliGenPythia():
 AliGenPythia::AliGenPythia(Int_t npart)
     :AliGenMC(npart),
      fProcess(kPyCharm),          
+     fItune(-1),
      fStrucFunc(kCTEQ5L), 
      fKineBias(0.),
      fTrials(0),
@@ -373,7 +375,7 @@ void AliGenPythia::Init()
        fRL = 0x0;
     }
  //
-    fPythia->ProcInit(fProcess,fEnergyCMS,fStrucFunc);
+    fPythia->ProcInit(fProcess,fEnergyCMS,fStrucFunc, fItune);
     //  Forward Paramters to the AliPythia object
     fDecayer->SetForceDecay(fForceDecay);    
 // Switch off Heavy Flavors on request  
index ba751f8..ca2f5a6 100644 (file)
@@ -195,6 +195,7 @@ class AliGenPythia : public AliGenMC
     void     MakeHeader();    
     void     GeneratePileup();
     Process_t   fProcess;           //Process type
+    Int_t       fItune;             // Pythia tune > 6.4
     StrucFunc_t fStrucFunc;         //Structure Function
     Float_t     fKineBias;          //!Bias from kinematic selection
     Int_t       fTrials;            //!Number of trials for current event
index 83f001f..af6a66b 100644 (file)
@@ -96,7 +96,7 @@ AliPythia::AliPythia(const AliPythia& pythia):
     pythia.Copy(*this);
 }
 
-void AliPythia::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfunc)
+void AliPythia::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfunc, Int_t itune)
 {
 // Initialise the process to generate 
     if (!AliPythiaRndm::GetPythiaRandom()) 
@@ -579,6 +579,11 @@ void AliPythia::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfun
     }
 //
 //  Initialize PYTHIA
+//
+//  Select the tune
+    if (itune > -1) Pytune(itune);
+    
+//  
     SetMSTP(41,1);   // all resonance decays switched on
     Initialize("CMS","p","p",fEcms);
     
index 0f13bdf..50198d7 100644 (file)
@@ -22,7 +22,7 @@ class AliPythia : public TPythia6, public AliRndm
     virtual Int_t CheckedLuComp(Int_t kf);
     // Pythia initialisation for selected processes
     virtual void ProcInit
-       (Process_t process, Float_t energy, StrucFunc_t strucfunc);
+       (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);
     // Print particle properties