PWGGA/GammaConv/macros/AddTask_GammaConvV1_PbPb.C

   1 void AddTask_GammaConvV1_PbPb(  Int_t trainConfig = 1,  //change different set of cuts
   2                               Bool_t isMC   = kFALSE, //run MC
   3                               Int_t enableQAMesonTask = 0, //enable QA in AliAnalysisTaskGammaConvV1
   4                               Int_t enableQAPhotonTask = 0, // enable additional QA task
   5                               TString fileNameInputForWeighting = "MCSpectraInput.root", // path to file for weigting input
   6                               Int_t headerSelectionInt = 0,  // 1 pi0 header, 2 eta header, 3 both (only for "named" boxes)
   7                               TString cutnumberAODBranch = "1000000060084000001500000",
   8                               TString periodName = "LHC13d2",  //name of the period for added signals and weighting
   9                                                           Bool_t doWeighting = kFALSE  //enable Weighting
  10                            ) {
  11
  12         // ================= Load Librariers =================================
  13         gSystem->Load("libCore.so");
  14         gSystem->Load("libTree.so");
  15         gSystem->Load("libGeom.so");
  16         gSystem->Load("libVMC.so");
  17         gSystem->Load("libPhysics.so");
  18         gSystem->Load("libMinuit");
  19         gSystem->Load("libSTEERBase");
  20         gSystem->Load("libESD");
  21         gSystem->Load("libAOD");
  22         gSystem->Load("libANALYSIS");
  23         gSystem->Load("libANALYSISalice");
  24         gSystem->Load("libPWGGAGammaConv.so");
  25         gSystem->Load("libCDB.so");
  26         gSystem->Load("libSTEER.so");
  27         gSystem->Load("libSTEERBase.so");
  28         gSystem->Load("libTENDER.so");
  29         gSystem->Load("libTENDERSupplies.so");
  30
  31         Int_t isHeavyIon = 1;
  32
  33         // ================== GetAnalysisManager ===============================
  34         AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
  35         if (!mgr) {
  36                 Error(Form("AddTask_GammaConvV1_%i",trainConfig), "No analysis manager found.");
  37                 return ;
  38         }
  39
  40         // ================== GetInputEventHandler =============================
  41         AliVEventHandler *inputHandler=mgr->GetInputEventHandler();
  42
  43         //========= Add PID Reponse to ANALYSIS manager ====
  44         if(!(AliPIDResponse*)mgr->GetTask("PIDResponseTask")){
  45                 gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
  46                 AddTaskPIDResponse(isMC);
  47         }
  48
  49         //=========  Set Cutnumber for V0Reader ================================
  50         TString cutnumberPhoton = "000084001001500000000";
  51         TString cutnumberEvent = "1000000";
  52         AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
  53         //========= Add V0 Reader to  ANALYSIS manager if not yet existent =====
  54         if( !(AliV0ReaderV1*)mgr->GetTask("V0ReaderV1") ){
  55                 AliV0ReaderV1 *fV0ReaderV1 = new AliV0ReaderV1("V0ReaderV1");
  56
  57                 fV0ReaderV1->SetUseOwnXYZCalculation(kTRUE);
  58                 fV0ReaderV1->SetCreateAODs(kFALSE);// AOD Output
  59                 fV0ReaderV1->SetUseAODConversionPhoton(kTRUE);
  60
  61                 if (!mgr) {
  62                         Error("AddTask_V0ReaderV1", "No analysis manager found.");
  63                         return;
  64                 }
  65
  66                 AliConvEventCuts *fEventCuts=NULL;
  67                 if(cutnumberEvent!=""){
  68                         fEventCuts= new AliConvEventCuts(cutnumberEvent.Data(),cutnumberEvent.Data());
  69                         fEventCuts->SetPreSelectionCutFlag(kTRUE);
  70                         if(fEventCuts->InitializeCutsFromCutString(cutnumberEvent.Data())){
  71                                 fV0ReaderV1->SetEventCuts(fEventCuts);
  72                                 fEventCuts->SetFillCutHistograms("",kTRUE);
  73                         }
  74                 }
  75
  76
  77                 // Set AnalysisCut Number
  78                 AliConversionPhotonCuts *fCuts=NULL;
  79                 if(cutnumberPhoton!=""){
  80                         fCuts= new AliConversionPhotonCuts(cutnumberPhoton.Data(),cutnumberPhoton.Data());
  81                         fCuts->SetPreSelectionCutFlag(kTRUE);
  82                         fCuts->SetIsHeavyIon(isHeavyIon);
  83                         if(fCuts->InitializeCutsFromCutString(cutnumberPhoton.Data())){
  84                                 fV0ReaderV1->SetConversionCuts(fCuts);
  85                                 fCuts->SetFillCutHistograms("",kTRUE);
  86                         }
  87                 }
  88
  89                 if(inputHandler->IsA()==AliAODInputHandler::Class()){
  90                 // AOD mode
  91                         fV0ReaderV1->SetDeltaAODBranchName(Form("GammaConv_%s_gamma",cutnumberAODBranch.Data()));
  92                 }
  93                 fV0ReaderV1->Init();
  94
  95                 AliLog::SetGlobalLogLevel(AliLog::kInfo);
  96
  97                 //connect input V0Reader
  98                 mgr->AddTask(fV0ReaderV1);
  99                 mgr->ConnectInput(fV0ReaderV1,0,cinput);
 100
 101         }
 102
 103         //================================================
 104         //========= Add task to the ANALYSIS manager =====
 105         //================================================
 106         AliAnalysisTaskGammaConvV1 *task=NULL;
 107         task= new AliAnalysisTaskGammaConvV1(Form("GammaConvV1_%i",trainConfig));
 108         task->SetIsHeavyIon(isHeavyIon);
 109         task->SetIsMC(isMC);
 110         // Cut Numbers to use in Analysis
 111         Int_t numberOfCuts = 5;
 112
 113         TString *eventCutArray = new TString[numberOfCuts];
 114         TString *photonCutArray = new TString[numberOfCuts];
 115         TString *mesonCutArray = new TString[numberOfCuts];
 116
 117         if (trainConfig == 1){ // Standard cuts
 118                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "042092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
 119                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "042092970023220000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
 120                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "042092970023220000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
 121                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "042092970023220000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
 122                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "042092970023220000000"; mesonCutArray[ 4] = "01522045009000"; // 0-20%
 123         } else if (trainConfig == 2) { // Standard cuts
 124                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "042092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
 125                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "042092970023220000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
 126                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "042092970023220000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
 127                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "042092970023220000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
 128                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "042092970023220000000"; mesonCutArray[ 4] = "01522065009000"; // 40-90%
 129         } else if (trainConfig == 3) { // Standard cuts only added signals
 130                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "042092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
 131                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "042092970023220000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
 132                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "042092970023220000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
 133                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "042092970023220000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
 134                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "042092970023220000000"; mesonCutArray[ 4] = "01522045009000"; // 0-20%
 135         } else if (trainConfig == 4) { // Standard cuts only added signals
 136                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "042092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
 137                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "042092970023220000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
 138                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "042092970023220000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
 139                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "042092970023220000000"; mesonCutArray[ 3] = "01522065009000"; // 20-40%
 140                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "042092970023220000000"; mesonCutArray[ 4] = "01522065009000"; // 40-90%
 141         } else if (trainConfig == 5){ // R-minCut 7.5 cm
 142                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "049092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
 143                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "049092970023220000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
 144                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "049092970023220000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
 145                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "049092970023220000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
 146                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "049092970023220000000"; mesonCutArray[ 4] = "01522045009000"; // 0-20%
 147         } else if (trainConfig == 6) { // R-minCut 7.5 cm
 148                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "049092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
 149                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "049092970023220000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
 150                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "049092970023220000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
 151                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "049092970023220000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
 152                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "049092970023220000000"; mesonCutArray[ 4] = "01522065009000"; // 40-90%
 153         } else if (trainConfig == 7) {// R-minCut 7.5 cm
 154                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "049092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
 155                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "049092970023220000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
 156                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "049092970023220000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
 157                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "049092970023220000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
 158                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "049092970023220000000"; mesonCutArray[ 4] = "01522045009000"; // 0-20%
 159         } else if (trainConfig == 8) { // R-minCut 7.5 cm
 160                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "040092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
 161                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "049092970023220000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
 162                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "049092970023220000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
 163                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "049092970023220000000"; mesonCutArray[ 3] = "01522065009000"; // 20-40%
 164                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "049092970023220000000"; mesonCutArray[ 4] = "01522065009000"; // 40-90%
 165         } else if (trainConfig == 9){ // R-minCut 12.5 cm
 166                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "048092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
 167                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "048092970023220000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
 168                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "048092970023220000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
 169                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "048092970023220000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
 170                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "048092970023220000000"; mesonCutArray[ 4] = "01522045009000"; // 0-20%
 171         } else if (trainConfig == 10) { // R-minCut 12.5 cm
 172                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "048092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
 173                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "048092970023220000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
 174                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "048092970023220000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
 175                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "048092970023220000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
 176                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "048092970023220000000"; mesonCutArray[ 4] = "01522065009000"; // 40-90%
 177         } else if (trainConfig == 11) {// R-minCut 12.5 cm
 178                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "048092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
 179                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "048092970023220000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
 180                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "048092970023220000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
 181                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "048092970023220000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
 182                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "048092970023220000000"; mesonCutArray[ 4] = "01522045009000"; // 0-20%
 183         } else if (trainConfig == 12) { // R-minCut 12.5 cm
 184                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "040092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
 185                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "048092970023220000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
 186                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "048092970023220000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
 187                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "048092970023220000000"; mesonCutArray[ 3] = "01522065009000"; // 20-40%
 188                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "048092970023220000000"; mesonCutArray[ 4] = "01522065009000"; // 40-90%
 189         } else  if (trainConfig == 13){ // LHC10h standard, eta 0.65, y = 0.6
 190                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "032092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
 191                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "032092970023220000000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
 192                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "032092970023220000000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
 193                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "032092970023220000000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
 194                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "032092970023220000000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
 195         } else if (trainConfig == 14) {  // LHC10h standard, eta 0.65, y = 0.6
 196                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "032092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
 197                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "032092970023220000000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
 198                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "032092970023220000000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
 199                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "032092970023220000000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
 200                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "032092970023220000000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
 201         } else if (trainConfig == 15) { // LHC10h standard, eta 0.65, y = 0.6  added signals
 202                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "032092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
 203                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "032092970023220000000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
 204                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "032092970023220000000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
 205                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "032092970023220000000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
 206                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "032092970023220000000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
 207         } else if (trainConfig == 16) { // LHC10h standard, eta 0.65, y = 0.6  added signals
 208                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "032092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
 209                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "032092970023220000000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
 210                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "032092970023220000000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
 211                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "032092970023220000000"; mesonCutArray[ 3] = "01523065009000"; // 20-40%
 212                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "032092970023220000000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
 213         } else  if (trainConfig == 17){ // LHC10h standard, eta 0.65, y = 0.6, photon quality 1
 214                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "032092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
 215                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "032092970023220020000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
 216                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "032092970023220020000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
 217                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "032092970023220020000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
 218                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "032092970023220020000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
 219         } else if (trainConfig == 18) {  // LHC10h standard, eta 0.65, y = 0.6, photon quality 1
 220                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "032092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
 221                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "032092970023220020000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
 222                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "032092970023220020000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
 223                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "032092970023220020000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
 224                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "032092970023220020000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
 225         } else if (trainConfig == 19) { // LHC10h standard, eta 0.65, y = 0.6, photon quality 1  added signal
 226                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "032092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
 227                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "032092970023220020000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
 228                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "032092970023220020000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
 229                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "032092970023220020000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
 230                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "032092970023220020000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
 231         } else if (trainConfig == 20) { // LHC10h standard, eta 0.65, y = 0.6, photon quality 1 added signal
 232                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "032092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
 233                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "032092970023220020000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
 234                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "032092970023220020000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
 235                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "032092970023220020000"; mesonCutArray[ 3] = "01523065009000"; // 20-40%
 236                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "032092970023220020000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
 237         } else  if (trainConfig == 21){ // LHC10h standard, eta 0.65, y = 0.6, photon quality 2
 238                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "032092970023220030000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
 239                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "032092970023220030000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
 240                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "032092970023220030000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
 241                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "032092970023220030000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
 242                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "032092970023220030000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
 243         } else if (trainConfig == 22) {  // LHC10h standard, eta 0.65, y = 0.6, photon quality 2
 244                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "032092970023220030000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
 245                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "032092970023220030000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
 246                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "032092970023220030000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
 247                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "032092970023220030000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
 248                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "032092970023220030000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
 249         } else if (trainConfig == 23) { // LHC10h standard, eta 0.65, y = 0.6, photon quality 2  added signal
 250                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "032092970023220030000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
 251                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "032092970023220030000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
 252                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "032092970023220030000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
 253                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "032092970023220030000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
 254                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "032092970023220030000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
 255         } else if (trainConfig == 24) { // LHC10h standard, eta 0.65, y = 0.6, photon quality 2 added signal
 256                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "032092970023220030000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
 257                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "032092970023220030000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
 258                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "032092970023220030000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
 259                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "032092970023220030000"; mesonCutArray[ 3] = "01523065009000"; // 20-40%
 260                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "032092970023220030000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
 261         } else  if (trainConfig == 25){ // LHC10h standard, eta 0.65, y = 0.6, photon quality 3
 262                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "032092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
 263                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "032092970023220040000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
 264                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "032092970023220040000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
 265                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "032092970023220040000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
 266                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "032092970023220040000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
 267         } else if (trainConfig == 26) {  // LHC10h standard, eta 0.65, y = 0.6, photon quality 3
 268                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "032092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
 269                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "032092970023220040000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
 270                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "032092970023220040000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
 271                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "032092970023220040000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
 272                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "032092970023220040000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
 273         } else if (trainConfig == 27) { // LHC10h standard, eta 0.65, y = 0.6, photon quality 3  added signal
 274                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "032092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
 275                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "032092970023220040000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
 276                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "032092970023220040000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
 277                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "032092970023220040000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
 278                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "032092970023220040000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
 279         } else if (trainConfig == 28) { // LHC10h standard, eta 0.65, y = 0.6, photon quality 3 added signal
 280                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "032092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
 281                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "032092970023220040000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
 282                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "032092970023220040000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
 283                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "032092970023220040000"; mesonCutArray[ 3] = "01523065009000"; // 20-40%
 284                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "032092970023220040000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
 285         } else  if (trainConfig == 29){ // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm
 286                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "037092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
 287                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "037092970023220000000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
 288                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "037092970023220000000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
 289                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "037092970023220000000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
 290                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "037092970023220000000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
 291         } else if (trainConfig == 30) {  // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm
 292                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "037092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
 293                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "037092970023220000000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
 294                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "037092970023220000000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
 295                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "037092970023220000000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
 296                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "037092970023220000000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
 297         } else if (trainConfig == 31) { // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm  added signals
 298                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "037092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
 299                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "037092970023220000000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
 300                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "037092970023220000000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
 301                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "037092970023220000000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
 302                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "037092970023220000000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
 303         } else if (trainConfig == 32) { // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm  added signals
 304                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "037092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
 305                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "037092970023220000000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
 306                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "037092970023220000000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
 307                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "037092970023220000000"; mesonCutArray[ 3] = "01523065009000"; // 20-40%
 308                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "037092970023220000000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
 309         } else  if (trainConfig == 33){ // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 1
 310                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "037092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
 311                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "037092970023220020000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
 312                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "037092970023220020000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
 313                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "037092970023220020000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
 314                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "037092970023220020000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
 315         } else if (trainConfig == 34) {  // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 1
 316                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "037092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
 317                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "037092970023220020000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
 318                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "037092970023220020000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
 319                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "037092970023220020000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
 320                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "037092970023220020000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
 321         } else if (trainConfig == 35) { // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 1  added signals
 322                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "037092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
 323                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "037092970023220020000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
 324                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "037092970023220020000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
 325                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "037092970023220020000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
 326                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "037092970023220020000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
 327         } else if (trainConfig == 36) { // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 1  added signals
 328                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "037092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
 329                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "037092970023220020000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
 330                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "037092970023220020000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
 331                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "037092970023220020000"; mesonCutArray[ 3] = "01523065009000"; // 20-40%
 332                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "037092970023220020000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
 333         } else  if (trainConfig == 37){ // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 3
 334                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "037092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
 335                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "037092970023220040000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
 336                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "037092970023220040000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
 337                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "037092970023220040000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
 338                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "037092970023220040000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
 339         } else if (trainConfig == 38) {  // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 3
 340                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "037092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
 341                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "037092970023220040000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
 342                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "037092970023220040000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
 343                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "037092970023220040000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
 344                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "037092970023220040000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
 345         } else if (trainConfig == 39) { // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 3  added signals
 346                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "037092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
 347                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "037092970023220040000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
 348                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "037092970023220040000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
 349                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "037092970023220040000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
 350                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "037092970023220040000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
 351         } else if (trainConfig == 40) { // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 3  added signals
 352                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "037092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
 353                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "037092970023220040000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
 354                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "037092970023220040000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
 355                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "037092970023220040000"; mesonCutArray[ 3] = "01523065009000"; // 20-40%
 356                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "037092970023220040000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
 357         } else if (trainConfig == 41){ // Standard cuts, eta 0.9, only to be run on data :kSemiCentral
 358                 eventCutArray[ 0] = "6016001"; photonCutArray[ 0] = "002092970023220000000"; mesonCutArray[ 0] = "01525045009000"; // 0-5%
 359                 eventCutArray[ 1] = "6126001"; photonCutArray[ 1] = "002092970023220000000"; mesonCutArray[ 1] = "01525045009000"; // 5-10%
 360                 eventCutArray[ 2] = "5016001"; photonCutArray[ 2] = "002092970023220000000"; mesonCutArray[ 2] = "01525045009000"; // 0-10%
 361                 eventCutArray[ 3] = "6126001"; photonCutArray[ 3] = "002092970023220000000"; mesonCutArray[ 3] = "01525045009000"; // 10-20%
 362                 eventCutArray[ 4] = "5026001"; photonCutArray[ 4] = "002092970023220000000"; mesonCutArray[ 4] = "01525045009000"; // 0-20%
 363         } else if (trainConfig == 42) { // Standard cuts, eta 0.9, only to be run on data :kSemiCentral
 364                 eventCutArray[ 0] = "5236001"; photonCutArray[ 0] = "002092970023220000000"; mesonCutArray[ 0] = "01525065009000";
 365                 eventCutArray[ 1] = "5346001"; photonCutArray[ 1] = "002092970023220000000"; mesonCutArray[ 1] = "01525065009000";
 366                 eventCutArray[ 2] = "5456001"; photonCutArray[ 2] = "002092970023220000000"; mesonCutArray[ 2] = "01525065009000";
 367                 eventCutArray[ 3] = "5566001"; photonCutArray[ 3] = "002092970023220000000"; mesonCutArray[ 3] = "01525065009000";
 368                 eventCutArray[ 4] = "5676001"; photonCutArray[ 4] = "002092970023220000000"; mesonCutArray[ 4] = "01525065009000";
 369         } else if (trainConfig == 43){ // Standard cuts, eta 0.9, only to be run on data
 370                 eventCutArray[ 0] = "5230001"; photonCutArray[ 0] = "002092970023220000000"; mesonCutArray[ 0] = "01525045009000";
 371                 eventCutArray[ 1] = "5340001"; photonCutArray[ 1] = "002092970023220000000"; mesonCutArray[ 1] = "01525065009000";
 372                 eventCutArray[ 2] = "5450001"; photonCutArray[ 2] = "002092970023220000000"; mesonCutArray[ 2] = "01525065009000";
 373                 eventCutArray[ 3] = "5560001"; photonCutArray[ 3] = "002092970023220000000"; mesonCutArray[ 3] = "01525065009000";
 374                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970023220000000"; mesonCutArray[ 4] = "01525065009000";
 375         } else if ( trainConfig == 44){ // qt elipse cut 0.05
 376                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028220000000"; mesonCutArray[ 0] = "01525065009000"; // 0-5%
 377                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028220000000"; mesonCutArray[ 1] = "01525065009000"; // 5-10%
 378                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028220000000"; mesonCutArray[ 2] = "01525065009000"; // 0-10%
 379                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002092970028220000000"; mesonCutArray[ 3] = "01525065009000"; // 10-20%
 380                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "002092970028220000000"; mesonCutArray[ 4] = "01525065009000"; // 0-20%
 381         } else if ( trainConfig == 45) { // qt elipse cut 0.05
 382                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "002092970028220000000"; mesonCutArray[ 0] = "01525065009000"; // 20-40%
 383                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "002092970028220000000"; mesonCutArray[ 1] = "01525065009000"; // 40-60%
 384                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "002092970028220000000"; mesonCutArray[ 2] = "01525065009000"; // 60-80%
 385                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002092970028220000000"; mesonCutArray[ 3] = "01525065009000"; // 40-80%
 386                 eventCutArray[ 4] = "5350001"; photonCutArray[ 4] = "002092970028220000000"; mesonCutArray[ 4] = "01525065009000"; // 30-50%
 387         } else if ( trainConfig == 46){ // qt elipse cut 0.05
 388                 eventCutArray[ 0] = "5230001"; photonCutArray[ 0] = "002092970028220000000"; mesonCutArray[ 0] = "01525065009000"; // 20-30%
 389                 eventCutArray[ 1] = "5340001"; photonCutArray[ 1] = "002092970028220000000"; mesonCutArray[ 1] = "01525065009000"; // 30-40%
 390                 eventCutArray[ 2] = "5450001"; photonCutArray[ 2] = "002092970028220000000"; mesonCutArray[ 2] = "01525065009000"; // 40-50%
 391                 eventCutArray[ 3] = "5560001"; photonCutArray[ 3] = "002092970028220000000"; mesonCutArray[ 3] = "01525065009000"; // 50-60%
 392                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028220000000"; mesonCutArray[ 4] = "01525065009000"; // 60-70%
 393         } else if ( trainConfig == 47){ // cos(theta_point) cut 0.85
 394                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970023220400000"; mesonCutArray[ 0] = "01525065009000"; // 0-5%
 395                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970023220400000"; mesonCutArray[ 1] = "01525065009000"; // 5-10%
 396                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970023220400000"; mesonCutArray[ 2] = "01525065009000"; // 0-10%
 397                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002092970023220400000"; mesonCutArray[ 3] = "01525065009000"; // 10-20%
 398                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "002092970023220400000"; mesonCutArray[ 4] = "01525065009000"; // 0-20%
 399         } else if ( trainConfig == 48) { // cos(theta_point) cut 0.85
 400                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "002092970023220400000"; mesonCutArray[ 0] = "01525065009000"; // 20-40%
 401                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "002092970023220400000"; mesonCutArray[ 1] = "01525065009000"; // 40-60%
 402                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "002092970023220400000"; mesonCutArray[ 2] = "01525065009000"; // 60-80%
 403                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002092970023220400000"; mesonCutArray[ 3] = "01525065009000"; // 40-80%
 404                 eventCutArray[ 4] = "5350001"; photonCutArray[ 4] = "002092970023220400000"; mesonCutArray[ 4] = "01525065009000"; // 30-50%
 405         } else if ( trainConfig == 49){ // cos(theta_point) cut 0.85
 406                 eventCutArray[ 0] = "5230001"; photonCutArray[ 0] = "002092970023220400000"; mesonCutArray[ 0] = "01525065009000"; // 20-30%
 407                 eventCutArray[ 1] = "5340001"; photonCutArray[ 1] = "002092970023220400000"; mesonCutArray[ 1] = "01525065009000"; // 30-40%
 408                 eventCutArray[ 2] = "5450001"; photonCutArray[ 2] = "002092970023220400000"; mesonCutArray[ 2] = "01525065009000"; // 40-50%
 409                 eventCutArray[ 3] = "5560001"; photonCutArray[ 3] = "002092970023220400000"; mesonCutArray[ 3] = "01525065009000"; // 50-60%
 410                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970023220400000"; mesonCutArray[ 4] = "01525065009000"; // 60-70%
 411         } else if ( trainConfig == 50){ // psi pair 2D 0.05
 412                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970023260000000"; mesonCutArray[ 0] = "01525065009000"; // 0-5%
 413                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970023260000000"; mesonCutArray[ 1] = "01525065009000"; // 5-10%
 414                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970023260000000"; mesonCutArray[ 2] = "01525065009000"; // 0-10%
 415                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002092970023260000000"; mesonCutArray[ 3] = "01525065009000"; // 10-20%
 416                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "002092970023260000000"; mesonCutArray[ 4] = "01525065009000"; // 0-20%
 417         } else if ( trainConfig == 51) { // psi pair 2D 0.05
 418                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "002092970023260000000"; mesonCutArray[ 0] = "01525065009000"; // 20-40%
 419                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "002092970023260000000"; mesonCutArray[ 1] = "01525065009000"; // 40-60%
 420                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "002092970023260000000"; mesonCutArray[ 2] = "01525065009000"; // 60-80%
 421                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002092970023260000000"; mesonCutArray[ 3] = "01525065009000"; // 40-80%
 422                 eventCutArray[ 4] = "5350001"; photonCutArray[ 4] = "002092970023260000000"; mesonCutArray[ 4] = "01525065009000"; // 30-50%
 423         } else if ( trainConfig == 52){ // psi pair 2D 0.05
 424                 eventCutArray[ 0] = "5230001"; photonCutArray[ 0] = "002092970023260000000"; mesonCutArray[ 0] = "01525065009000"; // 20-30%
 425                 eventCutArray[ 1] = "5340001"; photonCutArray[ 1] = "002092970023260000000"; mesonCutArray[ 1] = "01525065009000"; // 30-40%
 426                 eventCutArray[ 2] = "5450001"; photonCutArray[ 2] = "002092970023260000000"; mesonCutArray[ 2] = "01525065009000"; // 40-50%
 427                 eventCutArray[ 3] = "5560001"; photonCutArray[ 3] = "002092970023260000000"; mesonCutArray[ 3] = "01525065009000"; // 50-60%
 428                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970023260000000"; mesonCutArray[ 4] = "01525065009000"; // 60-70%
 429         } else if ( trainConfig == 53){ // psi pair 2D 0.1
 430                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970023250000000"; mesonCutArray[ 0] = "01525065009000"; // 0-5%
 431                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970023250000000"; mesonCutArray[ 1] = "01525065009000"; // 5-10%
 432                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970023250000000"; mesonCutArray[ 2] = "01525065009000"; // 0-10%
 433                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002092970023250000000"; mesonCutArray[ 3] = "01525065009000"; // 10-20%
 434                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "002092970023250000000"; mesonCutArray[ 4] = "01525065009000"; // 0-20%
 435         } else if ( trainConfig == 54) { // psi pair 2D 0.1
 436                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "002092970023250000000"; mesonCutArray[ 0] = "01525065009000"; // 20-40%
 437                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "002092970023250000000"; mesonCutArray[ 1] = "01525065009000"; // 40-60%
 438                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "002092970023250000000"; mesonCutArray[ 2] = "01525065009000"; // 60-80%
 439                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002092970023250000000"; mesonCutArray[ 3] = "01525065009000"; // 40-80%
 440                 eventCutArray[ 4] = "5350001"; photonCutArray[ 4] = "002092970023250000000"; mesonCutArray[ 4] = "01525065009000"; // 30-50%
 441         } else if ( trainConfig == 55){ // psi pair 2D 0.1
 442                 eventCutArray[ 0] = "5230001"; photonCutArray[ 0] = "002092970023250000000"; mesonCutArray[ 0] = "01525065009000"; // 20-30%
 443                 eventCutArray[ 1] = "5340001"; photonCutArray[ 1] = "002092970023250000000"; mesonCutArray[ 1] = "01525065009000"; // 30-40%
 444                 eventCutArray[ 2] = "5450001"; photonCutArray[ 2] = "002092970023250000000"; mesonCutArray[ 2] = "01525065009000"; // 40-50%
 445                 eventCutArray[ 3] = "5560001"; photonCutArray[ 3] = "002092970023250000000"; mesonCutArray[ 3] = "01525065009000"; // 50-60%
 446                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970023250000000"; mesonCutArray[ 4] = "01525065009000"; // 60-70%
 447         } else if ( trainConfig == 56){ // cleaner cuts
 448                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 449                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 450                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 451                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 452                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 453         } else if ( trainConfig == 57){ // cleaner cuts added signal
 454                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 455                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 456                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 457                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 458                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 459         } else if ( trainConfig == 58) { // cleaner cuts
 460                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 20-40%
 461                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 40-60%
 462                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 60-80%
 463                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 40-80%
 464                 eventCutArray[ 4] = "5350001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 30-50%
 465         } else if ( trainConfig == 59) { // cleaner cuts added signal
 466                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 20-40%
 467                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 40-60%
 468                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 60-80%
 469                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 40-80%
 470                 eventCutArray[ 4] = "5350002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 30-50%
 471         } else if ( trainConfig == 60){ // cleaner cuts
 472                 eventCutArray[ 0] = "5230001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 20-30%
 473                 eventCutArray[ 1] = "5340001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 30-40%
 474                 eventCutArray[ 2] = "5450001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 40-50%
 475                 eventCutArray[ 3] = "5560001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 50-60%
 476                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 60-70%
 477         } else if ( trainConfig == 61){ // cleaner cuts added signal
 478                 eventCutArray[ 0] = "5230002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 20-30%
 479                 eventCutArray[ 1] = "5340002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 30-40%
 480                 eventCutArray[ 2] = "5450002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 40-50%
 481                 eventCutArray[ 3] = "5560002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 50-60%
 482                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 60-70%
 483         } else if ( trainConfig == 62){ // cleaner cuts
 484                 eventCutArray[ 0] = "6230001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 485                 eventCutArray[ 1] = "6340001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 486                 eventCutArray[ 2] = "6450001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 487                 eventCutArray[ 3] = "6560001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 488                 eventCutArray[ 4] = "6670001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 489         } else if ( trainConfig == 63){ // cleaner cuts added signal
 490                 eventCutArray[ 0] = "6230002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 491                 eventCutArray[ 1] = "6340002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 492                 eventCutArray[ 2] = "6450002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 493                 eventCutArray[ 3] = "6560002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 494                 eventCutArray[ 4] = "6670002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 495         } else if ( trainConfig == 64){ // cleaner cuts
 496                 eventCutArray[ 0] = "6780001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 497                 eventCutArray[ 1] = "6890001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 498                 eventCutArray[ 2] = "5670001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 499                 eventCutArray[ 3] = "5780001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 500                 eventCutArray[ 4] = "5890001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 501         } else if ( trainConfig == 65){ // cleaner cuts added signal
 502                 eventCutArray[ 0] = "6780002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 503                 eventCutArray[ 1] = "6890002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 504                 eventCutArray[ 2] = "5670002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 505                 eventCutArray[ 3] = "5780002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 506                 eventCutArray[ 4] = "5890002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 507         } else if ( trainConfig == 66){ // cleaner cuts
 508                 eventCutArray[ 0] = "7010001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 509                 eventCutArray[ 1] = "7120001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 510                 eventCutArray[ 2] = "7230001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 511                 eventCutArray[ 3] = "7340001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 512                 eventCutArray[ 4] = "7450001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 513         } else if ( trainConfig == 67){ // cleaner cuts added signal
 514                 eventCutArray[ 0] = "7010002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 515                 eventCutArray[ 1] = "7120002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 516                 eventCutArray[ 2] = "7230002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 517                 eventCutArray[ 3] = "7340002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 518                 eventCutArray[ 4] = "7450002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 519         } else if ( trainConfig == 68){ // cleaner cuts
 520                 eventCutArray[ 0] = "7560001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 521                 eventCutArray[ 1] = "7670001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 522                 eventCutArray[ 2] = "7780001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 523                 eventCutArray[ 3] = "7890001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 10-20%
 524                 eventCutArray[ 4] = "7090001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 0-20%
 525         } else if ( trainConfig == 69){ // cleaner cuts added signal
 526                 eventCutArray[ 0] = "7560002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 527                 eventCutArray[ 1] = "7670002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 528                 eventCutArray[ 2] = "7780002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 529                 eventCutArray[ 3] = "7890002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 10-20%
 530                 eventCutArray[ 4] = "7090002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 0-20%
 531         } else if ( trainConfig == 70){ // variation eta  0.65
 532                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "032092970028250400000"; mesonCutArray[ 0]= "01523065000000"; // 0-5%
 533                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "032092970028250400000"; mesonCutArray[ 1]= "01523065000000"; // 5-10%
 534                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "032092970028250400000"; mesonCutArray[ 2]= "01523065000000"; // 0-10%
 535                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "032092970028250400000"; mesonCutArray[ 3]= "01523065000000"; // 20-40%
 536                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "032092970028250400000"; mesonCutArray[ 4]= "01523065000000"; // 20-50%
 537         } else if ( trainConfig == 71){ // variation eta  0.65 added signal
 538                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "032092970028250400000"; mesonCutArray[ 0]= "01523065000000"; // 0-5%
 539                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "032092970028250400000"; mesonCutArray[ 1]= "01523065000000"; // 5-10%
 540                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "032092970028250400000"; mesonCutArray[ 2]= "01523065000000"; // 0-10%
 541                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "032092970028250400000"; mesonCutArray[ 3]= "01523065000000"; // 20-40%
 542                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "032092970028250400000"; mesonCutArray[ 4]= "01523065000000"; // 20-50%
 543         } else if ( trainConfig == 72){ // variation eta  0.75
 544                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "042092970028250400000"; mesonCutArray[ 0]= "01522065009000"; // 0-5%
 545                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "042092970028250400000"; mesonCutArray[ 1]= "01522065009000"; // 5-10%
 546                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "042092970028250400000"; mesonCutArray[ 2]= "01522065009000"; // 0-10%
 547                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "042092970028250400000"; mesonCutArray[ 3]= "01522065009000"; // 20-40%
 548                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "042092970028250400000"; mesonCutArray[ 4]= "01522065009000"; // 20-50%
 549         } else if ( trainConfig == 73){ // variation eta  0.75 added signal
 550                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "042092970028250400000"; mesonCutArray[ 0]= "01522065009000"; // 0-5%
 551                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "042092970028250400000"; mesonCutArray[ 1]= "01522065009000"; // 5-10%
 552                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "042092970028250400000"; mesonCutArray[ 2]= "01522065009000"; // 0-10%
 553                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "042092970028250400000"; mesonCutArray[ 3]= "01522065009000"; // 20-40%
 554                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "042092970028250400000"; mesonCutArray[ 4]= "01522065009000"; // 20-50%
 555         } else if ( trainConfig == 74){ // single pt 0.075
 556                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002492970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 557                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002492970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 558                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002492970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 559                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002492970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 560                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002492970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 561         } else if ( trainConfig == 75){ // single pt 0.075 added signal
 562                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002492970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 563                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002492970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 564                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002492970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 565                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002492970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 566                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002492970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 567         } else if ( trainConfig == 76){ // single pt 0.1
 568                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002192970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 569                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002192970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 570                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002192970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 571                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002192970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 572                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002192970028250400000"; mesonCutArray[ 4]= "01525065000000"; //20-50%
 573         } else if ( trainConfig == 77){ // single pt 0.1 added signal
 574                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002192970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 575                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002192970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 576                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002192970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 577                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002192970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 578                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002192970028250400000"; mesonCutArray[ 4]= "01525065000000"; //20-50%
 579         } else if ( trainConfig == 78){ // variation TPC cls 0.7
 580                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002062970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 581                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002062970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 582                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002062970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 583                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002062970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 584                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002062970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 585         } else if ( trainConfig == 79){ // variation TPC cls 0.7 added signal
 586                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002062970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 587                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002062970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 588                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002062970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 589                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002062970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 590                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002062970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 591         } else if ( trainConfig == 80){ // variation TPC cls 0.35
 592                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002082970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 593                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002082970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 594                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002082970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 595                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002082970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 596                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002082970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 597         } else if ( trainConfig == 81){ // variation TPC cls 0.35 added signal
 598                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002082970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 599                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002082970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 600                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002082970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 601                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002082970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 602                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002082970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 603         } else if ( trainConfig == 82){ // variation edEdx  -4,5
 604                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002093970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 605                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002093970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 606                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002093970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 607                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002093970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 608                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002093970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 609         } else if ( trainConfig == 83){ // variation edEdx  -4,5  added signal
 610                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002093970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 611                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002093970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 612                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002093970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 613                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002093970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 614                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002093970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 615         } else if ( trainConfig == 84){ // variation edEdx  -2.5,4
 616                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002096970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 617                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002096970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 618                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002096970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 619                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002096970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 620                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002096970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 621         } else if ( trainConfig == 85){ // variation edEdx  -2.5,4  added signal
 622                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002096970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 623                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002096970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 624                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002096970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 625                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002096970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 626                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002096970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 627         } else if ( trainConfig == 86){ //variation pion p dEdx 0.3-5.
 628                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092951028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 629                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092951028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 630                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092951028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 631                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092951028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 632                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092951028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 633         } else if ( trainConfig == 87){ //variation pion p dEdx 0.3-5.  added signal
 634                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092951028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 635                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092951028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 636                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092951028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 637                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092951028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 638                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092951028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 639         } else if ( trainConfig == 88){ // TOF el. PID -3,5
 640                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970038250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 641                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970038250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 642                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970038250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 643                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970038250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 644                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970038250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 645         } else if ( trainConfig == 89){ // TOF el. PID -3,5  added signal
 646                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970038250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 647                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970038250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 648                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970038250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 649                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970038250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 650                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970038250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 651         } else if ( trainConfig == 90){ // TOF el. PID -2,3
 652                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970048250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 653                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970048250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 654                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970048250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 655                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970048250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 656                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970048250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 657         } else if ( trainConfig == 91){ // TOF el. PID -2,3  added signal
 658                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970048250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 659                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970048250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 660                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970048250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 661                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970048250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 662                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970048250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 663         } else if ( trainConfig == 92){ // qt 0.03
 664                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970029250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 665                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970029250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 666                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970029250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 667                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970029250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 668                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970029250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 669         } else if ( trainConfig == 93){ // qt 0.03  added signal
 670                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970029250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 671                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970029250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 672                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970029250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 673                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970029250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 674                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970029250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 675         } else if ( trainConfig == 94){ // qt 0.07 no2D
 676                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970022250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 677                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970022250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 678                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970022250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 679                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970022250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 680                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970022250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 681         } else if ( trainConfig         == 95){ // qt 0.07 no2D  added signal
 682                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970022250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 683                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970022250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 684                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970022250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 685                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970022250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 686                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970022250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 687         } else if ( trainConfig == 96){ // chi2  50.
 688                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028150400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 689                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028150400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 690                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028150400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 691                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028150400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 692                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028150400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 693         } else if ( trainConfig == 97){ // chi2  50.  added signal
 694                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028150400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 695                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028150400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 696                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028150400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 697                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028150400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 698                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028150400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 699         } else if ( trainConfig == 98){ // chi2  20.
 700                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028850400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 701                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028850400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 702                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028850400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 703                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028850400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 704                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028850400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 705         } else if ( trainConfig == 99){ // chi2  20.  added signal
 706                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028850400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 707                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028850400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 708                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028850400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 709                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028850400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 710                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028850400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 711         } else if ( trainConfig == 100){ // psi pair 0.05
 712                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028260400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 713                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028260400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 714                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028260400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 715                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028260400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 716                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028260400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 717         } else if ( trainConfig == 101){ // psi pair 0.05  added signal
 718                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028260400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 719                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028260400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 720                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028260400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 721                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028260400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 722                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028260400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 723         } else if ( trainConfig == 102){ // cosPA -1
 724                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250000000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 725                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250000000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 726                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250000000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 727                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250000000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 728                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250000000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 729         } else if ( trainConfig == 103){ // cosPA -1  added signal
 730                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028250000000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 731                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028250000000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 732                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028250000000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 733                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028250000000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 734                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028250000000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 735         } else if ( trainConfig == 104){ // variation alpha 0.75
 736                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0]= "01525055000000"; // 0-5%
 737                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1]= "01525055000000"; // 5-10%
 738                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2]= "01525055000000"; // 0-10%
 739                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3]= "01525055000000"; // 20-40%
 740                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4]= "01525055000000"; // 20-50%
 741         } else if ( trainConfig == 105){ // variation alpha 0.75  added signal
 742                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0]= "01525055000000"; // 0-5%
 743                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1]= "01525055000000"; // 5-10%
 744                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2]= "01525055000000"; // 0-10%
 745                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3]= "01525055000000"; // 20-40%
 746                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4]= "01525055000000"; // 20-50%
 747         } else if ( trainConfig == 106){ // variation alpha 0.85
 748                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0]= "01525075000000"; // 0-5%
 749                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1]= "01525075000000"; // 5-10%
 750                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2]= "01525075000000"; // 0-10%
 751                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3]= "01525075000000"; // 20-40%
 752                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4]= "01525075000000"; // 20-50%
 753         } else if ( trainConfig == 107){ // variation alpha 0.85  added signal
 754                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0]= "01525075000000"; // 0-5%
 755                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1]= "01525075000000"; // 5-10%
 756                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2]= "01525075000000"; // 0-10%
 757                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3]= "01525075000000"; // 20-40%
 758                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4]= "01525075000000"; // 20-50%
 759         } else if ( trainConfig == 108){ // psi pair 0.2
 760                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028280400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 761                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028280400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 762                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028280400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 763                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028280400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 764                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028280400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 765         } else if ( trainConfig == 109){ // psi pair 0.2  added signal
 766                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028280400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
 767                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028280400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
 768                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028280400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
 769                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028280400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
 770                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028280400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
 771         } else if ( trainConfig == 110){ // cleaner cuts
 772                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 773                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 774                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 775                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 776                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 777         } else if ( trainConfig == 111){ // cleaner cuts added signal
 778                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 779                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 780                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 781                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 782                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 783         } else if ( trainConfig == 112){ // cleaner cuts photon Quality 1
 784                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250420000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 785                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250420000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 786                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250420000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 787                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250420000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 788                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250420000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 789         } else if ( trainConfig == 113){ // cleaner cuts added signal photon Quality 1
 790                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028250420000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 791                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028250420000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 792                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028250420000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 793                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028250420000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 794                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028250420000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 795         } else if ( trainConfig == 114){ // cleaner cuts photon Quality 2
 796                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250430000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 797                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250430000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 798                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250430000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 799                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250430000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 800                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250430000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 801         } else if ( trainConfig == 115){ // cleaner cuts added signal photon Quality 2
 802                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028250430000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 803                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028250430000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 804                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028250430000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 805                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028250430000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 806                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028250430000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 807         } else if ( trainConfig == 116){ // cleaner cuts photon Quality 3
 808                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250440000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 809                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250440000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 810                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250440000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 811                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250440000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 812                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250440000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 813         } else if ( trainConfig == 117){ // cleaner cuts added signal photon Quality 3
 814                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028250440000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 815                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028250440000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 816                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028250440000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 817                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028250440000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 818                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028250440000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 819         } else if ( trainConfig == 118){ // cleaner cuts, min R = 35 cm
 820                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "007092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 821                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "007092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 822                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "007092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 823                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "007092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 824                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "007092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 825         } else if ( trainConfig == 119){ // cleaner cuts, min R = 35 cm added signal
 826                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "007092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 827                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "007092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 828                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "007092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 829                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "007092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 830                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "007092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 831         } else if ( trainConfig == 120){ // cleaner cuts, photon Quality 1, min R = 35 cm
 832                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "007092970028250420000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 833                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "007092970028250420000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 834                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "007092970028250420000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 835                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "007092970028250420000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 836                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "007092970028250420000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 837         } else if ( trainConfig == 121){ // cleaner cuts added signal, photon Quality 1, min R = 35 cm
 838                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "007092970028250420000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 839                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "007092970028250420000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 840                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "007092970028250420000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 841                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "007092970028250420000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 842                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "007092970028250420000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 843         } else if ( trainConfig == 122){ // cleaner cuts, photon Quality 3, min R = 35 cm
 844                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "007092970028250440000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 845                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "007092970028250440000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 846                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "007092970028250440000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 847                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "007092970028250440000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 848                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "007092970028250440000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 849         } else if ( trainConfig == 123){ // cleaner cuts added signal, photon Quality 3, min R = 35 cm
 850                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "007092970028250440000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
 851                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "007092970028250440000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
 852                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "007092970028250440000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
 853                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "007092970028250440000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
 854                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "007092970028250440000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
 855         } else if ( trainConfig == 124){ // cleaner cuts, specific centrality trigger selection
 856                 eventCutArray[ 0] = "5080001"; photonCutArray[ 0] = "007092970028250440000"; mesonCutArray[ 0] = "01525065000000"; // normal centrality
 857                 eventCutArray[ 1] = "5083101"; photonCutArray[ 1] = "007092970028250440000"; mesonCutArray[ 1] = "01525065000000"; // kCentral open
 858                 eventCutArray[ 2] = "5083201"; photonCutArray[ 2] = "007092970028250440000"; mesonCutArray[ 2] = "01525065000000"; // kCentral vertex
 859                 eventCutArray[ 3] = "5083301"; photonCutArray[ 3] = "007092970028250440000"; mesonCutArray[ 3] = "01525065000000"; // kCentral both
 860                 eventCutArray[ 4] = "5083401"; photonCutArray[ 4] = "007092970028250440000"; mesonCutArray[ 4] = "01525065000000"; // kSemiCentral open
 861         } else if ( trainConfig == 125){ // cleaner cuts, specific centrality trigger selection
 862                 eventCutArray[ 0] = "5080001"; photonCutArray[ 0] = "007092970028250440000"; mesonCutArray[ 0] = "01525065000000"; // normal centrality
 863                 eventCutArray[ 1] = "5083501"; photonCutArray[ 1] = "007092970028250440000"; mesonCutArray[ 1] = "01525065000000"; // kSemiCentral vertex
 864                 eventCutArray[ 2] = "5083601"; photonCutArray[ 2] = "007092970028250440000"; mesonCutArray[ 2] = "01525065000000"; // kSemiCentral both
 865                 eventCutArray[ 3] = "5083701"; photonCutArray[ 3] = "007092970028250440000"; mesonCutArray[ 3] = "01525065000000"; // kMB 1 only
 866                 eventCutArray[ 4] = "5083801"; photonCutArray[ 4] = "007092970028250440000"; mesonCutArray[ 4] = "01525065000000"; // kMB 2 only
 867         } else {
 868                 Error(Form("GammaConvV1_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
 869                 return;
 870         }
 871
 872         TList *EventCutList = new TList();
 873         TList *ConvCutList = new TList();
 874         TList *MesonCutList = new TList();
 875
 876         TList *HeaderList = new TList();
 877         if (periodName.CompareTo("LHC13d2")==0){
 878                 TObjString *Header1 = new TObjString("pi0_1");
 879                 HeaderList->Add(Header1);
 880         //    TObjString *Header3 = new TObjString("eta_2");
 881         //    HeaderList->Add(Header3);
 882
 883         } else if (periodName.CompareTo("LHC12a17x_fix")==0){
 884                 TObjString *Header1 = new TObjString("PARAM");
 885                 HeaderList->Add(Header1);
 886         } else if (periodName.CompareTo("LHC14a1a")==0){
 887                 if (headerSelectionInt == 1){
 888                         TObjString *Header1 = new TObjString("pi0_1");
 889                         HeaderList->Add(Header1);
 890                 } else if (headerSelectionInt == 2){
 891                         TObjString *Header1 = new TObjString("eta_2");
 892                         HeaderList->Add(Header1);
 893                 } else {
 894                         TObjString *Header1 = new TObjString("pi0_1");
 895                         HeaderList->Add(Header1);
 896                         TObjString *Header2 = new TObjString("eta_2");
 897                         HeaderList->Add(Header2);
 898                 }
 899         } else if (periodName.CompareTo("LHC14a1b")==0 || periodName.CompareTo("LHC14a1c")==0){
 900                 TObjString *Header1 = new TObjString("BOX");
 901                 HeaderList->Add(Header1);
 902         }
 903
 904         EventCutList->SetOwner(kTRUE);
 905         AliConvEventCuts **analysisEventCuts = new AliConvEventCuts*[numberOfCuts];
 906         ConvCutList->SetOwner(kTRUE);
 907         AliConversionPhotonCuts **analysisCuts = new AliConversionPhotonCuts*[numberOfCuts];
 908         MesonCutList->SetOwner(kTRUE);
 909         AliConversionMesonCuts **analysisMesonCuts = new AliConversionMesonCuts*[numberOfCuts];
 910
 911         for(Int_t i = 0; i<numberOfCuts; i++){
 912                 analysisEventCuts[i] = new AliConvEventCuts();
 913                 if (trainConfig == 1 ||trainConfig == 5 || trainConfig == 9 || trainConfig == 13 || trainConfig == 17 || trainConfig == 21 || trainConfig == 25 || trainConfig == 29 || trainConfig == 33 || trainConfig == 37){ // || trainConfig == 41
 914                         if (i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_0005TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0005V0M");
 915                         if (i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_0510TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0510V0M");
 916                         if (i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_0010TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0010V0M");
 917                         if (i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE,fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_1020TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_1020V0M");
 918                         if (i == 4 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE,fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_0020TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0020V0M");
 919                 } else if (trainConfig == 2 ||trainConfig == 6 || trainConfig == 10 || trainConfig == 14 || trainConfig == 18 || trainConfig == 22 || trainConfig == 26 || trainConfig == 30  || trainConfig == 34 || trainConfig == 38){ // || trainConfig == 42
 920                         if (i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_2040TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_2040V0M");
 921                         if (i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_4060TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_4060V0M");
 922                         if (i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_6080TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_6080V0M");
 923                         if (i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_4080TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_4080V0M");
 924                 } else if (trainConfig == 3 ||trainConfig == 7 || trainConfig == 11 || trainConfig == 15 || trainConfig == 19 || trainConfig == 23 || trainConfig == 27 || trainConfig == 31 || trainConfig == 35 || trainConfig == 39 ){ //|| trainConfig == 43
 925                         if (i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_0005TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0005V0M");
 926                         if (i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_0510TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0510V0M");
 927                         if (i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_0010TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0010V0M");
 928                         if (i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE,fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_1020TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_1020V0M");
 929                         if (i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_0020TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0020V0M");
 930                 } else if (trainConfig == 4 ||trainConfig == 8 || trainConfig == 12 || trainConfig == 16 || trainConfig == 20 || trainConfig == 24 || trainConfig == 28 || trainConfig == 32 || trainConfig == 36 || trainConfig == 40){ // || trainConfig == 44
 931                         if (i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_2040TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_2040V0M");
 932                         if (i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_4060TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_4060V0M");
 933                         if (i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_6080TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_6080V0M");
 934                         if (i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_4080TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_4080V0M");
 935                 }
 936
 937                 if (trainConfig == 56 ){
 938                         if (periodName.CompareTo("LHC14a1a") ==0 || periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
 939                                 if ( i == 0 && doWeighting)  analysisEventCuts[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");
 940                                 if ( i == 1 && doWeighting)  analysisEventCuts[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");
 941                                 if ( i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_0010TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_0010TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0010V0M","Eta_Fit_Data_PbPb_2760GeV_0010V0M");
 942                                 if ( i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_1020TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_1020TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_1020V0M","Eta_Fit_Data_PbPb_2760GeV_1020V0M");
 943                                 if ( i == 4 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_0020TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_0020TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0020V0M","Eta_Fit_Data_PbPb_2760GeV_0020V0M");
 944                         }
 945                 }
 946                 if (trainConfig == 57 ){
 947                         if (periodName.CompareTo("LHC14a1a") ==0 || periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
 948                                 if ( i == 0 && doWeighting)  analysisEventCuts[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");
 949                                 if ( i == 1 && doWeighting)  analysisEventCuts[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");
 950                                 if ( i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_0010TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_0010TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0010V0M","Eta_Fit_Data_PbPb_2760GeV_0010V0M");
 951                                 if ( i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_1020TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_1020TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_1020V0M","Eta_Fit_Data_PbPb_2760GeV_1020V0M");
 952                                 if ( i == 4 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_0020TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_0020TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0020V0M","Eta_Fit_Data_PbPb_2760GeV_0020V0M");
 953                         }
 954                 }
 955                 if (trainConfig == 58 ){
 956                         if (periodName.CompareTo("LHC14a1a") ==0 || periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
 957                                 if ( i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_2040TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_2040TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_2040V0M","Eta_Fit_Data_PbPb_2760GeV_2040V0M");
 958                                 if ( i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_4060TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_4060TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_4060V0M","Eta_Fit_Data_PbPb_2760GeV_4060V0M");
 959                                 if ( i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_6080TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_6080TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_6080V0M","Eta_Fit_Data_PbPb_2760GeV_6080V0M");
 960                                 if ( i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_4080TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_4080TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_1020V0M","Eta_Fit_Data_PbPb_2760GeV_1020V0M");
 961                                 if ( i == 4 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_3050TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_3050TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_3050V0M","Eta_Fit_Data_PbPb_2760GeV_3050V0M");
 962                         }
 963                 }
 964                 if (trainConfig == 59 ){
 965                         if (periodName.CompareTo("LHC14a1a") ==0 || periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
 966                                 if ( i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_2040TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_2040TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_2040V0M","Eta_Fit_Data_PbPb_2760GeV_2040V0M");
 967                                 if ( i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_4060TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_4060TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_4060V0M","Eta_Fit_Data_PbPb_2760GeV_4060V0M");
 968                                 if ( i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_6080TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_6080TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_6080V0M","Eta_Fit_Data_PbPb_2760GeV_6080V0M");
 969                                 if ( i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_4080TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_4080TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_1020V0M","Eta_Fit_Data_PbPb_2760GeV_1020V0M");
 970                                 if ( i == 4 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_3050TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_3050TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_3050V0M","Eta_Fit_Data_PbPb_2760GeV_3050V0M");
 971                         }
 972                 }
 973
 974                 if (trainConfig == 60 ){
 975                         if (periodName.CompareTo("LHC14a1a") ==0 || periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
 976                                 if ( i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_2030TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_2030TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_2030V0M","Eta_Fit_Data_PbPb_2760GeV_2030V0M");
 977                                 if ( i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_3040TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_3040TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_3040V0M","Eta_Fit_Data_PbPb_2760GeV_3040V0M");
 978                                 if ( i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_4050TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_4050TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_4050V0M","Eta_Fit_Data_PbPb_2760GeV_4050V0M");
 979                                 if ( i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_5060TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_5060TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_5060V0M","Eta_Fit_Data_PbPb_2760GeV_5060V0M");
 980                                 if ( i == 4 && doWeighting)  analysisEventCuts[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");
 981                         }
 982                 }
 983                 if (trainConfig == 61 ){
 984                         if (periodName.CompareTo("LHC14a1a") ==0 || periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
 985                                 if ( i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_2030TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_2030TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_2030V0M","Eta_Fit_Data_PbPb_2760GeV_2030V0M");
 986                                 if ( i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_3040TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_3040TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_3040V0M","Eta_Fit_Data_PbPb_2760GeV_3040V0M");
 987                                 if ( i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_4050TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_4050TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_4050V0M","Eta_Fit_Data_PbPb_2760GeV_4050V0M");
 988                                 if ( i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_5060TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_5060TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_5060V0M","Eta_Fit_Data_PbPb_2760GeV_5060V0M");
 989                                 if ( i == 4 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_2050TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_2050TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_2050V0M","Eta_Fit_Data_PbPb_2760GeV_2050V0M");
 990                         }
 991                 }
 992
 993                 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){
 994                         if (periodName.CompareTo("LHC14a1a") ==0 || periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
 995                                 if ( i == 0 && doWeighting)  analysisEventCuts[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");
 996                                 if ( i == 1 && doWeighting)  analysisEventCuts[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");
 997                                 if ( i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_0010TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_0010TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0010V0M","Eta_Fit_Data_PbPb_2760GeV_0010V0M");
 998                                 if ( i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_2040TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_2040TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_2040V0M","Eta_Fit_Data_PbPb_2760GeV_2040V0M");
 999                                 if ( i == 4 && doWeighting)  analysisEventCuts[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");
1000                         }
1001                 }
1002                 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){
1003                         if (periodName.CompareTo("LHC14a1a") ==0 || periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
1004                                 if ( i == 0 && doWeighting)  analysisEventCuts[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");
1005                                 if ( i == 1 && doWeighting)  analysisEventCuts[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");
1006                                 if ( i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_0010TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_0010TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0010V0M","Eta_Fit_Data_PbPb_2760GeV_0010V0M");
1007                                 if ( i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_2040TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_2040TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_2040V0M","Eta_Fit_Data_PbPb_2760GeV_2040V0M");
1008                                 if ( i == 4 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_2050TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_2050TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_2050V0M","Eta_Fit_Data_PbPb_2760GeV_2050V0M");
1009                         }
1010                 }
1011
1012                 analysisEventCuts[i]->InitializeCutsFromCutString(eventCutArray[i].Data());
1013                 if (periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
1014                         if (headerSelectionInt == 1) analysisEventCuts[i]->SetAddedSignalPDGCode(111);
1015                         if (headerSelectionInt == 2) analysisEventCuts[i]->SetAddedSignalPDGCode(221);
1016                 }
1017                 EventCutList->Add(analysisEventCuts[i]);
1018                 if (trainConfig == 37 || trainConfig == 38){
1019                         analysisEventCuts[i]->SelectSpecialTrigger(AliVEvent::kMB, "AliVEvent::kMB" );
1020                 }
1021                 if (trainConfig == 39 || trainConfig == 40){
1022                         analysisEventCuts[i]->SelectSpecialTrigger(AliVEvent::kCentral,"AliVEvent::kCentral" );
1023                 }
1024                 if (trainConfig == 41 || trainConfig == 42){
1025                         analysisEventCuts[i]->SelectSpecialTrigger(AliVEvent::kSemiCentral,"AliVEvent::kSemiCentral" );
1026                 }
1027                 analysisEventCuts[i]->SetFillCutHistograms("",kFALSE);
1028
1029                 analysisCuts[i] = new AliConversionPhotonCuts();
1030                 analysisCuts[i]->InitializeCutsFromCutString(photonCutArray[i].Data());
1031                 ConvCutList->Add(analysisCuts[i]);
1032                 analysisCuts[i]->SetFillCutHistograms("",kFALSE);
1033
1034                 analysisMesonCuts[i] = new AliConversionMesonCuts();
1035                 analysisMesonCuts[i]->InitializeCutsFromCutString(mesonCutArray[i].Data());
1036                 MesonCutList->Add(analysisMesonCuts[i]);
1037                 analysisMesonCuts[i]->SetFillCutHistograms("");
1038
1039                 analysisEventCuts[i]->SetAcceptedHeader(HeaderList);
1040         }
1041
1042         task->SetEventCutList(numberOfCuts,EventCutList);
1043         task->SetConversionCutList(numberOfCuts,ConvCutList);
1044         task->SetMesonCutList(numberOfCuts,MesonCutList);
1045         task->SetMoveParticleAccordingToVertex(kTRUE);
1046         task->SetDoMesonAnalysis(kTRUE);
1047         task->SetDoMesonQA(enableQAMesonTask); //Attention new switch for Pi0 QA
1048         task->SetDoPhotonQA(enableQAPhotonTask);  //Attention new switch small for Photon QA
1049
1050         //connect containers
1051         AliAnalysisDataContainer *coutput =
1052                 mgr->CreateContainer(Form("GammaConvV1_%i",trainConfig), TList::Class(),
1053                                                         AliAnalysisManager::kOutputContainer,Form("GammaConvV1_%i.root",trainConfig));
1054
1055         mgr->AddTask(task);
1056         mgr->ConnectInput(task,0,cinput);
1057         mgr->ConnectOutput(task,1,coutput);
1058
1059         return;
1060
1061 }