Skip to content

Commit dfe81e2

Browse files
dsekihatdsekihat
authored
[PWGEM/Dilepton] update resolution map vs. centrality (#9146)
1 parent c84fe86 commit dfe81e2

File tree

3 files changed

+132
-85
lines changed

3 files changed

+132
-85
lines changed

‎PWGEM/Dilepton/Tasks/createResolutionMap.cxx‎

Lines changed: 31 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -18,6 +18,7 @@
1818
#include <utility>
1919
#include <set>
2020
#include <vector>
21+
#include <array>
2122

2223
#include "Framework/runDataProcessing.h"
2324
#include "Framework/AnalysisTask.h"
@@ -28,6 +29,7 @@
2829
#include "Framework/HistogramRegistry.h"
2930
#include "Common/DataModel/TrackSelectionTables.h"
3031
#include "Common/DataModel/EventSelection.h"
32+
#include "Common/DataModel/Centrality.h"
3133

3234
#include "CCDB/BasicCCDBManager.h"
3335
#include "DataFormatsParameters/GRPMagField.h"
@@ -49,6 +51,9 @@ using namespace o2::soa;
4951
using MyCollisions = Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels>;
5052
using MyCollision = MyCollisions::iterator;
5153

54+
using MyCollisionsCent = soa::Join<MyCollisions, aod::CentFT0Ms, aod::CentFT0As, aod::CentFT0Cs>;
55+
using MyCollisionCent = MyCollisionsCent::iterator;
56+
5257
using MyMCTracks = soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::TracksCov, aod::McTrackLabels>;
5358
using MyMCTrack = MyMCTracks::iterator;
5459

@@ -70,8 +75,10 @@ struct CreateResolutionMap {
7075
Configurable<std::string> grpmagPath{"grpmagPath", "GLO/Config/GRPMagField", "CCDB path of the GRPMagField object"};
7176
Configurable<std::string> geoPath{"geoPath", "GLO/Config/GeometryAligned", "Path of the geometry file"};
7277
Configurable<int> cfgEventGeneratorType{"cfgEventGeneratorType", -1, "if positive, select event generator type. i.e. gap or signal"};
78+
Configurable<int> cfgCentEstimator{"cfgCentEstimator", 2, "FT0M:0, FT0A:1, FT0C:2"};
7379

7480
ConfigurableAxis ConfPtGenBins{"ConfPtGenBins", {VARIABLE_WIDTH, 0.00, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.07, 0.08, 0.09, 0.10, 0.11, 0.12, 0.13, 0.14, 0.15, 0.16, 0.17, 0.18, 0.19, 0.20, 0.21, 0.22, 0.23, 0.24, 0.25, 0.26, 0.27, 0.28, 0.29, 0.30, 0.31, 0.32, 0.33, 0.34, 0.35, 0.36, 0.37, 0.38, 0.39, 0.40, 0.41, 0.42, 0.43, 0.44, 0.45, 0.46, 0.47, 0.48, 0.49, 0.50, 0.51, 0.52, 0.53, 0.54, 0.55, 0.56, 0.57, 0.58, 0.59, 0.60, 0.61, 0.62, 0.63, 0.64, 0.65, 0.66, 0.67, 0.68, 0.69, 0.70, 0.71, 0.72, 0.73, 0.74, 0.75, 0.76, 0.77, 0.78, 0.79, 0.80, 0.81, 0.82, 0.83, 0.84, 0.85, 0.86, 0.87, 0.88, 0.89, 0.90, 0.91, 0.92, 0.93, 0.94, 0.95, 0.96, 0.97, 0.98, 0.99, 1.00, 1.10, 1.20, 1.30, 1.40, 1.50, 1.60, 1.70, 1.80, 1.90, 2.00, 2.10, 2.20, 2.30, 2.40, 2.50, 2.60, 2.70, 2.80, 2.90, 3.00, 3.10, 3.20, 3.30, 3.40, 3.50, 3.60, 3.70, 3.80, 3.90, 4.00, 4.10, 4.20, 4.30, 4.40, 4.50, 4.60, 4.70, 4.80, 4.90, 5.00, 5.50, 6.00, 6.50, 7.00, 7.50, 8.00, 8.50, 9.00, 9.50, 10.00, 11.00, 12.00, 13.00, 14.00, 15.00, 16.00, 17.00, 18.00, 19.00, 20.00}, "gen. pT bins for output histograms"};
81+
ConfigurableAxis ConfCentBins{"ConfCentBins", {VARIABLE_WIDTH, 0, 10, 30, 50, 110}, "centrality (%) bins for output histograms"};
7582

7683
ConfigurableAxis ConfEtaCBGenBins{"ConfEtaCBGenBins", {30, -1.5, +1.5}, "gen. eta bins at midrapidity for output histograms"};
7784
ConfigurableAxis ConfEtaFWDGenBins{"ConfEtaFWDGenBins", {40, -5.5, -1.5}, "gen. eta bins at forward rapidity for output histograms"};
@@ -88,6 +95,7 @@ struct CreateResolutionMap {
8895
Configurable<float> cfg_max_eta_track{"cfg_max_eta_track", +1.5, "max eta for single track"};
8996
Configurable<int> cfg_min_ncluster_tpc{"cfg_min_ncluster_tpc", 0, "min ncluster tpc"};
9097
Configurable<int> cfg_min_ncluster_its{"cfg_min_ncluster_its", 4, "min ncluster its"};
98+
Configurable<int> cfg_min_ncluster_itsib{"cfg_min_ncluster_itsib", 1, "min ncluster itsib"};
9199
Configurable<int> cfg_min_ncrossedrows{"cfg_min_ncrossedrows", 80, "min ncrossed rows"};
92100
Configurable<float> cfg_max_chi2tpc{"cfg_max_chi2tpc", 4.0, "max chi2/NclsTPC"};
93101
Configurable<float> cfg_max_chi2its{"cfg_max_chi2its", 5.0, "max chi2/NclsITS"};
@@ -129,6 +137,7 @@ struct CreateResolutionMap {
129137
ccdb->setFatalWhenNull(false);
130138
ccdbApi.init(ccdburl);
131139

140+
const AxisSpec axis_cent{ConfCentBins, "centrality (%)"};
132141
const AxisSpec axis_pt_gen{ConfPtGenBins, "p_{T,l}^{gen} (GeV/c)"};
133142
const AxisSpec axis_eta_cb_gen{ConfEtaCBGenBins, "#eta_{l}^{gen}"};
134143
const AxisSpec axis_eta_fwd_gen{ConfEtaFWDGenBins, "#eta_{l}^{gen}"};
@@ -145,9 +154,9 @@ struct CreateResolutionMap {
145154
registry.addClone("Electron/", "StandaloneMuon/");
146155
registry.addClone("Electron/", "GlobalMuon/");
147156

148-
registry.add("Electron/hs_reso", "7D resolution positive", kTHnSparseF, {axis_pt_gen, axis_eta_cb_gen, axis_phi_gen, axis_charge_gen, axis_dpt, axis_deta, axis_dphi}, true);
149-
registry.add("StandaloneMuon/hs_reso", "7D resolution positive", kTHnSparseF, {axis_pt_gen, axis_eta_fwd_gen, axis_phi_gen, axis_charge_gen, axis_dpt, axis_deta, axis_dphi}, true);
150-
registry.add("GlobalMuon/hs_reso", "7D resolution positive", kTHnSparseF, {axis_pt_gen, axis_eta_fwd_gen, axis_phi_gen, axis_charge_gen, axis_dpt, axis_deta, axis_dphi}, true);
157+
registry.add("Electron/hs_reso", "8D resolution positive", kTHnSparseF, {axis_cent, axis_pt_gen, axis_eta_cb_gen, axis_phi_gen, axis_charge_gen, axis_dpt, axis_deta, axis_dphi}, true);
158+
registry.add("StandaloneMuon/hs_reso", "8D resolution positive", kTHnSparseF, {axis_cent, axis_pt_gen, axis_eta_fwd_gen, axis_phi_gen, axis_charge_gen, axis_dpt, axis_deta, axis_dphi}, true);
159+
registry.add("GlobalMuon/hs_reso", "8D resolution positive", kTHnSparseF, {axis_cent, axis_pt_gen, axis_eta_fwd_gen, axis_phi_gen, axis_charge_gen, axis_dpt, axis_deta, axis_dphi}, true);
151160
}
152161

153162
void initCCDB(aod::BCsWithTimestamps::iterator const& bc)
@@ -185,6 +194,9 @@ struct CreateResolutionMap {
185194
if (track.itsNCls() < electroncuts.cfg_min_ncluster_its) {
186195
return false;
187196
}
197+
if (track.itsNClsInnerBarrel() < electroncuts.cfg_min_ncluster_itsib) {
198+
return false;
199+
}
188200

189201
auto hits = std::count_if(itsRequirement_ibany.second.begin(), itsRequirement_ibany.second.end(), [&](auto&& requiredLayer) { return track.itsClusterMap() & (1 << requiredLayer); });
190202
if (hits < itsRequirement_ibany.first) {
@@ -268,7 +280,7 @@ struct CreateResolutionMap {
268280
}
269281

270282
template <typename TMuon, typename TCollision>
271-
bool checkFwdTrack(TMuon const& muon, TCollision const& collision)
283+
bool checkFwdTrack(TMuon const& muon, TCollision const& collision, const float centrality)
272284
{
273285
o2::dataformats::GlobalFwdTrack propmuonAtPV = PropagateMuon(muon, collision, CreateResolutionMap::MuonExtrapolation::kToVertex);
274286
float pt = propmuonAtPV.getPt();
@@ -316,7 +328,7 @@ struct CreateResolutionMap {
316328

317329
auto mctrack = muon.template mcParticle_as<aod::McParticles>();
318330
if (muon.trackType() == static_cast<uint8_t>(o2::aod::fwdtrack::ForwardTrackTypeEnum::MuonStandaloneTrack)) {
319-
registry.fill(HIST("StandaloneMuon/hs_reso"), mctrack.pt(), mctrack.eta(), mctrack.phi(), -mctrack.pdgCode() / 13, (mctrack.pt() - pt) / mctrack.pt(), mctrack.eta() - eta, mctrack.phi() - phi);
331+
registry.fill(HIST("StandaloneMuon/hs_reso"), centrality, mctrack.pt(), mctrack.eta(), mctrack.phi(), -mctrack.pdgCode() / 13, (mctrack.pt() - pt) / mctrack.pt(), mctrack.eta() - eta, mctrack.phi() - phi);
320332
registry.fill(HIST("StandaloneMuon/Ptgen_RelDeltaPt"), mctrack.pt(), (mctrack.pt() - pt) / mctrack.pt());
321333
registry.fill(HIST("StandaloneMuon/Ptgen_DeltaEta"), mctrack.pt(), mctrack.eta() - eta);
322334
if (mctrack.pdgCode() == -13) { // positive muon
@@ -325,7 +337,7 @@ struct CreateResolutionMap {
325337
registry.fill(HIST("StandaloneMuon/Ptgen_DeltaPhi_Neg"), mctrack.pt(), mctrack.phi() - phi);
326338
}
327339
} else if (muon.trackType() == static_cast<uint8_t>(o2::aod::fwdtrack::ForwardTrackTypeEnum::GlobalMuonTrack)) {
328-
registry.fill(HIST("GlobalMuon/hs_reso"), mctrack.pt(), mctrack.eta(), mctrack.phi(), -mctrack.pdgCode() / 13, (mctrack.pt() - pt) / mctrack.pt(), mctrack.eta() - eta, mctrack.phi() - phi);
340+
registry.fill(HIST("GlobalMuon/hs_reso"), centrality, mctrack.pt(), mctrack.eta(), mctrack.phi(), -mctrack.pdgCode() / 13, (mctrack.pt() - pt) / mctrack.pt(), mctrack.eta() - eta, mctrack.phi() - phi);
329341
registry.fill(HIST("GlobalMuon/Ptgen_RelDeltaPt"), mctrack.pt(), (mctrack.pt() - pt) / mctrack.pt());
330342
registry.fill(HIST("GlobalMuon/Ptgen_DeltaEta"), mctrack.pt(), mctrack.eta() - eta);
331343
if (mctrack.pdgCode() == -13) { // positive muon
@@ -343,14 +355,16 @@ struct CreateResolutionMap {
343355

344356
Filter collisionFilter = o2::aod::evsel::sel8 == true && nabs(o2::aod::collision::posZ) < 10.f;
345357
using MyFilteredCollisions = soa::Filtered<MyCollisions>;
358+
using MyFilteredCollisionsCent = soa::Filtered<MyCollisionsCent>;
346359

347360
Filter trackFilter_mid = o2::aod::track::pt > electroncuts.cfg_min_pt_track&& electroncuts.cfg_min_eta_track < o2::aod::track::eta&& o2::aod::track::eta < electroncuts.cfg_max_eta_track&& o2::aod::track::tpcChi2NCl < electroncuts.cfg_max_chi2tpc&& o2::aod::track::itsChi2NCl < electroncuts.cfg_max_chi2its&& ncheckbit(aod::track::v001::detectorMap, (uint8_t)o2::aod::track::ITS) == true && ncheckbit(aod::track::v001::detectorMap, (uint8_t)o2::aod::track::TPC) == true && nabs(o2::aod::track::dcaXY) < electroncuts.cfg_max_dcaxy&& nabs(o2::aod::track::dcaZ) < electroncuts.cfg_max_dcaz;
348361
using MyFilteredMCTracks = soa::Filtered<MyMCTracks>;
349362

350363
Partition<MyMCFwdTracks> sa_muons = o2::aod::fwdtrack::trackType == uint8_t(o2::aod::fwdtrack::ForwardTrackTypeEnum::MuonStandaloneTrack); // MCH-MID
351364
Partition<MyMCFwdTracks> global_muons = o2::aod::fwdtrack::trackType == uint8_t(o2::aod::fwdtrack::ForwardTrackTypeEnum::GlobalMuonTrack); // MFT-MCH-MID
352365

353-
void process(MyFilteredCollisions const& collisions, aod::BCsWithTimestamps const&, MyFilteredMCTracks const& tracks, MyMCFwdTracks const&, aod::McCollisions const&, aod::McParticles const&)
366+
template <typename TCollisions>
367+
void process(TCollisions const& collisions, aod::BCsWithTimestamps const&, MyFilteredMCTracks const& tracks, MyMCFwdTracks const&, aod::McCollisions const&, aod::McParticles const&)
354368
{
355369
for (auto& collision : collisions) {
356370
auto bc = collision.template foundBC_as<aod::BCsWithTimestamps>();
@@ -365,6 +379,11 @@ struct CreateResolutionMap {
365379
continue;
366380
}
367381

382+
float centrality = 105.f;
383+
if constexpr (std::is_same_v<std::decay_t<TCollisions>, MyFilteredCollisionsCent>) {
384+
centrality = std::array{collision.centFT0M(), collision.centFT0A(), collision.centFT0C()}[cfgCentEstimator];
385+
}
386+
368387
auto tracks_per_coll = tracks.sliceBy(perCollision_mid, collision.globalIndex());
369388
for (auto& track : tracks_per_coll) {
370389
if (!track.has_mcParticle()) {
@@ -381,7 +400,7 @@ struct CreateResolutionMap {
381400
continue;
382401
}
383402

384-
registry.fill(HIST("Electron/hs_reso"), mctrack.pt(), mctrack.eta(), mctrack.phi(), -mctrack.pdgCode() / 11, (mctrack.pt() - track.pt()) / mctrack.pt(), mctrack.eta() - track.eta(), mctrack.phi() - track.phi());
403+
registry.fill(HIST("Electron/hs_reso"), centrality, mctrack.pt(), mctrack.eta(), mctrack.phi(), -mctrack.pdgCode() / 11, (mctrack.pt() - track.pt()) / mctrack.pt(), mctrack.eta() - track.eta(), mctrack.phi() - track.phi());
385404
registry.fill(HIST("Electron/Ptgen_RelDeltaPt"), mctrack.pt(), (mctrack.pt() - track.pt()) / mctrack.pt());
386405
registry.fill(HIST("Electron/Ptgen_DeltaEta"), mctrack.pt(), mctrack.eta() - track.eta());
387406
if (mctrack.pdgCode() == -11) { // positron
@@ -406,7 +425,7 @@ struct CreateResolutionMap {
406425
if (abs(mctrack.pdgCode()) != 13 || !(mctrack.isPhysicalPrimary() || mctrack.producedByGenerator())) {
407426
continue;
408427
}
409-
if (!checkFwdTrack(muon, collision)) {
428+
if (!checkFwdTrack(muon, collision, centrality)) {
410429
continue;
411430
}
412431
} // end of standalone muon loop
@@ -422,14 +441,16 @@ struct CreateResolutionMap {
422441
if (abs(mctrack.pdgCode()) != 13 || !(mctrack.isPhysicalPrimary() || mctrack.producedByGenerator())) {
423442
continue;
424443
}
425-
if (!checkFwdTrack(muon, collision)) {
444+
if (!checkFwdTrack(muon, collision, centrality)) {
426445
continue;
427446
}
428447

429448
} // end of global muon loop
430449

431450
} // end of collision loop
432451
}
452+
PROCESS_SWITCH_FULL(CreateResolutionMap, process<MyFilteredCollisionsCent>, processWithCent, "create resolution map wit centrality", true);
453+
PROCESS_SWITCH_FULL(CreateResolutionMap, process<MyFilteredCollisions>, processWithoutCent, "create resolution map without centrality", false);
433454
};
434455
WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)
435456
{

‎PWGEM/Dilepton/Tasks/smearing.cxx‎

Lines changed: 37 additions & 25 deletions
Original file line numberDiff line numberDiff line change
@@ -24,7 +24,7 @@
2424
#include "Framework/ASoA.h"
2525
#include "Framework/DataTypes.h"
2626
#include "Framework/HistogramRegistry.h"
27-
#include "PWGDQ/DataModel/ReducedInfoTables.h" // remove this later, because 2 data tables (covariant matrix) in this header confilict against EM tables.
27+
// #include "PWGDQ/DataModel/ReducedInfoTables.h" // remove this later, because 2 data tables (covariant matrix) in this header confilict against EM tables.
2828
#include "PWGEM/Dilepton/Utils/MomentumSmearer.h"
2929
#include "PWGEM/Dilepton/DataModel/dileptonTables.h"
3030

@@ -34,6 +34,11 @@ using namespace o2::framework::expressions;
3434
using namespace o2::aod;
3535

3636
struct ApplySmearing {
37+
enum class EMAnaType : int {
38+
kEfficiency = 0,
39+
kCocktail = 1,
40+
};
41+
3742
Produces<aod::SmearedElectrons> smearedelectron;
3843
Produces<aod::SmearedMuons> smearedmuon;
3944

@@ -167,8 +172,8 @@ struct ApplySmearing {
167172
smearer_GlobalMuon.init();
168173
}
169174

170-
template <typename TTracksMC>
171-
void applySmearing(TTracksMC const& tracksMC)
175+
template <EMAnaType type, typename TTracksMC, typename TCollisions, typename TMCCollisions>
176+
void applySmearing(TTracksMC const& tracksMC, TCollisions const&, TMCCollisions const&)
172177
{
173178
for (auto& mctrack : tracksMC) {
174179
float ptgen = mctrack.pt();
@@ -177,6 +182,13 @@ struct ApplySmearing {
177182
float efficiency = 1.;
178183
float dca = 0.;
179184

185+
float centrality = 105.f;
186+
if constexpr (type == EMAnaType::kEfficiency) {
187+
centrality = 0; // to be implemented later.
188+
} else {
189+
centrality = 0;
190+
}
191+
180192
int pdgCode = mctrack.pdgCode();
181193
if (abs(pdgCode) == 11) {
182194
int ch = -1;
@@ -185,7 +197,7 @@ struct ApplySmearing {
185197
}
186198
// apply smearing for electrons or muons.
187199
float ptsmeared, etasmeared, phismeared;
188-
smearer_Electron.applySmearing(ch, ptgen, etagen, phigen, ptsmeared, etasmeared, phismeared);
200+
smearer_Electron.applySmearing(centrality, ch, ptgen, etagen, phigen, ptsmeared, etasmeared, phismeared);
189201
// get the efficiency
190202
efficiency = smearer_Electron.getEfficiency(ptgen, etagen, phigen);
191203
// get DCA
@@ -200,13 +212,13 @@ struct ApplySmearing {
200212
}
201213
// apply smearing for muons based on resolution map of standalone muons
202214
float ptsmeared_sa = 0.f, etasmeared_sa = 0.f, phismeared_sa = 0.f, efficiency_sa = 1.f, dca_sa = 0.f;
203-
smearer_StandaloneMuon.applySmearing(ch, ptgen, etagen, phigen, ptsmeared_sa, etasmeared_sa, phismeared_sa);
215+
smearer_StandaloneMuon.applySmearing(centrality, ch, ptgen, etagen, phigen, ptsmeared_sa, etasmeared_sa, phismeared_sa);
204216
efficiency_sa = smearer_StandaloneMuon.getEfficiency(ptgen, etagen, phigen);
205217
dca_sa = smearer_StandaloneMuon.getDCA(ptsmeared_sa);
206218

207219
float ptsmeared_gl = 0.f, etasmeared_gl = 0.f, phismeared_gl = 0.f, efficiency_gl = 1.f, dca_gl = 0.f;
208220
// apply smearing for muons based on resolution map of global muons
209-
smearer_GlobalMuon.applySmearing(ch, ptgen, etagen, phigen, ptsmeared_gl, etasmeared_gl, phismeared_gl);
221+
smearer_GlobalMuon.applySmearing(centrality, ch, ptgen, etagen, phigen, ptsmeared_gl, etasmeared_gl, phismeared_gl);
210222
efficiency_gl = smearer_GlobalMuon.getEfficiency(ptgen, etagen, phigen);
211223
dca_gl = smearer_GlobalMuon.getDCA(ptsmeared_gl);
212224
smearedmuon(ptsmeared_sa, etasmeared_sa, phismeared_sa, efficiency_sa, dca_sa, ptsmeared_gl, etasmeared_gl, phismeared_gl, efficiency_gl, dca_gl);
@@ -220,19 +232,19 @@ struct ApplySmearing {
220232
}
221233
}
222234

223-
void processMCanalysisEM(aod::EMMCParticles const& tracksMC)
235+
void processMCanalysisEM(aod::EMMCParticles const& tracksMC, aod::EMEvents const& collisions, aod::EMMCEvents const& mccollisions)
224236
{
225-
applySmearing(tracksMC);
237+
applySmearing<EMAnaType::kEfficiency>(tracksMC, collisions, mccollisions);
226238
}
227239

228-
void processMCanalysisDQ(ReducedMCTracks const& tracksMC)
229-
{
230-
applySmearing(tracksMC);
231-
}
240+
// void processMCanalysisDQ(ReducedMCTracks const& tracksMC)
241+
// {
242+
// applySmearing<EMAnaType::kEfficiency>(tracksMC);
243+
// }
232244

233245
void processCocktail(aod::McParticles const& tracksMC)
234246
{
235-
applySmearing(tracksMC);
247+
applySmearing<EMAnaType::kCocktail>(tracksMC, nullptr, nullptr);
236248
}
237249

238250
void processDummyCocktail(aod::McParticles const& tracksMC)
@@ -252,20 +264,20 @@ struct ApplySmearing {
252264
}
253265

254266
void processDummyMCanalysisEM(aod::EMMCParticles const&) {}
255-
void processDummyMCanalysisDQ(ReducedMCTracks const&) {}
267+
// void processDummyMCanalysisDQ(ReducedMCTracks const&) {}
256268

257269
PROCESS_SWITCH(ApplySmearing, processMCanalysisEM, "Run for MC analysis which uses skimmed EM data format", false);
258-
PROCESS_SWITCH(ApplySmearing, processMCanalysisDQ, "Run for MC analysis which uses skimmed DQ data format", false);
270+
// PROCESS_SWITCH(ApplySmearing, processMCanalysisDQ, "Run for MC analysis which uses skimmed DQ data format", false);
259271
PROCESS_SWITCH(ApplySmearing, processCocktail, "Run for cocktail analysis", false);
260272
PROCESS_SWITCH(ApplySmearing, processDummyMCanalysisEM, "Dummy process function", false);
261-
PROCESS_SWITCH(ApplySmearing, processDummyMCanalysisDQ, "Dummy process function", false);
273+
// PROCESS_SWITCH(ApplySmearing, processDummyMCanalysisDQ, "Dummy process function", false);
262274
PROCESS_SWITCH(ApplySmearing, processDummyCocktail, "Dummy process function", true);
263275
};
264276

265277
struct CheckSmearing {
266278
using EMMCParticlesWithSmearing = soa::Join<aod::EMMCParticles, aod::SmearedElectrons>; // this is only for electrons
267-
using MyReducedTracks = soa::Join<ReducedMCTracks, aod::SmearedElectrons>; // this is only for electrons
268-
using MyCocktailTracks = soa::Join<aod::McParticles, aod::SmearedElectrons>; // this is only for electrons
279+
// using MyReducedTracks = soa::Join<ReducedMCTracks, aod::SmearedElectrons>; // this is only for electrons
280+
using MyCocktailTracks = soa::Join<aod::McParticles, aod::SmearedElectrons>; // this is only for electrons
269281

270282
// Run for electrons or muons
271283
Configurable<int> fPdgCode{"cfgPdgCode", 11, "Set the type of particle to be checked"};
@@ -336,25 +348,25 @@ struct CheckSmearing {
336348
Check(tracksMC);
337349
}
338350

339-
void processCheckMCanalysisDQ(MyReducedTracks const& tracksMC)
340-
{
341-
Check(tracksMC);
342-
}
351+
// void processCheckMCanalysisDQ(MyReducedTracks const& tracksMC)
352+
// {
353+
// Check(tracksMC);
354+
// }
343355

344356
void processCheckCocktail(MyCocktailTracks const& tracksMC)
345357
{
346358
Check(tracksMC);
347359
}
348360

349361
void processDummyMCanalysisEM(aod::EMMCParticles const&) {}
350-
void processDummyMCanalysisDQ(ReducedMCTracks const&) {}
362+
// void processDummyMCanalysisDQ(ReducedMCTracks const&) {}
351363
void processDummyCocktail(aod::McParticles const&) {}
352364

353365
PROCESS_SWITCH(CheckSmearing, processCheckMCanalysisEM, "Run for MC analysis", false);
354-
PROCESS_SWITCH(CheckSmearing, processCheckMCanalysisDQ, "Run for MC analysis", false);
366+
// PROCESS_SWITCH(CheckSmearing, processCheckMCanalysisDQ, "Run for MC analysis", false);
355367
PROCESS_SWITCH(CheckSmearing, processCheckCocktail, "Run for cocktail analysis", false);
356368
PROCESS_SWITCH(CheckSmearing, processDummyMCanalysisEM, "Dummy process function", false);
357-
PROCESS_SWITCH(CheckSmearing, processDummyMCanalysisDQ, "Dummy process function", false);
369+
// PROCESS_SWITCH(CheckSmearing, processDummyMCanalysisDQ, "Dummy process function", false);
358370
PROCESS_SWITCH(CheckSmearing, processDummyCocktail, "Dummy process function", true);
359371
};
360372
WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)

0 commit comments

Comments
 (0)