Skip to content

Commit 533ef62

Browse files
SuJeong-JiSuJeong-Ji
authored
Merge pull request #22 from alibuild/alibot-cleanup-13364
Please consider the following formatting changes to #13364
2 parents 17f4852 + 0a25f15 commit 533ef62

File tree

1 file changed

+129
-129
lines changed

1 file changed

+129
-129
lines changed

‎PWGLF/Tasks/Resonances/chk892pp.cxx‎

Lines changed: 129 additions & 129 deletions
Original file line numberDiff line numberDiff line change
@@ -109,7 +109,7 @@ struct Chk892pp {
109109

110110
SliceCache cache;
111111
Preslice<aod::Tracks> perCollision = aod::track::collisionId;
112-
Preslice<aod::V0s> perCollisionV0 = aod::v0data::collisionId;
112+
Preslice<aod::V0s> perCollisionV0 = aod::v0data::collisionId;
113113

114114
using EventCandidates = soa::Join<aod::Collisions, aod::EvSels, aod::FT0Mults, aod::FV0Mults, aod::TPCMults, aod::CentFV0As, aod::CentFT0Ms, aod::CentFT0Cs, aod::CentFT0As, aod::Mults, aod::PVMults>;
115115
using TrackCandidates = soa::Join<aod::FullTracks, aod::TracksExtra, aod::TracksDCA, aod::TrackSelection, aod::TrackSelectionExtension, aod::pidTPCFullPi, aod::pidTOFFullPi>;
@@ -172,7 +172,7 @@ struct Chk892pp {
172172
} EventCuts;
173173
RCTFlagsChecker rctChecker;
174174

175-
Configurable<bool> cfgFillQAPlots{"cfgFillQAPlots", false, "Fill QA plots"};
175+
Configurable<bool> cfgFillQAPlots{"cfgFillQAPlots", false, "Fill QA plots"};
176176
Configurable<int> cfgCentEst{"cfgCentEst", 2, "Centrality estimator, 1: FT0C, 2: FT0M"};
177177

178178
/// PID Selections, pion
@@ -923,7 +923,7 @@ struct Chk892pp {
923923
if (std::abs(yreco) > KstarCuts.cfgKstarMaxRap)
924924
continue;
925925
histos.fill(HIST("EffKstar/recoKstar"), ptreco, lCentrality);
926-
}
926+
}
927927
}
928928
} // effKstarProcessReco
929929

@@ -947,37 +947,37 @@ struct Chk892pp {
947947
auto trkbNSigmaPiTOF = (istrkbhasTOF) ? bTrack.tofNSigmaPi() : -999.;
948948

949949
if constexpr (!IsMix) {
950-
if (cfgFillQAPlots) {
951-
// Bachelor pion QA plots
952-
histos.fill(HIST("QA/before/trkbpionTPCPID"), trkbpt, trkbNSigmaPiTPC);
953-
if (istrkbhasTOF) {
954-
histos.fill(HIST("QA/before/trkbpionTOFPID"), trkbpt, trkbNSigmaPiTOF);
955-
histos.fill(HIST("QA/before/trkbpionTPCTOFPID"), trkbNSigmaPiTPC, trkbNSigmaPiTOF);
956-
}
957-
histos.fill(HIST("QA/before/trkbpionpT"), trkbpt);
958-
histos.fill(HIST("QA/before/trkbpionDCAxy"), bTrack.dcaXY());
959-
histos.fill(HIST("QA/before/trkbpionDCAz"), bTrack.dcaZ());
960-
}
961-
}
950+
if (cfgFillQAPlots) {
951+
// Bachelor pion QA plots
952+
histos.fill(HIST("QA/before/trkbpionTPCPID"), trkbpt, trkbNSigmaPiTPC);
953+
if (istrkbhasTOF) {
954+
histos.fill(HIST("QA/before/trkbpionTOFPID"), trkbpt, trkbNSigmaPiTOF);
955+
histos.fill(HIST("QA/before/trkbpionTPCTOFPID"), trkbNSigmaPiTPC, trkbNSigmaPiTOF);
956+
}
957+
histos.fill(HIST("QA/before/trkbpionpT"), trkbpt);
958+
histos.fill(HIST("QA/before/trkbpionDCAxy"), bTrack.dcaXY());
959+
histos.fill(HIST("QA/before/trkbpionDCAz"), bTrack.dcaZ());
960+
}
961+
}
962962

