Skip to content

Commit 1034dba

Browse files
author
jimun_lee
committed
[PWGLF] Added the option of cfgJetdR
1 parent 2f8682c commit 1034dba

1 file changed

Lines changed: 30 additions & 14 deletions

File tree

‎PWGLF/Tasks/Resonances/kstarInOO.cxx‎

Lines changed: 30 additions & 14 deletions
Original file line numberDiff line numberDiff line change
@@ -149,7 +149,8 @@ struct kstarInOO {
149149
//|
150150
//======================
151151
Configurable<float> cfgJetpT{"cfgJetpT", 8.0, "Set Jet pT minimum"};
152-
Configurable<float> cfgJetR{"cfgJetR", 0.4, "Set Jet radius parameter"};
152+
Configurable<float> cfgJetR{"cfgJetR", 0.4, "Anti-kT Radius"};
153+
Configurable<float> cfgJetdR{"cfgJetdR", 0.4, "Set Jet radius parameter"};
153154
Configurable<bool> cfgSingleJet{"cfgSingleJet", false, "Enforces strict phi-jet correspondance"};
154155
Configurable<bool> cfgReqJets{"cfgReqJets", false, "False: MB, True: Inside Jets"};
155156
Configurable<std::string> cfgRealTriggerMasks{"cfgRealTriggerMasks", "", "possible JE Trigger masks: fJetChLowPt,fJetChHighPt,fTrackLowPt,fTrackHighPt,fJetD0ChLowPt,fJetD0ChHighPt,fJetLcChLowPt,fJetLcChHighPt,fEMCALReadout,fJetFullHighPt,fJetFullLowPt,fJetNeutralHighPt,fJetNeutralLowPt,fGammaVeryHighPtEMCAL,fGammaVeryHighPtDCAL,fGammaHighPtEMCAL,fGammaHighPtDCAL,fGammaLowPtEMCAL,fGammaLowPtDCAL,fGammaVeryLowPtEMCAL,fGammaVeryLowPtDCAL"};
@@ -278,8 +279,14 @@ struct kstarInOO {
278279
histos.add("missed_kpi_INJets_8_infinite", "missed kpi In Jets with 8 < jetPt < infinite", {HistType::kTH2F, {{120, 0.0, 1.2}, {100, 0., 20.}}});
279280

280281
histos.add("recoveredJetpT_6_8to8_10", "recovered Jet pT", kTH1F, {{2000, 0., 100.}});
282+
histos.add("recoveredJetpT_6_8to8_10_kstarSpectra", "Kstar pT within the recovered Jet pT", kTH1F, {{2000, 0., 100.}});
281283

282-
histos.add("JetMigration", "bin to bin migration", {HistType::kTH2F, {{100, 0.0, 50.0, "True jet pT (GeV/c)"}, {100, 0., 50., "Recovered jet pT (GeV/c)"}}});
284+
histos.add("normalJetpT_8_kstarSpectra", "kstar pT in Jet > 8GeV/c", kTH1F, {{2000, 0., 100.}});
285+
histos.add("normalJetpT_6_8_kstarSpectra", "6 GeV/c < kstar pT in Jet < 8 GeV/c", kTH1F, {{2000, 0., 100.}});
286+
histos.add("normalJetpT_8_10_kstarSpectra", "8 GeV/c < kstar pT in Jet < 10 GeV/c", kTH1F, {{2000, 0., 100.}});
287+
histos.add("normalJetpT_10_12_kstarSpectra", "10 GeV/c < kstar pT in Jet < 12 GeV/c", kTH1F, {{2000, 0., 100.}});
288+
289+
histos.add("JetMigration", "bin to bin migration", {HistType::kTH2F, {{150, 0.0, 15.0, "True jet pT (GeV/c)"}, {150, 0., 15., "Recovered jet pT (GeV/c)"}}});
283290
}
284291

285292
////////////////////////////////////
@@ -1807,7 +1814,7 @@ struct kstarInOO {
18071814
bestJetEta = mcpjet.eta();
18081815
}
18091816
} // mcpJets
1810-
if (bestR > cfgJetR)
1817+
if (bestR > cfgJetdR)
18111818
continue;
18121819

