Skip to content

Commit 1b41ed7

Browse files
KotliarovArKotliarovAr
committed
[PWGJE]: add counting of TT for restricted phi range
1 parent 01d19d9 commit 1b41ed7

1 file changed

Lines changed: 72 additions & 1 deletion

File tree

‎PWGJE/Tasks/recoilJets.cxx‎

Lines changed: 72 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -355,7 +355,7 @@ struct RecoilJets {
355355

356356
auto tmpHistPointer = spectra.add<TH2>(Form("hScaled%s_Ntrig", eaAxis.label),
357357
Form("Scaled %s vs Total number of selected triggers per class", eaAxis.label),
358-
kTH2F, {{eaAxis.axis, eaAxis.axisName}, {2, 0.0, 2.}});
358+
kTH2F, {{eaAxis.axis, eaAxis.axisName}, {2, 0.0, 2.}}, hist.sumw2);
359359
tmpHistPointer->GetYaxis()->SetBinLabel(1, "TT_{Ref}");
360360
tmpHistPointer->GetYaxis()->SetBinLabel(2, "TT_{Sig}");
361361

@@ -411,6 +411,19 @@ struct RecoilJets {
411411
Form("Events w. TT_{Sig}: scaled %s & #rho", eaAxis.label),
412412
kTH2F, {{eaAxis.axis, eaAxis.axisName}, rho}, hist.sumw2);
413413

414+
// Rectricted phi range for TT selection
415+
spectra.add(Form("hScaled%s_Ntrig_RestrictedPhi", eaAxis.label),
416+
Form("Scaled %s vs Total number of selected triggers per class #in #varphi (%.2f, %.2f)", eaAxis.label, phiMin, phiMax),
417+
kTH2F, {{eaAxis.axis, eaAxis.axisName}, {2, 0.0, 2.}}, hist.sumw2);
418+
419+
spectra.add(Form("hScaled%s_Recoil_JetPt_Corr_TTRef_RestrictedPhi", eaAxis.label),
420+
Form("Events w. TT_{Ref}: scaled %s & #it{p}_{T} of recoil jets", eaAxis.label),
421+
kTH2F, {{eaAxis.axis, eaAxis.axisName}, jetPTcorr}, hist.sumw2);
422+
423+
spectra.add(Form("hScaled%s_Recoil_JetPt_Corr_TTSig_RestrictedPhi", eaAxis.label),
424+
Form("Events w. TT_{Sig}: scaled %s & #it{p}_{T} of recoil jets", eaAxis.label),
425+
kTH2F, {{eaAxis.axis, eaAxis.axisName}, jetPTcorr}, hist.sumw2);
426+
414427
spectra.add(Form("hScaled%s_DPhi_JetPt_Corr_TTRef_RestrictedPhi", eaAxis.label),
415428
Form("Events w. TT_{Ref} #in #varphi (%.2f, %.2f): scaled %s & #Delta#varphi & #it{p}_{T, jet}^{ch}", phiMin, phiMax, eaAxis.label),
416429
kTH3F, {{eaAxis.axis, eaAxis.axisName}, deltaPhiAngle, jetPTcorr}, hist.sumw2);
@@ -580,6 +593,18 @@ struct RecoilJets {
580593
kTH2F, {{eaAxis.axis, eaAxis.axisName}, rho}, hist.sumw2);
581594