963963
if (!trackCut(bTrack))
964964
continue;
965965
if (!selectionPIDPion(bTrack))
966966
continue;
967967

968968
if constexpr (!IsMix) {
969-
if (cfgFillQAPlots) {
970-
// Bachelor pion QA plots after applying cuts
971-
histos.fill(HIST("QA/after/trkbpionTPCPID"), trkbpt, trkbNSigmaPiTPC);
972-
if (istrkbhasTOF) {
973-
histos.fill(HIST("QA/after/trkbpionTOFPID"), trkbpt, trkbNSigmaPiTOF);
974-
histos.fill(HIST("QA/after/trkbpionTPCTOFPID"), trkbNSigmaPiTPC, trkbNSigmaPiTOF);
975-
}
976-
histos.fill(HIST("QA/after/trkbpionpT"), trkbpt);
977-
histos.fill(HIST("QA/after/trkbpionDCAxy"), bTrack.dcaXY());
978-
histos.fill(HIST("QA/after/trkbpionDCAz"), bTrack.dcaZ());
979-
}
980-
}
969+
if (cfgFillQAPlots) {
970+
// Bachelor pion QA plots after applying cuts
971+
histos.fill(HIST("QA/after/trkbpionTPCPID"), trkbpt, trkbNSigmaPiTPC);
972+
if (istrkbhasTOF) {
973+
histos.fill(HIST("QA/after/trkbpionTOFPID"), trkbpt, trkbNSigmaPiTOF);
974+
histos.fill(HIST("QA/after/trkbpionTPCTOFPID"), trkbNSigmaPiTPC, trkbNSigmaPiTOF);
975+
}
976+
histos.fill(HIST("QA/after/trkbpionpT"), trkbpt);
977+
histos.fill(HIST("QA/after/trkbpionDCAxy"), bTrack.dcaXY());
978+
histos.fill(HIST("QA/after/trkbpionDCAz"), bTrack.dcaZ());
979+
}
980+
}
981981
trackIndicies.push_back(bTrack.index());
982982
}
983983

