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 a5e77276ac03c3209cdf18b2e7dd71347d72d1a3..947bf64ba7c7f9b85f52484c4652873caed3ffa0 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 ba751f8977340211ff1cfb68a8debec812861ae6..ca2f5a6e7bc25929b6df0378ad446548b3eecce0 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 83f001fcb1a7ed7e79eee271a7efe4344c6c97a3..af6a66b958cdf396c6afb61576fb39d3c51b957e 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 0f13bdf6186ba527fdff933c892bf652019c4e8f..50198d7cb046b620046e59c40a859af5118adba8 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