Skip to content

Commit dc840db

Browse files
YazhenLinYazhenLin
authored
Add files via upload
1 parent 47c97e8 commit dc840db

File tree

4 files changed

+419
-20
lines changed

4 files changed

+419
-20
lines changed

‎PWGDQ/Core/HistogramsLibrary.cxx‎

Lines changed: 78 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -11,12 +11,15 @@
1111
//
1212
// Contact: iarsene@cern.ch, i.c.arsene@fys.uio.no
1313
//
14-
#include <vector>
15-
#include <algorithm>
1614
#include "PWGDQ/Core/HistogramsLibrary.h"
15+
1716
#include "VarManager.h"
17+
1818
#include "CommonConstants/MathConstants.h"
1919

20+
#include <algorithm>
21+
#include <vector>
22+
2023
void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* histClass, const char* groupName, const char* subGroupName)
2124
{
2225
//
@@ -959,7 +962,54 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
959962
hm->AddHistogram(histClass, "Weight", "", false, 50, 0.0, 5.0, VarManager::kMCParticleWeight);
960963
hm->AddHistogram(histClass, "MCImpPar_CentFT0CMC", "MC impact param vs MC Cent. FT0C", false, 20, 0.0, 20.0, VarManager::kMCEventImpParam, 100, 0.0, 100.0, VarManager::kMCEventCentrFT0C);
961964
}
965+
if (!groupStr.CompareTo("energy_correlator_gen")) {
966+
hm->AddHistogram(histClass, "MCCostheta", "Cos#theta", false, 40, -1.0, 1.0, VarManager::kMCCosTheta, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kMCWeight_before); ////
967+
hm->AddHistogram(histClass, "MCHadronPdgCode", "MCHadronPdgCode", false, 6000, -3000, 3000, VarManager::kMCHadronPdgCode); ////
968+
hm->AddHistogram(histClass, "MCMotherPdgCode", "MCMotherPdgCode", false, 6000, -3000, 3000, VarManager::kMCMotherPdgCode); ////
969+
hm->AddHistogram(histClass, "MCPdgCode", "MCPdgCode", false, 1000, -1000, 1000, VarManager::kMCPdgCode); ////
970+
hm->AddHistogram(histClass, "Coschi_pion", "", false, 40, -1.0, 1.0, VarManager::kMCCosChi_pion, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kMCWeight_pion);
971+
hm->AddHistogram(histClass, "Coschi_hadron", "", false, 40, -1.0, 1.0, VarManager::kMCCosChi_hadron, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kMCWeight_hadron);
972+
hm->AddHistogram(histClass, "Pt_Hadron", "", false, 120, 0.0, 30.0, VarManager::kMCHadronPt);
973+
hm->AddHistogram(histClass, "Eta_Hadron", "", false, 120, -2.0, 2.0, VarManager::kMCHadronEta);
974+
hm->AddHistogram(histClass, "DeltaEta", "", false, 20, -2.0, 2.0, VarManager::kMCdeltaeta);
975+
hm->AddHistogram(histClass, "DeltaPhi", "", false, 50, -8.0, 8.0, VarManager::kMCdeltaphi);
976+
hm->AddHistogram(histClass, "DeltaEta_DeltaPhi", "", false, 20, -2.0, 2.0, VarManager::kMCdeltaeta, 50, -8.0, 8.0, VarManager::kMCdeltaphi);
977+
hm->AddHistogram(histClass, "Coschi_hadron_kMCWeight_hadron_DeltaEta", "", false, 40, -1.0, 1.0, VarManager::kMCCosChi_hadron, 50, -8.0, 8.0, VarManager::kMCdeltaeta, 120, 0.0, 50, VarManager::kMCWeight_hadron);
978+
hm->AddHistogram(histClass, "Costheta_hadron_kMCWeight_before_DeltaEta", "", false, 40, -1.0, 1.0, VarManager::kMCCosTheta, 50, -8.0, 8.0, VarManager::kMCdeltaeta, 120, 0.0, 50, VarManager::kMCWeight_before);
979+
// for bkg
980+
981+
hm->AddHistogram(histClass, "DeltaPhi_minus", "", false, 50, -8.0, 8.0, VarManager::kMCdeltaphi_minus);
982+
hm->AddHistogram(histClass, "DeltaPhi_toward_minus", "", false, 50, -8.0, 8.0, VarManager::kMCdeltaphi_toward_minus);
983+
hm->AddHistogram(histClass, "DeltaPhi_away_minus", "", false, 50, -8.0, 8.0, VarManager::kMCdeltaphi_away_minus);
984+
hm->AddHistogram(histClass, "DeltaPhi_trans_minus", "", false, 50, -8.0, 8.0, VarManager::kMCdeltaphi_trans_minus);
985+
986+
hm->AddHistogram(histClass, "Coschi_minus", "", false, 40, -1.0, 1.0, VarManager::kMCCosChi_minus, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kMCWeight_minus);
987+
hm->AddHistogram(histClass, "Coschi_toward_minus", "", false, 40, -1.0, 1.0, VarManager::kMCCosChi_toward_minus, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kMCWeight_toward_minus);
988+
hm->AddHistogram(histClass, "Coschi_away_minus", "", false, 40, -1.0, 1.0, VarManager::kMCCosChi_away_minus, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kMCWeight_away_minus);
989+
hm->AddHistogram(histClass, "Coschi_trans_minus", "", false, 40, -1.0, 1.0, VarManager::kMCCosChi_trans_minus, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kMCWeight_trans_minus);
990+
991+
hm->AddHistogram(histClass, "Costheta_minus", "", false, 40, -1.0, 1.0, VarManager::kMCCosTheta_minus, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kMCWeight_before);
992+
hm->AddHistogram(histClass, "Costheta_toward_minus", "", false, 40, -1.0, 1.0, VarManager::kMCCosTheta_toward_minus, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kMCWeight_before);
993+
hm->AddHistogram(histClass, "Costheta_away_minus", "", false, 40, -1.0, 1.0, VarManager::kMCCosTheta_away_minus, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kMCWeight_before);
994+
hm->AddHistogram(histClass, "Costheta_trans_minus", "", false, 40, -1.0, 1.0, VarManager::kMCCosTheta_trans_minus, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kMCWeight_before);
995+
996+
hm->AddHistogram(histClass, "DeltaPhi_plus", "", false, 50, -8.0, 8.0, VarManager::kMCdeltaphi_plus);
997+
hm->AddHistogram(histClass, "DeltaPhi_toward_plus", "", false, 50, -8.0, 8.0, VarManager::kMCdeltaphi_toward_plus);
998+
hm->AddHistogram(histClass, "DeltaPhi_away_plus", "", false, 50, -8.0, 8.0, VarManager::kMCdeltaphi_away_plus);
999+
hm->AddHistogram(histClass, "DeltaPhi_trans_plus", "", false, 50, -8.0, 8.0, VarManager::kMCdeltaphi_trans_plus);
1000+
1001+
hm->AddHistogram(histClass, "Coschi_plus", "", false, 40, -1.0, 1.0, VarManager::kMCCosChi_plus, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kMCWeight_plus);
1002+
hm->AddHistogram(histClass, "Coschi_toward_plus", "", false, 40, -1.0, 1.0, VarManager::kMCCosChi_toward_plus, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kMCWeight_toward_plus);
1003+
hm->AddHistogram(histClass, "Coschi_away_plus", "", false, 40, -1.0, 1.0, VarManager::kMCCosChi_away_plus, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kMCWeight_away_plus);
1004+
hm->AddHistogram(histClass, "Coschi_trans_plus", "", false, 40, -1.0, 1.0, VarManager::kMCCosChi_trans_plus, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kMCWeight_trans_plus);
1005+
1006+
hm->AddHistogram(histClass, "Costheta_plus", "", false, 40, -1.0, 1.0, VarManager::kMCCosTheta_plus, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kMCWeight_before);
1007+
hm->AddHistogram(histClass, "Costheta_toward_plus", "", false, 40, -1.0, 1.0, VarManager::kMCCosTheta_toward_plus, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kMCWeight_before);
1008+
hm->AddHistogram(histClass, "Costheta_away_plus", "", false, 40, -1.0, 1.0, VarManager::kMCCosTheta_away_plus, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kMCWeight_before);
1009+
hm->AddHistogram(histClass, "Costheta_trans_plus", "", false, 40, -1.0, 1.0, VarManager::kMCCosTheta_trans_plus, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kMCWeight_before);
1010+
}
9621011

