From f6c249db2d880ae650438fe36397993ba7a3d50d Mon Sep 17 00:00:00 2001 From: fbock Date: Wed, 9 Apr 2014 15:58:48 +0200 Subject: [PATCH] - changed boundaries for MC part for PbPb --- PWGGA/GammaConv/AliConversionCuts.cxx | 50 +++--- .../macros/AddTask_GammaConvV1_PbPb.C | 161 ++++++++++++++++-- 2 files changed, 169 insertions(+), 42 deletions(-) diff --git a/PWGGA/GammaConv/AliConversionCuts.cxx b/PWGGA/GammaConv/AliConversionCuts.cxx index 8284d2a8c48..c51bfbd59e9 100644 --- a/PWGGA/GammaConv/AliConversionCuts.cxx +++ b/PWGGA/GammaConv/AliConversionCuts.cxx @@ -3244,40 +3244,40 @@ Bool_t AliConversionCuts::IsCentralitySelected(AliVEvent *event, AliVEvent *fMCE Int_t PrimaryTracks10[10][2] = { {9999,9999}, // 0 - {1210,2067}, // 10 - { 817,1450}, // 20 - { 536, 921}, // 30 - { 337, 572}, // 40 - { 197, 332}, // 50 - { 106, 173}, // 60 - { 51, 81}, // 70 - { 21, 34}, // 80 + {1210,928}, // 10 + { 817,683}, // 20 + { 536, 437}, // 30 + { 337, 277}, // 40 + { 197, 162}, // 50 + { 106, 100}, // 60 + { 51, 50}, // 70 + { 21, 10}, // 80 { 0, 0} // 90 }; Int_t PrimaryTracks5a[10][2] = { {9999,9999}, // 0 - {1485,2562}, // 5 - {1210,2067}, // 10 - { 995,1760}, // 15 - { 817,1450}, // 20 + {1485,1168}, // 5 + {1210,928}, // 10 + { 995,800}, // 15 + { 817,683}, // 20 { 666,1160}, // 25 - { 536, 921}, // 30 - { 428, 731}, // 35 - { 337, 572}, // 40 - { 260, 436} // 45 + { 536, 437}, // 30 + { 428, 350}, // 35 + { 337, 277}, // 40 + { 260, 200} // 45 }; Int_t PrimaryTracks5b[10][2] = { - { 260, 436}, // 45 - { 197, 327}, // 50 - { 147, 239}, // 55 - { 106, 173}, // 60 - { 75, 120}, // 65 - { 51, 81}, // 70 - { 34, 53}, // 75 - { 21, 34}, // 80 - { 13, 19}, // 85 + { 260, 200}, // 45 + { 197, 162}, // 50 + { 147, 125}, // 55 + { 106, 100}, // 60 + { 75, 75}, // 65 + { 51, 50}, // 70 + { 34, 30}, // 75 + { 21, 10}, // 80 + { 13, 5}, // 85 { 0, 0} // 90 }; diff --git a/PWGGA/GammaConv/macros/AddTask_GammaConvV1_PbPb.C b/PWGGA/GammaConv/macros/AddTask_GammaConvV1_PbPb.C index 08e0d52b0d1..96ea5eda0c9 100644 --- a/PWGGA/GammaConv/macros/AddTask_GammaConvV1_PbPb.C +++ b/PWGGA/GammaConv/macros/AddTask_GammaConvV1_PbPb.C @@ -3,7 +3,7 @@ void AddTask_GammaConvV1_PbPb( Int_t trainConfig = 1, //change different set o Int_t enableQAMesonTask = 0, //enable QA in AliAnalysisTaskGammaConvV1 Int_t enableQAPhotonTask = 0, // enable additional QA task TString fileNameInputForWeighting = "MCSpectraInput.root", // path to file for weigting input - Bool_t doWeighting = kFALSE, //enable Weighting + Int_t doWeightingInt = 0, //enable Weighting TString cutnumberAODBranch = "1000000060084000001500000", TString periodName = "LHC13d2" ) { @@ -26,6 +26,9 @@ void AddTask_GammaConvV1_PbPb( Int_t trainConfig = 1, //change different set o gSystem->Load("libSTEERBase.so"); gSystem->Load("libTENDER.so"); gSystem->Load("libTENDERSupplies.so"); + + Bool_t doWeighting = kFALSE; + if (doWeightingInt > 0) doWeighting = kTRUE; // ================== GetAnalysisManager =============================== AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager(); @@ -486,29 +489,153 @@ void AddTask_GammaConvV1_PbPb( Int_t trainConfig = 1, //change different set o cutarray[ 1] = "7670001002092970028250400000"; mesonCutArray[ 1] = "01525065009000"; // 5-10% cutarray[ 2] = "7780001002092970028250400000"; mesonCutArray[ 2] = "01525065009000"; // 0-10% cutarray[ 3] = "7890001002092970028250400000"; mesonCutArray[ 3] = "01525065009000"; // 10-20% - cutarray[ 4] = "7090001002092970028250400000"; mesonCutArray[ 4] = "01525065009000"; // 0-20% + cutarray[ 4] = "7090001002092970028250400000"; mesonCutArray[ 4] = "01525065009000"; // 0-20% + } else if ( trainConfig == 66){ // variation eta 0.65 + cutarray[ 0] = "6010001032092970028250400000"; mesonCutArray[0]= "01523065009000"; // 0-5% + cutarray[ 1] = "6120001032092970028250400000"; mesonCutArray[1]= "01523065009000"; // 5-10% + cutarray[ 2] = "5010001032092970028250400000"; mesonCutArray[2]= "01523065009000"; // 0-10% + cutarray[ 3] = "5240001032092970028250400000"; mesonCutArray[3]= "01523065009000"; // 20-40% + cutarray[ 4] = "5250001032092970028250400000"; mesonCutArray[4]= "01523065009000"; // 20-50% + } else if ( trainConfig == 67){ // variation eta 0.75 + cutarray[ 0] = "6010001072092970028250400000"; mesonCutArray[0]= "01522065009000"; // 0-5% + cutarray[ 1] = "6120001072092970028250400000"; mesonCutArray[1]= "01522065009000"; // 5-10% + cutarray[ 2] = "5010001072092970028250400000"; mesonCutArray[2]= "01522065009000"; // 0-10% + cutarray[ 3] = "5240001072092970028250400000"; mesonCutArray[3]= "01522065009000"; // 20-40% + cutarray[ 4] = "5250001072092970028250400000"; mesonCutArray[4]= "01522065009000"; // 20-50% + } else if ( trainConfig == 68){ // single pt 0.075 + cutarray[ 0] = "6010001002492970028250400000"; mesonCutArray[0]= "01525065009000"; // 0-5% + cutarray[ 1] = "6120001002492970028250400000"; mesonCutArray[1]= "01525065009000"; // 5-10% + cutarray[ 2] = "5010001002492970028250400000"; mesonCutArray[2]= "01525065009000"; // 0-10% + cutarray[ 3] = "5240001002492970028250400000"; mesonCutArray[3]= "01525065009000"; // 20-40% + cutarray[ 4] = "5250001002492970028250400000"; mesonCutArray[4]= "01525065009000"; // 20-50% + } else if ( trainConfig == 69){ // single pt 0.1 + cutarray[ 0] = "6010001002192970028250400000"; mesonCutArray[0]= "01525065009000"; // 0-5% + cutarray[ 1] = "6120001002192970028250400000"; mesonCutArray[1]= "01525065009000"; // 5-10% + cutarray[ 2] = "5010001002192970028250400000"; mesonCutArray[2]= "01525065009000"; // 0-10% + cutarray[ 3] = "5240001002192970028250400000"; mesonCutArray[3]= "01525065009000"; // 20-40% + cutarray[ 4] = "5250001002192970028250400000"; mesonCutArray[4]= "01525065009000"; //20-50% + } else if ( trainConfig == 70){ // variation TPC cls 0.7 + cutarray[ 0] = "6010001002062970028250400000"; mesonCutArray[0]= "01525065009000"; // 0-5% + cutarray[ 1] = "6120001002062970028250400000"; mesonCutArray[1]= "01525065009000"; // 5-10% + cutarray[ 2] = "5010001002062970028250400000"; mesonCutArray[2]= "01525065009000"; // 0-10% + cutarray[ 3] = "5240001002062970028250400000"; mesonCutArray[3]= "01525065009000"; // 20-40% + cutarray[ 4] = "5250001002062970028250400000"; mesonCutArray[4]= "01525065009000"; // 20-50% + } else if ( trainConfig == 71){ // variation TPC cls 0.35 + cutarray[ 0] = "6010001002082970028250400000"; mesonCutArray[0]= "01525065009000"; // 0-5% + cutarray[ 1] = "6120001002082970028250400000"; mesonCutArray[1]= "01525065009000"; // 5-10% + cutarray[ 2] = "5010001002082970028250400000"; mesonCutArray[2]= "01525065009000"; // 0-10% + cutarray[ 0] = "5240001002082970028250400000"; mesonCutArray[3]= "01525065009000"; // 20-40% + cutarray[ 4] = "5250001002082970028250400000"; mesonCutArray[4]= "01525065009000"; // 20-50% + } else if ( trainConfig == 72){ // variation edEdx -4,5 + cutarray[ 0] = "6010001002093970028250400000"; mesonCutArray[0]= "01525065009000"; // 0-5% + cutarray[ 1] = "6120001002093970028250400000"; mesonCutArray[1]= "01525065009000"; // 5-10% + cutarray[ 2] = "5010001002093970028250400000"; mesonCutArray[2]= "01525065009000"; // 0-10% + cutarray[ 3] = "5240001002093970028250400000"; mesonCutArray[3]= "01525065009000"; // 20-40% + cutarray[ 4] = "5250001002093970028250400000"; mesonCutArray[4]= "01525065009000"; // 20-50% + } else if ( trainConfig == 73){ // variation edEdx -2.5,4 + cutarray[ 0] = "6010001002096970028250400000"; mesonCutArray[0]= "01525065009000"; // 0-5% + cutarray[ 1] = "6120001002096970028250400000"; mesonCutArray[1]= "01525065009000"; // 5-10% + cutarray[ 2] = "5010001002096970028250400000"; mesonCutArray[2]= "01525065009000"; // 0-10% + cutarray[ 3] = "5240001002096970028250400000"; mesonCutArray[3]= "01525065009000"; // 20-40% + cutarray[ 4] = "5250001002096970028250400000"; mesonCutArray[4]= "01525065009000"; // 20-50% + } else if ( trainConfig == 74){ //variation pion p dEdx 0.3-5. + cutarray[ 0] = "6010001002092951028250400000"; mesonCutArray[0]= "01525065009000"; // 0-5% + cutarray[ 1] = "6120001002092951028250400000"; mesonCutArray[1]= "01525065009000"; // 5-10% + cutarray[ 2] = "5010001002092951028250400000"; mesonCutArray[2]= "01525065009000"; // 0-10% + cutarray[ 3] = "5240001002092951028250400000"; mesonCutArray[3]= "01525065009000"; // 20-40% + cutarray[ 4] = "5250001002092951028250400000"; mesonCutArray[4]= "01525065009000"; // 20-50% + } else if ( trainConfig == 75){ // TOF el. PID -3,5 + cutarray[ 0] = "6010001002092970038250400000"; mesonCutArray[0]= "01525065009000"; // 0-5% + cutarray[ 1] = "6120001002092970038250400000"; mesonCutArray[1]= "01525065009000"; // 5-10% + cutarray[ 2] = "5010001002092970038250400000"; mesonCutArray[2]= "01525065009000"; // 0-10% + cutarray[ 3] = "5240001002092970038250400000"; mesonCutArray[3]= "01525065009000"; // 20-40% + cutarray[ 4] = "5250001002092970038250400000"; mesonCutArray[4]= "01525065009000"; // 20-50% + } else if ( trainConfig == 76){ // TOF el. PID -2,3 + cutarray[ 0] = "6010001002092970048250400000"; mesonCutArray[0]= "01525065009000"; // 0-5% + cutarray[ 1] = "6120001002092970048250400000"; mesonCutArray[1]= "01525065009000"; // 5-10% + cutarray[ 2] = "5010001002092970048250400000"; mesonCutArray[2]= "01525065009000"; // 0-10% + cutarray[ 3] = "5240001002092970048250400000"; mesonCutArray[3]= "01525065009000"; // 20-40% + cutarray[ 4] = "5250001002092970048250400000"; mesonCutArray[4]= "01525065009000"; // 20-50% + } else if ( trainConfig == 77){ // qt 0.03 + cutarray[ 0] = "6010001002092970029250400000"; mesonCutArray[0]= "01525065009000"; // 0-5% + cutarray[ 1] = "6120001002092970029250400000"; mesonCutArray[1]= "01525065009000"; // 5-10% + cutarray[ 2] = "5010001002092970029250400000"; mesonCutArray[2]= "01525065009000"; // 0-10% + cutarray[ 3] = "5240001002092970029250400000"; mesonCutArray[3]= "01525065009000"; // 20-40% + cutarray[ 4] = "5250001002092970029250400000"; mesonCutArray[4]= "01525065009000"; // 20-50% + } else if ( trainConfig == 78){ // qt 0.07 no2D + cutarray[ 0] = "6010001002092970022250400000"; mesonCutArray[0]= "01525065009000"; // 0-5% + cutarray[ 1] = "6120001002092970022250400000"; mesonCutArray[1]= "01525065009000"; // 5-10% + cutarray[ 2] = "5010001002092970022250400000"; mesonCutArray[2]= "01525065009000"; // 0-10% + cutarray[ 3] = "5240001002092970022250400000"; mesonCutArray[3]= "01525065009000"; // 20-40% + cutarray[ 4] = "5250001002092970022250400000"; mesonCutArray[4]= "01525065009000"; // 20-50% + } else if ( trainConfig == 79){ // chi2 50. + cutarray[ 0] = "6010001002092970028150400000"; mesonCutArray[0]= "01525065009000"; // 0-5% + cutarray[ 1] = "6120001002092970028150400000"; mesonCutArray[1]= "01525065009000"; // 5-10% + cutarray[ 2] = "5010001002092970028150400000"; mesonCutArray[2]= "01525065009000"; // 0-10% + cutarray[ 3] = "5240001002092970028150400000"; mesonCutArray[3]= "01525065009000"; // 20-40% + cutarray[ 4] = "5250001002092970028150400000"; mesonCutArray[4]= "01525065009000"; // 20-50% + } else if ( trainConfig == 80){ // chi2 20. + cutarray[ 0] = "6010001002092970028850400000"; mesonCutArray[0]= "01525065009000"; // 0-5% + cutarray[ 1] = "6120001002092970028850400000"; mesonCutArray[1]= "01525065009000"; // 5-10% + cutarray[ 2] = "5010001002092970028850400000"; mesonCutArray[2]= "01525065009000"; // 0-10% + cutarray[ 3] = "5240001002092970028850400000"; mesonCutArray[3]= "01525065009000"; // 20-40% + cutarray[ 4] = "5250001002092970028850400000"; mesonCutArray[4]= "01525065009000"; // 20-50% + } else if ( trainConfig == 81){ // psi pair 0.05 + cutarray[ 0] = "6010001002092970028260400000"; mesonCutArray[0]= "01525065009000"; // 0-5% + cutarray[ 1] = "6120001002092970028260400000"; mesonCutArray[1]= "01525065009000"; // 5-10% + cutarray[ 2] = "5010001002092970028260400000"; mesonCutArray[2]= "01525065009000"; // 0-10% + cutarray[ 3] = "5240001002092970028260400000"; mesonCutArray[3]= "01525065009000"; // 20-40% + cutarray[ 4] = "5250001002092970028260400000"; mesonCutArray[4]= "01525065009000"; // 20-50% + } else if ( trainConfig == 82){ // cosPA -1 + cutarray[ 0] = "6010001002092970028250000000"; mesonCutArray[0]= "01525065009000"; // 0-5% + cutarray[ 1] = "6120001002092970028250000000"; mesonCutArray[1]= "01525065009000"; // 5-10% + cutarray[ 2] = "5010001002092970028250000000"; mesonCutArray[2]= "01525065009000"; // 0-10% + cutarray[ 3] = "5240001002092970028250000000"; mesonCutArray[3]= "01525065009000"; // 20-40% + cutarray[ 4] = "5250001002092970028250000000"; mesonCutArray[4]= "01525065009000"; // 20-50% + } else if ( trainConfig == 83){ // variation alpha 0.6&0.8 + cutarray[ 0] = "6010001002092970028250400000"; mesonCutArray[0]= "01525085009000"; // 0-5% + cutarray[ 1] = "6120001002092970028250400000"; mesonCutArray[1]= "01525085009000"; // 5-10% + cutarray[ 2] = "5010001002092970028250400000"; mesonCutArray[2]= "01525085009000"; // 0-10% + cutarray[ 3] = "5240001002092970028250400000"; mesonCutArray[3]= "01525065009000"; // 20-40% + cutarray[ 4] = "5250001002092970028250400000"; mesonCutArray[4]= "01525065009000"; // 20-50% + } else if ( trainConfig == 84){ // variation alpha 0.65&0.75 + cutarray[ 0] = "6010001002092970028250400000"; mesonCutArray[0]= "01525045009000"; // 0-5% + cutarray[ 1] = "6120001002092970028250400000"; mesonCutArray[1]= "01525045009000"; // 5-10% + cutarray[ 2] = "5010001002092970028250400000"; mesonCutArray[2]= "01525045009000"; // 0-10% + cutarray[ 3] = "5240001002092970028250400000"; mesonCutArray[3]= "01525055009000"; // 20-40% + cutarray[ 4] = "5250001002092970028250400000"; mesonCutArray[4]= "01525055009000"; // 20-50% } else { Error(Form("GammaConvV1_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration"); return; } - TList *ConvCutList = new TList(); - TList *MesonCutList = new TList(); + TList *ConvCutList = new TList(); + TList *MesonCutList = new TList(); - TList *HeaderList = new TList(); - if (periodName.CompareTo("LHC13d2")==0){ - TObjString *Header1 = new TObjString("pi0_1"); - HeaderList->Add(Header1); - // TObjString *Header3 = new TObjString("eta_2"); - // HeaderList->Add(Header3); + TList *HeaderList = new TList(); + if (periodName.CompareTo("LHC13d2")==0){ + TObjString *Header1 = new TObjString("pi0_1"); + HeaderList->Add(Header1); + // TObjString *Header3 = new TObjString("eta_2"); + // HeaderList->Add(Header3); - } else if (periodName.CompareTo("LHC12a17x_fix")==0){ - TObjString *Header1 = new TObjString("PARAM"); - HeaderList->Add(Header1); - } else if (periodName.CompareTo("LHC14a1x")==0){ - TObjString *Header1 = new TObjString("BOX"); - HeaderList->Add(Header1); - } + } else if (periodName.CompareTo("LHC12a17x_fix")==0){ + TObjString *Header1 = new TObjString("PARAM"); + HeaderList->Add(Header1); + } else if (periodName.CompareTo("LHC14a1x")==0){ + if (doWeightingInt == 1){ + TObjString *Header1 = new TObjString("pi0_1"); + HeaderList->Add(Header1); + } else if (doWeightingInt == 2){ + TObjString *Header1 = new TObjString("eta_2"); + HeaderList->Add(Header1); + } else { + TObjString *Header1 = new TObjString("pi0_1"); + HeaderList->Add(Header1); + TObjString *Header2 = new TObjString("eta_2"); + HeaderList->Add(Header2); + } + } ConvCutList->SetOwner(kTRUE); AliConversionCuts **analysisCuts = new AliConversionCuts*[numberOfCuts]; -- 2.43.0