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         0000011_01631031009000
 114
 115         if(trainConfig == 1){
 116                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , only Minbias MC
 117                 eventCutArray[ 1] = "0001011"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD, V0AND
 118                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002093260003800000000"; mesonCutArray[2] = "01631031009000"; //standard cut Gamma pp 2-76TeV
 119                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002093660003800000000"; mesonCutArray[3] = "01631031009000"; //standard cut Gamma pp 2-76TeV
 120         } else if (trainConfig == 2) {
 121                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , only boxes
 122                 eventCutArray[ 1] = "0001012"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD, V0AND , only boxes
 123                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002093260003800000000"; mesonCutArray[2] = "01631031009000"; //standard cut Gamma pp 2-76TeV , only boxes
 124                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002093660003800000000"; mesonCutArray[3] = "01631031009000"; //standard cut Gamma pp 2-76TeV
 125         } else if (trainConfig == 3) {
 126                 eventCutArray[ 0] = "0002011"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , only Minbias MC
 127                 eventCutArray[ 1] = "0003011"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD, V0AND , only Minbias MC
 128                 eventCutArray[ 2] = "0002012"; photonCutArray[ 2] = "002093663003800000000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , only Boxes MC
 129                 eventCutArray[ 3] = "0003012"; photonCutArray[ 3] = "002093663003800000000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD, V0AND, only Boxes MC
 130         } else if (trainConfig == 4) {
 131                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , all photon qualities
 132                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002093663003800020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 1
 133                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002093663003800030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 2
 134                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002093663003800040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 3
 135         } else if (trainConfig == 5) {
 136                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "007093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , all photon qualities, min R = 35 cm
 137                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "007093663003800020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 1, min R = 35 cm
 138                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "007093663003800030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 2, min R = 35 cm
 139                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "007093663003800040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 3, min R = 35 cm
 140         } else if (trainConfig == 6) {
 141                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002083663003200000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, all photon qualities
 142                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002083663003200020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 1
 143                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002083663003200030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 2
 144                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002083663003200040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 3
 145         } else if (trainConfig == 7) {
 146                 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
 147                 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
 148                 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
 149                 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
 150         } else if (trainConfig == 8) {
 151                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002083663000200000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 7TeV, all photon qualities
 152                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002083663000200020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 1
 153                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002083663000200030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 2
 154                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002083663000200040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 3
 155         } else if (trainConfig == 9) {
 156                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "007083663000200000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 7TeV, all photon qualities, min R = 35 cm
 157                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "007083663000200020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 1, min R = 35 cm
 158                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "007083663000200030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 2, min R = 35 cm
 159                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "007083663000200040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 3, min R = 35 cm
 160         } else if (trainConfig == 10) {
 161                 eventCutArray[ 0] = "0002011"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , all photon qualities
 162                 eventCutArray[ 1] = "0002011"; photonCutArray[ 1] = "002093663003800020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 1
 163                 eventCutArray[ 2] = "0002011"; photonCutArray[ 2] = "002093663003800030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 2
 164                 eventCutArray[ 3] = "0002011"; photonCutArray[ 3] = "002093663003800040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 3
 165         } else if (trainConfig == 11) {
 166                 eventCutArray[ 0] = "0002011"; photonCutArray[ 0] = "007093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , all photon qualities, min R = 35 cm
 167                 eventCutArray[ 1] = "0002011"; photonCutArray[ 1] = "007093663003800020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 1, min R = 35 cm
 168                 eventCutArray[ 2] = "0002011"; photonCutArray[ 2] = "007093663003800030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 2, min R = 35 cm
 169                 eventCutArray[ 3] = "0002011"; photonCutArray[ 3] = "007093663003800040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 3, min R = 35 cm
 170         } else if (trainConfig == 12) {
 171                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; //standard cut LHC11h pp 2.76TeV
 172                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "032092970028250400000"; mesonCutArray[1] = "01525065000000"; //variation eta 0.65
 173                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "042092970028250400000"; mesonCutArray[2] = "01525065000000"; //variation eta 0.75
 174                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002092950028250400000"; mesonCutArray[3] = "01525065000000"; //variation pion p dEdx 0.3-5.
 175         } else if (trainConfig == 13) { //added signals
 176                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; //standard cut LHC11h pp 2.76TeV
 177                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "032092970028250400000"; mesonCutArray[1] = "01525065000000"; //variation eta 0.65
 178                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "042092970028250400000"; mesonCutArray[2] = "01525065000000"; //variation eta 0.75
 179                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002092950028250400000"; mesonCutArray[3] = "01525065000000"; //variation pion p dEdx 0.3-5.
 180         } else if (trainConfig == 14) {
 181                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002492970028250400000"; mesonCutArray[0] = "01525065000000"; //variation pt 0.075
 182                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002192970028250400000"; mesonCutArray[1] = "01525065000000"; //variation pt 0.1
 183                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002062970028250400000"; mesonCutArray[2] = "01525065000000"; //variation TPC cls 0.7
 184                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002082970028250400000"; mesonCutArray[3] = "01525065000000"; //variation TPC cls 0.35
 185         } else if (trainConfig == 15) { //added signals
 186                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002492970028250400000"; mesonCutArray[0] = "01525065000000"; //variation pt 0.075
 187                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002192970028250400000"; mesonCutArray[1] = "01525065000000"; //variation pt 0.1
 188                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002062970028250400000"; mesonCutArray[2] = "01525065000000"; //variation TPC cls 0.7
 189                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002082970028250400000"; mesonCutArray[3] = "01525065000000"; //variation TPC cls 0.35
 190         } else if (trainConfig == 16) {
 191                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002093970028250400000"; mesonCutArray[0] = "01525065000000"; //variation edEdx -4,5
 192                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002096970028250400000"; mesonCutArray[1] = "01525065000000"; //variation edEdx -2.5,4
 193                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002092970038250400000"; mesonCutArray[2] = "01525065000000"; //variation TOF el. PID -3,5
 194                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002092970048250400000"; mesonCutArray[3] = "01525065000000"; //variation TOF el. PID -2,3
 195         } else if (trainConfig == 17) { //added signals
 196                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002093970028250400000"; mesonCutArray[0] = "01525065000000"; //variation edEdx -4,5
 197                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002096970028250400000"; mesonCutArray[1] = "01525065000000"; //variation edEdx -2.5,4
 198                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002092970038250400000"; mesonCutArray[2] = "01525065000000"; //variation TOF el. PID -3,5
 199                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002092970048250400000"; mesonCutArray[3] = "01525065000000"; //variation TOF el. PID -2,3
 200         } else if (trainConfig == 18) {
 201                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002092970029250400000"; mesonCutArray[0] = "01525065000000"; //variation qt 0.03
 202                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002092970022250400000"; mesonCutArray[1] = "01525065000000"; //variation qt 0.07 no2D
 203                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002092970028150400000"; mesonCutArray[2] = "01525065000000"; //variation chi2 50.
 204                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002092970028850400000"; mesonCutArray[3] = "01525065000000"; //variation chi2 20.
 205         } else if (trainConfig == 19) { //added signals
 206                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002092970029250400000"; mesonCutArray[0] = "01525065000000"; //variation qt 0.03
 207                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002092970022250400000"; mesonCutArray[1] = "01525065000000"; //variation qt 0.07 no2D
 208                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002092970028150400000"; mesonCutArray[2] = "01525065000000"; //variation chi2 50.
 209                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002092970028850400000"; mesonCutArray[3] = "01525065000000"; //variation chi2 20.
 210         } else if (trainConfig == 20) {
 211                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002092970028260400000"; mesonCutArray[0] = "01525065000000"; //variation psi pair 0.05
 212                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002092970028280400000"; mesonCutArray[1] = "01525065000000"; //variation psi pair 0.2
 213                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002092970028250000000"; mesonCutArray[2] = "01525065000000"; //variation cosPA -1
 214                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525055000000"; //variation alpha 0.75
 215         } else if (trainConfig == 21) { //added signals
 216                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002092970028260400000"; mesonCutArray[0] = "01525065000000"; //variation psi pair 0.05
 217                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002092970028280400000"; mesonCutArray[1] = "01525065000000"; //variation psi pair 0.2
 218                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002092970028250000000"; mesonCutArray[2] = "01525065000000"; //variation cosPA -1
 219                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525055000000"; //variation alpha 0.75
 220         } else if (trainConfig == 22) {
 221                 eventCutArray[ 0] = "0004011"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; // trigger kTRD
 222                 eventCutArray[ 1] = "0005011"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[1] = "01525065000000"; // trigger kEMC
 223                 eventCutArray[ 2] = "0006011"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01525065000000"; // trigger kPHI
 224                 eventCutArray[ 3] = "0007011"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525065000000"; // trigger kHighMult
 225         } else if (trainConfig == 23) {
 226                 eventCutArray[ 0] = "0008011"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; // trigger kEMCEGA
 227                 eventCutArray[ 1] = "0009011"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[1] = "01525065000000"; // trigger kEMCEJE
 228                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01525065000000"; // minimum bias
 229                 eventCutArray[ 3] = "0001111"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525065000000"; // trigger kINT8
 230         } else if (trainConfig == 24) {
 231                 eventCutArray[ 0] = "0004211"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; // trigger kTRD CINT8 HEE
 232                 eventCutArray[ 1] = "0004411"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[1] = "01525065000000"; // trigger kTRD CINT8 HSE
 233                 eventCutArray[ 2] = "0004611"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01525065000000"; // trigger kTRD CINT8 HJE
 234                 eventCutArray[ 3] = "0004811"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525065000000"; // trigger kTRD CINT8 HQU
 235         } else if (trainConfig == 25) {
 236                 eventCutArray[ 0] = "0004111"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; // trigger kTRD CINT7 HEE
 237                 eventCutArray[ 1] = "0004311"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[1] = "01525065000000"; // trigger kTRD CINT7 HSE
 238                 eventCutArray[ 2] = "0004511"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01525065000000"; // trigger kTRD CINT7 HJE
 239                 eventCutArray[ 3] = "0004711"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525065000000"; // trigger kTRD CINT7 HQU
 240         } else if (trainConfig == 26) {
 241                 eventCutArray[ 0] = "0005211"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; // trigger kEMC7
 242                 eventCutArray[ 1] = "0005311"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[1] = "01525065000000"; // trigger kEMC8
 243                 eventCutArray[ 2] = "0006211"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01525065000000"; // trigger kPHI7
 244                 eventCutArray[ 3] = "0006311"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525065000000"; // trigger kPHI8
 245         } else if (trainConfig == 27) {
 246                 eventCutArray[ 0] = "0005111"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; // trigger kEMC1
 247                 eventCutArray[ 1] = "0007111"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[1] = "01525065000000"; // trigger kSHM1
 248                 eventCutArray[ 2] = "0007211"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01525065000000"; // trigger kSHM7
 249                 eventCutArray[ 3] = "0007311"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525065000000"; // trigger kSHM8
 250         } else if (trainConfig == 28) {
 251                 eventCutArray[ 0] = "0008111"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; // trigger kEMCEGA + CINT7
 252                 eventCutArray[ 1] = "0008211"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[1] = "01525065000000"; // trigger kEMCEGA + CINT8
 253                 eventCutArray[ 2] = "0008311"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01525065000000"; // trigger kEMCEG1 + CINT7
 254                 eventCutArray[ 3] = "0008411"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525065000000"; // trigger kEMCEG1 + CINT8
 255         } else if (trainConfig == 29) {
 256                 eventCutArray[ 0] = "0008511"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; // trigger kEMCEG2 + CINT7
 257                 eventCutArray[ 1] = "0008611"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[1] = "01525065000000"; // trigger kEMCEG2 + CINT8
 258                 eventCutArray[ 2] = "0009111"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01525065000000"; // trigger kEMCEJE + CINT7
 259                 eventCutArray[ 3] = "0009211"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525065000000"; // trigger kEMCEJE + CINT8
 260         } else if (trainConfig == 30) {
 261                 eventCutArray[ 0] = "0009311"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[0] = "01525065000000"; // trigger kEMCEJ1 + CINT7
 262                 eventCutArray[ 1] = "0009411"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[1] = "01525065000000"; // trigger kEMCEJ1 + CINT8
 263                 eventCutArray[ 2] = "0009511"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01525065000000"; // trigger kEMCEJ2 + CINT7
 264                 eventCutArray[ 3] = "0009611"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[3] = "01525065000000"; // trigger kEMCEJ2 + CINT8
 265         } else if (trainConfig == 31) {
 266                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; //new standard cut for pp 8 TeV
 267                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002093570028250400000"; mesonCutArray[1] = "01521065000000"; //variation edEdx -4,5
 268                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002096570028250400000"; mesonCutArray[2] = "01521065000000"; //variation edEdx -2.5,4
 269                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002092550028250400000"; mesonCutArray[3] = "01521065000000"; //variation pion p dEdx 0.3-5.
 270         } else if (trainConfig == 32) { //added signals
 271                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; //new standard cut for pp 8 TeV
 272                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002093570028250400000"; mesonCutArray[1] = "01521065000000"; //variation edEdx -4,5
 273                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002096570028250400000"; mesonCutArray[2] = "01521065000000"; //variation edEdx -2.5,4
 274                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002092550028250400000"; mesonCutArray[3] = "01521065000000"; //variation pion p dEdx 0.3-5.
 275         } else if (trainConfig == 33) {
 276                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002492570028250400000"; mesonCutArray[0] = "01521065000000"; //variation pt 0.075
 277                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002192570028250400000"; mesonCutArray[1] = "01521065000000"; //variation pt 0.1
 278                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002062570028250400000"; mesonCutArray[2] = "01521065000000"; //variation TPC cls 0.7
 279                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002082570028250400000"; mesonCutArray[3] = "01521065000000"; //variation TPC cls 0.35
 280         } else if (trainConfig == 34) { //added signals
 281                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002492570028250400000"; mesonCutArray[0] = "01521065000000"; //variation pt 0.075
 282                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002192570028250400000"; mesonCutArray[1] = "01521065000000"; //variation pt 0.1
 283                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002062570028250400000"; mesonCutArray[2] = "01521065000000"; //variation TPC cls 0.7
 284                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002082570028250400000"; mesonCutArray[3] = "01521065000000"; //variation TPC cls 0.35
 285         } else if (trainConfig == 35) {
 286                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002092270028250400000"; mesonCutArray[0] = "01521065000000"; //variation pidEdx 1,-10
 287                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002092370028250400000"; mesonCutArray[1] = "01521065000000"; //variation pidEdx 2.5,-10
 288                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01521065000000"; //variation pidEdx 3,-10
 289                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002092500028250400000"; mesonCutArray[3] = "01521065000000"; //variation pion p dEdx 0.5-5
 290         } else if (trainConfig == 36) { //added signals
 291                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002092270028250400000"; mesonCutArray[0] = "01521065000000"; //variation pidEdx 1,-10
 292                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002092370028250400000"; mesonCutArray[1] = "01521065000000"; //variation pidEdx 2.5,-10
 293                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[2] = "01521065000000"; //variation pidEdx 3,-10
 294                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002092500028250400000"; mesonCutArray[3] = "01521065000000"; //variation pion p dEdx 0.5-5
 295         } else if (trainConfig == 37) {
 296                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002092570029250400000"; mesonCutArray[0] = "01521065000000"; //variation qt 0.03
 297                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002092570022250400000"; mesonCutArray[1] = "01521065000000"; //variation qt 0.07 no2D
 298                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002092570028150400000"; mesonCutArray[2] = "01521065000000"; //variation chi2 50.
 299                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002092570028850400000"; mesonCutArray[3] = "01521065000000"; //variation chi2 20.
 300         } else if (trainConfig == 38) { //added signals
 301                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002092570029250400000"; mesonCutArray[0] = "01521065000000"; //variation qt 0.03
 302                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002092570022250400000"; mesonCutArray[1] = "01521065000000"; //variation qt 0.07 no2D
 303                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002092570028150400000"; mesonCutArray[2] = "01521065000000"; //variation chi2 50.
 304                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002092570028850400000"; mesonCutArray[3] = "01521065000000"; //variation chi2 20.
 305         } else if (trainConfig == 39) {
 306                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002092570028260400000"; mesonCutArray[0] = "01521065000000"; //variation psi pair 0.05
 307                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002092570028280400000"; mesonCutArray[1] = "01521065000000"; //variation psi pair 0.2
 308                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002092570028250000000"; mesonCutArray[2] = "01521065000000"; //variation cosPA -1
 309                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002092570028250600000"; mesonCutArray[3] = "01521065000000"; //variation cosPA 0.9
 310         } else if (trainConfig == 40) { //added signals
 311                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002092570028260400000"; mesonCutArray[0] = "01521065000000"; //variation psi pair 0.05
 312                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002092570028280400000"; mesonCutArray[1] = "01521065000000"; //variation psi pair 0.2
 313                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002092570028250000000"; mesonCutArray[2] = "01521065000000"; //variation cosPA -1
 314                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002092570028250600000"; mesonCutArray[3] = "01521065000000"; //variation cosPA 0.9
 315         } else if (trainConfig == 41) {
 316                 eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002092570028950400000"; mesonCutArray[0] = "01521065000000"; //variation chi2 15
 317                 eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002092570028230400000"; mesonCutArray[1] = "01521065000000"; //variation psi pair 0.035
 318                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521055000000"; //variation alpha 0.75
 319                 eventCutArray[ 3] = "0000011"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521075000000"; //variation alpha 0.85
 320         } else if (trainConfig == 42) { //added signals
 321                 eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002092570028950400000"; mesonCutArray[0] = "01521065000000"; //variation chi2 15
 322                 eventCutArray[ 1] = "0000012"; photonCutArray[ 1] = "002092570028230400000"; mesonCutArray[1] = "01521065000000"; //variation psi pair 0.035
 323                 eventCutArray[ 2] = "0000012"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521055000000"; //variation alpha 0.75
 324                 eventCutArray[ 3] = "0000012"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521075000000"; //variation alpha 0.85
 325         } else if (trainConfig == 43) {
 326                 eventCutArray[ 0] = "0004011"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; // trigger kTRD with y 0.8
 327                 eventCutArray[ 1] = "0005011"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[1] = "01521065000000"; // trigger kEMC with y 0.8
 328                 eventCutArray[ 2] = "0006011"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521065000000"; // trigger kPHI with y 0.8
 329                 eventCutArray[ 3] = "0007011"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521065000000"; // trigger kHighMult with y 0.8
 330         } else if (trainConfig == 44) {
 331                 eventCutArray[ 0] = "0008011"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; // trigger kEMCEGA with y 0.8
 332                 eventCutArray[ 1] = "0009011"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[1] = "01521065000000"; // trigger kEMCEJE with y 0.8
 333                 eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521065000000"; // minimum bias with y 0.8
 334                 eventCutArray[ 3] = "0001111"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521065000000"; // trigger kINT8 with y 0.8
 335         } else if (trainConfig == 45) {
 336                 eventCutArray[ 0] = "0004211"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; // trigger kTRD CINT8 HEE
 337                 eventCutArray[ 1] = "0004411"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[1] = "01521065000000"; // trigger kTRD CINT8 HSE
 338                 eventCutArray[ 2] = "0004611"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521065000000"; // trigger kTRD CINT8 HJE
 339                 eventCutArray[ 3] = "0004811"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521065000000"; // trigger kTRD CINT8 HQU
 340         } else if (trainConfig == 46) {
 341                 eventCutArray[ 0] = "0004111"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; // trigger kTRD CINT7 HEE
 342                 eventCutArray[ 1] = "0004311"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[1] = "01521065000000"; // trigger kTRD CINT7 HSE
 343                 eventCutArray[ 2] = "0004511"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521065000000"; // trigger kTRD CINT7 HJE
 344                 eventCutArray[ 3] = "0004711"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521065000000"; // trigger kTRD CINT7 HQU
 345         } else if (trainConfig == 47) {
 346                 eventCutArray[ 0] = "0005211"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; // trigger kEMC7
 347                 eventCutArray[ 1] = "0005311"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[1] = "01521065000000"; // trigger kEMC8
 348                 eventCutArray[ 2] = "0006211"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521065000000"; // trigger kPHI7
 349                 eventCutArray[ 3] = "0006311"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521065000000"; // trigger kPHI8
 350         } else if (trainConfig == 48) {
 351                 eventCutArray[ 0] = "0005111"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; // trigger kEMC1
 352                 eventCutArray[ 1] = "0007111"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[1] = "01521065000000"; // trigger kSHM1
 353                 eventCutArray[ 2] = "0007211"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521065000000"; // trigger kSHM7
 354                 eventCutArray[ 3] = "0007311"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521065000000"; // trigger kSHM8
 355         } else if (trainConfig == 49) {
 356                 eventCutArray[ 0] = "0008111"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; // trigger kEMCEGA + CINT7
 357                 eventCutArray[ 1] = "0008211"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[1] = "01521065000000"; // trigger kEMCEGA + CINT8
 358                 eventCutArray[ 2] = "0008311"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521065000000"; // trigger kEMCEG1 + CINT7
 359                 eventCutArray[ 3] = "0008411"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521065000000"; // trigger kEMCEG1 + CINT8
 360         } else if (trainConfig == 50) {
 361                 eventCutArray[ 0] = "0008511"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; // trigger kEMCEG2 + CINT7
 362                 eventCutArray[ 1] = "0008611"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[1] = "01521065000000"; // trigger kEMCEG2 + CINT8
 363                 eventCutArray[ 2] = "0009111"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521065000000"; // trigger kEMCEJE + CINT7
 364                 eventCutArray[ 3] = "0009211"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521065000000"; // trigger kEMCEJE + CINT8
 365         } else if (trainConfig == 51) {
 366                 eventCutArray[ 0] = "0009311"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[0] = "01521065000000"; // trigger kEMCEJ1 + CINT7
 367                 eventCutArray[ 1] = "0009411"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[1] = "01521065000000"; // trigger kEMCEJ1 + CINT8
 368                 eventCutArray[ 2] = "0009511"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[2] = "01521065000000"; // trigger kEMCEJ2 + CINT7
 369                 eventCutArray[ 3] = "0009611"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[3] = "01521065000000"; // trigger kEMCEJ2 + CINT8
 370         } else if (trainConfig == 52) { //pp 2.76TeV cuts
 371                 eventCutArray[ 0] = "0004211"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; // trigger kTRD CINT8 HEE
 372                 eventCutArray[ 1] = "0004411"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; // trigger kTRD CINT8 HSE
 373                 eventCutArray[ 2] = "0004611"; photonCutArray[ 2] = "002093663003800000000"; mesonCutArray[2] = "01631031009000"; // trigger kTRD CINT8 HJE
 374                 eventCutArray[ 3] = "0004811"; photonCutArray[ 3] = "002093663003800000000"; mesonCutArray[3] = "01631031009000"; // trigger kTRD CINT8 HQU
 375         } else if (trainConfig == 53) { //pp 2.76TeV cuts
 376                 eventCutArray[ 0] = "0004111"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; // trigger kTRD CINT7 HEE
 377                 eventCutArray[ 1] = "0004311"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; // trigger kTRD CINT7 HSE
 378                 eventCutArray[ 2] = "0004511"; photonCutArray[ 2] = "002093663003800000000"; mesonCutArray[2] = "01631031009000"; // trigger kTRD CINT7 HJE
 379                 eventCutArray[ 3] = "0004711"; photonCutArray[ 3] = "002093663003800000000"; mesonCutArray[3] = "01631031009000"; // trigger kTRD CINT7 HQU
 380         } else if (trainConfig == 54) { //pp 2.76TeV cuts
 381                 eventCutArray[ 0] = "0005211"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; // trigger kEMC7
 382                 eventCutArray[ 1] = "0005311"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; // trigger kEMC8
 383                 eventCutArray[ 2] = "0006211"; photonCutArray[ 2] = "002093663003800000000"; mesonCutArray[2] = "01631031009000"; // trigger kPHI7
 384                 eventCutArray[ 3] = "0006311"; photonCutArray[ 3] = "002093663003800000000"; mesonCutArray[3] = "01631031009000"; // trigger kPHI8
 385         } else if (trainConfig == 55) { //pp 2.76TeV cuts
 386                 eventCutArray[ 0] = "0005111"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; // trigger kEMC1
 387                 eventCutArray[ 1] = "0007111"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; // trigger kSHM1
 388                 eventCutArray[ 2] = "0007211"; photonCutArray[ 2] = "002093663003800000000"; mesonCutArray[2] = "01631031009000"; // trigger kSHM7
 389                 eventCutArray[ 3] = "0007311"; photonCutArray[ 3] = "002093663003800000000"; mesonCutArray[3] = "01631031009000"; // trigger kSHM8
 390         } else if (trainConfig == 56) { //pp 2.76TeV cuts
 391                 eventCutArray[ 0] = "0008111"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; // trigger kEMCEGA + CINT7
 392                 eventCutArray[ 1] = "0008211"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; // trigger kEMCEGA + CINT8
 393                 eventCutArray[ 2] = "0008311"; photonCutArray[ 2] = "002093663003800000000"; mesonCutArray[2] = "01631031009000"; // trigger kEMCEG1 + CINT7
 394                 eventCutArray[ 3] = "0008411"; photonCutArray[ 3] = "002093663003800000000"; mesonCutArray[3] = "01631031009000"; // trigger kEMCEG1 + CINT8
 395         } else if (trainConfig == 57) { //pp 2.76TeV cuts
 396                 eventCutArray[ 0] = "0008511"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; // trigger kEMCEG2 + CINT7
 397                 eventCutArray[ 1] = "0008611"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; // trigger kEMCEG2 + CINT8
 398                 eventCutArray[ 2] = "0009111"; photonCutArray[ 2] = "002093663003800000000"; mesonCutArray[2] = "01631031009000"; // trigger kEMCEJE + CINT7
 399                 eventCutArray[ 3] = "0009211"; photonCutArray[ 3] = "002093663003800000000"; mesonCutArray[3] = "01631031009000"; // trigger kEMCEJE + CINT8
 400         } else if (trainConfig == 58) { //pp 2.76TeV cuts
 401                 eventCutArray[ 0] = "0009311"; photonCutArray[ 0] = "002093663003800000000"; mesonCutArray[0] = "01631031009000"; // trigger kEMCEJ1 + CINT7
 402                 eventCutArray[ 1] = "0009411"; photonCutArray[ 1] = "002093663003800000000"; mesonCutArray[1] = "01631031009000"; // trigger kEMCEJ1 + CINT8
 403                 eventCutArray[ 2] = "0009511"; photonCutArray[ 2] = "002093663003800000000"; mesonCutArray[2] = "01631031009000"; // trigger kEMCEJ2 + CINT7
 404                 eventCutArray[ 3] = "0009611"; photonCutArray[ 3] = "002093663003800000000"; mesonCutArray[3] = "01631031009000"; // trigger kEMCEJ2 + CINT8
 405
 406
 407         } else {
 408                         Error(Form("GammaConvV1_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
 409                         return;
 410         }
 411
 412         TList *EventCutList = new TList();
 413         TList *ConvCutList = new TList();
 414         TList *MesonCutList = new TList();
 415
 416         TList *HeaderList = new TList();
 417         if (periodname.CompareTo("LHC12i3") == 0){
 418                 TObjString *Header2 = new TObjString("BOX");
 419                 HeaderList->Add(Header2);
 420         } else if (periodname.CompareTo("LHC14e2b")==0){
 421                 TObjString *Header2 = new TObjString("pi0_1");
 422                 HeaderList->Add(Header2);
 423                 TObjString *Header3 = new TObjString("eta_2");
 424                 HeaderList->Add(Header3);
 425         }
 426
 427         EventCutList->SetOwner(kTRUE);
 428         AliConvEventCuts **analysisEventCuts = new AliConvEventCuts*[numberOfCuts];
 429         ConvCutList->SetOwner(kTRUE);
 430         AliConversionPhotonCuts **analysisCuts = new AliConversionPhotonCuts*[numberOfCuts];
 431         MesonCutList->SetOwner(kTRUE);
 432         AliConversionMesonCuts **analysisMesonCuts = new AliConversionMesonCuts*[numberOfCuts];
 433
 434
 435         for(Int_t i = 0; i<numberOfCuts; i++){
 436                 analysisEventCuts[i] = new AliConvEventCuts();
 437                 analysisEventCuts[i]->InitializeCutsFromCutString(eventCutArray[i].Data());
 438                 EventCutList->Add(analysisEventCuts[i]);
 439                 analysisEventCuts[i]->SetFillCutHistograms("",kFALSE);
 440
 441                 analysisCuts[i] = new AliConversionPhotonCuts();
 442                 analysisCuts[i]->InitializeCutsFromCutString(photonCutArray[i].Data());
 443                 ConvCutList->Add(analysisCuts[i]);
 444                 analysisCuts[i]->SetFillCutHistograms("",kFALSE);
 445
 446                 analysisMesonCuts[i] = new AliConversionMesonCuts();
 447                 analysisMesonCuts[i]->InitializeCutsFromCutString(mesonCutArray[i].Data());
 448                 MesonCutList->Add(analysisMesonCuts[i]);
 449                 analysisMesonCuts[i]->SetFillCutHistograms("");
 450                 analysisEventCuts[i]->SetAcceptedHeader(HeaderList);
 451         }
 452
 453         task->SetEventCutList(numberOfCuts,EventCutList);
 454         task->SetConversionCutList(numberOfCuts,ConvCutList);
 455         task->SetMesonCutList(numberOfCuts,MesonCutList);
 456         task->SetMoveParticleAccordingToVertex(kTRUE);
 457         task->SetDoMesonAnalysis(kTRUE);
 458         task->SetDoMesonQA(enableQAMesonTask); //Attention new switch for Pi0 QA
 459         task->SetDoPhotonQA(enableQAPhotonTask);  //Attention new switch small for Photon QA
 460
 461         //connect containers
 462         AliAnalysisDataContainer *coutput =
 463                 mgr->CreateContainer(Form("GammaConvV1_%i",trainConfig), TList::Class(),
 464                                                         AliAnalysisManager::kOutputContainer,Form("GammaConvV1_%i.root",trainConfig));
 465
 466         mgr->AddTask(task);
 467         mgr->ConnectInput(task,0,cinput);
 468         mgr->ConnectOutput(task,1,coutput);
 469
 470         return;
 471
 472 }