- changes by pedro for dalitz

- changes in AddTask for pPb

diff --git a/PWGGA/GammaConv/AliAnalysisTaskGammaConvDalitzV1.cxx b/PWGGA/GammaConv/AliAnalysisTaskGammaConvDalitzV1.cxx
index fe699d4e382879e05aae3174a68423ddc74a7cc1..f17542bed8adcacafe0fd8cc6ddd3b0c1fc6b3ae 100644 (file)
--- a/PWGGA/GammaConv/AliAnalysisTaskGammaConvDalitzV1.cxx
+++ b/PWGGA/GammaConv/AliAnalysisTaskGammaConvDalitzV1.cxx
@@ -517,6 +517,8 @@ void AliAnalysisTaskGammaConvDalitzV1::UserCreateOutputObjects()
    hESDDalitzPositronAfterTPCdEdxVsPhi = new TH2F*[fnCuts];
    hESDMotherPhi                   = new TH1F*[fnCuts];
    hESDEposEnegPsiPairDPhi         = new TH2F*[fnCuts];
+   hESDEposEnegInvMassPt           = new TH2F*[fnCuts];
+   hESDEposEnegLikeSignBackInvMassPt = new TH2F*[fnCuts];
    hESDConvGammaZR                 = new TH2F*[fnCuts];
    
    }
@@ -526,8 +528,6 @@ void AliAnalysisTaskGammaConvDalitzV1::UserCreateOutputObjects()
    hESDMotherInvMassPt             = new TH2F*[fnCuts];
    
    if(fDoChicAnalysis) {
-   hESDEposEnegInvMassPt             = new TH2F*[fnCuts];
-   hESDEposEnegLikeSignBackInvMassPt = new TH2F*[fnCuts];  
    hESDPi0MotherInvMassPt          = new TH2F*[fnCuts];
    hESDPi0MotherDiffInvMassPt      = new TH2F*[fnCuts];
    hESDPi0MotherDiffLimInvMassPt   = new TH2F*[fnCuts];
@@ -633,7 +633,7 @@ void AliAnalysisTaskGammaConvDalitzV1::UserCreateOutputObjects()
         
       
         
-      hESDConvGammaZR[iCut]= new TH2F("ESD_ConvGamma_ConversionPoint_ZR","ESD_ConvGamma_ConversionPoint_ZR",600,-150,150,240,0,120);
+      hESDConvGammaZR[iCut]= new TH2F("ESD_ConvGamma_ConversionPoint_ZR","ESD_ConvGamma_ConversionPoint_ZR",1200,-150,150,480,0,120);
       fQAFolder[iCut]->Add(hESDConvGammaZR[iCut]);
      
       hESDDalitzElectronAfterPt[iCut] = new TH1F("ESD_DalitzElectron_After_Pt","ESD_DalitzElectron_After_Pt",1000,0,25);
@@ -713,11 +713,11 @@ void AliAnalysisTaskGammaConvDalitzV1::UserCreateOutputObjects()
       hESDEposEnegPsiPairDPhi[iCut] = new TH2F("ESD_EposEneg_PsiPair_DPhi","ESD_EposEneg_PsiPair_DPhi", 100, -1.0,1.0,100,-1.0,1.0 );
       fQAFolder[iCut]->Add(hESDEposEnegPsiPairDPhi[iCut]);
 
-      //hESDEposEnegInvMassPt[iCut] = new TH2F("ESD_EposEneg_InvMassPt","ESD_EposEneg_InvMassPt",5000,0.,5.,100,0.,10.);
-      //fQAFolder[iCut]->Add(hESDEposEnegInvMassPt[iCut]);
+      hESDEposEnegInvMassPt[iCut] = new TH2F("ESD_EposEneg_InvMassPt","ESD_EposEneg_InvMassPt",5000,0.,5.,100,0.,10.);
+      fQAFolder[iCut]->Add(hESDEposEnegInvMassPt[iCut]);
       
-      //hESDEposEnegLikeSignBackInvMassPt[iCut]  = new TH2F("ESD_EposEneg_LikeSignBack_InvMassPt","ESD_EposEneg_LikeSignBack_InvMassPt",5000,0.,5.,100,0.,10.);
-      //fQAFolder[iCut]->Add(hESDEposEnegLikeSignBackInvMassPt[iCut]);
+      hESDEposEnegLikeSignBackInvMassPt[iCut]  = new TH2F("ESD_EposEneg_LikeSignBack_InvMassPt","ESD_EposEneg_LikeSignBack_InvMassPt",5000,0.,5.,100,0.,10.);
+      fQAFolder[iCut]->Add(hESDEposEnegLikeSignBackInvMassPt[iCut]);
       
       
       TAxis *AxisAfter = hESDDalitzElectronAfterTPCdEdxVsP[iCut]->GetXaxis(); 
@@ -768,14 +768,6 @@ void AliAnalysisTaskGammaConvDalitzV1::UserCreateOutputObjects()
       hESDPi0MotherDiffLimInvMassPt[iCut] = new TH2F("ESD_Pi0Mother_DiffLimInvMass_Pt","ESD_Pi0Mother_DiffLimInvMass_Pt",2000,0,2,250,0,25);
       fESDList[iCut]->Add(hESDPi0MotherDiffLimInvMassPt[iCut]);
       
-      hESDEposEnegInvMassPt[iCut] = new TH2F("ESD_EposEneg_InvMassPt","ESD_EposEneg_InvMassPt",5000,0.,5.,100,0.,10.);
-      fESDList[iCut]->Add(hESDEposEnegInvMassPt[iCut]);
-      
-      hESDEposEnegLikeSignBackInvMassPt[iCut]  = new TH2F("ESD_EposEneg_LikeSignBack_InvMassPt","ESD_EposEneg_LikeSignBack_InvMassPt",5000,0.,5.,100,0.,10.);
-      fESDList[iCut]->Add(hESDEposEnegLikeSignBackInvMassPt[iCut]);
-      
-      
-      
       }
 
 
