option to switch CR off
authormorsch <andreas.morsch@cern.ch>
Fri, 25 Apr 2014 13:16:42 +0000 (15:16 +0200)
committermorsch <andreas.morsch@cern.ch>
Fri, 25 Apr 2014 13:16:42 +0000 (15:16 +0200)
PYTHIA6/AliGenPythia.cxx
PYTHIA6/AliGenPythia.h

index f50efaa..3de870a 100644 (file)
@@ -76,6 +76,7 @@ AliGenPythia::AliGenPythia():
     fYHardMax(1.e10),
     fGinit(1),
     fGfinal(1),
+    fCRoff(0),
     fHadronisation(1),
     fPatchOmegaDalitz(0), 
     fNpartons(0),
@@ -192,6 +193,7 @@ AliGenPythia::AliGenPythia(Int_t npart)
      fYHardMax(1.e10),
      fGinit(kTRUE),
      fGfinal(kTRUE),
+     fCRoff(kFALSE),
      fHadronisation(kTRUE),
      fPatchOmegaDalitz(0), 
      fNpartons(0),
@@ -402,6 +404,8 @@ void AliGenPythia::Init()
     fPythia->SetMSTP(61,fGinit);
 //  final state radiation
     fPythia->SetMSTP(71,fGfinal);
+    //color reconnection strength
+    if(fCRoff==1)fPythia->SetMSTP(95,0);
 //  pt - kick
     if (fPtKick > 0.) {
        fPythia->SetMSTP(91,1);
index 0c82e8c..e888bba 100644 (file)
@@ -57,7 +57,8 @@ class AliGenPythia : public AliGenMC
     // Set initial and final state gluon radiation
     virtual void    SetGluonRadiation(Int_t iIn, Int_t iFin)
        {fGinit = iIn; fGfinal = iFin;}
-    // Intrinsic kT
+    virtual void SetColorReconnectionOff(Int_t iflag=0){fCRoff=iflag;} 
+   // Intrinsic kT
     virtual void    SetPtKick(Float_t kt = 1.)
        {fPtKick = kt;}
     // Use the Pythia 6.3 new multiple interations scenario
@@ -280,6 +281,7 @@ class AliGenPythia : public AliGenMC
     Float_t     fYHardMax;          //higher y-hard cut
     Int_t       fGinit;             //initial state gluon radiation
     Int_t       fGfinal;            //final state gluon radiation
+    Int_t       fCRoff;             //color reconnection off in the pythia6 annealying model
     Int_t       fHadronisation;     //hadronisation
     Bool_t      fPatchOmegaDalitz;  //flag for omega dalitz decay patch
     Int_t       fNpartons;          //Number of partons before hadronisation