corrected calo binning / changed track matching cuts
authordmuhlhei <d.muehlheim@cern.ch>
Fri, 19 Dec 2014 13:42:26 +0000 (14:42 +0100)
committerdmuhlhei <d.muehlheim@cern.ch>
Tue, 6 Jan 2015 15:57:20 +0000 (16:57 +0100)
PWGGA/GammaConv/AliCaloPhotonCuts.cxx
PWGGA/GammaConv/AliCaloPhotonCuts.h
PWGGA/GammaConv/macros/AddTask_GammaConvCalo_PbPb.C
PWGGA/GammaConv/macros/AddTask_GammaConvCalo_pPb.C
PWGGA/GammaConv/macros/AddTask_GammaConvCalo_pp.C

index 327fb92..360bcd3 100644 (file)
@@ -146,12 +146,12 @@ AliCaloPhotonCuts::AliCaloPhotonCuts(const char *name,const char *title) :
     fHistClusterdEtadPhiNegTracksBeforeQA(NULL),
     fHistClusterdEtadPhiPosTracksAfterQA(NULL),
     fHistClusterdEtadPhiNegTracksAfterQA(NULL),
-    fHistClusterdEtadPhiPosTracks_P1_BeforeQA(NULL),
-    fHistClusterdEtadPhiPosTracks_P14_BeforeQA(NULL),
-    fHistClusterdEtadPhiPosTracks_P4_BeforeQA(NULL),
-    fHistClusterdEtadPhiNegTracks_P1_BeforeQA(NULL),
-    fHistClusterdEtadPhiNegTracks_P14_BeforeQA(NULL),
-    fHistClusterdEtadPhiNegTracks_P4_BeforeQA(NULL),
+    fHistClusterdEtadPhiPosTracksP_00_05BeforeQA(NULL),
+    fHistClusterdEtadPhiPosTracksP_05_15BeforeQA(NULL),
+    fHistClusterdEtadPhiPosTracksP_15_99BeforeQA(NULL),
+    fHistClusterdEtadPhiNegTracksP_00_05BeforeQA(NULL),
+    fHistClusterdEtadPhiNegTracksP_05_15BeforeQA(NULL),
+    fHistClusterdEtadPhiNegTracksP_15_99BeforeQA(NULL),
     fHistClusterdEtadPtBeforeQA(NULL),
     fHistClusterdPhidPtBeforeQA(NULL),
     fHistClusterM20Pt_dPhiBeforeQA(NULL),
@@ -236,12 +236,12 @@ AliCaloPhotonCuts::AliCaloPhotonCuts(const AliCaloPhotonCuts &ref) :
     fHistClusterdEtadPhiNegTracksBeforeQA(NULL),
     fHistClusterdEtadPhiPosTracksAfterQA(NULL),
     fHistClusterdEtadPhiNegTracksAfterQA(NULL),
-    fHistClusterdEtadPhiPosTracks_P1_BeforeQA(NULL),
-    fHistClusterdEtadPhiPosTracks_P14_BeforeQA(NULL),
-    fHistClusterdEtadPhiPosTracks_P4_BeforeQA(NULL),
-    fHistClusterdEtadPhiNegTracks_P1_BeforeQA(NULL),
-    fHistClusterdEtadPhiNegTracks_P14_BeforeQA(NULL),
-    fHistClusterdEtadPhiNegTracks_P4_BeforeQA(NULL),
+    fHistClusterdEtadPhiPosTracksP_00_05BeforeQA(NULL),
+    fHistClusterdEtadPhiPosTracksP_05_15BeforeQA(NULL),
+    fHistClusterdEtadPhiPosTracksP_15_99BeforeQA(NULL),
+    fHistClusterdEtadPhiNegTracksP_00_05BeforeQA(NULL),
+    fHistClusterdEtadPhiNegTracksP_05_15BeforeQA(NULL),
+    fHistClusterdEtadPhiNegTracksP_15_99BeforeQA(NULL),
     fHistClusterdEtadPtBeforeQA(NULL),
     fHistClusterdPhidPtBeforeQA(NULL),
     fHistClusterM20Pt_dPhiBeforeQA(NULL),
