Skip to content

Commit 969adb6

Browse files
sejeong8alibuild
sejeong8
and
alibuild
authored
[PWGHF] Added dcaHadron Histogram in electron analysis (#14765)
Co-authored-by: ALICE Action Bot <alibuild@cern.ch>
1 parent 88ecd82 commit 969adb6

File tree

1 file changed

+110
-103
lines changed

1 file changed

+110
-103
lines changed

‎PWGHF/HFL/Tasks/taskSingleElectron.cxx‎

Lines changed: 110 additions & 103 deletions
Original file line numberDiff line numberDiff line change
@@ -11,7 +11,7 @@
1111
///
1212
/// \file taskSingleElectron.cxx
1313
/// \brief task for electrons from heavy-flavour hadron decays
14-
/// \author Jonghan Park (Jeonbuk National University)
14+
/// \author Jonghan Park (Jeonbuk National University), Seul I Jeong (Pusan National University)
1515

1616
#include "Common/DataModel/EventSelection.h"
1717
#include "Common/DataModel/PIDResponseTOF.h"
@@ -41,12 +41,6 @@ using namespace o2::constants::physics;
4141
using namespace o2::framework;
4242
using namespace o2::framework::expressions;
4343

44-
enum PdgCode {
45-
kEta = 221,
46-
kOmega = 223,
47-
kEtaPrime = 331
48-
};
49-
5044
enum SourceType {
5145
NotElec = 0, // not electron
5246
DirectCharm = 1, // electrons from prompt charm hadrons
@@ -88,20 +82,23 @@ struct HfTaskSingleElectron {
8882
Configurable<float> ptTrackMax{"ptTrackMax", 10., "max pt cut"};
8983
Configurable<float> ptTrackMin{"ptTrackMin", 0.5, "min pt cut"};
9084
Configurable<float> etaTrackMax{"etaTrackMax", 0.8, "eta cut"};
91-
Configurable<int> tpcNCrossedRowMin{"tpcNCrossedRowMin", 70, "max of TPC n cluster crossed rows"};
92-
Configurable<float> tpcNClsFoundOverFindableMin{"tpcNClsFoundOverFindableMin", 0.8, "min # of TPC found/findable clusters"};
93-
Configurable<float> tpcChi2perNClMax{"tpcChi2perNClMax", 4., "min # of tpc chi2 per clusters"};
94-
Configurable<int> itsIBClsMin{"itsIBClsMin", 3, "min # of its clusters in IB"};
95-
Configurable<float> itsChi2perNClMax{"itsChi2perNClMax", 6., "min # of tpc chi2 per clusters"};
85+
Configurable<int> nCrossedRowTPCMin{"nCrossedRowTPCMin", 70, "max of TPC n cluster crossed rows"};
86+
Configurable<float> nClsFoundOverFindableTpcMin{"nClsFoundOverFindableTpcMin", 0.8, "min # of TPC found/findable clusters"};
87+
Configurable<float> Chi2pernClTpcMax{"Chi2pernClTpcMax", 4., "min # of tpc chi2 per clusters"};
88+
Configurable<int> ClsIBItsMin{"ClsIBItsMin", 3, "min # of its clusters in IB"};
89+
Configurable<float> Chi2pernClItsMax{"Chi2pernClItsMax", 6., "min # of tpc chi2 per clusters"};
9690
Configurable<float> dcaxyMax{"dcaxyMax", 1., "max of track dca in xy"};
9791
Configurable<float> dcazMax{"dcazMax", 2., "max of track dca in z"};
98-
Configurable<float> tofNSigmaMax{"tofNSigmaMax", 3., "max of tof nsigma"};
99-
Configurable<float> tpcNSigmaMin{"tpcNSigmaMin", -1., "min of tpc nsigma"};
100-
Configurable<float> tpcNSigmaMax{"tpcNSigmaMax", 3., "max of tpc nsigma"};
92+
Configurable<float> nSigmaTofMax{"nSigmaTofMax", 3., "max of tof nsigma"};
93+
Configurable<float> nSigmaTpcMin{"nSigmaTpcMin", -1., "min of tpc nsigma"};
94+
Configurable<float> nSigmaTPCMax{"nSigmaTPCMax", 3., "max of tpc nsigma"};
10195

10296
Configurable<int> nBinsP{"nBinsP", 1500, "number of bins of particle momentum"};
10397
Configurable<int> nBinsPt{"nBinsPt", 100, "N bins in pT histo"};
10498

99+
Configurable<int> nSigmaTpcHadronMax{"nSigmaTpcHadronMax", -3, "max of tpc hadron nsigma"};
100+
Configurable<int> nSigmaTpcHadronMin{"nSigmaTpcHadronMin", -5, "min of tpc hadron nsigma"};
101+
105102
// SliceCache
106103
SliceCache cache;
107104

@@ -119,7 +116,7 @@ struct HfTaskSingleElectron {
119116
ConfigurableAxis axisPtEl{"axisPtEl", {VARIABLE_WIDTH, 0.5f, 0.6f, 0.7f, 0.8f, 0.9f, 1.f, 1.1f, 1.2f, 1.3f, 1.4f, 1.5f, 1.75f, 2.0f, 2.25f, 2.5f, 2.75f, 3.f, 3.5f, 4.0f, 5.0f, 6.0f, 8.0f, 10.0f}, "electron pt bins"};
120117

121118
// Histogram registry
122-
HistogramRegistry histos{"histos", {}, OutputObjHandlingPolicy::AnalysisObject};
119+
HistogramRegistry histos{"histos"};
123120

124121
void init(InitContext const&)
125122
{
@@ -134,35 +131,36 @@ struct HfTaskSingleElectron {
134131
const AxisSpec axisTrackIp{4000, -0.2, 0.2, "dca"};
135132

136133
// create histograms
137-
histos.add("nEvents", "Number of events", kTH1D, {{1, 0., 1.}});
138-
histos.add("VtxZ", "VtxZ; cm; entries", kTH1D, {axisPosZ});
139-
histos.add("etaTrack", "etaTrack; #eta; entries", kTH1D, {axisEta});
140-
histos.add("ptTrack", "#it{p}_{T} distribution of selected tracks; #it{p}_{T} (GeV/#it{c}); entries", kTH1D, {axisPt});
134+
histos.add("hNEvents", "Number of events", kTH1D, {{1, 0., 1.}});
135+
histos.add("hVtxZ", "hVtxZ; cm; entries", kTH1D, {axisPosZ});
136+
histos.add("hEtaTrack", "hEtaTrack; #eta; entries", kTH1D, {axisEta});
137+
histos.add("hPtTrack", "#it{p}_{T} distribution of selected tracks; #it{p}_{T} (GeV/#it{c}); entries", kTH1D, {axisPt});
141138

142139
// QA plots for trigger track selection
143-
histos.add("tpcNClsTrack", "tpcNClsTrack", kTH1D, {{200, 0, 200}});
144-
histos.add("tpcFoundFindableTrack", "", kTH1D, {{10, 0, 1}});
145-
histos.add("tpcChi2Track", "", kTH1D, {{100, 0, 10}});
146-
histos.add("itsIBClsTrack", "", kTH1D, {{10, 0, 10}});
147-
histos.add("itsChi2Track", "", kTH1D, {{50, 0, 50}});
148-
histos.add("dcaXYTrack", "", kTH1D, {{600, -3, 3}});
149-
histos.add("dcaZTrack", "", kTH1D, {{600, -3, 3}});
140+
histos.add("hNClsTpcTrack", "hNClsTpcTrack", kTH1D, {{200, 0, 200}});
141+
histos.add("hNClsFoundFindableTpcTrack", "", kTH1D, {{10, 0, 1}});
142+
histos.add("hChi2TpcTrack", "", kTH1D, {{100, 0, 10}});
143+
histos.add("hibClsItsTrack", "", kTH1D, {{10, 0, 10}});
144+
histos.add("hChi2ItsTrack", "", kTH1D, {{50, 0, 50}});
145+
histos.add("hDcaXYTrack", "", kTH1D, {{600, -3, 3}});
146+
histos.add("hDcaZTrack", "", kTH1D, {{600, -3, 3}});
150147

151148
// pid
152-
histos.add("tofNSigPt", "", kTH2D, {{axisPtEl}, {axisNsig}});
153-
histos.add("tofNSigPtQA", "", kTH2D, {{axisPtEl}, {axisNsig}});
154-
histos.add("tpcNSigP", "", kTH2D, {{axisP}, {axisNsig}});
155-
histos.add("tpcNSigPt", "", kTH2D, {{axisPtEl}, {axisNsig}});
156-
histos.add("tpcNSigPAfterTofCut", "", kTH2D, {{axisP}, {axisNsig}});
157-
histos.add("tpcNSigPtAfterTofCut", "", kTH2D, {{axisPtEl}, {axisNsig}});
158-
histos.add("tpcNSigPtQA", "", kTH2D, {{axisPtEl}, {axisNsig}});
149+
histos.add("hTofNSigPt", "", kTH2D, {{axisPtEl}, {axisNsig}});
150+
histos.add("hTofNSigPtQA", "", kTH2D, {{axisPtEl}, {axisNsig}});
151+
histos.add("hTPCNSigP", "", kTH2D, {{axisP}, {axisNsig}});
152+
histos.add("hTPCNSigPt", "", kTH2D, {{axisPtEl}, {axisNsig}});
153+
histos.add("hTPCNSigPAfterTofCut", "", kTH2D, {{axisP}, {axisNsig}});
154+
histos.add("hTPCNSigPtAfterTofCut", "", kTH2D, {{axisPtEl}, {axisNsig}});
155+
histos.add("hTPCNSigPtQA", "", kTH2D, {{axisPtEl}, {axisNsig}});
159156

160157
// track impact parameter
161-
histos.add("dcaTrack", "", kTH2D, {{axisPtEl}, {axisTrackIp}});
162-
histos.add("dcaBeauty", "", kTH2D, {{axisPtEl}, {axisTrackIp}});
163-
histos.add("dcaCharm", "", kTH2D, {{axisPtEl}, {axisTrackIp}});
164-
histos.add("dcaDalitz", "", kTH2D, {{axisPtEl}, {axisTrackIp}});
165-
histos.add("dcaConv", "", kTH2D, {{axisPtEl}, {axisTrackIp}});
158+
histos.add("hDcaTrack", "", kTH2D, {{axisPtEl}, {axisTrackIp}});
159+
histos.add("hDcaBeauty", "", kTH2D, {{axisPtEl}, {axisTrackIp}});
160+
histos.add("hDcaChram", "", kTH2D, {{axisPtEl}, {axisTrackIp}});
161+
histos.add("hDcaDalitz", "", kTH2D, {{axisPtEl}, {axisTrackIp}});
162+
histos.add("hDcaConv", "", kTH2D, {{axisPtEl}, {axisTrackIp}});
163+
histos.add("hDcaHadron", "", kTH2D, {{axisPtEl}, {axisTrackIp}});
166164

167165
// QA plots for MC
168166
histos.add("hPdgC", "", kTH1D, {{10001, -0.5, 10000.5}});
@@ -181,23 +179,23 @@ struct HfTaskSingleElectron {
181179
return false;
182180
}
183181

184-
if (track.tpcNClsCrossedRows() < tpcNCrossedRowMin) {
182+
if (track.tpcNClsCrossedRows() < nCrossedRowTPCMin) {
185183
return false;
186184
}
187185

188-
if (track.tpcCrossedRowsOverFindableCls() < tpcNClsFoundOverFindableMin) {
186+
if (track.tpcCrossedRowsOverFindableCls() < nClsFoundOverFindableTpcMin) {
189187
return false;
190188
}
191189

192-
if (track.tpcChi2NCl() > tpcChi2perNClMax) {
190+
if (track.tpcChi2NCl() > Chi2pernClTpcMax) {
193191
return false;
194192
}
195193

196-
if (!(track.itsNClsInnerBarrel() == itsIBClsMin)) {
194+
if (!(track.itsNClsInnerBarrel() == ClsIBItsMin)) {
197195
return false;
198196
}
199197

200-
if (track.itsChi2NCl() > itsChi2perNClMax) {
198+
if (track.itsChi2NCl() > Chi2pernClItsMax) {
201199
return false;
202200
}
203201

@@ -279,16 +277,16 @@ struct HfTaskSingleElectron {
279277
if (grmotherPdg == kPi0) {
280278
return GammaPi0;
281279
}
282-
if (grmotherPdg == PdgCode::kEta) {
280+
if (grmotherPdg == Pdg::kEta) {
283281
return GammaEta;
284282
}
285-
if (grmotherPdg == PdgCode::kOmega) {
283+
if (grmotherPdg == Pdg::kOmega) {
286284
return GammaOmega;
287285
}
288-
if (grmotherPdg == kPhi) {
286+
if (grmotherPdg == Pdg::kPhi) {
289287
return GammaPhi;
290288
}
291-
if (grmotherPdg == PdgCode::kEtaPrime) {
289+
if (grmotherPdg == Pdg::kEtaPrime) {
292290
return GammaEtaPrime;
293291
}
294292
if (grmotherPdg == kRho770_0) {
@@ -320,22 +318,22 @@ struct HfTaskSingleElectron {
320318
mpt = grmotherPt;
321319
return GammaPi0;
322320
}
323-
if (grmotherPdg == PdgCode::kEta) {
321+
if (grmotherPdg == Pdg::kEta) {
324322
mpdg = grmotherPdg;
325323
mpt = grmotherPt;
326324
return GammaEta;
327325
}
328-
if (grmotherPdg == PdgCode::kOmega) {
326+
if (grmotherPdg == Pdg::kOmega) {
329327
mpdg = grmotherPdg;
330328
mpt = grmotherPt;
331329
return GammaOmega;
332330
}
333-
if (grmotherPdg == kPhi) {
331+
if (grmotherPdg == Pdg::kPhi) {
334332
mpdg = grmotherPdg;
335333
mpt = grmotherPt;
336334
return GammaPhi;
337335
}
338-
if (grmotherPdg == PdgCode::kEtaPrime) {
336+
if (grmotherPdg == Pdg::kEtaPrime) {
339337
mpdg = grmotherPdg;
340338
mpt = grmotherPt;
341339
return GammaEtaPrime;
@@ -355,10 +353,10 @@ struct HfTaskSingleElectron {
355353
if (grmothersIdsVec.empty()) {
356354
static const std::map<int, SourceType> pdgToSource = {
357355
{kPi0, Pi0},
358-
{PdgCode::kEta, Eta},
359-
{PdgCode::kOmega, Omega},
360-
{kPhi, Phi},
361-
{PdgCode::kEtaPrime, EtaPrime},
356+
{Pdg::kEta, Eta},
357+
{Pdg::kOmega, Omega},
358+
{Pdg::kPhi, Phi},
359+
{Pdg::kEtaPrime, EtaPrime},
362360
{kRho770_0, Rho0},
363361
{kKPlus, Ke3},
364362
{kK0Long, K0l}};
@@ -393,16 +391,16 @@ struct HfTaskSingleElectron {
393391
mpdg = motherPdg;
394392
return Pi0;
395393
}
396-
if (motherPdg == PdgCode::kEta) {
394+
if (motherPdg == Pdg::kEta) {
397395
return Eta;
398396
}
399-
if (motherPdg == PdgCode::kOmega) {
397+
if (motherPdg == Pdg::kOmega) {
400398
return Omega;
401399
}
402-
if (motherPdg == kPhi) {
400+
if (motherPdg == Pdg::kPhi) {
403401
return Phi;
404402
}
405-
if (motherPdg == PdgCode::kEtaPrime) {
403+
if (motherPdg == Pdg::kEtaPrime) {
406404
return EtaPrime;
407405
}
408406
if (motherPdg == kRho770_0) {
@@ -434,8 +432,8 @@ struct HfTaskSingleElectron {
434432
return;
435433
}
436434

437-
histos.fill(HIST("VtxZ"), collision.posZ());
438-
histos.fill(HIST("nEvents"), flagAnalysedEvt);
435+
histos.fill(HIST("hVtxZ"), collision.posZ());
436+
histos.fill(HIST("hNEvents"), flagAnalysedEvt);
439437

440438
for (const auto& track : tracks) {
441439

@@ -447,34 +445,40 @@ struct HfTaskSingleElectron {
447445
continue;
448446
}
449447

450-
histos.fill(HIST("etaTrack"), track.eta());
451-
histos.fill(HIST("ptTrack"), track.pt());
448+
histos.fill(HIST("hEtaTrack"), track.eta());
449+
histos.fill(HIST("hPtTrack"), track.pt());
452450

453-
histos.fill(HIST("tpcNClsTrack"), track.tpcNClsCrossedRows());
454-
histos.fill(HIST("tpcFoundFindableTrack"), track.tpcCrossedRowsOverFindableCls());
455-
histos.fill(HIST("tpcChi2Track"), track.tpcChi2NCl());
456-
histos.fill(HIST("itsIBClsTrack"), track.itsNClsInnerBarrel());
457-
histos.fill(HIST("itsChi2Track"), track.itsChi2NCl());
458-
histos.fill(HIST("dcaXYTrack"), track.dcaXY());
459-
histos.fill(HIST("dcaZTrack"), track.dcaZ());
451+
histos.fill(HIST("hNClsTpcTrack"), track.tpcNClsCrossedRows());
452+
histos.fill(HIST("hNClsFoundFindableTpcTrack"), track.tpcCrossedRowsOverFindableCls());
453+
histos.fill(HIST("hChi2TpcTrack"), track.tpcChi2NCl());
454+
histos.fill(HIST("hibClsItsTrack"), track.itsNClsInnerBarrel());
455+
histos.fill(HIST("hChi2ItsTrack"), track.itsChi2NCl());
456+
histos.fill(HIST("hDcaXYTrack"), track.dcaXY());
457+
histos.fill(HIST("hDcaZTrack"), track.dcaZ());
460458

461-
histos.fill(HIST("tofNSigPt"), track.pt(), track.tofNSigmaEl());
462-
histos.fill(HIST("tpcNSigP"), track.p(), track.tpcNSigmaEl());
463-
histos.fill(HIST("tpcNSigPt"), track.pt(), track.tpcNSigmaEl());
459+
histos.fill(HIST("hTofNSigPt"), track.pt(), track.tofNSigmaEl());
460+
histos.fill(HIST("hTPCNSigP"), track.p(), track.tpcNSigmaEl());
461+
histos.fill(HIST("hTPCNSigPt"), track.pt(), track.tpcNSigmaEl());
464462

465-
if (std::abs(track.tofNSigmaEl()) > tofNSigmaMax) {
463+
if (std::abs(track.tofNSigmaEl()) > nSigmaTofMax) {
466464
continue;
467465
}
468-
histos.fill(HIST("tofNSigPtQA"), track.pt(), track.tofNSigmaEl());
469-
histos.fill(HIST("tpcNSigPAfterTofCut"), track.p(), track.tpcNSigmaEl());
470-
histos.fill(HIST("tpcNSigPtAfterTofCut"), track.pt(), track.tpcNSigmaEl());
466+
histos.fill(HIST("hTofNSigPtQA"), track.pt(), track.tofNSigmaEl());
467+
histos.fill(HIST("hTPCNSigPAfterTofCut"), track.p(), track.tpcNSigmaEl());
468+
histos.fill(HIST("hTPCNSigPtAfterTofCut"), track.pt(), track.tpcNSigmaEl());
471469

472-
if (track.tpcNSigmaEl() < tpcNSigmaMin || track.tpcNSigmaEl() > tpcNSigmaMax) {
470+
if (track.tpcNSigmaEl() < nSigmaTpcMin || track.tpcNSigmaEl() > nSigmaTPCMax) {
473471
continue;
474472
}
475-
histos.fill(HIST("tpcNSigPtQA"), track.pt(), track.tpcNSigmaEl());
476473

477-
histos.fill(HIST("dcaTrack"), track.pt(), track.dcaXY());
474+
if (track.tpcNSigmaEl() < nSigmaTpcHadronMax && track.tpcNSigmaEl() > nSigmaTpcHadronMin) {
475+
476+
histos.fill(HIST("hDcaHadron"), track.pt(), track.dcaXY());
477+
}
478+
479+
histos.fill(HIST("hTPCNSigPtQA"), track.pt(), track.tpcNSigmaEl());
480+
481+
histos.fill(HIST("hDcaTrack"), track.pt(), track.dcaXY());
478482
}
479483
}
480484
PROCESS_SWITCH(HfTaskSingleElectron, processData, "For real data", true);
@@ -493,64 +497,67 @@ struct HfTaskSingleElectron {
493497
return;
494498
}
495499

496-
histos.fill(HIST("VtxZ"), collision.posZ());
497-
histos.fill(HIST("nEvents"), flagAnalysedEvt);
500+
histos.fill(HIST("hVtxZ"), collision.posZ());
501+
histos.fill(HIST("hNEvents"), flagAnalysedEvt);
498502

499503
for (const auto& track : tracks) {
500504

501505
if (!trackSel(track)) {
502506
continue;
503507
}
504508

505-
histos.fill(HIST("etaTrack"), track.eta());
506-
histos.fill(HIST("ptTrack"), track.pt());
509+
histos.fill(HIST("hEtaTrack"), track.eta());
510+
histos.fill(HIST("hPtTrack"), track.pt());
507511

508-
histos.fill(HIST("tpcNClsTrack"), track.tpcNClsCrossedRows());
509-
histos.fill(HIST("tpcFoundFindableTrack"), track.tpcCrossedRowsOverFindableCls());
510-
histos.fill(HIST("tpcChi2Track"), track.tpcChi2NCl());
511-
histos.fill(HIST("itsIBClsTrack"), track.itsNClsInnerBarrel());
512-
histos.fill(HIST("dcaXYTrack"), track.dcaXY());
513-
histos.fill(HIST("dcaZTrack"), track.dcaZ());
512+
histos.fill(HIST("hNClsTpcTrack"), track.tpcNClsCrossedRows());
513+
histos.fill(HIST("hNClsFoundFindableTpcTrack"), track.tpcCrossedRowsOverFindableCls());
514+
histos.fill(HIST("hChi2TpcTrack"), track.tpcChi2NCl());
515+
histos.fill(HIST("hibClsItsTrack"), track.itsNClsInnerBarrel());
516+
histos.fill(HIST("hDcaXYTrack"), track.dcaXY());
517+
histos.fill(HIST("hDcaZTrack"), track.dcaZ());
514518

515-
histos.fill(HIST("tofNSigPt"), track.pt(), track.tofNSigmaEl());
516-
histos.fill(HIST("tpcNSigPt"), track.pt(), track.tpcNSigmaEl());
519+
histos.fill(HIST("hTofNSigPt"), track.pt(), track.tofNSigmaEl());
520+
histos.fill(HIST("hTPCNSigPt"), track.pt(), track.tpcNSigmaEl());
517521

518522
int mpdg; // electron source pdg code
519523
double mpt; // electron source pt
520524
int const source = getElecSource(track, mpt, mpdg);
521525

522526
if (source == DirectBeauty || source == BeautyCharm) {
523527
histos.fill(HIST("hPdgB"), mpdg);
524-
histos.fill(HIST("dcaBeauty"), track.pt(), track.dcaXY());
528+
histos.fill(HIST("hDcaBeauty"), track.pt(), track.dcaXY());
525529
}
526530

527531
if (source == DirectCharm) {
528532
histos.fill(HIST("hPdgC"), mpdg);
529-
histos.fill(HIST("dcaCharm"), track.pt(), track.dcaXY());
533+
histos.fill(HIST("hDcaChram"), track.pt(), track.dcaXY());
530534
}
531535

532536
if (source >= GammaPi0 && source <= GammaSigma) {
533537
histos.fill(HIST("hPdgCo"), mpdg);
534-
histos.fill(HIST("dcaConv"), track.pt(), track.dcaXY());
538+
histos.fill(HIST("hDcaConv"), track.pt(), track.dcaXY());
535539
}
536540

537541
if (source >= Pi0 && source <= Sigma) {
538542
histos.fill(HIST("hPdgDa"), mpdg);
539-
histos.fill(HIST("dcaDalitz"), track.pt(), track.dcaXY());
543+
histos.fill(HIST("hDcaDalitz"), track.pt(), track.dcaXY());
540544
}
541545

542-
if (std::abs(track.tofNSigmaEl()) > tofNSigmaMax) {
546+
if (track.tpcNSigmaEl() < nSigmaTpcHadronMax && track.tpcNSigmaEl() > nSigmaTpcHadronMin)
547+
histos.fill(HIST("hDcaHadron"), track.pt(), track.dcaXY());
548+
549+
if (std::abs(track.tofNSigmaEl()) > nSigmaTofMax) {
543550
continue;
544551
}
545-
histos.fill(HIST("tofNSigPtQA"), track.pt(), track.tofNSigmaEl());
546-
histos.fill(HIST("tpcNSigPtAfterTofCut"), track.pt(), track.tpcNSigmaEl());
552+
histos.fill(HIST("hTofNSigPtQA"), track.pt(), track.tofNSigmaEl());
553+
histos.fill(HIST("hTPCNSigPtAfterTofCut"), track.pt(), track.tpcNSigmaEl());
547554

548-
if (track.tpcNSigmaEl() < tpcNSigmaMin || track.tpcNSigmaEl() > tpcNSigmaMax) {
555+
if (track.tpcNSigmaEl() < nSigmaTpcMin || track.tpcNSigmaEl() > nSigmaTPCMax) {
549556
continue;
550557
}
551-
histos.fill(HIST("tpcNSigPtQA"), track.pt(), track.tpcNSigmaEl());
558+
histos.fill(HIST("hTPCNSigPtQA"), track.pt(), track.tpcNSigmaEl());
552559

553-
histos.fill(HIST("dcaTrack"), track.pt(), track.dcaXY());
560+
histos.fill(HIST("hDcaTrack"), track.pt(), track.dcaXY());
554561
}
555562
}
556563
PROCESS_SWITCH(HfTaskSingleElectron, processMc, "For real data", false);

0 commit comments

Comments
 (0)