1012+
// }
9631013
if (!groupStr.CompareTo("pair")) {
9641014
if (subGroupStr.Contains("cepf")) {
9651015
hm->AddHistogram(histClass, "Mass", "", false, 300, 0.0, 12.0, VarManager::kMass);
@@ -999,6 +1049,8 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
9991049
hm->AddHistogram(histClass, "Eta_Pt", "", false, 40, -2.0, 2.0, VarManager::kEta, 40, 0.0, 20.0, VarManager::kPt);
10001050
hm->AddHistogram(histClass, "Y_Pt", "", false, 40, -2.0, 2.0, VarManager::kRap, 40, 0.0, 20.0, VarManager::kPt);
10011051
hm->AddHistogram(histClass, "Mass_VtxZ", "", true, 30, -15.0, 15.0, VarManager::kVtxZ, 500, 0.0, 5.0, VarManager::kMass);
1052+
hm->AddHistogram(histClass, "Mass_Y_Pt", "", false, 125, 0.0, 5.0, VarManager::kMass, 40, -2.0, 2.0, VarManager::kRap, 40, 0.0, 20.0, VarManager::kPt);
1053+
hm->AddHistogram(histClass, "Y", ";y", false, 40, -2.0, 2.0, VarManager::kRap);
10021054
if (subGroupStr.Contains("pbpb")) {
10031055
hm->AddHistogram(histClass, "Mass_CentFT0C", "", false, 125, 0.0, 5.0, VarManager::kMass, 20, 0.0, 100.0, VarManager::kCentFT0C);
10041056
hm->AddHistogram(histClass, "Pt_CentFT0C", "", false, 100, 0.0, 10.0, VarManager::kPt, 20, 0.0, 100.0, VarManager::kCentFT0C);
@@ -1799,6 +1851,18 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
17991851
hm->AddHistogram(histClass, "Mass", "", false, 750, 0.0, 30.0, VarManager::kPairMass);
18001852
hm->AddHistogram(histClass, "Pt", "", false, 750, 0.0, 30.0, VarManager::kPairPt);
18011853
}
1854+
if (subGroupStr.Contains("mixedevent_energy_correlator")) { // for mixed event energy_correlator
1855+
hm->AddHistogram(histClass, "Coschi", "", false, 40, -1.0, 1.0, VarManager::kCosChi, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kECWeight);
1856+
hm->AddHistogram(histClass, "CosTheta", "", false, 40, -1.0, 1.0, VarManager::kCosTheta, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kEWeight_before);
1857+
hm->AddHistogram(histClass, "Pt_Hadron", ";P_{T}", false, 120, 0.0, 30.0, VarManager::kPtDau);
1858+
hm->AddHistogram(histClass, "Eta_Hadron", ";#eta", false, 120, -2.0, 2.0, VarManager::kEtaDau);
1859+
hm->AddHistogram(histClass, "Phi_Hadron", ";#phi", false, 120, -8, 8, VarManager::kPhiDau);
1860+
hm->AddHistogram(histClass, "DeltaEta_DeltaPhi", "", false, 20, -2.0, 2.0, VarManager::kDeltaEta, 50, -8.0, 8.0, VarManager::kDeltaPhi);
1861+
hm->AddHistogram(histClass, "DeltaEta", "", false, 20, -2.0, 2.0, VarManager::kDeltaEta);
1862+
hm->AddHistogram(histClass, "DeltaPhi", "", false, 50, -8.0, 8.0, VarManager::kDeltaPhi);
1863+
hm->AddHistogram(histClass, "Coschi_hadron_kMCWeight_hadron_DeltaEta", "", false, 40, -1.0, 1.0, VarManager::kCosChi, 50, -8.0, 8.0, VarManager::kDeltaEta, 120, 0.0, 50, VarManager::kECWeight);
1864+
hm->AddHistogram(histClass, "Costheta_hadron_kMCWeight_before_DeltaEta", "", false, 40, -1.0, 1.0, VarManager::kCosTheta, 50, -8.0, 8.0, VarManager::kDeltaEta, 120, 0.0, 50, VarManager::kEWeight_before);
1865+
}
18021866
if (subGroupStr.Contains("invmass")) {
18031867
hm->AddHistogram(histClass, "Mass_Dilepton", "", false, 125, 0.0, 5.0, VarManager::kPairMassDau);
18041868
hm->AddHistogram(histClass, "Mass_Hadron", "", false, 125, 0.0, 5.0, VarManager::kMassDau);
@@ -1898,6 +1962,18 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
18981962
if (subGroupStr.Contains("opencharm")) {
18991963
hm->AddHistogram(histClass, "Delta_Mass_DstarD0region", "", false, 50, 0.14, 0.16, VarManager::kDeltaMass);
19001964
}
1965+
if (subGroupStr.Contains("energy_correlator")) {
1966+
hm->AddHistogram(histClass, "Coschi", "", false, 40, -1.0, 1.0, VarManager::kCosChi, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kECWeight);
1967+
hm->AddHistogram(histClass, "CosTheta", "", false, 40, -1.0, 1.0, VarManager::kCosTheta, 0, 0, 0, -1, 0, 0, 0, -1, "", "", "", -1, VarManager::kEWeight_before);
1968+
hm->AddHistogram(histClass, "Pt_Hadron", ";P_{T}", false, 120, 0.0, 30.0, VarManager::kPtDau);
1969+
hm->AddHistogram(histClass, "Eta_Hadron", ";#eta", false, 120, -2.0, 2.0, VarManager::kEtaDau);
1970+
hm->AddHistogram(histClass, "Phi_Hadron", ";#phi", false, 120, -8, 8, VarManager::kPhiDau);
1971+
hm->AddHistogram(histClass, "DeltaEta_DeltaPhi", "", false, 20, -2.0, 2.0, VarManager::kDeltaEta, 50, -8.0, 8.0, VarManager::kDeltaPhi);
1972+
hm->AddHistogram(histClass, "DeltaEta", "", false, 20, -2.0, 2.0, VarManager::kDeltaEta);
1973+
hm->AddHistogram(histClass, "DeltaPhi", "", false, 50, -8.0, 8.0, VarManager::kDeltaPhi);
1974+
hm->AddHistogram(histClass, "Coschi_hadron_kMCWeight_hadron_DeltaEta", "", false, 40, -1.0, 1.0, VarManager::kCosChi, 50, -8.0, 8.0, VarManager::kDeltaEta, 120, 0.0, 50, VarManager::kECWeight);
1975+
hm->AddHistogram(histClass, "Costheta_hadron_kMCWeight_before_DeltaEta", "", false, 40, -1.0, 1.0, VarManager::kCosTheta, 50, -8.0, 8.0, VarManager::kDeltaEta, 120, 0.0, 50, VarManager::kEWeight_before);
1976+
}
19011977
}
19021978

19031979
if (!groupStr.CompareTo("dilepton-charmhadron")) {

‎PWGDQ/Core/MCSignalLibrary.cxx‎

Lines changed: 56 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -15,11 +15,13 @@
1515
#include <vector>
1616
// #include <iostream>
1717

18-
#include <TPDGCode.h>
19-
#include "CommonConstants/PhysicsConstants.h"
2018
#include "PWGDQ/Core/MCSignalLibrary.h"
19+
20+
#include "CommonConstants/PhysicsConstants.h"
2121
#include "Framework/Logger.h"
2222

23+
#include <TPDGCode.h>
24+
2325
using namespace o2::constants::physics;
2426
// using std::cout;
2527
// using std::endl;
@@ -882,6 +884,22 @@ MCSignal* o2::aod::dqmcsignals::GetMCSignal(const char* name)
882884
signal = new MCSignal(name, "ee pairs from j/psi decays", {prong, prong}, {1, 1}); // signal at pair level
883885
return signal;
884886
}
887+
/////
888+
if (!nameStr.compare("eeFrompromptJpsi")) {
889+
MCProng prong(2, {11, 443}, {true, true}, {false, false}, {0, 0}, {0, 0}, {false, false}, false, {503}, {true});
890+
prong.SetSourceBit(0, MCProng::kPhysicalPrimary);
891+
signal = new MCSignal(name, "ee pairs from prompt j/psi decays", {prong, prong}, {1, 1}); // signal at pair level
892+
return signal;
893+
}
894+
895+
if (!nameStr.compare("eeFromnonpromptJpsi")) {
896+
MCProng prong(3, {11, 443, 503}, {true, true, true}, {false, false, false}, {0, 0, 0}, {0, 0, 0}, {false, false, false});
897+
prong.SetSourceBit(0, MCProng::kPhysicalPrimary);
898+
signal = new MCSignal(name, "ee pairs from nonprompt j/psi decays", {prong, prong}, {2, 2}); // signal at pair level
899+
return signal;
900+
}
901+
902+
//////
885903
if (!nameStr.compare("eeFromJpsiExclusive")) {
886904
MCProng prong(2, {11, 443}, {true, true}, {false, false}, {0, 0}, {0, 0}, {false, false});
887905
prong.SetSourceBit(0, MCProng::kPhysicalPrimary);
@@ -1885,6 +1903,42 @@ MCSignal* o2::aod::dqmcsignals::GetMCSignal(const char* name)
18851903
signal = new MCSignal(name, "Electron pair and pion pair from Psi2S", {pronge, pronge, prongPi, prongPi}, {2, 2, 1, 1});
18861904
return signal;
18871905
}
1906+
1907+
if (!nameStr.compare("eeFromPromptJpsiAnyPrimary")) {
1908+
MCProng pronge(2, {11, 443}, {true, true}, {false, false}, {0, 0}, {0, 0}, {false, false}, false, {503}, {true});
1909+
pronge.SetSourceBit(0, MCProng::kPhysicalPrimary);
1910+
MCProng prongPrimary(1);
1911+
prongPrimary.SetSourceBit(0, MCProng::kPhysicalPrimary);
1912+
signal = new MCSignal(name, "anyprimary and electron pair from prompt jpsi", {pronge, pronge, prongPrimary}, {1, 1, -1});
1913+
return signal;
1914+
}
1915+
1916+
if (!nameStr.compare("eeFromJpsiAnyPrimary")) {
1917+
MCProng pronge(2, {11, 443}, {true, true}, {false, false}, {0, 0}, {0, 0}, {false, false});
1918+
pronge.SetSourceBit(0, MCProng::kPhysicalPrimary);
1919+
MCProng prongPrimary(1);
1920+
prongPrimary.SetSourceBit(0, MCProng::kPhysicalPrimary);
1921+
signal = new MCSignal(name, "anyprimary and electron pair from prompt jpsi", {pronge, pronge, prongPrimary}, {1, 1, -1});
1922+
return signal;
1923+
}
1924+
1925+
if (!nameStr.compare("eeFromNonPromptJpsiAnyPrimary")) {
1926+
MCProng pronge(3, {11, 443, 503}, {true, true, true}, {false, false, false}, {0, 0, 0}, {0, 0, 0}, {false, false, false});
1927+
pronge.SetSourceBit(0, MCProng::kPhysicalPrimary);
1928+
MCProng prongPrimary(1);
1929+
prongPrimary.SetSourceBit(0, MCProng::kPhysicalPrimary);
1930+
signal = new MCSignal(name, "anyprimary and electron pair from non-prompt jpsi", {pronge, pronge, prongPrimary}, {2, 2, -1});
1931+
return signal;
1932+
}
1933+
1934+
// if (!nameStr.compare("eeFromNonPromptJpsiAnyPrimary")) {
1935+
// MCProng pronge(2, {11, 443}, {true, true}, {false, false}, {0, 0}, {0, 0}, {false, false}, false, {503}, {false});
1936+
// pronge.SetSourceBit(0, MCProng::kPhysicalPrimary);
1937+
// MCProng prongPrimary(1);
1938+
// prongPrimary.SetSourceBit(0, MCProng::kPhysicalPrimary);
1939+
// signal = new MCSignal(name, "anyprimary and electron pair from non-prompt jpsi", {pronge, pronge, prongPrimary}, {-1, -1, -1});
1940+
// return signal;
1941+
//}
18881942
return nullptr;
18891943
}
18901944

‎PWGDQ/Core/VarManager.cxx‎

Lines changed: 55 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -8,12 +8,14 @@
88
// In applying this license CERN does not waive the privileges and immunities
99
// granted to it by virtue of its status as an Intergovernmental Organization
1010
// or submit itself to any jurisdiction.
11+
#include "PWGDQ/Core/VarManager.h"
12+
13+
#include "Tools/KFparticle/KFUtilities.h"
14+
1115
#include <cmath>
1216
#include <iostream>
13-
#include <vector>
1417
#include <map>
15-
#include "PWGDQ/Core/VarManager.h"
16-
#include "Tools/KFparticle/KFUtilities.h"
18+
#include <vector>
1719

1820
using std::cout;
1921
using std::endl;
@@ -29,7 +31,7 @@ bool VarManager::fgUsedKF = false;
2931
float VarManager::fgMagField = 0.5;
3032
float VarManager::fgzMatching = -77.5;
3133
float VarManager::fgValues[VarManager::kNVars] = {0.0f};
32-
float VarManager::fgTPCInterSectorBoundary = 1.0; // cm
34+
float VarManager::fgTPCInterSectorBoundary = 1.0; // cm
3335
int VarManager::fgITSROFbias = 0;
3436
int VarManager::fgITSROFlength = 100;
3537
int VarManager::fgITSROFBorderMarginLow = 0;
@@ -728,6 +730,28 @@ void VarManager::SetDefaultVarNames()
728730
fgVariableUnits[kMuonTimeRes] = "ns";
729731
fgVariableNames[kMCPdgCode] = "MC PDG code";
730732
fgVariableUnits[kMCPdgCode] = "";
733+
/////
734+
fgVariableNames[kMCCosTheta] = "Cos#theta";
735+
fgVariableUnits[kMCCosTheta] = "";
736+
fgVariableNames[kMCHadronPdgCode] = "HadronPdgCode";
737+
fgVariableUnits[kMCHadronPdgCode] = "";
738+
fgVariableNames[kMCCosChi_pion] = "Cos#chi";
739+
fgVariableUnits[kMCCosChi_pion] = "";
740+
fgVariableNames[kMCCosChi_hadron] = "Cos#chi";
741+
fgVariableUnits[kMCCosChi_hadron] = "";
742+
fgVariableNames[kMCJpsiPt] = "Jpsi p_{T}";
743+
fgVariableUnits[kMCJpsiPt] = "GeV/c";
744+
fgVariableNames[kMCHadronPt] = "Hadron p_{T}";
745+
fgVariableUnits[kMCHadronPt] = "GeV/c";
746+
fgVariableNames[kMCHadronEta] = "Hadron #eta";
747+
fgVariableUnits[kMCHadronEta] = "";
748+
fgVariableNames[kMCdeltaphi] = "#Delta#phi";
749+
fgVariableUnits[kMCdeltaphi] = "";
750+
fgVariableNames[kMCdeltaeta] = "#Delta#eta";
751+
fgVariableUnits[kMCdeltaeta] = "";
752+
fgVariableNames[kNhadron] = "N_{hadron}";
753+
fgVariableUnits[kNhadron] = "";
754+
/////
731755
fgVariableNames[kMCParticleWeight] = "MC particle weight";
732756
fgVariableUnits[kMCParticleWeight] = "";
733757
fgVariableNames[kMCPx] = "MC px";
@@ -1156,6 +1180,16 @@ void VarManager::SetDefaultVarNames()
11561180
fgVariableUnits[kDeltaPhi] = "rad.";
11571181
fgVariableNames[kDeltaPhiSym] = "#Delta#phi";
11581182
fgVariableUnits[kDeltaPhiSym] = "rad.";
1183+
fgVariableNames[kCosChi] = "Cos#chi";
1184+
fgVariableUnits[kCosChi] = "";
1185+
fgVariableNames[kCosTheta] = "Cos#theta";
1186+
fgVariableUnits[kCosTheta] = "";
1187+
fgVariableNames[kPtDau] = "hadron P_{T}";
1188+
fgVariableUnits[kPtDau] = "GeV/c";
1189+
fgVariableNames[kEtaDau] = "hadron #eta";
1190+
fgVariableUnits[kEtaDau] = "";
1191+
fgVariableNames[kPhiDau] = "hadron #phi";
1192+
fgVariableUnits[kPhiDau] = "";
11591193
fgVariableNames[kCosThetaHE] = "cos#it{#theta}";
11601194
fgVariableUnits[kCosThetaHE] = "";
11611195
fgVariableNames[kPhiHE] = "#varphi_{HE}";
@@ -1685,6 +1719,16 @@ void VarManager::SetDefaultVarNames()
16851719
fgVarNamesMap["kMuonNAssocsOutOfBunch"] = kMuonNAssocsOutOfBunch;
16861720
fgVarNamesMap["kNMuonTrackVariables"] = kNMuonTrackVariables;
16871721
fgVarNamesMap["kMCPdgCode"] = kMCPdgCode;
1722+
//////
1723+
fgVarNamesMap["kMCCosTheta"] = kMCCosTheta;
1724+
fgVarNamesMap["kMCHadronPdgCode"] = kMCHadronPdgCode;
1725+
fgVarNamesMap["kMCCosChi_pion"] = kMCCosChi_pion;
1726+
fgVarNamesMap["kMCCosChi_hadron"] = kMCCosChi_hadron;
1727+
fgVarNamesMap["kMCHadronPt"] = kMCHadronPt;
1728+
fgVarNamesMap["kMCCosChi_minus"] = kMCCosChi_minus;
1729+
fgVarNamesMap["kMCWeight_before"] = kMCWeight_before;
1730+
1731+
/////
16881732
fgVarNamesMap["kMCParticleWeight"] = kMCParticleWeight;
16891733
fgVarNamesMap["kMCPx"] = kMCPx;
16901734
fgVarNamesMap["kMCPy"] = kMCPy;
@@ -1890,6 +1934,13 @@ void VarManager::SetDefaultVarNames()
18901934
fgVarNamesMap["kDeltaEta"] = kDeltaEta;
18911935
fgVarNamesMap["kDeltaPhi"] = kDeltaPhi;
18921936
fgVarNamesMap["kDeltaPhiSym"] = kDeltaPhiSym;
1937+
fgVarNamesMap["kCosTheta"] = kCosTheta;
1938+
fgVarNamesMap["kCosChi"] = kCosChi;
1939+
fgVarNamesMap["kECWeight"] = kECWeight;
1940+
fgVarNamesMap["kEWeight_before"] = kEWeight_before;
1941+
fgVarNamesMap["kPtDau"] = kPtDau;
1942+
fgVarNamesMap["kEtaDau"] = kEtaDau;
1943+
fgVarNamesMap["kPhiDau"] = kPhiDau;
18931944
fgVarNamesMap["kNCorrelationVariables"] = kNCorrelationVariables;
18941945
fgVarNamesMap["kQuadMass"] = kQuadMass;
18951946
fgVarNamesMap["kQuadDefaultDileptonMass"] = kQuadDefaultDileptonMass;

0 commit comments

Comments
 (0)