@@ -1013,39 +1013,39 @@ struct Chk892pp {
10131013
lResoSecondary = LorentzVectorSetXYZM(k0sCand.px(), k0sCand.py(), k0sCand.pz(), trkkMass);
10141014

10151015
if constexpr (!IsMix) {
1016-
if (cfgFillQAPlots) {
1017-
// Seconddary QA plots
1018-
histos.fill(HIST("QA/before/trkppionTPCPID"), trkppt, trkpNSigmaPiTPC);
1019-
if (istrkphasTOF) {
1020-
histos.fill(HIST("QA/before/trkppionTOFPID"), trkppt, trkpNSigmaPiTOF);
1021-
histos.fill(HIST("QA/before/trkppionTPCTOFPID"), trkpNSigmaPiTPC, trkpNSigmaPiTOF);
1022-
}
1023-
histos.fill(HIST("QA/before/trkppionpT"), trkppt);
1024-
histos.fill(HIST("QA/before/trkppionDCAxy"), posDauTrack.dcaXY());
1025-
histos.fill(HIST("QA/before/trkppionDCAz"), posDauTrack.dcaZ());
1026-
1027-
histos.fill(HIST("QA/before/trknpionTPCPID"), trknpt, trknNSigmaPiTPC);
1028-
if (istrknhasTOF) {
1029-
histos.fill(HIST("QA/before/trknpionTOFPID"), trknpt, trknNSigmaPiTOF);
1030-
histos.fill(HIST("QA/before/trknpionTPCTOFPID"), trknNSigmaPiTPC, trknNSigmaPiTOF);
1031-
}
1032-
histos.fill(HIST("QA/before/trknpionpT"), trknpt);
1033-
histos.fill(HIST("QA/before/trknpionDCAxy"), negDauTrack.dcaXY());
1034-
histos.fill(HIST("QA/before/trknpionDCAz"), negDauTrack.dcaZ());
1035-
1036-
histos.fill(HIST("QA/before/hDauDCASecondary"), trkkDauDCA);
1037-
histos.fill(HIST("QA/before/hDauPosDCAtoPVSecondary"), trkkDauDCAPostoPV);
1038-
histos.fill(HIST("QA/before/hDauNegDCAtoPVSecondary"), trkkDauDCANegtoPV);
1039-
1040-
histos.fill(HIST("QA/before/hpT_Secondary"), trkkpt);
1041-
histos.fill(HIST("QA/before/hy_Secondary"), trkky);
1042-
histos.fill(HIST("QA/before/hRadiusSecondary"), trkkRadius);
1043-
histos.fill(HIST("QA/before/hDCAtoPVSecondary"), trkkDCAtoPV);
1044-
histos.fill(HIST("QA/before/hCPASecondary"), trkkCPA);
1045-
histos.fill(HIST("QA/before/hPropTauSecondary"), trkkPropTau);
1046-
histos.fill(HIST("QA/before/hInvmassSecondary"), trkkMass);
1047-
}
1048-
}
1016+
if (cfgFillQAPlots) {
1017+
// Seconddary QA plots
1018+
histos.fill(HIST("QA/before/trkppionTPCPID"), trkppt, trkpNSigmaPiTPC);
1019+
if (istrkphasTOF) {
1020+
histos.fill(HIST("QA/before/trkppionTOFPID"), trkppt, trkpNSigmaPiTOF);
1021+
histos.fill(HIST("QA/before/trkppionTPCTOFPID"), trkpNSigmaPiTPC, trkpNSigmaPiTOF);
1022+
}
1023+
histos.fill(HIST("QA/before/trkppionpT"), trkppt);
1024+
histos.fill(HIST("QA/before/trkppionDCAxy"), posDauTrack.dcaXY());
1025+
histos.fill(HIST("QA/before/trkppionDCAz"), posDauTrack.dcaZ());
1026+
1027+
histos.fill(HIST("QA/before/trknpionTPCPID"), trknpt, trknNSigmaPiTPC);
1028+
if (istrknhasTOF) {
1029+
histos.fill(HIST("QA/before/trknpionTOFPID"), trknpt, trknNSigmaPiTOF);
1030+
histos.fill(HIST("QA/before/trknpionTPCTOFPID"), trknNSigmaPiTPC, trknNSigmaPiTOF);
1031+
}
1032+
histos.fill(HIST("QA/before/trknpionpT"), trknpt);
1033+
histos.fill(HIST("QA/before/trknpionDCAxy"), negDauTrack.dcaXY());
1034+
histos.fill(HIST("QA/before/trknpionDCAz"), negDauTrack.dcaZ());
1035+
1036+
histos.fill(HIST("QA/before/hDauDCASecondary"), trkkDauDCA);
1037+
histos.fill(HIST("QA/before/hDauPosDCAtoPVSecondary"), trkkDauDCAPostoPV);
1038+
histos.fill(HIST("QA/before/hDauNegDCAtoPVSecondary"), trkkDauDCANegtoPV);
1039+
1040+
histos.fill(HIST("QA/before/hpT_Secondary"), trkkpt);
1041+
histos.fill(HIST("QA/before/hy_Secondary"), trkky);
1042+
histos.fill(HIST("QA/before/hRadiusSecondary"), trkkRadius);
1043+
histos.fill(HIST("QA/before/hDCAtoPVSecondary"), trkkDCAtoPV);
1044+
histos.fill(HIST("QA/before/hCPASecondary"), trkkCPA);
1045+
histos.fill(HIST("QA/before/hPropTauSecondary"), trkkPropTau);
1046+
histos.fill(HIST("QA/before/hInvmassSecondary"), trkkMass);
1047+
}
1048+
}
10491049

