Skip to content

Commit b6ffa56

Browse files
YubiaoWangYubiaoWang
authored
Add files via upload
1 parent af4f88b commit b6ffa56

1 file changed

Lines changed: 68 additions & 2 deletions

File tree

‎PWGJE/Tasks/jetChargedV2.cxx‎

Lines changed: 68 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -307,12 +307,19 @@ struct JetChargedV2 {
307307
histosQA.add(Form("histEvtPlTwistV%d", cfgnMods->at(i)), "", {HistType::kTH2F, {axisEvtPl, axisCent}});
308308
histosQA.add(Form("histEvtPlFinalV%d", cfgnMods->at(i)), "", {HistType::kTH2F, {axisEvtPl, axisCent}});
309309

310+
histosQA.add(Form("h_ep2_FT0CV%d", cfgnMods->at(i)), "", {HistType::kTH2F, {axisEvtPl, axisCent}});
311+
histosQA.add(Form("h_ep2_FT0MV%d", cfgnMods->at(i)), "", {HistType::kTH2F, {axisEvtPl, axisCent}});
312+
310313
histosQA.add(Form("histEvtPlRes_SigRefAV%d", cfgnMods->at(i)), "", {HistType::kTH2F, {axisEvtPl, axisCent}});
311314
histosQA.add(Form("histEvtPlRes_SigRefBV%d", cfgnMods->at(i)), "", {HistType::kTH2F, {axisEvtPl, axisCent}});
312315
histosQA.add(Form("histEvtPlRes_RefARefBV%d", cfgnMods->at(i)), "", {HistType::kTH2F, {axisEvtPl, axisCent}});
313316
}
314317
histosQA.add("histCent", "Centrality TrkProcess", HistType::kTH1F, {axisCent});
315318

