]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWGDQ/dielectron/AliDielectronVarManager.cxx
including switch to set on/off iso-track core removal, cleaning and bug fix
[u/mrichter/AliRoot.git] / PWGDQ / dielectron / AliDielectronVarManager.cxx
index 9bb48c0093ff3831be388918785d233fd1f900b5..627eb05e597cf231741460825be9fffc4b3488ce 100644 (file)
@@ -31,6 +31,7 @@ const char* AliDielectronVarManager::fgkParticleNames[AliDielectronVarManager::k
   {"Py",                     "#it{p}_{y}",                                         "(GeV/#it{c})"},
   {"Pz",                     "#it{p}_{z}",                                         "(GeV/#it{c})"},
   {"Pt",                     "#it{p}_{T}",                                         "(GeV/#it{c})"},
+  {"PtSq",                   "#it{p}_{T}^{2}",                                     "((GeV/#it{c})^{2})"},
   {"P",                      "#it{p}",                                             "(GeV/#it{c})"},
   {"Xv",                     "x_{vtx}",                                            "(cm)"},
   {"Yv",                     "y_{vtx}",                                            "(cm)"},
@@ -58,6 +59,8 @@ const char* AliDielectronVarManager::fgkParticleNames[AliDielectronVarManager::k
   {"TPCchi2PerCluster",      "#chi^{2}/^{}N_{cls}^{TPC}",                          ""},
   {"TPCclsDiff",             "N_{d#it{E}/d#it{x} points}^{TPC} - N_{cls}^{TPC}",   ""},
   {"TPCclsSegments",         "N_{segments}^{TPC}",                                 ""},
+  {"TPCclsIRO",              "N_{clsIRO}^{TPC}",                                   ""},
+  {"TPCclsORO",              "N_{clsORO}^{TPC}",                                   ""},
   {"TrackStatus",            "TrackStatus",                                        ""},
   {"FilterBit",              "AOD filter bit",                                     ""},
     
@@ -70,7 +73,7 @@ const char* AliDielectronVarManager::fgkParticleNames[AliDielectronVarManager::k
   {"TRDpidProb2D_Electrons", "P(2dim. PID_{e}^{TRD})",                             ""},
   {"TRDpidProb2D_Pions",     "P(2dim. PID_{#pi}^{TRD})",                           ""},
   {"TRDphi",                 "#phi^{TRD}",                                         ""},
-  {"TRDpidEffLeg",           "#epsilon^{TRD legs}(PID)",                           ""},
+  {"TRDpidEffLeg",           "#varepsilon^{TRD legs}(PID)",                           ""},
   {"TRD_signal",             "TRD d#it{E}/d#it{x} + TR",                           "(a.u.)"},
 
   {"ImpactParXY",            "#it{dca}_{xy}",                                      "(cm)"},
@@ -110,8 +113,9 @@ const char* AliDielectronVarManager::fgkParticleNames[AliDielectronVarManager::k
   {"TOF_signal",             "TOF signal",                                         "(ps)"},
   {"TOF_beta",               "TOF #beta",                                          "(#it{c})"},
   {"TOF_PIDbit",             "TOF PID bit",                                        ""},
-  {"TOF_mismProb",           "TOF mismatch probability",                               ""},
-  
+  {"TOF_mismProb",           "TOF mismatch probability"                            ""},
+
+  {"TPC_nSigma_ElectronsRaw","n#sigma_{e}^{TPC} (raw)",                            ""},
   {"TPC_nSigma_Electrons",   "n#sigma_{e}^{TPC}",                                  ""},
   {"TPC_nSigma_Pions",       "n#sigma_{#pi}^{TPC}",                                ""},
   {"TPC_nSigma_Muons",       "n#sigma_{#mu}^{TPC}",                                ""},
@@ -132,6 +136,8 @@ const char* AliDielectronVarManager::fgkParticleNames[AliDielectronVarManager::k
   {"EMCAL_M20",              "M20 EMCAL showershape param.",                       ""},
   {"EMCAL_Dispersion",       "EMCAL dispersion param.",                            ""},
 
+  {"EffLeg",                 "A#times#varepsilon",                                    ""},
+  {"OneOverEffLeg",          "(A#times#varepsilon)^{-1}",                             ""},
   {"V0Index0",               "V0Index0",                                           ""},
   {"KinkIndex0",             "KinkIndex0",                                         ""},
   //
@@ -152,14 +158,31 @@ const char* AliDielectronVarManager::fgkParticleNames[AliDielectronVarManager::k
   {"ThetaSqCS",              "cos^{2}(#theta_{CS})",                               ""},
   {"PsiPair",                "#Psi^{pair}",                                        "(rad.)"},
   {"PhivPair",               "#Phi_{v}^{pair}",                                    "(rad.)"},
-  {"PairPlanev0rpH2Angle",   "PairPlanev0rpH2Angle",                               "(rad.)"},
-  {"PairPlaneMagAngle",      "PairPlaneMagAngle",                                  "(rad.)"},
-  {"PairPlaneAngle",         "PairPlaneAngle",                                     "(rad.)"},
-  {"kRotPairx",                 "kRotPairx",                                          "(rad.)"},
-  {"kRotPairy",                 "kRotPairy",                                          "(rad.)"},
-  {"kRotPairz",                 "kRotPairz",                                          "(rad.)"},
+  {"PairPlaneAngle1A",       "#Phi_{1A}",                                          "(rad.)"},
+  {"PairPlaneAngle2A",       "#Phi_{2A}",                                          "(rad.)"},
+  {"PairPlaneAngle3A",       "#Phi_{3A}",                                          "(rad.)"},
+  {"PairPlaneAngle4A",       "#Phi_{4A}",                                          "(rad.)"},
+  {"PairPlaneAngle1C",       "#Phi_{1C}",                                          "(rad.)"},
+  {"PairPlaneAngle2C",       "#Phi_{2C}",                                          "(rad.)"},
+  {"PairPlaneAngle3C",       "#Phi_{3C}",                                          "(rad.)"},
+  {"PairPlaneAngle4C",       "#Phi_{4C}",                                          "(rad.)"},
+  {"PairPlaneAngle1AC",      "#Phi_{1AC}",                                         "(rad.)"},
+  {"PairPlaneAngle2AC",      "#Phi_{2AC}",                                         "(rad.)"},
+  {"PairPlaneAngle3AC",      "#Phi_{3AC}",                                         "(rad.)"},
+  {"PairPlaneAngle4AC",      "#Phi_{4AC}",                                         "(rad.)"},
+  {"PairPlaneAngle1Ran",     "#Phi_{1Ran}",                                        "(rad.)"},
+  {"PairPlaneAngle2Ran",     "#Phi_{2Ran}",                                        "(rad.)"},
+  {"PairPlaneAngle3Ran",     "#Phi_{3Ran}",                                        "(rad.)"},
+  {"PairPlaneAngle4Ran",     "#Phi_{4Ran}",                                        "(rad.)"},
+  {"RandomRP",               "#Phi_{RanRP}",                                       "(rad.)"},
+  {"DeltaPhiRandomRP",       "#Delta #Phi_{RanRP}",                                ""},
+
+  {"PairPlaneMagInPro",      "cos(#Phi)"                                           ""},
+
   {"Cos2PhiCS",              "cos(2#phi_{CS})",                                    ""},
   {"CosTilPhiCS",            "cos(#phi_{CS})",                                     ""},
+  {"CosPhiH2",               "cos(2#phi)",                                         ""},
+  {"SinPhiH2",               "sin(2#phi)",                                         ""},
   {"DeltaPhiV0ArpH2",        "#phi^{pair}-#Psi^{V0A}",                             ""},
   {"DeltaPhiV0CrpH2",        "#phi^{pair}-#Psi^{V0C}",                             ""},
   {"DeltaPhiV0ACrpH2",       "#phi^{pair}-#Psi^{V0AC}",                            ""},
@@ -169,9 +192,12 @@ const char* AliDielectronVarManager::fgkParticleNames[AliDielectronVarManager::k
   {"DeltaPhiv0ArpH2",        "#phi^{pair}-#Psi^{V0A}",                             ""},
   {"DeltaPhiv0CrpH2",        "#phi^{pair}-#Psi^{V0C}",                             ""},
   {"DeltaPhiv0ACrpH2",       "#phi^{pair}-#Psi^{V0AC}",                            ""},
+  {"DeltaPhiTPCrpH2",        "#phi^{pair}-#Psi^{TPC}",                             ""},
   {"v0ArpH2FlowV2",          "cos(2(#phi^{pair}-#Psi^{V0A}))",                     ""},
   {"v0CrpH2FlowV2",          "cos(2(#phi^{pair}-#Psi^{V0C}))",                     ""},
   {"v0ACrpH2FlowV2",         "cos(2(#phi^{pair}-#Psi^{V0AC}))",                    ""},
+  {"TPCrpH2FlowV2",          "cos(2(#phi^{pair}-#Psi^{TPC}))",                     ""},
+  {"TPCrpH2FlowV2Sin",       "sin(2(#phi^{pair}-#Psi^{TPC}))",                     ""},
   {"LegDistance",            "d^{legs}",                                           "(cm)"},
   {"LegDistanceXY",          "d^{legs}_{xy}",                                      "(cm)"},
   {"DeltaEta",               "#Delta #eta",                                        ""},
@@ -183,12 +209,17 @@ const char* AliDielectronVarManager::fgkParticleNames[AliDielectronVarManager::k
   {"PseudoProperTimeErr",    "#tau error",                                         "(#mus)"},
   {"PseudoProperTimeResolution", "(#tau-#tau^{MC truth})",                         "(#mus)"},
   {"PseudoProperTimePull",   "#frac{(#tau-#tau^{MC truth})}{#tau error}",          ""},
-  {"TRDpidEffPair",          "#epsilon^{TRD pair}(PID)",                           ""},
+  {"TRDpidEffPair",          "#varepsilon^{TRD pair}(PID)",                           ""},
   {"MomAsymDau1",            "#it{p}^{leg1}/#it{p}^{pair}",                        ""},
   {"MomAsymDau2",            "#it{p}^{leg2}/#it{p}^{pair}",                        ""},
+  {"EffPair",                "A#times#varepsilon",                                 ""},
+  {"OneOverEffPair",         "(A#times#varepsilon)^{-1}",                          ""},
+  {"OneOverEffPairSq",       "(A#times#varepsilon)^{-2}",                          ""},
+  {"RndmPair",               "P",                                                  ""},
+  {"Pairs",                  "pairs/event",                                        ""},
   //
-  {"X",                      "x_{prim.vtx}",                                       "(cm)"},
-  {"Y",                      "y_{prim.vtx}",                                       "(cm)"},
+  {"Xprim",                  "x_{prim.vtx}",                                       "(cm)"},
+  {"Yprim",                  "y_{prim.vtx}",                                       "(cm)"},
   {"Z",                      "z_{prim.vtx}",                                       "(cm)"},
   {"XRes",                   "#Delta x_{prim.vtx}",                                "(cm)"},
   {"YRes",                   "#Delta y_{prim.vtx}",                                "(cm)"},
@@ -225,10 +256,13 @@ const char* AliDielectronVarManager::fgkParticleNames[AliDielectronVarManager::k
   {"MultV0A",                "multiplicity V0A",                                   ""},
   {"MultV0C",                "multiplicity V0C",                                   ""},
   {"MultV0",                 "multiplicity V0",                                    ""},
-  {"AdcV0A",                   "AdcV0A",                                   ""},
-  {"AdcV0C",                   "AdcV0C",                                   ""},
-  {"AdcV0",    "AdcV0",                                     ""},
-  {"VZERO_ch0",   "VZERO_ch0",                                    ""},
+  {"EqMultV0A",              "equalized multiplicity V0A",                         ""},
+  {"EqMultV0C",              "equalized multiplicity V0C",                         ""},
+  {"EqMultV0",               "equalized multiplicity V0",                          ""},
+  {"AdcV0A",                 "AdcV0A",                                             ""},
+  {"AdcV0C",                 "AdcV0C",                                             ""},
+  {"AdcV0",                  "AdcV0",                                              ""},
+  {"VZERO_ch0",   "VZERO_ch0", ""},
   {"VZERO_ch1",   "VZERO_ch1", ""},
   {"VZERO_ch2",   "VZERO_ch2", ""},
   {"VZERO_ch3",   "VZERO_ch3", ""},
@@ -316,8 +350,10 @@ const char* AliDielectronVarManager::fgkParticleNames[AliDielectronVarManager::k
   {"V0AV0CDiffH2",           "cos(2(#Psi^{V0A}-#Psi^{V0C}))",                      ""},
   {"TPCxH2",                 "Q_{x}^{TPC}",                                        ""},
   {"TPCyH2",                 "Q_{y}^{TPC}",                                        ""},
-  {"TPCmagH2",               "|#vec{Q}^{TPC}|",                                   ""},
+  {"TPCmagH2",               "|#vec{Q}^{TPC}|",                                    ""},
   {"TPCrpH2",                "#Psi^{TPC}",                                         ""},
+  {"CosTPCrpH2",             "cos(2#Psi^{TPC})",                                   ""},
+  {"SinTPCrpH2",             "sin(2#Psi^{TPC})",                                   ""},
   {"TPCsub1xH2",             "Q_{x}^{TPCsub1}",                                    ""},
   {"TPCsub1yH2",             "Q_{y}^{TPCsub1}",                                    ""},
   {"TPCsub1rpH2",            "#Psi^{TPCsub1}",                                     ""},
@@ -339,8 +375,14 @@ const char* AliDielectronVarManager::fgkParticleNames[AliDielectronVarManager::k
   {"TPCsub2rpH2uc",          "#Psi^{TPCsub2} (uncorr.)",                           ""},
   {"TPCsub12DiffH2uc",       "cos(2(#Psi^{TPCsub1}-#Psi^{TPCsub2})) (uncorr.)",    ""},
 
+  {"ZDCArpH1",               "#Psi^{ZDCA}"                                         ""},
+  {"ZDCCrpH1",               "#Psi^{ZDCA}"                                         ""},
+  {"ZDCACrpH1",              "#Psi^{ZDCA}"                                         ""},
+  {"kZDCrpResH1",            ""                                                    ""},
+  {"kv0ZDCrpRes",            ""                                                    ""},
+
   {"NTrk",                   "N_{trk}",                                            ""},
-  {"Tracks",                 "tracks",                                             ""},
+  {"Tracks",                 "tracks/per event",                                             ""},
   {"NVtxContrib",            "N_{vtx. contrib.}",                                  ""},
   {"NVtxContribTPC",         "N_{vtx. contrib.}^{TPC}",                            ""},
   {"Nacc",                   "N_{acc} #cbar_{#||{#eta}<0.9}",                      ""},
@@ -374,6 +416,9 @@ const char* AliDielectronVarManager::fgkParticleNames[AliDielectronVarManager::k
   {"Nch10",                  "N_{ch} #cbar_{#||{#eta}<1.0}",                       ""},
   {"Centrality",             "centrality",                                         "(%)"},
   {"CentralitySPD",          "centrality_{SPD}",                                   "(%)"},
+  {"TriggerInclONL",         "online trigger bit (inclusive)",                     ""},
+  {"TriggerInclOFF",         "offline trigger bit (inclusive)",                    ""},
+  {"TriggerExclOFF",         "offline trigger bit (exclusive)",                    ""},
   {"Nevents",                "N_{evt}",                                            ""},
   {"RunNumber",              "run",                                                ""},
   {"MixingBin",              "mixing bin",                                         ""}
@@ -385,11 +430,16 @@ AliEventplane*  AliDielectronVarManager::fgTPCEventPlane    = 0x0;
 AliKFVertex*    AliDielectronVarManager::fgKFVertex         = 0x0;
 TProfile*       AliDielectronVarManager::fgMultEstimatorAvg[4][9] = {{0x0}};
 TH3D*           AliDielectronVarManager::fgTRDpidEff[10][4] = {{0x0}};
+TObject*        AliDielectronVarManager::fgLegEffMap           = 0x0;
+TObject*        AliDielectronVarManager::fgPairEffMap          = 0x0;
+TBits*          AliDielectronVarManager::fgFillMap          = 0x0;
 Double_t        AliDielectronVarManager::fgTRDpidEffCentRanges[10][4] = {{0.0}};
 TString         AliDielectronVarManager::fgVZEROCalibrationFile = "";
 TString         AliDielectronVarManager::fgVZERORecenteringFile = "";
+TString         AliDielectronVarManager::fgZDCRecenteringFile = "";
 TProfile2D*     AliDielectronVarManager::fgVZEROCalib[64] = {0x0};
 TProfile2D*     AliDielectronVarManager::fgVZERORecentering[2][2] = {{0x0,0x0},{0x0,0x0}};
+TProfile3D*     AliDielectronVarManager::fgZDCRecentering[3][2] = {{0x0,0x0},{0x0,0x0},{0x0,0x0}};
 Int_t           AliDielectronVarManager::fgCurrentRun = -1;
 Double_t        AliDielectronVarManager::fgData[AliDielectronVarManager::kNMaxValues] = {0.};
 //________________________________________________________________
@@ -409,6 +459,10 @@ AliDielectronVarManager::AliDielectronVarManager() :
   for(Int_t i=0; i<2; ++i) {
     for(Int_t j=0; j<2; ++j) fgVZERORecentering[i][j] = 0x0;
   }
+  for(Int_t i=0; i<3; ++i)
+    for(Int_t j=0; j<2; ++j) fgZDCRecentering[i][j] = 0x0;
+
+  gRandom->SetSeed();
 }
 
 //________________________________________________________________
@@ -428,6 +482,11 @@ AliDielectronVarManager::AliDielectronVarManager(const char* name, const char* t
   for(Int_t i=0; i<2; ++i)
     for(Int_t j=0; j<2; ++j) 
       fgVZERORecentering[i][j] = 0x0;
+  for(Int_t i=0; i<3; ++i)
+    for(Int_t j=0; j<2; ++j) fgZDCRecentering[i][j] = 0x0;
+
+
+  gRandom->SetSeed();
 }
 
 //________________________________________________________________
@@ -447,6 +506,10 @@ AliDielectronVarManager::~AliDielectronVarManager()
   for(Int_t i=0; i<2; ++i)
     for(Int_t j=0; j<2; ++j) 
       if(fgVZERORecentering[i][j]) delete fgVZERORecentering[i][j]; 
+  for(Int_t i=0; i<3; ++i)
+    for(Int_t j=0; j<2; ++j)
+      if(fgZDCRecentering[i][j]) delete fgZDCRecentering[i][j];
+
 }
 
 //________________________________________________________________