10501050
if (!SecondaryCuts.cfgByPassDauPIDSelection && !selectionPIDPion(posDauTrack))
10511051
continue;
@@ -1055,46 +1055,46 @@ struct Chk892pp {
10551055
continue;
10561056

10571057
if constexpr (!IsMix) {
1058-
if (cfgFillQAPlots) {
1059-
// Seconddary QA plots after applying cuts
1060-
1061-
histos.fill(HIST("QA/after/trkppionTPCPID"), trkppt, trkpNSigmaPiTPC);
1062-
if (istrkphasTOF) {
1063-
histos.fill(HIST("QA/after/trkppionTOFPID"), trkppt, trkpNSigmaPiTOF);
1064-
histos.fill(HIST("QA/after/trkppionTPCTOFPID"), trkpNSigmaPiTPC, trkpNSigmaPiTOF);
1065-
}
1066-
histos.fill(HIST("QA/after/trkppionpT"), trkppt);
1067-
histos.fill(HIST("QA/after/trkppionDCAxy"), posDauTrack.dcaXY());
1068-
histos.fill(HIST("QA/after/trkppionDCAz"), posDauTrack.dcaZ());
1069-
1070-
histos.fill(HIST("QA/after/trknpionTPCPID"), trknpt, trknNSigmaPiTPC);
1071-
if (istrknhasTOF) {
1072-
histos.fill(HIST("QA/after/trknpionTOFPID"), trknpt, trknNSigmaPiTOF);
1073-
histos.fill(HIST("QA/after/trknpionTPCTOFPID"), trknNSigmaPiTPC, trknNSigmaPiTOF);
1074-
}
1075-
histos.fill(HIST("QA/after/trknpionpT"), trknpt);
1076-
histos.fill(HIST("QA/after/trknpionDCAxy"), negDauTrack.dcaXY());
1077-
histos.fill(HIST("QA/after/trknpionDCAz"), negDauTrack.dcaZ());
1078-
1079-
histos.fill(HIST("QA/after/hDauDCASecondary"), trkkDauDCA);
1080-
histos.fill(HIST("QA/after/hDauPosDCAtoPVSecondary"), trkkDauDCAPostoPV);
1081-
histos.fill(HIST("QA/after/hDauNegDCAtoPVSecondary"), trkkDauDCANegtoPV);
1082-
1083-
histos.fill(HIST("QA/after/hpT_Secondary"), trkkpt);
1084-
histos.fill(HIST("QA/after/hy_Secondary"), trkky);
1085-
histos.fill(HIST("QA/after/hRadiusSecondary"), trkkRadius);
1086-
histos.fill(HIST("QA/after/hDCAtoPVSecondary"), trkkDCAtoPV);
1087-
histos.fill(HIST("QA/after/hCPASecondary"), trkkCPA);
1088-
histos.fill(HIST("QA/after/hPropTauSecondary"), trkkPropTau);
1089-
histos.fill(HIST("QA/after/hInvmassSecondary"), trkkMass);
1090-
}
1091-
histos.fill(HIST("hInvmass_K0s"), lCentrality, lResoSecondary.Pt(), lResoSecondary.M());
1092-
}
1093-
k0sIndicies.push_back(k0sCand.index());
1094-
}
1058+
if (cfgFillQAPlots) {
1059+
// Seconddary QA plots after applying cuts
1060+
1061+
histos.fill(HIST("QA/after/trkppionTPCPID"), trkppt, trkpNSigmaPiTPC);
1062+
if (istrkphasTOF) {
1063+
histos.fill(HIST("QA/after/trkppionTOFPID"), trkppt, trkpNSigmaPiTOF);
1064+
histos.fill(HIST("QA/after/trkppionTPCTOFPID"), trkpNSigmaPiTPC, trkpNSigmaPiTOF);
1065+
}
1066+
histos.fill(HIST("QA/after/trkppionpT"), trkppt);
1067+
histos.fill(HIST("QA/after/trkppionDCAxy"), posDauTrack.dcaXY());
1068+
histos.fill(HIST("QA/after/trkppionDCAz"), posDauTrack.dcaZ());
1069+
1070+
histos.fill(HIST("QA/after/trknpionTPCPID"), trknpt, trknNSigmaPiTPC);
1071+
if (istrknhasTOF) {
1072+
histos.fill(HIST("QA/after/trknpionTOFPID"), trknpt, trknNSigmaPiTOF);
1073+
histos.fill(HIST("QA/after/trknpionTPCTOFPID"), trknNSigmaPiTPC, trknNSigmaPiTOF);
1074+
}
1075+
histos.fill(HIST("QA/after/trknpionpT"), trknpt);
1076+
histos.fill(HIST("QA/after/trknpionDCAxy"), negDauTrack.dcaXY());
1077+
histos.fill(HIST("QA/after/trknpionDCAz"), negDauTrack.dcaZ());
1078+
1079+
histos.fill(HIST("QA/after/hDauDCASecondary"), trkkDauDCA);
1080+
histos.fill(HIST("QA/after/hDauPosDCAtoPVSecondary"), trkkDauDCAPostoPV);
1081+
histos.fill(HIST("QA/after/hDauNegDCAtoPVSecondary"), trkkDauDCANegtoPV);
1082+
1083+
histos.fill(HIST("QA/after/hpT_Secondary"), trkkpt);
1084+
histos.fill(HIST("QA/after/hy_Secondary"), trkky);
1085+
histos.fill(HIST("QA/after/hRadiusSecondary"), trkkRadius);
1086+
histos.fill(HIST("QA/after/hDCAtoPVSecondary"), trkkDCAtoPV);
1087+
histos.fill(HIST("QA/after/hCPASecondary"), trkkCPA);
1088+
histos.fill(HIST("QA/after/hPropTauSecondary"), trkkPropTau);
1089+
histos.fill(HIST("QA/after/hInvmassSecondary"), trkkMass);
1090+
}
1091+
histos.fill(HIST("hInvmass_K0s"), lCentrality, lResoSecondary.Pt(), lResoSecondary.M());
1092+
}
1093+
k0sIndicies.push_back(k0sCand.index());
1094+
}
10951095