@@ -320,12 +320,42 @@ void AliCaloPhotonCuts::InitCutHistograms(TString name){
        fHistograms->Add(fHistClusterIdentificationCuts);
 
        // Acceptance related histogramms
-    fHistClusterEtavsPhiBeforeAcc=new TH2F(Form("EtaPhi_beforeAcceptance %s",GetCutNumber().Data()),"EtaPhi_beforeAcceptance",430,0,2*TMath::Pi(),96,-0.7,0.7); // 462->430, 110->96
-       fHistograms->Add(fHistClusterEtavsPhiBeforeAcc);
-    fHistClusterEtavsPhiAfterAcc=new TH2F(Form("EtaPhi_afterAcceptance %s",GetCutNumber().Data()),"EtaPhi_afterAcceptance",430,0,2*TMath::Pi(),96,-0.7,0.7); // 462->430, 110->96
-       fHistograms->Add(fHistClusterEtavsPhiAfterAcc);
-    fHistClusterEtavsPhiAfterQA=new TH2F(Form("EtaPhi_afterClusterQA %s",GetCutNumber().Data()),"EtaPhi_afterClusterQA",430,0,2*TMath::Pi(),96,-0.7,0.7); // 462->430, 110->96
-       fHistograms->Add(fHistClusterEtavsPhiAfterQA);
+    const Int_t nEmcalEtaBins = 96;
+    const Int_t nEmcalPhiBins = 124;
+    Float_t EmcalEtaBins[nEmcalEtaBins+1] = {-0.66687,-0.653,-0.63913,-0.62526,-0.61139,-0.59752,-0.58365,-0.56978,-0.55591,-0.54204,-0.52817,-0.5143,-0.50043,-0.48656,-0.47269,-0.45882,-0.44495,-0.43108,-0.41721,-0.40334,-0.38947,-0.3756,-0.36173,-0.34786,-0.33399,-0.32012,-0.30625,-0.29238,-0.27851,-0.26464,-0.25077,-0.2369,-0.22303,-0.20916,-0.19529,-0.18142,-0.16755,-0.15368,-0.13981,-0.12594,-0.11207,-0.0982,-0.08433,-0.07046,-0.05659,-0.04272,-0.02885,-0.01498,-0.00111,0.01276,0.02663,0.0405,0.05437,0.06824,0.08211,0.09598,0.10985,0.12372,0.13759,0.15146,0.16533,0.1792,0.19307,0.20694,0.22081,0.23468,0.24855,0.26242,0.27629,0.29016,0.30403,0.3179,0.33177,0.34564,0.35951,0.37338,0.38725,0.40112,0.41499,0.42886,0.44273,0.4566,0.47047,0.48434,0.49821,0.51208,0.52595,0.53982,0.55369,0.56756,0.58143,0.5953,0.60917,0.62304,0.63691,0.65078,0.66465};
+    Float_t EmcalPhiBins[nEmcalPhiBins+1] = {1.408,1.4215,1.435,1.4485,1.462,1.4755,1.489,1.5025,1.516,1.5295,1.543,1.5565,1.57,1.5835,1.597,1.6105,1.624,1.6375,1.651,1.6645,1.678,1.6915,1.705,1.7185,1.732,1.758,1.7715,1.785,1.7985,1.812,1.8255,1.839,1.8525,1.866,1.8795,1.893,1.9065,1.92,1.9335,1.947,1.9605,1.974,1.9875,2.001,2.0145,2.028,2.0415,2.055,2.0685,2.082,2.108,2.1215,2.135,2.1485,2.162,2.1755,2.189,2.2025,2.216,2.2295,2.243,2.2565,2.27,2.2835,2.297,2.3105,2.324,2.3375,2.351,2.3645,2.378,2.3915,2.405,2.4185,2.432,2.456,2.4695,2.483,2.4965,2.51,2.5235,2.537,2.5505,2.564,2.5775,2.591,2.6045,2.618,2.6315,2.645,2.6585,2.672,2.6855,2.699,2.7125,2.726,2.7395,2.753,2.7665,2.78,2.804,2.8175,2.831,2.8445,2.858,2.8715,2.885,2.8985,2.912,2.9255,2.939,2.9525,2.966,2.9795,2.993,3.0065,3.02,3.0335,3.047,3.0605,3.074,3.0875,3.101,3.1145,3.128};
+
+    const Int_t nPhosEtaBins = 56;
+    const Int_t nPhosPhiBins = 192;
+    const Float_t PhosEtaRange[2] = {-0.16, 0.16};
+    const Float_t PhosPhiRange[2] = {4.5, 5.6};
+
+    if( GetClusterType() == 1 ){ //EMCAL
+        fHistClusterEtavsPhiBeforeAcc=new TH2F(Form("EtaPhi_beforeAcceptance %s",GetCutNumber().Data()),"EtaPhi_beforeAcceptance",nEmcalPhiBins,EmcalPhiBins,nEmcalEtaBins,EmcalEtaBins);
+        fHistograms->Add(fHistClusterEtavsPhiBeforeAcc);
+        fHistClusterEtavsPhiAfterAcc=new TH2F(Form("EtaPhi_afterAcceptance %s",GetCutNumber().Data()),"EtaPhi_afterAcceptance",nEmcalPhiBins,EmcalPhiBins,nEmcalEtaBins,EmcalEtaBins);
+        fHistograms->Add(fHistClusterEtavsPhiAfterAcc);
+        fHistClusterEtavsPhiAfterQA=new TH2F(Form("EtaPhi_afterClusterQA %s",GetCutNumber().Data()),"EtaPhi_afterClusterQA",nEmcalPhiBins,EmcalPhiBins,nEmcalEtaBins,EmcalEtaBins);
+        fHistograms->Add(fHistClusterEtavsPhiAfterQA);
+    }
+    else if( GetClusterType() == 2 ){ //PHOS
+        fHistClusterEtavsPhiBeforeAcc=new TH2F(Form("EtaPhi_beforeAcceptance %s",GetCutNumber().Data()),"EtaPhi_beforeAcceptance",nPhosPhiBins,PhosPhiRange[0],PhosPhiRange[1],nPhosEtaBins,PhosEtaRange[0],PhosEtaRange[1]);
+        fHistograms->Add(fHistClusterEtavsPhiBeforeAcc);
+        fHistClusterEtavsPhiAfterAcc=new TH2F(Form("EtaPhi_afterAcceptance %s",GetCutNumber().Data()),"EtaPhi_afterAcceptance",nPhosPhiBins,PhosPhiRange[0],PhosPhiRange[1],nPhosEtaBins,PhosEtaRange[0],PhosEtaRange[1]);
+        fHistograms->Add(fHistClusterEtavsPhiAfterAcc);
+        fHistClusterEtavsPhiAfterQA=new TH2F(Form("EtaPhi_afterClusterQA %s",GetCutNumber().Data()),"EtaPhi_afterClusterQA",nPhosPhiBins,PhosPhiRange[0],PhosPhiRange[1],nPhosEtaBins,PhosEtaRange[0],PhosEtaRange[1]);
+        fHistograms->Add(fHistClusterEtavsPhiAfterQA);
+    }
+    else if( GetClusterType() == 0 ){ //all
+        fHistClusterEtavsPhiBeforeAcc=new TH2F(Form("EtaPhi_beforeAcceptance %s",GetCutNumber().Data()),"EtaPhi_beforeAcceptance",462,0,2*TMath::Pi(),110,-0.7,0.7);
+        fHistograms->Add(fHistClusterEtavsPhiBeforeAcc);
+        fHistClusterEtavsPhiAfterAcc=new TH2F(Form("EtaPhi_afterAcceptance %s",GetCutNumber().Data()),"EtaPhi_afterAcceptance",462,0,2*TMath::Pi(),110,-0.7,0.7);
+        fHistograms->Add(fHistClusterEtavsPhiAfterAcc);
+        fHistClusterEtavsPhiAfterQA=new TH2F(Form("EtaPhi_afterClusterQA %s",GetCutNumber().Data()),"EtaPhi_afterClusterQA",462,0,2*TMath::Pi(),110,-0.7,0.7);
+        fHistograms->Add(fHistClusterEtavsPhiAfterQA);
+    }
+    else{AliError(Form("Cluster Type is not EMCAL nor PHOS: %i",GetClusterType()));}
+
     //fHistDistanceToBadChannelBeforeAcc = new TH1F(Form("DistanceToBadChannel_beforeAcceptance %s",GetCutNumber().Data()),"DistanceToBadChannel_beforeAcceptance",200,0,40);
     //fHistograms->Add(fHistDistanceToBadChannelBeforeAcc);
     //fHistDistanceToBadChannelAfterAcc = new TH1F(Form("DistanceToBadChannel_afterAcceptance %s",GetCutNumber().Data()),"DistanceToBadChannel_afterAcceptance",200,0,40);
@@ -334,10 +364,11 @@ void AliCaloPhotonCuts::InitCutHistograms(TString name){
        // Cluster quality related histograms
     Double_t timeMin = -2e-6;
     Double_t timeMax = 8e-6;
-       if (GetClusterType() == 1){
+    if( GetClusterType() == 1 ){
                timeMin = -2e-7;
                timeMax = 12e-7;
-       }       
+    }
+
     fHistClusterTimevsEBeforeQA=new TH2F(Form("ClusterTimeVsE_beforeClusterQA %s",GetCutNumber().Data()),"ClusterTimeVsE_beforeClusterQA",800,timeMin,timeMax,100,0,40);
        fHistograms->Add(fHistClusterTimevsEBeforeQA);
     fHistClusterTimevsEAfterQA=new TH2F(Form("ClusterTimeVsE_afterClusterQA %s",GetCutNumber().Data()),"ClusterTimeVsE_afterClusterQA",800,timeMin,timeMax,100,0,40);
@@ -375,13 +406,18 @@ void AliCaloPhotonCuts::InitCutHistograms(TString name){
 
     //TrackMatching histograms
     if(fUseDistTrackToCluster) {
+        const Int_t nEtaBins = 300;
+        const Int_t nPhiBins = 300;
+        const Float_t EtaRange[2] = {-0.3, 0.3};
+        const Float_t PhiRange[2] = {-0.3, 0.3};
+
         fHistClusterRBeforeQA = new TH1F(Form("R_Cluster_beforeClusterQA %s",GetCutNumber().Data()),"R of cluster",200,400,500);
         fHistograms->Add(fHistClusterRBeforeQA);
         fHistClusterRAfterQA = new TH1F(Form("R_Cluster_afterClusterQA %s",GetCutNumber().Data()),"R of cluster_matched",200,400,500);
         fHistograms->Add(fHistClusterRAfterQA);
-        fHistClusterdEtadPhiBeforeQA=new TH2F(Form("dEtaVsdPhi_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_beforeClusterQA",240,-0.3,0.3,240,-0.3,0.3);
+        fHistClusterdEtadPhiBeforeQA=new TH2F(Form("dEtaVsdPhi_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_beforeClusterQA",nEtaBins,EtaRange[0],EtaRange[1],nPhiBins,PhiRange[0],PhiRange[1]);
         fHistograms->Add(fHistClusterdEtadPhiBeforeQA);
-        fHistClusterdEtadPhiAfterQA=new TH2F(Form("dEtaVsdPhi_afterClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_afterClusterQA",240,-0.3,0.3,240,-0.3,0.3);
+        fHistClusterdEtadPhiAfterQA=new TH2F(Form("dEtaVsdPhi_afterClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_afterClusterQA",nEtaBins,EtaRange[0],EtaRange[1],nPhiBins,PhiRange[0],PhiRange[1]);
         fHistograms->Add(fHistClusterdEtadPhiAfterQA);
         fHistDistanceTrackToClusterBeforeQA = new TH1F(Form("DistanceToTrack_beforeClusterQA %s",GetCutNumber().Data()),"DistanceToTrack_beforeClusterQA",200,0,2);
         fHistograms->Add(fHistDistanceTrackToClusterBeforeQA);
@@ -389,29 +425,29 @@ void AliCaloPhotonCuts::InitCutHistograms(TString name){
         fHistograms->Add(fHistDistanceTrackToClusterAfterQA);
        
         if(fExtendedMatching){
-            fHistClusterdEtadPhiPosTracksBeforeQA = new TH2F(Form("dEtaVsdPhi_posTracks_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_posTracks_beforeClusterQA",240,-0.3,0.3,240,-0.3,0.3);
+            fHistClusterdEtadPhiPosTracksBeforeQA = new TH2F(Form("dEtaVsdPhi_posTracks_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_posTracks_beforeClusterQA",nEtaBins,EtaRange[0],EtaRange[1],nPhiBins,PhiRange[0],PhiRange[1]);
             fHistograms->Add(fHistClusterdEtadPhiPosTracksBeforeQA);
-            fHistClusterdEtadPhiNegTracksBeforeQA = new TH2F(Form("dEtaVsdPhi_negTracks_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_negTracks_beforeClusterQA",240,-0.3,0.3,240,-0.3,0.3);
+            fHistClusterdEtadPhiNegTracksBeforeQA = new TH2F(Form("dEtaVsdPhi_negTracks_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_negTracks_beforeClusterQA",nEtaBins,EtaRange[0],EtaRange[1],nPhiBins,PhiRange[0],PhiRange[1]);
             fHistograms->Add(fHistClusterdEtadPhiNegTracksBeforeQA);
-            fHistClusterdEtadPhiPosTracksAfterQA = new TH2F(Form("dEtaVsdPhi_posTracks_afterClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_posTracks_afterClusterQA",240,-0.3,0.3,240,-0.3,0.3);
+            fHistClusterdEtadPhiPosTracksAfterQA = new TH2F(Form("dEtaVsdPhi_posTracks_afterClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_posTracks_afterClusterQA",nEtaBins,EtaRange[0],EtaRange[1],nPhiBins,PhiRange[0],PhiRange[1]);
             fHistograms->Add(fHistClusterdEtadPhiPosTracksAfterQA);
-            fHistClusterdEtadPhiNegTracksAfterQA = new TH2F(Form("dEtaVsdPhi_negTracks_afterClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_negTracks_afterClusterQA",240,-0.3,0.3,240,-0.3,0.3);
+            fHistClusterdEtadPhiNegTracksAfterQA = new TH2F(Form("dEtaVsdPhi_negTracks_afterClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_negTracks_afterClusterQA",nEtaBins,EtaRange[0],EtaRange[1],nPhiBins,PhiRange[0],PhiRange[1]);
             fHistograms->Add(fHistClusterdEtadPhiNegTracksAfterQA);
-            fHistClusterdEtadPhiPosTracks_P1_BeforeQA = new TH2F(Form("dEtaVsdPhi_posTracks_P<1_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_posTracks_P<1_beforeClusterQA",240,-0.3,0.3,240,-0.3,0.3);
-            fHistograms->Add(fHistClusterdEtadPhiPosTracks_P1_BeforeQA);
-            fHistClusterdEtadPhiPosTracks_P14_BeforeQA = new TH2F(Form("dEtaVsdPhi_posTracks_1<P<4_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_posTracks_1<P<4_beforeClusterQA",240,-0.3,0.3,240,-0.3,0.3);
-            fHistograms->Add(fHistClusterdEtadPhiPosTracks_P14_BeforeQA);
-            fHistClusterdEtadPhiPosTracks_P4_BeforeQA = new TH2F(Form("dEtaVsdPhi_posTracks_P>4_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_posTracks_P>4_beforeClusterQA",240,-0.3,0.3,240,-0.3,0.3);
-            fHistograms->Add(fHistClusterdEtadPhiPosTracks_P4_BeforeQA);
-            fHistClusterdEtadPhiNegTracks_P1_BeforeQA= new TH2F(Form("dEtaVsdPhi_negTrack_P<1s_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_negTrack_P<1s_beforeClusterQA",240,-0.3,0.3,240,-0.3,0.3);
-            fHistograms->Add(fHistClusterdEtadPhiNegTracks_P1_BeforeQA);
-            fHistClusterdEtadPhiNegTracks_P14_BeforeQA = new TH2F(Form("dEtaVsdPhi_negTracks_1<P<4_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_negTracks_1<P<4_beforeClusterQA",240,-0.3,0.3,240,-0.3,0.3);
-            fHistograms->Add(fHistClusterdEtadPhiNegTracks_P14_BeforeQA);
-            fHistClusterdEtadPhiNegTracks_P4_BeforeQA = new TH2F(Form("dEtaVsdPhi_negTracks_P>4_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_negTracks_P>4_beforeClusterQA",240,-0.3,0.3,240,-0.3,0.3);
-            fHistograms->Add(fHistClusterdEtadPhiNegTracks_P4_BeforeQA);
-            fHistClusterdEtadPtBeforeQA = new TH2F(Form("dEtaVsPt_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsPt_beforeClusterQA",240,-0.3,0.3,250,0,25);
+            fHistClusterdEtadPhiPosTracksP_00_05BeforeQA = new TH2F(Form("dEtaVsdPhi_posTracks_P<0.5_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_posTracks_P<0.5_beforeClusterQA",nEtaBins,EtaRange[0],EtaRange[1],nPhiBins,PhiRange[0],PhiRange[1]);
+            fHistograms->Add(fHistClusterdEtadPhiPosTracksP_00_05BeforeQA);
+            fHistClusterdEtadPhiPosTracksP_05_15BeforeQA = new TH2F(Form("dEtaVsdPhi_posTracks_0.5<P<1.5_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_posTracks_0.5<P<1.5_beforeClusterQA",nEtaBins,EtaRange[0],EtaRange[1],nPhiBins,PhiRange[0],PhiRange[1]);
+            fHistograms->Add(fHistClusterdEtadPhiPosTracksP_05_15BeforeQA);
+            fHistClusterdEtadPhiPosTracksP_15_99BeforeQA = new TH2F(Form("dEtaVsdPhi_posTracks_P>1.5_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_posTracks_P>1.5_beforeClusterQA",nEtaBins,EtaRange[0],EtaRange[1],nPhiBins,PhiRange[0],PhiRange[1]);
+            fHistograms->Add(fHistClusterdEtadPhiPosTracksP_15_99BeforeQA);
+            fHistClusterdEtadPhiNegTracksP_00_05BeforeQA= new TH2F(Form("dEtaVsdPhi_negTrack_P<0.5_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_negTrack_P<0.5_beforeClusterQA",nEtaBins,EtaRange[0],EtaRange[1],nPhiBins,PhiRange[0],PhiRange[1]);
+            fHistograms->Add(fHistClusterdEtadPhiNegTracksP_00_05BeforeQA);
+            fHistClusterdEtadPhiNegTracksP_05_15BeforeQA = new TH2F(Form("dEtaVsdPhi_negTracks_0.5<P<1.5_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_negTracks_0.5<P<1.5_beforeClusterQA",nEtaBins,EtaRange[0],EtaRange[1],nPhiBins,PhiRange[0],PhiRange[1]);
+            fHistograms->Add(fHistClusterdEtadPhiNegTracksP_05_15BeforeQA);
+            fHistClusterdEtadPhiNegTracksP_15_99BeforeQA = new TH2F(Form("dEtaVsdPhi_negTracks_P>1.5_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsdPhi_negTracks_P>1.5_beforeClusterQA",nEtaBins,EtaRange[0],EtaRange[1],nPhiBins,PhiRange[0],PhiRange[1]);
+            fHistograms->Add(fHistClusterdEtadPhiNegTracksP_15_99BeforeQA);
+            fHistClusterdEtadPtBeforeQA = new TH2F(Form("dEtaVsPt_beforeClusterQA %s",GetCutNumber().Data()),"dEtaVsPt_beforeClusterQA",nEtaBins,EtaRange[0],EtaRange[1],250,0,25);
             fHistograms->Add(fHistClusterdEtadPtBeforeQA);
-            fHistClusterdPhidPtBeforeQA = new TH2F(Form("dPhiVsPt_beforeClusterQA %s",GetCutNumber().Data()),"dPhiVsPt_beforeClusterQA",480,-0.6,0.6,250,0,25);
+            fHistClusterdPhidPtBeforeQA = new TH2F(Form("dPhiVsPt_beforeClusterQA %s",GetCutNumber().Data()),"dPhiVsPt_beforeClusterQA",2*nPhiBins,2*PhiRange[0],2*PhiRange[1],250,0,25);
             fHistograms->Add(fHistClusterdPhidPtBeforeQA);
             fHistClusterM20Pt_dPhiBeforeQA = new TH2F(Form("M20VsPt_dPhi_beforeClusterQA %s",GetCutNumber().Data()),"M20VsPt_dPhi_beforeClusterQA",200,0,2.5,250,0,25);
             fHistograms->Add(fHistClusterM20Pt_dPhiBeforeQA);
@@ -882,15 +918,15 @@ Bool_t AliCaloPhotonCuts::MatchConvPhotonToCluster(AliAODConversionPhoton* convP
             if(fExtendedMatching){
                 if(inTrack->Charge() > 0) {
                     fHistClusterdEtadPhiPosTracksBeforeQA->Fill(dEta, dPhi);
-                    if(inTrack->P() < 1) fHistClusterdEtadPhiPosTracks_P1_BeforeQA->Fill(dEta, dPhi);
-                    else if(inTrack->P() < 4) fHistClusterdEtadPhiPosTracks_P14_BeforeQA->Fill(dEta, dPhi);
-                    else fHistClusterdEtadPhiPosTracks_P4_BeforeQA->Fill(dEta, dPhi);
+                    if(inTrack->P() < 0.5) fHistClusterdEtadPhiPosTracksP_00_05BeforeQA->Fill(dEta, dPhi);
+                    else if(inTrack->P() < 1.5) fHistClusterdEtadPhiPosTracksP_05_15BeforeQA->Fill(dEta, dPhi);
+                    else fHistClusterdEtadPhiPosTracksP_15_99BeforeQA->Fill(dEta, dPhi);
                 }
                 else{
                     fHistClusterdEtadPhiNegTracksBeforeQA->Fill(dEta, dPhi);
-                    if(inTrack->P() < 1) fHistClusterdEtadPhiNegTracks_P1_BeforeQA->Fill(dEta, dPhi);
-                    else if(inTrack->P() < 4) fHistClusterdEtadPhiNegTracks_P14_BeforeQA->Fill(dEta, dPhi);
-                    else fHistClusterdEtadPhiNegTracks_P4_BeforeQA->Fill(dEta, dPhi);
+                    if(inTrack->P() < 0.5) fHistClusterdEtadPhiNegTracksP_00_05BeforeQA->Fill(dEta, dPhi);
+                    else if(inTrack->P() < 1.5) fHistClusterdEtadPhiNegTracksP_05_15BeforeQA->Fill(dEta, dPhi);
+                    else fHistClusterdEtadPhiNegTracksP_15_99BeforeQA->Fill(dEta, dPhi);
                 }
                 fHistClusterdEtadPtBeforeQA->Fill(dEta, inTrack->Pt());
                 fHistClusterdPhidPtBeforeQA->Fill(dPhi, inTrack->Pt());
@@ -1346,57 +1382,57 @@ Bool_t AliCaloPhotonCuts::SetTrackMatchingCut(Int_t trackMatching)
                break;
        case 1: 
                if (!fUseDistTrackToCluster) fUseDistTrackToCluster=1;
-        fMaxDistTrackToClusterEta = 0.015;
-        fMinDistTrackToClusterPhi = -0.01;
-        fMaxDistTrackToClusterPhi = 0.03;      //0.04;
+        fMaxDistTrackToClusterEta = 0.008;      //0.015;
+        fMinDistTrackToClusterPhi = -0.03;      //-0.01;
+        fMaxDistTrackToClusterPhi = 0.03;       //0.03;        //0.04;
                break;
        case 2: 
                if (!fUseDistTrackToCluster) fUseDistTrackToCluster=1;
-        fMaxDistTrackToClusterEta = 0.015;
-        fMinDistTrackToClusterPhi = -0.01;
-        fMaxDistTrackToClusterPhi = 0.035;     //0.05;
+        fMaxDistTrackToClusterEta = 0.012;      //0.015;
+        fMinDistTrackToClusterPhi = -0.05;      //-0.01;
+        fMaxDistTrackToClusterPhi = 0.04;       //0.035;       //0.05;
                break;
        case 3: 
                if (!fUseDistTrackToCluster) fUseDistTrackToCluster=1;
-        fMaxDistTrackToClusterEta = 0.015;
-        fMinDistTrackToClusterPhi = -0.015;
-        fMaxDistTrackToClusterPhi = 0.04;      //0.1;
+        fMaxDistTrackToClusterEta = 0.016;      //0.015;
+        fMinDistTrackToClusterPhi = -0.09;      //-0.015;
+        fMaxDistTrackToClusterPhi = 0.06;       //0.04;        //0.1;
                break;
        case 4: 
                if (!fUseDistTrackToCluster) fUseDistTrackToCluster=1;
-        fMaxDistTrackToClusterEta = 0.015;
-        fMinDistTrackToClusterPhi = -0.015;
-        fMaxDistTrackToClusterPhi = 0.045;     //0.13;
+        fMaxDistTrackToClusterEta = 0.018;      //0.015;
+        fMinDistTrackToClusterPhi = -0.11;      //-0.015;
+        fMaxDistTrackToClusterPhi = 0.07;       //0.045;       //0.13;
                break;
        case 5: 
                if (!fUseDistTrackToCluster) fUseDistTrackToCluster=1;
-        fMaxDistTrackToClusterEta = 0.015;
-        fMinDistTrackToClusterPhi = -0.02;
-        fMaxDistTrackToClusterPhi = 0.05;      //0.15
+        fMaxDistTrackToClusterEta = 0.02;       //0.015;
+        fMinDistTrackToClusterPhi = -0.13;      //-0.02;
+        fMaxDistTrackToClusterPhi = 0.08;       //0.05;        //0.15
                break;
        case 6: 
                if (!fUseDistTrackToCluster) fUseDistTrackToCluster=1;
-        fMaxDistTrackToClusterEta = 0.015;
-        fMinDistTrackToClusterPhi = -0.02;
-        fMaxDistTrackToClusterPhi = 0.055;     //0.2;
+        fMaxDistTrackToClusterEta = 0.022;      //0.015;
+        fMinDistTrackToClusterPhi = -0.15;      //-0.02;
+        fMaxDistTrackToClusterPhi = 0.10;       //0.055;       //0.2;
                break;
-       case 7: 
+    case 7: //PHOS
                if (!fUseDistTrackToCluster) fUseDistTrackToCluster=1;
-        fMaxDistTrackToClusterEta = 0.015;
-        fMinDistTrackToClusterPhi = -0.025;
-        fMaxDistTrackToClusterPhi = 0.06;      //0.3;
+        fMaxDistTrackToClusterEta = 0.005;      //0.015;
+        fMinDistTrackToClusterPhi = -0.03;      //-0.025;
+        fMaxDistTrackToClusterPhi = 0.03;       //0.06;        //0.3;
                break;
-       case 8: 
+    case 8: //PHOS
                if (!fUseDistTrackToCluster) fUseDistTrackToCluster=1;
-        fMaxDistTrackToClusterEta = 0.015;
-        fMinDistTrackToClusterPhi = -0.025;
-        fMaxDistTrackToClusterPhi = 0.07;      //0.4;
+        fMaxDistTrackToClusterEta = 0.01;       //0.015;
+        fMinDistTrackToClusterPhi = -0.09;      //-0.025;
+        fMaxDistTrackToClusterPhi = 0.07;       //0.07;        //0.4;
                break;
-       case 9: 
+    case 9: //PHOS
                if (!fUseDistTrackToCluster) fUseDistTrackToCluster=1;
-        fMaxDistTrackToClusterEta = 0.02;
-        fMinDistTrackToClusterPhi = -0.03;
-        fMaxDistTrackToClusterPhi = 0.1;       //0.5;
+        fMaxDistTrackToClusterEta = 0.015;      //0.02;
+        fMinDistTrackToClusterPhi = -0.15;      //-0.03;
+        fMaxDistTrackToClusterPhi = 0.11;       //0.1; //0.5;
                break;
 
        default:
index 0938091..d295be7 100644 (file)
@@ -207,12 +207,12 @@ class AliCaloPhotonCuts : public AliAnalysisCuts {
         TH2F*          fHistClusterdEtadPhiNegTracksBeforeQA;  // 2-dim plot dEta vs. dPhi
         TH2F*          fHistClusterdEtadPhiPosTracksAfterQA;   // 2-dim plot dEta vs. dPhi
         TH2F*          fHistClusterdEtadPhiNegTracksAfterQA;   // 2-dim plot dEta vs. dPhi
-        TH2F*       fHistClusterdEtadPhiPosTracks_P1_BeforeQA;   // 2-dim plot dEta vs. dPhi, positive Tracks, P < 1
-        TH2F*       fHistClusterdEtadPhiPosTracks_P14_BeforeQA;  // 2-dim plot dEta vs. dPhi, positive Tracks, 1 < P < 4
-        TH2F*       fHistClusterdEtadPhiPosTracks_P4_BeforeQA;   // 2-dim plot dEta vs. dPhi, positive Tracks, P > 4
-        TH2F*       fHistClusterdEtadPhiNegTracks_P1_BeforeQA;   // 2-dim plot dEta vs. dPhi, negative Tracks, P < 1
-        TH2F*       fHistClusterdEtadPhiNegTracks_P14_BeforeQA;  // 2-dim plot dEta vs. dPhi, negative Tracks, 1 < P < 4
-        TH2F*       fHistClusterdEtadPhiNegTracks_P4_BeforeQA;   // 2-dim plot dEta vs. dPhi, negative Tracks, P > 4
+        TH2F*       fHistClusterdEtadPhiPosTracksP_00_05BeforeQA;// 2-dim plot dEta vs. dPhi, positive Tracks, P < 0.5
+        TH2F*       fHistClusterdEtadPhiPosTracksP_05_15BeforeQA;// 2-dim plot dEta vs. dPhi, positive Tracks, 0.5 < P < 1.5
+        TH2F*       fHistClusterdEtadPhiPosTracksP_15_99BeforeQA;// 2-dim plot dEta vs. dPhi, positive Tracks, P > 1.5
+        TH2F*       fHistClusterdEtadPhiNegTracksP_00_05BeforeQA;// 2-dim plot dEta vs. dPhi, negative Tracks, P < 0.5
+        TH2F*       fHistClusterdEtadPhiNegTracksP_05_15BeforeQA;// 2-dim plot dEta vs. dPhi, negative Tracks, 0.5 < P < 1.5
+        TH2F*       fHistClusterdEtadPhiNegTracksP_15_99BeforeQA;// 2-dim plot dEta vs. dPhi, negative Tracks, P > 1.5
         TH2F*          fHistClusterdEtadPtBeforeQA;                    // 2-dim plot dEta vs. Pt
         TH2F*          fHistClusterdPhidPtBeforeQA;                    // 2-dim plot dEta vs. Pt
         TH2F*          fHistClusterM20Pt_dPhiBeforeQA;                 // 2-dim plot M20 vs. Pt for given dPhi>0.05
index fbe92f1..73f5a29 100644 (file)
@@ -156,41 +156,41 @@ void AddTask_GammaConvCalo_PbPb(          Int_t trainConfig = 1,  //change different s
                eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002000092970028250400000"; clusterCutArray[3] = "10000043032030000"; mesonCutArray[ 3] = "01631031000000"; // 40-80%
                eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002000092970028250400000"; clusterCutArray[4] = "10000043032030000"; mesonCutArray[ 4] = "01631031000000"; // 20-50%                               
        } else if (trainConfig == 31){ // PHOS clusters
-               eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002000092970028250400000"; clusterCutArray[0] = "20000042033200000"; mesonCutArray[ 0] = "01631031000000"; // 0-5%
-               eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002000092970028250400000"; clusterCutArray[1] = "20000042033200000"; mesonCutArray[ 1] = "01631031000000"; // 5-10%
-               eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002000092970028250400000"; clusterCutArray[2] = "20000042033200000"; mesonCutArray[ 2] = "01631031000000"; // 0-10%
-               eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002000092970028250400000"; clusterCutArray[3] = "20000042033200000"; mesonCutArray[ 3] = "01631031000000"; // 20-40%
-               eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002000092970028250400000"; clusterCutArray[4] = "20000042033200000"; mesonCutArray[ 4] = "01631031000000"; // 20-50%
+        eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002000092970028250400000"; clusterCutArray[0] = "20000048033200000"; mesonCutArray[ 0] = "01631031000000"; // 0-5%
+        eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002000092970028250400000"; clusterCutArray[1] = "20000048033200000"; mesonCutArray[ 1] = "01631031000000"; // 5-10%
+        eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002000092970028250400000"; clusterCutArray[2] = "20000048033200000"; mesonCutArray[ 2] = "01631031000000"; // 0-10%
+        eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002000092970028250400000"; clusterCutArray[3] = "20000048033200000"; mesonCutArray[ 3] = "01631031000000"; // 20-40%
+        eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002000092970028250400000"; clusterCutArray[4] = "20000048033200000"; mesonCutArray[ 4] = "01631031000000"; // 20-50%
        } else if (trainConfig == 32){ // PHOS clusters
-               eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002000092970028250400000"; clusterCutArray[0] = "20000042033200000"; mesonCutArray[ 0] = "01631031000000"; // 0-5%
-               eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002000092970028250400000"; clusterCutArray[1] = "20000042033200000"; mesonCutArray[ 1] = "01631031000000"; // 5-10%
-               eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002000092970028250400000"; clusterCutArray[2] = "20000042033200000"; mesonCutArray[ 2] = "01631031000000"; // 0-10%
-               eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002000092970028250400000"; clusterCutArray[3] = "20000042033200000"; mesonCutArray[ 3] = "01631031000000"; // 10-20%
-               eventCutArray[ 4] = "5240001"; photonCutArray[ 4] = "002000092970028250400000"; clusterCutArray[4] = "20000042033200000"; mesonCutArray[ 4] = "01631031000000"; // 20-40%               
+        eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002000092970028250400000"; clusterCutArray[0] = "20000048033200000"; mesonCutArray[ 0] = "01631031000000"; // 0-5%
+        eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002000092970028250400000"; clusterCutArray[1] = "20000048033200000"; mesonCutArray[ 1] = "01631031000000"; // 5-10%
+        eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002000092970028250400000"; clusterCutArray[2] = "20000048033200000"; mesonCutArray[ 2] = "01631031000000"; // 0-10%
+        eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002000092970028250400000"; clusterCutArray[3] = "20000048033200000"; mesonCutArray[ 3] = "01631031000000"; // 10-20%
+        eventCutArray[ 4] = "5240001"; photonCutArray[ 4] = "002000092970028250400000"; clusterCutArray[4] = "20000048033200000"; mesonCutArray[ 4] = "01631031000000"; // 20-40%
        } else if (trainConfig == 33){ // PHOS clusters
-               eventCutArray[ 0] = "5460001"; photonCutArray[ 0] = "002000092970028250400000"; clusterCutArray[0] = "20000042033200000"; mesonCutArray[ 0] = "01631031000000"; // 40-60%
-               eventCutArray[ 1] = "5680001"; photonCutArray[ 1] = "002000092970028250400000"; clusterCutArray[1] = "20000042033200000"; mesonCutArray[ 1] = "01631031000000"; // 60-80%
-               eventCutArray[ 2] = "5260001"; photonCutArray[ 2] = "002000092970028250400000"; clusterCutArray[2] = "20000042033200000"; mesonCutArray[ 2] = "01631031000000"; // 20-60%
-               eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002000092970028250400000"; clusterCutArray[3] = "20000042033200000"; mesonCutArray[ 3] = "01631031000000"; // 40-80%
-               eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002000092970028250400000"; clusterCutArray[4] = "20000042033200000"; mesonCutArray[ 4] = "01631031000000"; // 20-50%                                               
+        eventCutArray[ 0] = "5460001"; photonCutArray[ 0] = "002000092970028250400000"; clusterCutArray[0] = "20000048033200000"; mesonCutArray[ 0] = "01631031000000"; // 40-60%
+        eventCutArray[ 1] = "5680001"; photonCutArray[ 1] = "002000092970028250400000"; clusterCutArray[1] = "20000048033200000"; mesonCutArray[ 1] = "01631031000000"; // 60-80%
+        eventCutArray[ 2] = "5260001"; photonCutArray[ 2] = "002000092970028250400000"; clusterCutArray[2] = "20000048033200000"; mesonCutArray[ 2] = "01631031000000"; // 20-60%
+        eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002000092970028250400000"; clusterCutArray[3] = "20000048033200000"; mesonCutArray[ 3] = "01631031000000"; // 40-80%
+        eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002000092970028250400000"; clusterCutArray[4] = "20000048033200000"; mesonCutArray[ 4] = "01631031000000"; // 20-50%
        } else if (trainConfig == 34){ // PHOS clusters
-               eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002000092970028250400000"; clusterCutArray[0] = "20000043033200000"; mesonCutArray[ 0] = "01631031000000"; // 0-5%
-               eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002000092970028250400000"; clusterCutArray[1] = "20000043033200000"; mesonCutArray[ 1] = "01631031000000"; // 5-10%
-               eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002000092970028250400000"; clusterCutArray[2] = "20000043033200000"; mesonCutArray[ 2] = "01631031000000"; // 0-10%
-               eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002000092970028250400000"; clusterCutArray[3] = "20000043033200000"; mesonCutArray[ 3] = "01631031000000"; // 20-40%
-               eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002000092970028250400000"; clusterCutArray[4] = "20000043033200000"; mesonCutArray[ 4] = "01631031000000"; // 20-50%
+        eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002000092970028250400000"; clusterCutArray[0] = "20000048033200000"; mesonCutArray[ 0] = "01631031000000"; // 0-5%
+        eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002000092970028250400000"; clusterCutArray[1] = "20000048033200000"; mesonCutArray[ 1] = "01631031000000"; // 5-10%
+        eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002000092970028250400000"; clusterCutArray[2] = "20000048033200000"; mesonCutArray[ 2] = "01631031000000"; // 0-10%
+        eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002000092970028250400000"; clusterCutArray[3] = "20000048033200000"; mesonCutArray[ 3] = "01631031000000"; // 20-40%
+        eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002000092970028250400000"; clusterCutArray[4] = "20000048033200000"; mesonCutArray[ 4] = "01631031000000"; // 20-50%
        } else if (trainConfig == 35){ // PHOS clusters
-               eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002000092970028250400000"; clusterCutArray[0] = "20000043033200000"; mesonCutArray[ 0] = "01631031000000"; // 0-5%
-               eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002000092970028250400000"; clusterCutArray[1] = "20000043033200000"; mesonCutArray[ 1] = "01631031000000"; // 5-10%
-               eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002000092970028250400000"; clusterCutArray[2] = "20000043033200000"; mesonCutArray[ 2] = "01631031000000"; // 0-10%
-               eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002000092970028250400000"; clusterCutArray[3] = "20000043033200000"; mesonCutArray[ 3] = "01631031000000"; // 10-20%
-               eventCutArray[ 4] = "5240001"; photonCutArray[ 4] = "002000092970028250400000"; clusterCutArray[4] = "20000043033200000"; mesonCutArray[ 4] = "01631031000000"; // 20-40%               
+        eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002000092970028250400000"; clusterCutArray[0] = "20000048033200000"; mesonCutArray[ 0] = "01631031000000"; // 0-5%
+        eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002000092970028250400000"; clusterCutArray[1] = "20000048033200000"; mesonCutArray[ 1] = "01631031000000"; // 5-10%
+        eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002000092970028250400000"; clusterCutArray[2] = "20000048033200000"; mesonCutArray[ 2] = "01631031000000"; // 0-10%
+        eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002000092970028250400000"; clusterCutArray[3] = "20000048033200000"; mesonCutArray[ 3] = "01631031000000"; // 10-20%
+        eventCutArray[ 4] = "5240001"; photonCutArray[ 4] = "002000092970028250400000"; clusterCutArray[4] = "20000048033200000"; mesonCutArray[ 4] = "01631031000000"; // 20-40%
        } else if (trainConfig == 36){ // PHOS clusters
-               eventCutArray[ 0] = "5460001"; photonCutArray[ 0] = "002000092970028250400000"; clusterCutArray[0] = "20000043033200000"; mesonCutArray[ 0] = "01631031000000"; // 40-60%
-               eventCutArray[ 1] = "5680001"; photonCutArray[ 1] = "002000092970028250400000"; clusterCutArray[1] = "20000043033200000"; mesonCutArray[ 1] = "01631031000000"; // 60-80%
-               eventCutArray[ 2] = "5260001"; photonCutArray[ 2] = "002000092970028250400000"; clusterCutArray[2] = "20000043033200000"; mesonCutArray[ 2] = "01631031000000"; // 20-60%
-               eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002000092970028250400000"; clusterCutArray[3] = "20000043033200000"; mesonCutArray[ 3] = "01631031000000"; // 40-80%
-               eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002000092970028250400000"; clusterCutArray[4] = "20000043033200000"; mesonCutArray[ 4] = "01631031000000"; // 20-50%                                                               
+        eventCutArray[ 0] = "5460001"; photonCutArray[ 0] = "002000092970028250400000"; clusterCutArray[0] = "20000048033200000"; mesonCutArray[ 0] = "01631031000000"; // 40-60%
+        eventCutArray[ 1] = "5680001"; photonCutArray[ 1] = "002000092970028250400000"; clusterCutArray[1] = "20000048033200000"; mesonCutArray[ 1] = "01631031000000"; // 60-80%
+        eventCutArray[ 2] = "5260001"; photonCutArray[ 2] = "002000092970028250400000"; clusterCutArray[2] = "20000048033200000"; mesonCutArray[ 2] = "01631031000000"; // 20-60%
+        eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002000092970028250400000"; clusterCutArray[3] = "20000048033200000"; mesonCutArray[ 3] = "01631031000000"; // 40-80%
+        eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002000092970028250400000"; clusterCutArray[4] = "20000048033200000"; mesonCutArray[ 4] = "01631031000000"; // 20-50%
        } else {
                Error(Form("GammaConvCalo_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
                return;
index 2c59fd3..ce98f90 100644 (file)
@@ -112,11 +112,14 @@ void AddTask_GammaConvCalo_pPb(  Int_t trainConfig = 1,  //change different set
        task= new AliAnalysisTaskGammaConvCalo(Form("GammaConvCalo_%i",trainConfig));
        task->SetIsHeavyIon(isHeavyIon);
        task->SetIsMC(isMC);
+
        // Cut Numbers to use in Analysis
        Int_t numberOfCuts = 2;
-       if (trainConfig==10 ){ numberOfCuts =4;}
-       if (trainConfig==2 || trainConfig==4  ){ numberOfCuts =6;}
-       if (trainConfig==5 || trainConfig==6 || trainConfig==7 || trainConfig==8 || trainConfig==32 || trainConfig==33){ numberOfCuts =5;}
+
+    if ( trainConfig==32 || trainConfig==33 ) {numberOfCuts = 3;}
+    if ( trainConfig==10 ) {numberOfCuts = 4;}
+    if ( trainConfig==7 || trainConfig==8 ) {numberOfCuts = 5;}
+    if ( trainConfig==2 || trainConfig==4 || trainConfig==5 || trainConfig==6 ) {numberOfCuts = 6;}
        
        TString *eventCutArray = new TString[numberOfCuts];
        TString *photonCutArray = new TString[numberOfCuts];
@@ -129,74 +132,72 @@ void AddTask_GammaConvCalo_pPb(  Int_t trainConfig = 1,  //change different set
        
        //************************************************ EMCAL clusters **********************************************************
        if (trainConfig == 1){ // min energy = 0.3 GeV/c
-               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "10000042022030000"; mesonCutArray[0] = "01631031000000"; //standart cut, kINT7 
-               eventCutArray[ 1] = "8005201"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "10000042022030000"; mesonCutArray[1] = "01631031000000"; //standard cut, kEMC7 
+        eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "10000043022030000"; mesonCutArray[0] = "01631031000000"; //standart cut, kINT7
+        eventCutArray[ 1] = "8005201"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "10000043022030000"; mesonCutArray[1] = "01631031000000"; //standard cut, kEMC7
        } else if (trainConfig == 2){  // min energy = 0.3 GeV/c
-               eventCutArray[ 0] = "8008101"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "10000042022030000"; mesonCutArray[0] = "01631031000000"; //standart cut, kEMCEGA based on INT7 
-               eventCutArray[ 1] = "8008301"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "10000042022030000"; mesonCutArray[1] = "01631031000000"; //standard cut, kEMCEG1 based on INT7 
-               eventCutArray[ 2] = "8008501"; photonCutArray[ 2] = "002000093270028250400000"; clusterCutArray[2] = "10000042022030000"; mesonCutArray[2] = "01631031000000"; //standard cut, kEMCEG2 based on INT7 
-               eventCutArray[ 3] = "8009101"; photonCutArray[ 3] = "002000093270028250400000"; clusterCutArray[3] = "10000042022030000"; mesonCutArray[3] = "01631031000000"; //standard cut, kEMCEJE based on INT7 
-               eventCutArray[ 4] = "8009301"; photonCutArray[ 4] = "002000093270028250400000"; clusterCutArray[4] = "10000042022030000"; mesonCutArray[4] = "01631031000000"; //standard cut, kEMCEJ1 based on INT7 
-               eventCutArray[ 5] = "8009501"; photonCutArray[ 5] = "002000093270028250400000"; clusterCutArray[5] = "10000042022030000"; mesonCutArray[5] = "01631031000000"; //standard cut, kEMCEG2 based on INT7 
+        eventCutArray[ 0] = "8008101"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "10000043022030000"; mesonCutArray[0] = "01631031000000"; //standart cut, kEMCEGA based on INT7
+        eventCutArray[ 1] = "8008301"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "10000043022030000"; mesonCutArray[1] = "01631031000000"; //standard cut, kEMCEG1 based on INT7
+        eventCutArray[ 2] = "8008501"; photonCutArray[ 2] = "002000093270028250400000"; clusterCutArray[2] = "10000043022030000"; mesonCutArray[2] = "01631031000000"; //standard cut, kEMCEG2 based on INT7
+        eventCutArray[ 3] = "8009101"; photonCutArray[ 3] = "002000093270028250400000"; clusterCutArray[3] = "10000043022030000"; mesonCutArray[3] = "01631031000000"; //standard cut, kEMCEJE based on INT7
+        eventCutArray[ 4] = "8009301"; photonCutArray[ 4] = "002000093270028250400000"; clusterCutArray[4] = "10000043022030000"; mesonCutArray[4] = "01631031000000"; //standard cut, kEMCEJ1 based on INT7
+        eventCutArray[ 5] = "8009501"; photonCutArray[ 5] = "002000093270028250400000"; clusterCutArray[5] = "10000043022030000"; mesonCutArray[5] = "01631031000000"; //standard cut, kEMCEG2 based on INT7
        } else if (trainConfig == 3){ // min energy = 0.4 GeV/c
-               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "10000042032030000"; mesonCutArray[0] = "01631031000000"; //standart cut, kINT7 
-               eventCutArray[ 1] = "8005201"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01631031000000"; //standard cut, kEMC7    
+        eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "10000043032030000"; mesonCutArray[0] = "01631031000000"; //standart cut, kINT7
+        eventCutArray[ 1] = "8005201"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "10000043032030000"; mesonCutArray[1] = "01631031000000"; //standard cut, kEMC7
        } else if (trainConfig == 4){ // min energy = 0.4 GeV/c
-               eventCutArray[ 0] = "8008101"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "10000042032030000"; mesonCutArray[0] = "01631031000000"; //standart cut, kEMCEGA based on INT7 
-               eventCutArray[ 1] = "8008301"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01631031000000"; //standard cut, kEMCEG1 based on INT7 
-               eventCutArray[ 2] = "8008501"; photonCutArray[ 2] = "002000093270028250400000"; clusterCutArray[2] = "10000042032030000"; mesonCutArray[2] = "01631031000000"; //standard cut, kEMCEG2 based on INT7 
-               eventCutArray[ 3] = "8009101"; photonCutArray[ 3] = "002000093270028250400000"; clusterCutArray[3] = "10000042032030000"; mesonCutArray[3] = "01631031000000"; //standard cut, kEMCEJE based on INT7 
-               eventCutArray[ 4] = "8009301"; photonCutArray[ 4] = "002000093270028250400000"; clusterCutArray[4] = "10000042032030000"; mesonCutArray[4] = "01631031000000"; //standard cut, kEMCEJ1 based on INT7 
-               eventCutArray[ 5] = "8009501"; photonCutArray[ 5] = "002000093270028250400000"; clusterCutArray[5] = "10000042032030000"; mesonCutArray[5] = "01631031000000"; //standard cut, kEMCEG2 based on INT7
+        eventCutArray[ 0] = "8008101"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "10000043032030000"; mesonCutArray[0] = "01631031000000"; //standart cut, kEMCEGA based on INT7
+        eventCutArray[ 1] = "8008301"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "10000043032030000"; mesonCutArray[1] = "01631031000000"; //standard cut, kEMCEG1 based on INT7
+        eventCutArray[ 2] = "8008501"; photonCutArray[ 2] = "002000093270028250400000"; clusterCutArray[2] = "10000043032030000"; mesonCutArray[2] = "01631031000000"; //standard cut, kEMCEG2 based on INT7
+        eventCutArray[ 3] = "8009101"; photonCutArray[ 3] = "002000093270028250400000"; clusterCutArray[3] = "10000043032030000"; mesonCutArray[3] = "01631031000000"; //standard cut, kEMCEJE based on INT7
+        eventCutArray[ 4] = "8009301"; photonCutArray[ 4] = "002000093270028250400000"; clusterCutArray[4] = "10000043032030000"; mesonCutArray[4] = "01631031000000"; //standard cut, kEMCEJ1 based on INT7
+        eventCutArray[ 5] = "8009501"; photonCutArray[ 5] = "002000093270028250400000"; clusterCutArray[5] = "10000043032030000"; mesonCutArray[5] = "01631031000000"; //standard cut, kEMCEG2 based on INT7
        } else if (trainConfig == 5){ //EMCAL variation of track matching
                eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "10000041032030000"; mesonCutArray[0] = "01631031000000"; // 
                eventCutArray[ 1] = "8000001"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01631031000000"; 
                eventCutArray[ 2] = "8000001"; photonCutArray[ 2] = "002000093270028250400000"; clusterCutArray[2] = "10000043032030000"; mesonCutArray[2] = "01631031000000"; 
                eventCutArray[ 3] = "8000001"; photonCutArray[ 3] = "002000093270028250400000"; clusterCutArray[3] = "10000044032030000"; mesonCutArray[3] = "01631031000000"; 
                eventCutArray[ 4] = "8000001"; photonCutArray[ 4] = "002000093270028250400000"; clusterCutArray[4] = "10000045032030000"; mesonCutArray[4] = "01631031000000"; 
+        eventCutArray[ 5] = "8000001"; photonCutArray[ 5] = "002000093270028250400000"; clusterCutArray[5] = "10000046032030000"; mesonCutArray[5] = "01631031000000";
        } else if (trainConfig == 6){ //EMCAL added signal
                eventCutArray[ 0] = "8000002"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "10000041032030000"; mesonCutArray[0] = "01631031000000"; 
                eventCutArray[ 1] = "8000002"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01631031000000"; 
                eventCutArray[ 2] = "8000002"; photonCutArray[ 2] = "002000093270028250400000"; clusterCutArray[2] = "10000043032030000"; mesonCutArray[2] = "01631031000000"; 
                eventCutArray[ 3] = "8000002"; photonCutArray[ 3] = "002000093270028250400000"; clusterCutArray[3] = "10000044032030000"; mesonCutArray[3] = "01631031000000"; 
                eventCutArray[ 4] = "8000002"; photonCutArray[ 4] = "002000093270028250400000"; clusterCutArray[4] = "10000045032030000"; mesonCutArray[4] = "01631031000000"; 
+        eventCutArray[ 5] = "8000002"; photonCutArray[ 5] = "002000093270028250400000"; clusterCutArray[5] = "10000046032030000"; mesonCutArray[5] = "01631031000000";
        } else if (trainConfig == 7){ //EMCAL minEnergy variation
-               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "10000042012030000"; mesonCutArray[0] = "01631031000000"; //0.2 GeV/c
-               eventCutArray[ 1] = "8000001"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "10000042022030000"; mesonCutArray[1] = "01631031000000"; //0.3 GeV/c
-               eventCutArray[ 2] = "8000001"; photonCutArray[ 2] = "002000093270028250400000"; clusterCutArray[2] = "10000042032030000"; mesonCutArray[2] = "01631031000000"; //0.4 GeV/c default
-               eventCutArray[ 3] = "8000001"; photonCutArray[ 3] = "002000093270028250400000"; clusterCutArray[3] = "10000042042030000"; mesonCutArray[3] = "01631031000000"; //0.5 GeV/c
-               eventCutArray[ 4] = "8000001"; photonCutArray[ 4] = "002000093270028250400000"; clusterCutArray[4] = "10000042052030000"; mesonCutArray[4] = "01631031000000"; //0.6 GeV/c
+        eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "10000043012030000"; mesonCutArray[0] = "01631031000000"; //0.2 GeV/c
+        eventCutArray[ 1] = "8000001"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "10000043022030000"; mesonCutArray[1] = "01631031000000"; //0.3 GeV/c
+        eventCutArray[ 2] = "8000001"; photonCutArray[ 2] = "002000093270028250400000"; clusterCutArray[2] = "10000043032030000"; mesonCutArray[2] = "01631031000000"; //0.4 GeV/c default
+        eventCutArray[ 3] = "8000001"; photonCutArray[ 3] = "002000093270028250400000"; clusterCutArray[3] = "10000043042030000"; mesonCutArray[3] = "01631031000000"; //0.5 GeV/c
+        eventCutArray[ 4] = "8000001"; photonCutArray[ 4] = "002000093270028250400000"; clusterCutArray[4] = "10000043052030000"; mesonCutArray[4] = "01631031000000"; //0.6 GeV/c
        } else if (trainConfig == 8){ //EMCAL minNCells variation
-               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "10000042031030000"; mesonCutArray[0] = "01631031000000"; //n cells >= 1
-               eventCutArray[ 1] = "8000001"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "10000042033030000"; mesonCutArray[1] = "01631031000000"; //n cells >= 3
-               eventCutArray[ 2] = "8000001"; photonCutArray[ 2] = "002000093270028250400000"; clusterCutArray[2] = "10000042032000000"; mesonCutArray[2] = "01631031000000"; //no M02 cut
-               eventCutArray[ 3] = "8000001"; photonCutArray[ 3] = "002000093270028250400000"; clusterCutArray[3] = "10021042032030000"; mesonCutArray[3] = "01631031000000"; //only modules with TRD infront
-               eventCutArray[ 4] = "8000001"; photonCutArray[ 4] = "002000093270028250400000"; clusterCutArray[4] = "10013042032030000"; mesonCutArray[4] = "01631031000000"; //no modules with TRD infront            
+        eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "10000043031030000"; mesonCutArray[0] = "01631031000000"; //n cells >= 1
+        eventCutArray[ 1] = "8000001"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "10000043033030000"; mesonCutArray[1] = "01631031000000"; //n cells >= 3
+        eventCutArray[ 2] = "8000001"; photonCutArray[ 2] = "002000093270028250400000"; clusterCutArray[2] = "10000043032000000"; mesonCutArray[2] = "01631031000000"; //no M02 cut
+        eventCutArray[ 3] = "8000001"; photonCutArray[ 3] = "002000093270028250400000"; clusterCutArray[3] = "10021043032030000"; mesonCutArray[3] = "01631031000000"; //only modules with TRD infront
+        eventCutArray[ 4] = "8000001"; photonCutArray[ 4] = "002000093270028250400000"; clusterCutArray[4] = "10013043032030000"; mesonCutArray[4] = "01631031000000"; //no modules with TRD infront
        } else if (trainConfig == 9){ //PCM restriction in acceptance 
-               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002022093270028250400000"; clusterCutArray[0] = "10000042032030000"; mesonCutArray[0] = "01631031000000"; // PCM photons pointing to EMCAL loose
-               eventCutArray[ 1] = "8000001"; photonCutArray[ 1] = "002044093270028250400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01631031000000"; // PCM photons pointing to EMCAL tight
+        eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002022093270028250400000"; clusterCutArray[0] = "10000043032030000"; mesonCutArray[0] = "01631031000000"; // PCM photons pointing to EMCAL loose
+        eventCutArray[ 1] = "8000001"; photonCutArray[ 1] = "002044093270028250400000"; clusterCutArray[1] = "10000043032030000"; mesonCutArray[1] = "01631031000000"; // PCM photons pointing to EMCAL tight
        } else if (trainConfig == 10){ 
-               eventCutArray[ 0] = "8005201"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "10000042062030000"; mesonCutArray[0] = "01631031000000"; 
-               eventCutArray[ 1] = "8005201"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "10000042072030000"; mesonCutArray[1] = "01631031000000"; 
-               eventCutArray[ 2] = "8005201"; photonCutArray[ 2] = "002000093270028250400000"; clusterCutArray[2] = "10000042082030000"; mesonCutArray[2] = "01631031000000"; 
-               eventCutArray[ 3] = "8005201"; photonCutArray[ 3] = "002000093270028250400000"; clusterCutArray[3] = "10000042092030000"; mesonCutArray[3] = "01631031000000"; 
+        eventCutArray[ 0] = "8005201"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "10000043062030000"; mesonCutArray[0] = "01631031000000";
+        eventCutArray[ 1] = "8005201"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "10000043072030000"; mesonCutArray[1] = "01631031000000";
+        eventCutArray[ 2] = "8005201"; photonCutArray[ 2] = "002000093270028250400000"; clusterCutArray[2] = "10000043082030000"; mesonCutArray[2] = "01631031000000";
+        eventCutArray[ 3] = "8005201"; photonCutArray[ 3] = "002000093270028250400000"; clusterCutArray[3] = "10000043092030000"; mesonCutArray[3] = "01631031000000";
                
        //************************************************ PHOS clusters **********************************************************     
        } else if (trainConfig == 31) { // min energy = 0.3 GeV/c
-               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "20000042033200000"; mesonCutArray[0] = "01631031000000"; //standart cut, kINT7 
-               eventCutArray[ 1] = "8006201"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "20000042033200000"; mesonCutArray[1] = "01631031000000"; //standard cut, kPHI7            
+        eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "20000048033200000"; mesonCutArray[0] = "01631031000000"; //standart cut, kINT7
+        eventCutArray[ 1] = "8006201"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "20000048033200000"; mesonCutArray[1] = "01631031000000"; //standard cut, kPHI7
        } else if (trainConfig == 32) { //PHOS
-               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "20000041033200000"; mesonCutArray[0] = "01631031000000"; 
-               eventCutArray[ 1] = "8000001"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "20000042033200000"; mesonCutArray[1] = "01631031000000"; 
-               eventCutArray[ 2] = "8000001"; photonCutArray[ 2] = "002000093270028250400000"; clusterCutArray[2] = "20000043033200000"; mesonCutArray[2] = "01631031000000"; 
-               eventCutArray[ 3] = "8000001"; photonCutArray[ 3] = "002000093270028250400000"; clusterCutArray[3] = "20000044033200000"; mesonCutArray[3] = "01631031000000"; 
-               eventCutArray[ 4] = "8000001"; photonCutArray[ 4] = "002000093270028250400000"; clusterCutArray[4] = "20000045033200000"; mesonCutArray[4] = "01631031000000"; 
+        eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "20000047033200000"; mesonCutArray[0] = "01631031000000";
+        eventCutArray[ 1] = "8000001"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "20000048033200000"; mesonCutArray[1] = "01631031000000";
+        eventCutArray[ 2] = "8000001"; photonCutArray[ 2] = "002000093270028250400000"; clusterCutArray[2] = "20000049033200000"; mesonCutArray[2] = "01631031000000";
        } else if (trainConfig == 33) { //PHOS
-               eventCutArray[ 0] = "8000002"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "20000041033200000"; mesonCutArray[0] = "01631031000000"; 
-               eventCutArray[ 1] = "8000002"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "20000042033200000"; mesonCutArray[1] = "01631031000000"; 
-               eventCutArray[ 2] = "8000002"; photonCutArray[ 2] = "002000093270028250400000"; clusterCutArray[2] = "20000043033200000"; mesonCutArray[2] = "01631031000000"; 
-               eventCutArray[ 3] = "8000002"; photonCutArray[ 3] = "002000093270028250400000"; clusterCutArray[3] = "20000044033200000"; mesonCutArray[3] = "01631031000000"; 
-               eventCutArray[ 4] = "8000002"; photonCutArray[ 4] = "002000093270028250400000"; clusterCutArray[4] = "20000045033200000"; mesonCutArray[4] = "01631031000000"; 
+        eventCutArray[ 0] = "8000002"; photonCutArray[ 0] = "002000093270028250400000"; clusterCutArray[0] = "20000047033200000"; mesonCutArray[0] = "01631031000000";
+        eventCutArray[ 1] = "8000002"; photonCutArray[ 1] = "002000093270028250400000"; clusterCutArray[1] = "20000048033200000"; mesonCutArray[1] = "01631031000000";
+        eventCutArray[ 2] = "8000002"; photonCutArray[ 2] = "002000093270028250400000"; clusterCutArray[2] = "20000049033200000"; mesonCutArray[2] = "01631031000000";
        } else {
                Error(Form("GammaConvCalo_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
                return;
index 6937f5c..bfb5b1a 100644 (file)
@@ -112,9 +112,9 @@ void AddTask_GammaConvCalo_pp(  Int_t trainConfig = 1,  //change different set o
        task->SetIsMC(isMC);
        // Cut Numbers to use in Analysis
        Int_t numberOfCuts = 2;
-       if (trainConfig==2 || trainConfig==3 || trainConfig==4 || trainConfig==5 || trainConfig==6 || trainConfig==7 ){ numberOfCuts =5;}
-       if (trainConfig==8 || trainConfig==10){ numberOfCuts =4;}
-       if (trainConfig==11 || trainConfig==31 || trainConfig==32 ){ numberOfCuts =6;}
+    if (trainConfig==8 || trainConfig==10) {numberOfCuts = 4;}
+    if (trainConfig==2 || trainConfig==3 || trainConfig==5 || trainConfig==6 || trainConfig==7 ) {numberOfCuts = 5;}
+    if (trainConfig==4 || trainConfig==11 || trainConfig==31 || trainConfig==32 ) {numberOfCuts = 6;}
 
        TString *eventCutArray = new TString[numberOfCuts];
        TString *photonCutArray = new TString[numberOfCuts];
@@ -128,83 +128,84 @@ void AddTask_GammaConvCalo_pp(  Int_t trainConfig = 1,  //change different set o
        // ************************************* EMCAL cuts ****************************************************
        // LHC11a
        if (trainConfig == 1){ // EMCAL clusters 2.76 TeV LHC11a, with SDD (0), kEMC1 (1)
-               eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002000093270008250400000"; clusterCutArray[0] = "10000042032030000"; mesonCutArray[0] = "01631031000000"; // 400 MeV cluster min energy
-               eventCutArray[ 1] = "0005111"; photonCutArray[ 1] = "002000093270008250400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01631031000000"; // 400 MeV cluster min energy
+        eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002000093270008250400000"; clusterCutArray[0] = "10000043032030000"; mesonCutArray[0] = "01631031000000"; // 400 MeV cluster min energy
+        eventCutArray[ 1] = "0005111"; photonCutArray[ 1] = "002000093270008250400000"; clusterCutArray[1] = "10000043032030000"; mesonCutArray[1] = "01631031000000"; // 400 MeV cluster min energy
        } else if (trainConfig == 2){ //EMCAL minEnergy variation
-               eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002000093270008250400000"; clusterCutArray[0] = "10000042012030000"; mesonCutArray[0] = "01631031000000"; //0.2 GeV/c
-               eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002000093270008250400000"; clusterCutArray[1] = "10000042022030000"; mesonCutArray[1] = "01631031000000"; //0.3 GeV/c
-               eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002000093270008250400000"; clusterCutArray[2] = "10000042032030000"; mesonCutArray[2] = "01631031000000"; //0.4 GeV/c default
-               eventCutArray[ 3] = "0000311"; photonCutArray[ 3] = "002000093270008250400000"; clusterCutArray[3] = "10000042042030000"; mesonCutArray[3] = "01631031000000"; //0.5 GeV/c
-               eventCutArray[ 4] = "0000311"; photonCutArray[ 4] = "002000093270008250400000"; clusterCutArray[4] = "10000042052030000"; mesonCutArray[4] = "01631031000000"; //0.6 GeV/c
+        eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002000093270008250400000"; clusterCutArray[0] = "10000043012030000"; mesonCutArray[0] = "01631031000000"; //0.2 GeV/c
+        eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002000093270008250400000"; clusterCutArray[1] = "10000043022030000"; mesonCutArray[1] = "01631031000000"; //0.3 GeV/c
+        eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002000093270008250400000"; clusterCutArray[2] = "10000043032030000"; mesonCutArray[2] = "01631031000000"; //0.4 GeV/c default
+        eventCutArray[ 3] = "0000311"; photonCutArray[ 3] = "002000093270008250400000"; clusterCutArray[3] = "10000043042030000"; mesonCutArray[3] = "01631031000000"; //0.5 GeV/c
+        eventCutArray[ 4] = "0000311"; photonCutArray[ 4] = "002000093270008250400000"; clusterCutArray[4] = "10000043052030000"; mesonCutArray[4] = "01631031000000"; //0.6 GeV/c
        } else if (trainConfig == 3){ //EMCAL minNCells variation
-               eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002000093270008250400000"; clusterCutArray[0] = "10000042031030000"; mesonCutArray[0] = "01631031000000"; //n cells >= 1
-               eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002000093270008250400000"; clusterCutArray[1] = "10000042033030000"; mesonCutArray[1] = "01631031000000"; //n cells >= 3
-               eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002000093270008250400000"; clusterCutArray[2] = "10000042032000000"; mesonCutArray[2] = "01631031000000"; //no M02 cut
-               eventCutArray[ 3] = "0000311"; photonCutArray[ 3] = "002000093270008250400000"; clusterCutArray[3] = "10031042032030000"; mesonCutArray[3] = "01631031000000"; //only modules with TRD infront
-               eventCutArray[ 4] = "0000311"; photonCutArray[ 4] = "002000093270008250400000"; clusterCutArray[4] = "10012042032030000"; mesonCutArray[4] = "01631031000000"; //no modules with TRD infront            
+        eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002000093270008250400000"; clusterCutArray[0] = "10000043031030000"; mesonCutArray[0] = "01631031000000"; //n cells >= 1
+        eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002000093270008250400000"; clusterCutArray[1] = "10000043033030000"; mesonCutArray[1] = "01631031000000"; //n cells >= 3
+        eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002000093270008250400000"; clusterCutArray[2] = "10000043032000000"; mesonCutArray[2] = "01631031000000"; //no M02 cut
+        eventCutArray[ 3] = "0000311"; photonCutArray[ 3] = "002000093270008250400000"; clusterCutArray[3] = "10031043032030000"; mesonCutArray[3] = "01631031000000"; //only modules with TRD infront
+        eventCutArray[ 4] = "0000311"; photonCutArray[ 4] = "002000093270008250400000"; clusterCutArray[4] = "10012043032030000"; mesonCutArray[4] = "01631031000000"; //no modules with TRD infront
        } else if (trainConfig == 4){ // EMCAL track matching variations 
                eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002000093270008250400000"; clusterCutArray[0] = "10000041032030000"; mesonCutArray[0] = "01631031000000"; // 
                eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002000093270008250400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01631031000000"; // 
                eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002000093270008250400000"; clusterCutArray[2] = "10000043032030000"; mesonCutArray[2] = "01631031000000"; // 
                eventCutArray[ 3] = "0000311"; photonCutArray[ 3] = "002000093270008250400000"; clusterCutArray[3] = "10000044032030000"; mesonCutArray[3] = "01631031000000"; // 
                eventCutArray[ 4] = "0000311"; photonCutArray[ 4] = "002000093270008250400000"; clusterCutArray[4] = "10000045032030000"; mesonCutArray[4] = "01631031000000"; // 
+        eventCutArray[ 5] = "0000311"; photonCutArray[ 5] = "002000093270008250400000"; clusterCutArray[5] = "10000046032030000"; mesonCutArray[5] = "01631031000000"; //
        } else if (trainConfig == 5){ // PCM variations
-               eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002000092270008250400000"; clusterCutArray[0] = "10000042032030000"; mesonCutArray[0] = "01631031000000"; // dEdx e -3, 5
-               eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002000091270008250400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01631031000000"; // dEdx e -5, 5
-               eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002000093570008250400000"; clusterCutArray[2] = "10000042032030000"; mesonCutArray[2] = "01631031000000"; // dEdx pi 2
-               eventCutArray[ 3] = "0000311"; photonCutArray[ 3] = "002000093170008250400000"; clusterCutArray[3] = "10000042032030000"; mesonCutArray[3] = "01631031000000"; // dEdx pi 0
-               eventCutArray[ 4] = "0000311"; photonCutArray[ 4] = "002000093873008250400000"; clusterCutArray[4] = "10000042032030000"; mesonCutArray[4] = "01631031000000"; // dEdx pi 2 high 1 (> 3.5 GeV)
+        eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002000092270008250400000"; clusterCutArray[0] = "10000043032030000"; mesonCutArray[0] = "01631031000000"; // dEdx e -3, 5
+        eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002000091270008250400000"; clusterCutArray[1] = "10000043032030000"; mesonCutArray[1] = "01631031000000"; // dEdx e -5, 5
+        eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002000093570008250400000"; clusterCutArray[2] = "10000043032030000"; mesonCutArray[2] = "01631031000000"; // dEdx pi 2
+        eventCutArray[ 3] = "0000311"; photonCutArray[ 3] = "002000093170008250400000"; clusterCutArray[3] = "10000043032030000"; mesonCutArray[3] = "01631031000000"; // dEdx pi 0
+        eventCutArray[ 4] = "0000311"; photonCutArray[ 4] = "002000093873008250400000"; clusterCutArray[4] = "10000043032030000"; mesonCutArray[4] = "01631031000000"; // dEdx pi 2 high 1 (> 3.5 GeV)
        } else if (trainConfig == 6){ // PCM variations
-               eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002000093270009250400000"; clusterCutArray[0] = "10000042032030000"; mesonCutArray[0] = "01631031000000"; // qt 2D 0.03
-               eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002000093270003250400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01631031000000"; // qt 1D 0.05
-               eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002000093270002250400000"; clusterCutArray[2] = "10000042032030000"; mesonCutArray[2] = "01631031000000"; // qt 1D 0.07
-               eventCutArray[ 3] = "0000311"; photonCutArray[ 3] = "002000493270008250400000"; clusterCutArray[3] = "10000042032030000"; mesonCutArray[3] = "01631031000000"; // single pt > 0.075
-               eventCutArray[ 4] = "0000311"; photonCutArray[ 4] = "002000193270008250400000"; clusterCutArray[4] = "10000042032030000"; mesonCutArray[4] = "01631031000000"; // single pt > 0.1
+        eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002000093270009250400000"; clusterCutArray[0] = "10000043032030000"; mesonCutArray[0] = "01631031000000"; // qt 2D 0.03
+        eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002000093270003250400000"; clusterCutArray[1] = "10000043032030000"; mesonCutArray[1] = "01631031000000"; // qt 1D 0.05
+        eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002000093270002250400000"; clusterCutArray[2] = "10000043032030000"; mesonCutArray[2] = "01631031000000"; // qt 1D 0.07
+        eventCutArray[ 3] = "0000311"; photonCutArray[ 3] = "002000493270008250400000"; clusterCutArray[3] = "10000043032030000"; mesonCutArray[3] = "01631031000000"; // single pt > 0.075
+        eventCutArray[ 4] = "0000311"; photonCutArray[ 4] = "002000193270008250400000"; clusterCutArray[4] = "10000043032030000"; mesonCutArray[4] = "01631031000000"; // single pt > 0.1
        } else if (trainConfig == 7){ // PCM variations
-               eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002000093270008850400000"; clusterCutArray[0] = "10000042032030000"; mesonCutArray[0] = "01631031000000"; // 2D psi pair chi2 var
-               eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002000093270008260400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01631031000000"; // 2D psi pair chi2 var
-               eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002000093270008860400000"; clusterCutArray[2] = "10000042032030000"; mesonCutArray[2] = "01631031000000"; // 2D psi pair chi2 var
-               eventCutArray[ 3] = "0000311"; photonCutArray[ 3] = "002000093270008280400000"; clusterCutArray[3] = "10000042032030000"; mesonCutArray[3] = "01631031000000"; // 2D psi pair chi2 var
-               eventCutArray[ 4] = "0000311"; photonCutArray[ 4] = "002000093270008880400000"; clusterCutArray[4] = "10000042032030000"; mesonCutArray[4] = "01631031000000"; // 2D psi pair chi2 var
+        eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002000093270008850400000"; clusterCutArray[0] = "10000043032030000"; mesonCutArray[0] = "01631031000000"; // 2D psi pair chi2 var
+        eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002000093270008260400000"; clusterCutArray[1] = "10000043032030000"; mesonCutArray[1] = "01631031000000"; // 2D psi pair chi2 var
+        eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002000093270008860400000"; clusterCutArray[2] = "10000043032030000"; mesonCutArray[2] = "01631031000000"; // 2D psi pair chi2 var
+        eventCutArray[ 3] = "0000311"; photonCutArray[ 3] = "002000093270008280400000"; clusterCutArray[3] = "10000043032030000"; mesonCutArray[3] = "01631031000000"; // 2D psi pair chi2 var
+        eventCutArray[ 4] = "0000311"; photonCutArray[ 4] = "002000093270008880400000"; clusterCutArray[4] = "10000043032030000"; mesonCutArray[4] = "01631031000000"; // 2D psi pair chi2 var
        } else if (trainConfig == 8){ // PCM variations
-               eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002000063270008250400000"; clusterCutArray[0] = "10000042032030000"; mesonCutArray[0] = "01631031000000"; // min TPC cl > 0.7
-               eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002000083270008250400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01631031000000"; // min TPC cl > 0.35
-               eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002000093270008250400000"; clusterCutArray[2] = "10000042032030000"; mesonCutArray[2] = "01631061000000"; // alpha < 0.8
-               eventCutArray[ 3] = "0000311"; photonCutArray[ 3] = "002000093270008250400000"; clusterCutArray[3] = "10000042032030000"; mesonCutArray[3] = "01631051000000"; // alpha < 0.75
+        eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002000063270008250400000"; clusterCutArray[0] = "10000043032030000"; mesonCutArray[0] = "01631031000000"; // min TPC cl > 0.7
+        eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002000083270008250400000"; clusterCutArray[1] = "10000043032030000"; mesonCutArray[1] = "01631031000000"; // min TPC cl > 0.35
+        eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002000093270008250400000"; clusterCutArray[2] = "10000043032030000"; mesonCutArray[2] = "01631061000000"; // alpha < 0.8
+        eventCutArray[ 3] = "0000311"; photonCutArray[ 3] = "002000093270008250400000"; clusterCutArray[3] = "10000043032030000"; mesonCutArray[3] = "01631051000000"; // alpha < 0.75
        } else if (trainConfig == 9){ // PCM variations
-               eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002022093270008250400000"; clusterCutArray[0] = "10000042032030000"; mesonCutArray[0] = "01631031000000"; // restrict acceptance to EMCAL loose
-               eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002044093270008250400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01631031000000"; // restrict acceptance to EMCAL tight
+        eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002022093270008250400000"; clusterCutArray[0] = "10000043032030000"; mesonCutArray[0] = "01631031000000"; // restrict acceptance to EMCAL loose
+        eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002044093270008250400000"; clusterCutArray[1] = "10000043032030000"; mesonCutArray[1] = "01631031000000"; // restrict acceptance to EMCAL tight
        } else if (trainConfig == 10){ // EMCAL clusters 2.76 TeV LHC11a, with SDD (0), kEMC1 (1)
-               eventCutArray[ 0] = "0005111"; photonCutArray[ 0] = "002000093270008250400000"; clusterCutArray[0] = "10000042062030000"; mesonCutArray[0] = "01631031000000"; // min Energy cluster = 4.5 GeV
-               eventCutArray[ 1] = "0005111"; photonCutArray[ 1] = "002000093270008250400000"; clusterCutArray[1] = "10000042072030000"; mesonCutArray[1] = "01631031000000"; // min Energy cluster = 5.0 GeV
-               eventCutArray[ 2] = "0005111"; photonCutArray[ 2] = "002000093270008250400000"; clusterCutArray[2] = "10000042082030000"; mesonCutArray[2] = "01631031000000"; // min Energy cluster = 5.5 GeV
-               eventCutArray[ 3] = "0005111"; photonCutArray[ 3] = "002000093270008250400000"; clusterCutArray[3] = "10000042092030000"; mesonCutArray[3] = "01631031000000"; // min Energy cluster = 6.0 GeV
+        eventCutArray[ 0] = "0005111"; photonCutArray[ 0] = "002000093270008250400000"; clusterCutArray[0] = "10000043062030000"; mesonCutArray[0] = "01631031000000"; // min Energy cluster = 4.5 GeV
+        eventCutArray[ 1] = "0005111"; photonCutArray[ 1] = "002000093270008250400000"; clusterCutArray[1] = "10000043072030000"; mesonCutArray[1] = "01631031000000"; // min Energy cluster = 5.0 GeV
+        eventCutArray[ 2] = "0005111"; photonCutArray[ 2] = "002000093270008250400000"; clusterCutArray[2] = "10000043082030000"; mesonCutArray[2] = "01631031000000"; // min Energy cluster = 5.5 GeV
+        eventCutArray[ 3] = "0005111"; photonCutArray[ 3] = "002000093270008250400000"; clusterCutArray[3] = "10000043092030000"; mesonCutArray[3] = "01631031000000"; // min Energy cluster = 6.0 GeV
                // LHC13g       
        } else if (trainConfig == 11){  // EMCAL clusters, EMCEGA triggers, track matching 0.035
-               eventCutArray[ 0] = "0008311"; photonCutArray[ 0] = "002000093270008250400000"; clusterCutArray[0] = "10000042032030000"; mesonCutArray[0] = "01631031000000"; // EMCEG1, 
-               eventCutArray[ 1] = "0008511"; photonCutArray[ 1] = "002000093270008250400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01631031000000"; // EMCEG2, 
-               eventCutArray[ 2] = "0009311"; photonCutArray[ 2] = "002000093270008250400000"; clusterCutArray[2] = "10000042032030000"; mesonCutArray[2] = "01631031000000"; // EMCEJ1, 
-               eventCutArray[ 3] = "0009511"; photonCutArray[ 3] = "002000093270008250400000"; clusterCutArray[3] = "10000042032030000"; mesonCutArray[3] = "01631031000000"; // EMCEJ2, 
-               eventCutArray[ 4] = "0000011"; photonCutArray[ 4] = "002000093270008250400000"; clusterCutArray[4] = "10000042032030000"; mesonCutArray[4] = "01631031000000"; // INT7
-               eventCutArray[ 5] = "0005211"; photonCutArray[ 5] = "002000093270008250400000"; clusterCutArray[5] = "10000042032030000"; mesonCutArray[5] = "01631031000000"; // EMC7
+        eventCutArray[ 0] = "0008311"; photonCutArray[ 0] = "002000093270008250400000"; clusterCutArray[0] = "10000043032030000"; mesonCutArray[0] = "01631031000000"; // EMCEG1,
+        eventCutArray[ 1] = "0008511"; photonCutArray[ 1] = "002000093270008250400000"; clusterCutArray[1] = "10000043032030000"; mesonCutArray[1] = "01631031000000"; // EMCEG2,
+        eventCutArray[ 2] = "0009311"; photonCutArray[ 2] = "002000093270008250400000"; clusterCutArray[2] = "10000043032030000"; mesonCutArray[2] = "01631031000000"; // EMCEJ1,
+        eventCutArray[ 3] = "0009511"; photonCutArray[ 3] = "002000093270008250400000"; clusterCutArray[3] = "10000043032030000"; mesonCutArray[3] = "01631031000000"; // EMCEJ2,
+        eventCutArray[ 4] = "0000011"; photonCutArray[ 4] = "002000093270008250400000"; clusterCutArray[4] = "10000043032030000"; mesonCutArray[4] = "01631031000000"; // INT7
+        eventCutArray[ 5] = "0005211"; photonCutArray[ 5] = "002000093270008250400000"; clusterCutArray[5] = "10000043032030000"; mesonCutArray[5] = "01631031000000"; // EMC7
                
        // ************************************* PHOS cuts ****************************************************
        // LHC11a       
        } else if (trainConfig == 31) { //PHOS clusters
-               eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002000093270008250400000"; clusterCutArray[0] = "20000041033200000"; mesonCutArray[0] = "01631031000000"; 
-               eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002000093270008250400000"; clusterCutArray[1] = "20000042033200000"; mesonCutArray[1] = "01631031000000"; 
-               eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002000093270008250400000"; clusterCutArray[2] = "20000043033200000"; mesonCutArray[2] = "01631031000000"; 
-               eventCutArray[ 3] = "0006111"; photonCutArray[ 3] = "002000093270008250400000"; clusterCutArray[3] = "20000041033200000"; mesonCutArray[3] = "01631031000000"; 
-               eventCutArray[ 4] = "0006111"; photonCutArray[ 4] = "002000093270008250400000"; clusterCutArray[4] = "20000042033200000"; mesonCutArray[4] = "01631031000000"; 
-               eventCutArray[ 5] = "0006111"; photonCutArray[ 5] = "002000093270008250400000"; clusterCutArray[5] = "20000043033200000"; mesonCutArray[5] = "01631031000000"; 
+        eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002000093270008250400000"; clusterCutArray[0] = "20000047033200000"; mesonCutArray[0] = "01631031000000";
+        eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002000093270008250400000"; clusterCutArray[1] = "20000048033200000"; mesonCutArray[1] = "01631031000000";
+        eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002000093270008250400000"; clusterCutArray[2] = "20000049033200000"; mesonCutArray[2] = "01631031000000";
+        eventCutArray[ 3] = "0006111"; photonCutArray[ 3] = "002000093270008250400000"; clusterCutArray[3] = "20000047033200000"; mesonCutArray[3] = "01631031000000";
+        eventCutArray[ 4] = "0006111"; photonCutArray[ 4] = "002000093270008250400000"; clusterCutArray[4] = "20000048033200000"; mesonCutArray[4] = "01631031000000";
+        eventCutArray[ 5] = "0006111"; photonCutArray[ 5] = "002000093270008250400000"; clusterCutArray[5] = "20000049033200000"; mesonCutArray[5] = "01631031000000";
        // LHC13g & LHC12x
        } else if (trainConfig == 32) { //PHOS clusters
-               eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002000093270008250400000"; clusterCutArray[0] = "20000041033200000"; mesonCutArray[0] = "01631031000000"; 
-               eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002000093270008250400000"; clusterCutArray[1] = "20000042033200000"; mesonCutArray[1] = "01631031000000"; 
-               eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002000093270008250400000"; clusterCutArray[2] = "20000043033200000"; mesonCutArray[2] = "01631031000000"; 
-               eventCutArray[ 3] = "0006211"; photonCutArray[ 3] = "002000093270008250400000"; clusterCutArray[3] = "20000041033200000"; mesonCutArray[3] = "01631031000000"; 
-               eventCutArray[ 4] = "0006211"; photonCutArray[ 4] = "002000093270008250400000"; clusterCutArray[4] = "20000042033200000"; mesonCutArray[4] = "01631031000000"; 
-               eventCutArray[ 5] = "0006211"; photonCutArray[ 5] = "002000093270008250400000"; clusterCutArray[5] = "20000043033200000"; mesonCutArray[5] = "01631031000000"; 
+        eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002000093270008250400000"; clusterCutArray[0] = "20000047033200000"; mesonCutArray[0] = "01631031000000";
+        eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002000093270008250400000"; clusterCutArray[1] = "20000048033200000"; mesonCutArray[1] = "01631031000000";
+        eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002000093270008250400000"; clusterCutArray[2] = "20000049033200000"; mesonCutArray[2] = "01631031000000";
+        eventCutArray[ 3] = "0006211"; photonCutArray[ 3] = "002000093270008250400000"; clusterCutArray[3] = "20000047033200000"; mesonCutArray[3] = "01631031000000";
+        eventCutArray[ 4] = "0006211"; photonCutArray[ 4] = "002000093270008250400000"; clusterCutArray[4] = "20000048033200000"; mesonCutArray[4] = "01631031000000";
+        eventCutArray[ 5] = "0006211"; photonCutArray[ 5] = "002000093270008250400000"; clusterCutArray[5] = "20000049033200000"; mesonCutArray[5] = "01631031000000";
        } else {
                Error(Form("GammaConvCalo_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
                return;