Skip to content

Commit d3aa8cf

Browse files
sashingosashingo
committed
implemeqted comments and suggestions
1 parent ec18869 commit d3aa8cf

File tree

1 file changed

+44
-23
lines changed

1 file changed

+44
-23
lines changed

‎PWGHF/HFL/Tasks/taskElectronWeakBoson.cxx‎

Lines changed: 44 additions & 23 deletions
Original file line numberDiff line numberDiff line change
@@ -69,8 +69,29 @@ struct HfTaskElectronWeakBoson {
6969
// configurable parameters
7070
Configurable<int> nBinsPt{"nBinsPt", 100, "N bins in pt registry"};
7171
Configurable<float> binPtmax{"binPtmax", 100.0, "maximum pt registry"};
72+
Configurable<int> nBinsPtZee{"nBinsPtZee", 60, "N bins in pt Zee registry"};
73+
Configurable<float> binPtZeemin{"binPtZeemin", 20.0, "minimum pt Zee registry"};
74+
Configurable<float> binPtZeemax{"binPtZeemax", 80.0, "maximum pt Zee registry"};
7275
Configurable<int> nBinsE{"nBinsE", 100, "N bins in E registry"};
7376
Configurable<float> binEmax{"binEmax", 100.0, "maximum E registry"};
77+
Configurable<int> nBinsEop{"nBinsEop", 300, "N bins in E/p registry"};
78+
Configurable<float> binEopmin{"binEopmin", -0.01, "minimum E/p registry"};
79+
Configurable<float> binEopmax{"binEopmax", 1.49, "maximum E/p registry"};
80+
Configurable<int> nBinsEiso{"nBinsEiso", 255, "N bins in Eiso registry"};
81+
Configurable<float> binEisomin{"binEisomin", 0, "minimum Eiso registry"};
82+
Configurable<float> binEisomax{"binEisomax", 2.0, "maximum Eiso registry"};
83+
Configurable<int> nBinsPiso{"nBinsPiso", 255, "N bins in mom. isolation registry"};
84+
Configurable<float> binPisomin{"binPisomin", 0, "minimum mom. isolation registry"};
85+
Configurable<float> binPisomax{"binPisomax", 1.5, "maximum mom. isolation registry"};
86+
Configurable<int> nBinsTriso{"nBinsTriso", 25, "N bins in Track isolation registry"};
87+
Configurable<float> binTrisomin{"binTrisomin", -0.5, "minimum Track isolation registry"};
88+
Configurable<float> binTrisomax{"binTrisomax", 24.5, "maximum Track isolation registry"};
89+
Configurable<int> nBinsMassZg{"nBinsMassZg", 150, "N bins in Mass Zgamma registry"};
90+
Configurable<float> binMassZgmax{"binMassZgmax", 150, "maximum Mass Zgamma registry"};
91+
Configurable<int> nBinsMassZ{"nBinsMassZ", 130, "N bins in Mass Z registry"};
92+
Configurable<float> binMassZmin{"binMassZmin", 20, "minimum Mass Z registry"};
93+
Configurable<int> nBinsZfrag{"nBinsZfrag", 200, "N bins in Z frag registry"};
94+
Configurable<float> binZfragmin{"binZfragmin", 2, "maximum Z frag registry"};
7495

7596
Configurable<float> vtxZ{"vtxZ", 10.f, ""};
7697

@@ -214,9 +235,9 @@ struct HfTaskElectronWeakBoson {
214235
const AxisSpec axisCounter{1, 0, 1, "events"};
215236
const AxisSpec axisEta{20, -1.0, 1.0, "#eta"};
216237
const AxisSpec axisPt{nBinsPt, 0, binPtmax, "p_{T}"};
217-
const AxisSpec axisPtZee{60, 20, 80, "p_{T}"};
218-
const AxisSpec axisPtZneg{60, 20, 80, "p_{T,neg} (GeV/c)"};
219-
const AxisSpec axisPtZpos{60, 20, 80, "p_{T,pos} (GeV/c)"};
238+
const AxisSpec axisPtZee{nBinsPtZee, binPtZeemin, binPtZeemax, "p_{T}"};
239+
const AxisSpec axisPtZneg{nBinsPtZee, binPtZeemin, binPtZeemax, "p_{T,neg} (GeV/c)"};
240+
const AxisSpec axisPtZpos{nBinsPtZee, binPtZeemin, binPtZeemax, "p_{T,pos} (GeV/c)"};
220241
const AxisSpec axisDCAxyneg{150, 0, 0.3, "DCAxy_{neg}"};
221242
const AxisSpec axisDCAxypos{150, 0, 0.3, "DCAxy_{pos}"};
222243
const AxisSpec axisDCAzneg{150, 0, 0.3, "DCAz_{neg}"};
@@ -236,31 +257,31 @@ struct HfTaskElectronWeakBoson {
236257
const AxisSpec axisdR{20, 0.0, 0.2, "dR"};
237258
const AxisSpec axisNcell{50, 0.0, 50.0, "Ncell"};
238259
const AxisSpec axisPhi{350, 0, 7, "Phi"};
239-
const AxisSpec axisEop{200, 0, 2, "E/p"};
240-
const AxisSpec axisEopZneg{300, -0.01, 1.49, "E/p neg"};
241-
const AxisSpec axisEopZpos{300, -0.01, 1.49, "E/p pos"};
260+
const AxisSpec axisEop{nBinsEop, binEopmin, binEopmax, "E/p"};
261+
const AxisSpec axisEopZneg{nBinsEop, binEopmin, binEopmax, "E/p"};
262+
const AxisSpec axisEopZpos{nBinsEop, binEopmin, binEopmax, "E/p"};
242263
const AxisSpec axisChi2{250, 0.0, 25.0, "#chi^{2}"};
243264
const AxisSpec axisCluster{100, 0.0, 200.0, "counts"};
244265
const AxisSpec axisITSNCls{10, 0.0, 10, "counts"};
245266
const AxisSpec axisEMCtime{100, -50.0, 50, "EMC time"};
246-
const AxisSpec axisIsoEnergy{100, 0, 1.0, "E_{iso}"};
247-
const AxisSpec axisIsoEnergyZneg{255, 0, 1.5, "E_{iso,neg}"};
248-
const AxisSpec axisIsoEnergyZpos{255, 0, 1.5, "E_{iso,pos}"};
249-
const AxisSpec axisIsoMomentum{100, 0, 10.0, "Isolation momentum(GeV/C)"};
250-
const AxisSpec axisIsoMomentumZneg{255, 0, 2.0, "p_{iso,neg}"};
251-
const AxisSpec axisIsoMomentumZpos{255, 0, 2.0, "p_{iso,pos}"};
252-
const AxisSpec axisIsoTrack{25, -0.5, 24.5, "Isolation Track"};
253-
const AxisSpec axisIsoTrackZneg{25, -0.5, 24.5, "N_{isotrk,neg}"};
254-
const AxisSpec axisIsoTrackZpos{25, -0.5, 24.5, "N_{isotrk,pos}"};
255-
const AxisSpec axisInvMassZgamma{150, 0, 150, "M_{ee} (GeV/c^{2})"};
256-
const AxisSpec axisInvMassZ{130, 20, 150, "M_{ee} (GeV/c^{2})"};
267+
const AxisSpec axisIsoEnergy{nBinsEiso, binEisomin, binEisomax, "E_{iso}"};
268+
const AxisSpec axisIsoEnergyZneg{nBinsEiso, binEisomin, binEisomax, "E_{iso}"};
269+
const AxisSpec axisIsoEnergyZpos{nBinsEiso, binEisomin, binEisomax, "E_{iso}"};
270+
const AxisSpec axisIsoMomentum{nBinsPiso, binPisomin, binPisomax, "Isolation momentum(GeV/C)"};
271+
const AxisSpec axisIsoMomentumZneg{nBinsPiso, binPisomin, binPisomax, "Isolation momentum(GeV/C)"};
272+
const AxisSpec axisIsoMomentumZpos{nBinsPiso, binPisomin, binPisomax, "Isolation momentum(GeV/C)"};
273+
const AxisSpec axisIsoTrack{nBinsTriso, binTrisomin, binTrisomax, "Isolation Track"};
274+
const AxisSpec axisIsoTrackZneg{nBinsTriso, binTrisomin, binTrisomax, "N_{isotrk,neg}"};
275+
const AxisSpec axisIsoTrackZpos{nBinsTriso, binTrisomin, binTrisomax, "N_{isotrk,pos}"};
276+
const AxisSpec axisInvMassZgamma{nBinsMassZg, 0, binMassZgmax, "M_{ee} (GeV/c^{2})"};
277+
const AxisSpec axisInvMassZ{nBinsMassZ, binMassZmin, binMassZgmax, "M_{ee} (GeV/c^{2})"};
257278
const AxisSpec axisTrigger{3, -0.5, 2.5, "Trigger status of zorro"};
258279
const AxisSpec axisDPhiZh{64, -o2::constants::math::PIHalf, 3 * o2::constants::math::PIHalf, "#Delta#phi(Z-h)"};
259-
const AxisSpec axisPtHadron{50, 0, 50, "p_{T,hadron} (GeV/c)"};
260-
const AxisSpec axisPtZ{150, 0, 150, "p_{T,Z} (GeV/c)"};
280+
const AxisSpec axisPtHadron{nBinsPt, 0, binPtmax, "p_{T,hadron} (GeV/c)"};
281+
const AxisSpec axisPtZ{nBinsPt, 0, binPtmax, "p_{T,Z} (GeV/c)"};
261282
const AxisSpec axisSign{2, -2, 2, "charge sign"};
262283
const AxisSpec axisCentrality{centralityBins};
263-
const AxisSpec axisPtRatio{200, 0, 2.0, "pt ratio for h and Z"};
284+
const AxisSpec axisZfrag{nBinsZfrag, 0, binZfragmin, "p_{T,h}/p_{T,Z}"};
264285

265286
// create registrygrams
266287
registry.add("hZvtx", "Z vertex", kTH1D, {axisZvtx});
@@ -295,7 +316,7 @@ struct HfTaskElectronWeakBoson {
295316
registry.add("hTHnTrMatch", "Track EMC Match", HistType::kTHnSparseF, {axisPt, axisdPhi, axisdEta});
296317

297318
// Z-hadron correlation histograms
298-
registry.add("hZHadronDphi", "Z-hadron #Delta#phi correlation", HistType::kTHnSparseF, {axisCentrality, axisSign, axisPtZ, axisDPhiZh, axisPtRatio, axisPtHadron});
319+
registry.add("hZHadronDphi", "Z-hadron #Delta#phi correlation", HistType::kTHnSparseF, {axisCentrality, axisSign, axisPtZ, axisDPhiZh, axisZfrag, axisPtHadron});
299320
registry.add("hZptSpectrum", "Z boson p_{T} spectrum", kTH2F, {{axisSign}, {axisPtZ}});
300321

301322
// hisotgram for EMCal trigger
@@ -751,7 +772,7 @@ struct HfTaskElectronWeakBoson {
751772
// bg e-e-
752773
for (const auto& trackEle : selectedElectronsIso) {
753774
for (const auto& trackEle2 : selectedElectronsIso) {
754-
if (trackEle.pt == trackEle2.pt) {
775+
if (&trackEle.pt == &trackEle2.pt) {
755776
continue;
756777
}
757778
auto child1 = RecoDecayPtEtaPhi::pVector(trackEle.pt, trackEle.eta, trackEle.phi);
@@ -767,7 +788,7 @@ struct HfTaskElectronWeakBoson {
767788
// bg e+e+
768789
for (const auto& trackPos : selectedPositronsIso) {
769790
for (const auto& trackPos2 : selectedPositronsIso) {
770-
if (trackPos.pt == trackPos2.pt) {
791+
if (&trackPos.pt == &trackPos2.pt) {
771792
continue;
772793
}
773794
auto child1 = RecoDecayPtEtaPhi::pVector(trackPos.pt, trackPos.eta, trackPos.phi);

0 commit comments

Comments
 (0)