Skip to content

Commit 5937038

Browse files
alibuildalibuild
committed
Please consider the following formatting changes
1 parent f314069 commit 5937038

File tree

3 files changed

+269
-263
lines changed

3 files changed

+269
-263
lines changed

‎PWGDQ/Core/VarManager.h‎

Lines changed: 35 additions & 32 deletions
Original file line numberDiff line numberDiff line change
@@ -1104,9 +1104,10 @@ class VarManager : public TObject
11041104
return false;
11051105
}
11061106

1107-
// Flag to set PV recalculation via KF
1108-
static void SetPVrecalculationKF(const bool pvRecalKF) {
1109-
fgPVrecalKF = pvRecalKF;
1107+
// Flag to set PV recalculation via KF
1108+
static void SetPVrecalculationKF(const bool pvRecalKF)
1109+
{
1110+
fgPVrecalKF = pvRecalKF;
11101111
}
11111112

11121113
// Setup the collision system
@@ -1311,7 +1312,7 @@ class VarManager : public TObject
13111312
static void FillQuadMC(T1 const& t1, T2 const& t2, T2 const& t3, float* values = nullptr);
13121313
template <int pairType, uint32_t collFillMap, uint32_t fillMap, typename C, typename T>
13131314
static void FillPairVertexing(C const& collision, T const& t1, T const& t2, bool propToSV = false, float* values = nullptr);
1314-
template <int pairType, uint32_t collFillMap, uint32_t fillMap, typename C, typename T>
1315+
template <int pairType, uint32_t collFillMap, uint32_t fillMap, typename C, typename T>
13151316
static void FillPairVertexingRecomputePV(C const& /*collision*/, T const& t1, T const& t2, o2::dataformats::VertexBase pvRefitted, float* values = nullptr);
13161317
template <uint32_t collFillMap, uint32_t fillMap, typename C, typename T>
13171318
static void FillTripletVertexing(C const& collision, T const& t1, T const& t2, T const& t3, PairCandidateType tripletType, float* values = nullptr);
@@ -4202,7 +4203,7 @@ void VarManager::FillPairVertexing(C const& collision, T const& t1, T const& t2,
42024203

42034204
if constexpr ((pairType == kDecayToEE || pairType == kDecayToKPi) && trackHasCov) {
42044205
secondaryVertex = fgFitterTwoProngBarrel.getPCACandidate();
4205-
// printf("secVtx (first) %f %f %f \n",secondaryVertex[0],secondaryVertex[1],secondaryVertex[2]);
4206+
// printf("secVtx (first) %f %f %f \n",secondaryVertex[0],secondaryVertex[1],secondaryVertex[2]);
42064207
covMatrixPCA = fgFitterTwoProngBarrel.calcPCACovMatrixFlat();
42074208
auto chi2PCA = fgFitterTwoProngBarrel.getChi2AtPCACandidate();
42084209
auto trackParVar0 = fgFitterTwoProngBarrel.getTrack(0);
@@ -4211,7 +4212,8 @@ void VarManager::FillPairVertexing(C const& collision, T const& t1, T const& t2,
42114212
v1 = {trackParVar0.getPt(), trackParVar0.getEta(), trackParVar0.getPhi(), m1};
42124213
v2 = {trackParVar1.getPt(), trackParVar1.getEta(), trackParVar1.getPhi(), m2};
42134214
v12 = v1 + v2;
4214-
if(fgPVrecalKF) primaryVertexNew = RecalculatePrimaryVertex(t1, t2, collision);
4215+
if (fgPVrecalKF)
4216+
primaryVertexNew = RecalculatePrimaryVertex(t1, t2, collision);
42154217

42164218
} else if constexpr (pairType == kDecayToMuMu && muonHasCov) {
42174219
// Get pca candidate from forward DCA fitter
@@ -4269,14 +4271,15 @@ void VarManager::FillPairVertexing(C const& collision, T const& t1, T const& t2,
42694271
values[kVertexingLxyProjected] = values[kVertexingLxyProjected] / TMath::Sqrt((v12.Px() * v12.Px()) + (v12.Py() * v12.Py()));
42704272
values[kVertexingLxyzProjected] = ((secondaryVertex[0] - collision.posX()) * v12.Px()) + ((secondaryVertex[1] - collision.posY()) * v12.Py()) + ((secondaryVertex[2] - collision.posZ()) * v12.Pz());
42714273
values[kVertexingLxyzProjected] = values[kVertexingLxyzProjected] / TMath::Sqrt((v12.Px() * v12.Px()) + (v12.Py() * v12.Py()) + (v12.Pz() * v12.Pz()));
4272-
if(fgPVrecalKF){
4273-
values[kVertexingLxyProjectedRecalculatePV] = (secondaryVertex[0] - primaryVertexNew.getX()) * v12.Px() + (secondaryVertex[1] - primaryVertexNew.getY()) * v12.Py();
4274-
values[kVertexingLxyProjectedRecalculatePV] = values[kVertexingLxyProjectedRecalculatePV] / v12.Pt();
4274+
if (fgPVrecalKF) {
4275+
values[kVertexingLxyProjectedRecalculatePV] = (secondaryVertex[0] - primaryVertexNew.getX()) * v12.Px() + (secondaryVertex[1] - primaryVertexNew.getY()) * v12.Py();
4276+
values[kVertexingLxyProjectedRecalculatePV] = values[kVertexingLxyProjectedRecalculatePV] / v12.Pt();
42754277
}
42764278
values[kVertexingTauxyProjected] = values[kVertexingLxyProjected] * v12.M() / (v12.Pt());
42774279
values[kVertexingTauxyProjectedPoleJPsiMass] = values[kVertexingLxyProjected] * o2::constants::physics::MassJPsi / (v12.Pt());
42784280
values[kVertexingTauxyProjectedNs] = values[kVertexingTauxyProjected] / o2::constants::physics::LightSpeedCm2NS;
4279-
if(fgPVrecalKF) values[kVertexingTauxyProjectedPoleJPsiMassRecalculatePV] = values[kVertexingLxyProjectedRecalculatePV] * o2::constants::physics::MassJPsi / (v12.Pt());
4281+
if (fgPVrecalKF)
4282+
values[kVertexingTauxyProjectedPoleJPsiMassRecalculatePV] = values[kVertexingLxyProjectedRecalculatePV] * o2::constants::physics::MassJPsi / (v12.Pt());
42804283
values[kVertexingTauzProjected] = values[kVertexingLzProjected] * v12.M() / TMath::Abs(v12.Pz());
42814284
values[kVertexingTauxyzProjected] = values[kVertexingLxyzProjected] * v12.M() / (v12.P());
42824285
}
@@ -4495,37 +4498,38 @@ void VarManager::FillPairVertexing(C const& collision, T const& t1, T const& t2,
44954498
}
44964499

44974500
template <int pairType, uint32_t collFillMap, uint32_t fillMap, typename C, typename T>
4498-
void VarManager::FillPairVertexingRecomputePV(C const& /*collision*/, T const& t1, T const& t2, o2::dataformats::VertexBase pvRefitted, float* values )
4499-
{
4500-
// recompute decay lenght variables using updated primary vertex
4501+
void VarManager::FillPairVertexingRecomputePV(C const& /*collision*/, T const& t1, T const& t2, o2::dataformats::VertexBase pvRefitted, float* values)
4502+
{
4503+
// recompute decay lenght variables using updated primary vertex
45014504

4502-
// check at compile time that the event and cov matrix have the cov matrix
4503-
constexpr bool eventHasVtxCov = ((collFillMap & Collision) > 0 || (collFillMap & ReducedEventVtxCov) > 0);
4504-
constexpr bool trackHasCov = ((fillMap & TrackCov) > 0 || (fillMap & ReducedTrackBarrelCov) > 0);
4505-
constexpr bool muonHasCov = ((fillMap & MuonCov) > 0 || (fillMap & ReducedMuonCov) > 0);
4505+
// check at compile time that the event and cov matrix have the cov matrix
4506+
constexpr bool eventHasVtxCov = ((collFillMap & Collision) > 0 || (collFillMap & ReducedEventVtxCov) > 0);
4507+
constexpr bool trackHasCov = ((fillMap & TrackCov) > 0 || (fillMap & ReducedTrackBarrelCov) > 0);
4508+
constexpr bool muonHasCov = ((fillMap & MuonCov) > 0 || (fillMap & ReducedMuonCov) > 0);
45064509

4507-
if (!values) {
4510+
if (!values) {
45084511
values = fgValues;
4509-
}
4512+
}
45104513

4511-
float m1 = o2::constants::physics::MassElectron;
4512-
float m2 = o2::constants::physics::MassElectron;
4513-
if constexpr (pairType == kDecayToKPi) {
4514+
float m1 = o2::constants::physics::MassElectron;
4515+
float m2 = o2::constants::physics::MassElectron;
4516+
if constexpr (pairType == kDecayToKPi) {
45144517
m1 = o2::constants::physics::MassKaonCharged;
45154518
m2 = o2::constants::physics::MassPionCharged;
4516-
}
4517-
if constexpr (pairType == kDecayToMuMu && muonHasCov) {
4519+
}
4520+
if constexpr (pairType == kDecayToMuMu && muonHasCov) {
45184521
m1 = o2::constants::physics::MassMuon;
45194522
m2 = o2::constants::physics::MassMuon;
4520-
}
4521-
ROOT::Math::PtEtaPhiMVector v1(t1.pt(), t1.eta(), t1.phi(), m1);
4522-
ROOT::Math::PtEtaPhiMVector v2(t2.pt(), t2.eta(), t2.phi(), m2);
4523-
ROOT::Math::PtEtaPhiMVector v12 = v1 + v2;
4523+
}
4524+
ROOT::Math::PtEtaPhiMVector v1(t1.pt(), t1.eta(), t1.phi(), m1);
4525+
ROOT::Math::PtEtaPhiMVector v2(t2.pt(), t2.eta(), t2.phi(), m2);
4526+
ROOT::Math::PtEtaPhiMVector v12 = v1 + v2;
45244527

4525-
if(fgFitterTwoProngBarrel.getNCandidates() == 0) return;
4526-
Vec3D secondaryVertex;
4528+
if (fgFitterTwoProngBarrel.getNCandidates() == 0)
4529+
return;
4530+
Vec3D secondaryVertex;
45274531

4528-
if (!fgUsedKF) { // to be updated when seconday vertex is computed with KF
4532+
if (!fgUsedKF) { // to be updated when seconday vertex is computed with KF
45294533
if constexpr (eventHasVtxCov) {
45304534

45314535
if constexpr ((pairType == kDecayToEE || pairType == kDecayToKPi) && trackHasCov) {
@@ -4554,7 +4558,6 @@ void VarManager::FillPairVertexingRecomputePV(C const& /*collision*/, T const& t
45544558
values[kVertexingTauxyProjectedPoleJPsiMassRecalculatePV] = values[kVertexingLxyProjectedRecalculatePV] * o2::constants::physics::MassJPsi / (v12.Pt());
45554559
}
45564560
}
4557-
45584561
}
45594562

45604563
template <uint32_t collFillMap, uint32_t fillMap, typename C, typename T>

‎PWGDQ/Tasks/dqEfficiency_withAssoc_direct.cxx‎

Lines changed: 62 additions & 60 deletions
Original file line numberDiff line numberDiff line change
@@ -23,9 +23,6 @@
2323
#include "PWGDQ/Core/MixingLibrary.h"
2424
#include "PWGDQ/Core/VarManager.h"
2525
#include "PWGDQ/DataModel/ReducedInfoTables.h"
26-
#include <DetectorsVertexing/PVertexer.h>
27-
#include <ReconstructionDataFormats/Track.h>
28-
#include <Common/Core/trackUtilities.h>
2926

3027
#include "Common/Core/PID/PIDTOFParamService.h"
3128
#include "Common/Core/TableHelper.h"
@@ -34,6 +31,7 @@
3431
#include "Common/DataModel/McCollisionExtra.h"
3532
#include "Common/DataModel/Multiplicity.h"
3633
#include "Common/DataModel/TrackSelectionTables.h"
34+
#include <Common/Core/trackUtilities.h>
3735

3836
#include "CCDB/BasicCCDBManager.h"
3937
#include "DataFormatsParameters/GRPMagField.h"
@@ -46,6 +44,8 @@
4644
#include "Framework/AnalysisHelpers.h"
4745
#include "Framework/AnalysisTask.h"
4846
#include "Framework/runDataProcessing.h"
47+
#include <DetectorsVertexing/PVertexer.h>
48+
#include <ReconstructionDataFormats/Track.h>
4949

5050
#include "TGeoGlobalMagField.h"
5151
#include <TH1F.h>
@@ -1163,7 +1163,6 @@ struct AnalysisSameEventPairing {
11631163
std::vector<o2::track::TrackParCov> pvContribTrackPars;
11641164
std::vector<bool> vec_useTrk_PVrefit;
11651165

1166-
11671166
// keep histogram class names in maps, so we don't have to buld their names in the pair loops
11681167
std::map<int, std::vector<TString>> fTrackHistNames;
11691168
std::map<int, std::vector<TString>> fBarrelHistNamesMCmatched;
@@ -1509,17 +1508,16 @@ struct AnalysisSameEventPairing {
15091508
cout << "AnalysisSameEventPairing::init() completed" << endl;
15101509
}
15111510

1512-
15131511
void initParamsFromCCDB(uint64_t timestamp, bool withTwoProngFitter = true)
15141512
{
15151513
cout << "AnalysisSameEventPairing::initParamsFromCCDB() called for timestamp " << timestamp << endl;
15161514
if (fConfigOptions.useRemoteField.value) {
15171515
o2::parameters::GRPMagField* grpmag = fCCDB->getForTimeStamp<o2::parameters::GRPMagField>(fConfigCCDB.grpMagPath, timestamp);
1518-
o2::base::MatLayerCylSet* lut = o2::base::MatLayerCylSet::rectifyPtrFromFile(fCCDB->get<o2::base::MatLayerCylSet>(fConfigCCDB.lutPath));
1516+
o2::base::MatLayerCylSet* lut = o2::base::MatLayerCylSet::rectifyPtrFromFile(fCCDB->get<o2::base::MatLayerCylSet>(fConfigCCDB.lutPath));
15191517
float magField = 0.0;
15201518
if (grpmag != nullptr) {
15211519
magField = grpmag->getNominalL3Field();
1522-
o2::base::Propagator::initFieldFromGRP(grpmag);
1520+
o2::base::Propagator::initFieldFromGRP(grpmag);
15231521
o2::base::Propagator::Instance()->setMatLUT(lut);
15241522
} else {
15251523
LOGF(fatal, "GRP object is not available in CCDB at timestamp=%llu", timestamp);
@@ -1549,53 +1547,56 @@ struct AnalysisSameEventPairing {
15491547
cout << "AnalysisSameEventPairing::initParamsFromCCDB() completed" << endl;
15501548
}
15511549

1552-
template <typename Events, typename TTracks, typename Tracks>
1553-
bool refitPVWithPVertexer(Events const& collision, TTracks const& tracks, Tracks const& t1, Tracks const& t2, o2::dataformats::VertexBase& pvRefitted)
1554-
{
1555-
// --- build PV contributor list ---
1556-
pvContribGlobIDs.clear();
1557-
pvContribTrackPars.clear();
1558-
// int nMyPVContrib = 0; int nMyPVContribOrig = 0;
1559-
for (auto const& trk : tracks) {
1560-
// check if it is PV contributor
1561-
if (!trk.isPVContributor()) continue;
1562-
// check if it contributes to the vtx of this collision
1563-
if (trk.collisionId() != collision.globalIndex()) continue;
1564-
// nMyPVContribOrig++;
1565-
// --- remove t1 and t2 if they are PV contributors ---
1566-
if (trk.globalIndex() == t1.globalIndex() || trk.globalIndex() == t2.globalIndex()) continue;
1567-
// add tracks and parameters to the list
1568-
pvContribGlobIDs.push_back(trk.globalIndex());
1569-
pvContribTrackPars.push_back(getTrackParCov(trk));
1570-
// nMyPVContrib++;
1571-
}
1550+
template <typename Events, typename TTracks, typename Tracks>
1551+
bool refitPVWithPVertexer(Events const& collision, TTracks const& tracks, Tracks const& t1, Tracks const& t2, o2::dataformats::VertexBase& pvRefitted)
1552+
{
1553+
// --- build PV contributor list ---
1554+
pvContribGlobIDs.clear();
1555+
pvContribTrackPars.clear();
1556+
// int nMyPVContrib = 0; int nMyPVContribOrig = 0;
1557+
for (auto const& trk : tracks) {
1558+
// check if it is PV contributor
1559+
if (!trk.isPVContributor())
1560+
continue;
1561+
// check if it contributes to the vtx of this collision
1562+
if (trk.collisionId() != collision.globalIndex())
1563+
continue;
1564+
// nMyPVContribOrig++;
1565+
// --- remove t1 and t2 if they are PV contributors ---
1566+
if (trk.globalIndex() == t1.globalIndex() || trk.globalIndex() == t2.globalIndex())
1567+
continue;
1568+
// add tracks and parameters to the list
1569+
pvContribGlobIDs.push_back(trk.globalIndex());
1570+
pvContribTrackPars.push_back(getTrackParCov(trk));
1571+
// nMyPVContrib++;
1572+
}
15721573

1573-
// cout << "contributors from collision: " << collision.numContrib() << " - from refitting: before -> " << nMyPVContribOrig << " after -> " << nMyPVContrib << endl;
1574-
vec_useTrk_PVrefit.assign(pvContribGlobIDs.size(), true);
1575-
// --- build VertexBase from event collision ---
1576-
o2::dataformats::VertexBase Pvtx;
1577-
Pvtx.setX(collision.posX());
1578-
Pvtx.setY(collision.posY());
1579-
Pvtx.setZ(collision.posZ());
1580-
Pvtx.setCov(collision.covXX(), collision.covXY(), collision.covYY(),
1581-
collision.covXZ(), collision.covYZ(), collision.covZZ());
1582-
1583-
// --- configure vertexer ---
1584-
o2::vertexing::PVertexer vertexer;
1585-
if (fConfigOptions.removeDiamondConstrPV) {
1586-
o2::conf::ConfigurableParam::updateFromString("pvertexer.useMeanVertexConstraint=false");
1587-
}
1588-
vertexer.init();
1574+
// cout << "contributors from collision: " << collision.numContrib() << " - from refitting: before -> " << nMyPVContribOrig << " after -> " << nMyPVContrib << endl;
1575+
vec_useTrk_PVrefit.assign(pvContribGlobIDs.size(), true);
1576+
// --- build VertexBase from event collision ---
1577+
o2::dataformats::VertexBase Pvtx;
1578+
Pvtx.setX(collision.posX());
1579+
Pvtx.setY(collision.posY());
1580+
Pvtx.setZ(collision.posZ());
1581+
Pvtx.setCov(collision.covXX(), collision.covXY(), collision.covYY(),
1582+
collision.covXZ(), collision.covYZ(), collision.covZZ());
15891583

1590-
bool PVrefit_doable = vertexer.prepareVertexRefit(pvContribTrackPars, Pvtx);
1591-
if (!PVrefit_doable) return false;
1584+
// --- configure vertexer ---
1585+
o2::vertexing::PVertexer vertexer;
1586+
if (fConfigOptions.removeDiamondConstrPV) {
1587+
o2::conf::ConfigurableParam::updateFromString("pvertexer.useMeanVertexConstraint=false");
1588+
}
1589+
vertexer.init();
15921590

1593-
// --- do the refit ---
1594-
pvRefitted = vertexer.refitVertex(vec_useTrk_PVrefit, Pvtx);
1591+
bool PVrefit_doable = vertexer.prepareVertexRefit(pvContribTrackPars, Pvtx);
1592+
if (!PVrefit_doable)
1593+
return false;
15951594

1596-
return true;
1597-
}
1595+
// --- do the refit ---
1596+
pvRefitted = vertexer.refitVertex(vec_useTrk_PVrefit, Pvtx);
15981597

1598+
return true;
1599+
}
15991600

16001601
// Template function to run same event pairing (barrel-barrel, muon-muon, barrel-muon)
16011602
template <bool TTwoProngFitter, int TPairType, uint32_t TEventFillMap, uint32_t TTrackFillMap, typename TEvents, typename TTracks>
@@ -1710,17 +1711,18 @@ bool refitPVWithPVertexer(Events const& collision, TTracks const& tracks, Track
17101711
}
17111712
if constexpr (TTwoProngFitter) {
17121713
VarManager::FillPairVertexing<TPairType, TEventFillMap, TTrackFillMap>(event, t1, t2, fConfigOptions.propToPCA);
1713-
if(fConfigOptions.recomputePV){
1714-
VarManager::SetPVrecalculationKF(false);
1715-
VarManager::ResetValues(VarManager::kVertexingLxyProjectedRecalculatePV, VarManager::kVertexingLxyProjectedRecalculatePV+1);
1716-
VarManager::ResetValues(VarManager::kVertexingTauxyProjectedPoleJPsiMassRecalculatePV, VarManager::kVertexingTauxyProjectedPoleJPsiMassRecalculatePV+1);
1717-
// cout << "primary vertex (before): x -> " << event.posX() << " y -> " << event.posY() << " z -> " << event.posZ() << endl;
1718-
o2::dataformats::VertexBase pvRefit;
1719-
bool ok = refitPVWithPVertexer(event, tracks, t1, t2, pvRefit);
1720-
if(ok) VarManager::FillPairVertexingRecomputePV<TPairType, TEventFillMap, TTrackFillMap>(event, t1, t2, pvRefit);
1721-
// cout << "primary vertex (after): ok -> " << ok << " x -> " << pvRefit.getX() << " y -> " << pvRefit.getY() << " z -> " << pvRefit.getZ() << endl;
1714+
if (fConfigOptions.recomputePV) {
1715+
VarManager::SetPVrecalculationKF(false);
1716+
VarManager::ResetValues(VarManager::kVertexingLxyProjectedRecalculatePV, VarManager::kVertexingLxyProjectedRecalculatePV + 1);
1717+
VarManager::ResetValues(VarManager::kVertexingTauxyProjectedPoleJPsiMassRecalculatePV, VarManager::kVertexingTauxyProjectedPoleJPsiMassRecalculatePV + 1);
1718+
// cout << "primary vertex (before): x -> " << event.posX() << " y -> " << event.posY() << " z -> " << event.posZ() << endl;
1719+
o2::dataformats::VertexBase pvRefit;
1720+
bool ok = refitPVWithPVertexer(event, tracks, t1, t2, pvRefit);
1721+
if (ok)
1722+
VarManager::FillPairVertexingRecomputePV<TPairType, TEventFillMap, TTrackFillMap>(event, t1, t2, pvRefit);
1723+
// cout << "primary vertex (after): ok -> " << ok << " x -> " << pvRefit.getX() << " y -> " << pvRefit.getY() << " z -> " << pvRefit.getZ() << endl;
17221724
}
1723-
}
1725+
}
17241726
if constexpr (eventHasQvector) {
17251727
VarManager::FillPairVn<TPairType>(t1, t2);
17261728
}
@@ -1858,7 +1860,7 @@ bool refitPVWithPVertexer(Events const& collision, TTracks const& tracks, Track
18581860
fHistMan->FillHistClass(histNames[icut][0].Data(), VarManager::fgValues); // reconstructed, unmatched
18591861
for (unsigned int isig = 0; isig < fRecMCSignals.size(); isig++) { // loop over MC signals
18601862
if (mcDecision & (static_cast<uint32_t>(1) << isig)) {
1861-
PromptNonPromptSepTable(VarManager::fgValues[VarManager::kMass], VarManager::fgValues[VarManager::kPt], VarManager::fgValues[VarManager::kEta], VarManager::fgValues[VarManager::kRap], VarManager::fgValues[VarManager::kPhi], VarManager::fgValues[VarManager::kVertexingTauxyProjected], VarManager::fgValues[VarManager::kVertexingTauxyProjectedPoleJPsiMass], VarManager::fgValues[VarManager::kVertexingTauzProjected], VarManager::fgValues[VarManager::kVertexingTauxyProjectedPoleJPsiMassRecalculatePV], isAmbiInBunch, isAmbiOutOfBunch, isCorrect_pair, VarManager::fgValues[VarManager::kMultFT0A], VarManager::fgValues[VarManager::kMultFT0C], VarManager::fgValues[VarManager::kCentFT0M], VarManager::fgValues[VarManager::kVtxNcontribReal]);
1863+
PromptNonPromptSepTable(VarManager::fgValues[VarManager::kMass], VarManager::fgValues[VarManager::kPt], VarManager::fgValues[VarManager::kEta], VarManager::fgValues[VarManager::kRap], VarManager::fgValues[VarManager::kPhi], VarManager::fgValues[VarManager::kVertexingTauxyProjected], VarManager::fgValues[VarManager::kVertexingTauxyProjectedPoleJPsiMass], VarManager::fgValues[VarManager::kVertexingTauzProjected], VarManager::fgValues[VarManager::kVertexingTauxyProjectedPoleJPsiMassRecalculatePV], isAmbiInBunch, isAmbiOutOfBunch, isCorrect_pair, VarManager::fgValues[VarManager::kMultFT0A], VarManager::fgValues[VarManager::kMultFT0C], VarManager::fgValues[VarManager::kCentFT0M], VarManager::fgValues[VarManager::kVtxNcontribReal]);
18621864
fHistMan->FillHistClass(histNamesMC[icut * fRecMCSignals.size() + isig][0].Data(), VarManager::fgValues); // matched signal
18631865
/*if (fConfigOptions.fConfigMiniTree) {
18641866
if constexpr (TPairType == VarManager::kDecayToMuMu) {

0 commit comments

Comments
 (0)