18131820
//==================
@@ -1824,12 +1831,12 @@ struct kstarInOO {
18241831
double deta_kaon = bestJetEta - daughter.eta();
18251832
dR_kaon = TMath::Sqrt((dphi_kaon * dphi_kaon) + (deta_kaon * deta_kaon));
18261833

1827-
if (bestR < cfgJetR) {
1834+
if (bestR < cfgJetdR) {
18281835
if (cfgJetdRHistos) {
18291836
histos.fill(HIST("dR_taggedjet_kaon"), dR_kaon, lResonance.Pt());
18301837
histos.fill(HIST("dR_taggedjet_all"), dR_kaon, lResonance.Pt());
18311838
}
1832-
if (dR_kaon > cfgJetR) {
1839+
if (dR_kaon > cfgJetdR) {
18331840
kaon_out = true;
18341841
missing_pt += daughter.pt();
18351842
}
@@ -1841,15 +1848,15 @@ struct kstarInOO {
18411848
double deta_pion = bestJetEta - daughter.eta();
18421849
dR_pion = TMath::Sqrt((dphi_pion * dphi_pion) + (deta_pion * deta_pion));
18431850

1844-
if (bestR < cfgJetR) {
1851+
if (bestR < cfgJetdR) {
18451852
if (cfgJetdRHistos) {
18461853
histos.fill(HIST("dR_taggedjet_pion"), dR_pion, lResonance.Pt());
18471854
histos.fill(HIST("dR_taggedjet_all"), dR_pion, lResonance.Pt());
18481855

18491856
if (bestJetpT > 6.0 && bestJetpT < 8.0)
18501857
histos.fill(HIST("dR_taggedjet_all_6_8"), dR_pion, lResonance.Pt());
18511858
}
1852-
if (dR_pion > cfgJetR) {
1859+
if (dR_pion > cfgJetdR) {
18531860
pion_out = true;
18541861
missing_pt += daughter.pt();
18551862
}
@@ -1861,7 +1868,7 @@ struct kstarInOO {
18611868
double deta_kaon = bestJetEta - daughter.eta();
18621869
dR_kaon = TMath::Sqrt((dphi_kaon * dphi_kaon) + (deta_kaon * deta_kaon));
18631870

1864-
if (bestR < cfgJetR) {
1871+
if (bestR < cfgJetdR) {
18651872
if (cfgJetdRHistos) {
18661873
histos.fill(HIST("dR_taggedjet_kaon"), dR_kaon, lResonance.Pt());
18671874
histos.fill(HIST("dR_taggedjet_all"), dR_kaon, lResonance.Pt());
@@ -1870,7 +1877,7 @@ struct kstarInOO {
18701877
histos.fill(HIST("dR_taggedjet_all_6_8"), dR_kaon, lResonance.Pt());
18711878
}
18721879

1873-
if (dR_kaon > cfgJetR) {
1880+
if (dR_kaon > cfgJetdR) {
18741881
kaon_out = true;
18751882
missing_pt = daughter.pt();
18761883
}
@@ -1883,24 +1890,33 @@ struct kstarInOO {
18831890
double recoveredJetpT = bestJetpT + missing_pt;
18841891
if (cfgJetdRHistos) {
18851892
if (bestJetpT > 6.0 && bestJetpT < 8.0) {
1893+
histos.fill(HIST("normalJetpT_6_8_kstarSpectra"), lResonance.Pt());
18861894
histos.fill(HIST("missed_kpi_INJets_6_8"), (bestJetpT - missing_pt) / bestJetpT, lResonance.Pt());
1887-
if (recoveredJetpT > 8.0)
1895+
if (recoveredJetpT > 8.0) {
18881896
histos.fill(HIST("recoveredJetpT_6_8to8_10"), recoveredJetpT);
1897+
histos.fill(HIST("recoveredJetpT_6_8to8_10_kstarSpectra"), lResonance.Pt());
1898+
}
18891899
}
1890-
if (bestJetpT > 8.0 && bestJetpT < 10.0)
1900+
1901+
if (bestJetpT > 8.0 && bestJetpT < 10.0) {
18911902
histos.fill(HIST("missed_kpi_INJets_8_10"), (bestJetpT - missing_pt) / bestJetpT, lResonance.Pt());
1892-
if (bestJetpT > 10.0 && bestJetpT < 12.0)
1903+
histos.fill(HIST("normalJetpT_8_10_kstarSpectra"), lResonance.Pt());
1904+
}
1905+
if (bestJetpT > 10.0 && bestJetpT < 12.0) {
18931906
histos.fill(HIST("missed_kpi_INJets_10_12"), (bestJetpT - missing_pt) / bestJetpT, lResonance.Pt());
1907+
histos.fill(HIST("normalJetpT_10_12_kstarSpectra"), lResonance.Pt());
1908+
}
18941909
if (bestJetpT > 12.0 && bestJetpT < 15.0)
18951910
histos.fill(HIST("missed_kpi_INJets_12_15"), (bestJetpT - missing_pt) / bestJetpT, lResonance.Pt());
18961911
if (bestJetpT > 15.0 && bestJetpT < 25.0)
18971912
histos.fill(HIST("missed_kpi_INJets_15_25"), (bestJetpT - missing_pt) / bestJetpT, lResonance.Pt());
18981913
if (bestJetpT > 25.0)
18991914
histos.fill(HIST("missed_kpi_INJets_25_infinite"), (bestJetpT - missing_pt) / bestJetpT, lResonance.Pt());
19001915

1901-
if (bestJetpT > 8.0)
1916+
if (bestJetpT > 8.0) {
19021917
histos.fill(HIST("missed_kpi_INJets_8_infinite"), (bestJetpT - missing_pt) / bestJetpT, lResonance.Pt());
1903-
1918+
histos.fill(HIST("normalJetpT_8_kstarSpectra"), lResonance.Pt());
1919+
}
19041920
histos.fill(HIST("JetMigration"), bestJetpT, recoveredJetpT);
19051921
} // cfgJetdRHistos
19061922
} // kaon_out || pion_out

0 commit comments

Comments
 (0)