Skip to content

Commit 37cfc2b

Browse files
YubiaoWangYubiaoWang
authored
Add files via upload
1 parent 0a33ba2 commit 37cfc2b

1 file changed

Lines changed: 20 additions & 39 deletions

File tree

‎PWGJE/Tasks/jetChargedV2.cxx‎

Lines changed: 20 additions & 39 deletions
Original file line numberDiff line numberDiff line change
@@ -276,6 +276,8 @@ struct JetChargedV2 {
276276

277277
//< RC test plots >//
278278
registry.add("h3_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}, {40, 0., o2::constants::math::TwoPI}}});
279+
registry.add("h3_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}, {40, 0., o2::constants::math::TwoPI}}});
280+
279281
registry.add("h3_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}, {40, 0., o2::constants::math::TwoPI}}});
280282
registry.add("h3_centrality_deltapT_RandomCornPhi_localrhovsphiwithoutoneleadingjet", "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}, {40, 0., o2::constants::math::TwoPI}}});
281283
registry.add("h3_centrality_deltapT_RandomCornPhi_localrhovsphiwithouttwoleadingjet", "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}, {40, 0., o2::constants::math::TwoPI}}});
@@ -617,6 +619,10 @@ struct JetChargedV2 {
617619
registry.fill(HIST("h_v3obs_centrality"), collision.centrality(), temppara[3]);
618620
registry.fill(HIST("h_evtnum_centrlity"), evtnum, collision.centrality());
619621

622+
if (temppara[0] == 0) {
623+
return;
624+
}
625+
620626
int nDF = 1;
621627
int numOfFreePara = 2;
622628
nDF = static_cast<int>(fFitModulationV2v3->GetXaxis()->GetNbins()) - numOfFreePara;
@@ -646,27 +652,6 @@ struct JetChargedV2 {
646652
} else if (evtcent >= 30 && evtcent <= 50) {
647653
registry.fill(HIST("h2_PChi2_CombinFitB"), cDF, chiSqr / (static_cast<float>(nDF)));
648654
}
649-
650-
// integral local rho
651-
double integralValue = 0.;
652-
for (auto const& jet : jets) {
653-
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
654-
continue;
655-
}
656-
if (!isAcceptedJet<aod::JetTracks>(jet)) {
657-
continue;
658-
}
659-
if (jet.r() != round(selectedJetsRadius * 100.0f)) {
660-
continue;
661-
}
662-
integralValue = fFitModulationV2v3->Integral(jet.phi() - jetRadius, jet.phi() + jetRadius);
663-
}
664-
double rholocal = collision.rho() / (2 * jetRadius * temppara[0]) * integralValue;
665-
registry.fill(HIST("h2_phi_rholocal_cent"), collision.centrality(), rholocal, 1.0);
666-
registry.fill(HIST("h2_phi_rhomedian_cent"), collision.centrality(), collision.rho(), 1.0);
667-
668-
// float evtPl2 = 0.;
669-
// float evtPl3 = 0.;
670655
for (uint i = 0; i < cfgnMods->size(); i++) {
671656
int nmode = cfgnMods->at(i);
672657
int detInd = detId * 4 + cfgnTotalSystem * 4 * (nmode - 2);
@@ -681,14 +666,19 @@ struct JetChargedV2 {
681666
if (jet.r() != round(selectedJetsRadius * 100.0f)) {
682667
continue;
683668
}
669+
670+
double integralValue = fFitModulationV2v3->Integral(jet.phi() - jetRadius, jet.phi() + jetRadius);
671+
double rholocal = collision.rho() / (2 * jetRadius * temppara[0]) * integralValue;
672+
registry.fill(HIST("h2_phi_rholocal_cent"), collision.centrality(), rholocal, 1.0);
673+
registry.fill(HIST("h2_phi_rhomedian_cent"), collision.centrality(), collision.rho(), 1.0);
674+
684675
if (nmode == 2) {
685676
registry.fill(HIST("h_jet_pt_rholocal"), jet.pt() - (rholocal * jet.area()), 1.0);
686677

687678
double phiMinusPsi2;
688679
if (collision.qvecAmp()[detId] < 1e-8) {
689680
continue;
690681
}
691-
// evtPl2 = helperEP.GetEventPlane(collision.qvecRe()[detInd], collision.qvecIm()[detInd], nmode);
692682
phiMinusPsi2 = jet.phi() - ep2;
693683

694684
registry.fill(HIST("h2_phi_rholocal"), jet.phi() - ep2, rholocal, 1.0);
@@ -727,6 +717,9 @@ struct JetChargedV2 {
727717
float randomConePhi = randomNumber.Uniform(0.0, o2::constants::math::TwoPI);
728718
float randomConePt = 0;
729719

720+
double integralValueRC = fFitModulationV2v3->Integral(randomConePhi - randomConeR, randomConePhi + randomConeR);
721+
double rholocalRC = collision.rho() / (2 * randomConeR * temppara[0]) * integralValueRC;
722+
730723
// int nmode = cfgnMods->at(i);
731724
// int detInd = detId * 4 + cfgnTotalSystem * 4 * (nmode - 2);
732725

@@ -750,17 +743,8 @@ struct JetChargedV2 {
750743
}
751744
}
752745
}
753-
// double integralValue = 0;
754-
// double rholocal = 0;
755-
// for (auto const& jet : jets) {
756-
// if (temppara[0] == 0) {
757-
// break;
758-
// }
759-
// integralValue = fFitModulationV2v3->Integral(jet.phi() - jetRadius, jet.phi() + jetRadius);
760-
// rholocal = collision.rho() / (2 * jetRadius * temppara[0]) * integralValue;
761-
// registry.fill(HIST("h2_phi_rholocal_RC"), jet.phi() - evtPl2, rholocal, 1.0);
762-
// registry.fill(HIST("h2_phi_rhomedian_RC"), jet.phi() - evtPl2, collision.rho(), 1.0);
763-
// }
746+
registry.fill(HIST("h3_centrality_deltapT_RandomCornPhi_localrhovsphi"), collision.centrality(), randomConePt - o2::constants::math::PI * randomConeR * randomConeR * rholocalRC, rcPhiPsi2, 1.0);
747+
764748
// removing the leading jet from the random cone
765749
if (jets.size() > 0) { // if there are no jets in the acceptance (from the jetfinder cuts) then there can be no leading jet
766750
float dPhiLeadingJet = RecoDecay::constrainAngle(jets.iteratorAt(0).phi() - randomConePhi, static_cast<float>(-o2::constants::math::PI));
@@ -787,10 +771,7 @@ struct JetChargedV2 {
787771
}
788772
}
789773
}
790-
if (temppara[0] == 0) {
791-
break;
792-
}
793-
registry.fill(HIST("h3_centrality_deltapT_RandomCornPhi_localrhovsphiwithoutleadingjet"), collision.centrality(), randomConePt - o2::constants::math::PI * randomConeR * randomConeR * rholocal, rcPhiPsi2, 1.0);
774+
registry.fill(HIST("h3_centrality_deltapT_RandomCornPhi_localrhovsphiwithoutleadingjet"), collision.centrality(), randomConePt - o2::constants::math::PI * randomConeR * randomConeR * rholocalRC, rcPhiPsi2, 1.0);
794775
registry.fill(HIST("h3_centrality_deltapT_RandomCornPhi_rhorandomconewithoutleadingjet"), collision.centrality(), randomConePt - o2::constants::math::PI * randomConeR * randomConeR * collision.rho(), rcPhiPsi2, 1.0);
795776

796777
// randomised eta,phi for tracks, to assess part of fluctuations coming from statistically independently emitted particles, removing tracks from 2 leading jets
@@ -810,8 +791,8 @@ struct JetChargedV2 {
810791
}
811792
}
812793
}
813-
registry.fill(HIST("h3_centrality_deltapT_RandomCornPhi_localrhovsphiwithoutoneleadingjet"), collision.centrality(), randomConePtWithoutOneLeadJet - o2::constants::math::PI * randomConeR * randomConeR * rholocal, rcPhiPsi2, 1.0);
814-
registry.fill(HIST("h3_centrality_deltapT_RandomCornPhi_localrhovsphiwithouttwoleadingjet"), collision.centrality(), randomConePtWithoutTwoLeadJet - o2::constants::math::PI * randomConeR * randomConeR * rholocal, rcPhiPsi2, 1.0);
794+
registry.fill(HIST("h3_centrality_deltapT_RandomCornPhi_localrhovsphiwithoutoneleadingjet"), collision.centrality(), randomConePtWithoutOneLeadJet - o2::constants::math::PI * randomConeR * randomConeR * rholocalRC, rcPhiPsi2, 1.0);
795+
registry.fill(HIST("h3_centrality_deltapT_RandomCornPhi_localrhovsphiwithouttwoleadingjet"), collision.centrality(), randomConePtWithoutTwoLeadJet - o2::constants::math::PI * randomConeR * randomConeR * rholocalRC, rcPhiPsi2, 1.0);
815796
// } else if (nmode == 3) {
816797
// continue;
817798
// }

0 commit comments

Comments
 (0)