Skip to content

Commit 23a71b6

Browse files
gbencedigbencedi
committed
Added explicit track selection for reassociation
1 parent aedbd51 commit 23a71b6

File tree

1 file changed

+34
-35
lines changed

1 file changed

+34
-35
lines changed

‎PWGMM/Mult/Tasks/dndetaMFTPbPb.cxx‎

Lines changed: 34 additions & 35 deletions
Original file line numberDiff line numberDiff line change
@@ -106,7 +106,8 @@ struct DndetaMFTPbPb {
106106
Configurable<bool> cfgUseIRCut{"cfgUseIRCut", false, "Flag to cut on IR rate"};
107107
Configurable<bool> cfgIRCrashOnNull{"cfgIRCrashOnNull", false, "Flag to avoid CTP RateFetcher crash"};
108108
Configurable<std::string> cfgIRSource{"cfgIRSource", "ZNC hadronic", "Estimator of the interaction rate (Pb-Pb: ZNC hadronic)"};
109-
Configurable<bool> cfgUseTrackSel{"cfgUseTrackSel", false, "Flag to apply track selection"}; Configurable<bool> cfgUseParticleSel{"cfgUseParticleSel", false, "Flag to apply particle selection"};
109+
Configurable<bool> cfgUseTrackSel{"cfgUseTrackSel", false, "Flag to apply track selection"};
110+
Configurable<bool> cfgUseParticleSel{"cfgUseParticleSel", false, "Flag to apply particle selection"};
110111

111112
struct : ConfigurableGroup {
112113
ConfigurableAxis interactionRateBins{"interactionRateBins", {500, 0, 50}, "Binning for the interaction rate (kHz)"};
@@ -819,10 +820,10 @@ struct DndetaMFTPbPb {
819820
qaregistry.add({"Tracks/hDCAxyBestGenSecMat", "; DCA_{xy}^{gen} (cm); counts", {HistType::kTH1F, {dcaxyAxis}}});
820821
qaregistry.add({"Tracks/hDCAzBestGenSecMat", "; DCA_{z}^{gen} (cm); counts", {HistType::kTH1F, {dcazAxis}}});
821822
//
822-
qaregistry.add({"Tracks/hDCAxyBestPtRec","; p_{T} (GeV/c) ; DCA_{XY} (cm)", {HistType::kTH2F, {ptAxis, dcaxyAxis}}});
823-
qaregistry.add({"Tracks/hDCAzBestPtRec","; p_{T} (GeV/c) ; DCA_{Z} (cm)", {HistType::kTH2F, {ptAxis, dcazAxis}}});
824-
qaregistry.add({"Tracks/hDCAxyBestRec","; DCA_{XY} (cm)", {HistType::kTH1F, {dcaxyAxis}}});
825-
qaregistry.add({"Tracks/hDCAzBestRec","; DCA_{Z} (cm)", {HistType::kTH1F, {dcazAxis}}});
823+
qaregistry.add({"Tracks/hDCAxyBestPtRec", "; p_{T} (GeV/c) ; DCA_{XY} (cm)", {HistType::kTH2F, {ptAxis, dcaxyAxis}}});
824+
qaregistry.add({"Tracks/hDCAzBestPtRec", "; p_{T} (GeV/c) ; DCA_{Z} (cm)", {HistType::kTH2F, {ptAxis, dcazAxis}}});
825+
qaregistry.add({"Tracks/hDCAxyBestRec", "; DCA_{XY} (cm)", {HistType::kTH1F, {dcaxyAxis}}});
826+
qaregistry.add({"Tracks/hDCAzBestRec", "; DCA_{Z} (cm)", {HistType::kTH1F, {dcazAxis}}});
826827
}
827828
if (doprocessDCAMcCentFT0C) {
828829
qaregistry.add({"Tracks/Centrality/hDCAxBestGenPrim", "; DCA_{x}^{gen} (cm); counts", {HistType::kTH1F, {dcaxyAxis}}});
@@ -840,10 +841,10 @@ struct DndetaMFTPbPb {
840841
qaregistry.add({"Tracks/Centrality/hDCAxyBestGenSecMat", "; DCA_{xy}^{gen} (cm); counts", {HistType::kTH1F, {dcaxyAxis}}});
841842
qaregistry.add({"Tracks/Centrality/hDCAzBestGenSecMat", "; DCA_{z}^{gen} (cm); counts", {HistType::kTH1F, {dcazAxis}}});
842843
//
843-
qaregistry.add({"Tracks/Centrality/hDCAxyBestPtRec","; p_{T} (GeV/c) ; DCA_{XY} (cm); centrality", {HistType::kTHnSparseF, {ptAxis, dcaxyAxis, centralityAxis}}});
844-
qaregistry.add({"Tracks/Centrality/hDCAzBestPtRec","; p_{T} (GeV/c) ; DCA_{Z} (cm); centrality", {HistType::kTHnSparseF, {ptAxis, dcazAxis, centralityAxis}}});
845-
qaregistry.add({"Tracks/Centrality/hDCAxyBestRec","; DCA_{XY} (cm); centrality", {HistType::kTH2F, {dcaxyAxis, centralityAxis}}});
846-
qaregistry.add({"Tracks/Centrality/hDCAzBestRec","; DCA_{Z} (cm); centrality", {HistType::kTH2F, {dcazAxis, centralityAxis}}});
844+
qaregistry.add({"Tracks/Centrality/hDCAxyBestPtRec", "; p_{T} (GeV/c) ; DCA_{XY} (cm); centrality", {HistType::kTHnSparseF, {ptAxis, dcaxyAxis, centralityAxis}}});
845+
qaregistry.add({"Tracks/Centrality/hDCAzBestPtRec", "; p_{T} (GeV/c) ; DCA_{Z} (cm); centrality", {HistType::kTHnSparseF, {ptAxis, dcazAxis, centralityAxis}}});
846+
qaregistry.add({"Tracks/Centrality/hDCAxyBestRec", "; DCA_{XY} (cm); centrality", {HistType::kTH2F, {dcaxyAxis, centralityAxis}}});
847+
qaregistry.add({"Tracks/Centrality/hDCAzBestRec", "; DCA_{Z} (cm); centrality", {HistType::kTH2F, {dcazAxis, centralityAxis}}});
847848
}
848849
}
849850

@@ -860,6 +861,7 @@ struct DndetaMFTPbPb {
860861
/// Joined tables
861862
using FullBCs = soa::Join<aod::BCsWithTimestamps, aod::BcSels>;
862863
using CollBCs = soa::Join<aod::BCsWithTimestamps, aod::Run3MatchedToBCSparse>;
864+
// Collisions
863865
using Colls = soa::Join<aod::Collisions, aod::EvSels>;
864866
using Coll = Colls::iterator;
865867
using CollsCentFT0C = soa::Join<aod::Collisions, aod::CentFT0Cs, aod::EvSels>;
@@ -874,7 +876,7 @@ struct DndetaMFTPbPb {
874876
aod::CentFT0Cs, aod::EvSels>;
875877
using CollGenCent = CollsGenCentFT0C::iterator;
876878
using CollsCorr = soa::Join<aod::Collisions, aod::EvSels, aod::Mults, aod::PVMults, aod::CentFT0Cs, aod::CentFV0As, aod::CentFT0CVariant1s, aod::CentFT0Ms, aod::CentNGlobals, aod::CentMFTs>;
877-
879+
// Tracks
878880
using MFTTracksLabeled = soa::Join<aod::MFTTracks, aod::McMFTTrackLabels>;
879881
using MftTracksWColls = soa::Join<aod::MFTTracks, aod::MFTTrkCompColls>;
880882

@@ -2201,12 +2203,12 @@ struct DndetaMFTPbPb {
22012203
if (itrack.eta() > trackCuts.minEta && itrack.eta() < trackCuts.maxEta) {
22022204
if constexpr (has_reco_cent<C>) {
22032205
qaregistry.fill(HIST("Tracks/Centrality/hPtPhiEtaZvtxEffBestRec"),
2204-
particle.pt(), particle.phi(), particle.eta(),
2205-
mcCollision.posZ(), crec, occrec);
2206+
particle.pt(), particle.phi(), particle.eta(),
2207+
mcCollision.posZ(), crec, occrec);
22062208
} else {
22072209
qaregistry.fill(HIST("Tracks/hPtPhiEtaZvtxEffBestRec"), particle.pt(),
2208-
particle.phi(), particle.eta(), mcCollision.posZ(),
2209-
occrec);
2210+
particle.phi(), particle.eta(), mcCollision.posZ(),
2211+
occrec);
22102212
}
22112213
}
22122214
} else {
@@ -2233,12 +2235,12 @@ struct DndetaMFTPbPb {
22332235
if (track.eta() > trackCuts.minEta && track.eta() < trackCuts.maxEta) {
22342236
if constexpr (has_reco_cent<C>) {
22352237
qaregistry.fill(HIST("Tracks/Centrality/hPtPhiEtaZvtxEffBestRec"),
2236-
particle.pt(), particle.phi(), particle.eta(),
2237-
mcCollision.posZ(), crec, occrec);
2238+
particle.pt(), particle.phi(), particle.eta(),
2239+
mcCollision.posZ(), crec, occrec);
22382240
} else {
22392241
qaregistry.fill(HIST("Tracks/hPtPhiEtaZvtxEffBestRec"), particle.pt(),
2240-
particle.phi(), particle.eta(), mcCollision.posZ(),
2241-
occrec);
2242+
particle.phi(), particle.eta(), mcCollision.posZ(),
2243+
occrec);
22422244
}
22432245
}
22442246
} else {
@@ -2379,7 +2381,7 @@ struct DndetaMFTPbPb {
23792381
if constexpr (has_reco_cent<C>) {
23802382
qaregistry.fill(HIST("Tracks/Centrality/hMftTracksAmbDegreeWoTrivial"), track.compatibleColl().size(), c);
23812383
} else {
2382-
qaregistry.fill(HIST("Tracks/hMftTracksAmbDegreeWoTrivial"), track.compatibleColl().size());
2384+
qaregistry.fill(HIST("Tracks/hMftTracksAmbDegreeWoTrivial"), track.compatibleColl().size());
23832385
}
23842386
continue;
23852387
}
@@ -2457,10 +2459,10 @@ struct DndetaMFTPbPb {
24572459
/// @brief process template function for MC QA checks
24582460
template <typename C>
24592461
void processMcQA(typename soa::Join<C, aod::McCollisionLabels> const& collisions,
2460-
MFTTracksLabeled const& tracks,
2461-
aod::AmbiguousMFTTracks const& atracks,
2462-
aod::McCollisions const& mcCollisions,
2463-
FiltParticles const& /*particles*/)
2462+
MFTTracksLabeled const& tracks,
2463+
aod::AmbiguousMFTTracks const& atracks,
2464+
aod::McCollisions const& mcCollisions,
2465+
FiltParticles const& /*particles*/)
24642466
{
24652467
for (const auto& collision : collisions) {
24662468
float crec = getRecoCent(collision);
@@ -2529,15 +2531,14 @@ struct DndetaMFTPbPb {
25292531
PROCESS_SWITCH(DndetaMFTPbPb, processMcQACentFT0C,
25302532
"Process MC QA checks (in FT0 centrality bins)", false);
25312533

2532-
25332534
/// @brief process template function for DCA MC checks
25342535
template <typename C, typename MC>
25352536
void processDCAMc(typename soa::Join<C, aod::McCollisionLabels>::iterator const& collision,
25362537
FiltMcMftTracks const& /*tracks*/,
25372538
soa::SmallGroups<aod::BestCollisionsFwd3d> const& besttracks,
25382539
MC const& /*mcCollisions*/,
25392540
aod::McParticles const& /*particles*/
2540-
)
2541+
)
25412542
{
25422543
if (!isGoodEvent<false>(collision)) {
25432544
return;
@@ -2596,8 +2597,8 @@ struct DndetaMFTPbPb {
25962597
const auto dcaZtruth(particle.vz() - particle.mcCollision().posZ());
25972598
auto dcaXYtruth = std::sqrt(dcaXtruth * dcaXtruth + dcaYtruth * dcaYtruth);
25982599

2599-
if (!particle.isPhysicalPrimary()) { // Secondaries (weak decays and material)
2600-
if (particle.getProcess() == kNumDecay) { // Particles from decay
2600+
if (!particle.isPhysicalPrimary()) { // Secondaries (weak decays and material)
2601+
if (particle.getProcess() == kNumDecay) { // Particles from decay
26012602
if constexpr (has_reco_cent<C>) {
26022603
qaregistry.fill(HIST("Tracks/Centrality/hDCAxBestGenSecWeak"), dcaXtruth, crec);
26032604
qaregistry.fill(HIST("Tracks/Centrality/hDCAyBestGenSecWeak"), dcaYtruth, crec);
@@ -2622,7 +2623,7 @@ struct DndetaMFTPbPb {
26222623
qaregistry.fill(HIST("Tracks/hDCAzBestGenSecMat"), dcaZtruth);
26232624
}
26242625
}
2625-
} else { // Primaries
2626+
} else { // Primaries
26262627
if constexpr (has_reco_cent<C>) {
26272628
qaregistry.fill(HIST("Tracks/Centrality/hDCAxBestGenPrim"), dcaXtruth, crec);
26282629
qaregistry.fill(HIST("Tracks/Centrality/hDCAyBestGenPrim"), dcaYtruth, crec);
@@ -2639,11 +2640,10 @@ struct DndetaMFTPbPb {
26392640
}
26402641

26412642
void processDCAMcInlcusive(soa::Join<Colls, aod::McCollisionLabels>::iterator const& collision,
2642-
FiltMcMftTracks const& tracks,
2643-
soa::SmallGroups<aod::BestCollisionsFwd3d> const& besttracks,
2644-
aod::McCollisions const& mccollisions,
2645-
aod::McParticles const& particles
2646-
)
2643+
FiltMcMftTracks const& tracks,
2644+
soa::SmallGroups<aod::BestCollisionsFwd3d> const& besttracks,
2645+
aod::McCollisions const& mccollisions,
2646+
aod::McParticles const& particles)
26472647
{
26482648
processDCAMc<Colls, aod::McCollisions>(collision, tracks, besttracks, mccollisions, particles);
26492649
}
@@ -2654,8 +2654,7 @@ struct DndetaMFTPbPb {
26542654
FiltMcMftTracks const& tracks,
26552655
soa::SmallGroups<aod::BestCollisionsFwd3d> const& besttracks,
26562656
aod::McCollisions const& mccollisions,
2657-
aod::McParticles const& particles
2658-
)
2657+
aod::McParticles const& particles)
26592658
{
26602659
processDCAMc<CollsCentFT0C, aod::McCollisions>(collision, tracks, besttracks, mccollisions, particles);
26612660
}

0 commit comments

Comments
 (0)