@@ -1681,9 +1673,12 @@ void AliAnalysisTaskGammaConvDalitzV1::ProcessElectronCandidates(){
         
       
         
-        if( fDoMesonQA ) hESDEposEnegPsiPairDPhi[fiCut]->Fill(deltaPhi,psiPair);        
-        if( fDoChicAnalysis ) hESDEposEnegInvMassPt[fiCut]->Fill(vphoton->GetMass(),vphoton->Pt());
+        if ( fDoMesonQA ) {
           
+         hESDEposEnegPsiPairDPhi[fiCut]->Fill(deltaPhi,psiPair);   
+         hESDEposEnegInvMassPt[fiCut]->Fill(vphoton->GetMass(),vphoton->Pt());
+                        
+        }
                                 
         if( ! fDoChicAnalysis ) {
        
@@ -1697,7 +1692,10 @@ void AliAnalysisTaskGammaConvDalitzV1::ProcessElectronCandidates(){
                       MassCutMax = ((AliDalitzElectronCuts*) fCutElectronArray->At(fiCut))->GetMassCutHighPt();
                     }
                     
-                    if( vphoton->GetMass() > MassCutMax ) {
+                    Bool_t DoMassMinCut = ((AliDalitzElectronCuts*)fCutElectronArray->At(fiCut))->DoMassMinCut();
+                    Double_t MassMinCut = ((AliDalitzElectronCuts*)fCutElectronArray->At(fiCut))->GetMassMinCut();
+                    
+                    if( vphoton->GetMass() > MassCutMax || (  DoMassMinCut && vphoton->GetMass() < MassMinCut )  ) {
                       
                       
                       delete vphoton;
@@ -1722,7 +1720,7 @@ void AliAnalysisTaskGammaConvDalitzV1::ProcessElectronCandidates(){
  
    //Computing mixing event
    
-   if(  fDoChicAnalysis ) {
+   if(  fDoMesonQA ) {
 
       for(UInt_t i = 0; i < lGoodElectronIndex.size(); i++){
          
@@ -1794,12 +1792,11 @@ void AliAnalysisTaskGammaConvDalitzV1::ProcessElectronCandidates(){
                virtualPhoton = 0x0;
           
          }
-   }   
+      }   
 
    }
 }
 
-
 //________________________________________________________________________
 void AliAnalysisTaskGammaConvDalitzV1::CalculatePi0DalitzCandidates(){
 

diff --git a/PWGGA/GammaConv/AliDalitzElectronCuts.cxx b/PWGGA/GammaConv/AliDalitzElectronCuts.cxx
index 19b2a8f681a4d64b86a4ac01b9e75f0b39329068..8a6cf46ee245902b1238345148b2d390a3132d94 100644 (file)
--- a/PWGGA/GammaConv/AliDalitzElectronCuts.cxx
+++ b/PWGGA/GammaConv/AliDalitzElectronCuts.cxx
@@ -115,9 +115,11 @@ AliDalitzElectronCuts::AliDalitzElectronCuts(const char *name,const char *title)
     fBKGMethod(0),
     fnumberOfRotationEventsForBG(0),
     fDoMassCut(kFALSE),
+    fDoMassMinCut(kFALSE),
     fMassCutLowPt(999.),
     fMassCutHighPt(999.),
     fMassCutPtMin(-100.0),
+    fMassMinCut(-999.),
     fDoWeights(kFALSE),
     fCutString(NULL),
     hCutIndex(NULL),
@@ -168,6 +170,8 @@ void AliDalitzElectronCuts::InitCutHistograms(TString name, Bool_t preCut,TStrin
 
     // Initialize Cut Histograms for QA (only initialized and filled if function is called)
 
+     TH1::AddDirectory(kFALSE);
+
      TString cutName = "";
     
      if( cutNumber==""){
@@ -332,7 +336,8 @@ void AliDalitzElectronCuts::InitCutHistograms(TString name, Bool_t preCut,TStrin
     }
     delete [] newBins;
 
-        
+    TH1::AddDirectory(kTRUE);        
+
     // Event Cuts and Info
 }
 
@@ -1764,63 +1769,88 @@ Bool_t AliDalitzElectronCuts::SetMassCut(Int_t massCut)
                         fMassCutPtMin  = -999.; //GeV
                         fMassCutLowPt  =  999.; //GeV/c^2
                         fMassCutHighPt =  999.; //GeV/c^2
+                        fMassMinCut = -999;
                         fDoMassCut = kFALSE;   
+                       fDoMassMinCut = kFALSE;
                         break;
     case 1:
                         //fMassCut = 0.135;             //GeV/c^2
                         fMassCutPtMin  = -999.; //GeV
                         fMassCutLowPt  = 0.135; //GeV/c^2
                         fMassCutHighPt = 0.135; //GeV/c^2
+                        fMassMinCut = -999;
                         fDoMassCut = kTRUE;
+                       fDoMassMinCut = kFALSE;
                         break; 
     case 2:
                         //fMassCut = 0.100;     //GeV/c^2
                         fMassCutPtMin  = -999.; //GeV
                         fMassCutLowPt  = 0.100; //GeV/c^2
                         fMassCutHighPt = 0.100; //GeV/c^2
+                        fMassMinCut = -999;
                         fDoMassCut = kTRUE;
+                       fDoMassMinCut = kFALSE;
                         break;
     case 3:
-                        //fMassCut = 0.075;     //GeV/c^2
+                        /*fMassCut = 0.075;     //GeV/c^2 Changed from Feb 25
                         fMassCutPtMin  = -999.; //GeV
                         fMassCutLowPt  = 0.075; //GeV/c^2
                         fMassCutHighPt = 0.075; //GeV/c^2
                         fDoMassCut = kTRUE;
+                       fDoMassMinCut = kFALSE;*/
+                       fMassCutPtMin  = 1.0;   //GeV
+                        fMassCutLowPt  = 0.015; //GeV/c^2
+                        fMassCutHighPt = 0.035; //GeV/c^2
+                        fMassMinCut    = 0.002;
+                        fDoMassCut = kTRUE;
+                       fDoMassMinCut = kTRUE;
                         break;
     case 4:
                         //fMassCut = 0.050;     //GeV/c^2
                         fMassCutPtMin  = -999.; //GeV
                         fMassCutLowPt  = 0.050; //GeV/c^2
                         fMassCutHighPt = 0.050; //GeV/c^2
+                        fMassMinCut = -999;
                         fDoMassCut = kTRUE;
+                       fDoMassMinCut = kFALSE;
                         break;
     case 5:
                         
                         fMassCutPtMin  = -999.; //GeV
                         fMassCutLowPt  = 0.035; //GeV/c^2
                         fMassCutHighPt = 0.035; //GeV/c^2
+                        fMassMinCut = -999;
                         fDoMassCut = kTRUE;
+                       fDoMassMinCut = kFALSE;
                         break;
     case 6:
                         fMassCutPtMin  = -999.; //GeV
                         fMassCutLowPt  = 0.015; //GeV/c^2
                         fMassCutHighPt = 0.015; //GeV/c^2
+                        fMassMinCut = -999;
                         fDoMassCut = kTRUE;
+                       fDoMassMinCut = kFALSE;
                         break;
     case 7:             fMassCutPtMin  = 1.0;   //GeV
                         fMassCutLowPt  = 0.015; //GeV/c^2
                         fMassCutHighPt = 0.035; //GeV/c^2
+                        fMassMinCut = -999;
                         fDoMassCut = kTRUE;
+                       fDoMassMinCut = kFALSE;
                         break;
     case 8:             fMassCutPtMin  = 1.0;   //GeV
                         fMassCutLowPt  = 0.015; //GeV/c^2
                         fMassCutHighPt = 0.050; //GeV/c^2
+                        fMassMinCut = -999;
                         fDoMassCut = kTRUE;
+                       fDoMassMinCut = kFALSE;
                        break;
     case 9:             fMassCutPtMin  = 1.0;   //GeV
                         fMassCutLowPt  = 0.025; //GeV/c^2
                         fMassCutHighPt = 0.035; //GeV/c^2
+                        fMassMinCut = -999;
                         fDoMassCut = kTRUE;
+                       fDoMassMinCut = kFALSE;
                        break;
     default:
                         cout<<"Warning: MassCut not defined "<<massCut<<endl;

diff --git a/PWGGA/GammaConv/AliDalitzElectronCuts.h b/PWGGA/GammaConv/AliDalitzElectronCuts.h
index 5c7ce1f2c5d55aecadc62a4799fbe333d42566ad..a1913b430bf1a0dba34dc54662aa5bf85f6741a5 100644 (file)
--- a/PWGGA/GammaConv/AliDalitzElectronCuts.h
+++ b/PWGGA/GammaConv/AliDalitzElectronCuts.h
@@ -156,9 +156,11 @@ class AliDalitzElectronCuts : public AliAnalysisCuts {
   Int_t    GetBKGMethod(){ return fBKGMethod; }
   Int_t    NumberOfRotationEvents(){return fnumberOfRotationEventsForBG;}
   Bool_t   DoMassCut(){return  fDoMassCut;}
+  Bool_t   DoMassMinCut(){return fDoMassMinCut;}
   Double_t GetMassCutLowPt(){return fMassCutLowPt;}
   Double_t GetMassCutHighPt(){return fMassCutHighPt;}
   Double_t GetPtMinMassCut(){return fMassCutPtMin;}
+  Double_t GetMassMinCut(){return fMassMinCut;}
   Bool_t   DoWeights(){return fDoWeights;}
   
 
@@ -212,9 +214,11 @@ class AliDalitzElectronCuts : public AliAnalysisCuts {
   Int_t    fBKGMethod;
   Int_t    fnumberOfRotationEventsForBG;
   Bool_t   fDoMassCut;
+  Bool_t   fDoMassMinCut;
   Double_t fMassCutLowPt;
   Double_t fMassCutHighPt;
   Double_t fMassCutPtMin;
+  Double_t fMassMinCut;
   Bool_t   fDoWeights;
 
 

diff --git a/PWGGA/GammaConv/macros/AddTask_GammaConvDalitzQAV1_pPb.C b/PWGGA/GammaConv/macros/AddTask_GammaConvDalitzQAV1_pPb.C
index 8d0fdc45d58a1163335969a7280bde4a10f24a37..746700ec0cbdfca658c58dffb716d85df39aecbc 100644 (file)
--- a/PWGGA/GammaConv/macros/AddTask_GammaConvDalitzQAV1_pPb.C
+++ b/PWGGA/GammaConv/macros/AddTask_GammaConvDalitzQAV1_pPb.C
@@ -188,17 +188,22 @@ if( trainConfig == 1 ) {  // No eta shift |Y| < 0.8
        
 }  else if( trainConfig == 2 ) {  // No eta shift |Y| < 0.8
   
-       ConvCutarray[0] = "8000012032093603007200000000"; ElecCutarray[0] = "9047540023910262371"; MesonCutarray[0] = "01033035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 +  |Y| < 0.6 and |Gamma_eta| < 0.65 and |e+_eta| < 0.65 and |e-_eta| < 0.65 
+       ConvCutarray[0] = "8000011032093603007200000000"; ElecCutarray[0] = "9047540023910262371"; MesonCutarray[0] = "01033035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 +  |Y| < 0.6 and |Gamma_eta| < 0.65 and |e+_eta| < 0.65 and |e-_eta| < 0.65 
        
 }  else if( trainConfig == 3 ) {  // No eta shift |Y| < 0.8
   
-       ConvCutarray[0] = "8000012042093603007200000000"; ElecCutarray[0] = "9047540023510262371"; MesonCutarray[0] = "01032035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 +  |Y| < 0.7 and |Gamma_eta| < 0.75 and |e+_eta| < 0.75 and |e-_eta| < 0.75
+       ConvCutarray[0] = "8000011042093603007200000000"; ElecCutarray[0] = "9047540023510262371"; MesonCutarray[0] = "01032035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 +  |Y| < 0.7 and |Gamma_eta| < 0.75 and |e+_eta| < 0.75 and |e-_eta| < 0.75
 
 }  else if( trainConfig == 4 ) {  // No eta shift  |Y| < 0.8
   
-       ConvCutarray[0] = "8000012012093603007200000000"; ElecCutarray[0] = "9047540023610262371"; MesonCutarray[0] = "01034035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 +  |Y| < 0.5 and |Gamma_eta| < 0.60 and |e+_eta| < 0.60 and |e-_eta| < 0.60  
+       ConvCutarray[0] = "8000011012093603007200000000"; ElecCutarray[0] = "9047540023610262371"; MesonCutarray[0] = "01034035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 +  |Y| < 0.5 and |Gamma_eta| < 0.60 and |e+_eta| < 0.60 and |e-_eta| < 0.60  
     
-}  
+} else if ( trainConfig == 5 ) {
+
+       ConvCutarray[0] = "8000011002093603007200000000"; ElecCutarray[0] = "9047540023310262331"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011
+
+
+} 
 
 
 
@@ -230,7 +235,7 @@ if( trainConfig == 1 ) {  // No eta shift |Y| < 0.8
             return 0;
       } else {
 
-         if (  ( trainConfig >= 1 && trainConfig <= 4 )   ){
+         if (  ( trainConfig >= 1 && trainConfig <= 5 )   ){
            
            if (doWeighting){
              if (generatorName.CompareTo("DPMJET")==0){

diff --git a/PWGGA/GammaConv/macros/AddTask_GammaConvV1_pPb.C b/PWGGA/GammaConv/macros/AddTask_GammaConvV1_pPb.C
index 79a698f4c627ee4752becc49c33df10328761dd6..38e0f87af804c0267247d5c6098dbca7b3d894db 100644 (file)
--- a/PWGGA/GammaConv/macros/AddTask_GammaConvV1_pPb.C
+++ b/PWGGA/GammaConv/macros/AddTask_GammaConvV1_pPb.C
@@ -102,14 +102,14 @@ void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of
    Bool_t doEtaShiftIndCuts = kFALSE;
    TString stringShift = "";
        if(trainConfig == 1){
-               cutarray[ 0] = "8000011003092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
-               cutarray[ 1] = "8000011004092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
-               cutarray[ 2] = "8000011001092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
+               cutarray[ 0] = "8000011032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
+               cutarray[ 1] = "8000011042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
+               cutarray[ 2] = "8000011012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
                cutarray[ 3] = "8000011002092370023220000000"; mesonCutArray[ 3] = "01621035009000";                       
        } else if (trainConfig == 2) {
-               cutarray[ 0] = "8000012003092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
-               cutarray[ 1] = "8000012004092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
-               cutarray[ 2] = "8000012001092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
+               cutarray[ 0] = "8000012032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
+               cutarray[ 1] = "8000012042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
+               cutarray[ 2] = "8000012012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
                cutarray[ 3] = "8000012002092370023220000000"; mesonCutArray[ 3] = "01621035009000";                       
        } else if (trainConfig == 3) {
                cutarray[ 0] = "8000011002092770023220000000"; mesonCutArray[ 0] = "01621035009000";                       
@@ -212,14 +212,14 @@ void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of
                cutarray[ 2] = "8000012002093172003290000000"; mesonCutArray[ 2] = "01621035000000";  //old standard eta=0.9 y=0.8 
                cutarray[ 3] = "8000012002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
        } else if (trainConfig == 23){
-               cutarray[ 0] = "8020011003092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
-               cutarray[ 1] = "8020011004092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
-               cutarray[ 2] = "8020011001092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
+               cutarray[ 0] = "8020011032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
+               cutarray[ 1] = "8020011042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
+               cutarray[ 2] = "8020011012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
                cutarray[ 3] = "8020011002092370023220000000"; mesonCutArray[ 3] = "01621035009000";                       
        } else if (trainConfig == 24) {
-               cutarray[ 0] = "8020012003092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
-               cutarray[ 1] = "8020012004092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
-               cutarray[ 2] = "8020012001092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
+               cutarray[ 0] = "8020012032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
+               cutarray[ 1] = "8020012042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
+               cutarray[ 2] = "8020012012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
                cutarray[ 3] = "8020012002092370023220000000"; mesonCutArray[ 3] = "01621035009000";                       
        } else if (trainConfig == 25) {
                cutarray[ 0] = "8020011002092770023220000000"; mesonCutArray[ 0] = "01621035009000";                       
@@ -322,14 +322,14 @@ void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of
                cutarray[ 2] = "8020012002093172003290000000"; mesonCutArray[ 2] = "01621035000000";  //old standard eta=0.9 y=0.8 
                cutarray[ 3] = "8020012002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
        } else if (trainConfig == 45){
-               cutarray[ 0] = "8240011003092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
-               cutarray[ 1] = "8240011004092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
-               cutarray[ 2] = "8240011001092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
+               cutarray[ 0] = "8240011032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
+               cutarray[ 1] = "8240011042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
+               cutarray[ 2] = "8240011012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
                cutarray[ 3] = "8240011002092370023220000000"; mesonCutArray[ 3] = "01621035009000";                       
        } else if (trainConfig == 46) {
-               cutarray[ 0] = "8240012003092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
-               cutarray[ 1] = "8240012004092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
-               cutarray[ 2] = "8240012001092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
+               cutarray[ 0] = "8240012032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
+               cutarray[ 1] = "8240012042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
+               cutarray[ 2] = "8240012012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
                cutarray[ 3] = "8240012002092370023220000000"; mesonCutArray[ 3] = "01621035009000";                       
        } else if (trainConfig == 47) {
                cutarray[ 0] = "8240011002092770023220000000"; mesonCutArray[ 0] = "01621035009000";                       
@@ -432,14 +432,14 @@ void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of
                cutarray[ 2] = "8240012002093172003290000000"; mesonCutArray[ 2] = "01621035000000";  //old standard eta=0.9 y=0.8 
                cutarray[ 3] = "8240012002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
        } else if (trainConfig == 67){
-               cutarray[ 0] = "8460011003092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
-               cutarray[ 1] = "8460011004092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
-               cutarray[ 2] = "8460011001092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
+               cutarray[ 0] = "8460011032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
+               cutarray[ 1] = "8460011042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
+               cutarray[ 2] = "8460011012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
                cutarray[ 3] = "8460011002092370023220000000"; mesonCutArray[ 3] = "01621035009000";                       
        } else if (trainConfig == 68) {
-               cutarray[ 0] = "8460012003092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
-               cutarray[ 1] = "8460012004092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
-               cutarray[ 2] = "8460012001092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
+               cutarray[ 0] = "8460012032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
+               cutarray[ 1] = "8460012042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
+               cutarray[ 2] = "8460012012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
                cutarray[ 3] = "8460012002092370023220000000"; mesonCutArray[ 3] = "01621035009000";                       
        } else if (trainConfig == 69) {
                cutarray[ 0] = "8460011002092770023220000000"; mesonCutArray[ 0] = "01621035009000";                       
@@ -542,14 +542,14 @@ void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of
                cutarray[ 2] = "8460012002093172003290000000"; mesonCutArray[ 2] = "01621035000000";  //old standard eta=0.9 y=0.8 
                cutarray[ 3] = "8460012002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;   
        } else if (trainConfig == 89){
-               cutarray[ 0] = "8680011003092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
-               cutarray[ 1] = "8680011004092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
-               cutarray[ 2] = "8680011001092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
+               cutarray[ 0] = "8680011032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
+               cutarray[ 1] = "8680011042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
+               cutarray[ 2] = "8680011012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
                cutarray[ 3] = "8680011002092370023220000000"; mesonCutArray[ 3] = "01621035009000";                       
        } else if (trainConfig == 90) {
-               cutarray[ 0] = "8680012003092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
-               cutarray[ 1] = "8680012004092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
-               cutarray[ 2] = "8680012001092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
+               cutarray[ 0] = "8680012032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
+               cutarray[ 1] = "8680012042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
+               cutarray[ 2] = "8680012012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
                cutarray[ 3] = "8680012002092370023220000000"; mesonCutArray[ 3] = "01621035009000";                       
        } else if (trainConfig == 91) {
                cutarray[ 0] = "8680011002092770023220000000"; mesonCutArray[ 0] = "01621035009000";                       
@@ -652,14 +652,14 @@ void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of
                cutarray[ 2] = "8680012002093172003290000000"; mesonCutArray[ 2] = "01621035000000";  //old standard eta=0.9 y=0.8 
                cutarray[ 3] = "8680012002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;   
        } else if (trainConfig == 111){
-               cutarray[ 0] = "8600011003092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
-               cutarray[ 1] = "8600011004092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
-               cutarray[ 2] = "8600011001092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
+               cutarray[ 0] = "8600011032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
+               cutarray[ 1] = "8600011042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
+               cutarray[ 2] = "8600011012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
                cutarray[ 3] = "8600011002092370023220000000"; mesonCutArray[ 3] = "01621035009000";                       
        } else if (trainConfig == 112) {
-               cutarray[ 0] = "8600012003092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
-               cutarray[ 1] = "8600012004092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
-               cutarray[ 2] = "8600012001092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
+               cutarray[ 0] = "8600012032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
+               cutarray[ 1] = "8600012042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
+               cutarray[ 2] = "8600012012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
                cutarray[ 3] = "8600012002092370023220000000"; mesonCutArray[ 3] = "01621035009000";                       
        } else if (trainConfig == 113) {
                cutarray[ 0] = "8600011002092770023220000000"; mesonCutArray[ 0] = "01621035009000";                       
@@ -762,14 +762,14 @@ void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of
                cutarray[ 2] = "8600012002093172003290000000"; mesonCutArray[ 2] = "01621035000000";  //old standard eta=0.9 y=0.8 
                cutarray[ 3] = "8600012002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;   
        } else if (trainConfig == 133) {
-               cutarray[ 0] = "8000011003092170008260400000"; mesonCutArray[ 0] = "01623035000000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
-               cutarray[ 1] = "8000011004092170008260400000"; mesonCutArray[ 1] = "01622035000000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
-               cutarray[ 2] = "8000011001092170008260400000"; mesonCutArray[ 2] = "01624035000000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
+               cutarray[ 0] = "8000011032092170008260400000"; mesonCutArray[ 0] = "01623035000000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
+               cutarray[ 1] = "8000011042092170008260400000"; mesonCutArray[ 1] = "01622035000000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
+               cutarray[ 2] = "8000011012092170008260400000"; mesonCutArray[ 2] = "01624035000000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
                cutarray[ 3] = "8000011002092170003260000000"; mesonCutArray[ 3] = "01621035000000"; //tighten Psi pair and chi2 in 2D
        } else if (trainConfig == 134) {
-               cutarray[ 0] = "8000012003092170008260400000"; mesonCutArray[ 0] = "01623035000000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
-               cutarray[ 1] = "8000012004092170008260400000"; mesonCutArray[ 1] = "01622035000000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
-               cutarray[ 2] = "8000012001092170008260400000"; mesonCutArray[ 2] = "01624035000000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
+               cutarray[ 0] = "8000012032092170008260400000"; mesonCutArray[ 0] = "01623035000000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
+               cutarray[ 1] = "8000012042092170008260400000"; mesonCutArray[ 1] = "01622035000000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
+               cutarray[ 2] = "8000012012092170008260400000"; mesonCutArray[ 2] = "01624035000000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
                cutarray[ 3] = "8000012002092170003260000000"; mesonCutArray[ 3] = "01621035000000"; //tighten Psi pair and chi2 in 2D
        } else if (trainConfig == 135) {          
                cutarray[ 0] = "8000011002092170008220000000"; mesonCutArray[ 0] = "01621035000000"; //tighten psi pair and qt in 2D
@@ -862,14 +862,14 @@ void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of
                cutarray[ 2] = "8000012002093172003290000000"; mesonCutArray[ 2] = "01621035000000";  //old standard eta=0.9 y=0.8 
                cutarray[ 3] = "8000012002092170008260400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
        } else if (trainConfig == 153) {
-               cutarray[ 0] = "8020011003092170008260400000"; mesonCutArray[ 0] = "01623035000000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
-               cutarray[ 1] = "8020011004092170008260400000"; mesonCutArray[ 1] = "01622035000000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
-               cutarray[ 2] = "8020011001092170008260400000"; mesonCutArray[ 2] = "01624035000000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
+               cutarray[ 0] = "8020011032092170008260400000"; mesonCutArray[ 0] = "01623035000000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
+               cutarray[ 1] = "8020011042092170008260400000"; mesonCutArray[ 1] = "01622035000000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
+               cutarray[ 2] = "8020011012092170008260400000"; mesonCutArray[ 2] = "01624035000000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
                cutarray[ 3] = "8020011002092170003260000000"; mesonCutArray[ 3] = "01621035000000"; //tighten Psi pair and chi2 in 2D  
        } else if (trainConfig == 154) {
-               cutarray[ 0] = "8020012003092170008260400000"; mesonCutArray[ 0] = "01623035000000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
-               cutarray[ 1] = "8020012004092170008260400000"; mesonCutArray[ 1] = "01622035000000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
-               cutarray[ 2] = "8020012001092170008260400000"; mesonCutArray[ 2] = "01624035000000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
+               cutarray[ 0] = "8020012032092170008260400000"; mesonCutArray[ 0] = "01623035000000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
+               cutarray[ 1] = "8020012042092170008260400000"; mesonCutArray[ 1] = "01622035000000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
+               cutarray[ 2] = "8020012012092170008260400000"; mesonCutArray[ 2] = "01624035000000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
                cutarray[ 3] = "8020012002092170003260000000"; mesonCutArray[ 3] = "01621035000000"; //tighten Psi pair and chi2 in 2D
        } else if (trainConfig == 155) {          
                cutarray[ 0] = "8020011002092170008220000000"; mesonCutArray[ 0] = "01621035000000"; //tighten psi pair and qt in 2D