Skip to content

[PWGCF - FemtoDream]: Including MC collision table #5799

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
victor-gonzalez merged 5 commits into from
Apr 23, 2024
Merged

[PWGCF - FemtoDream]: Including MC collision table #5799

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
7071259
[PWGCF - FemtoDream]: Including MC collision table
Apr 13, 2024
f26074a
clang format
Apr 15, 2024
a03e60c
Fixing MC configurable for track-V0 task
Apr 22, 2024
dfe6db2
Fix the name of the newly added histogram
Apr 23, 2024
224cf74
clang format
Apr 23, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
16 changes: 16 additions & 0 deletions PWGCF/DataModel/FemtoDerived.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -66,6 +66,22 @@ DECLARE_SOA_TABLE(FDColMasks, "AOD", "FDCOLMASK",
DECLARE_SOA_TABLE(FDDownSample, "AOD", "FDDOWNSAMPLE",
femtodreamcollision::Downsample);

namespace femtodreamMCcollision
{
DECLARE_SOA_COLUMN(MultMCgenPartEta08, multMCgenPartEta08, int); //! Multiplicity of the event as given by the generator in |eta|<0.8
}

DECLARE_SOA_TABLE(FDMCCollisions, "AOD", "FDMCCOLLISION",
o2::soa::Index<>,
femtodreamMCcollision::MultMCgenPartEta08);
using FDMCCollision = FDMCCollisions::iterator;

namespace mcfdcolllabel
{
DECLARE_SOA_INDEX_COLUMN(FDMCCollision, fdMCCollision); //! MC collision for femtodreamcollision
}
DECLARE_SOA_TABLE(FDMCCollLabels, "AOD", "FDMCCollLabel", mcfdcolllabel::FDMCCollisionId);

/// FemtoDreamTrack
namespace femtodreamparticle
{
Expand Down
14 changes: 12 additions & 2 deletions PWGCF/FemtoDream/Core/femtoDreamEventHisto.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -31,7 +31,7 @@ class FemtoDreamEventHisto
virtual ~FemtoDreamEventHisto() = default;
/// Initializes histograms for the task
/// \param registry Histogram registry to be passed
void init(HistogramRegistry* registry)
void init(HistogramRegistry* registry, bool isMC)
{
mHistogramRegistry = registry;
mHistogramRegistry->add("Event/hZvtx", "; vtx_{z} (cm); Entries", kTH1F, {{300, -12.5, 12.5}});
Expand All @@ -40,12 +40,16 @@ class FemtoDreamEventHisto
mHistogramRegistry->add("Event/hMultNTrVsZvtx", "; Multiplicity (MultNtr); vtx_{z} (cm)", kTH2F, {{200, 0, 200}, {300, -12.5, 12.5}});
mHistogramRegistry->add("Event/hMultNTrVsMultPercentile", "; Multiplicity (MultNtr); Multiplicity Percentile (FT0M)", kTH2F, {{200, 0, 200}, {110, 0, 110}});
mHistogramRegistry->add("Event/hMultPercentileVsZvtx", "; Multiplicity Percentile (FT0M); vtx_{z} (cm)", kTH2F, {{110, 0, 110}, {300, -12.5, 12.5}});

if (isMC) {
mHistogramRegistry->add("Event_MC/hGenMult08VsMultPercentile", "; generated MC multiplicity (#eta<0.8); Multiplicity Percentile (FT0M)", kTH2F, {{200, 0, 200}, {110, 0, 110}});
}
}

/// Some basic QA of the event
/// \tparam T type of the collision
/// \param col Collision
template <typename T>
template <bool isMC, typename T>
void fillQA(T const& col)
{
if (mHistogramRegistry) {
Expand All @@ -55,6 +59,12 @@ class FemtoDreamEventHisto
mHistogramRegistry->fill(HIST("Event/hMultNTrVsZvtx"), col.multNtr(), col.posZ());
mHistogramRegistry->fill(HIST("Event/hMultNTrVsMultPercentile"), col.multNtr(), col.multV0M());
mHistogramRegistry->fill(HIST("Event/hMultPercentileVsZvtx"), col.multV0M(), col.posZ());

if constexpr (isMC) {
if (col.has_fdMCCollision()) {
mHistogramRegistry->fill(HIST("Event_MC/hGenMult08VsMultPercentile"), col.fdMCCollision().multMCgenPartEta08(), col.multV0M());
}
}
}
}

Expand Down
25 changes: 21 additions & 4 deletions PWGCF/FemtoDream/TableProducer/femtoDreamProducerTask.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -45,10 +45,11 @@ using namespace o2::analysis::femtoDream;
namespace o2::aod
{

using FemtoFullCollision =
soa::Join<aod::Collisions, aod::EvSels, aod::Mults, aod::CentFT0Ms>::iterator;
using FemtoFullCollision = soa::Join<aod::Collisions, aod::EvSels, aod::Mults, aod::CentFT0Ms>::iterator;
using FemtoFullCollisionMC = soa::Join<aod::Collisions, aod::EvSels, aod::Mults, aod::CentFT0Ms, aod::McCollisionLabels>::iterator;
using FemtoFullCollision_noCent_MC = soa::Join<aod::Collisions, aod::EvSels, aod::Mults, aod::McCollisionLabels>::iterator;
using FemtoFullMCgenCollisions = soa::Join<aod::McCollisions, MultsExtraMC>;
using FemtoFullMCgenCollision = FemtoFullMCgenCollisions::iterator;

using FemtoFullTracks =
soa::Join<aod::FullTracks, aod::TracksDCA,
Expand All @@ -72,6 +73,8 @@ int getRowDaughters(int daughID, T const& vecID)
struct femtoDreamProducerTask {

Produces<aod::FDCollisions> outputCollision;
Produces<aod::FDMCCollisions> outputMCCollision;
Produces<aod::FDMCCollLabels> outputCollsMCLabels;
Produces<aod::FDParticles> outputParts;
Produces<aod::FDMCParticles> outputPartsMC;
Produces<aod::FDExtParticles> outputDebugParts;
Expand Down Expand Up @@ -355,6 +358,17 @@ struct femtoDreamProducerTask {
}
}

template <typename CollisionType>
void fillMCCollision(CollisionType const& col)
{
if (col.has_mcCollision()) {
auto genMCcol = col.template mcCollision_as<aod::FemtoFullMCgenCollisions>();
outputMCCollision(genMCcol.multMCNParticlesEta08());
outputCollsMCLabels(outputMCCollision.lastIndex());
} else {
outputCollsMCLabels(-1);
}
}
template <bool isMC, bool useCentrality, typename V0Type, typename TrackType, typename CollisionType>
void fillCollisionsAndTracksAndV0(CollisionType const& col, TrackType const& tracks, V0Type const& fullV0s)
{
Expand Down Expand Up @@ -392,6 +406,9 @@ struct femtoDreamProducerTask {
}

outputCollision(vtxZ, mult, multNtr, spher, mMagField);
if constexpr (isMC) {
fillMCCollision(col);
}

std::vector<int> childIDs = {0, 0}; // these IDs are necessary to keep track of the children
std::vector<int> tmpIDtrack; // this vector keeps track of the matching of the primary track table row <-> aod::track table global index
Expand Down Expand Up @@ -530,7 +547,7 @@ struct femtoDreamProducerTask {
void processMC(aod::FemtoFullCollisionMC const& col,
aod::BCsWithTimestamps const&,
soa::Join<aod::FemtoFullTracks, aod::McTrackLabels> const& tracks,
aod::McCollisions const&,
aod::FemtoFullMCgenCollisions const&,
aod::McParticles const&,
soa::Join<o2::aod::V0Datas, aod::McV0Labels> const& fullV0s) /// \todo with FilteredFullV0s
{
Expand All @@ -544,7 +561,7 @@ struct femtoDreamProducerTask {
void processMC_noCentrality(aod::FemtoFullCollision_noCent_MC const& col,
aod::BCsWithTimestamps const&,
soa::Join<aod::FemtoFullTracks, aod::McTrackLabels> const& tracks,
aod::McCollisions const&,
aod::FemtoFullMCgenCollisions const&,
aod::McParticles const&,
soa::Join<o2::aod::V0Datas, aod::McV0Labels> const& fullV0s) /// \todo with FilteredFullV0s
{
Expand Down
13 changes: 8 additions & 5 deletions PWGCF/FemtoDream/Tasks/femtoDreamDebugTrack.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -62,6 +62,9 @@ struct femtoDreamDebugTrack {
Configurable<int> ConfTempFitVarMomentum{"ConfTempFitVarMomentum", 0, "Momentum used for binning: 0 -> pt; 1 -> preco; 2 -> ptpc"};
ConfigurableAxis ConfDummy{"ConfDummy", {1, 0, 1}, "Dummy axis for inv mass"};

using FemtoMCCollisions = Join<aod::FDCollisions, aod::FDMCCollLabels>;
using FemtoMCCollision = FemtoMCCollisions::iterator;

using FemtoFullParticles = soa::Join<aod::FDParticles, aod::FDExtParticles>;

Partition<FemtoFullParticles> partsOne = (aod::femtodreamparticle::partType == uint8_t(aod::femtodreamparticle::ParticleType::kTrack)) &&
Expand Down Expand Up @@ -90,15 +93,15 @@ struct femtoDreamDebugTrack {

void init(InitContext&)
{
eventHisto.init(&qaRegistry);
eventHisto.init(&qaRegistry, ConfIsMC);
trackHisto.init(&qaRegistry, ConfBinmult, ConfBinmultPercentile, ConfBinpT, ConfBineta, ConfBinphi, ConfTempFitVarBins, ConfNsigmaTPCBins, ConfNsigmaTOFBins, ConfNsigmaTPCTOFBins, ConfTPCclustersBins, ConfDummy, ConfIsMC, ConfTrk1_PDGCode.value, true, ConfOptCorrelatedPlots);
}

/// Porduce QA plots for sigle track selection in FemtoDream framework
template <bool isMC, typename PartitionType>
void FillDebugHistos(o2::aod::FDCollision& col, PartitionType& groupPartsOne)
template <bool isMC, typename CollisionType, typename PartitionType>
void FillDebugHistos(CollisionType& col, PartitionType& groupPartsOne)
{
eventHisto.fillQA(col);
eventHisto.fillQA<isMC>(col);
for (auto& part : groupPartsOne) {
trackHisto.fillQA<isMC, true>(part, static_cast<aod::femtodreamparticle::MomentumType>(ConfTempFitVarMomentum.value), col.multNtr(), col.multV0M(), ConfOptCorrelatedPlots);
}
Expand All @@ -119,7 +122,7 @@ struct femtoDreamDebugTrack {
/// \param col subscribe to FemtoDreamCollision table
/// \param parts subscribe to the joined table of FemtoDreamParticles and FemtoDreamMCLabels table
/// \param FemtoDramMCParticles subscribe to the table containing the Monte Carlo Truth information
void processMC(o2::aod::FDCollision& col, FemtoFullParticlesMC&, aod::FDMCParticles&, aod::FDExtMCParticles&)
void processMC(FemtoMCCollision& col, o2::aod::FDMCCollisions&, FemtoFullParticlesMC& parts, aod::FDMCParticles&, aod::FDExtMCParticles&)
{
auto groupPartsOne = partsOneMC->sliceByCached(aod::femtodreamparticle::fdCollisionId, col.globalIndex(), cache);
FillDebugHistos<true>(col, groupPartsOne);
Expand Down
4 changes: 2 additions & 2 deletions PWGCF/FemtoDream/Tasks/femtoDreamDebugV0.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -79,7 +79,7 @@ struct femtoDreamDebugV0 {

void init(InitContext&)
{
eventHisto.init(&EventRegistry);
eventHisto.init(&EventRegistry, false);
posChildHistos.init(&V0Registry, ConfBinmult, ConfDummy, ConfV0ChildTempFitVarMomentumBins, ConfDummy, ConfDummy, ConfChildTempFitVarBins, ConfV0ChildNsigmaTPCBins, ConfV0ChildNsigmaTOFBins, ConfV0ChildNsigmaTPCTOFBins, ConfDummy, ConfV0InvMassBins, false, ConfV01_ChildPos_PDGCode.value, true);
negChildHistos.init(&V0Registry, ConfBinmult, ConfDummy, ConfV0ChildTempFitVarMomentumBins, ConfDummy, ConfDummy, ConfChildTempFitVarBins, ConfV0ChildNsigmaTPCBins, ConfV0ChildNsigmaTOFBins, ConfV0ChildNsigmaTPCTOFBins, ConfDummy, ConfV0InvMassBins, false, ConfV01_ChildNeg_PDGCode, true);
V0Histos.init(&V0Registry, ConfBinmult, ConfDummy, ConfV0TempFitVarMomentumBins, ConfDummy, ConfDummy, ConfV0TempFitVarBins, ConfV0ChildNsigmaTPCBins, ConfV0ChildNsigmaTOFBins, ConfV0ChildNsigmaTPCTOFBins, ConfDummy, ConfV0InvMassBins, false, ConfV01_PDGCode.value, true);
Expand All @@ -89,7 +89,7 @@ struct femtoDreamDebugV0 {
void process(o2::aod::FDCollision const& col, FemtoFullParticles const& parts)
{
auto groupPartsOne = partsOne->sliceByCached(aod::femtodreamparticle::fdCollisionId, col.globalIndex(), cache);
eventHisto.fillQA(col);
eventHisto.fillQA<false>(col);
for (auto& part : groupPartsOne) {
if (!part.has_children()) {
continue;
Expand Down
30 changes: 19 additions & 11 deletions PWGCF/FemtoDream/Tasks/femtoDreamPairTaskTrackTrack.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -79,8 +79,14 @@ struct femtoDreamPairTaskTrackTrack {

using FilteredCollisions = soa::Filtered<FDCollisions>;
using FilteredCollision = FilteredCollisions::iterator;
using FilteredMCCollisions = soa::Filtered<soa::Join<aod::FDCollisions, aod::FDMCCollLabels>>;
using FilteredMCCollision = FilteredMCCollisions::iterator;

using FilteredMaskedCollisions = soa::Filtered<soa::Join<FDCollisions, FDColMasks, FDDownSample>>;
using FilteredMaskedCollision = FilteredMaskedCollisions::iterator;
using FilteredMaskedMCCollisions = soa::Filtered<soa::Join<FDCollisions, aod::FDMCCollLabels, FDColMasks, FDDownSample>>;
using FilteredMaskedMCCollision = FilteredMaskedMCCollisions::iterator;

femtodreamcollision::BitMaskType BitMask = -1;

/// Track 1
Expand Down Expand Up @@ -216,7 +222,7 @@ struct femtoDreamPairTaskTrackTrack {
if (Option.RandomizePair.value) {
random = new TRandom3(0);
}
eventHisto.init(&qaRegistry);
eventHisto.init(&qaRegistry, Option.IsMC);
trackHistoPartOne.init(&qaRegistry, Binning.multTempFit, Option.Dummy, Binning.TrackpT, Option.Dummy, Option.Dummy, Binning.TempFitVar, Option.Dummy, Option.Dummy, Option.Dummy, Option.Dummy, Option.Dummy, Option.IsMC, Track1.PDGCode);
if (!Option.SameSpecies) {
trackHistoPartTwo.init(&qaRegistry, Binning.multTempFit, Option.Dummy, Binning.TrackpT, Option.Dummy, Option.Dummy, Binning.TempFitVar, Option.Dummy, Option.Dummy, Option.Dummy, Option.Dummy, Option.Dummy, Option.IsMC, Track2.PDGCode);
Expand Down Expand Up @@ -286,10 +292,10 @@ struct femtoDreamPairTaskTrackTrack {
}
};

template <typename CollisionType>
template <bool isMC, typename CollisionType>
void fillCollision(CollisionType col)
{
eventHisto.fillQA(col);
eventHisto.fillQA<isMC>(col);
}

/// This function processes the same event and takes care of all the histogramming
Expand Down Expand Up @@ -358,7 +364,7 @@ struct femtoDreamPairTaskTrackTrack {
/// \param parts subscribe to the femtoDreamParticleTable
void processSameEvent(FilteredCollision& col, o2::aod::FDParticles& parts)
{
fillCollision(col);
fillCollision<false>(col);
auto SliceTrk1 = PartitionTrk1->sliceByCached(aod::femtodreamparticle::fdCollisionId, col.globalIndex(), cache);
auto SliceTrk2 = PartitionTrk2->sliceByCached(aod::femtodreamparticle::fdCollisionId, col.globalIndex(), cache);
if (SliceTrk1.size() == 0 && SliceTrk2.size() == 0) {
Expand All @@ -379,7 +385,7 @@ struct femtoDreamPairTaskTrackTrack {
return;
}
}
fillCollision(col);
fillCollision<false>(col);
auto SliceTrk1 = PartitionTrk1->sliceByCached(aod::femtodreamparticle::fdCollisionId, col.globalIndex(), cache);
auto SliceTrk2 = PartitionTrk2->sliceByCached(aod::femtodreamparticle::fdCollisionId, col.globalIndex(), cache);
doSameEvent<false>(SliceTrk1, SliceTrk2, parts, col);
Expand All @@ -390,10 +396,12 @@ struct femtoDreamPairTaskTrackTrack {
/// \param col subscribe to the collision table (Monte Carlo Reconstructed reconstructed)
/// \param parts subscribe to joined table FemtoDreamParticles and FemtoDreamMCLables to access Monte Carlo truth
/// \param FemtoDreamMCParticles subscribe to the Monte Carlo truth table
void processSameEventMC(FilteredCollision& col, soa::Join<o2::aod::FDParticles, o2::aod::FDMCLabels>& parts,
void processSameEventMC(FilteredMCCollision& col,
o2::aod::FDMCCollisions&,
soa::Join<o2::aod::FDParticles, o2::aod::FDMCLabels>& parts,
o2::aod::FDMCParticles&)
{
fillCollision(col);
fillCollision<true>(col);
auto SliceMCTrk1 = PartitionMCTrk1->sliceByCached(aod::femtodreamparticle::fdCollisionId, col.globalIndex(), cache);
auto SliceMCTrk2 = PartitionMCTrk2->sliceByCached(aod::femtodreamparticle::fdCollisionId, col.globalIndex(), cache);
if (SliceMCTrk1.size() == 0 && SliceMCTrk2.size() == 0) {
Expand All @@ -403,7 +411,7 @@ struct femtoDreamPairTaskTrackTrack {
}
PROCESS_SWITCH(femtoDreamPairTaskTrackTrack, processSameEventMC, "Enable processing same event for Monte Carlo", false);

void processSameEventMCMasked(FilteredMaskedCollision& col, soa::Join<o2::aod::FDParticles, o2::aod::FDMCLabels>& parts,
void processSameEventMCMasked(FilteredMaskedMCCollision& col, o2::aod::FDMCCollisions&, soa::Join<o2::aod::FDParticles, o2::aod::FDMCLabels>& parts,
o2::aod::FDMCParticles&)
{
if (Option.SameSpecies.value) {
Expand All @@ -415,7 +423,7 @@ struct femtoDreamPairTaskTrackTrack {
return;
}
}
fillCollision(col);
fillCollision<true>(col);
auto SliceMCTrk1 = PartitionMCTrk1->sliceByCached(aod::femtodreamparticle::fdCollisionId, col.globalIndex(), cache);
auto SliceMCTrk2 = PartitionMCTrk2->sliceByCached(aod::femtodreamparticle::fdCollisionId, col.globalIndex(), cache);
doSameEvent<true>(SliceMCTrk1, SliceMCTrk2, parts, col);
Expand Down Expand Up @@ -539,7 +547,7 @@ struct femtoDreamPairTaskTrackTrack {
/// @param cols subscribe to the collisions table (Monte Carlo Reconstructed reconstructed)
/// @param parts subscribe to joined table FemtoDreamParticles and FemtoDreamMCLables to access Monte Carlo truth
/// @param FemtoDreamMCParticles subscribe to the Monte Carlo truth table
void processMixedEventMC(FilteredCollisions& cols, soa::Join<o2::aod::FDParticles, o2::aod::FDMCLabels>& parts, o2::aod::FDMCParticles&)
void processMixedEventMC(FilteredMCCollisions& cols, o2::aod::FDMCCollisions&, soa::Join<o2::aod::FDParticles, o2::aod::FDMCLabels>& parts, o2::aod::FDMCParticles&)
{
switch (Mixing.Policy.value) {
case femtodreamcollision::kMult:
Expand All @@ -557,7 +565,7 @@ struct femtoDreamPairTaskTrackTrack {
}
PROCESS_SWITCH(femtoDreamPairTaskTrackTrack, processMixedEventMC, "Enable processing mixed events MC", false);

void processMixedEventMCMasked(FilteredMaskedCollisions& cols, soa::Join<o2::aod::FDParticles, o2::aod::FDMCLabels>& parts, o2::aod::FDMCParticles&)
void processMixedEventMCMasked(FilteredMaskedMCCollisions& cols, o2::aod::FDMCCollisions&, soa::Join<o2::aod::FDParticles, o2::aod::FDMCLabels>& parts, o2::aod::FDMCParticles&)
{
switch (Mixing.Policy.value) {
case femtodreamcollision::kMult:
Expand Down
15 changes: 8 additions & 7 deletions PWGCF/FemtoDream/Tasks/femtoDreamPairTaskTrackV0.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -69,10 +69,11 @@ struct femtoDreamPairTaskTrackV0 {
Filter EventMultiplicity = aod::femtodreamcollision::multNtr >= EventSel.MultMin && aod::femtodreamcollision::multNtr <= EventSel.MultMax;
Filter EventMultiplicityPercentile = aod::femtodreamcollision::multV0M >= EventSel.MultPercentileMin && aod::femtodreamcollision::multV0M <= EventSel.MultPercentileMax;

using FilteredCollisions = soa::Filtered<FDCollisions>;
using FilteredCollision = FilteredCollisions::iterator;
using FilteredMaskedCollisions = soa::Filtered<soa::Join<FDCollisions, FDColMasks, FDDownSample>>;
using FilteredMaskedCollision = FilteredMaskedCollisions::iterator;
using FilteredMaskedMCCollisions = soa::Filtered<soa::Join<FDCollisions, aod::FDMCCollLabels, FDColMasks, FDDownSample>>;
using FilteredMaskedMCCollision = FilteredMaskedMCCollisions::iterator;

femtodreamcollision::BitMaskType BitMask = -1;

/// Particle 1 (track)
Expand Down Expand Up @@ -205,7 +206,7 @@ struct femtoDreamPairTaskTrackV0 {

void init(InitContext& context)
{
eventHisto.init(&qaRegistry);
eventHisto.init(&qaRegistry, Option.IsMC);
// void init(HistogramRegistry* registry,
// T& MomentumBins, T& tempFitVarBins, T& NsigmaTPCBins, T& NsigmaTOFBins, T& NsigmaTPCTOFBins, T& InvMassBins,
// bool isMC, int pdgCode, bool isDebug = false)
Expand Down Expand Up @@ -332,19 +333,19 @@ struct femtoDreamPairTaskTrackV0 {
if ((col.bitmaskTrackOne() & BitMask) != BitMask && (col.bitmaskTrackTwo() & BitMask) != BitMask) {
return;
}
eventHisto.fillQA(col);
eventHisto.fillQA<false>(col);
auto SliceTrk1 = PartitionTrk1->sliceByCached(aod::femtodreamparticle::fdCollisionId, col.globalIndex(), cache);
auto SliceV02 = PartitionV02->sliceByCached(aod::femtodreamparticle::fdCollisionId, col.globalIndex(), cache);
doSameEvent<false>(SliceTrk1, SliceV02, parts, col);
}
PROCESS_SWITCH(femtoDreamPairTaskTrackV0, processSameEventMasked, "Enable processing same event with masks", true);

void processSameEventMCMasked(FilteredMaskedCollision const& col, FilteredFDMCParts const& parts, o2::aod::FDMCParticles const&)
void processSameEventMCMasked(FilteredMaskedMCCollision const& col, o2::aod::FDMCCollisions&, FilteredFDMCParts const& parts, o2::aod::FDMCParticles const&)
{
if ((col.bitmaskTrackOne() & BitMask) != BitMask && (col.bitmaskTrackTwo() & BitMask) != BitMask) {
return;
}
eventHisto.fillQA(col);
eventHisto.fillQA<true>(col);
auto SliceMCTrk1 = PartitionMCTrk1->sliceByCached(aod::femtodreamparticle::fdCollisionId, col.globalIndex(), cache);
auto SliceMCV02 = PartitionMCV02->sliceByCached(aod::femtodreamparticle::fdCollisionId, col.globalIndex(), cache);
doSameEvent<true>(SliceMCTrk1, SliceMCV02, parts, col);
Expand Down Expand Up @@ -410,7 +411,7 @@ struct femtoDreamPairTaskTrackV0 {
}
PROCESS_SWITCH(femtoDreamPairTaskTrackV0, processMixedEventMasked, "Enable processing mixed events with masks", true);

void processMixedEventMCMasked(FilteredMaskedCollisions const& cols, FilteredFDMCParts const& parts, o2::aod::FDMCParticles const&)
void processMixedEventMCMasked(FilteredMaskedMCCollisions const& cols, o2::aod::FDMCCollisions&, FilteredFDMCParts const& parts, o2::aod::FDMCParticles const&)
{
switch (Mixing.Policy.value) {
case femtodreamcollision::kMult:
Expand Down
Loading