PWGGA/GammaConv/macros/AddTask_GammaConvV1_pp.C

   1 void AddTask_GammaConvV1_pp(  Int_t trainConfig = 1,                                                                            // change different set of cuts
   2                               Bool_t isMC   = kFALSE,                                                                           // run MC
   3                               Int_t enableQAMesonTask = 0,                                                                      // enable meson QA in AliAnalysisTaskGammaConvV1
   4                               Int_t enableQAPhotonTask = 0,                                                             // enable photon QA in AliAnalysisTaskGammaConvV1
   5                               TString fileNameInputForWeighting = "MCSpectraInput.root",        // path to file for weigting input
   6                               TString cutnumberAODBranch = "0000000060084001001500000",         // cutnumber for AOD branch
   7                                                           TString periodname = "LHC12f1x"                                                               // period name
   8                            ) {
   9
  10         // ================= Load Librariers =================================
  11         gSystem->Load("libCore.so");
  12         gSystem->Load("libTree.so");
  13         gSystem->Load("libGeom.so");
  14         gSystem->Load("libVMC.so");
  15         gSystem->Load("libPhysics.so");
  16         gSystem->Load("libMinuit");
  17         gSystem->Load("libSTEERBase");
  18         gSystem->Load("libESD");
  19         gSystem->Load("libAOD");
  20         gSystem->Load("libANALYSIS");
  21         gSystem->Load("libANALYSISalice");
  22         gSystem->Load("libPWGGAGammaConv.so");
  23         gSystem->Load("libCDB.so");
  24         gSystem->Load("libSTEER.so");
  25         gSystem->Load("libSTEERBase.so");
  26         gSystem->Load("libTENDER.so");
  27         gSystem->Load("libTENDERSupplies.so");
  28
  29         Int_t isHeavyIon = 0;
  30         // ================== GetAnalysisManager ===============================
  31         AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
  32         if (!mgr) {
  33                 Error(Form("AddTask_GammaConvV1_%i",trainConfig), "No analysis manager found.");
  34                 return ;
  35         }
  36
  37         // ================== GetInputEventHandler =============================
  38         AliVEventHandler *inputHandler=mgr->GetInputEventHandler();
  39
  40         //========= Add PID Reponse to ANALYSIS manager ====
  41         if(!(AliPIDResponse*)mgr->GetTask("PIDResponseTask")){
  42                 gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
  43                 AddTaskPIDResponse(isMC);
  44         }
  45
  46         //=========  Set Cutnumber for V0Reader ================================
  47         TString cutnumberPhoton = "002084000002200000000";
  48         TString cutnumberEvent = "0000000";
  49         AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
  50
  51         //========= Add V0 Reader to  ANALYSIS manager if not yet existent =====
  52         if( !(AliV0ReaderV1*)mgr->GetTask("V0ReaderV1") ){
  53                 AliV0ReaderV1 *fV0ReaderV1 = new AliV0ReaderV1("V0ReaderV1");
  54
  55                 fV0ReaderV1->SetUseOwnXYZCalculation(kTRUE);
  56                 fV0ReaderV1->SetCreateAODs(kFALSE);// AOD Output
  57                 fV0ReaderV1->SetUseAODConversionPhoton(kTRUE);
  58
  59                 if (!mgr) {
  60                         Error("AddTask_V0ReaderV1", "No analysis manager found.");
  61                         return;
  62                 }
  63
  64                 AliConvEventCuts *fEventCuts=NULL;
  65                 if(cutnumberEvent!=""){
  66                         fEventCuts= new AliConvEventCuts(cutnumberEvent.Data(),cutnumberEvent.Data());
  67                         fEventCuts->SetPreSelectionCutFlag(kTRUE);
  68                         if(fEventCuts->InitializeCutsFromCutString(cutnumberEvent.Data())){
  69                                 fV0ReaderV1->SetEventCuts(fEventCuts);
  70                                 fEventCuts->SetFillCutHistograms("",kTRUE);
  71                         }
  72                 }
  73
  74                 // Set AnalysisCut Number
  75                 AliConversionPhotonCuts *fCuts=NULL;
  76                 if(cutnumberPhoton!=""){
  77                         fCuts= new AliConversionPhotonCuts(cutnumberPhoton.Data(),cutnumberPhoton.Data());
  78                         fCuts->SetPreSelectionCutFlag(kTRUE);
  79                         fCuts->SetIsHeavyIon(isHeavyIon);
  80                         if(fCuts->InitializeCutsFromCutString(cutnumberPhoton.Data())){
  81                                 fV0ReaderV1->SetConversionCuts(fCuts);
  82                                 fCuts->SetFillCutHistograms("",kTRUE);
  83                         }
  84                 }
  85                 if(inputHandler->IsA()==AliAODInputHandler::Class()){
  86                 // AOD mode
  87                         fV0ReaderV1->SetDeltaAODBranchName(Form("GammaConv_%s_gamma",cutnumberAODBranch.Data()));
  88                 }
  89                 fV0ReaderV1->Init();
  90
  91                 AliLog::SetGlobalLogLevel(AliLog::kInfo);
  92
  93                 //connect input V0Reader
  94                 mgr->AddTask(fV0ReaderV1);
  95                 mgr->ConnectInput(fV0ReaderV1,0,cinput);
  96
  97         }
  98
  99         //================================================
 100         //========= Add task to the ANALYSIS manager =====
 101         //            find input container
 102         AliAnalysisTaskGammaConvV1 *task=NULL;
 103         task= new AliAnalysisTaskGammaConvV1(Form("GammaConvV1_%i",trainConfig));
 104         task->SetIsHeavyIon(isHeavyIon);
 105         task->SetIsMC(isMC);
 106         // Cut Numbers to use in Analysis
 107         Int_t numberOfCuts = 4;
 108
 109         TString *eventCutArray = new TString[numberOfCuts];
 110         TString *photonCutArray = new TString[numberOfCuts];
 111         TString *mesonCutArray = new TString[numberOfCuts];
 112
 113         if(trainConfig == 1){
 114                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , only Minbias MC
 115                 eventCutArray[ 1] = "0001011"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD, V0AND
 116                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002093260003800000000"; mesonCutArray[2] = "01631031009000"; //standard cut Gamma pp 2-76TeV
 117                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002093660003800000000"; mesonCutArray[3] = "01631031009000"; //standard cut Gamma pp 2-76TeV
 118         } else if (trainConfig == 2) {
 119                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , only boxes
 120                 eventCutArray[ 1] = "0001012"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD, V0AND , only boxes
 121                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002093260003800000000"; mesonCutArray[2] = "01631031009000"; //standard cut Gamma pp 2-76TeV , only boxes
 122                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002093660003800000000"; mesonCutArray[3] = "01631031009000"; //standard cut Gamma pp 2-76TeV
 123         } else if (trainConfig == 3) {
 124                 eventCutArray[ 0] = "0002011"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , only Minbias MC
 125                 eventCutArray[ 1] = "0003011"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD, V0AND , only Minbias MC
 126                 eventCutArray[ 2] = "0002012"; photonCutArray[ 2] = "002093663003800000000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , only Boxes MC
 127                 eventCutArray[ 3] = "0003012"; photonCutArray[ 3] = "002093663003800000000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD, V0AND, only Boxes MC
 128         } else if (trainConfig == 4) {
 129                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , all photon qualities
 130                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002093663003800020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 1
 131                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002093663003800030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 2
 132                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002093663003800040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 3
 133         } else if (trainConfig == 5) {
 134                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "007093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , all photon qualities, min R = 35 cm
 135                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "007093663003800020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 1, min R = 35 cm
 136                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "007093663003800030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 2, min R = 35 cm
 137                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "007093663003800040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 3, min R = 35 cm
 138         } else if (trainConfig == 6) {
 139                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002083663003200000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, all photon qualities
 140                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002083663003200020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 1
 141                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002083663003200030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 2
 142                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002083663003200040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 3
 143         } else if (trainConfig == 7) {
 144                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "007083663003200000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, all photon qualities, min R = 35 cm
 145                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "007083663003200020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 1, min R = 35 cm
 146                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "007083663003200030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 2, min R = 35 cm
 147                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "007083663003200040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 3, min R = 35 cm
 148         } else if (trainConfig == 8) {
 149                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002083663000200000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 7TeV, all photon qualities
 150                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002083663000200020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 1
 151                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002083663000200030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 2
 152                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002083663000200040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 3
 153         } else if (trainConfig == 9) {
 154                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "007083663000200000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 7TeV, all photon qualities, min R = 35 cm
 155                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "007083663000200020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 1, min R = 35 cm
 156                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "007083663000200030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 2, min R = 35 cm
 157                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "007083663000200040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 3, min R = 35 cm
 158         } else if (trainConfig == 10) {
 159                 eventCutArray[ 0] = "0002011"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , all photon qualities
 160                 eventCutArray[ 1] = "0002011"; photonCutArray[ 1] = "002093663003800020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 1
 161                 eventCutArray[ 2] = "0002011"; photonCutArray[ 2] = "002093663003800030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 2
 162                 eventCutArray[ 3] = "0002011"; photonCutArray[ 3] = "002093663003800040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 3
 163         } else if (trainConfig == 11) {
 164                 eventCutArray[ 0] = "0002011"; photonCutArray[ 0] = "007093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , all photon qualities, min R = 35 cm
 165                 eventCutArray[ 1] = "0002011"; photonCutArray[ 1] = "007093663003800020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 1, min R = 35 cm
 166                 eventCutArray[ 2] = "0002011"; photonCutArray[ 2] = "007093663003800030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 2, min R = 35 cm
 167                 eventCutArray[ 3] = "0002011"; photonCutArray[ 3] = "007093663003800040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 3, min R = 35 cm
 168         } else if (trainConfig == 12) {
 169                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; //standard cut LHC11h pp 2.76TeV
 170                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "032092970028250400000"; mesonCutArray[1] = "01525065000000"; //variation eta 0.65
 171                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "042092970028250400000"; mesonCutArray[2] = "01525065000000"; //variation eta 0.75
 172                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002092950028250400000"; mesonCutArray[3] = "01525065000000"; //variation pion p dEdx 0.3-5.
 173         } else if (trainConfig == 13) { //added signals
 174                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; //standard cut LHC11h pp 2.76TeV
 175                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "032092970028250400000"; mesonCutArray[1] = "01525065000000"; //variation eta 0.65
 176                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "042092970028250400000"; mesonCutArray[2] = "01525065000000"; //variation eta 0.75
 177                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002092950028250400000"; mesonCutArray[3] = "01525065000000"; //variation pion p dEdx 0.3-5.
 178         } else if (trainConfig == 14) {
 179                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002492970028250400000"; mesonCutArray[0] = "01525065000000"; //variation pt 0.075
 180                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002192970028250400000"; mesonCutArray[1] = "01525065000000"; //variation pt 0.1
 181                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002062970028250400000"; mesonCutArray[2] = "01525065000000"; //variation TPC cls 0.7
 182                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002082970028250400000"; mesonCutArray[3] = "01525065000000"; //variation TPC cls 0.35
 183         } else if (trainConfig == 15) { //added signals
 184                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002492970028250400000"; mesonCutArray[0] = "01525065000000"; //variation pt 0.075
 185                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002192970028250400000"; mesonCutArray[1] = "01525065000000"; //variation pt 0.1
 186                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002062970028250400000"; mesonCutArray[2] = "01525065000000"; //variation TPC cls 0.7
 187                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002082970028250400000"; mesonCutArray[3] = "01525065000000"; //variation TPC cls 0.35
 188         } else if (trainConfig == 16) {
 189                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002093970028250400000"; mesonCutArray[0] = "01525065000000"; //variation edEdx -4,5
 190                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002096970028250400000"; mesonCutArray[1] = "01525065000000"; //variation edEdx -2.5,4
 191                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002092970038250400000"; mesonCutArray[2] = "01525065000000"; //variation TOF el. PID -3,5
 192                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002092970048250400000"; mesonCutArray[3] = "01525065000000"; //variation TOF el. PID -2,3
 193         } else if (trainConfig == 17) { //added signals
 194                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002093970028250400000"; mesonCutArray[0] = "01525065000000"; //variation edEdx -4,5
 195                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002096970028250400000"; mesonCutArray[1] = "01525065000000"; //variation edEdx -2.5,4
 196                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002092970038250400000"; mesonCutArray[2] = "01525065000000"; //variation TOF el. PID -3,5
 197                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002092970048250400000"; mesonCutArray[3] = "01525065000000"; //variation TOF el. PID -2,3
 198         } else if (trainConfig == 18) {
 199                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002092970029250400000"; mesonCutArray[0] = "01525065000000"; //variation qt 0.03
 200                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002092970022250400000"; mesonCutArray[1] = "01525065000000"; //variation qt 0.07 no2D
 201                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002092970028150400000"; mesonCutArray[2] = "01525065000000"; //variation chi2 50.
 202                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002092970028850400000"; mesonCutArray[3] = "01525065000000"; //variation chi2 20.
 203         } else if (trainConfig == 19) { //added signals
 204                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002092970029250400000"; mesonCutArray[0] = "01525065000000"; //variation qt 0.03
 205                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002092970022250400000"; mesonCutArray[1] = "01525065000000"; //variation qt 0.07 no2D
 206                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002092970028150400000"; mesonCutArray[2] = "01525065000000"; //variation chi2 50.
 207                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002092970028850400000"; mesonCutArray[3] = "01525065000000"; //variation chi2 20.
 208         } else if (trainConfig == 20) {
 209                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002092970028260400000"; mesonCutArray[0] = "01525065000000"; //variation psi pair 0.05
 210                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002092970028280400000"; mesonCutArray[1] = "01525065000000"; //variation psi pair 0.2
 211                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002092970028250000000"; mesonCutArray[2] = "01525065000000"; //variation cosPA -1
 212                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525055000000"; //variation alpha 0.75
 213         } else if (trainConfig == 21) { //added signals
 214                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002092970028260400000"; mesonCutArray[0] = "01525065000000"; //variation psi pair 0.05
 215                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002092970028280400000"; mesonCutArray[1] = "01525065000000"; //variation psi pair 0.2
 216                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002092970028250000000"; mesonCutArray[2] = "01525065000000"; //variation cosPA -1
 217                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525055000000"; //variation alpha 0.75
 218         } else if (trainConfig == 22) {
 219                 eventCutArray[ 0] = "0004011"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; // trigger kTRD
 220                 eventCutArray[ 1] = "0005011"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[1] = "01525065000000"; // trigger kEMC
 221                 eventCutArray[ 2] = "0006011"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01525065000000"; // trigger kPHI
 222                 eventCutArray[ 3] = "0007011"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525065000000"; // trigger kHighMult
 223         } else if (trainConfig == 23) {
 224                 eventCutArray[ 0] = "0008011"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; // trigger kEMCEGA
 225                 eventCutArray[ 1] = "0009011"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[1] = "01525065000000"; // trigger kEMCEJE
 226                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01525065000000"; // minimum bias
 227                 eventCutArray[ 3] = "0001111"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525065000000"; // trigger kINT8
 228         } else if (trainConfig == 24) {
 229                 eventCutArray[ 0] = "0004211"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; // trigger kTRD CINT8 HEE
 230                 eventCutArray[ 1] = "0004411"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[1] = "01525065000000"; // trigger kTRD CINT8 HSE
 231                 eventCutArray[ 2] = "0004611"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01525065000000"; // trigger kTRD CINT8 HJE
 232                 eventCutArray[ 3] = "0004811"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525065000000"; // trigger kTRD CINT8 HQU
 233         } else if (trainConfig == 25) {
 234                 eventCutArray[ 0] = "0004111"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; // trigger kTRD CINT7 HEE
 235                 eventCutArray[ 1] = "0004311"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[1] = "01525065000000"; // trigger kTRD CINT7 HSE
 236                 eventCutArray[ 2] = "0004511"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01525065000000"; // trigger kTRD CINT7 HJE
 237                 eventCutArray[ 3] = "0004711"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525065000000"; // trigger kTRD CINT7 HQU
 238         } else if (trainConfig == 26) {
 239                 eventCutArray[ 0] = "0005211"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; // trigger kEMC7
 240                 eventCutArray[ 1] = "0005311"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[1] = "01525065000000"; // trigger kEMC8
 241                 eventCutArray[ 2] = "0006211"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01525065000000"; // trigger kPHI7
 242                 eventCutArray[ 3] = "0006311"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525065000000"; // trigger kPHI8
 243         } else if (trainConfig == 27) {
 244                 eventCutArray[ 0] = "0005111"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; // trigger kEMC1
 245                 eventCutArray[ 1] = "0007111"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[1] = "01525065000000"; // trigger kSHM1
 246                 eventCutArray[ 2] = "0007211"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01525065000000"; // trigger kSHM7
 247                 eventCutArray[ 3] = "0007311"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525065000000"; // trigger kSHM8
 248         } else if (trainConfig == 28) {
 249                 eventCutArray[ 0] = "0008111"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; // trigger kEMCEGA + CINT7
 250                 eventCutArray[ 1] = "0008211"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[1] = "01525065000000"; // trigger kEMCEGA + CINT8
 251                 eventCutArray[ 2] = "0008311"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01525065000000"; // trigger kEMCEG1 + CINT7
 252                 eventCutArray[ 3] = "0008411"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525065000000"; // trigger kEMCEG1 + CINT8
 253         } else if (trainConfig == 29) {
 254                 eventCutArray[ 0] = "0008511"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; // trigger kEMCEG2 + CINT7
 255                 eventCutArray[ 1] = "0008611"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[1] = "01525065000000"; // trigger kEMCEG2 + CINT8
 256                 eventCutArray[ 2] = "0009111"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01525065000000"; // trigger kEMCEJE + CINT7
 257                 eventCutArray[ 3] = "0009211"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525065000000"; // trigger kEMCEJE + CINT8
 258         } else if (trainConfig == 30) {
 259                 eventCutArray[ 0] = "0009311"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; // trigger kEMCEJ1 + CINT7
 260                 eventCutArray[ 1] = "0009411"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[1] = "01525065000000"; // trigger kEMCEJ1 + CINT8
 261                 eventCutArray[ 2] = "0009511"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01525065000000"; // trigger kEMCEJ2 + CINT7
 262                 eventCutArray[ 3] = "0009611"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525065000000"; // trigger kEMCEJ2 + CINT8
 263         } else if (trainConfig == 31) {
 264                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; //new standard cut for pp 8 TeV
 265                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002093570028250400000"; mesonCutArray[1] = "01521065000000"; //variation edEdx -4,5
 266                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002096570028250400000"; mesonCutArray[2] = "01521065000000"; //variation edEdx -2.5,4
 267                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002092550028250400000"; mesonCutArray[3] = "01521065000000"; //variation pion p dEdx 0.3-5.
 268         } else if (trainConfig == 32) { //added signals
 269                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; //new standard cut for pp 8 TeV
 270                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002093570028250400000"; mesonCutArray[1] = "01521065000000"; //variation edEdx -4,5
 271                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002096570028250400000"; mesonCutArray[2] = "01521065000000"; //variation edEdx -2.5,4
 272                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002092550028250400000"; mesonCutArray[3] = "01521065000000"; //variation pion p dEdx 0.3-5.
 273         } else if (trainConfig == 33) {
 274                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002492570028250400000"; mesonCutArray[0] = "01521065000000"; //variation pt 0.075
 275                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002192570028250400000"; mesonCutArray[1] = "01521065000000"; //variation pt 0.1
 276                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002062570028250400000"; mesonCutArray[2] = "01521065000000"; //variation TPC cls 0.7
 277                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002082570028250400000"; mesonCutArray[3] = "01521065000000"; //variation TPC cls 0.35
 278         } else if (trainConfig == 34) { //added signals
 279                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002492570028250400000"; mesonCutArray[0] = "01521065000000"; //variation pt 0.075
 280                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002192570028250400000"; mesonCutArray[1] = "01521065000000"; //variation pt 0.1
 281                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002062570028250400000"; mesonCutArray[2] = "01521065000000"; //variation TPC cls 0.7
 282                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002082570028250400000"; mesonCutArray[3] = "01521065000000"; //variation TPC cls 0.35
 283         } else if (trainConfig == 35) {
 284                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002092270028250400000"; mesonCutArray[0] = "01521065000000"; //variation pidEdx 1,-10
 285                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002092370028250400000"; mesonCutArray[1] = "01521065000000"; //variation pidEdx 2.5,-10
 286                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01521065000000"; //variation pidEdx 3,-10
 287                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002092500028250400000"; mesonCutArray[3] = "01521065000000"; //variation pion p dEdx 0.5-5
 288         } else if (trainConfig == 36) { //added signals
 289                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002092270028250400000"; mesonCutArray[0] = "01521065000000"; //variation pidEdx 1,-10
 290                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002092370028250400000"; mesonCutArray[1] = "01521065000000"; //variation pidEdx 2.5,-10
 291                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01521065000000"; //variation pidEdx 3,-10
 292                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002092500028250400000"; mesonCutArray[3] = "01521065000000"; //variation pion p dEdx 0.5-5
 293         } else if (trainConfig == 37) {
 294                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002092570029250400000"; mesonCutArray[0] = "01521065000000"; //variation qt 0.03
 295                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002092570022250400000"; mesonCutArray[1] = "01521065000000"; //variation qt 0.07 no2D
 296                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002092570028150400000"; mesonCutArray[2] = "01521065000000"; //variation chi2 50.
 297                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002092570028850400000"; mesonCutArray[3] = "01521065000000"; //variation chi2 20.
 298         } else if (trainConfig == 38) { //added signals
 299                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002092570029250400000"; mesonCutArray[0] = "01521065000000"; //variation qt 0.03
 300                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002092570022250400000"; mesonCutArray[1] = "01521065000000"; //variation qt 0.07 no2D
 301                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002092570028150400000"; mesonCutArray[2] = "01521065000000"; //variation chi2 50.
 302                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002092570028850400000"; mesonCutArray[3] = "01521065000000"; //variation chi2 20.
 303         } else if (trainConfig == 39) {
 304                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002092570028260400000"; mesonCutArray[0] = "01521065000000"; //variation psi pair 0.05
 305                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002092570028280400000"; mesonCutArray[1] = "01521065000000"; //variation psi pair 0.2
 306                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002092570028250000000"; mesonCutArray[2] = "01521065000000"; //variation cosPA -1
 307                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002092570028250600000"; mesonCutArray[3] = "01521065000000"; //variation cosPA 0.9
 308         } else if (trainConfig == 40) { //added signals
 309                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002092570028260400000"; mesonCutArray[0] = "01521065000000"; //variation psi pair 0.05
 310                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002092570028280400000"; mesonCutArray[1] = "01521065000000"; //variation psi pair 0.2
 311                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002092570028250000000"; mesonCutArray[2] = "01521065000000"; //variation cosPA -1
 312                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002092570028250600000"; mesonCutArray[3] = "01521065000000"; //variation cosPA 0.9
 313         } else if (trainConfig == 41) {
 314                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002092570028950400000"; mesonCutArray[0] = "01521065000000"; //variation chi2 15
 315                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002092570028230400000"; mesonCutArray[1] = "01521065000000"; //variation psi pair 0.035
 316                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521055000000"; //variation alpha 0.75
 317                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521075000000"; //variation alpha 0.85
 318         } else if (trainConfig == 42) { //added signals
 319                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002092570028950400000"; mesonCutArray[0] = "01521065000000"; //variation chi2 15
 320                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002092570028230400000"; mesonCutArray[1] = "01521065000000"; //variation psi pair 0.035
 321                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521055000000"; //variation alpha 0.75
 322                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521075000000"; //variation alpha 0.85
 323         } else if (trainConfig == 43) {
 324                 eventCutArray[ 0] = "0004011"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; // trigger kTRD with y 0.8
 325                 eventCutArray[ 1] = "0005011"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[1] = "01521065000000"; // trigger kEMC with y 0.8
 326                 eventCutArray[ 2] = "0006011"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521065000000"; // trigger kPHI with y 0.8
 327                 eventCutArray[ 3] = "0007011"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521065000000"; // trigger kHighMult with y 0.8
 328         } else if (trainConfig == 44) {
 329                 eventCutArray[ 0] = "0008011"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; // trigger kEMCEGA with y 0.8
 330                 eventCutArray[ 1] = "0009011"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[1] = "01521065000000"; // trigger kEMCEJE with y 0.8
 331                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521065000000"; // minimum bias with y 0.8
 332                 eventCutArray[ 3] = "0001111"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521065000000"; // trigger kINT8 with y 0.8
 333         } else if (trainConfig == 45) {
 334                 eventCutArray[ 0] = "0004211"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; // trigger kTRD CINT8 HEE
 335                 eventCutArray[ 1] = "0004411"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[1] = "01521065000000"; // trigger kTRD CINT8 HSE
 336                 eventCutArray[ 2] = "0004611"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521065000000"; // trigger kTRD CINT8 HJE
 337                 eventCutArray[ 3] = "0004811"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521065000000"; // trigger kTRD CINT8 HQU
 338         } else if (trainConfig == 46) {
 339                 eventCutArray[ 0] = "0004111"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; // trigger kTRD CINT7 HEE
 340                 eventCutArray[ 1] = "0004311"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[1] = "01521065000000"; // trigger kTRD CINT7 HSE
 341                 eventCutArray[ 2] = "0004511"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521065000000"; // trigger kTRD CINT7 HJE
 342                 eventCutArray[ 3] = "0004711"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521065000000"; // trigger kTRD CINT7 HQU
 343         } else if (trainConfig == 47) {
 344                 eventCutArray[ 0] = "0005211"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; // trigger kEMC7
 345                 eventCutArray[ 1] = "0005311"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[1] = "01521065000000"; // trigger kEMC8
 346                 eventCutArray[ 2] = "0006211"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521065000000"; // trigger kPHI7
 347                 eventCutArray[ 3] = "0006311"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521065000000"; // trigger kPHI8
 348         } else if (trainConfig == 48) {
 349                 eventCutArray[ 0] = "0005111"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; // trigger kEMC1
 350                 eventCutArray[ 1] = "0007111"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[1] = "01521065000000"; // trigger kSHM1
 351                 eventCutArray[ 2] = "0007211"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521065000000"; // trigger kSHM7
 352                 eventCutArray[ 3] = "0007311"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521065000000"; // trigger kSHM8
 353         } else if (trainConfig == 49) {
 354                 eventCutArray[ 0] = "0008111"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; // trigger kEMCEGA + CINT7
 355                 eventCutArray[ 1] = "0008211"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[1] = "01521065000000"; // trigger kEMCEGA + CINT8
 356                 eventCutArray[ 2] = "0008311"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521065000000"; // trigger kEMCEG1 + CINT7
 357                 eventCutArray[ 3] = "0008411"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521065000000"; // trigger kEMCEG1 + CINT8
 358         } else if (trainConfig == 50) {
 359                 eventCutArray[ 0] = "0008511"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; // trigger kEMCEG2 + CINT7
 360                 eventCutArray[ 1] = "0008611"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[1] = "01521065000000"; // trigger kEMCEG2 + CINT8
 361                 eventCutArray[ 2] = "0009111"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521065000000"; // trigger kEMCEJE + CINT7
 362                 eventCutArray[ 3] = "0009211"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521065000000"; // trigger kEMCEJE + CINT8
 363         } else if (trainConfig == 51) {
 364                 eventCutArray[ 0] = "0009311"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; // trigger kEMCEJ1 + CINT7
 365                 eventCutArray[ 1] = "0009411"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[1] = "01521065000000"; // trigger kEMCEJ1 + CINT8
 366                 eventCutArray[ 2] = "0009511"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521065000000"; // trigger kEMCEJ2 + CINT7
 367                 eventCutArray[ 3] = "0009611"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521065000000"; // trigger kEMCEJ2 + CINT8
 368         } else if (trainConfig == 52) { //pp 2.76TeV cuts
 369                 eventCutArray[ 0] = "0004211"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; // trigger kTRD CINT8 HEE
 370                 eventCutArray[ 1] = "0004411"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; // trigger kTRD CINT8 HSE
 371                 eventCutArray[ 2] = "0004611"; photonCutArray[ 2] = "002093663003800000000"; mesonCutArray[2] = "01631031009000"; // trigger kTRD CINT8 HJE
 372                 eventCutArray[ 3] = "0004811"; photonCutArray[ 3] = "002093663003800000000"; mesonCutArray[3] = "01631031009000"; // trigger kTRD CINT8 HQU
 373         } else if (trainConfig == 53) { //pp 2.76TeV cuts
 374                 eventCutArray[ 0] = "0004111"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; // trigger kTRD CINT7 HEE
 375                 eventCutArray[ 1] = "0004311"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; // trigger kTRD CINT7 HSE
 376                 eventCutArray[ 2] = "0004511"; photonCutArray[ 2] = "002093663003800000000"; mesonCutArray[2] = "01631031009000"; // trigger kTRD CINT7 HJE
 377                 eventCutArray[ 3] = "0004711"; photonCutArray[ 3] = "002093663003800000000"; mesonCutArray[3] = "01631031009000"; // trigger kTRD CINT7 HQU
 378         } else if (trainConfig == 54) { //pp 2.76TeV cuts
 379                 eventCutArray[ 0] = "0005211"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; // trigger kEMC7
 380                 eventCutArray[ 1] = "0005311"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; // trigger kEMC8
 381                 eventCutArray[ 2] = "0006211"; photonCutArray[ 2] = "002093663003800000000"; mesonCutArray[2] = "01631031009000"; // trigger kPHI7
 382                 eventCutArray[ 3] = "0006311"; photonCutArray[ 3] = "002093663003800000000"; mesonCutArray[3] = "01631031009000"; // trigger kPHI8
 383         } else if (trainConfig == 55) { //pp 2.76TeV cuts
 384                 eventCutArray[ 0] = "0005111"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; // trigger kEMC1
 385                 eventCutArray[ 1] = "0007111"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; // trigger kSHM1
 386                 eventCutArray[ 2] = "0007211"; photonCutArray[ 2] = "002093663003800000000"; mesonCutArray[2] = "01631031009000"; // trigger kSHM7
 387                 eventCutArray[ 3] = "0007311"; photonCutArray[ 3] = "002093663003800000000"; mesonCutArray[3] = "01631031009000"; // trigger kSHM8
 388         } else if (trainConfig == 56) { //pp 2.76TeV cuts
 389                 eventCutArray[ 0] = "0008111"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; // trigger kEMCEGA + CINT7
 390                 eventCutArray[ 1] = "0008211"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; // trigger kEMCEGA + CINT8
 391                 eventCutArray[ 2] = "0008311"; photonCutArray[ 2] = "002093663003800000000"; mesonCutArray[2] = "01631031009000"; // trigger kEMCEG1 + CINT7
 392                 eventCutArray[ 3] = "0008411"; photonCutArray[ 3] = "002093663003800000000"; mesonCutArray[3] = "01631031009000"; // trigger kEMCEG1 + CINT8
 393         } else if (trainConfig == 57) { //pp 2.76TeV cuts
 394                 eventCutArray[ 0] = "0008511"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; // trigger kEMCEG2 + CINT7
 395                 eventCutArray[ 1] = "0008611"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; // trigger kEMCEG2 + CINT8
 396                 eventCutArray[ 2] = "0009111"; photonCutArray[ 2] = "002093663003800000000"; mesonCutArray[2] = "01631031009000"; // trigger kEMCEJE + CINT7
 397                 eventCutArray[ 3] = "0009211"; photonCutArray[ 3] = "002093663003800000000"; mesonCutArray[3] = "01631031009000"; // trigger kEMCEJE + CINT8
 398         } else if (trainConfig == 58) { //pp 2.76TeV cuts
 399                 eventCutArray[ 0] = "0009311"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; // trigger kEMCEJ1 + CINT7
 400                 eventCutArray[ 1] = "0009411"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; // trigger kEMCEJ1 + CINT8
 401                 eventCutArray[ 2] = "0009511"; photonCutArray[ 2] = "002093663003800000000"; mesonCutArray[2] = "01631031009000"; // trigger kEMCEJ2 + CINT7
 402                 eventCutArray[ 3] = "0009611"; photonCutArray[ 3] = "002093663003800000000"; mesonCutArray[3] = "01631031009000"; // trigger kEMCEJ2 + CINT8
 403
 404
 405         } else {
 406                         Error(Form("GammaConvV1_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
 407                         return;
 408         }
 409
 410         TList *EventCutList = new TList();
 411         TList *ConvCutList = new TList();
 412         TList *MesonCutList = new TList();
 413
 414         TList *HeaderList = new TList();
 415         if (periodname.CompareTo("LHC12i3") == 0){
 416                 TObjString *Header2 = new TObjString("BOX");
 417                 HeaderList->Add(Header2);
 418         } else if (periodname.CompareTo("LHC14e2b")==0){
 419                 TObjString *Header2 = new TObjString("pi0_1");
 420                 HeaderList->Add(Header2);
 421                 TObjString *Header3 = new TObjString("eta_2");
 422                 HeaderList->Add(Header3);
 423         }
 424
 425         EventCutList->SetOwner(kTRUE);
 426         AliConvEventCuts **analysisEventCuts = new AliConvEventCuts*[numberOfCuts];
 427         ConvCutList->SetOwner(kTRUE);
 428         AliConversionPhotonCuts **analysisCuts = new AliConversionPhotonCuts*[numberOfCuts];
 429         MesonCutList->SetOwner(kTRUE);
 430         AliConversionMesonCuts **analysisMesonCuts = new AliConversionMesonCuts*[numberOfCuts];
 431
 432
 433         for(Int_t i = 0; i<numberOfCuts; i++){
 434                 analysisEventCuts[i] = new AliConvEventCuts();
 435                 analysisEventCuts[i]->InitializeCutsFromCutString(eventCutArray[i].Data());
 436                 EventCutList->Add(analysisEventCuts[i]);
 437                 analysisEventCuts[i]->SetFillCutHistograms("",kFALSE);
 438
 439                 analysisCuts[i] = new AliConversionPhotonCuts();
 440                 analysisCuts[i]->InitializeCutsFromCutString(photonCutArray[i].Data());
 441                 ConvCutList->Add(analysisCuts[i]);
 442                 analysisCuts[i]->SetFillCutHistograms("",kFALSE);
 443
 444                 analysisMesonCuts[i] = new AliConversionMesonCuts();
 445                 analysisMesonCuts[i]->InitializeCutsFromCutString(mesonCutArray[i].Data());
 446                 MesonCutList->Add(analysisMesonCuts[i]);
 447                 analysisMesonCuts[i]->SetFillCutHistograms("");
 448                 analysisEventCuts[i]->SetAcceptedHeader(HeaderList);
 449         }
 450
 451         task->SetEventCutList(numberOfCuts,EventCutList);
 452         task->SetConversionCutList(numberOfCuts,ConvCutList);
 453         task->SetMesonCutList(numberOfCuts,MesonCutList);
 454         task->SetMoveParticleAccordingToVertex(kTRUE);
 455         task->SetDoMesonAnalysis(kTRUE);
 456         task->SetDoMesonQA(enableQAMesonTask); //Attention new switch for Pi0 QA
 457         task->SetDoPhotonQA(enableQAPhotonTask);  //Attention new switch small for Photon QA
 458
 459         //connect containers
 460         AliAnalysisDataContainer *coutput =
 461                 mgr->CreateContainer(Form("GammaConvV1_%i",trainConfig), TList::Class(),
 462                                                         AliAnalysisManager::kOutputContainer,Form("GammaConvV1_%i.root",trainConfig));
 463
 464         mgr->AddTask(task);
 465         mgr->ConnectInput(task,0,cinput);
 466         mgr->ConnectOutput(task,1,coutput);
 467
 468         return;
 469
 470 }