10961096
for (const auto& trackIndex : trackIndicies) {
1097-
auto bTrack = dTracks1.rawIteratorAt(trackIndex);
1097+
auto bTrack = dTracks1.rawIteratorAt(trackIndex);
10981098
for (const auto& k0sIndex : k0sIndicies) {
10991099
auto k0sCand = dTracks2.rawIteratorAt(k0sIndex);
11001100

@@ -1106,30 +1106,30 @@ struct Chk892pp {
11061106

11071107
// QA plots
11081108
if constexpr (!IsMix) {
1109-
if (cfgFillQAPlots) {
1110-
histos.fill(HIST("QA/before/KstarRapidity"), lResoKstar.Rapidity());
1111-
histos.fill(HIST("QA/before/kstarinvmass"), lResoKstar.M());
1112-
}
1113-
}
1109+
if (cfgFillQAPlots) {
1110+
histos.fill(HIST("QA/before/KstarRapidity"), lResoKstar.Rapidity());
1111+
histos.fill(HIST("QA/before/kstarinvmass"), lResoKstar.M());
1112+
}
1113+
}
11141114

11151115
if (lResoKstar.Rapidity() > KstarCuts.cfgKstarMaxRap || lResoKstar.Rapidity() < KstarCuts.cfgKstarMinRap)
11161116
continue;
11171117

11181118
if constexpr (!IsMix) {
11191119
unsigned int typeKstar = bTrack.sign() > 0 ? BinType::kKstarP : BinType::kKstarN;
1120-
if (cfgFillQAPlots) {
1120+
if (cfgFillQAPlots) {
11211121

1122-
histos.fill(HIST("QA/after/KstarRapidity"), lResoKstar.Rapidity());
1123-
histos.fill(HIST("QA/after/kstarinvmass"), lResoKstar.M());
1124-
}
1125-
histos.fill(HIST("hInvmass_Kstar"), typeKstar, lCentrality, lResoKstar.Pt(), lResoKstar.M());
1122+
histos.fill(HIST("QA/after/KstarRapidity"), lResoKstar.Rapidity());
1123+
histos.fill(HIST("QA/after/kstarinvmass"), lResoKstar.M());
1124+
}
1125+
histos.fill(HIST("hInvmass_Kstar"), typeKstar, lCentrality, lResoKstar.Pt(), lResoKstar.M());
11261126

11271127
if (BkgEstimationConfig.cfgFillRotBkg) {
11281128
for (int i = 0; i < BkgEstimationConfig.cfgNrotBkg; i++) {
11291129
auto lRotAngle = BkgEstimationConfig.cfgMinRot + i * ((BkgEstimationConfig.cfgMaxRot - BkgEstimationConfig.cfgMinRot) / (BkgEstimationConfig.cfgNrotBkg - 1));
1130-
if (cfgFillQAPlots) {
1131-
histos.fill(HIST("QA/RotBkg/hRotBkg"), lRotAngle);
1132-
}
1130+
if (cfgFillQAPlots) {
1131+
histos.fill(HIST("QA/RotBkg/hRotBkg"), lRotAngle);
1132+
}
11331133
if (BkgEstimationConfig.cfgRotPion) {
11341134
lDaughterRot = lDecayDaughter_bach;
11351135
// lDaughterRot.RotateZ(lRotAngle);
@@ -1150,7 +1150,7 @@ struct Chk892pp {
11501150
}
11511151
}
11521152
} // IsMix
1153-
} // K0scand
1153+
} // K0scand
11541154
} // bTrack
11551155

11561156
count++;
@@ -1186,32 +1186,32 @@ struct Chk892pp {
11861186
MCEventCandidates const& events)
11871187
{
11881188
buildAllowedMcIds(events);
1189-
effK0sProcessGen(mcpart);
1190-
effK0sProcessReco(v0s);
1191-
effKstarProcessGen(mcpart);
1192-
effKstarProcessReco(v0s, tracks);
1193-
}
1189+
effK0sProcessGen(mcpart);
1190+
effK0sProcessReco(v0s);
1191+
effKstarProcessGen(mcpart);
1192+
effKstarProcessReco(v0s, tracks);
1193+
}
11941194
PROCESS_SWITCH(Chk892pp, processMC, "Process Event for MC", false);
11951195

