PWGGA/GammaConv/macros/AddTask_GammaConvV1_pPb.C

   1 void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of cuts
   2                               Bool_t isMC   = kFALSE, //run MC
   3                               Int_t enableQAMesonTask = 0, //enable QA in AliAnalysisTaskGammaConvV1
   4                               Int_t enableQAPhotonTask = 0, // enable additional QA task
   5                               TString fileNameInputForWeighting = "MCSpectraInput.root", // path to file for weigting input
   6                               Int_t doWeightingPart = 0,  //enable Weighting
   7                               TString generatorName = "DPMJET",
   8                               TString cutnumberAODBranch = "8000000060084000001500000" // cutnumber for AOD branch
   9                            ) {
  10
  11    // ================= Load Librariers =================================
  12    gSystem->Load("libCore.so");
  13    gSystem->Load("libTree.so");
  14    gSystem->Load("libGeom.so");
  15    gSystem->Load("libVMC.so");
  16    gSystem->Load("libPhysics.so");
  17    gSystem->Load("libMinuit");
  18    gSystem->Load("libSTEERBase");
  19    gSystem->Load("libESD");
  20    gSystem->Load("libAOD");
  21    gSystem->Load("libANALYSIS");
  22    gSystem->Load("libANALYSISalice");
  23    gSystem->Load("libPWGGAGammaConv.so");
  24    gSystem->Load("libCDB.so");
  25    gSystem->Load("libSTEER.so");
  26    gSystem->Load("libSTEERBase.so");
  27    gSystem->Load("libTENDER.so");
  28    gSystem->Load("libTENDERSupplies.so");
  29
  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 cutnumber = "8000000060084001001500000000";
  48    Bool_t doEtaShift = kFALSE;
  49    AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
  50    //========= Add V0 Reader to  ANALYSIS manager if not yet existent =====
  51    if( !(AliV0ReaderV1*)mgr->GetTask("V0ReaderV1") ){
  52       AliV0ReaderV1 *fV0ReaderV1 = new AliV0ReaderV1("V0ReaderV1");
  53
  54       fV0ReaderV1->SetUseOwnXYZCalculation(kTRUE);
  55       fV0ReaderV1->SetCreateAODs(kFALSE);// AOD Output
  56       fV0ReaderV1->SetUseAODConversionPhoton(kTRUE);
  57
  58       if (!mgr) {
  59          Error("AddTask_V0ReaderV1", "No analysis manager found.");
  60          return;
  61       }
  62
  63       // Set AnalysisCut Number
  64       AliConversionCuts *fCuts=NULL;
  65       if(cutnumber!=""){
  66          fCuts= new AliConversionCuts(cutnumber.Data(),cutnumber.Data());
  67          fCuts->SetPreSelectionCutFlag(kTRUE);
  68          if(fCuts->InitializeCutsFromCutString(cutnumber.Data())){
  69             fCuts->DoEtaShift(doEtaShift);
  70             fV0ReaderV1->SetConversionCuts(fCuts);
  71             fCuts->SetFillCutHistograms("",kTRUE);
  72          }
  73       }
  74       if(inputHandler->IsA()==AliAODInputHandler::Class()){
  75       // AOD mode
  76          cout << "AOD handler: adding " << cutnumberAODBranch.Data() << " as conversion branch" << endl;
  77          fV0ReaderV1->SetDeltaAODBranchName(Form("GammaConv_%s_gamma",cutnumberAODBranch.Data()));
  78       }
  79       fV0ReaderV1->Init();
  80
  81       AliLog::SetGlobalLogLevel(AliLog::kInfo);
  82
  83       //connect input V0Reader
  84       mgr->AddTask(fV0ReaderV1);
  85       mgr->ConnectInput(fV0ReaderV1,0,cinput);
  86
  87    }
  88
  89    //================================================
  90    //========= Add task to the ANALYSIS manager =====
  91    //================================================
  92    //            find input container
  93    AliAnalysisTaskGammaConvV1 *task=NULL;
  94    task= new AliAnalysisTaskGammaConvV1(Form("GammaConvV1_%i",trainConfig));
  95    task->SetIsHeavyIon(2);
  96    task->SetIsMC(isMC);
  97    // Cut Numbers to use in Analysis
  98    Int_t numberOfCuts = 4;
  99
 100    TString *cutarray = new TString[numberOfCuts];
 101    TString *mesonCutArray = new TString[numberOfCuts];
 102    Bool_t doEtaShiftIndCuts = kFALSE;
 103    TString stringShift = "";
 104         if(trainConfig == 1){
 105                 cutarray[ 0] = "8000011032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
 106                 cutarray[ 1] = "8000011042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
 107                 cutarray[ 2] = "8000011012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
 108                 cutarray[ 3] = "8000011002092370023220000000"; mesonCutArray[ 3] = "01621035009000";
 109         } else if (trainConfig == 2) {
 110                 cutarray[ 0] = "8000012032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
 111                 cutarray[ 1] = "8000012042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
 112                 cutarray[ 2] = "8000012012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
 113                 cutarray[ 3] = "8000012002092370023220000000"; mesonCutArray[ 3] = "01621035009000";
 114         } else if (trainConfig == 3) {
 115                 cutarray[ 0] = "8000011002092770023220000000"; mesonCutArray[ 0] = "01621035009000";
 116                 cutarray[ 1] = "8000011002092551023220000000"; mesonCutArray[ 1] = "01621035009000";
 117                 cutarray[ 2] = "8000011002092170003220000000"; mesonCutArray[ 2] = "01621035009000"; //just tighten Psi pair
 118                 cutarray[ 3] = "8000011002092170003260000000"; mesonCutArray[ 3] = "01621035009000"; //tighten Psi pair and chi2 in 2D
 119         } else if (trainConfig == 4) {
 120                 cutarray[ 0] = "8000012002092770023220000000"; mesonCutArray[ 0] = "01621035009000";
 121                 cutarray[ 1] = "8000012002092551023220000000"; mesonCutArray[ 1] = "01621035009000";
 122                 cutarray[ 2] = "8000012002092170003220000000"; mesonCutArray[ 2] = "01621035009000"; //just tighten Psi pair
 123                 cutarray[ 3] = "8000012002092170003260000000"; mesonCutArray[ 3] = "01621035009000"; //tighten Psi pair and chi2 in 2D
 124         } else if (trainConfig == 5) {
 125                 cutarray[ 0] = "8000011002092170008220000000"; mesonCutArray[ 0] = "01621035009000"; //tighten psi pair and qt in 2D
 126                 cutarray[ 1] = "8000011002092170008260000000"; mesonCutArray[ 1] = "01621035009000"; //tighten psi pair and chi2 in 2D and qt in 2D
 127                 cutarray[ 2] = "8000011002092170008220400000"; mesonCutArray[ 2] = "01621035009000"; //clean cuts
 128                 cutarray[ 3] = "8000011002092170008260400000"; mesonCutArray[ 3] = "01621035009000"; //clean cuts ///New STANDARD CUT
 129         } else if (trainConfig == 6) {
 130                 cutarray[ 0] = "8000012002092170008220000000"; mesonCutArray[ 0] = "01621035009000"; //tighten psi pair and qt in 2D
 131                 cutarray[ 1] = "8000012002092170008260000000"; mesonCutArray[ 1] = "01621035009000"; //tighten psi pair and chi2 in 2D and qt in 2D
 132                 cutarray[ 2] = "8000012002092170008220400000"; mesonCutArray[ 2] = "01621035009000"; //clean cuts
 133                 cutarray[ 3] = "8000012002092170008260400000"; mesonCutArray[ 3] = "01621035009000"; //clean cuts ///New STANDARD CUT
 134         } else if (trainConfig == 7) {
 135                 cutarray[ 0] = "8000011001092170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8  // minR 2.8
 136                 cutarray[ 1] = "8000011009092170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8  //minR 7.5
 137                 cutarray[ 2] = "8000011002792170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0. GeV/c
 138                 cutarray[ 3] = "8000011002192170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0.1 GeV/c
 139         } else if (trainConfig == 8) {
 140                 cutarray[ 0] = "8000012001092170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8  // minR 2.8
 141                 cutarray[ 1] = "8000012009092170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8  //minR 7.5
 142                 cutarray[ 2] = "8000012002792170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0. GeV/c
 143                 cutarray[ 3] = "8000012002192170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0.1 GeV/c
 144         } else if (trainConfig == 9) {
 145                 cutarray[ 0] = "8000011002082170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.35
 146                 cutarray[ 1] = "8000011002062170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.7
 147                 cutarray[ 2] = "8000011002093170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -4,5
 148                 cutarray[ 3] = "8000011002096170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -2.5,4
 149         } else if (trainConfig == 10) {
 150                 cutarray[ 0] = "8000012002082170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.35
 151                 cutarray[ 1] = "8000012002062170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.7
 152                 cutarray[ 2] = "8000012002093170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -4,5
 153                 cutarray[ 3] = "8000012002096170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -2.5,4
 154         } else if (trainConfig == 11) {
 155                 cutarray[ 0] = "8000011002092270008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 1,-10
 156                 cutarray[ 1] = "8000011002092570008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 2,-10
 157                 cutarray[ 2] = "8000011002092160008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi minMom 0.25
 158                 cutarray[ 3] = "8000011002092150008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi min Mom 0.3
 159         } else if (trainConfig == 12) {
 160                 cutarray[ 0] = "8000012002092270008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 1,-10
 161                 cutarray[ 1] = "8000012002092570008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 2,-10
 162                 cutarray[ 2] = "8000012002092160008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi minMom 0.25
 163                 cutarray[ 3] = "8000012002092150008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi min Mom 0.3
 164         } else if (trainConfig == 13) {
 165                 cutarray[ 0] = "8000011002092172008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi maxMom 4GeV
 166                 cutarray[ 1] = "8000011002092162008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi minMom 0.25GeV maxMom 4GeV
 167                 cutarray[ 2] = "8000011002092260008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV
 168                 cutarray[ 3] = "8000011002092262008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV maxMom 4GeV
 169         } else if (trainConfig == 14) {
 170                 cutarray[ 0] = "8000012002092172008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi maxMom 4GeV
 171                 cutarray[ 1] = "8000012002092162008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi minMom 0.25GeV maxMom 4GeV
 172                 cutarray[ 2] = "8000012002092260008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV
 173                 cutarray[ 3] = "8000012002092262008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV maxMom 4GeV
 174         } else if (trainConfig == 15) {
 175                 cutarray[ 0] = "8000011002092170003260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.05 1D
 176                 cutarray[ 1] = "8000011002092170009260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.03 2D
 177                 cutarray[ 2] = "8000011002092170002260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.07 1D
 178                 cutarray[ 3] = "8000011002092170008220400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 and psi Pair 1D
 179         } else if (trainConfig == 16) {
 180                 cutarray[ 0] = "8000012002092170003260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.05 1D
 181                 cutarray[ 1] = "8000012002092170009260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.03 2D
 182                 cutarray[ 2] = "8000012002092170002260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.07 1D
 183                 cutarray[ 3] = "8000012002092170008220400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 and psi Pair 1D
 184         } else if (trainConfig == 17) {
 185                 cutarray[ 0] = "8000011002092170008160400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 50  2D
 186                 cutarray[ 1] = "8000011002092170008860400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 20  2D
 187                 cutarray[ 2] = "8000011002092170008250400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.1
 188                 cutarray[ 3] = "8000011002092170008270400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.035
 189         } else if (trainConfig == 18) {
 190                 cutarray[ 0] = "8000012002092170008160400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 50  2D
 191                 cutarray[ 1] = "8000012002092170008860400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 20  2D
 192                 cutarray[ 2] = "8000012002092170008250400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.1
 193                 cutarray[ 3] = "8000012002092170008270400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.035
 194         } else if (trainConfig == 19) {
 195                 cutarray[ 0] = "8000011002092170008260300000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // cos pointing angle 0.75
 196                 cutarray[ 1] = "8000011002092170008260600000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //cos pointing angle 0.9
 197                 cutarray[ 2] = "8000011002092170008260400000"; mesonCutArray[ 2] = "01621065009000";  //new standard eta=0.9 y=0.8 // alpha meson cut 0.8
 198                 cutarray[ 3] = "8000011002092170008260400000"; mesonCutArray[ 3] = "01621034009000";  //new standard eta=0.9 y=0.8 //chi2 meson 500
 199         } else if (trainConfig == 20) {
 200                 cutarray[ 0] = "8000012002092170008260300000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // cos pointing angle 0.75
 201                 cutarray[ 1] = "8000012002092170008260600000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //cos pointing angle 0.9
 202                 cutarray[ 2] = "8000012002092170008260400000"; mesonCutArray[ 2] = "01621065009000";  //new standard eta=0.9 y=0.8 // alpha meson cut 0.8
 203                 cutarray[ 3] = "8000012002092170008260400000"; mesonCutArray[ 3] = "01621034009000";  //new standard eta=0.9 y=0.8 //chi2 meson 500
 204         } else if (trainConfig == 21) {
 205                 cutarray[ 0] = "8000011002092170008260400000"; mesonCutArray[ 0] = "02621035009000";  //new standard eta=0.9 y=0.8 // BG track multiplicity
 206                 cutarray[ 1] = "8000011002092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //no MCP smearing
 207                 cutarray[ 2] = "8000011002093172003290000000"; mesonCutArray[ 2] = "01621035009000";  //old standard eta=0.9 y=0.8
 208                 cutarray[ 3] = "8000011002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
 209         } else if (trainConfig == 22) {
 210                 cutarray[ 0] = "8000012002092170008260400000"; mesonCutArray[ 0] = "02621035009000";  //new standard eta=0.9 y=0.8 // BG track multiplicity
 211                 cutarray[ 1] = "8000012002092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //no MCP smearing
 212                 cutarray[ 2] = "8000012002093172003290000000"; mesonCutArray[ 2] = "01621035009000";  //old standard eta=0.9 y=0.8
 213                 cutarray[ 3] = "8000012002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
 214         } else if (trainConfig == 23){
 215                 cutarray[ 0] = "8020011032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
 216                 cutarray[ 1] = "8020011042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
 217                 cutarray[ 2] = "8020011012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
 218                 cutarray[ 3] = "8020011002092370023220000000"; mesonCutArray[ 3] = "01621035009000";
 219         } else if (trainConfig == 24) {
 220                 cutarray[ 0] = "8020012032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
 221                 cutarray[ 1] = "8020012042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
 222                 cutarray[ 2] = "8020012012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
 223                 cutarray[ 3] = "8020012002092370023220000000"; mesonCutArray[ 3] = "01621035009000";
 224         } else if (trainConfig == 25) {
 225                 cutarray[ 0] = "8020011002092770023220000000"; mesonCutArray[ 0] = "01621035009000";
 226                 cutarray[ 1] = "8020011002092551023220000000"; mesonCutArray[ 1] = "01621035009000";
 227                 cutarray[ 2] = "8020011002092170003220000000"; mesonCutArray[ 2] = "01621035009000"; //just tighten Psi pair
 228                 cutarray[ 3] = "8020011002092170003260000000"; mesonCutArray[ 3] = "01621035009000"; //tighten Psi pair and chi2 in 2D
 229         } else if (trainConfig == 26) {
 230                 cutarray[ 0] = "8020012002092770023220000000"; mesonCutArray[ 0] = "01621035009000";
 231                 cutarray[ 1] = "8020012002092551023220000000"; mesonCutArray[ 1] = "01621035009000";
 232                 cutarray[ 2] = "8020012002092170003220000000"; mesonCutArray[ 2] = "01621035009000"; //just tighten Psi pair
 233                 cutarray[ 3] = "8020012002092170003260000000"; mesonCutArray[ 3] = "01621035009000"; //tighten Psi pair and chi2 in 2D
 234         } else if (trainConfig == 27) {
 235                 cutarray[ 0] = "8020011002092170008220000000"; mesonCutArray[ 0] = "01621035009000"; //tighten psi pair and qt in 2D
 236                 cutarray[ 1] = "8020011002092170008260000000"; mesonCutArray[ 1] = "01621035009000"; //tighten psi pair and chi2 in 2D and qt in 2D
 237                 cutarray[ 2] = "8020011002092170008220400000"; mesonCutArray[ 2] = "01621035009000"; //clean cuts
 238                 cutarray[ 3] = "8020011002092170008260400000"; mesonCutArray[ 3] = "01621035009000"; //clean cuts
 239         } else if (trainConfig == 28) {
 240                 cutarray[ 0] = "8020012002092170008220000000"; mesonCutArray[ 0] = "01621035009000"; //tighten psi pair and qt in 2D
 241                 cutarray[ 1] = "8020012002092170008260000000"; mesonCutArray[ 1] = "01621035009000"; //tighten psi pair and chi2 in 2D and qt in 2D
 242                 cutarray[ 2] = "8020012002092170008220400000"; mesonCutArray[ 2] = "01621035009000"; //clean cuts
 243                 cutarray[ 3] = "8020012002092170008260400000"; mesonCutArray[ 3] = "01621035009000"; //clean cuts
 244         } else if (trainConfig == 29) {
 245                 cutarray[ 0] = "8020011001092170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8  // minR 2.8
 246                 cutarray[ 1] = "8020011009092170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8  //minR 7.5
 247                 cutarray[ 2] = "8020011002792170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0. GeV/c
 248                 cutarray[ 3] = "8020011002192170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0.1 GeV/c
 249         } else if (trainConfig == 30) {
 250                 cutarray[ 0] = "8020012001092170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8  // minR 2.8
 251                 cutarray[ 1] = "8020012009092170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8  //minR 7.5
 252                 cutarray[ 2] = "8020012002792170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0. GeV/c
 253                 cutarray[ 3] = "8020012002192170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0.1 GeV/c
 254         } else if (trainConfig == 31) {
 255                 cutarray[ 0] = "8020011002082170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.35
 256                 cutarray[ 1] = "8020011002062170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.7
 257                 cutarray[ 2] = "8020011002093170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -4,5
 258                 cutarray[ 3] = "8020011002096170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -2.5,4
 259         } else if (trainConfig == 32) {
 260                 cutarray[ 0] = "8020012002082170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.35
 261                 cutarray[ 1] = "8020012002062170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.7
 262                 cutarray[ 2] = "8020012002093170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -4,5
 263                 cutarray[ 3] = "8020012002096170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -2.5,4
 264         } else if (trainConfig == 33) {
 265                 cutarray[ 0] = "8020011002092270008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 1,-10
 266                 cutarray[ 1] = "8020011002092570008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 2,-10
 267                 cutarray[ 2] = "8020011002092160008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi minMom 0.25
 268                 cutarray[ 3] = "8020011002092150008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi min Mom 0.3
 269         } else if (trainConfig == 34) {
 270                 cutarray[ 0] = "8020012002092270008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 1,-10
 271                 cutarray[ 1] = "8020012002092570008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 2,-10
 272                 cutarray[ 2] = "8020012002092160008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi minMom 0.25
 273                 cutarray[ 3] = "8020012002092150008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi min Mom 0.3
 274         } else if (trainConfig == 35) {
 275                 cutarray[ 0] = "8020011002092172008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi maxMom 4GeV
 276                 cutarray[ 1] = "8020011002092162008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi minMom 0.25GeV maxMom 4GeV
 277                 cutarray[ 2] = "8020011002092260008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV
 278                 cutarray[ 3] = "8020011002092262008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV maxMom 4GeV
 279         } else if (trainConfig == 36) {
 280                 cutarray[ 0] = "8020012002092172008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi maxMom 4GeV
 281                 cutarray[ 1] = "8020012002092162008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi minMom 0.25GeV maxMom 4GeV
 282                 cutarray[ 2] = "8020012002092260008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV
 283                 cutarray[ 3] = "8020012002092262008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV maxMom 4GeV
 284         } else if (trainConfig == 37) {
 285                 cutarray[ 0] = "8020011002092170003260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.05 1D
 286                 cutarray[ 1] = "8020011002092170009260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.03 2D
 287                 cutarray[ 2] = "8020011002092170002260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.07 1D
 288                 cutarray[ 3] = "8020011002092170008220400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 and psi Pair 1D
 289         } else if (trainConfig == 38) {
 290                 cutarray[ 0] = "8020012002092170003260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.05 1D
 291                 cutarray[ 1] = "8020012002092170009260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.03 2D
 292                 cutarray[ 2] = "8020012002092170002260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.07 1D
 293                 cutarray[ 3] = "8020012002092170008220400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 and psi Pair 1D
 294         } else if (trainConfig == 39) {
 295                 cutarray[ 0] = "8020011002092170008160400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 50  2D
 296                 cutarray[ 1] = "8020011002092170008860400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 20  2D
 297                 cutarray[ 2] = "8020011002092170008250400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.1
 298                 cutarray[ 3] = "8020011002092170008270400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.035
 299         } else if (trainConfig == 40) {
 300                 cutarray[ 0] = "8020012002092170008160400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 50  2D
 301                 cutarray[ 1] = "8020012002092170008860400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 20  2D
 302                 cutarray[ 2] = "8020012002092170008250400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.1
 303                 cutarray[ 3] = "8020012002092170008270400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.035
 304         } else if (trainConfig == 41) {
 305                 cutarray[ 0] = "8020011002092170008260300000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // cos pointing angle 0.75
 306                 cutarray[ 1] = "8020011002092170008260600000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //cos pointing angle 0.9
 307                 cutarray[ 2] = "8020011002092170008260400000"; mesonCutArray[ 2] = "01621065009000";  //new standard eta=0.9 y=0.8 // alpha meson cut 0.8
 308                 cutarray[ 3] = "8020011002092170008260400000"; mesonCutArray[ 3] = "01621034009000";  //new standard eta=0.9 y=0.8 //chi2 meson 500
 309         } else if (trainConfig == 42) {
 310                 cutarray[ 0] = "8020012002092170008260300000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // cos pointing angle 0.75
 311                 cutarray[ 1] = "8020012002092170008260600000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //cos pointing angle 0.9
 312                 cutarray[ 2] = "8020012002092170008260400000"; mesonCutArray[ 2] = "01621065009000";  //new standard eta=0.9 y=0.8 // alpha meson cut 0.8
 313                 cutarray[ 3] = "8020012002092170008260400000"; mesonCutArray[ 3] = "01621034009000";  //new standard eta=0.9 y=0.8 //chi2 meson 500
 314         } else if (trainConfig == 43) {
 315                 cutarray[ 0] = "8020011002092170008260400000"; mesonCutArray[ 0] = "02621035009000";  //new standard eta=0.9 y=0.8 // BG track multiplicity
 316                 cutarray[ 1] = "8020011002092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //no MCP smearing
 317                 cutarray[ 2] = "8020011002093172003290000000"; mesonCutArray[ 2] = "01621035009000";  //old standard eta=0.9 y=0.8
 318                 cutarray[ 3] = "8020011002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
 319         } else if (trainConfig == 44) {
 320                 cutarray[ 0] = "8020012002092170008260400000"; mesonCutArray[ 0] = "02621035009000";  //new standard eta=0.9 y=0.8 // BG track multiplicity
 321                 cutarray[ 1] = "8020012002092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //no MCP smearing
 322                 cutarray[ 2] = "8020012002093172003290000000"; mesonCutArray[ 2] = "01621035009000";  //old standard eta=0.9 y=0.8
 323                 cutarray[ 3] = "8020012002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
 324         } else if (trainConfig == 45){
 325                 cutarray[ 0] = "8240011032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
 326                 cutarray[ 1] = "8240011042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
 327                 cutarray[ 2] = "8240011012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
 328                 cutarray[ 3] = "8240011002092370023220000000"; mesonCutArray[ 3] = "01621035009000";
 329         } else if (trainConfig == 46) {
 330                 cutarray[ 0] = "8240012032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
 331                 cutarray[ 1] = "8240012042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
 332                 cutarray[ 2] = "8240012012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
 333                 cutarray[ 3] = "8240012002092370023220000000"; mesonCutArray[ 3] = "01621035009000";
 334         } else if (trainConfig == 47) {
 335                 cutarray[ 0] = "8240011002092770023220000000"; mesonCutArray[ 0] = "01621035009000";
 336                 cutarray[ 1] = "8240011002092551023220000000"; mesonCutArray[ 1] = "01621035009000";
 337                 cutarray[ 2] = "8240011002092170003220000000"; mesonCutArray[ 2] = "01621035009000"; //just tighten Psi pair
 338                 cutarray[ 3] = "8240011002092170003260000000"; mesonCutArray[ 3] = "01621035009000"; //tighten Psi pair and chi2 in 2D
 339         } else if (trainConfig == 48) {
 340                 cutarray[ 0] = "8240012002092770023220000000"; mesonCutArray[ 0] = "01621035009000";
 341                 cutarray[ 1] = "8240012002092551023220000000"; mesonCutArray[ 1] = "01621035009000";
 342                 cutarray[ 2] = "8240012002092170003220000000"; mesonCutArray[ 2] = "01621035009000"; //just tighten Psi pair
 343                 cutarray[ 3] = "8240012002092170003260000000"; mesonCutArray[ 3] = "01621035009000"; //tighten Psi pair and chi2 in 2D
 344         } else if (trainConfig == 49) {
 345                 cutarray[ 0] = "8240011002092170008220000000"; mesonCutArray[ 0] = "01621035009000"; //tighten psi pair and qt in 2D
 346                 cutarray[ 1] = "8240011002092170008260000000"; mesonCutArray[ 1] = "01621035009000"; //tighten psi pair and chi2 in 2D and qt in 2D
 347                 cutarray[ 2] = "8240011002092170008220400000"; mesonCutArray[ 2] = "01621035009000"; //clean cuts
 348                 cutarray[ 3] = "8240011002092170008260400000"; mesonCutArray[ 3] = "01621035009000"; //clean cuts
 349         } else if (trainConfig == 50) {
 350                 cutarray[ 0] = "8240012002092170008220000000"; mesonCutArray[ 0] = "01621035009000"; //tighten psi pair and qt in 2D
 351                 cutarray[ 1] = "8240012002092170008260000000"; mesonCutArray[ 1] = "01621035009000"; //tighten psi pair and chi2 in 2D and qt in 2D
 352                 cutarray[ 2] = "8240012002092170008220400000"; mesonCutArray[ 2] = "01621035009000"; //clean cuts
 353                 cutarray[ 3] = "8240012002092170008260400000"; mesonCutArray[ 3] = "01621035009000"; //clean cuts
 354         } else if (trainConfig == 51) {
 355                 cutarray[ 0] = "8240011001092170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8  // minR 2.8
 356                 cutarray[ 1] = "8240011009092170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8  //minR 7.5
 357                 cutarray[ 2] = "8240011002792170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0. GeV/c
 358                 cutarray[ 3] = "8240011002192170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0.1 GeV/c
 359         } else if (trainConfig == 52) {
 360                 cutarray[ 0] = "8240012001092170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8  // minR 2.8
 361                 cutarray[ 1] = "8240012009092170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8  //minR 7.5
 362                 cutarray[ 2] = "8240012002792170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0. GeV/c
 363                 cutarray[ 3] = "8240012002192170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0.1 GeV/c
 364         } else if (trainConfig == 53) {
 365                 cutarray[ 0] = "8240011002082170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.35
 366                 cutarray[ 1] = "8240011002062170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.7
 367                 cutarray[ 2] = "8240011002093170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -4,5
 368                 cutarray[ 3] = "8240011002096170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -2.5,4
 369         } else if (trainConfig == 54) {
 370                 cutarray[ 0] = "8240012002082170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.35
 371                 cutarray[ 1] = "8240012002062170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.7
 372                 cutarray[ 2] = "8240012002093170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -4,5
 373                 cutarray[ 3] = "8240012002096170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -2.5,4
 374         } else if (trainConfig == 55) {
 375                 cutarray[ 0] = "8240011002092270008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 1,-10
 376                 cutarray[ 1] = "8240011002092570008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 2,-10
 377                 cutarray[ 2] = "8240011002092160008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi minMom 0.25
 378                 cutarray[ 3] = "8240011002092150008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi min Mom 0.3
 379         } else if (trainConfig == 56) {
 380                 cutarray[ 0] = "8240012002092270008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 1,-10
 381                 cutarray[ 1] = "8240012002092570008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 2,-10
 382                 cutarray[ 2] = "8240012002092160008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi minMom 0.25
 383                 cutarray[ 3] = "8240012002092150008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi min Mom 0.3
 384         } else if (trainConfig == 57) {
 385                 cutarray[ 0] = "8240011002092172008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi maxMom 4GeV
 386                 cutarray[ 1] = "8240011002092162008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi minMom 0.25GeV maxMom 4GeV
 387                 cutarray[ 2] = "8240011002092260008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV
 388                 cutarray[ 3] = "8240011002092262008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV maxMom 4GeV
 389         } else if (trainConfig == 58) {
 390                 cutarray[ 0] = "8240012002092172008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi maxMom 4GeV
 391                 cutarray[ 1] = "8240012002092162008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi minMom 0.25GeV maxMom 4GeV
 392                 cutarray[ 2] = "8240012002092260008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV
 393                 cutarray[ 3] = "8240012002092262008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV maxMom 4GeV
 394         } else if (trainConfig == 59) {
 395                 cutarray[ 0] = "8240011002092170003260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.05 1D
 396                 cutarray[ 1] = "8240011002092170009260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.03 2D
 397                 cutarray[ 2] = "8240011002092170002260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.07 1D
 398                 cutarray[ 3] = "8240011002092170008220400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 and psi Pair 1D
 399         } else if (trainConfig == 60) {
 400                 cutarray[ 0] = "8240012002092170003260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.05 1D
 401                 cutarray[ 1] = "8240012002092170009260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.03 2D
 402                 cutarray[ 2] = "8240012002092170002260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.07 1D
 403                 cutarray[ 3] = "8240012002092170008220400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 and psi Pair 1D
 404         } else if (trainConfig == 61) {
 405                 cutarray[ 0] = "8240011002092170008160400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 50  2D
 406                 cutarray[ 1] = "8240011002092170008860400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 20  2D
 407                 cutarray[ 2] = "8240011002092170008250400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.1
 408                 cutarray[ 3] = "8240011002092170008270400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.035
 409         } else if (trainConfig == 62) {
 410                 cutarray[ 0] = "8240012002092170008160400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 50  2D
 411                 cutarray[ 1] = "8240012002092170008860400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 20  2D
 412                 cutarray[ 2] = "8240012002092170008250400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.1
 413                 cutarray[ 3] = "8240012002092170008270400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.035
 414         } else if (trainConfig == 63) {
 415                 cutarray[ 0] = "8240011002092170008260300000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // cos pointing angle 0.75
 416                 cutarray[ 1] = "8240011002092170008260600000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //cos pointing angle 0.9
 417                 cutarray[ 2] = "8240011002092170008260400000"; mesonCutArray[ 2] = "01621065009000";  //new standard eta=0.9 y=0.8 // alpha meson cut 0.8
 418                 cutarray[ 3] = "8240011002092170008260400000"; mesonCutArray[ 3] = "01621034009000";  //new standard eta=0.9 y=0.8 //chi2 meson 500
 419         } else if (trainConfig == 64) {
 420                 cutarray[ 0] = "8240012002092170008260300000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // cos pointing angle 0.75
 421                 cutarray[ 1] = "8240012002092170008260600000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //cos pointing angle 0.9
 422                 cutarray[ 2] = "8240012002092170008260400000"; mesonCutArray[ 2] = "01621065009000";  //new standard eta=0.9 y=0.8 // alpha meson cut 0.8
 423                 cutarray[ 3] = "8240012002092170008260400000"; mesonCutArray[ 3] = "01621034009000";  //new standard eta=0.9 y=0.8 //chi2 meson 500
 424         } else if (trainConfig == 65) {
 425                 cutarray[ 0] = "8240011002092170008260400000"; mesonCutArray[ 0] = "02621035009000";  //new standard eta=0.9 y=0.8 // BG track multiplicity
 426                 cutarray[ 1] = "8240011002092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //no MCP smearing
 427                 cutarray[ 2] = "8240011002093172003290000000"; mesonCutArray[ 2] = "01621035009000";  //old standard eta=0.9 y=0.8
 428                 cutarray[ 3] = "8240011002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
 429         } else if (trainConfig == 66) {
 430                 cutarray[ 0] = "8240012002092170008260400000"; mesonCutArray[ 0] = "02621035009000";  //new standard eta=0.9 y=0.8 // BG track multiplicity
 431                 cutarray[ 1] = "8240012002092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //no MCP smearing
 432                 cutarray[ 2] = "8240012002093172003290000000"; mesonCutArray[ 2] = "01621035009000";  //old standard eta=0.9 y=0.8
 433                 cutarray[ 3] = "8240012002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
 434         } else if (trainConfig == 67){
 435                 cutarray[ 0] = "8460011032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
 436                 cutarray[ 1] = "8460011042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
 437                 cutarray[ 2] = "8460011012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
 438                 cutarray[ 3] = "8460011002092370023220000000"; mesonCutArray[ 3] = "01621035009000";
 439         } else if (trainConfig == 68) {
 440                 cutarray[ 0] = "8460012032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
 441                 cutarray[ 1] = "8460012042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
 442                 cutarray[ 2] = "8460012012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
 443                 cutarray[ 3] = "8460012002092370023220000000"; mesonCutArray[ 3] = "01621035009000";
 444         } else if (trainConfig == 69) {
 445                 cutarray[ 0] = "8460011002092770023220000000"; mesonCutArray[ 0] = "01621035009000";
 446                 cutarray[ 1] = "8460011002092551023220000000"; mesonCutArray[ 1] = "01621035009000";
 447                 cutarray[ 2] = "8460011002092170003220000000"; mesonCutArray[ 2] = "01621035009000"; //just tighten Psi pair
 448                 cutarray[ 3] = "8460011002092170003260000000"; mesonCutArray[ 3] = "01621035009000"; //tighten Psi pair and chi2 in 2D
 449         } else if (trainConfig == 70) {
 450                 cutarray[ 0] = "8460012002092770023220000000"; mesonCutArray[ 0] = "01621035009000";
 451                 cutarray[ 1] = "8460012002092551023220000000"; mesonCutArray[ 1] = "01621035009000";
 452                 cutarray[ 2] = "8460012002092170003220000000"; mesonCutArray[ 2] = "01621035009000"; //just tighten Psi pair
 453                 cutarray[ 3] = "8460012002092170003260000000"; mesonCutArray[ 3] = "01621035009000"; //tighten Psi pair and chi2 in 2D
 454         } else if (trainConfig == 71) {
 455                 cutarray[ 0] = "8460011002092170008220000000"; mesonCutArray[ 0] = "01621035009000"; //tighten psi pair and qt in 2D
 456                 cutarray[ 1] = "8460011002092170008260000000"; mesonCutArray[ 1] = "01621035009000"; //tighten psi pair and chi2 in 2D and qt in 2D
 457                 cutarray[ 2] = "8460011002092170008220400000"; mesonCutArray[ 2] = "01621035009000"; //clean cuts
 458                 cutarray[ 3] = "8460011002092170008260400000"; mesonCutArray[ 3] = "01621035009000"; //clean cuts
 459         } else if (trainConfig == 72) {
 460                 cutarray[ 0] = "8460012002092170008220000000"; mesonCutArray[ 0] = "01621035009000"; //tighten psi pair and qt in 2D
 461                 cutarray[ 1] = "8460012002092170008260000000"; mesonCutArray[ 1] = "01621035009000"; //tighten psi pair and chi2 in 2D and qt in 2D
 462                 cutarray[ 2] = "8460012002092170008220400000"; mesonCutArray[ 2] = "01621035009000"; //clean cuts
 463                 cutarray[ 3] = "8460012002092170008260400000"; mesonCutArray[ 3] = "01621035009000"; //clean cuts
 464         } else if (trainConfig == 73) {
 465                 cutarray[ 0] = "8460011001092170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8  // minR 2.8
 466                 cutarray[ 1] = "8460011009092170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8  //minR 7.5
 467                 cutarray[ 2] = "8460011002792170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0. GeV/c
 468                 cutarray[ 3] = "8460011002192170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0.1 GeV/c
 469         } else if (trainConfig == 74) {
 470                 cutarray[ 0] = "8460012001092170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8  // minR 2.8
 471                 cutarray[ 1] = "8460012009092170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8  //minR 7.5
 472                 cutarray[ 2] = "8460012002792170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0. GeV/c
 473                 cutarray[ 3] = "8460012002192170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0.1 GeV/c
 474         } else if (trainConfig == 75) {
 475                 cutarray[ 0] = "8460011002082170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.35
 476                 cutarray[ 1] = "8460011002062170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.7
 477                 cutarray[ 2] = "8460011002093170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -4,5
 478                 cutarray[ 3] = "8460011002096170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -2.5,4
 479         } else if (trainConfig == 76) {
 480                 cutarray[ 0] = "8460012002082170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.35
 481                 cutarray[ 1] = "8460012002062170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.7
 482                 cutarray[ 2] = "8460012002093170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -4,5
 483                 cutarray[ 3] = "8460012002096170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -2.5,4
 484         } else if (trainConfig == 77) {
 485                 cutarray[ 0] = "8460011002092270008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 1,-10
 486                 cutarray[ 1] = "8460011002092570008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 2,-10
 487                 cutarray[ 2] = "8460011002092160008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi minMom 0.25
 488                 cutarray[ 3] = "8460011002092150008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi min Mom 0.3
 489         } else if (trainConfig == 78) {
 490                 cutarray[ 0] = "8460012002092270008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 1,-10
 491                 cutarray[ 1] = "8460012002092570008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 2,-10
 492                 cutarray[ 2] = "8460012002092160008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi minMom 0.25
 493                 cutarray[ 3] = "8460012002092150008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi min Mom 0.3
 494         } else if (trainConfig == 79) {
 495                 cutarray[ 0] = "8460011002092172008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi maxMom 4GeV
 496                 cutarray[ 1] = "8460011002092162008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi minMom 0.25GeV maxMom 4GeV
 497                 cutarray[ 2] = "8460011002092260008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV
 498                 cutarray[ 3] = "8460011002092262008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV maxMom 4GeV
 499         } else if (trainConfig == 80) {
 500                 cutarray[ 0] = "8460012002092172008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi maxMom 4GeV
 501                 cutarray[ 1] = "8460012002092162008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi minMom 0.25GeV maxMom 4GeV
 502                 cutarray[ 2] = "8460012002092260008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV
 503                 cutarray[ 3] = "8460012002092262008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV maxMom 4GeV
 504         } else if (trainConfig == 81) {
 505                 cutarray[ 0] = "8460011002092170003260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.05 1D
 506                 cutarray[ 1] = "8460011002092170009260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.03 2D
 507                 cutarray[ 2] = "8460011002092170002260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.07 1D
 508                 cutarray[ 3] = "8460011002092170008220400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 and psi Pair 1D
 509         } else if (trainConfig == 82) {
 510                 cutarray[ 0] = "8460012002092170003260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.05 1D
 511                 cutarray[ 1] = "8460012002092170009260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.03 2D
 512                 cutarray[ 2] = "8460012002092170002260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.07 1D
 513                 cutarray[ 3] = "8460012002092170008220400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 and psi Pair 1D
 514         } else if (trainConfig == 83) {
 515                 cutarray[ 0] = "8460011002092170008160400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 50  2D
 516                 cutarray[ 1] = "8460011002092170008860400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 20  2D
 517                 cutarray[ 2] = "8460011002092170008250400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.1
 518                 cutarray[ 3] = "8460011002092170008270400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.035
 519         } else if (trainConfig == 84) {
 520                 cutarray[ 0] = "8460012002092170008160400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 50  2D
 521                 cutarray[ 1] = "8460012002092170008860400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 20  2D
 522                 cutarray[ 2] = "8460012002092170008250400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.1
 523                 cutarray[ 3] = "8460012002092170008270400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.035
 524         } else if (trainConfig == 85) {
 525                 cutarray[ 0] = "8460011002092170008260300000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // cos pointing angle 0.75
 526                 cutarray[ 1] = "8460011002092170008260600000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //cos pointing angle 0.9
 527                 cutarray[ 2] = "8460011002092170008260400000"; mesonCutArray[ 2] = "01621065009000";  //new standard eta=0.9 y=0.8 // alpha meson cut 0.8
 528                 cutarray[ 3] = "8460011002092170008260400000"; mesonCutArray[ 3] = "01621034009000";  //new standard eta=0.9 y=0.8 //chi2 meson 500
 529         } else if (trainConfig == 86) {
 530                 cutarray[ 0] = "8460012002092170008260300000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // cos pointing angle 0.75
 531                 cutarray[ 1] = "8460012002092170008260600000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //cos pointing angle 0.9
 532                 cutarray[ 2] = "8460012002092170008260400000"; mesonCutArray[ 2] = "01621065009000";  //new standard eta=0.9 y=0.8 // alpha meson cut 0.8
 533                 cutarray[ 3] = "8460012002092170008260400000"; mesonCutArray[ 3] = "01621034009000";  //new standard eta=0.9 y=0.8 //chi2 meson 500
 534         } else if (trainConfig == 87) {
 535                 cutarray[ 0] = "8460011002092170008260400000"; mesonCutArray[ 0] = "02621035009000";  //new standard eta=0.9 y=0.8 // BG track multiplicity
 536                 cutarray[ 1] = "8460011002092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //no MCP smearing
 537                 cutarray[ 2] = "8460011002093172003290000000"; mesonCutArray[ 2] = "01621035009000";  //old standard eta=0.9 y=0.8
 538                 cutarray[ 3] = "8460011002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
 539         } else if (trainConfig == 88) {
 540                 cutarray[ 0] = "8460012002092170008260400000"; mesonCutArray[ 0] = "02621035009000";  //new standard eta=0.9 y=0.8 // BG track multiplicity
 541                 cutarray[ 1] = "8460012002092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //no MCP smearing
 542                 cutarray[ 2] = "8460012002093172003290000000"; mesonCutArray[ 2] = "01621035009000";  //old standard eta=0.9 y=0.8
 543                 cutarray[ 3] = "8460012002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
 544         } else if (trainConfig == 89){
 545                 cutarray[ 0] = "8680011032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
 546                 cutarray[ 1] = "8680011042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
 547                 cutarray[ 2] = "8680011012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
 548                 cutarray[ 3] = "8680011002092370023220000000"; mesonCutArray[ 3] = "01621035009000";
 549         } else if (trainConfig == 90) {
 550                 cutarray[ 0] = "8680012032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
 551                 cutarray[ 1] = "8680012042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
 552                 cutarray[ 2] = "8680012012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
 553                 cutarray[ 3] = "8680012002092370023220000000"; mesonCutArray[ 3] = "01621035009000";
 554         } else if (trainConfig == 91) {
 555                 cutarray[ 0] = "8680011002092770023220000000"; mesonCutArray[ 0] = "01621035009000";
 556                 cutarray[ 1] = "8680011002092551023220000000"; mesonCutArray[ 1] = "01621035009000";
 557                 cutarray[ 2] = "8680011002092170003220000000"; mesonCutArray[ 2] = "01621035009000"; //just tighten Psi pair
 558                 cutarray[ 3] = "8680011002092170003260000000"; mesonCutArray[ 3] = "01621035009000"; //tighten Psi pair and chi2 in 2D
 559         } else if (trainConfig == 92) {
 560                 cutarray[ 0] = "8680012002092770023220000000"; mesonCutArray[ 0] = "01621035009000";
 561                 cutarray[ 1] = "8680012002092551023220000000"; mesonCutArray[ 1] = "01621035009000";
 562                 cutarray[ 2] = "8680012002092170003220000000"; mesonCutArray[ 2] = "01621035009000"; //just tighten Psi pair
 563                 cutarray[ 3] = "8680012002092170003260000000"; mesonCutArray[ 3] = "01621035009000"; //tighten Psi pair and chi2 in 2D
 564         } else if (trainConfig == 93) {
 565                 cutarray[ 0] = "8680011002092170008220000000"; mesonCutArray[ 0] = "01621035009000"; //tighten psi pair and qt in 2D
 566                 cutarray[ 1] = "8680011002092170008260000000"; mesonCutArray[ 1] = "01621035009000"; //tighten psi pair and chi2 in 2D and qt in 2D
 567                 cutarray[ 2] = "8680011002092170008220400000"; mesonCutArray[ 2] = "01621035009000"; //clean cuts
 568                 cutarray[ 3] = "8680011002092170008260400000"; mesonCutArray[ 3] = "01621035009000"; //clean cuts
 569         } else if (trainConfig == 94) {
 570                 cutarray[ 0] = "8680012002092170008220000000"; mesonCutArray[ 0] = "01621035009000"; //tighten psi pair and qt in 2D
 571                 cutarray[ 1] = "8680012002092170008260000000"; mesonCutArray[ 1] = "01621035009000"; //tighten psi pair and chi2 in 2D and qt in 2D
 572                 cutarray[ 2] = "8680012002092170008220400000"; mesonCutArray[ 2] = "01621035009000"; //clean cuts
 573                 cutarray[ 3] = "8680012002092170008260400000"; mesonCutArray[ 3] = "01621035009000"; //clean cuts
 574         } else if (trainConfig == 95) {
 575                 cutarray[ 0] = "8680011001092170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8  // minR 2.8
 576                 cutarray[ 1] = "8680011009092170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8  //minR 7.5
 577                 cutarray[ 2] = "8680011002792170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0. GeV/c
 578                 cutarray[ 3] = "8680011002192170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0.1 GeV/c
 579         } else if (trainConfig == 96) {
 580                 cutarray[ 0] = "8680012001092170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8  // minR 2.8
 581                 cutarray[ 1] = "8680012009092170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8  //minR 7.5
 582                 cutarray[ 2] = "8680012002792170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0. GeV/c
 583                 cutarray[ 3] = "8680012002192170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0.1 GeV/c
 584         } else if (trainConfig == 97) {
 585                 cutarray[ 0] = "8680011002082170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.35
 586                 cutarray[ 1] = "8680011002062170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.7
 587                 cutarray[ 2] = "8680011002093170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -4,5
 588                 cutarray[ 3] = "8680011002096170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -2.5,4
 589         } else if (trainConfig == 98) {
 590                 cutarray[ 0] = "8680012002082170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.35
 591                 cutarray[ 1] = "8680012002062170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.7
 592                 cutarray[ 2] = "8680012002093170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -4,5
 593                 cutarray[ 3] = "8680012002096170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -2.5,4
 594         } else if (trainConfig == 99) {
 595                 cutarray[ 0] = "8680011002092270008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 1,-10
 596                 cutarray[ 1] = "8680011002092570008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 2,-10
 597                 cutarray[ 2] = "8680011002092160008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi minMom 0.25
 598                 cutarray[ 3] = "8680011002092150008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi min Mom 0.3
 599         } else if (trainConfig == 100) {
 600                 cutarray[ 0] = "8680012002092270008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 1,-10
 601                 cutarray[ 1] = "8680012002092570008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 2,-10
 602                 cutarray[ 2] = "8680012002092160008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi minMom 0.25
 603                 cutarray[ 3] = "8680012002092150008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi min Mom 0.3
 604         } else if (trainConfig == 101) {
 605                 cutarray[ 0] = "8680011002092172008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi maxMom 4GeV
 606                 cutarray[ 1] = "8680011002092162008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi minMom 0.25GeV maxMom 4GeV
 607                 cutarray[ 2] = "8680011002092260008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV
 608                 cutarray[ 3] = "8680011002092262008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV maxMom 4GeV
 609         } else if (trainConfig == 102) {
 610                 cutarray[ 0] = "8680012002092172008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi maxMom 4GeV
 611                 cutarray[ 1] = "8680012002092162008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi minMom 0.25GeV maxMom 4GeV
 612                 cutarray[ 2] = "8680012002092260008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV
 613                 cutarray[ 3] = "8680012002092262008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV maxMom 4GeV
 614         } else if (trainConfig == 103) {
 615                 cutarray[ 0] = "8680011002092170003260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.05 1D
 616                 cutarray[ 1] = "8680011002092170009260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.03 2D
 617                 cutarray[ 2] = "8680011002092170002260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.07 1D
 618                 cutarray[ 3] = "8680011002092170008220400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 and psi Pair 1D
 619         } else if (trainConfig == 104) {
 620                 cutarray[ 0] = "8680012002092170003260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.05 1D
 621                 cutarray[ 1] = "8680012002092170009260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.03 2D
 622                 cutarray[ 2] = "8680012002092170002260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.07 1D
 623                 cutarray[ 3] = "8680012002092170008220400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 and psi Pair 1D
 624         } else if (trainConfig == 105) {
 625                 cutarray[ 0] = "8680011002092170008160400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 50  2D
 626                 cutarray[ 1] = "8680011002092170008860400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 20  2D
 627                 cutarray[ 2] = "8680011002092170008250400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.1
 628                 cutarray[ 3] = "8680011002092170008270400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.035
 629         } else if (trainConfig == 106) {
 630                 cutarray[ 0] = "8680012002092170008160400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 50  2D
 631                 cutarray[ 1] = "8680012002092170008860400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 20  2D
 632                 cutarray[ 2] = "8680012002092170008250400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.1
 633                 cutarray[ 3] = "8680012002092170008270400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.035
 634         } else if (trainConfig == 107) {
 635                 cutarray[ 0] = "8680011002092170008260300000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // cos pointing angle 0.75
 636                 cutarray[ 1] = "8680011002092170008260600000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //cos pointing angle 0.9
 637                 cutarray[ 2] = "8680011002092170008260400000"; mesonCutArray[ 2] = "01621065009000";  //new standard eta=0.9 y=0.8 // alpha meson cut 0.8
 638                 cutarray[ 3] = "8680011002092170008260400000"; mesonCutArray[ 3] = "01621034009000";  //new standard eta=0.9 y=0.8 //chi2 meson 500
 639         } else if (trainConfig == 108) {
 640                 cutarray[ 0] = "8680012002092170008260300000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // cos pointing angle 0.75
 641                 cutarray[ 1] = "8680012002092170008260600000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //cos pointing angle 0.9
 642                 cutarray[ 2] = "8680012002092170008260400000"; mesonCutArray[ 2] = "01621065009000";  //new standard eta=0.9 y=0.8 // alpha meson cut 0.8
 643                 cutarray[ 3] = "8680012002092170008260400000"; mesonCutArray[ 3] = "01621034009000";  //new standard eta=0.9 y=0.8 //chi2 meson 500
 644         } else if (trainConfig == 109) {
 645                 cutarray[ 0] = "8680011002092170008260400000"; mesonCutArray[ 0] = "02621035009000";  //new standard eta=0.9 y=0.8 // BG track multiplicity
 646                 cutarray[ 1] = "8680011002092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //no MCP smearing
 647                 cutarray[ 2] = "8680011002093172003290000000"; mesonCutArray[ 2] = "01621035009000";  //old standard eta=0.9 y=0.8
 648                 cutarray[ 3] = "8680011002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
 649         } else if (trainConfig == 110) {
 650                 cutarray[ 0] = "8680012002092170008260400000"; mesonCutArray[ 0] = "02621035009000";  //new standard eta=0.9 y=0.8 // BG track multiplicity
 651                 cutarray[ 1] = "8680012002092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //no MCP smearing
 652                 cutarray[ 2] = "8680012002093172003290000000"; mesonCutArray[ 2] = "01621035009000";  //old standard eta=0.9 y=0.8
 653                 cutarray[ 3] = "8680012002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
 654         } else if (trainConfig == 111){
 655                 cutarray[ 0] = "8600011032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
 656                 cutarray[ 1] = "8600011042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
 657                 cutarray[ 2] = "8600011012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
 658                 cutarray[ 3] = "8600011002092370023220000000"; mesonCutArray[ 3] = "01621035009000";
 659         } else if (trainConfig == 112) {
 660                 cutarray[ 0] = "8600012032092170008260400000"; mesonCutArray[ 0] = "01623035009000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
 661                 cutarray[ 1] = "8600012042092170008260400000"; mesonCutArray[ 1] = "01622035009000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
 662                 cutarray[ 2] = "8600012012092170008260400000"; mesonCutArray[ 2] = "01624035009000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
 663                 cutarray[ 3] = "8600012002092370023220000000"; mesonCutArray[ 3] = "01621035009000";
 664         } else if (trainConfig == 113) {
 665                 cutarray[ 0] = "8600011002092770023220000000"; mesonCutArray[ 0] = "01621035009000";
 666                 cutarray[ 1] = "8600011002092551023220000000"; mesonCutArray[ 1] = "01621035009000";
 667                 cutarray[ 2] = "8600011002092170003220000000"; mesonCutArray[ 2] = "01621035009000"; //just tighten Psi pair
 668                 cutarray[ 3] = "8600011002092170003260000000"; mesonCutArray[ 3] = "01621035009000"; //tighten Psi pair and chi2 in 2D
 669         } else if (trainConfig == 114) {
 670                 cutarray[ 0] = "8600012002092770023220000000"; mesonCutArray[ 0] = "01621035009000";
 671                 cutarray[ 1] = "8600012002092551023220000000"; mesonCutArray[ 1] = "01621035009000";
 672                 cutarray[ 2] = "8600012002092170003220000000"; mesonCutArray[ 2] = "01621035009000"; //just tighten Psi pair
 673                 cutarray[ 3] = "8600012002092170003260000000"; mesonCutArray[ 3] = "01621035009000"; //tighten Psi pair and chi2 in 2D
 674         } else if (trainConfig == 115) {
 675                 cutarray[ 0] = "8600011002092170008220000000"; mesonCutArray[ 0] = "01621035009000"; //tighten psi pair and qt in 2D
 676                 cutarray[ 1] = "8600011002092170008260000000"; mesonCutArray[ 1] = "01621035009000"; //tighten psi pair and chi2 in 2D and qt in 2D
 677                 cutarray[ 2] = "8600011002092170008220400000"; mesonCutArray[ 2] = "01621035009000"; //clean cuts
 678                 cutarray[ 3] = "8600011002092170008260400000"; mesonCutArray[ 3] = "01621035009000"; //clean cuts
 679         } else if (trainConfig == 116) {
 680                 cutarray[ 0] = "8600012002092170008220000000"; mesonCutArray[ 0] = "01621035009000"; //tighten psi pair and qt in 2D
 681                 cutarray[ 1] = "8600012002092170008260000000"; mesonCutArray[ 1] = "01621035009000"; //tighten psi pair and chi2 in 2D and qt in 2D
 682                 cutarray[ 2] = "8600012002092170008220400000"; mesonCutArray[ 2] = "01621035009000"; //clean cuts
 683                 cutarray[ 3] = "8600012002092170008260400000"; mesonCutArray[ 3] = "01621035009000"; //clean cuts
 684         } else if (trainConfig == 117) {
 685                 cutarray[ 0] = "8600011001092170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8  // minR 2.8
 686                 cutarray[ 1] = "8600011009092170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8  //minR 7.5
 687                 cutarray[ 2] = "8600011002792170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0. GeV/c
 688                 cutarray[ 3] = "8600011002192170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0.1 GeV/c
 689         } else if (trainConfig == 118) {
 690                 cutarray[ 0] = "8600012001092170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8  // minR 2.8
 691                 cutarray[ 1] = "8600012009092170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8  //minR 7.5
 692                 cutarray[ 2] = "8600012002792170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0. GeV/c
 693                 cutarray[ 3] = "8600012002192170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8  //single pT 0.1 GeV/c
 694         } else if (trainConfig == 119) {
 695                 cutarray[ 0] = "8600011002082170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.35
 696                 cutarray[ 1] = "8600011002062170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.7
 697                 cutarray[ 2] = "8600011002093170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -4,5
 698                 cutarray[ 3] = "8600011002096170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -2.5,4
 699         } else if (trainConfig == 120) {
 700                 cutarray[ 0] = "8600012002082170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.35
 701                 cutarray[ 1] = "8600012002062170008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.7
 702                 cutarray[ 2] = "8600012002093170008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -4,5
 703                 cutarray[ 3] = "8600012002096170008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //edEdx -2.5,4
 704         } else if (trainConfig == 121) {
 705                 cutarray[ 0] = "8600011002092270008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 1,-10
 706                 cutarray[ 1] = "8600011002092570008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 2,-10
 707                 cutarray[ 2] = "8600011002092160008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi minMom 0.25
 708                 cutarray[ 3] = "8600011002092150008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi min Mom 0.3
 709         } else if (trainConfig == 122) {
 710                 cutarray[ 0] = "8600012002092270008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 1,-10
 711                 cutarray[ 1] = "8600012002092570008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi dEdx 2,-10
 712                 cutarray[ 2] = "8600012002092160008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi minMom 0.25
 713                 cutarray[ 3] = "8600012002092150008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi min Mom 0.3
 714         } else if (trainConfig == 123) {
 715                 cutarray[ 0] = "8600011002092172008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi maxMom 4GeV
 716                 cutarray[ 1] = "8600011002092162008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi minMom 0.25GeV maxMom 4GeV
 717                 cutarray[ 2] = "8600011002092260008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV
 718                 cutarray[ 3] = "8600011002092262008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV maxMom 4GeV
 719         } else if (trainConfig == 124) {
 720                 cutarray[ 0] = "8600012002092172008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi maxMom 4GeV
 721                 cutarray[ 1] = "8600012002092162008260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 // pi minMom 0.25GeV maxMom 4GeV
 722                 cutarray[ 2] = "8600012002092260008260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV
 723                 cutarray[ 3] = "8600012002092262008260400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV maxMom 4GeV
 724         } else if (trainConfig == 125) {
 725                 cutarray[ 0] = "8600011002092170003260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.05 1D
 726                 cutarray[ 1] = "8600011002092170009260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.03 2D
 727                 cutarray[ 2] = "8600011002092170002260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.07 1D
 728                 cutarray[ 3] = "8600011002092170008220400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 and psi Pair 1D
 729         } else if (trainConfig == 126) {
 730                 cutarray[ 0] = "8600012002092170003260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.05 1D
 731                 cutarray[ 1] = "8600012002092170009260400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.03 2D
 732                 cutarray[ 2] = "8600012002092170002260400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //qT 0.07 1D
 733                 cutarray[ 3] = "8600012002092170008220400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 and psi Pair 1D
 734         } else if (trainConfig == 127) {
 735                 cutarray[ 0] = "8600011002092170008160400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 50  2D
 736                 cutarray[ 1] = "8600011002092170008860400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 20  2D
 737                 cutarray[ 2] = "8600011002092170008250400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.1
 738                 cutarray[ 3] = "8600011002092170008270400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.035
 739         } else if (trainConfig == 128) {
 740                 cutarray[ 0] = "8600012002092170008160400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 50  2D
 741                 cutarray[ 1] = "8600012002092170008860400000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //chi2 20  2D
 742                 cutarray[ 2] = "8600012002092170008250400000"; mesonCutArray[ 2] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.1
 743                 cutarray[ 3] = "8600012002092170008270400000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 //psi pair 0.035
 744         } else if (trainConfig == 129) {
 745                 cutarray[ 0] = "8600011002092170008260300000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // cos pointing angle 0.75
 746                 cutarray[ 1] = "8600011002092170008260600000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //cos pointing angle 0.9
 747                 cutarray[ 2] = "8600011002092170008260400000"; mesonCutArray[ 2] = "01621065009000";  //new standard eta=0.9 y=0.8 // alpha meson cut 0.8
 748                 cutarray[ 3] = "8600011002092170008260400000"; mesonCutArray[ 3] = "01621034009000";  //new standard eta=0.9 y=0.8 //chi2 meson 500
 749         } else if (trainConfig == 130) {
 750                 cutarray[ 0] = "8600012002092170008260300000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // cos pointing angle 0.75
 751                 cutarray[ 1] = "8600012002092170008260600000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //cos pointing angle 0.9
 752                 cutarray[ 2] = "8600012002092170008260400000"; mesonCutArray[ 2] = "01621065009000";  //new standard eta=0.9 y=0.8 // alpha meson cut 0.8
 753                 cutarray[ 3] = "8600012002092170008260400000"; mesonCutArray[ 3] = "01621034009000";  //new standard eta=0.9 y=0.8 //chi2 meson 500
 754         } else if (trainConfig == 131) {
 755                 cutarray[ 0] = "8600011002092170008260400000"; mesonCutArray[ 0] = "02621035009000";  //new standard eta=0.9 y=0.8 // BG track multiplicity
 756                 cutarray[ 1] = "8600011002092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //no MCP smearing
 757                 cutarray[ 2] = "8600011002093172003290000000"; mesonCutArray[ 2] = "01621035009000";  //old standard eta=0.9 y=0.8
 758                 cutarray[ 3] = "8600011002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
 759         } else if (trainConfig == 132) {
 760                 cutarray[ 0] = "8600012002092170008260400000"; mesonCutArray[ 0] = "02621035009000";  //new standard eta=0.9 y=0.8 // BG track multiplicity
 761                 cutarray[ 1] = "8600012002092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //no MCP smearing
 762                 cutarray[ 2] = "8600012002093172003290000000"; mesonCutArray[ 2] = "01621035009000";  //old standard eta=0.9 y=0.8
 763                 cutarray[ 3] = "8600012002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
 764         } else if (trainConfig == 133) {
 765                 cutarray[ 0] = "8000011032092170008260400000"; mesonCutArray[ 0] = "01623035000000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
 766                 cutarray[ 1] = "8000011042092170008260400000"; mesonCutArray[ 1] = "01622035000000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
 767                 cutarray[ 2] = "8000011012092170008260400000"; mesonCutArray[ 2] = "01624035000000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
 768                 cutarray[ 3] = "8000011002092170003260000000"; mesonCutArray[ 3] = "01621035000000"; //tighten Psi pair and chi2 in 2D
 769         } else if (trainConfig == 134) {
 770                 cutarray[ 0] = "8000012032092170008260400000"; mesonCutArray[ 0] = "01623035000000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
 771                 cutarray[ 1] = "8000012042092170008260400000"; mesonCutArray[ 1] = "01622035000000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
 772                 cutarray[ 2] = "8000012012092170008260400000"; mesonCutArray[ 2] = "01624035000000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
 773                 cutarray[ 3] = "8000012002092170003260000000"; mesonCutArray[ 3] = "01621035000000"; //tighten Psi pair and chi2 in 2D
 774         } else if (trainConfig == 135) {
 775                 cutarray[ 0] = "8000011002092170008220000000"; mesonCutArray[ 0] = "01621035000000"; //tighten psi pair and qt in 2D
 776                 cutarray[ 1] = "8000011002092170008260000000"; mesonCutArray[ 1] = "01621035000000"; //tighten psi pair and chi2 in 2D and qt in 2D
 777                 cutarray[ 2] = "8000011002092170008220400000"; mesonCutArray[ 2] = "01621035000000"; //clean cuts
 778                 cutarray[ 3] = "8000011002092170008260400000"; mesonCutArray[ 3] = "01621035000000"; //clean cuts
 779         } else if (trainConfig == 136) {
 780                 cutarray[ 0] = "8000012002092170008220000000"; mesonCutArray[ 0] = "01621035000000"; //tighten psi pair and qt in 2D
 781                 cutarray[ 1] = "8000012002092170008260000000"; mesonCutArray[ 1] = "01621035000000"; //tighten psi pair and chi2 in 2D and qt in 2D
 782                 cutarray[ 2] = "8000012002092170008220400000"; mesonCutArray[ 2] = "01621035000000"; //clean cuts
 783                 cutarray[ 3] = "8000012002092170008260400000"; mesonCutArray[ 3] = "01621035000000"; //clean cuts
 784         } else if (trainConfig == 137) {
 785                 cutarray[ 0] = "8000011001092170008260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8  // minR 2.8
 786                 cutarray[ 1] = "8000011009092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8  //minR 7.5
 787                 cutarray[ 2] = "8000011002792170008260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8  //single pT 0. GeV/c
 788                 cutarray[ 3] = "8000011002192170008260400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8  //single pT 0.1 GeV/c
 789         } else if (trainConfig == 138) {
 790                 cutarray[ 0] = "8000012001092170008260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8  // minR 2.8
 791                 cutarray[ 1] = "8000012009092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8  //minR 7.5
 792                 cutarray[ 2] = "8000012002792170008260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8  //single pT 0. GeV/c
 793                 cutarray[ 3] = "8000012002192170008260400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8  //single pT 0.1 GeV/c
 794         } else if (trainConfig == 139) {
 795                 cutarray[ 0] = "8000011002082170008260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.35
 796                 cutarray[ 1] = "8000011002062170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.7
 797                 cutarray[ 2] = "8000011002093170008260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //edEdx -4,5
 798                 cutarray[ 3] = "8000011002096170008260400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //edEdx -2.5,4
 799         } else if (trainConfig == 140) {
 800                 cutarray[ 0] = "8000012002082170008260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.35
 801                 cutarray[ 1] = "8000012002062170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.7
 802                 cutarray[ 2] = "8000012002093170008260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //edEdx -4,5
 803                 cutarray[ 3] = "8000012002096170008260400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //edEdx -2.5,4
 804         } else if (trainConfig == 141) {
 805                 cutarray[ 0] = "8000011002092270008260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 // pi dEdx 1,-10
 806                 cutarray[ 1] = "8000011002092570008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 // pi dEdx 2,-10
 807                 cutarray[ 2] = "8000011002092160008260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi minMom 0.25
 808                 cutarray[ 3] = "8000011002092150008260400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi min Mom 0.3
 809         } else if (trainConfig == 142) {
 810                 cutarray[ 0] = "8000012002092270008260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 // pi dEdx 1,-10
 811                 cutarray[ 1] = "8000012002092570008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 // pi dEdx 2,-10
 812                 cutarray[ 2] = "8000012002092160008260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi minMom 0.25
 813                 cutarray[ 3] = "8000012002092150008260400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi min Mom 0.3
 814         } else if (trainConfig == 143) {
 815                 cutarray[ 0] = "8000011002092172008260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi maxMom 4GeV
 816                 cutarray[ 1] = "8000011002092162008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 // pi minMom 0.25GeV maxMom 4GeV
 817                 cutarray[ 2] = "8000011002092260008260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV
 818                 cutarray[ 3] = "8000011002092262008260400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV maxMom 4GeV
 819         } else if (trainConfig == 144) {
 820                 cutarray[ 0] = "8000012002092172008260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi maxMom 4GeV
 821                 cutarray[ 1] = "8000012002092162008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 // pi minMom 0.25GeV maxMom 4GeV
 822                 cutarray[ 2] = "8000012002092260008260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV
 823                 cutarray[ 3] = "8000012002092262008260400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV maxMom 4GeV
 824         } else if (trainConfig == 145) {
 825                 cutarray[ 0] = "8000011002092170003260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 //qT 0.05 1D
 826                 cutarray[ 1] = "8000011002092170009260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //qT 0.03 2D
 827                 cutarray[ 2] = "8000011002092170002260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //qT 0.07 1D
 828                 cutarray[ 3] = "8000011002092170008220400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //chi2 and psi Pair 1D
 829         } else if (trainConfig == 146) {
 830                 cutarray[ 0] = "8000012002092170003260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 //qT 0.05 1D
 831                 cutarray[ 1] = "8000012002092170009260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //qT 0.03 2D
 832                 cutarray[ 2] = "8000012002092170002260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //qT 0.07 1D
 833                 cutarray[ 3] = "8000012002092170008220400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //chi2 and psi Pair 1D
 834         } else if (trainConfig == 147) {
 835                 cutarray[ 0] = "8000011002092170008160400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 //chi2 50  2D
 836                 cutarray[ 1] = "8000011002092170008860400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //chi2 20  2D
 837                 cutarray[ 2] = "8000011002092170008250400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //psi pair 0.1
 838                 cutarray[ 3] = "8000011002092170008270400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //psi pair 0.035
 839         } else if (trainConfig == 148) {
 840                 cutarray[ 0] = "8000012002092170008160400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 //chi2 50  2D
 841                 cutarray[ 1] = "8000012002092170008860400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //chi2 20  2D
 842                 cutarray[ 2] = "8000012002092170008250400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //psi pair 0.1
 843                 cutarray[ 3] = "8000012002092170008270400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //psi pair 0.035
 844         } else if (trainConfig == 149) {
 845                 cutarray[ 0] = "8000011002092170008260300000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 // cos pointing angle 0.75
 846                 cutarray[ 1] = "8000011002092170008260600000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //cos pointing angle 0.9
 847                 cutarray[ 2] = "8000011002092170008260400000"; mesonCutArray[ 2] = "01621065000000";  //new standard eta=0.9 y=0.8 // alpha meson cut 0.8
 848                 cutarray[ 3] = "8000011002092170008260400000"; mesonCutArray[ 3] = "01621034000000";  //new standard eta=0.9 y=0.8 //chi2 meson 500
 849         } else if (trainConfig == 150) {
 850                 cutarray[ 0] = "8000012002092170008260300000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 // cos pointing angle 0.75
 851                 cutarray[ 1] = "8000012002092170008260600000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //cos pointing angle 0.9
 852                 cutarray[ 2] = "8000012002092170008260400000"; mesonCutArray[ 2] = "01621065000000";  //new standard eta=0.9 y=0.8 // alpha meson cut 0.8
 853                 cutarray[ 3] = "8000012002092170008260400000"; mesonCutArray[ 3] = "01621034000000";  //new standard eta=0.9 y=0.8 //chi2 meson 500
 854         } else if (trainConfig == 151) {
 855                 cutarray[ 0] = "8000011002092170008260400000"; mesonCutArray[ 0] = "02621035000000";  //new standard eta=0.9 y=0.8 // BG track multiplicity
 856                 cutarray[ 1] = "8000011002092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //no MCP smearing
 857                 cutarray[ 2] = "8000011002093172003290000000"; mesonCutArray[ 2] = "01621035000000";  //old standard eta=0.9 y=0.8
 858                 cutarray[ 3] = "8000011002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
 859         } else if (trainConfig == 152) {
 860                 cutarray[ 0] = "8000012002092170008260400000"; mesonCutArray[ 0] = "02621035000000";  //new standard eta=0.9 y=0.8 // BG track multiplicity
 861                 cutarray[ 1] = "8000012002092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //no MCP smearing
 862                 cutarray[ 2] = "8000012002093172003290000000"; mesonCutArray[ 2] = "01621035000000";  //old standard eta=0.9 y=0.8
 863                 cutarray[ 3] = "8000012002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
 864         } else if (trainConfig == 153) {
 865                 cutarray[ 0] = "8020011032092170008260400000"; mesonCutArray[ 0] = "01623035000000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
 866                 cutarray[ 1] = "8020011042092170008260400000"; mesonCutArray[ 1] = "01622035000000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
 867                 cutarray[ 2] = "8020011012092170008260400000"; mesonCutArray[ 2] = "01624035000000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
 868                 cutarray[ 3] = "8020011002092170003260000000"; mesonCutArray[ 3] = "01621035000000"; //tighten Psi pair and chi2 in 2D
 869         } else if (trainConfig == 154) {
 870                 cutarray[ 0] = "8020012032092170008260400000"; mesonCutArray[ 0] = "01623035000000"; //New STANDARD CUT |eta| < 0.65, |y| < 0.6
 871                 cutarray[ 1] = "8020012042092170008260400000"; mesonCutArray[ 1] = "01622035000000"; //New STANDARD CUT |eta| < 0.75, |y| < 0.7
 872                 cutarray[ 2] = "8020012012092170008260400000"; mesonCutArray[ 2] = "01624035000000"; //New STANDARD CUT |eta| < 0.6, |y| < 0.5
 873                 cutarray[ 3] = "8020012002092170003260000000"; mesonCutArray[ 3] = "01621035000000"; //tighten Psi pair and chi2 in 2D
 874         } else if (trainConfig == 155) {
 875                 cutarray[ 0] = "8020011002092170008220000000"; mesonCutArray[ 0] = "01621035000000"; //tighten psi pair and qt in 2D
 876                 cutarray[ 1] = "8020011002092170008260000000"; mesonCutArray[ 1] = "01621035000000"; //tighten psi pair and chi2 in 2D and qt in 2D
 877                 cutarray[ 2] = "8020011002092170008220400000"; mesonCutArray[ 2] = "01621035000000"; //clean cuts
 878                 cutarray[ 3] = "8020011002092170008260400000"; mesonCutArray[ 3] = "01621035000000"; //clean cuts
 879         } else if (trainConfig == 156) {
 880                 cutarray[ 0] = "8020012002092170008220000000"; mesonCutArray[ 0] = "01621035000000"; //tighten psi pair and qt in 2D
 881                 cutarray[ 1] = "8020012002092170008260000000"; mesonCutArray[ 1] = "01621035000000"; //tighten psi pair and chi2 in 2D and qt in 2D
 882                 cutarray[ 2] = "8020012002092170008220400000"; mesonCutArray[ 2] = "01621035000000"; //clean cuts
 883                 cutarray[ 3] = "8020012002092170008260400000"; mesonCutArray[ 3] = "01621035000000"; //clean cuts
 884         } else if (trainConfig == 157) {
 885                 cutarray[ 0] = "8020011001092170008260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8  // minR 2.8
 886                 cutarray[ 1] = "8020011009092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8  //minR 7.5
 887                 cutarray[ 2] = "8020011002792170008260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8  //single pT 0. GeV/c
 888                 cutarray[ 3] = "8020011002192170008260400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8  //single pT 0.1 GeV/c
 889         } else if (trainConfig == 158) {
 890                 cutarray[ 0] = "8020012001092170008260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8  // minR 2.8
 891                 cutarray[ 1] = "8020012009092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8  //minR 7.5
 892                 cutarray[ 2] = "8020012002792170008260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8  //single pT 0. GeV/c
 893                 cutarray[ 3] = "8020012002192170008260400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8  //single pT 0.1 GeV/c
 894         } else if (trainConfig == 159) {
 895                 cutarray[ 0] = "8020011002082170008260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.35
 896                 cutarray[ 1] = "8020011002062170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.7
 897                 cutarray[ 2] = "8020011002093170008260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //edEdx -4,5
 898                 cutarray[ 3] = "8020011002096170008260400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //edEdx -2.5,4
 899         } else if (trainConfig == 160) {
 900                 cutarray[ 0] = "8020012002082170008260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.35
 901                 cutarray[ 1] = "8020012002062170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //TPC Cluster 0.7
 902                 cutarray[ 2] = "8020012002093170008260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //edEdx -4,5
 903                 cutarray[ 3] = "8020012002096170008260400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //edEdx -2.5,4
 904         } else if (trainConfig == 161) {
 905                 cutarray[ 0] = "8020011002092270008260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 // pi dEdx 1,-10
 906                 cutarray[ 1] = "8020011002092570008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 // pi dEdx 2,-10
 907                 cutarray[ 2] = "8020011002092160008260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi minMom 0.25
 908                 cutarray[ 3] = "8020011002092150008260400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi min Mom 0.3
 909         } else if (trainConfig == 162) {
 910                 cutarray[ 0] = "8020012002092270008260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 // pi dEdx 1,-10
 911                 cutarray[ 1] = "8020012002092570008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 // pi dEdx 2,-10
 912                 cutarray[ 2] = "8020012002092160008260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi minMom 0.25
 913                 cutarray[ 3] = "8020012002092150008260400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi min Mom 0.3
 914         } else if (trainConfig == 163) {
 915                 cutarray[ 0] = "8020011002092172008260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi maxMom 4GeV
 916                 cutarray[ 1] = "8020011002092162008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 // pi minMom 0.25GeV maxMom 4GeV
 917                 cutarray[ 2] = "8020011002092260008260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV
 918                 cutarray[ 3] = "8020011002092262008260400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV maxMom 4GeV
 919         } else if (trainConfig == 164) {
 920                 cutarray[ 0] = "8020012002092172008260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi maxMom 4GeV
 921                 cutarray[ 1] = "8020012002092162008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 // pi minMom 0.25GeV maxMom 4GeV
 922                 cutarray[ 2] = "8020012002092260008260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV
 923                 cutarray[ 3] = "8020012002092262008260400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //pi dEdx 1,-10  minMom 0.25GeV maxMom 4GeV
 924         } else if (trainConfig == 165) {
 925                 cutarray[ 0] = "8020011002092170003260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 //qT 0.05 1D
 926                 cutarray[ 1] = "8020011002092170009260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //qT 0.03 2D
 927                 cutarray[ 2] = "8020011002092170002260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //qT 0.07 1D
 928                 cutarray[ 3] = "8020011002092170008220400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //chi2 and psi Pair 1D
 929         } else if (trainConfig == 166) {
 930                 cutarray[ 0] = "8020012002092170003260400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 //qT 0.05 1D
 931                 cutarray[ 1] = "8020012002092170009260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //qT 0.03 2D
 932                 cutarray[ 2] = "8020012002092170002260400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //qT 0.07 1D
 933                 cutarray[ 3] = "8020012002092170008220400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //chi2 and psi Pair 1D
 934         } else if (trainConfig == 167) {
 935                 cutarray[ 0] = "8020011002092170008160400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 //chi2 50  2D
 936                 cutarray[ 1] = "8020011002092170008860400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //chi2 20  2D
 937                 cutarray[ 2] = "8020011002092170008250400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //psi pair 0.1
 938                 cutarray[ 3] = "8020011002092170008270400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //psi pair 0.035
 939         } else if (trainConfig == 168) {
 940                 cutarray[ 0] = "8020012002092170008160400000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 //chi2 50  2D
 941                 cutarray[ 1] = "8020012002092170008860400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //chi2 20  2D
 942                 cutarray[ 2] = "8020012002092170008250400000"; mesonCutArray[ 2] = "01621035000000";  //new standard eta=0.9 y=0.8 //psi pair 0.1
 943                 cutarray[ 3] = "8020012002092170008270400000"; mesonCutArray[ 3] = "01621035000000";  //new standard eta=0.9 y=0.8 //psi pair 0.035
 944         } else if (trainConfig == 169) {
 945                 cutarray[ 0] = "8020011002092170008260300000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 // cos pointing angle 0.75
 946                 cutarray[ 1] = "8020011002092170008260600000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //cos pointing angle 0.9
 947                 cutarray[ 2] = "8020011002092170008260400000"; mesonCutArray[ 2] = "01621065000000";  //new standard eta=0.9 y=0.8 // alpha meson cut 0.8
 948                 cutarray[ 3] = "8020011002092170008260400000"; mesonCutArray[ 3] = "01621034000000";  //new standard eta=0.9 y=0.8 //chi2 meson 500
 949         } else if (trainConfig == 170) {
 950                 cutarray[ 0] = "8020012002092170008260300000"; mesonCutArray[ 0] = "01621035000000";  //new standard eta=0.9 y=0.8 // cos pointing angle 0.75
 951                 cutarray[ 1] = "8020012002092170008260600000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //cos pointing angle 0.9
 952                 cutarray[ 2] = "8020012002092170008260400000"; mesonCutArray[ 2] = "01621065000000";  //new standard eta=0.9 y=0.8 // alpha meson cut 0.8
 953                 cutarray[ 3] = "8020012002092170008260400000"; mesonCutArray[ 3] = "01621034000000";  //new standard eta=0.9 y=0.8 //chi2 meson 500
 954         } else if (trainConfig == 171) {
 955                 cutarray[ 0] = "8020011002092170008260400000"; mesonCutArray[ 0] = "02621035000000";  //new standard eta=0.9 y=0.8 // BG track multiplicity
 956                 cutarray[ 1] = "8020011002092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //no MCP smearing
 957                 cutarray[ 2] = "8020011002093172003290000000"; mesonCutArray[ 2] = "01621035000000";  //old standard eta=0.9 y=0.8
 958                 cutarray[ 3] = "8020011002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
 959         } else if (trainConfig == 172) {
 960                 cutarray[ 0] = "8020012002092170008260400000"; mesonCutArray[ 0] = "02621035000000";  //new standard eta=0.9 y=0.8 // BG track multiplicity
 961                 cutarray[ 1] = "8020012002092170008260400000"; mesonCutArray[ 1] = "01621035000000";  //new standard eta=0.9 y=0.8 //no MCP smearing
 962                 cutarray[ 2] = "8020012002093172003290000000"; mesonCutArray[ 2] = "01621035000000";  //old standard eta=0.9 y=0.8
 963                 cutarray[ 3] = "8020012002092170008260400000"; mesonCutArray[ 3] = "01621035008000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
 964         } else if (trainConfig == 173) {
 965                 cutarray[ 0] = "8000011002092170008260400000"; mesonCutArray[ 0] = "01621035009000"; //clean cuts ///New STANDARD CUT
 966                 cutarray[ 1] = "8000011002092170008260420000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //no MCP smearing
 967                 cutarray[ 2] = "8000011002092170008260430000"; mesonCutArray[ 2] = "01621035009000";  //old standard eta=0.9 y=0.8
 968                 cutarray[ 3] = "8000011002092170008260440000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
 969         } else if (trainConfig == 174) {
 970                 cutarray[ 0] = "8000012002092170008260400000"; mesonCutArray[ 0] = "01621035009000";  //new standard eta=0.9 y=0.8 // BG track multiplicity
 971                 cutarray[ 1] = "8000012002092170008260420000"; mesonCutArray[ 1] = "01621035009000";  //new standard eta=0.9 y=0.8 //no MCP smearing
 972                 cutarray[ 2] = "8000012002092170008260430000"; mesonCutArray[ 2] = "01621035009000";  //old standard eta=0.9 y=0.8
 973                 cutarray[ 3] = "8000012002092170008260440000"; mesonCutArray[ 3] = "01621035009000";  //new standard eta=0.9 y=0.8 // fPSigSmearingCte=0.014;
 974
 975         } else {
 976                 Error(Form("GammaConvV1_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
 977                 return;
 978         }
 979
 980         TList *ConvCutList = new TList();
 981         TList *MesonCutList = new TList();
 982
 983         TList *HeaderList = new TList();
 984         if (doWeightingPart==1) {
 985                 TObjString *Header1 = new TObjString("pi0_1");
 986                 HeaderList->Add(Header1);
 987         }
 988         if (doWeightingPart==2){
 989                 TObjString *Header3 = new TObjString("eta_2");
 990                 HeaderList->Add(Header3);
 991         }
 992         if (doWeightingPart==3) {
 993                 TObjString *Header1 = new TObjString("pi0_1");
 994                 HeaderList->Add(Header1);
 995                 TObjString *Header3 = new TObjString("eta_2");
 996                 HeaderList->Add(Header3);
 997         }
 998
 999
1000         Bool_t doWeighting = kFALSE;
1001         if (doWeightingPart == 1 || doWeightingPart == 2 || doWeightingPart == 3) doWeighting = kTRUE;
1002
1003         ConvCutList->SetOwner(kTRUE);
1004         AliConversionCuts **analysisCuts = new AliConversionCuts*[numberOfCuts];
1005         MesonCutList->SetOwner(kTRUE);
1006         AliConversionMesonCuts **analysisMesonCuts = new AliConversionMesonCuts*[numberOfCuts];
1007
1008
1009         for(Int_t i = 0; i<numberOfCuts; i++){
1010
1011                 analysisCuts[i] = new AliConversionCuts();
1012                 if ( trainConfig == 1 || trainConfig == 3 || trainConfig == 5 || trainConfig == 7 || trainConfig == 9 || trainConfig == 11 || trainConfig == 13 || trainConfig == 15|| trainConfig == 17|| trainConfig == 19 || trainConfig == 21 || trainConfig == 133 || trainConfig == 135 || trainConfig == 137 || trainConfig == 139 || trainConfig == 141 || trainConfig == 143 || trainConfig == 145 || trainConfig == 147 || trainConfig == 149 || trainConfig == 151){
1013                         if (doWeighting){
1014                                 if (generatorName.CompareTo("DPMJET")==0){
1015                                         analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_DPMJET_LHC13b2_efix_pPb_5023GeV_MBV0A", "Eta_DPMJET_LHC13b2_efix_pPb_5023GeV_MBV0A", "","Pi0_Fit_Data_pPb_5023GeV_MBV0A","Eta_Fit_Data_pPb_5023GeV_MBV0A");
1016                                 } else if (generatorName.CompareTo("HIJING")==0){
1017                                         analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13e7_pPb_5023GeV_MBV0A", "Eta_Hijing_LHC13e7_pPb_5023GeV_MBV0A", "","Pi0_Fit_Data_pPb_5023GeV_MBV0A","Eta_Fit_Data_pPb_5023GeV_MBV0A");
1018                                 }
1019                         }
1020                 }
1021                 if ( trainConfig == 2 || trainConfig == 4 || trainConfig == 6 || trainConfig == 8 || trainConfig == 10 || trainConfig == 12 || trainConfig == 14 || trainConfig == 16|| trainConfig == 18|| trainConfig == 20|| trainConfig == 22 || trainConfig == 134 || trainConfig == 136 || trainConfig == 138 || trainConfig == 140 || trainConfig == 142 || trainConfig == 144 || trainConfig == 146 || trainConfig == 148 || trainConfig == 150 || trainConfig == 152){
1022                         if (doWeighting){
1023                                 analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13e7_addSig_pPb_5023GeV_MBV0A", "Eta_Hijing_LHC13e7_addSig_pPb_5023GeV_MBV0A", "","Pi0_Fit_Data_pPb_5023GeV_MBV0A","Eta_Fit_Data_pPb_5023GeV_MBV0A");
1024                         }
1025
1026                 }
1027                 if ( trainConfig == 23 || trainConfig == 25 || trainConfig == 27 || trainConfig == 29 || trainConfig == 31 || trainConfig == 33 || trainConfig == 35 || trainConfig == 37|| trainConfig == 39|| trainConfig == 41 || trainConfig == 43 || trainConfig == 153 || trainConfig == 155 || trainConfig == 157 || trainConfig == 159 || trainConfig == 161 || trainConfig == 163 || trainConfig == 165 || trainConfig == 167 || trainConfig == 169 || trainConfig == 171){
1028                         if (doWeighting){
1029                                 if (generatorName.CompareTo("DPMJET")==0){
1030                                         analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_DPMJET_LHC13b2_efix_pPb_5023GeV_0020V0A", "Eta_DPMJET_LHC13b2_efix_pPb_5023GeV_0020V0A", "","Pi0_Fit_Data_pPb_5023GeV_0020V0A","Eta_Fit_Data_pPb_5023GeV_0020V0A");
1031                                 } else if (generatorName.CompareTo("HIJING")==0){
1032                                         analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13e7_pPb_5023GeV_0020V0A", "Eta_Hijing_LHC13e7_pPb_5023GeV_0020V0A", "","Pi0_Fit_Data_pPb_5023GeV_0020V0A","Eta_Fit_Data_pPb_5023GeV_0020V0A");
1033                                 }
1034                         }
1035                 }
1036                 if ( trainConfig == 24 || trainConfig == 26 || trainConfig == 28 || trainConfig == 30 || trainConfig == 32 || trainConfig == 34 || trainConfig == 36 || trainConfig == 38|| trainConfig == 40|| trainConfig == 42|| trainConfig == 44 || trainConfig == 154 || trainConfig == 156 || trainConfig == 158 || trainConfig == 160 || trainConfig == 162 || trainConfig == 164 || trainConfig == 166 || trainConfig == 168 || trainConfig == 170 || trainConfig == 172){
1037                         if (doWeighting){
1038                                 analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13e7_addSig_pPb_5023GeV_0020V0A", "Eta_Hijing_LHC13e7_addSig_pPb_5023GeV_0020V0A", "","Pi0_Fit_Data_pPb_5023GeV_0020V0A","Eta_Fit_Data_pPb_5023GeV_0020V0A");
1039                         }
1040                 }
1041                 if ( trainConfig == 45 || trainConfig == 47 || trainConfig == 49 || trainConfig == 51 || trainConfig == 53 || trainConfig == 55 || trainConfig == 57 || trainConfig == 59|| trainConfig == 61|| trainConfig == 63 || trainConfig == 65 ){
1042                         if (doWeighting){
1043                                 if (generatorName.CompareTo("DPMJET")==0){
1044                                         analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_DPMJET_LHC13b2_efix_pPb_5023GeV_2040V0A", "Eta_DPMJET_LHC13b2_efix_pPb_5023GeV_2040V0A", "","Pi0_Fit_Data_pPb_5023GeV_2040V0A","Eta_Fit_Data_pPb_5023GeV_2040V0A");
1045                                 } else if (generatorName.CompareTo("HIJING")==0){
1046                                         analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13e7_pPb_5023GeV_2040V0A", "Eta_Hijing_LHC13e7_pPb_5023GeV_2040V0A", "","Pi0_Fit_Data_pPb_5023GeV_2040V0A","Eta_Fit_Data_pPb_5023GeV_2040V0A");
1047                                 }
1048                         }
1049                 }
1050                 if ( trainConfig == 46 || trainConfig == 48 || trainConfig == 50 || trainConfig == 52 || trainConfig == 54 || trainConfig == 56 || trainConfig == 58 || trainConfig == 60|| trainConfig == 62|| trainConfig == 64|| trainConfig == 66){
1051                         if (doWeighting){
1052                                 analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13e7_addSig_pPb_5023GeV_2040V0A", "Eta_Hijing_LHC13e7_addSig_pPb_5023GeV_2040V0A", "","Pi0_Fit_Data_pPb_5023GeV_2040V0A","Eta_Fit_Data_pPb_5023GeV_2040V0A");
1053                         }
1054                 }
1055                 if ( trainConfig == 67 || trainConfig == 69 || trainConfig == 71 || trainConfig == 73 || trainConfig == 75 || trainConfig == 77 || trainConfig == 79 || trainConfig == 81 || trainConfig == 83 || trainConfig == 85 || trainConfig == 87 ){
1056                         if (doWeighting){
1057                                 if (generatorName.CompareTo("DPMJET")==0){
1058                                         analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_DPMJET_LHC13b2_efix_pPb_5023GeV_4060V0A", "Eta_DPMJET_LHC13b2_efix_pPb_5023GeV_4060V0A", "","Pi0_Fit_Data_pPb_5023GeV_4060V0A","Eta_Fit_Data_pPb_5023GeV_4060V0A");
1059                                 } else if (generatorName.CompareTo("HIJING")==0){
1060                                         analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13e7_pPb_5023GeV_4060V0A", "Eta_Hijing_LHC13e7_pPb_5023GeV_4060V0A", "","Pi0_Fit_Data_pPb_5023GeV_4060V0A","Eta_Fit_Data_pPb_5023GeV_4060V0A");
1061                                 }
1062                         }
1063                 }
1064                 if ( trainConfig == 68 || trainConfig == 70 || trainConfig == 72 || trainConfig == 74 || trainConfig == 76 || trainConfig == 78 || trainConfig == 80 || trainConfig == 82 || trainConfig == 84 || trainConfig == 86 || trainConfig == 88){
1065                         if (doWeighting){
1066                                 analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13e7_addSig_pPb_5023GeV_4060V0A", "Eta_Hijing_LHC13e7_addSig_pPb_5023GeV_4060V0A", "","Pi0_Fit_Data_pPb_5023GeV_4060V0A","Eta_Fit_Data_pPb_5023GeV_4060V0A");
1067                         }
1068                 }
1069                 if ( trainConfig == 89 || trainConfig == 91 || trainConfig == 93 || trainConfig == 95 || trainConfig == 97 || trainConfig == 99 || trainConfig == 101 || trainConfig == 103 || trainConfig == 105 || trainConfig == 107 || trainConfig == 109 ){
1070                         if (doWeighting){
1071                                 if (generatorName.CompareTo("DPMJET")==0){
1072                                         analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_DPMJET_LHC13b2_efix_pPb_5023GeV_6080V0A", "Eta_DPMJET_LHC13b2_efix_pPb_5023GeV_6080V0A", "","Pi0_Fit_Data_pPb_5023GeV_6080V0A","Eta_Fit_Data_pPb_5023GeV_6080V0A");
1073                                 } else if (generatorName.CompareTo("HIJING")==0){
1074                                         analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13e7_pPb_5023GeV_6080V0A", "Eta_Hijing_LHC13e7_pPb_5023GeV_6080V0A", "","Pi0_Fit_Data_pPb_5023GeV_6080V0A","Eta_Fit_Data_pPb_5023GeV_6080V0A");
1075                                 }
1076                         }
1077                 }
1078                 if ( trainConfig == 90 || trainConfig == 92 || trainConfig == 94 || trainConfig == 96 || trainConfig == 98 || trainConfig == 100 || trainConfig == 102 || trainConfig == 104 || trainConfig == 106 || trainConfig == 108 || trainConfig == 108 ){
1079                         if (doWeighting){
1080                                 analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13e7_addSig_pPb_5023GeV_6080V0A", "Eta_Hijing_LHC13e7_addSig_pPb_5023GeV_6080V0A", "","Pi0_Fit_Data_pPb_5023GeV_6080V0A","Eta_Fit_Data_pPb_5023GeV_6080V0A");
1081                         }
1082                 }
1083                 if ( trainConfig == 111 || trainConfig == 113 || trainConfig == 115 || trainConfig == 117 || trainConfig == 119 || trainConfig == 121 || trainConfig == 123 || trainConfig == 125 || trainConfig == 127 || trainConfig == 129 || trainConfig == 131 ){
1084                         if (doWeighting){
1085                                 if (generatorName.CompareTo("DPMJET")==0){
1086                                         analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_DPMJET_LHC13b2_efix_pPb_5023GeV_60100V0A", "Eta_DPMJET_LHC13b2_efix_pPb_5023GeV_60100V0A", "","Pi0_Fit_Data_pPb_5023GeV_60100V0A","Eta_Fit_Data_pPb_5023GeV_60100V0A");
1087                                 } else if (generatorName.CompareTo("HIJING")==0){
1088                                         analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13e7_pPb_5023GeV_60100V0A", "Eta_Hijing_LHC13e7_pPb_5023GeV_60100V0A", "","Pi0_Fit_Data_pPb_5023GeV_60100V0A","Eta_Fit_Data_pPb_5023GeV_60100V0A");
1089                                 }
1090                         }
1091                 }
1092                 if ( trainConfig == 112 || trainConfig == 114 || trainConfig == 116 || trainConfig == 118 || trainConfig == 120 || trainConfig == 122 || trainConfig == 124 || trainConfig == 126 || trainConfig == 128 || trainConfig == 130 || trainConfig == 130){
1093                         if (doWeighting){
1094                                 analysisCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13e7_addSig_pPb_5023GeV_60100V0A", "Eta_Hijing_LHC13e7_addSig_pPb_5023GeV_60100V0A", "","Pi0_Fit_Data_pPb_5023GeV_60100V0A","Eta_Fit_Data_pPb_5023GeV_60100V0A");
1095                         }
1096                 }
1097                 analysisCuts[i]->InitializeCutsFromCutString(cutarray[i].Data());
1098                 if (doEtaShiftIndCuts) {
1099                         analysisCuts[i]->DoEtaShift(doEtaShiftIndCuts);
1100                         analysisCuts[i]->SetEtaShift(stringShift);
1101                 }
1102
1103                 ConvCutList->Add(analysisCuts[i]);
1104
1105                 analysisCuts[i]->SetFillCutHistograms("",kFALSE);
1106                 analysisMesonCuts[i] = new AliConversionMesonCuts();
1107                 analysisMesonCuts[i]->InitializeCutsFromCutString(mesonCutArray[i].Data());
1108                 MesonCutList->Add(analysisMesonCuts[i]);
1109                 analysisMesonCuts[i]->SetFillCutHistograms("");
1110                 analysisCuts[i]->SetAcceptedHeader(HeaderList);
1111         }
1112
1113         task->SetConversionCutList(numberOfCuts,ConvCutList);
1114         task->SetMesonCutList(numberOfCuts,MesonCutList);
1115         task->SetMoveParticleAccordingToVertex(kTRUE);
1116         task->SetDoMesonAnalysis(kTRUE);
1117         task->SetDoMesonQA(enableQAMesonTask); //Attention new switch for Pi0 QA
1118         task->SetDoPhotonQA(enableQAPhotonTask);  //Attention new switch small for Photon QA
1119
1120         //connect containers
1121         AliAnalysisDataContainer *coutput =
1122                 mgr->CreateContainer(Form("GammaConvV1_%i",trainConfig), TList::Class(),
1123                                                         AliAnalysisManager::kOutputContainer,Form("GammaConvV1_%i.root",trainConfig));
1124
1125         mgr->AddTask(task);
1126         mgr->ConnectInput(task,0,cinput);
1127         mgr->ConnectOutput(task,1,coutput);
1128
1129         return;
1130
1131 }