@@ -53,9 +53,7 @@ namespace o2::aod
5353
5454namespace dqanalysisflags
5555{
56- // DECLARE_SOA_COLUMN(IsEventSelected, isEventSelected, int);
5756DECLARE_SOA_BITMAP_COLUMN (IsEventSelected, isEventSelected, 8 ); // ! Event decision
58- // DECLARE_SOA_COLUMN(IsMuonSelected, isMuonSelected, int);
5957DECLARE_SOA_BITMAP_COLUMN (IsMuonSelected, isMuonSelected, 32 ); // ! Muon track decisions (joinable to ReducedMuonsAssoc)
6058}
6159
@@ -64,16 +62,10 @@ DECLARE_SOA_TABLE(MuonTrackCuts, "AOD", "DQANAMUONCUTSA", dqanalysisflags::IsMuo
6462}
6563
6664// Declarations of various short names
67- using MyEvents = soa::Join<aod::ReducedEvents, aod::ReducedEventsExtended>;
6865using MyEventsSelected = soa::Join<aod::ReducedEvents, aod::ReducedEventsExtended, aod::EventCuts>;
69- using MyEventsVtxCov = soa::Join<aod::ReducedEvents, aod::ReducedEventsExtended, aod::ReducedEventsVtxCov>;
70- using MyEventsVtxCovSelected = soa::Join<aod::ReducedEvents, aod::ReducedEventsExtended, aod::ReducedEventsVtxCov, aod::EventCuts>;
7166
7267using MyPairCandidatesSelected = soa::Join<aod::Dimuons, aod::DimuonsExtra>;
73- using MyMuonTracks = soa::Join<aod::ReducedMuons, aod::ReducedMuonsExtra>;
74- using MyMuonTracksSelected = soa::Join<aod::ReducedMuons, aod::ReducedMuonsExtra, aod::MuonTrackCuts>;
75- using MyMuonTracksWithCov = soa::Join<aod::ReducedMuons, aod::ReducedMuonsExtra, aod::ReducedMuonsCov>;
76- using MyMuonTracksSelectedWithCov = soa::Join<aod::ReducedMuons, aod::ReducedMuonsExtra, aod::ReducedMuonsCov, aod::MuonTrackCuts>;
68+ using MyMuonTracksSelected = soa::Join<aod::ReducedMuonsAssoc, aod::MuonTrackCuts>;
7769
7870// bit maps used for the Fill functions of the VarManager
7971constexpr static uint32_t gkEventFillMap = VarManager::ObjTypes::ReducedEvent | VarManager::ObjTypes::ReducedEventExtended;
@@ -103,6 +95,9 @@ struct DqJPsiMuonCorrelations {
10395 float * fValuesMuon ;
10496
10597 HistogramRegistry registry{" registry" };
98+
99+ // int nMuons;
100+ // int nEvents;
106101
107102 void init (o2::framework::InitContext&)
108103 {
@@ -114,7 +109,8 @@ struct DqJPsiMuonCorrelations {
114109 fValuesMuon = new float [VarManager::kNVars ];
115110 VarManager::SetDefaultVarNames ();
116111
117- LOG (info) << " Initialized DqJPsiMuonCorrelations task" << std::endl;
112+ // nMuons = 0;
113+ // nEvents = 0;
118114 }
119115
120116 // Template function to run pair - muon combinations
@@ -127,16 +123,40 @@ struct DqJPsiMuonCorrelations {
127123 VarManager::FillEvent<TEventFillMap>(event, fValuesMuon );
128124 VarManager::FillEvent<TEventFillMap>(event, fValuesDilepton );
129125
130- if (dileptons.size () > 0 ) {
131- LOG (info) << " Dileptons in event: " << dileptons.size () << std::endl;
132- } else {
133- LOG (info) << " No dileptons in event" << std::endl;
126+ if (!event.isEventSelected_bit (0 )) {
127+ return ;
134128 }
129+
130+ // if (dileptons.size() > 0) {
131+ // LOG(info) << "Processing event " << event.globalIndex() << " with " << muons.size() << " muons and " << dileptons.size() << " dileptons" << std::endl;
132+ // LOG(info) << "Dilepton leg indexes: ";
133+ // for (auto& dilepton : dileptons) {
134+ // LOG(info) << dilepton.index0Id() << " and " << dilepton.index1Id() << " ";
135+ // }
136+ // LOG(info) << std::endl;
137+ // LOG(info) << "Muon indexes: ";
138+ // for (auto& muon : muons) {
139+ // LOG(info) << muon.reducedmuonId() << " ";
140+ // }
141+ // LOG(info) << std::endl;
142+ // }
143+
144+ // for (auto& muon : muons) {
145+ // // VarManager::FillTrack<TMuonFillMap>(muon, fValuesMuon);
146+ // if (!muon.isMuonSelected_raw()) {
147+ // continue;
148+ // }
149+ // nMuons++;
150+ // LOG(info) << "Total number of muons processed: " << nMuons << std::endl;
151+ // }
152+ // nEvents++;
153+ // LOG(info) << "Total number of events processed: " << nEvents << std::endl;
154+
155+
135156 }
136157
137158 void processSkimmedDimuon (MyEventsSelected::iterator const & event, MyMuonTracksSelected const & muons, soa::Filtered<MyPairCandidatesSelected> const & dileptons)
138159 {
139- LOG (info) << " Processing event " << event.globalIndex () << " with " << muons.size () << " muons and " << dileptons.size () << " dileptons" << std::endl;
140160 runDileptonMuon<VarManager::kDecayToMuMu , gkEventFillMap, gkMuonFillMap>(event, muons, dileptons);
141161 }
142162 void processDummy (MyEvents&)
0 commit comments