- added cuts for selection of photon quality for all systems
authorfbock <friederike.bock@cern.ch>
Thu, 17 Apr 2014 12:46:13 +0000 (14:46 +0200)
committerfbock <friederike.bock@cern.ch>
Thu, 17 Apr 2014 12:46:13 +0000 (14:46 +0200)
- changed minR cut variable  7 from minR = 26 to 35

PWGGA/GammaConv/AliConversionCuts.cxx
PWGGA/GammaConv/macros/AddTask_GammaConvV1_PbPb.C
PWGGA/GammaConv/macros/AddTask_GammaConvV1_pPb.C
PWGGA/GammaConv/macros/AddTask_GammaConvV1_pp.C

index 2d2d975..a8118d0 100644 (file)
@@ -2434,7 +2434,6 @@ Bool_t AliConversionCuts::SetRCut(Int_t RCut){
       fMaxR = 70.;
       fMinR = 5.;
       break;
-      // High purity cuts for PbPb (remove first layers of material)
    case 5:
       fMaxR = 180.;
       fMinR = 10.;
@@ -2445,7 +2444,7 @@ Bool_t AliConversionCuts::SetRCut(Int_t RCut){
       break;
    case 7:
       fMaxR = 180.;
-      fMinR = 26.;
+      fMinR = 35.; //old 26.
       break;
    case 8:
       fMaxR = 180.;
index 370370b..e22ac51 100644 (file)
@@ -170,174 +170,174 @@ void AddTask_GammaConvV1_PbPb(  Int_t trainConfig = 1,  //change different set o
       cutarray[ 2] = "5680002048092970023220000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
       cutarray[ 3] = "5480002048092970023220000000"; mesonCutArray[ 3] = "01522065009000"; // 20-40% 
       cutarray[ 4] = "5490002048092970023220000000"; mesonCutArray[ 4] = "01522065009000"; // 40-90%      
-   } else  if (trainConfig == 13){ // eta 0.65 (new standard), y = 0.6 (new Standard)
+   } else  if (trainConfig == 13){ // LHC10h standard, eta 0.65, y = 0.6 
       cutarray[ 0] = "6010001032092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
       cutarray[ 1] = "6120001032092970023220000000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
       cutarray[ 2] = "5010001032092970023220000000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
       cutarray[ 3] = "5120001032092970023220000000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
       cutarray[ 4] = "5020001032092970023220000000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
-   } else if (trainConfig == 14) {  // eta 0.65 (new standard), y = 0.6 (new Standard)
+   } else if (trainConfig == 14) {  // LHC10h standard, eta 0.65, y = 0.6 
       cutarray[ 0] = "5240001032092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
       cutarray[ 1] = "5460001032092970023220000000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
       cutarray[ 2] = "5680001032092970023220000000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
       cutarray[ 3] = "5480001032092970023220000000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
       cutarray[ 4] = "5490001032092970023220000000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
-   } else if (trainConfig == 15) { // eta 0.65 (new standard), y = 0.6 (new Standard) cuts only added signals
+   } else if (trainConfig == 15) { // LHC10h standard, eta 0.65, y = 0.6  added signals
       cutarray[ 0] = "6010002032092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
       cutarray[ 1] = "6120002032092970023220000000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
       cutarray[ 2] = "5010002032092970023220000000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
       cutarray[ 3] = "5120002032092970023220000000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
       cutarray[ 4] = "5020002032092970023220000000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
-   } else if (trainConfig == 16) { // eta 0.65 (new standard), y = 0.6 (new Standard) cuts only added signals
+   } else if (trainConfig == 16) { // LHC10h standard, eta 0.65, y = 0.6  added signals
       cutarray[ 0] = "5240002032092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
       cutarray[ 1] = "5460002032092970023220000000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
       cutarray[ 2] = "5680002032092970023220000000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
       cutarray[ 3] = "5480002032092970023220000000"; mesonCutArray[ 3] = "01523065009000"; // 20-40% 
       cutarray[ 4] = "5490002032092970023220000000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
-   } else  if (trainConfig == 17){ // eta 0.6, y = 0.6 (new Standard)
-      cutarray[ 0] = "6010001012092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
-      cutarray[ 1] = "6120001012092970023220000000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
-      cutarray[ 2] = "5010001012092970023220000000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
-      cutarray[ 3] = "5120001012092970023220000000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
-      cutarray[ 4] = "5020001012092970023220000000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
-   } else if (trainConfig == 18) {  // eta 0.6, y = 0.6 (new Standard)
-      cutarray[ 0] = "5240001012092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
-      cutarray[ 1] = "5460001012092970023220000000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
-      cutarray[ 2] = "5680001012092970023220000000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
-      cutarray[ 3] = "5480001012092970023220000000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
-      cutarray[ 4] = "5490001012092970023220000000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
-   } else if (trainConfig == 19) { // eta 0.6, y = 0.6 (new Standard) cuts only added signals
-      cutarray[ 0] = "6010002012092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
-      cutarray[ 1] = "6120002012092970023220000000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
-      cutarray[ 2] = "5010002012092970023220000000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
-      cutarray[ 3] = "5120002012092970023220000000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
-      cutarray[ 4] = "5020002012092970023220000000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
-   } else if (trainConfig == 20) { // eta 0.6, y = 0.6 (new Standard) cuts only added signals
-      cutarray[ 0] = "5240002012092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
-      cutarray[ 1] = "5460002012092970023220000000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
-      cutarray[ 2] = "5680002012092970023220000000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
-      cutarray[ 3] = "5480002012092970023220000000"; mesonCutArray[ 3] = "01523065009000"; // 20-40% 
-      cutarray[ 4] = "5490002012092970023220000000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
-   } else  if (trainConfig == 21){ // eta 0.7, y = 0.6 (new Standard)
-      cutarray[ 0] = "6010001072092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
-      cutarray[ 1] = "6120001072092970023220000000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
-      cutarray[ 2] = "5010001072092970023220000000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
-      cutarray[ 3] = "5120001072092970023220000000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
-      cutarray[ 4] = "5020001072092970023220000000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
-   } else if (trainConfig == 22) {  // eta 0.7, y = 0.6 (new Standard)
-      cutarray[ 0] = "5240001072092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
-      cutarray[ 1] = "5460001072092970023220000000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
-      cutarray[ 2] = "5680001072092970023220000000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
-      cutarray[ 3] = "5480001072092970023220000000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
-      cutarray[ 4] = "5490001072092970023220000000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
-   } else if (trainConfig == 23) { // eta 0.7, y = 0.6 (new Standard) cuts only added signals
-      cutarray[ 0] = "6010002072092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
-      cutarray[ 1] = "6120002072092970023220000000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
-      cutarray[ 2] = "5010002072092970023220000000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
-      cutarray[ 3] = "5120002072092970023220000000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
-      cutarray[ 4] = "5020002072092970023220000000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
-   } else if (trainConfig == 24) { // eta 0.7, y = 0.6 (new Standard) cuts only added signals
-      cutarray[ 0] = "5240002072092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
-      cutarray[ 1] = "5460002072092970023220000000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
-      cutarray[ 2] = "5680002072092970023220000000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
-      cutarray[ 3] = "5480002072092970023220000000"; mesonCutArray[ 3] = "01523065009000"; // 20-40% 
-      cutarray[ 4] = "5490002072092970023220000000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
-   } else  if (trainConfig == 25){ // eta 0.5, y = 0.6 (new Standard)
-      cutarray[ 0] = "6010001052092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
-      cutarray[ 1] = "6120001052092970023220000000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
-      cutarray[ 2] = "5010001052092970023220000000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
-      cutarray[ 3] = "5120001052092970023220000000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
-      cutarray[ 4] = "5020001052092970023220000000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
-   } else if (trainConfig == 26) {  // eta 0.5, y = 0.6 (new Standard)
-      cutarray[ 0] = "5240001052092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
-      cutarray[ 1] = "5460001052092970023220000000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
-      cutarray[ 2] = "5680001052092970023220000000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
-      cutarray[ 3] = "5480001052092970023220000000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
-      cutarray[ 4] = "5490001052092970023220000000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
-   } else if (trainConfig == 27) { // eta 0.5, y = 0.6 (new Standard) cuts only added signals
-      cutarray[ 0] = "6010002052092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
-      cutarray[ 1] = "6120002052092970023220000000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
-      cutarray[ 2] = "5010002052092970023220000000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
-      cutarray[ 3] = "5120002052092970023220000000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
-      cutarray[ 4] = "5020002052092970023220000000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
-   } else if (trainConfig == 28) { // eta 0.5, y = 0.6 (new Standard) cuts only added signals
-      cutarray[ 0] = "5240002052092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
-      cutarray[ 1] = "5460002052092970023220000000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
-      cutarray[ 2] = "5680002052092970023220000000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
-      cutarray[ 3] = "5480002052092970023220000000"; mesonCutArray[ 3] = "01523065009000"; // 20-40% 
-      cutarray[ 4] = "5490002052092970023220000000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
-   } else  if (trainConfig == 29){ // eta 0.65 (new standard), y = 0.6 (new Standard) pt dependent alpha
-      cutarray[ 0] = "6010001032092970023220000000"; mesonCutArray[ 0] = "01523015009000"; // 0-5%
-      cutarray[ 1] = "6120001032092970023220000000"; mesonCutArray[ 1] = "01523015009000"; // 5-10%
-      cutarray[ 2] = "5010001032092970023220000000"; mesonCutArray[ 2] = "01523015009000"; // 0-10%
-      cutarray[ 3] = "5120001032092970023220000000"; mesonCutArray[ 3] = "01523015009000"; // 10-20%
-      cutarray[ 4] = "5020001032092970023220000000"; mesonCutArray[ 4] = "01523015009000"; // 0-20%
-   } else if (trainConfig == 30) {  // eta 0.65 (new standard), y = 0.6 (new Standard) pt dependent alpha
-      cutarray[ 0] = "5240001032092970023220000000"; mesonCutArray[ 0] = "01523015009000"; // 20-40%
-      cutarray[ 1] = "5460001032092970023220000000"; mesonCutArray[ 1] = "01523025009000"; // 40-60%
-      cutarray[ 2] = "5680001032092970023220000000"; mesonCutArray[ 2] = "01523025009000"; // 60-80%
-      cutarray[ 3] = "5480001032092970023220000000"; mesonCutArray[ 3] = "01523025009000"; // 40-80%
-      cutarray[ 4] = "5490001032092970023220000000"; mesonCutArray[ 4] = "01523025009000"; // 40-90%
-   } else if (trainConfig == 31) { // eta 0.65 (new standard), y = 0.6 (new Standard) cuts only added signals, pt dependent alpha
-      cutarray[ 0] = "6010002032092970023220000000"; mesonCutArray[ 0] = "01523015009000"; // 0-5%
-      cutarray[ 1] = "6120002032092970023220000000"; mesonCutArray[ 1] = "01523015009000"; // 5-10%
-      cutarray[ 2] = "5010002032092970023220000000"; mesonCutArray[ 2] = "01523015009000"; // 0-10%
-      cutarray[ 3] = "5120002032092970023220000000"; mesonCutArray[ 3] = "01523015009000"; // 10-20%
-      cutarray[ 4] = "5020002032092970023220000000"; mesonCutArray[ 4] = "01523015009000"; // 0-20%
-   } else if (trainConfig == 32) { // eta 0.65 (new standard), y = 0.6 (new Standard) cuts only added signals, pt dependent alpha
-      cutarray[ 0] = "5240002032092970023220000000"; mesonCutArray[ 0] = "01523015009000"; // 20-40%
-      cutarray[ 1] = "5460002032092970023220000000"; mesonCutArray[ 1] = "01523025009000"; // 40-60%
-      cutarray[ 2] = "5680002032092970023220000000"; mesonCutArray[ 2] = "01523025009000"; // 60-80%
-      cutarray[ 3] = "5480002032092970023220000000"; mesonCutArray[ 3] = "01523025009000"; // 20-40% 
-      cutarray[ 4] = "5490002032092970023220000000"; mesonCutArray[ 4] = "01523025009000"; // 40-90%
-   } else if (trainConfig == 33){ // Standard cuts, eta 0.9, only to be run on data
-      cutarray[ 0] = "6010001002092970023220000000"; mesonCutArray[ 0] = "01525045009000"; // 0-5%
-      cutarray[ 1] = "6120001002092970023220000000"; mesonCutArray[ 1] = "01525045009000"; // 5-10%
-      cutarray[ 2] = "5010001002092970023220000000"; mesonCutArray[ 2] = "01525045009000"; // 0-10%
-      cutarray[ 3] = "5120001002092970023220000000"; mesonCutArray[ 3] = "01525045009000"; // 10-20%
-      cutarray[ 4] = "5020001002092970023220000000"; mesonCutArray[ 4] = "01525045009000"; // 0-20%
-   } else if (trainConfig == 34) { // Standard cuts, eta 0.9, only to be run on data
-      cutarray[ 0] = "5240001002092970023220000000"; mesonCutArray[ 0] = "01525045009000"; // 20-40%
-      cutarray[ 1] = "5460001002092970023220000000"; mesonCutArray[ 1] = "01525065009000"; // 40-60%
-      cutarray[ 2] = "5680001002092970023220000000"; mesonCutArray[ 2] = "01525065009000"; // 60-80%
-      cutarray[ 3] = "5480001002092970023220000000"; mesonCutArray[ 3] = "01525065009000"; // 40-80%
-      cutarray[ 4] = "5350001002092970023220000000"; mesonCutArray[ 4] = "01525065009000"; // 40-90%   
-   } else if (trainConfig == 35){ // Standard cuts, eta 1.4, only to be run on data
-      cutarray[ 0] = "6010001022092970023220000000"; mesonCutArray[ 0] = "01520045009000"; // 0-5%
-      cutarray[ 1] = "6120001022092970023220000000"; mesonCutArray[ 1] = "01520045009000"; // 5-10%
-      cutarray[ 2] = "5010001022092970023220000000"; mesonCutArray[ 2] = "01520045009000"; // 0-10%
-      cutarray[ 3] = "5120001022092970023220000000"; mesonCutArray[ 3] = "01520045009000"; // 10-20%
-      cutarray[ 4] = "5020001022092970023220000000"; mesonCutArray[ 4] = "01520045009000"; // 0-20%
-   } else if (trainConfig == 36) { // Standard cuts, eta 1.4, only to be run on data
-      cutarray[ 0] = "5240001022092970023220000000"; mesonCutArray[ 0] = "01520045009000"; // 20-40%
-      cutarray[ 1] = "5460001022092970023220000000"; mesonCutArray[ 1] = "01520065009000"; // 40-60%
-      cutarray[ 2] = "5680001022092970023220000000"; mesonCutArray[ 2] = "01520065009000"; // 60-80%
-      cutarray[ 3] = "5480001022092970023220000000"; mesonCutArray[ 3] = "01520065009000"; // 40-80%
-      cutarray[ 4] = "5350001022092970023220000000"; mesonCutArray[ 4] = "01520065009000"; // 40-90%   
-   } else if (trainConfig == 37){ // Standard cuts, eta 0.9, only to be run on data : kMB
-      cutarray[ 0] = "6014001002092970023220000000"; mesonCutArray[ 0] = "01525045009000"; // 0-5%
-      cutarray[ 1] = "6124001002092970023220000000"; mesonCutArray[ 1] = "01525045009000"; // 5-10%
-      cutarray[ 2] = "5014001002092970023220000000"; mesonCutArray[ 2] = "01525045009000"; // 0-10%
-      cutarray[ 3] = "5124001002092970023220000000"; mesonCutArray[ 3] = "01525045009000"; // 10-20%
-      cutarray[ 4] = "5024001002092970023220000000"; mesonCutArray[ 4] = "01525045009000"; // 0-20%
-   } else if (trainConfig == 38) { // Standard cuts, eta 0.9, only to be run on data : kMB
-      cutarray[ 0] = "5244001002092970023220000000"; mesonCutArray[ 0] = "01525045009000"; // 20-40%
-      cutarray[ 1] = "5464001002092970023220000000"; mesonCutArray[ 1] = "01525065009000"; // 40-60%
-      cutarray[ 2] = "5684001002092970023220000000"; mesonCutArray[ 2] = "01525065009000"; // 60-80%
-      cutarray[ 3] = "5484001002092970023220000000"; mesonCutArray[ 3] = "01525065009000"; // 40-80%
-      cutarray[ 4] = "5354001002092970023220000000"; mesonCutArray[ 4] = "01525065009000"; // 40-90%   
-   } else if (trainConfig == 39){ // Standard cuts, eta 0.9, only to be run on data :kCentral
-      cutarray[ 0] = "6015001002092970023220000000"; mesonCutArray[ 0] = "01525045009000"; // 0-5%
-      cutarray[ 1] = "6125001002092970023220000000"; mesonCutArray[ 1] = "01525045009000"; // 5-10%
-      cutarray[ 2] = "5015001002092970023220000000"; mesonCutArray[ 2] = "01525045009000"; // 0-10%
-      cutarray[ 3] = "5125001002092970023220000000"; mesonCutArray[ 3] = "01525045009000"; // 10-20%
-      cutarray[ 4] = "5025001002092970023220000000"; mesonCutArray[ 4] = "01525045009000"; // 0-20%
-   } else if (trainConfig == 40) { // Standard cuts, eta 0.9, only to be run on data : kCentral
-      cutarray[ 0] = "5245001002092970023220000000"; mesonCutArray[ 0] = "01525045009000"; // 20-40%
-      cutarray[ 1] = "5465001002092970023220000000"; mesonCutArray[ 1] = "01525065009000"; // 40-60%
-      cutarray[ 2] = "5685001002092970023220000000"; mesonCutArray[ 2] = "01525065009000"; // 60-80%
-      cutarray[ 3] = "5485001002092970023220000000"; mesonCutArray[ 3] = "01525065009000"; // 40-80%
-      cutarray[ 4] = "5355001002092970023220000000"; mesonCutArray[ 4] = "01525065009000"; // 40-90%   
+   } else  if (trainConfig == 17){ // LHC10h standard, eta 0.65, y = 0.6, photon quality 1 
+      cutarray[ 0] = "6010001032092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
+      cutarray[ 1] = "6120001032092970023220020000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
+      cutarray[ 2] = "5010001032092970023220020000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
+      cutarray[ 3] = "5120001032092970023220020000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
+      cutarray[ 4] = "5020001032092970023220020000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
+   } else if (trainConfig == 18) {  // LHC10h standard, eta 0.65, y = 0.6, photon quality 1 
+      cutarray[ 0] = "5240001032092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
+      cutarray[ 1] = "5460001032092970023220020000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
+      cutarray[ 2] = "5680001032092970023220020000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
+      cutarray[ 3] = "5480001032092970023220020000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
+      cutarray[ 4] = "5490001032092970023220020000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
+   } else if (trainConfig == 19) { // LHC10h standard, eta 0.65, y = 0.6, photon quality 1  added signal
+      cutarray[ 0] = "6010002032092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
+      cutarray[ 1] = "6120002032092970023220020000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
+      cutarray[ 2] = "5010002032092970023220020000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
+      cutarray[ 3] = "5120002032092970023220020000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
+      cutarray[ 4] = "5020002032092970023220020000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
+   } else if (trainConfig == 20) { // LHC10h standard, eta 0.65, y = 0.6, photon quality 1 added signal
+      cutarray[ 0] = "5240002032092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
+      cutarray[ 1] = "5460002032092970023220020000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
+      cutarray[ 2] = "5680002032092970023220020000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
+      cutarray[ 3] = "5480002032092970023220020000"; mesonCutArray[ 3] = "01523065009000"; // 20-40% 
+      cutarray[ 4] = "5490002032092970023220020000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
+   } else  if (trainConfig == 21){ // LHC10h standard, eta 0.65, y = 0.6, photon quality 2 
+      cutarray[ 0] = "6010001032092970023220030000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
+      cutarray[ 1] = "6120001032092970023220030000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
+      cutarray[ 2] = "5010001032092970023220030000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
+      cutarray[ 3] = "5120001032092970023220030000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
+      cutarray[ 4] = "5020001032092970023220030000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
+   } else if (trainConfig == 22) {  // LHC10h standard, eta 0.65, y = 0.6, photon quality 2 
+      cutarray[ 0] = "5240001032092970023220030000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
+      cutarray[ 1] = "5460001032092970023220030000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
+      cutarray[ 2] = "5680001032092970023220030000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
+      cutarray[ 3] = "5480001032092970023220030000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
+      cutarray[ 4] = "5490001032092970023220030000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
+   } else if (trainConfig == 23) { // LHC10h standard, eta 0.65, y = 0.6, photon quality 2  added signal
+      cutarray[ 0] = "6010002032092970023220030000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
+      cutarray[ 1] = "6120002032092970023220030000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
+      cutarray[ 2] = "5010002032092970023220030000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
+      cutarray[ 3] = "5120002032092970023220030000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
+      cutarray[ 4] = "5020002032092970023220030000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
+   } else if (trainConfig == 24) { // LHC10h standard, eta 0.65, y = 0.6, photon quality 2 added signal
+      cutarray[ 0] = "5240002032092970023220030000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
+      cutarray[ 1] = "5460002032092970023220030000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
+      cutarray[ 2] = "5680002032092970023220030000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
+      cutarray[ 3] = "5480002032092970023220030000"; mesonCutArray[ 3] = "01523065009000"; // 20-40% 
+      cutarray[ 4] = "5490002032092970023220030000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
+   } else  if (trainConfig == 25){ // LHC10h standard, eta 0.65, y = 0.6, photon quality 3 
+      cutarray[ 0] = "6010001032092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
+      cutarray[ 1] = "6120001032092970023220040000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
+      cutarray[ 2] = "5010001032092970023220040000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
+      cutarray[ 3] = "5120001032092970023220040000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
+      cutarray[ 4] = "5020001032092970023220040000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
+   } else if (trainConfig == 26) {  // LHC10h standard, eta 0.65, y = 0.6, photon quality 3 
+      cutarray[ 0] = "5240001032092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
+      cutarray[ 1] = "5460001032092970023220040000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
+      cutarray[ 2] = "5680001032092970023220040000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
+      cutarray[ 3] = "5480001032092970023220040000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
+      cutarray[ 4] = "5490001032092970023220040000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
+   } else if (trainConfig == 27) { // LHC10h standard, eta 0.65, y = 0.6, photon quality 3  added signal
+      cutarray[ 0] = "6010002032092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
+      cutarray[ 1] = "6120002032092970023220040000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
+      cutarray[ 2] = "5010002032092970023220040000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
+      cutarray[ 3] = "5120002032092970023220040000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
+      cutarray[ 4] = "5020002032092970023220040000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
+   } else if (trainConfig == 28) { // LHC10h standard, eta 0.65, y = 0.6, photon quality 3 added signal
+      cutarray[ 0] = "5240002032092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
+      cutarray[ 1] = "5460002032092970023220040000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
+      cutarray[ 2] = "5680002032092970023220040000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
+      cutarray[ 3] = "5480002032092970023220040000"; mesonCutArray[ 3] = "01523065009000"; // 20-40% 
+      cutarray[ 4] = "5490002032092970023220040000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
+   } else  if (trainConfig == 29){ // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm 
+      cutarray[ 0] = "6010001037092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
+      cutarray[ 1] = "6120001037092970023220000000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
+      cutarray[ 2] = "5010001037092970023220000000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
+      cutarray[ 3] = "5120001037092970023220000000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
+      cutarray[ 4] = "5020001037092970023220000000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
+   } else if (trainConfig == 30) {  // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm 
+      cutarray[ 0] = "5240001037092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
+      cutarray[ 1] = "5460001037092970023220000000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
+      cutarray[ 2] = "5680001037092970023220000000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
+      cutarray[ 3] = "5480001037092970023220000000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
+      cutarray[ 4] = "5490001037092970023220000000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
+   } else if (trainConfig == 31) { // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm  added signals
+      cutarray[ 0] = "6010002037092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
+      cutarray[ 1] = "6120002037092970023220000000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
+      cutarray[ 2] = "5010002037092970023220000000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
+      cutarray[ 3] = "5120002037092970023220000000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
+      cutarray[ 4] = "5020002037092970023220000000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
+   } else if (trainConfig == 32) { // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm  added signals
+      cutarray[ 0] = "5240002037092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
+      cutarray[ 1] = "5460002037092970023220000000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
+      cutarray[ 2] = "5680002037092970023220000000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
+      cutarray[ 3] = "5480002037092970023220000000"; mesonCutArray[ 3] = "01523065009000"; // 20-40% 
+      cutarray[ 4] = "5490002037092970023220000000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
+   } else  if (trainConfig == 33){ // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 1 
+      cutarray[ 0] = "6010001037092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
+      cutarray[ 1] = "6120001037092970023220020000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
+      cutarray[ 2] = "5010001037092970023220020000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
+      cutarray[ 3] = "5120001037092970023220020000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
+      cutarray[ 4] = "5020001037092970023220020000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
+   } else if (trainConfig == 34) {  // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 1  
+      cutarray[ 0] = "5240001037092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
+      cutarray[ 1] = "5460001037092970023220020000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
+      cutarray[ 2] = "5680001037092970023220020000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
+      cutarray[ 3] = "5480001037092970023220020000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
+      cutarray[ 4] = "5490001037092970023220020000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
+   } else if (trainConfig == 35) { // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 1  added signals
+      cutarray[ 0] = "6010002037092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
+      cutarray[ 1] = "6120002037092970023220020000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
+      cutarray[ 2] = "5010002037092970023220020000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
+      cutarray[ 3] = "5120002037092970023220020000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
+      cutarray[ 4] = "5020002037092970023220020000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
+   } else if (trainConfig == 36) { // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 1  added signals
+      cutarray[ 0] = "5240002037092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
+      cutarray[ 1] = "5460002037092970023220020000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
+      cutarray[ 2] = "5680002037092970023220020000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
+      cutarray[ 3] = "5480002037092970023220020000"; mesonCutArray[ 3] = "01523065009000"; // 20-40% 
+      cutarray[ 4] = "5490002037092970023220020000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
+   } else  if (trainConfig == 37){ // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 3 
+      cutarray[ 0] = "6010001037092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
+      cutarray[ 1] = "6120001037092970023220040000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
+      cutarray[ 2] = "5010001037092970023220040000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
+      cutarray[ 3] = "5120001037092970023220040000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
+      cutarray[ 4] = "5020001037092970023220040000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
+   } else if (trainConfig == 38) {  // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 3  
+      cutarray[ 0] = "5240001037092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
+      cutarray[ 1] = "5460001037092970023220040000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
+      cutarray[ 2] = "5680001037092970023220040000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
+      cutarray[ 3] = "5480001037092970023220040000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
+      cutarray[ 4] = "5490001037092970023220040000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
+   } else if (trainConfig == 39) { // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 3  added signals
+      cutarray[ 0] = "6010002037092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
+      cutarray[ 1] = "6120002037092970023220040000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
+      cutarray[ 2] = "5010002037092970023220040000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
+      cutarray[ 3] = "5120002037092970023220040000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
+      cutarray[ 4] = "5020002037092970023220040000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
+   } else if (trainConfig == 40) { // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 3  added signals
+      cutarray[ 0] = "5240002037092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
+      cutarray[ 1] = "5460002037092970023220040000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
+      cutarray[ 2] = "5680002037092970023220040000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
+      cutarray[ 3] = "5480002037092970023220040000"; mesonCutArray[ 3] = "01523065009000"; // 20-40% 
+      cutarray[ 4] = "5490002037092970023220040000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
    } else if (trainConfig == 41){ // Standard cuts, eta 0.9, only to be run on data :kSemiCentral
       cutarray[ 0] = "6016001002092970023220000000"; mesonCutArray[ 0] = "01525045009000"; // 0-5%
       cutarray[ 1] = "6126001002092970023220000000"; mesonCutArray[ 1] = "01525045009000"; // 5-10%
@@ -463,7 +463,7 @@ void AddTask_GammaConvV1_PbPb(  Int_t trainConfig = 1,  //change different set o
       cutarray[ 1] = "5340002002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 30-40% 
       cutarray[ 2] = "5450002002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 40-50% 
       cutarray[ 3] = "5560002002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 50-60% 
-      cutarray[ 4] = "5670002002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 60-70%                   
+      cutarray[ 4] = "5250002002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 60-70%                   
    } else if ( trainConfig == 62){ // cleaner cuts
       cutarray[ 0] = "6230001002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
       cutarray[ 1] = "6340001002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
@@ -764,6 +764,78 @@ void AddTask_GammaConvV1_PbPb(  Int_t trainConfig = 1,  //change different set o
                cutarray[ 2] = "5010002002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
                cutarray[ 3] = "5240002002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
                cutarray[ 4] = "5250002002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%                   
+   } else if ( trainConfig == 112){ // cleaner cuts photon Quality 1
+               cutarray[ 0] = "6010001002092970028250420000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
+               cutarray[ 1] = "6120001002092970028250420000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
+               cutarray[ 2] = "5010001002092970028250420000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
+               cutarray[ 3] = "5240001002092970028250420000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
+               cutarray[ 4] = "5250001002092970028250420000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
+   } else if ( trainConfig == 113){ // cleaner cuts added signal photon Quality 1
+               cutarray[ 0] = "6010002002092970028250420000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
+               cutarray[ 1] = "6120002002092970028250420000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
+               cutarray[ 2] = "5010002002092970028250420000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
+               cutarray[ 3] = "5240002002092970028250420000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
+               cutarray[ 4] = "5250002002092970028250420000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%                   
+   } else if ( trainConfig == 114){ // cleaner cuts photon Quality 2
+               cutarray[ 0] = "6010001002092970028250430000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
+               cutarray[ 1] = "6120001002092970028250430000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
+               cutarray[ 2] = "5010001002092970028250430000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
+               cutarray[ 3] = "5240001002092970028250430000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
+               cutarray[ 4] = "5250001002092970028250430000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
+   } else if ( trainConfig == 115){ // cleaner cuts added signal photon Quality 2
+               cutarray[ 0] = "6010002002092970028250430000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
+               cutarray[ 1] = "6120002002092970028250430000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
+               cutarray[ 2] = "5010002002092970028250430000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
+               cutarray[ 3] = "5240002002092970028250430000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
+               cutarray[ 4] = "5250002002092970028250430000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%                   
+   } else if ( trainConfig == 116){ // cleaner cuts photon Quality 3
+               cutarray[ 0] = "6010001002092970028250440000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
+               cutarray[ 1] = "6120001002092970028250440000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
+               cutarray[ 2] = "5010001002092970028250440000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
+               cutarray[ 3] = "5240001002092970028250440000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
+               cutarray[ 4] = "5250001002092970028250440000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
+   } else if ( trainConfig == 117){ // cleaner cuts added signal photon Quality 3
+               cutarray[ 0] = "6010002002092970028250440000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
+               cutarray[ 1] = "6120002002092970028250440000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
+               cutarray[ 2] = "5010002002092970028250440000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
+               cutarray[ 3] = "5240002002092970028250440000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
+               cutarray[ 4] = "5250002002092970028250440000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%                   
+   } else if ( trainConfig == 118){ // cleaner cuts, min R = 35 cm
+               cutarray[ 0] = "6010001007092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
+               cutarray[ 1] = "6120001007092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
+               cutarray[ 2] = "5010001007092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
+               cutarray[ 3] = "5240001007092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
+               cutarray[ 4] = "5250001007092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
+   } else if ( trainConfig == 119){ // cleaner cuts, min R = 35 cm added signal
+               cutarray[ 0] = "6010002007092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
+               cutarray[ 1] = "6120002007092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
+               cutarray[ 2] = "5010002007092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
+               cutarray[ 3] = "5240002007092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
+               cutarray[ 4] = "5250002007092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%                   
+   } else if ( trainConfig == 120){ // cleaner cuts, photon Quality 1, min R = 35 cm
+               cutarray[ 0] = "6010001007092970028250420000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
+               cutarray[ 1] = "6120001007092970028250420000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
+               cutarray[ 2] = "5010001007092970028250420000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
+               cutarray[ 3] = "5240001007092970028250420000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
+               cutarray[ 4] = "5250001007092970028250420000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
+   } else if ( trainConfig == 121){ // cleaner cuts added signal, photon Quality 1, min R = 35 cm
+               cutarray[ 0] = "6010002007092970028250420000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
+               cutarray[ 1] = "6120002007092970028250420000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
+               cutarray[ 2] = "5010002007092970028250420000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
+               cutarray[ 3] = "5240002007092970028250420000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
+               cutarray[ 4] = "5250002007092970028250420000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%                   
+   } else if ( trainConfig == 122){ // cleaner cuts, photon Quality 3, min R = 35 cm
+               cutarray[ 0] = "6010001007092970028250440000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
+               cutarray[ 1] = "6120001007092970028250440000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
+               cutarray[ 2] = "5010001007092970028250440000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
+               cutarray[ 3] = "5240001007092970028250440000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
+               cutarray[ 4] = "5250001007092970028250440000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
+   } else if ( trainConfig == 123){ // cleaner cuts added signal, photon Quality 3, min R = 35 cm
+               cutarray[ 0] = "6010002007092970028250440000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
+               cutarray[ 1] = "6120002007092970028250440000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
+               cutarray[ 2] = "5010002007092970028250440000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
+               cutarray[ 3] = "5240002007092970028250440000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
+               cutarray[ 4] = "5250002007092970028250440000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%                   
        } else {
                Error(Form("GammaConvV1_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
                return;
@@ -807,24 +879,24 @@ void AddTask_GammaConvV1_PbPb(  Int_t trainConfig = 1,  //change different set o
 
        for(Int_t i = 0; i<numberOfCuts; i++){
                analysisCuts[i] = new AliConversionCuts();
-               if (trainConfig == 1 ||trainConfig == 5 || trainConfig == 9 || trainConfig == 13 || trainConfig == 17 || trainConfig == 21 || trainConfig == 25 || trainConfig == 29 ){ //|| trainConfig == 33 || trainConfig == 37 || trainConfig == 41 
+               if (trainConfig == 1 ||trainConfig == 5 || trainConfig == 9 || trainConfig == 13 || trainConfig == 17 || trainConfig == 21 || trainConfig == 25 || trainConfig == 29 || trainConfig == 33 || trainConfig == 37){ // || trainConfig == 41 
                        if (i == 0 && doWeighting)  analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_0005TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0005V0M");
                        if (i == 1 && doWeighting)  analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_0510TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0510V0M");
                        if (i == 2 && doWeighting)  analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_0010TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0010V0M");
                        if (i == 3 && doWeighting)  analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE,fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_1020TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_1020V0M");
                        if (i == 4 && doWeighting)  analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE,fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_0020TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0020V0M");
-               } else if (trainConfig == 2 ||trainConfig == 6 || trainConfig == 10 || trainConfig == 14 || trainConfig == 18 || trainConfig == 22 || trainConfig == 26 || trainConfig == 30  ){ //|| trainConfig == 34 || trainConfig == 38 || trainConfig == 42
+               } else if (trainConfig == 2 ||trainConfig == 6 || trainConfig == 10 || trainConfig == 14 || trainConfig == 18 || trainConfig == 22 || trainConfig == 26 || trainConfig == 30  || trainConfig == 34 || trainConfig == 38){ // || trainConfig == 42
                        if (i == 0 && doWeighting)  analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_2040TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_2040V0M");
                        if (i == 1 && doWeighting) analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_4060TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_4060V0M");
                        if (i == 2 && doWeighting) analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_6080TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_6080V0M");
                        if (i == 3 && doWeighting)  analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_4080TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_4080V0M");
-               } else if (trainConfig == 3 ||trainConfig == 7 || trainConfig == 11 || trainConfig == 15 || trainConfig == 19 || trainConfig == 23 || trainConfig == 27 || trainConfig == 31 ){ //|| trainConfig == 35 || trainConfig == 39 || trainConfig == 43 
+               } else if (trainConfig == 3 ||trainConfig == 7 || trainConfig == 11 || trainConfig == 15 || trainConfig == 19 || trainConfig == 23 || trainConfig == 27 || trainConfig == 31 || trainConfig == 35 || trainConfig == 39 ){ //|| trainConfig == 43 
                        if (i == 0 && doWeighting)  analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_0005TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0005V0M");
                        if (i == 1 && doWeighting)  analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_0510TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0510V0M");
                        if (i == 2 && doWeighting)  analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_0010TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0010V0M");
                        if (i == 3 && doWeighting)  analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE,fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_1020TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_1020V0M");
                        if (i == 0 && doWeighting)  analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_0020TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0020V0M");
-               } else if (trainConfig == 4 ||trainConfig == 8 || trainConfig == 12 || trainConfig == 16 || trainConfig == 20 || trainConfig == 24 || trainConfig == 28 || trainConfig == 32 ){ //|| trainConfig == 36 || trainConfig == 40 || trainConfig == 44 
+               } else if (trainConfig == 4 ||trainConfig == 8 || trainConfig == 12 || trainConfig == 16 || trainConfig == 20 || trainConfig == 24 || trainConfig == 28 || trainConfig == 32 || trainConfig == 36 || trainConfig == 40){ // || trainConfig == 44 
                        if (i == 0 && doWeighting)  analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_2040TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_2040V0M");
                        if (i == 1 && doWeighting)  analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_4060TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_4060V0M");
                        if (i == 2  && doWeighting) analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_6080TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_6080V0M");
@@ -887,7 +959,7 @@ void AddTask_GammaConvV1_PbPb(  Int_t trainConfig = 1,  //change different set o
                        }       
                }         
                
-               if ( trainConfig == 70 || trainConfig == 72  || trainConfig == 74  || trainConfig == 76  || trainConfig == 78  || trainConfig == 80  || trainConfig == 82  || trainConfig == 84 || trainConfig == 86  || trainConfig == 88  || trainConfig == 90 || trainConfig == 92 || trainConfig == 94 || trainConfig == 96  || trainConfig == 98  || trainConfig == 100 || trainConfig == 102  || trainConfig == 104 || trainConfig == 106 || trainConfig == 108 || trainConfig == 110){
+               if ( trainConfig == 70 || trainConfig == 72  || trainConfig == 74  || trainConfig == 76  || trainConfig == 78  || trainConfig == 80  || trainConfig == 82  || trainConfig == 84 || trainConfig == 86  || trainConfig == 88  || trainConfig == 90 || trainConfig == 92 || trainConfig == 94 || trainConfig == 96  || trainConfig == 98  || trainConfig == 100 || trainConfig == 102  || trainConfig == 104 || trainConfig == 106 || trainConfig == 108 || trainConfig == 110 || trainConfig == 112 || trainConfig == 114 || trainConfig == 116 || trainConfig == 118 || trainConfig == 120 || trainConfig == 122){
                        if (periodName.CompareTo("LHC14a1a") ==0 || periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
                                if ( i == 0 && doWeighting)  analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_0005TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_0005TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0005V0M","Eta_Fit_Data_PbPb_2760GeV_0005V0M");
                                if ( i == 1 && doWeighting)  analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_0510TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_0510TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0510V0M","Eta_Fit_Data_PbPb_2760GeV_0510V0M");
@@ -896,7 +968,7 @@ void AddTask_GammaConvV1_PbPb(  Int_t trainConfig = 1,  //change different set o
                                if ( i == 4 && doWeighting)  analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_2050TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_2050TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_2050V0M","Eta_Fit_Data_PbPb_2760GeV_2050V0M");
                        }       
                } 
-               if ( trainConfig == 71 || trainConfig == 73  || trainConfig == 75  || trainConfig == 77  || trainConfig == 79  || trainConfig == 81  || trainConfig == 83  || trainConfig == 85 || trainConfig == 87  || trainConfig == 89  || trainConfig == 91 || trainConfig == 93 || trainConfig == 95 || trainConfig == 97  || trainConfig == 99  || trainConfig == 101 || trainConfig == 103  || trainConfig == 105 || trainConfig == 107 || trainConfig == 109 || trainConfig == 111){
+               if ( trainConfig == 71 || trainConfig == 73  || trainConfig == 75  || trainConfig == 77  || trainConfig == 79  || trainConfig == 81  || trainConfig == 83  || trainConfig == 85 || trainConfig == 87  || trainConfig == 89  || trainConfig == 91 || trainConfig == 93 || trainConfig == 95 || trainConfig == 97  || trainConfig == 99  || trainConfig == 101 || trainConfig == 103  || trainConfig == 105 || trainConfig == 107 || trainConfig == 109 || trainConfig == 111 || trainConfig == 113 || trainConfig == 115 || trainConfig == 117 || trainConfig == 119 || trainConfig == 121 || trainConfig == 123){
                        if (periodName.CompareTo("LHC14a1a") ==0 || periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
                                if ( i == 0 && doWeighting)  analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_0005TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_0005TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0005V0M","Eta_Fit_Data_PbPb_2760GeV_0005V0M");
                                if ( i == 1 && doWeighting)  analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_0510TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_0510TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0510V0M","Eta_Fit_Data_PbPb_2760GeV_0510V0M");
index f6a5450..6887e3e 100644 (file)
@@ -962,15 +962,15 @@ 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 == 173) {
-               cutarray[ 0] = "8000011002092170008260400000"; mesonCutArray[ 0] = "01621035009000"; //clean cuts ///New STANDARD CUT
-               cutarray[ 1] = "8000011002092170008260420000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //no MCP smearing
-               cutarray[ 2] = "8000011002092170008260430000"; mesonCutArray[ 2] = "01621035009000";  //old standard eta=0.9 y=0.8 
-               cutarray[ 3] = "8000011002092170008260440000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
+               cutarray[ 0] = "8000011002092170008260400000"; mesonCutArray[ 0] = "01621035009000";  // all Photon Qualities
+               cutarray[ 1] = "8000011002092170008260420000"; mesonCutArray[ 1] = "01621035009000";  // only photon quality 1
+               cutarray[ 2] = "8000011002092170008260430000"; mesonCutArray[ 2] = "01621035009000";  // only photon quality 2
+               cutarray[ 3] = "8000011002092170008260440000"; mesonCutArray[ 3] = "01621035009000";  // only photon quality 3
        } else if (trainConfig == 174) {
-               cutarray[ 0] = "8000012002092170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // BG track multiplicity
-               cutarray[ 1] = "8000012002092170008260420000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //no MCP smearing      
-               cutarray[ 2] = "8000012002092170008260430000"; mesonCutArray[ 2] = "01621035009000";  //old standard eta=0.9 y=0.8 
-               cutarray[ 3] = "8000012002092170008260440000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;           
+               cutarray[ 0] = "8000012002092170008260400000"; mesonCutArray[ 0] = "01621035009000";  // all Photon Qualities
+               cutarray[ 1] = "8000012002092170008260420000"; mesonCutArray[ 1] = "01621035009000";  // only photon quality 1
+               cutarray[ 2] = "8000012002092170008260430000"; mesonCutArray[ 2] = "01621035009000";  // only photon quality 2
+               cutarray[ 3] = "8000012002092170008260440000"; mesonCutArray[ 3] = "01621035009000";  // only photon quality 3
        } else if (trainConfig == 175) {
                cutarray[ 0] = "8000011002092170008260410000"; mesonCutArray[ 0] = "01621035009000";  //no shared electron cut on
                cutarray[ 1] = "8000011002192170008260400000"; mesonCutArray[ 1] = "01621035009000";  //single pt cut 0.1
@@ -981,6 +981,16 @@ void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of
                cutarray[ 1] = "8000012002192170008260400000"; mesonCutArray[ 1] = "01621035009000";  //single pt cut 0.1
                cutarray[ 2] = "8000012002292170008260400000"; mesonCutArray[ 2] = "01621035009000";  //single pt cut 0.15
                cutarray[ 3] = "8000012002492170008260400000"; mesonCutArray[ 3] = "01621035009000";  //single pt cut 0.075
+       } else if (trainConfig == 177) {
+               cutarray[ 0] = "8000011007092170008260400000"; mesonCutArray[ 0] = "01621035009000";  // min R =35 cm
+               cutarray[ 1] = "8000011007092170008260420000"; mesonCutArray[ 1] = "01621035009000";  // min R =35 cm & only photon quality 1
+               cutarray[ 2] = "8000011007092170008260430000"; mesonCutArray[ 2] = "01621035009000";  // min R =35 cm & only photon quality 2
+               cutarray[ 3] = "8000011007092170008260440000"; mesonCutArray[ 3] = "01621035009000";  // min R =35 cm & only photon quality 3
+       } else if (trainConfig == 178) {
+               cutarray[ 0] = "8000012007092170008260400000"; mesonCutArray[ 0] = "01621035009000";  // min R =35 cm
+               cutarray[ 1] = "8000012007092170008260420000"; mesonCutArray[ 1] = "01621035009000";  // min R =35 cm & only photon quality 1
+               cutarray[ 2] = "8000012007092170008260430000"; mesonCutArray[ 2] = "01621035009000";  // min R =35 cm & only photon quality 2
+               cutarray[ 3] = "8000012007092170008260440000"; mesonCutArray[ 3] = "01621035009000";  // min R =35 cm & only photon quality 3
        } else {
                Error(Form("GammaConvV1_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
                return;
@@ -1018,7 +1028,7 @@ void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of
        for(Int_t i = 0; i<numberOfCuts; i++){
                
                analysisCuts[i] = new AliConversionCuts();
-               if ( trainConfig == 1 || trainConfig == 3 || trainConfig == 5 || trainConfig == 7 || trainConfig == 9 || trainConfig == 11 || trainConfig == 13 || trainConfig == 15|| trainConfig == 17|| trainConfig == 19 || trainConfig == 21 || trainConfig == 133 || trainConfig == 135 || trainConfig == 137 || trainConfig == 139 || trainConfig == 141 || trainConfig == 143 || trainConfig == 145 || trainConfig == 147 || trainConfig == 149 || trainConfig == 151 || trainConfig == 173 || trainConfig == 175   ){
+               if ( trainConfig == 1 || trainConfig == 3 || trainConfig == 5 || trainConfig == 7 || trainConfig == 9 || trainConfig == 11 || trainConfig == 13 || trainConfig == 15|| trainConfig == 17|| trainConfig == 19 || trainConfig == 21 || trainConfig == 133 || trainConfig == 135 || trainConfig == 137 || trainConfig == 139 || trainConfig == 141 || trainConfig == 143 || trainConfig == 145 || trainConfig == 147 || trainConfig == 149 || trainConfig == 151 || trainConfig == 173 || trainConfig == 175 || trainConfig == 177  ){
                        if (doWeighting){
                                if (generatorName.CompareTo("DPMJET")==0){
                                        analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_DPMJET_LHC13b2_efix_pPb_5023GeV_MBV0A", "Eta_DPMJET_LHC13b2_efix_pPb_5023GeV_MBV0A", "","Pi0_Fit_Data_pPb_5023GeV_MBV0A","Eta_Fit_Data_pPb_5023GeV_MBV0A");
@@ -1027,7 +1037,7 @@ void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of
                                }
                        }
                }   
-               if ( trainConfig == 2 || trainConfig == 4 || trainConfig == 6 || trainConfig == 8 || trainConfig == 10 || trainConfig == 12 || trainConfig == 14 || trainConfig == 16|| trainConfig == 18|| trainConfig == 20|| trainConfig == 22 || trainConfig == 134 || trainConfig == 136 || trainConfig == 138 || trainConfig == 140 || trainConfig == 142 || trainConfig == 144 || trainConfig == 146 || trainConfig == 148 || trainConfig == 150 || trainConfig == 152 || trainConfig == 174 || trainConfig == 176){
+               if ( trainConfig == 2 || trainConfig == 4 || trainConfig == 6 || trainConfig == 8 || trainConfig == 10 || trainConfig == 12 || trainConfig == 14 || trainConfig == 16|| trainConfig == 18|| trainConfig == 20|| trainConfig == 22 || trainConfig == 134 || trainConfig == 136 || trainConfig == 138 || trainConfig == 140 || trainConfig == 142 || trainConfig == 144 || trainConfig == 146 || trainConfig == 148 || trainConfig == 150 || trainConfig == 152 || trainConfig == 174 || trainConfig == 176 || trainConfig == 178){
                        if (doWeighting){
                                analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13e7_addSig_pPb_5023GeV_MBV0A", "Eta_Hijing_LHC13e7_addSig_pPb_5023GeV_MBV0A", "","Pi0_Fit_Data_pPb_5023GeV_MBV0A","Eta_Fit_Data_pPb_5023GeV_MBV0A");
                        }
index 70463f4..03ec3e9 100644 (file)
@@ -104,15 +104,45 @@ void AddTask_GammaConvV1_pp(  Int_t trainConfig = 1,  //change different set of
       cutarray[ 0] = "0000011002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , only Minbias MC
       cutarray[ 1] = "0001011002093663003800000000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD, V0AND
       cutarray[ 2] = "0000011002093260003800000000"; mesonCutArray[2] = "01631031009000"; //standard cut Gamma pp 2-76TeV
-      cutarray[ 3] = "0000011002093260003800000000"; mesonCutArray[3] = "01631031009000"; //standard cut Gamma pp 2-76TeV , only boxes
+      cutarray[ 3] = "0000011002093260003800000000"; mesonCutArray[3] = "01631031009000"; //standard cut Gamma pp 2-76TeV
    } else if (trainConfig == 3) {
       cutarray[ 0] = "0002011002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , only Minbias MC
       cutarray[ 1] = "0003011002093663003800000000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD, V0AND , only Minbias MC
       cutarray[ 2] = "0002012002093663003800000000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , only Boxes MC
       cutarray[ 3] = "0003012002093663003800000000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD, V0AND, only Boxes MC
-   } else {
-      Error(Form("GammaConvV1_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
-      return;
+   } else if (trainConfig == 4) {
+      cutarray[ 0] = "0000011002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , all photon qualities
+      cutarray[ 1] = "0000011002093663003800020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 1
+      cutarray[ 2] = "0000011002093663003800030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 2
+         cutarray[ 3] = "0000011002093663003800040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 3
+   } else if (trainConfig == 5) {
+      cutarray[ 0] = "0000011007093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , all photon qualities, min R = 35 cm
+      cutarray[ 1] = "0000011007093663003800020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 1, min R = 35 cm
+      cutarray[ 2] = "0000011007093663003800030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 2, min R = 35 cm
+         cutarray[ 3] = "0000011007093663003800040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 3, min R = 35 cm
+   } else if (trainConfig == 4) {
+      cutarray[ 0] = "0000011002083663003200000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, all photon qualities
+      cutarray[ 1] = "0000011002083663003200020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 1
+      cutarray[ 2] = "0000011002083663003200030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 2
+         cutarray[ 3] = "0000011002083663003200040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 3   
+   } else if (trainConfig == 5) {
+      cutarray[ 0] = "0000011007083663003200000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, all photon qualities, min R = 35 cm
+      cutarray[ 1] = "0000011007083663003200020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 1, min R = 35 cm
+      cutarray[ 2] = "0000011007083663003200030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 2, min R = 35 cm
+         cutarray[ 3] = "0000011007083663003200040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 3, min R = 35 cm
+   } else if (trainConfig == 6) {
+      cutarray[ 0] = "0000011002083663000200000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 7TeV, all photon qualities
+      cutarray[ 1] = "0000011002083663000200020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 1
+      cutarray[ 2] = "0000011002083663000200030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 2
+         cutarray[ 3] = "0000011002083663000200040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 3   
+   } else if (trainConfig == 7) {
+      cutarray[ 0] = "0000011007083663000200000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 7TeV, all photon qualities, min R = 35 cm
+      cutarray[ 1] = "0000011007083663000200020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 1, min R = 35 cm
+      cutarray[ 2] = "0000011007083663000200030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 2, min R = 35 cm
+         cutarray[ 3] = "0000011007083663000200040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 3, min R = 35 cm          
+       } else {
+               Error(Form("GammaConvV1_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
+               return;
    }
 
    TList *ConvCutList = new TList();