1196-
void processMCQA(MCEventCandidates::iterator const& collision,
1197-
MCTrackCandidates const& tracks,
1198-
MCV0Candidates const& v0s,
1199-
aod::BCsWithTimestamps const&)
1200-
{
1201-
if (!colCuts.isSelected(collision))
1202-
return;
1203-
if (EventCuts.cfgEvtUseRCTFlagChecker && !rctChecker(collision))
1204-
return;
1205-
1206-
lCentrality = getCentrality(collision);
1196+
void processMCQA(MCEventCandidates::iterator const& collision,
1197+
MCTrackCandidates const& tracks,
1198+
MCV0Candidates const& v0s,
1199+
aod::BCsWithTimestamps const&)
1200+
{
1201+
if (!colCuts.isSelected(collision))
1202+
return;
1203+
if (EventCuts.cfgEvtUseRCTFlagChecker && !rctChecker(collision))
1204+
return;
1205+
1206+
lCentrality = getCentrality(collision);
12071207
if (lCentrality < EventCuts.cfgEventCentralityMin || lCentrality > EventCuts.cfgEventCentralityMax)
12081208
return;
12091209
if (!collision.isInelGt0())
12101210
return;
12111211
colCuts.fillQA(collision);
12121212

12131213
fillHistograms<true, false>(collision, tracks, v0s);
1214-
}
1214+
}
12151215
PROCESS_SWITCH(Chk892pp, processMCQA, "Process Event for MC and fill QA plots", true);
12161216
};
12171217

0 commit comments

Comments
 (0)