582595
// Rectricted phi range for TT selection
596+
spectra.add(Form("hScaled%s_Ntrig_RestrictedPhi_Part", eaAxis.label),
597+
Form("Scaled %s vs Total number of selected triggers per class #in #varphi (%.2f, %.2f)", eaAxis.label, phiMin, phiMax),
598+
kTH2F, {{eaAxis.axis, eaAxis.axisName}, {2, 0.0, 2.}}, hist.sumw2);
599+
600+
spectra.add(Form("hScaled%s_Recoil_JetPt_Corr_TTRef_RestrictedPhi_Part", eaAxis.label),
601+
Form("Events w. TT_{Ref}: scaled %s & #it{p}_{T} of recoil jets", eaAxis.label),
602+
kTH2F, {{eaAxis.axis, eaAxis.axisName}, jetPTcorr}, hist.sumw2);
603+
604+
spectra.add(Form("hScaled%s_Recoil_JetPt_Corr_TTSig_RestrictedPhi_Part", eaAxis.label),
605+
Form("Events w. TT_{Sig}: scaled %s & #it{p}_{T} of recoil jets", eaAxis.label),
606+
kTH2F, {{eaAxis.axis, eaAxis.axisName}, jetPTcorr}, hist.sumw2);
607+
583608
spectra.add(Form("hScaled%s_DPhi_JetPt_Corr_TTRef_RestrictedPhi_Part", eaAxis.label),
584609
Form("MC events w. TT_{Ref} #in #varphi (%.2f, %.2f): scaled %s & #Delta#varphi & #it{p}_{T, jet}^{ch}", phiMin, phiMax, eaAxis.label),
585610
kTH3F, {{eaAxis.axis, eaAxis.axisName}, deltaPhiAngle, jetPTcorr}, hist.sumw2);
@@ -1327,11 +1352,19 @@ struct RecoilJets {
13271352

13281353
phiTT = getPhiTT(vPhiOfTT);
13291354

1355+
const auto phiMin = tt.phiRestr->at(0);
1356+
const auto phiMax = tt.phiRestr->at(1);
1357+
13301358
if (bSigEv) {
13311359
// EA
13321360
spectra.fill(HIST("hScaledFT0C_Ntrig"), scaledFT0C, addCountToTTSig, weight);
13331361
spectra.fill(HIST("hScaledFT0M_Ntrig"), scaledFT0M, addCountToTTSig, weight);
13341362

1363+
if (phiTT > phiMin && phiTT < phiMax) {
1364+
spectra.fill(HIST("hScaledFT0C_Ntrig_RestrictedPhi"), scaledFT0C, addCountToTTSig, weight);
1365+
spectra.fill(HIST("hScaledFT0M_Ntrig_RestrictedPhi"), scaledFT0M, addCountToTTSig, weight);
1366+
}
1367+
13351368
spectra.fill(HIST("hScaledFT0C_TTSig_per_event"), scaledFT0C, nTT, weight);
13361369
spectra.fill(HIST("hScaledFT0M_TTSig_per_event"), scaledFT0M, nTT, weight);
13371370

@@ -1362,6 +1395,11 @@ struct RecoilJets {
13621395
spectra.fill(HIST("hScaledFT0C_Ntrig"), scaledFT0C, addCountToTTRef, weight);
13631396
spectra.fill(HIST("hScaledFT0M_Ntrig"), scaledFT0M, addCountToTTRef, weight);
13641397

1398+
if (phiTT > phiMin && phiTT < phiMax) {
1399+
spectra.fill(HIST("hScaledFT0C_Ntrig_RestrictedPhi"), scaledFT0C, addCountToTTRef, weight);
1400+
spectra.fill(HIST("hScaledFT0M_Ntrig_RestrictedPhi"), scaledFT0M, addCountToTTRef, weight);
1401+
}
1402+
13651403
spectra.fill(HIST("hScaledFT0C_TTRef_per_event"), scaledFT0C, nTT, weight);
13661404
spectra.fill(HIST("hScaledFT0M_TTRef_per_event"), scaledFT0M, nTT, weight);
13671405

@@ -1441,6 +1479,11 @@ struct RecoilJets {
14411479
spectra.fill(HIST("hCentFT0C_Recoil_JetPt_Corr_TTSig"), centFT0C, jetPtCorr, weight);
14421480
spectra.fill(HIST("hCentFT0M_Recoil_JetPt_Corr_TTSig"), centFT0M, jetPtCorr, weight);
14431481
spectra.fill(HIST("hCentFT0CVar1_Recoil_JetPt_Corr_TTSig"), centFT0CVar1, jetPtCorr, weight);
1482+
1483+
if (phiTT > phiMin && phiTT < phiMax) {
1484+
spectra.fill(HIST("hScaledFT0C_Recoil_JetPt_Corr_TTSig_RestrictedPhi"), scaledFT0C, jetPtCorr, weight);
1485+
spectra.fill(HIST("hScaledFT0M_Recoil_JetPt_Corr_TTSig_RestrictedPhi"), scaledFT0M, jetPtCorr, weight);
1486+
}
14441487
}
14451488
} else {
14461489

@@ -1476,6 +1519,11 @@ struct RecoilJets {
14761519
spectra.fill(HIST("hCentFT0C_Recoil_JetPt_Corr_TTRef"), centFT0C, jetPtCorr, weight);
14771520
spectra.fill(HIST("hCentFT0M_Recoil_JetPt_Corr_TTRef"), centFT0M, jetPtCorr, weight);
14781521
spectra.fill(HIST("hCentFT0CVar1_Recoil_JetPt_Corr_TTRef"), centFT0CVar1, jetPtCorr, weight);
1522+
1523+
if (phiTT > phiMin && phiTT < phiMax) {
1524+
spectra.fill(HIST("hScaledFT0C_Recoil_JetPt_Corr_TTRef_RestrictedPhi"), scaledFT0C, jetPtCorr, weight);
1525+
spectra.fill(HIST("hScaledFT0M_Recoil_JetPt_Corr_TTRef_RestrictedPhi"), scaledFT0M, jetPtCorr, weight);
1526+
}
14791527
}
14801528
}
14811529
}
@@ -1566,12 +1614,20 @@ struct RecoilJets {
15661614

15671615
phiTT = getPhiTT(vPhiOfTT);
15681616

1617+
const auto phiMin = tt.phiRestr->at(0);
1618+
const auto phiMax = tt.phiRestr->at(1);
1619+
15691620
if (bSigEv) {
15701621

15711622
// EA
15721623
spectra.fill(HIST("hScaledFT0C_Ntrig_Part"), scaledFT0C, addCountToTTSig, weight);
15731624
spectra.fill(HIST("hScaledFT0M_Ntrig_Part"), scaledFT0M, addCountToTTSig, weight);
15741625

1626+
if (phiTT > phiMin && phiTT < phiMax) {
1627+
spectra.fill(HIST("hScaledFT0C_Ntrig_RestrictedPhi_Part"), scaledFT0C, addCountToTTSig, weight);
1628+
spectra.fill(HIST("hScaledFT0M_Ntrig_RestrictedPhi_Part"), scaledFT0M, addCountToTTSig, weight);
1629+
}
1630+
15751631
spectra.fill(HIST("hScaledFT0C_TTSig_per_event_Part"), scaledFT0C, nTT, weight);
15761632
spectra.fill(HIST("hScaledFT0M_TTSig_per_event_Part"), scaledFT0M, nTT, weight);
15771633

@@ -1599,6 +1655,11 @@ struct RecoilJets {
15991655
spectra.fill(HIST("hScaledFT0C_Ntrig_Part"), scaledFT0C, addCountToTTRef, weight);
16001656
spectra.fill(HIST("hScaledFT0M_Ntrig_Part"), scaledFT0M, addCountToTTRef, weight);
16011657

1658+
if (phiTT > phiMin && phiTT < phiMax) {
1659+
spectra.fill(HIST("hScaledFT0C_Ntrig_RestrictedPhi_Part"), scaledFT0C, addCountToTTRef, weight);
1660+
spectra.fill(HIST("hScaledFT0M_Ntrig_RestrictedPhi_Part"), scaledFT0M, addCountToTTRef, weight);
1661+
}
1662+
16021663
spectra.fill(HIST("hScaledFT0C_TTRef_per_event_Part"), scaledFT0C, nTT, weight);
16031664
spectra.fill(HIST("hScaledFT0M_TTRef_per_event_Part"), scaledFT0M, nTT, weight);
16041665

@@ -1669,6 +1730,11 @@ struct RecoilJets {
16691730
spectra.fill(HIST("hCentFT0A_Recoil_JetPt_Corr_TTSig_Part"), centFT0A, jetPtCorr, weight);
16701731
spectra.fill(HIST("hCentFT0C_Recoil_JetPt_Corr_TTSig_Part"), centFT0C, jetPtCorr, weight);
16711732
spectra.fill(HIST("hCentFT0M_Recoil_JetPt_Corr_TTSig_Part"), centFT0M, jetPtCorr, weight);
1733+
1734+
if (phiTT > phiMin && phiTT < phiMax) {
1735+
spectra.fill(HIST("hScaledFT0C_Recoil_JetPt_Corr_TTSig_RestrictedPhi_Part"), scaledFT0C, jetPtCorr, weight);
1736+
spectra.fill(HIST("hScaledFT0M_Recoil_JetPt_Corr_TTSig_RestrictedPhi_Part"), scaledFT0M, jetPtCorr, weight);
1737+
}
16721738
}
16731739
} else {
16741740

@@ -1702,6 +1768,11 @@ struct RecoilJets {
17021768
spectra.fill(HIST("hCentFT0A_Recoil_JetPt_Corr_TTRef_Part"), centFT0A, jetPtCorr, weight);
17031769
spectra.fill(HIST("hCentFT0C_Recoil_JetPt_Corr_TTRef_Part"), centFT0C, jetPtCorr, weight);
17041770
spectra.fill(HIST("hCentFT0M_Recoil_JetPt_Corr_TTRef_Part"), centFT0M, jetPtCorr, weight);
1771+
1772+
if (phiTT > phiMin && phiTT < phiMax) {
1773+
spectra.fill(HIST("hScaledFT0C_Recoil_JetPt_Corr_TTRef_RestrictedPhi_Part"), scaledFT0C, jetPtCorr, weight);
1774+
spectra.fill(HIST("hScaledFT0M_Recoil_JetPt_Corr_TTRef_RestrictedPhi_Part"), scaledFT0M, jetPtCorr, weight);
1775+
}
17051776
}
17061777
}
17071778
}

0 commit comments

Comments
 (0)