319+
// histosQA.add("h2_ep2_FT0C_FT0M", "event_plane_FT0C_FT0M; FT0C; FT0M", {HistType::kTH3F, {{100, 0.0, 100.0}, {100, 0.0, 100.0}, {100, 0.0, 100.0}}});
320+
histosQA.add("h2_ep2_FT0C_FT0M", "event_plane_FT0C_FT0M; FT0C; FT0M; cent", {HistType::kTH3F, {{100, -o2::constants::math::PIHalf., o2::constants::math::PIHalf.}, {100, -o2::constants::math::PIHalf., o2::constants::math::PIHalf}, {100, 0.0, 100.0}}});
321+
histosQA.add("h2_ep2_FT0C_FT0M_bumpRegion", "event_plane_FT0C_FT0M_bumpRegion; FT0C; FT0M; cent", {HistType::kTH3F, {{100, -o2::constants::math::PIHalf, o2::constants::math::PIHalf}, {100, -o2::constants::math::PIHalf, o2::constants::math::PIHalf}, {100, 0.0, 100.0}}});
322+
316323
//< fit quality >//
317324
registry.add("h_PvalueCDF_CombinFit", "cDF #chi^{2}; entries", {HistType::kTH1F, {{50, 0, 1}}});
318325
registry.add("h2_PvalueCDFCent_CombinFit", "p-value cDF vs centrality; centrality; p-value", {HistType::kTH2F, {{100, 0, 100}, {40, 0, 1}}});
@@ -334,7 +341,11 @@ struct JetChargedV2 {
334341

335342
registry.add("h2_phi_rholocal", "#varphi vs #rho(#varphi); #varphi - #Psi_{EP,2}; #rho(#varphi) ", {HistType::kTH2F, {{40, 0., o2::constants::math::TwoPI}, {210, -10.0, 200.0}}});
336343
registry.add("h2_phi_rholocal_absDelta", "#varphi vs #rho(#varphi), absDelta; #varphi - #Psi_{EP,2}; #rho(#varphi) ", {HistType::kTH2F, {{40, 0., o2::constants::math::TwoPI}, {210, -10.0, 200.0}}});
344+
registry.add("h2_phi_mediumrho_absDelta", "#varphi vs #rho(#varphi), absDelta; #varphi - #Psi_{EP,2}; #rho(#varphi) ", {HistType::kTH2F, {{40, 0., o2::constants::math::TwoPI}, {210, -10.0, 200.0}}});
337345

346+
registry.add("h2_phi_rholocal_absDelta_low", "#varphi vs #rho(#varphi), absDelta; #varphi - #Psi_{EP,2}; #rho(#varphi) ", {HistType::kTH2F, {{40, 0., o2::constants::math::TwoPI}, {210, -10.0, 200.0}}});
347+
registry.add("h2_phi_rholocal_absDelta_medium", "#varphi vs #rho(#varphi), absDelta; #varphi - #Psi_{EP,2}; #rho(#varphi) ", {HistType::kTH2F, {{40, 0., o2::constants::math::TwoPI}, {210, -10.0, 200.0}}});
348+
registry.add("h2_phi_rholocal_absDelta_high", "#varphi vs #rho(#varphi), absDelta; #varphi - #Psi_{EP,2}; #rho(#varphi) ", {HistType::kTH2F, {{40, 0., o2::constants::math::TwoPI}, {210, -10.0, 200.0}}});
338349
registry.add("h2_rholocal_cent", "#varphi vs #rho(#varphi); #cent; #rho(#varphi) ", {HistType::kTH2F, {{100, 0., 100}, {210, -10.0, 200.0}}});
339350
//< \sigma p_T at local rho test plot | end >
340351

@@ -1239,6 +1250,7 @@ struct JetChargedV2 {
12391250
for (uint i = 0; i < cfgnMods->size(); i++) {
12401251
int nmode = cfgnMods->at(i);
12411252
int detInd = detId * 4 + cfgnTotalSystem * 4 * (nmode - 2);
1253+
int detIndFT0M = (detId + 2) * 4 + cfgnTotalSystem * 4 * (nmode - 2);
12421254
int refAInd = refAId * 4 + cfgnTotalSystem * 4 * (nmode - 2);
12431255
int refBInd = refBId * 4 + cfgnTotalSystem * 4 * (nmode - 2);
12441256

@@ -1254,6 +1266,11 @@ struct JetChargedV2 {
12541266
histosQA.fill(HIST("histEvtPlTwistV2"), helperEP.GetEventPlane(collision.qvecRe()[detInd + 2], collision.qvecIm()[detInd + 2], nmode), collision.cent());
12551267
histosQA.fill(HIST("histEvtPlFinalV2"), helperEP.GetEventPlane(collision.qvecRe()[detInd + 3], collision.qvecIm()[detInd + 3], nmode), collision.cent());
12561268

1269+
// mark
1270+
histosQA.fill(HIST("h_ep2_FT0CV2"), helperEP.GetEventPlane(collision.qvecRe()[detInd + 3], collision.qvecIm()[detInd + 3], nmode), collision.cent());
1271+
histosQA.fill(HIST("h_ep2_FT0MV2"), helperEP.GetEventPlane(collision.qvecRe()[detIndFT0M + 3], collision.qvecIm()[detIndFT0M + 3], nmode), collision.cent());
1272+
histosQA.fill(HIST("h2_ep2_FT0C_FT0M"), helperEP.GetEventPlane(collision.qvecRe()[detInd + 3], collision.qvecIm()[detInd + 3], nmode), helperEP.GetEventPlane(collision.qvecRe()[detIndFT0M + 3], collision.qvecIm()[detIndFT0M + 3], nmode), collision.cent());
1273+
12571274
histosQA.fill(HIST("histEvtPlRes_SigRefAV2"), helperEP.GetResolution(helperEP.GetEventPlane(collision.qvecRe()[detInd + 3], collision.qvecIm()[detInd + 3], nmode), helperEP.GetEventPlane(collision.qvecRe()[refAInd + 3], collision.qvecIm()[refAInd + 3], nmode), nmode), collision.cent());
12581275
histosQA.fill(HIST("histEvtPlRes_SigRefBV2"), helperEP.GetResolution(helperEP.GetEventPlane(collision.qvecRe()[detInd + 3], collision.qvecIm()[detInd + 3], nmode), helperEP.GetEventPlane(collision.qvecRe()[refBInd + 3], collision.qvecIm()[refBInd + 3], nmode), nmode), collision.cent());
12591276
histosQA.fill(HIST("histEvtPlRes_RefARefBV2"), helperEP.GetResolution(helperEP.GetEventPlane(collision.qvecRe()[refAInd + 3], collision.qvecIm()[refAInd + 3], nmode), helperEP.GetEventPlane(collision.qvecRe()[refBInd + 3], collision.qvecIm()[refBInd + 3], nmode), nmode), collision.cent());
@@ -1539,6 +1556,7 @@ struct JetChargedV2 {
15391556
for (uint i = 0; i < cfgnMods->size(); i++) {
15401557
int nmode = cfgnMods->at(i);
15411558
int detInd = detId * 4 + cfgnTotalSystem * 4 * (nmode - 2);
1559+
15421560
if (nmode == cfgNmodA) {
15431561
if (collision.qvecAmp()[detId] > collQvecAmpDetId) {
15441562
ep2 = helperEP.GetEventPlane(collision.qvecRe()[detInd + 3], collision.qvecIm()[detInd + 3], nmode);
@@ -1551,7 +1569,7 @@ struct JetChargedV2 {
15511569
}
15521570

15531571
const char* fitFunctionV2v3 = "[0] * (1. + 2. * ([1] * std::cos(2. * (x - [2])) + [3] * std::cos(3. * (x - [4]))))";
1554-
fFitModulationV2v3 = new TF1("fit_kV3", fitFunctionV2v3, 0, o2::constants::math::TwoPI);
1572+
fFitModulationV2v3 = new TF1("fit_kV3", fitFunctionV2v3, -2, o2::constants::math::TwoPI + 2);
15551573
//=========================< set parameter >=========================//
15561574
fFitModulationV2v3->SetParameter(0, 1.);
15571575
fFitModulationV2v3->SetParameter(1, 0.01);
@@ -1638,7 +1656,11 @@ struct JetChargedV2 {
16381656
for (uint i = 0; i < cfgnMods->size(); i++) {
16391657
int nmode = cfgnMods->at(i);
16401658
int detInd = detId * 4 + cfgnTotalSystem * 4 * (nmode - 2);
1659+
int detIndFT0M = 2 * 4 + cfgnTotalSystem * 4 * (nmode - 2);
1660+
int bumpLow = 75;
1661+
int bumpUp = 125;
16411662

1663+
bool hasBumpJet = false;
16421664
for (auto const& jet : jets) {
16431665
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
16441666
continue;
@@ -1649,8 +1671,32 @@ struct JetChargedV2 {
16491671
if (jet.r() != round(selectedJetsRadius * 100.0f)) {
16501672
continue;
16511673
}
1674+
if (jet.pt() > bumpLow && jet.pt() < bumpUp) {
1675+
hasBumpJet = true;
1676+
}
1677+
}
1678+
if (hasBumpJet && nmode == cfgNmodA) {
1679+
if (collision.qvecAmp()[detId] >= collQvecAmpDetId) {
1680+
histosQA.fill(HIST("h2_ep2_FT0C_FT0M_bumpRegion"), helperEP.GetEventPlane(collision.qvecRe()[detInd + 3], collision.qvecIm()[detInd + 3], nmode), helperEP.GetEventPlane(collision.qvecRe()[detIndFT0M + 3], collision.qvecIm()[detIndFT0M + 3], nmode), collision.cent());
1681+
}
1682+
}
16521683

1653-
double integralValue = fFitModulationV2v3->Integral(jet.phi() - selectedJetsRadius, jet.phi() + selectedJetsRadius);
1684+
for (auto const& jet : jets) {
1685+
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
1686+
continue;
1687+
}
1688+
if (!isAcceptedJet<aod::JetTracks>(jet)) {
1689+
continue;
1690+
}
1691+
if (jet.r() != round(selectedJetsRadius * 100.0f)) {
1692+
continue;
1693+
}
1694+
double twoPi = o2::constants::math::TwoPI;
1695+
double phi = std::fmod(jet.phi(), twoPi);
1696+
if (phi < 0) {
1697+
phi += twoPi;
1698+
}
1699+
double integralValue = fFitModulationV2v3->Integral(phi - selectedJetsRadius, phi + selectedJetsRadius);
16541700
double rholocal = collision.rho() / (2 * selectedJetsRadius * temppara[0]) * integralValue;
16551701
registry.fill(HIST("h2_rholocal_cent"), centrality, rholocal, 1.0);
16561702

@@ -1664,7 +1710,27 @@ struct JetChargedV2 {
16641710
float absDelta = std::abs(phiMinusPsi2);
16651711
registry.fill(HIST("h2_phi_rholocal"), jet.phi() - ep2, rholocal, 1.0);
16661712
registry.fill(HIST("h2_phi_rholocal_absDelta"), absDelta, rholocal, 1.0);
1713+
registry.fill(HIST("h2_phi_mediumrho_absDelta"), absDelta, collision.rho(), 1.0);
1714+
1715+
int lowPtCut = 20;
1716+
int mediumPtCut = 40;
1717+
int highPtCut = 70;
1718+
int highPtCutEnd = 100;
1719+
1720+
if (jet.pt() >= lowPtCut && jet.pt() < mediumPtCut) {
1721+
registry.fill(HIST("h2_phi_rholocal_absDelta_low"), absDelta, rholocal, 1.0);
1722+
} else if (jet.pt() >= mediumPtCut && jet.pt() < highPtCut) {
1723+
registry.fill(HIST("h2_phi_rholocal_absDelta_medium"), absDelta, rholocal, 1.0);
1724+
} else if (jet.pt() >= highPtCut && jet.pt() < highPtCutEnd) {
1725+
registry.fill(HIST("h2_phi_rholocal_absDelta_high"), absDelta, rholocal, 1.0);
1726+
}
1727+
16671728
registry.fill(HIST("h_jet_pt_inclusive_v2_rho"), jet.pt() - (rholocal * jet.area()), 1.0);
1729+
// mark
1730+
if (jet.pt() > bumpLow && jet.pt() < bumpUp) {
1731+
histosQA.fill(HIST("h2_ep2_FT0C_FT0M_bumpRegion"), helperEP.GetEventPlane(collision.qvecRe()[detInd + 3], collision.qvecIm()[detInd + 3], nmode), helperEP.GetEventPlane(collision.qvecRe()[detIndFT0M + 3], collision.qvecIm()[detIndFT0M + 3], nmode), collision.cent());
1732+
}
1733+
16681734
if ((absDelta < o2::constants::math::PIQuarter) || (absDelta >= evtPlnAngleA * o2::constants::math::PIQuarter) || (absDelta >= evtPlnAngleB * o2::constants::math::PIQuarter && absDelta < evtPlnAngleC * o2::constants::math::PIQuarter)) {
16691735
registry.fill(HIST("h_jet_pt_in_plane_v2_rho"), jet.pt() - (rholocal * jet.area()), 1.0);
16701736
registry.fill(HIST("h2_centrality_jet_pt_in_plane_v2_rho"), centrality, jet.pt() - (rholocal * jet.area()), 1.0);

0 commit comments

Comments
 (0)