Skip to content

Commit a7a4a60

Browse files
wefeng1110wefeng1110
authored
formatting change for the code
1 parent ea38636 commit a7a4a60

File tree

1 file changed

+55
-57
lines changed

1 file changed

+55
-57
lines changed

‎PWGJE/Tasks/jetSpectraCharged.cxx‎

Lines changed: 55 additions & 57 deletions
Original file line numberDiff line numberDiff line change
@@ -83,7 +83,6 @@ struct JetSpectraCharged {
8383
Configurable<bool> checkGeoPtMatched{"checkGeoPtMatched", false, "0: turn off geometry and pT matching, 1: do geometry and pT matching"};
8484
Configurable<int> acceptSplitCollisions{"acceptSplitCollisions", 0, "0: only look at mcCollisions that are not split; 1: accept split mcCollisions, 2: accept split mcCollisions but only look at the first reco collision associated with it"};
8585

86-
8786
std::vector<int> eventSelectionBits;
8887
int trackSelection = -1;
8988

@@ -111,7 +110,7 @@ struct JetSpectraCharged {
111110
}
112111
}
113112

114-
if (doprocessSpectraData || doprocessSpectraMCD || doprocessSpectraMCDWeighted){
113+
if (doprocessSpectraData || doprocessSpectraMCD || doprocessSpectraMCDWeighted) {
115114
registry.add("h_jet_pt", "jet pT;#it{p}_{T,jet} (GeV/#it{c}); counts", {HistType::kTH1F, {jetPtAxis}});
116115
registry.add("h_jet_eta", "jet eta;#eta; counts", {HistType::kTH1F, {jetEtaAxis}});
117116
registry.add("h_jet_phi", "jet phi;#phi; counts", {HistType::kTH1F, {phiAxis}});
@@ -182,7 +181,7 @@ struct JetSpectraCharged {
182181
registry.add("h2_jet_pt_part_jet_area_part_rhoareasubtracted", "part jet #it{p}_{T,jet} vs. Area_{jet}; #it{p}_{T,jet}^{part} (GeV/#it{c}); Area_{jet}^{part}", {HistType::kTH2F, {jetPtAxisRhoAreaSub, {150, 0., 1.5}}});
183182
registry.add("h2_jet_pt_part_jet_ntracks_part_rhoareasubtracted", "part jet #it{p}_{T,jet} vs. N_{jet tracks}; #it{p}_{T,jet}^{part} (GeV/#it{c}); N_{jet, tracks}{part}", {HistType::kTH2F, {jetPtAxisRhoAreaSub, {200, -0.5, 199.5}}});
184183
registry.add("h3_jet_pt_eta_phi_part_rhoareasubtracted", "part jet pt vs. eta vs.phi", {HistType::kTH3F, {jetPtAxisRhoAreaSub, jetEtaAxis, phiAxis}});
185-
}
184+
}
186185

187186
if (doprocessEvtWiseConstSubJetsData || doprocessEvtWiseConstSubJetsMCD) {
188187
registry.add("h2_centrality_jet_pt_eventwiseconstituentsubtracted", "centrality vs. jet pT;centrality;#it{p}_{T,jet} (GeV/#it{c});entries", {HistType::kTH2F, {centralityAxis, jetPtAxis}});
@@ -325,7 +324,7 @@ struct JetSpectraCharged {
325324
registry.fill(HIST("h3_jet_pt_eta_phi_part"), jet.pt(), jet.eta(), jet.phi(), weight);
326325
registry.fill(HIST("h2_jet_pt_part_jet_area_part"), jet.pt(), jet.area(), weight);
327326
registry.fill(HIST("h2_jet_pt_part_jet_ntracks_part"), jet.pt(), jet.tracksIds().size(), weight);
328-
}
327+
}
329328

330329
for (const auto& constituent : jet.template tracks_as<aod::JetParticles>()) {
331330
registry.fill(HIST("h2_jet_pt_part_track_pt_part"), jet.pt(), constituent.pt(), weight);
@@ -337,14 +336,14 @@ struct JetSpectraCharged {
337336
{
338337
if (jet.r() == round(selectedJetsRadius * 100.0f)) {
339338
// fill mcp jet histograms
340-
double jetcorrpt = jet.pt() - (rho * jet.area());
341-
registry.fill(HIST("h_jet_pt_part_rhoareasubtracted"), jetcorrpt, weight);
342-
registry.fill(HIST("h3_jet_pt_eta_phi_part_rhoareasubtracted"), jetcorrpt, jet.eta(), jet.phi(), weight);
343-
if (jetcorrpt > 0) {
344-
registry.fill(HIST("h_jet_eta_part_rhoareasubtracted"), jet.eta(), weight);
345-
registry.fill(HIST("h_jet_phi_part_rhoareasubtracted"), jet.phi(), weight);
346-
registry.fill(HIST("h2_jet_pt_part_jet_area_part_rhoareasubtracted"), jetcorrpt, jet.area(), weight);
347-
registry.fill(HIST("h2_jet_pt_part_jet_ntracks_part_rhoareasubtracted"), jetcorrpt, jet.tracksIds().size(), weight);
339+
double jetcorrpt = jet.pt() - (rho * jet.area());
340+
registry.fill(HIST("h_jet_pt_part_rhoareasubtracted"), jetcorrpt, weight);
341+
registry.fill(HIST("h3_jet_pt_eta_phi_part_rhoareasubtracted"), jetcorrpt, jet.eta(), jet.phi(), weight);
342+
if (jetcorrpt > 0) {
343+
registry.fill(HIST("h_jet_eta_part_rhoareasubtracted"), jet.eta(), weight);
344+
registry.fill(HIST("h_jet_phi_part_rhoareasubtracted"), jet.phi(), weight);
345+
registry.fill(HIST("h2_jet_pt_part_jet_area_part_rhoareasubtracted"), jetcorrpt, jet.area(), weight);
346+
registry.fill(HIST("h2_jet_pt_part_jet_ntracks_part_rhoareasubtracted"), jetcorrpt, jet.tracksIds().size(), weight);
348347
}
349348
}
350349
}
@@ -516,19 +515,19 @@ struct JetSpectraCharged {
516515
aod::JetTracks const&)
517516
{
518517
if (!jetderiveddatautilities::selectCollision(collision, eventSelectionBits)) {
519-
return;
518+
return;
520519
}
521520
if (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange()) {
522-
return;
521+
return;
523522
}
524523
for (auto const& jet : jets) {
525-
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
526-
continue;
527-
}
528-
if (!isAcceptedJet<aod::JetTracks>(jet)) {
529-
continue;
530-
}
531-
fillJetHistograms(jet, collision.centrality());
524+
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
525+
continue;
526+
}
527+
if (!isAcceptedJet<aod::JetTracks>(jet)) {
528+
continue;
529+
}
530+
fillJetHistograms(jet, collision.centrality());
532531
}
533532
}
534533
PROCESS_SWITCH(JetSpectraCharged, processSpectraData, "jet spectra for Data", false);
@@ -538,19 +537,19 @@ struct JetSpectraCharged {
538537
aod::JetTracks const&)
539538
{
540539
if (!jetderiveddatautilities::selectCollision(collision, eventSelectionBits)) {
541-
return;
540+
return;
542541
}
543542
if (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange()) {
544-
return;
543+
return;
545544
}
546545
for (auto const& jet : jets) {
547-
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
548-
continue;
549-
}
550-
if (!isAcceptedJet<aod::JetTracks>(jet)) {
551-
continue;
552-
}
553-
fillJetHistograms(jet, collision.centrality());
546+
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
547+
continue;
548+
}
549+
if (!isAcceptedJet<aod::JetTracks>(jet)) {
550+
continue;
551+
}
552+
fillJetHistograms(jet, collision.centrality());
554553
}
555554
}
556555
PROCESS_SWITCH(JetSpectraCharged, processSpectraMCD, "jet spectra for MCD", false);
@@ -633,61 +632,61 @@ struct JetSpectraCharged {
633632
soa::Join<aod::ChargedMCParticleLevelJets, aod::ChargedMCParticleLevelJetConstituents> const& jets,
634633
aod::JetParticles const&)
635634
{
636-
registry.fill(HIST("h_mcColl_counts_areasub"), 0.5); //all mc collisions
635+
registry.fill(HIST("h_mcColl_counts_areasub"), 0.5);
637636
if (std::abs(mccollision.posZ()) > vertexZCut) {
638637
return;
639638
}
640-
registry.fill(HIST("h_mcColl_counts_areasub"), 1.5); //cut on mccollision.posZ()
639+
registry.fill(HIST("h_mcColl_counts_areasub"), 1.5);
641640
if (collisions.size() < 1) {
642641
return;
643642
}
644-
registry.fill(HIST("h_mcColl_counts_areasub"), 2.5); //mcCollisions with at least one resonstructed collision
643+
registry.fill(HIST("h_mcColl_counts_areasub"), 2.5);
645644
if (acceptSplitCollisions == 0 && collisions.size() > 1) {
646645
return;
647646
}
648-
registry.fill(HIST("h_mcColl_counts_areasub"), 3.5); // split mcCollisions condition
647+
registry.fill(HIST("h_mcColl_counts_areasub"), 3.5);
649648

650649
bool hasSel8Coll = false;
651650
bool centralityIsGood = false;
652651
bool occupancyIsGood = false;
653-
if (acceptSplitCollisions == 2) {
654-
if (jetderiveddatautilities::selectCollision(collisions.begin(), eventSelectionBits)) { // Skipping MC events that have not a single selected reconstructed collision ; effect unclear if mcColl is split
652+
if (acceptSplitCollisions == 2) {
653+
if (jetderiveddatautilities::selectCollision(collisions.begin(), eventSelectionBits)) {
655654
hasSel8Coll = true;
656655
}
657-
if ((centralityMin < collisions.begin().centrality()) && (collisions.begin().centrality() < centralityMax)) { // effect unclear if mcColl is split
656+
if ((centralityMin < collisions.begin().centrality()) && (collisions.begin().centrality() < centralityMax)) {
658657
centralityIsGood = true;
659658
}
660659
if ((trackOccupancyInTimeRangeMin < collisions.begin().trackOccupancyInTimeRange()) && (collisions.begin().trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMax)) {
661660
occupancyIsGood = true;
662661
}
663662
}
664-
else { // check that at least one of the reconstructed collisions passes the checks
663+
else {
665664
for (auto& collision : collisions) {
666-
if (jetderiveddatautilities::selectCollision(collision, eventSelectionBits)) { // Skipping MC events that have not a single selected reconstructed collision ; effect unclear if mcColl is split
665+
if (jetderiveddatautilities::selectCollision(collision, eventSelectionBits)) {
667666
hasSel8Coll = true;
668667
}
669-
if ((centralityMin < collision.centrality()) && (collision.centrality() < centralityMax)) { // effect unclear if mcColl is split
668+
if ((centralityMin < collision.centrality()) && (collision.centrality() < centralityMax)) {
670669
centralityIsGood = true;
671670
}
672671
if ((trackOccupancyInTimeRangeMin < collision.trackOccupancyInTimeRange()) && (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMax)) {
673-
occupancyIsGood = true;
672+
occupancyIsGood = true;
674673
}
675674
}
676675
}
677676
if (!hasSel8Coll) {
678677
return;
679678
}
680-
registry.fill(HIST("h_mcColl_counts_areasub"), 4.5); // at least one of the reconstructed collisions associated with this mcCollision is selected
679+
registry.fill(HIST("h_mcColl_counts_areasub"), 4.5);
681680

682681
if (!centralityIsGood) {
683682
return;
684683
}
685-
registry.fill(HIST("h_mcColl_counts_areasub"), 5.5); // at least one of the reconstructed collisions associated with this mcCollision is selected with regard to centrality
684+
registry.fill(HIST("h_mcColl_counts_areasub"), 5.5);
686685

687686
if (!occupancyIsGood) {
688687
return;
689688
}
690-
registry.fill(HIST("h_mcColl_counts_areasub"), 6.5); // at least one of the reconstructed collisions associated with this mcCollision is selected with regard to centrality
689+
registry.fill(HIST("h_mcColl_counts_areasub"), 6.5);
691690
registry.fill(HIST("h_mcColl_rho"), mccollision.rho());
692691
for (auto const& jet : jets) {
693692
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
@@ -706,15 +705,15 @@ struct JetSpectraCharged {
706705
soa::Join<aod::ChargedMCParticleLevelJets, aod::ChargedMCParticleLevelJetConstituents> const& jets,
707706
aod::JetParticles const&)
708707
{
709-
registry.fill(HIST("h_mcColl_counts"), 0.5); //all mc collisions
708+
registry.fill(HIST("h_mcColl_counts"), 0.5);
710709
if (std::abs(mccollision.posZ()) > vertexZCut) {
711710
return;
712711
}
713-
registry.fill(HIST("h_mcColl_counts"), 1.5); //cut on mccollision.posZ()
712+
registry.fill(HIST("h_mcColl_counts"), 1.5);
714713
if (collisions.size() < 1) {
715714
return;
716715
}
717-
registry.fill(HIST("h_mcColl_counts"), 2.5); //mcCollisions with at least one resonstructed collision
716+
registry.fill(HIST("h_mcColl_counts"), 2.5);
718717

719718
bool hasSel8Coll = false;
720719
bool centralityIsGood = false;
@@ -733,44 +732,43 @@ struct JetSpectraCharged {
733732
if (!hasSel8Coll) {
734733
return;
735734
}
736-
registry.fill(HIST("h_mcColl_counts"), 3.5); // at least one of the reconstructed collisions associated with this mcCollision is selected
735+
registry.fill(HIST("h_mcColl_counts"), 3.5);
737736
if (!centralityIsGood) {
738737
return;
739738
}
740-
registry.fill(HIST("h_mcColl_counts"), 4.5); // at least one of the reconstructed collisions associated with this mcCollision is selected with regard to centrality
739+
registry.fill(HIST("h_mcColl_counts"), 4.5);
741740
if (!occupancyIsGood) {
742741
return;
743742
}
744-
registry.fill(HIST("h_mcColl_counts"), 5.5); // at least one of the reconstructed collisions associated with this mcCollision is selected with regard to centrality
743+
registry.fill(HIST("h_mcColl_counts"), 5.5);
745744
registry.fill(HIST("h_mc_zvertex"), mccollision.posZ());
746745

747746
for (auto const& jet : jets) {
748747
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
749-
continue;
748+
continue;
750749
}
751750
if (!isAcceptedJet<aod::JetParticles>(jet)) {
752-
continue;
751+
continue;
753752
}
754753
fillMCPHistograms(jet);
755754
}
756755
}
757756
PROCESS_SWITCH(JetSpectraCharged, processSpectraMCP, "jet spectra for MC particle level", false);
758757

759-
760758
void processSpectraMCPWeighted(aod::JetMcCollision const& mccollision,
761759
soa::SmallGroups<aod::JetCollisionsMCD> const& collisions,
762760
soa::Join<aod::ChargedMCParticleLevelJets, aod::ChargedMCParticleLevelJetConstituents, aod::ChargedMCParticleLevelJetEventWeights> const& jets,
763761
aod::JetParticles const&)
764762
{
765-
registry.fill(HIST("h_mcColl_counts"), 0.5); //all mc collisions
763+
registry.fill(HIST("h_mcColl_counts"), 0.5);
766764
if (std::abs(mccollision.posZ()) > vertexZCut) {
767765
return;
768766
}
769-
registry.fill(HIST("h_mcColl_counts"), 1.5); //cut on mccollision.posZ()
767+
registry.fill(HIST("h_mcColl_counts"), 1.5);
770768
if (collisions.size() < 1) {
771769
return;
772770
}
773-
registry.fill(HIST("h_mcColl_counts"), 2.5); //mcCollisions with at least one resonstructed collision
771+
registry.fill(HIST("h_mcColl_counts"), 2.5);
774772

775773
bool hasSel8Coll = false;
776774
for (auto& collision : collisions) {
@@ -781,7 +779,7 @@ struct JetSpectraCharged {
781779
if (!hasSel8Coll) {
782780
return;
783781
}
784-
registry.fill(HIST("h_mcColl_counts"), 3.5); // at least one of the reconstructed collisions associated with this mcCollision is selected
782+
registry.fill(HIST("h_mcColl_counts"), 3.5);
785783
registry.fill(HIST("h_mc_zvertex"), mccollision.posZ());
786784

787785
for (auto const& jet : jets) {

0 commit comments

Comments
 (0)