Skip to content

Commit b9e64fd

Browse files
YubiaoWangYubiaoWang
authored
jetChargedV2.cxx
1 parent ec3d34b commit b9e64fd

File tree

1 file changed

+21
-20
lines changed

1 file changed

+21
-20
lines changed

‎PWGJE/Tasks/jetChargedV2.cxx‎

Lines changed: 21 additions & 20 deletions
Original file line numberDiff line numberDiff line change
@@ -273,16 +273,16 @@ struct JetChargedV2 {
273273
registry.add("h_mcp_jet_pt_out_of_plane_v2_rho", "jet pT;#it{p}^{out-of-plane}_{T,jet} (GeV/#it{c});entries", {HistType::kTH1F, {jetPtAxisRhoAreaSub}});
274274
registry.add("h_mcp_jet_pt_in_plane_v3_rho", "jet pT;#it{p}^{in-plane}_{T,jet} (GeV/#it{c});entries", {HistType::kTH1F, {jetPtAxisRhoAreaSub}});
275275
registry.add("h_mcp_jet_pt_out_of_plane_v3_rho", "jet pT;#it{p}^{out-of-plane}_{T,jet} (GeV/#it{c});entries", {HistType::kTH1F, {jetPtAxisRhoAreaSub}});
276-
276+
277277
registry.add("h2_mcp_centrality_jet_pt_in_plane_v2_rho", "centrality vs #it{p}^{in-plane}_{T,jet}; centrality; #it{p}_{T,jet} (GeV/#it{c})", {HistType::kTH2F, {{120, -10.0, 110.0}, jetPtAxisRhoAreaSub}});
278278
registry.add("h2_mcp_centrality_jet_pt_out_of_plane_v2_rho", "centrality vs #it{p}^{out-of-plane}_{T,jet}; centrality; #it{p}_{T,jet} (GeV/#it{c})", {HistType::kTH2F, {{120, -10.0, 110.0}, jetPtAxisRhoAreaSub}});
279279
registry.add("h2_mcp_centrality_jet_pt_in_plane_v3_rho", "centrality vs #it{p}^{in-plane}_{T,jet}; centrality; #it{p}_{T,jet} (GeV/#it{c})", {HistType::kTH2F, {{120, -10.0, 110.0}, jetPtAxisRhoAreaSub}});
280280
registry.add("h2_mcp_centrality_jet_pt_out_of_plane_v3_rho", "centrality vs #it{p}^{out-of-plane}_{T,jet}; centrality; #it{p}_{T,jet} (GeV/#it{c})", {HistType::kTH2F, {{120, -10.0, 110.0}, jetPtAxisRhoAreaSub}});
281-
281+
282282
registry.add("h3_mcp_centrality_deltapT_RandomCornPhi_localrhovsphi", "centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho}; #Delta#varphi_{jet}", {HistType::kTH3F, {{100, 0.0, 100.0}, {400, -200.0, 200.0}, {100, 0., o2::constants::math::TwoPI}}});
283283
registry.add("h3_mcp_centrality_deltapT_RandomCornPhi_rhorandomconewithoutleadingjet", "centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho}; #Delta#varphi_{jet}", {HistType::kTH3F, {{100, 0.0, 100.0}, {400, -200.0, 200.0}, {100, 0., o2::constants::math::TwoPI}}});
284284
registry.add("h3_mcp_centrality_deltapT_RandomCornPhi_localrhovsphiwithoutleadingjet", "centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho}(#varphi); #Delta#varphi_{jet}", {HistType::kTH3F, {{100, 0.0, 100.0}, {400, -200.0, 200.0}, {100, 0., o2::constants::math::TwoPI}}});
285-
285+
286286
registry.add("h_mcColl_counts_areasub", " number of mc events; event status; entries", {HistType::kTH1F, {{10, 0, 10}}});
287287
registry.get<TH1>(HIST("h_mcColl_counts_areasub"))->GetXaxis()->SetBinLabel(1, "allMcColl");
288288
registry.get<TH1>(HIST("h_mcColl_counts_areasub"))->GetXaxis()->SetBinLabel(2, "vertexZ");
@@ -292,7 +292,7 @@ struct JetChargedV2 {
292292
registry.get<TH1>(HIST("h_mcColl_counts_areasub"))->GetXaxis()->SetBinLabel(6, "centralitycut");
293293
registry.get<TH1>(HIST("h_mcColl_counts_areasub"))->GetXaxis()->SetBinLabel(7, "occupancycut");
294294
registry.add("h_mcColl_rho", "mc collision rho;#rho (GeV/#it{c}); counts", {HistType::kTH1F, {{500, 0.0, 500.0}}});
295-
295+
296296
//< \sigma p_T at local rho test plot >
297297
registry.add("h_accept_Track", "all and accept track;Track;entries", {HistType::kTH1F, {{10, 0.0, 10.0}}});
298298
registry.get<TH1>(HIST("h_accept_Track"))->GetXaxis()->SetBinLabel(1, "acceptTrk");
@@ -498,7 +498,8 @@ struct JetChargedV2 {
498498

499499
template <typename T, typename J>
500500
void fillInOutjetV2EP(T const& collision, J const& jets)
501-
{
501+
{
502+
float collQvecAmpDetId = 1e-8;
502503
int cfgNmodA = 2;
503504
int cfgNmodB = 3;
504505
for (uint i = 0; i < cfgnMods->size(); i++) {
@@ -508,7 +509,7 @@ struct JetChargedV2 {
508509
int refBInd = refBId * 4 + cfgnTotalSystem * 4 * (nmode - 2);
509510

510511
if (nmode == cfgNmodA) {
511-
if (collision.qvecAmp()[detId] > 1e-8 || collision.qvecAmp()[refAId] < 1e-8 || collision.qvecAmp()[refBId] < 1e-8) {
512+
if (collision.qvecAmp()[detId] > collQvecAmpDetId || collision.qvecAmp()[refAId] < collQvecAmpDetId || collision.qvecAmp()[refBId] < collQvecAmpDetId) {
512513
histosQA.fill(HIST("histQvecUncorV2"), collision.qvecRe()[detInd], collision.qvecIm()[detInd], collision.cent());
513514
histosQA.fill(HIST("histQvecRectrV2"), collision.qvecRe()[detInd + 1], collision.qvecIm()[detInd + 1], collision.cent());
514515
histosQA.fill(HIST("histQvecTwistV2"), collision.qvecRe()[detInd + 2], collision.qvecIm()[detInd + 2], collision.cent());
@@ -549,7 +550,7 @@ struct JetChargedV2 {
549550

550551
if (nmode == cfgNmodA) {
551552
double phiMinusPsi2;
552-
if (collision.qvecAmp()[detId] < 1e-8) {
553+
if (collision.qvecAmp()[detId] < collQvecAmpDetId) {
553554
continue;
554555
}
555556
float ep2 = helperEP.GetEventPlane(collision.qvecRe()[detInd], collision.qvecIm()[detInd], nmode);
@@ -564,7 +565,7 @@ struct JetChargedV2 {
564565
continue;
565566
}
566567
registry.fill(HIST("h_jet_pt_rhoareasubtracted"), jet.pt() - (collision.rho() * jet.area()), 1.0);
567-
568+
568569
phiMinusPsi2 = jet.phi() - ep2;
569570
if ((phiMinusPsi2 < o2::constants::math::PIQuarter) || (phiMinusPsi2 >= evtPlnAngleA * o2::constants::math::PIQuarter) || (phiMinusPsi2 >= evtPlnAngleB * o2::constants::math::PIQuarter && phiMinusPsi2 < evtPlnAngleC * o2::constants::math::PIQuarter)) {
570571
registry.fill(HIST("h_jet_pt_in_plane_v2"), jet.pt() - (collision.rho() * jet.area()), 1.0);
@@ -588,7 +589,7 @@ struct JetChargedV2 {
588589
continue;
589590
}
590591
phiMinusPsi3 = jet.phi() - ep3;
591-
592+
592593
if ((phiMinusPsi3 < o2::constants::math::PIQuarter) || (phiMinusPsi3 >= evtPlnAngleA * o2::constants::math::PIQuarter) || (phiMinusPsi3 >= evtPlnAngleB * o2::constants::math::PIQuarter && phiMinusPsi3 < evtPlnAngleC * o2::constants::math::PIQuarter)) {
593594
registry.fill(HIST("h_jet_pt_in_plane_v3"), jet.pt() - (collision.rho() * jet.area()), 1.0);
594595
registry.fill(HIST("h2_centrality_jet_pt_in_plane_v3"), collision.centrality(), jet.pt() - (collision.rho() * jet.area()), 1.0);
@@ -648,7 +649,7 @@ struct JetChargedV2 {
648649
registry.fill(HIST("h_accept_Track"), 1.5);
649650
}
650651
}
651-
652+
652653
for (auto const& track : tracks) {
653654
registry.fill(HIST("h_accept_Track"), 2.5);
654655
if (jetderiveddatautilities::selectTrack(track, trackSelection) && (std::fabs(track.eta() - leadingJetEta) > jetRadius) && track.pt() >= localRhoFitPtMin && track.pt() <= localRhoFitPtMax) {
@@ -678,11 +679,11 @@ struct JetChargedV2 {
678679
int nmode = cfgnMods->at(i);
679680
int detInd = detId * 4 + cfgnTotalSystem * 4 * (nmode - 2);
680681
if (nmode == cfgNmodA) {
681-
if (collision.qvecAmp()[detId] > 1e-8) {
682+
if (collision.qvecAmp()[detId] > collQvecAmpDetId) {
682683
ep2 = helperEP.GetEventPlane(collision.qvecRe()[detInd + 3], collision.qvecIm()[detInd + 3], nmode);
683684
}
684685
} else if (nmode == cfgNmodB) {
685-
if (collision.qvecAmp()[detId] > 1e-8) {
686+
if (collision.qvecAmp()[detId] > collQvecAmpDetId) {
686687
ep3 = helperEP.GetEventPlane(collision.qvecRe()[detInd + 3], collision.qvecIm()[detInd + 3], nmode);
687688
}
688689
}
@@ -769,7 +770,7 @@ struct JetChargedV2 {
769770
int evtMidAreaMax = 50;
770771
if (evtcent >= evtCentAreaMin && evtcent <= evtCentAreaMax) {
771772
registry.fill(HIST("h2_PChi2_CombinFitA"), cDF, chiSqr / (static_cast<float>(nDF)));
772-
773+
773774
} else if (evtcent >= evtMidAreaMin && evtcent <= evtMidAreaMax) {
774775
registry.fill(HIST("h2_PChi2_CombinFitB"), cDF, chiSqr / (static_cast<float>(nDF)));
775776
}
@@ -800,7 +801,7 @@ struct JetChargedV2 {
800801
registry.fill(HIST("h_jet_pt_rholocal"), jet.pt() - (rholocal * jet.area()), 1.0);
801802

802803
double phiMinusPsi2;
803-
if (collision.qvecAmp()[detId] < 1e-8) {
804+
if (collision.qvecAmp()[detId] < collQvecAmpDetId) {
804805
continue;
805806
}
806807
phiMinusPsi2 = jet.phi() - ep2;
@@ -817,7 +818,7 @@ struct JetChargedV2 {
817818
}
818819
} else if (nmode == cfgNmodB) {
819820
double phiMinusPsi3;
820-
if (collision.qvecAmp()[detId] < 1e-8) {
821+
if (collision.qvecAmp()[detId] < collQvecAmpDetId) {
821822
continue;
822823
}
823824
ep3 = helperEP.GetEventPlane(collision.qvecRe()[detInd], collision.qvecIm()[detInd], nmode);
@@ -937,11 +938,11 @@ struct JetChargedV2 {
937938
int nmode = cfgnMods->at(i);
938939
int detInd = detId * 4 + cfgnTotalSystem * 4 * (nmode - 2);
939940
if (nmode == cfgNmodA) {
940-
if (collision.qvecAmp()[detId] > 1e-8) {
941+
if (collision.qvecAmp()[detId] > collQvecAmpDetId) {
941942
ep2 = helperEP.GetEventPlane(collision.qvecRe()[detInd + 3], collision.qvecIm()[detInd + 3], nmode);
942943
}
943944
} else if (nmode == cfgNmodB) {
944-
if (collision.qvecAmp()[detId] > 1e-8) {
945+
if (collision.qvecAmp()[detId] > collQvecAmpDetId) {
945946
ep3 = helperEP.GetEventPlane(collision.qvecRe()[detInd + 3], collision.qvecIm()[detInd + 3], nmode);
946947
}
947948
}
@@ -1014,7 +1015,7 @@ struct JetChargedV2 {
10141015
registry.fill(HIST("h_mcp_jet_pt_rholocal"), jet.pt() - (rholocal * jet.area()), 1.0);
10151016

10161017
double phiMinusPsi2;
1017-
if (collision.qvecAmp()[detId] < 1e-8) {
1018+
if (collision.qvecAmp()[detId] < collQvecAmpDetId) {
10181019
continue;
10191020
}
10201021
phiMinusPsi2 = jet.phi() - ep2;
@@ -1030,7 +1031,7 @@ struct JetChargedV2 {
10301031
}
10311032
} else if (nmode == cfgNmodB) {
10321033
double phiMinusPsi3;
1033-
if (collision.qvecAmp()[detId] < 1e-8) {
1034+
if (collision.qvecAmp()[detId] < collQvecAmpDetId) {
10341035
continue;
10351036
}
10361037
ep3 = helperEP.GetEventPlane(collision.qvecRe()[detInd], collision.qvecIm()[detInd], nmode);
@@ -1423,4 +1424,4 @@ struct JetChargedV2 {
14231424
WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)
14241425
{
14251426
return WorkflowSpec{adaptAnalysisTask<JetChargedV2>(cfgc)};
1426-
}
1427+
}

0 commit comments

Comments
 (0)