Skip to content

Commit 505a5c4

Browse files
alibuildalibuild
committed
Please consider the following formatting changes
1 parent c6e3541 commit 505a5c4

File tree

1 file changed

+75
-78
lines changed

1 file changed

+75
-78
lines changed

‎PWGLF/Tasks/Resonances/chk892pp.cxx‎

Lines changed: 75 additions & 78 deletions
Original file line numberDiff line numberDiff line change
@@ -84,7 +84,7 @@ using namespace o2::constants::physics;
8484
using namespace o2::aod::rctsel;
8585

8686
struct Chk892pp {
87-
enum BinType : unsigned int {
87+
enum BinType : unsigned int {
8888
kKstarP = 0,
8989
kKstarN,
9090
kKstarP_Mix,
@@ -115,21 +115,21 @@ struct Chk892pp {
115115
struct : ConfigurableGroup {
116116
Configurable<std::string> cfgURL{"cfgURL", "http://alice-ccdb.cern.ch", "Address of the CCDB to browse"};
117117
} CCDBConfig;
118-
// Configurable<int64_t> nolaterthan{"ccdb-no-later-than", std::chrono::duration_cast<std::chrono::milliseconds>(std::chrono::system_clock::now().time_since_epoch()).count(), "Latest acceptable timestamp of creation for the object"};
118+
// Configurable<int64_t> nolaterthan{"ccdb-no-later-than", std::chrono::duration_cast<std::chrono::milliseconds>(std::chrono::system_clock::now().time_since_epoch()).count(), "Latest acceptable timestamp of creation for the object"};
119119

120120
// Configurables
121121
struct : ConfigurableGroup {
122122
ConfigurableAxis cfgBinsPt{"cfgBinsPt", {VARIABLE_WIDTH, 0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.1, 2.2, 2.3, 2.4, 2.5, 2.6, 2.7, 2.8, 2.9, 3.0, 3.1, 3.2, 3.3, 3.4, 3.5, 3.6, 3.7, 3.8, 3.9, 4.0, 4.1, 4.2, 4.3, 4.4, 4.5, 4.6, 4.7, 4.8, 4.9, 5.0, 5.1, 5.2, 5.3, 5.4, 5.5, 5.6, 5.7, 5.8, 5.9, 6.0, 6.1, 6.2, 6.3, 6.4, 6.5, 6.6, 6.7, 6.8, 6.9, 7.0, 7.1, 7.2, 7.3, 7.4, 7.5, 7.6, 7.7, 7.8, 7.9, 8.0, 8.1, 8.2, 8.3, 8.4, 8.5, 8.6, 8.7, 8.8, 8.9, 9.0, 9.1, 9.2, 9.3, 9.4, 9.5, 9.6, 9.7, 9.8, 9.9, 10.0, 10.1, 10.2, 10.3, 10.4, 10.5, 10.6, 10.7, 10.8, 10.9, 11.0, 11.1, 11.2, 11.3, 11.4, 11.5, 11.6, 11.7, 11.8, 11.9, 12.0, 12.1, 12.2, 12.3, 12.4, 12.5, 12.6, 12.7, 12.8, 12.9, 13.0, 13.1, 13.2, 13.3, 13.4, 13.5, 13.6, 13.7, 13.8, 13.9, 14.0, 14.1, 14.2, 14.3, 14.4, 14.5, 14.6, 14.7, 14.8, 14.9, 15.0}, "Binning of the pT axis"};
123-
ConfigurableAxis cfgBinsPtQA{"cfgBinsPtQA", {VARIABLE_WIDTH, 0.0, 0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.4, 1.6, 1.8, 2.0, 2.2, 2.4, 2.6, 2.8, 3.0, 3.2, 3.4, 3.6, 3.8, 4.0, 4.2, 4.4, 4.6, 4.8, 5.0, 5.2, 5.4, 5.6, 5.8, 6.0, 6.2, 6.4, 6.6, 6.8, 7.0, 7.2, 7.4, 7.6, 7.8, 8.0, 8.2, 8.4, 8.6, 8.8, 9.0, 9.2, 9.4, 9.6, 9.8, 10.0}, "Binning of the pT axis"};
124-
ConfigurableAxis cfgBinsCent{"cfgBinsCent", {VARIABLE_WIDTH, 0.0, 1.0, 5.0, 10.0, 20.0, 30.0, 40.0, 50.0, 60.0, 70.0, 80.0, 90.0, 100.0, 110.0}, "Binning of the centrality axis"};
125-
ConfigurableAxis cfgBinsVtxZ{"cfgBinsVtxZ", {VARIABLE_WIDTH, -10.0, -9.0, -8.0, -7.0, -6.0, -5.0, -4.0, -3.0, -2.0, -1.0, 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0}, "Binning of the z-vertex axis"};
126-
Configurable<int> cNbinsDiv{"cNbinsDiv", 1, "Integer to divide the number of bins"};
123+
ConfigurableAxis cfgBinsPtQA{"cfgBinsPtQA", {VARIABLE_WIDTH, 0.0, 0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.4, 1.6, 1.8, 2.0, 2.2, 2.4, 2.6, 2.8, 3.0, 3.2, 3.4, 3.6, 3.8, 4.0, 4.2, 4.4, 4.6, 4.8, 5.0, 5.2, 5.4, 5.6, 5.8, 6.0, 6.2, 6.4, 6.6, 6.8, 7.0, 7.2, 7.4, 7.6, 7.8, 8.0, 8.2, 8.4, 8.6, 8.8, 9.0, 9.2, 9.4, 9.6, 9.8, 10.0}, "Binning of the pT axis"};
124+
ConfigurableAxis cfgBinsCent{"cfgBinsCent", {VARIABLE_WIDTH, 0.0, 1.0, 5.0, 10.0, 20.0, 30.0, 40.0, 50.0, 60.0, 70.0, 80.0, 90.0, 100.0, 110.0}, "Binning of the centrality axis"};
125+
ConfigurableAxis cfgBinsVtxZ{"cfgBinsVtxZ", {VARIABLE_WIDTH, -10.0, -9.0, -8.0, -7.0, -6.0, -5.0, -4.0, -3.0, -2.0, -1.0, 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0}, "Binning of the z-vertex axis"};
126+
Configurable<int> cNbinsDiv{"cNbinsDiv", 1, "Integer to divide the number of bins"};
127127
Configurable<int> cNbinsDivQA{"cNbinsDivQA", 1, "Integer to divide the number of bins for QA"};
128128
} AxisConfig;
129129

130130
/// Event cuts
131131
o2::analysis::CollisonCuts colCuts;
132-
struct : ConfigurableGroup {
132+
struct : ConfigurableGroup {
133133
Configurable<float> cfgEvtZvtx{"cfgEvtZvtx", 10.f, "Evt sel: Max. z-Vertex (cm)"};
134134
Configurable<int> cfgEvtOccupancyInTimeRangeMax{"cfgEvtOccupancyInTimeRangeMax", -1, "Evt sel: maximum track occupancy"};
135135
Configurable<int> cfgEvtOccupancyInTimeRangeMin{"cfgEvtOccupancyInTimeRangeMin", -1, "Evt sel: minimum track occupancy"};
@@ -143,7 +143,7 @@ struct Chk892pp {
143143
Configurable<bool> cfgEvtNoITSROBorderCut{"cfgEvtNoITSROBorderCut", false, "Evt sel: apply NoITSRO border cut"};
144144
Configurable<bool> cfgincludeCentralityMC{"cfgincludeCentralityMC", false, "Include centrality in MC"};
145145
Configurable<bool> cfgEvtCollInTimeRangeStandard{"cfgEvtCollInTimeRangeStandard", true, "Evt sel: apply NoCollInTimeRangeStandard"};
146-
Configurable<float> cfgEventCentralityMin{"cfgEventCentralityMin", 0.0f, "Event sel: minimum centrality"};
146+
Configurable<float> cfgEventCentralityMin{"cfgEventCentralityMin", 0.0f, "Event sel: minimum centrality"};
147147
Configurable<float> cfgEventCentralityMax{"cfgEventCentralityMax", 80.0f, "Event sel: maximum centrality"};
148148
Configurable<bool> cfgEvtUseRCTFlagChecker{"cfgEvtUseRCTFlagChecker", false, "Evt sel: use RCT flag checker"};
149149
Configurable<std::string> cfgEvtRCTFlagCheckerLabel{"cfgEvtRCTFlagCheckerLabel", "CBT_hadronPID", "Evt sel: RCT flag checker label"};
@@ -161,19 +161,17 @@ struct Chk892pp {
161161
*/
162162
Configurable<int> cfgCentEst{"cfgCentEst", 1, "Centrality estimator, 1: FT0C, 2: FT0M"};
163163

164-
165-
166164
/// PID Selections, pion
167165
struct : ConfigurableGroup {
168-
Configurable<bool> cfgTPConly{"cfgTPConly", true, "Use only TPC for PID"}; // bool
169-
Configurable<float> cfgMaxTPCnSigmaPion{"cfgMaxTPCnSigmaPion", 3.0, "TPC nSigma cut for Pion"}; // TPC
170-
Configurable<float> cfgMaxTOFnSigmaPion{"cfgMaxTOFnSigmaPion", 3.0, "TOF nSigma cut for Pion"}; // TOF
166+
Configurable<bool> cfgTPConly{"cfgTPConly", true, "Use only TPC for PID"}; // bool
167+
Configurable<float> cfgMaxTPCnSigmaPion{"cfgMaxTPCnSigmaPion", 3.0, "TPC nSigma cut for Pion"}; // TPC
168+
Configurable<float> cfgMaxTOFnSigmaPion{"cfgMaxTOFnSigmaPion", 3.0, "TOF nSigma cut for Pion"}; // TOF
171169
Configurable<float> cfgNsigmaCutCombinedPion{"cfgNsigmaCutCombinedPion", -999, "Combined nSigma cut for Pion"}; // Combined
172-
Configurable<bool> cfgTOFVeto{"cfgTOFVeto", true, "TOF Veto, if false, TOF is nessessary for PID selection"}; // TOF Veto
170+
Configurable<bool> cfgTOFVeto{"cfgTOFVeto", true, "TOF Veto, if false, TOF is nessessary for PID selection"}; // TOF Veto
173171
} PIDCuts;
174172

175173
// Track selections
176-
struct : ConfigurableGroup {
174+
struct : ConfigurableGroup {
177175
Configurable<float> cfgMinPtcut{"cfgMinPtcut", 0.15, "Track minium pt cut"};
178176
Configurable<float> cfgMaxEtacut{"cfgMaxEtacut", 0.8, "Track maximum eta cut"};
179177
Configurable<bool> cfgPrimaryTrack{"cfgPrimaryTrack", true, "Primary track selection"}; // kGoldenChi2 | kDCAxy | kDCAz
@@ -182,51 +180,51 @@ struct Chk892pp {
182180
Configurable<bool> cfgPVContributor{"cfgPVContributor", false, "PV contributor track selection"}; // PV Contriuibutor
183181

184182
Configurable<bool> cfgpTdepDCAxyCut{"cfgpTdepDCAxyCut", false, "pT-dependent DCAxy cut"};
185-
Configurable<int> cfgITScluster{"cfgITScluster", 0, "Number of ITS cluster"};
186-
Configurable<int> cfgTPCcluster{"cfgTPCcluster", 0, "Number of TPC cluster"};
187-
Configurable<float> cfgRatioTPCRowsOverFindableCls{"cfgRatioTPCRowsOverFindableCls", 0.0f, "TPC Crossed Rows to Findable Clusters"};
188-
Configurable<float> cfgITSChi2NCl{"cfgITSChi2NCl", 999.0, "ITS Chi2/NCl"};
189-
Configurable<float> cfgTPCChi2NCl{"cfgTPCChi2NCl", 999.0, "TPC Chi2/NCl"};
190-
Configurable<bool> cfgUseTPCRefit{"cfgUseTPCRefit", false, "Require TPC Refit"};
191-
Configurable<bool> cfgUseITSRefit{"cfgUseITSRefit", false, "Require ITS Refit"};
192-
Configurable<bool> cfgHasITS{"cfgHasITS", false, "Require ITS"};
193-
Configurable<bool> cfgHasTPC{"cfgHasTPC", false, "Require TPC"};
194-
Configurable<bool> cfgHasTOF{"cfgHasTOF", false, "Require TOF"};
195-
// DCAr to PV
196-
Configurable<float> cfgMaxbDCArToPVcut{"cfgMaxbDCArToPVcut", 0.1, "Track DCAr cut to PV Maximum"};
197-
// DCAz to PV
198-
Configurable<float> cfgMaxbDCAzToPVcut{"cfgMaxbDCAzToPVcut", 0.1, "Track DCAz cut to PV Maximum"};
199-
} TrackCuts;
183+
Configurable<int> cfgITScluster{"cfgITScluster", 0, "Number of ITS cluster"};
184+
Configurable<int> cfgTPCcluster{"cfgTPCcluster", 0, "Number of TPC cluster"};
185+
Configurable<float> cfgRatioTPCRowsOverFindableCls{"cfgRatioTPCRowsOverFindableCls", 0.0f, "TPC Crossed Rows to Findable Clusters"};
186+
Configurable<float> cfgITSChi2NCl{"cfgITSChi2NCl", 999.0, "ITS Chi2/NCl"};
187+
Configurable<float> cfgTPCChi2NCl{"cfgTPCChi2NCl", 999.0, "TPC Chi2/NCl"};
188+
Configurable<bool> cfgUseTPCRefit{"cfgUseTPCRefit", false, "Require TPC Refit"};
189+
Configurable<bool> cfgUseITSRefit{"cfgUseITSRefit", false, "Require ITS Refit"};
190+
Configurable<bool> cfgHasITS{"cfgHasITS", false, "Require ITS"};
191+
Configurable<bool> cfgHasTPC{"cfgHasTPC", false, "Require TPC"};
192+
Configurable<bool> cfgHasTOF{"cfgHasTOF", false, "Require TOF"};
193+
// DCAr to PV
194+
Configurable<float> cfgMaxbDCArToPVcut{"cfgMaxbDCArToPVcut", 0.1, "Track DCAr cut to PV Maximum"};
195+
// DCAz to PV
196+
Configurable<float> cfgMaxbDCAzToPVcut{"cfgMaxbDCAzToPVcut", 0.1, "Track DCAz cut to PV Maximum"};
197+
} TrackCuts;
200198

201199
// Secondary Selection
202-
struct : ConfigurableGroup {
200+
struct : ConfigurableGroup {
203201
Configurable<bool> cfgReturnFlag{"cfgReturnFlag", false, "Return Flag for debugging"};
204-
Configurable<bool> cfgSecondaryRequire{"cfgSecondaryRequire", true, "Secondary cuts on/off"};
205-
Configurable<bool> cfgSecondaryArmenterosCut{"cfgSecondaryArmenterosCut", true, "cut on Armenteros-Podolanski graph"};
202+
Configurable<bool> cfgSecondaryRequire{"cfgSecondaryRequire", true, "Secondary cuts on/off"};
203+
Configurable<bool> cfgSecondaryArmenterosCut{"cfgSecondaryArmenterosCut", true, "cut on Armenteros-Podolanski graph"};
206204
Configurable<bool> cfgSecondaryCrossMassHypothesisCut{"cfgSecondaryCrossMassHypothesisCut", false, "Apply cut based on the lambda mass hypothesis"};
207205

208-
Configurable<bool> cfgByPassDauPIDSelection{"cfgByPassDauPIDSelection", true, "Bypass Daughters PID selection"};
209-
Configurable<float> cfgSecondaryDauDCAMax{"cfgSecondaryDauDCAMax", 1., "Maximum DCA Secondary daughters to PV"};
210-
Configurable<float> cfgSecondaryDauPosDCAtoPVMin{"cfgSecondaryDauPosDCAtoPVMin", 0.0, "Minimum DCA Secondary positive daughters to PV"};
211-
Configurable<float> cfgSecondaryDauNegDCAtoPVMin{"cfgSecondaryDauNegDCAtoPVMin", 0.0, "Minimum DCA Secondary negative daughters to PV"};
206+
Configurable<bool> cfgByPassDauPIDSelection{"cfgByPassDauPIDSelection", true, "Bypass Daughters PID selection"};
207+
Configurable<float> cfgSecondaryDauDCAMax{"cfgSecondaryDauDCAMax", 1., "Maximum DCA Secondary daughters to PV"};
208+
Configurable<float> cfgSecondaryDauPosDCAtoPVMin{"cfgSecondaryDauPosDCAtoPVMin", 0.0, "Minimum DCA Secondary positive daughters to PV"};
209+
Configurable<float> cfgSecondaryDauNegDCAtoPVMin{"cfgSecondaryDauNegDCAtoPVMin", 0.0, "Minimum DCA Secondary negative daughters to PV"};
212210

213-
Configurable<float> cfgSecondaryPtMin{"cfgSecondaryPtMin", 0.f, "Minimum transverse momentum of Secondary"};
214-
Configurable<float> cfgSecondaryRapidityMax{"cfgSecondaryRapidityMax", 0.8, "Maximum rapidity of Secondary"};
215-
Configurable<float> cfgSecondaryRadiusMin{"cfgSecondaryRadiusMin", 1.2, "Minimum transverse radius of Secondary"};
211+
Configurable<float> cfgSecondaryPtMin{"cfgSecondaryPtMin", 0.f, "Minimum transverse momentum of Secondary"};
212+
Configurable<float> cfgSecondaryRapidityMax{"cfgSecondaryRapidityMax", 0.8, "Maximum rapidity of Secondary"};
213+
Configurable<float> cfgSecondaryRadiusMin{"cfgSecondaryRadiusMin", 1.2, "Minimum transverse radius of Secondary"};
216214
Configurable<float> cfgSecondaryRadiusMax{"cfgSecondaryRadiusMax", 999.9, "Maximum transverse radius of Secondary"};
217215
Configurable<float> cfgSecondaryCosPAMin{"cfgSecondaryCosPAMin", 0.995, "Mininum cosine pointing angle of Secondary"};
218-
Configurable<float> cfgSecondaryDCAtoPVMax{"cfgSecondaryDCAtoPVMax", 0.3, "Maximum DCA Secondary to PV"};
219-
Configurable<float> cfgSecondaryProperLifetimeMax{"cfgSecondaryProperLifetimeMax", 20, "Maximum Secondary Lifetime"};
220-
Configurable<float> cfgSecondaryparamArmenterosCut{"cfgSecondaryparamArmenterosCut", 0.2, "parameter for Armenteros Cut"};
221-
Configurable<float> cfgSecondaryMassWindow{"cfgSecondaryMassWindow", 0.03, "Secondary inv mass selciton window"};
216+
Configurable<float> cfgSecondaryDCAtoPVMax{"cfgSecondaryDCAtoPVMax", 0.3, "Maximum DCA Secondary to PV"};
217+
Configurable<float> cfgSecondaryProperLifetimeMax{"cfgSecondaryProperLifetimeMax", 20, "Maximum Secondary Lifetime"};
218+
Configurable<float> cfgSecondaryparamArmenterosCut{"cfgSecondaryparamArmenterosCut", 0.2, "parameter for Armenteros Cut"};
219+
Configurable<float> cfgSecondaryMassWindow{"cfgSecondaryMassWindow", 0.03, "Secondary inv mass selciton window"};
222220
Configurable<float> cfgSecondaryCrossMassCutWindow{"cfgSecondaryCrossMassCutWindow", 0.05, "Secondary inv mass selection window with (anti)lambda hypothesis"};
223221
} SecondaryCuts;
224222

225223
// K* selection
226-
struct : ConfigurableGroup {
224+
struct : ConfigurableGroup {
227225
Configurable<float> cfgKstarMaxRap{"cfgKstarMaxRap", 0.5, "Kstar maximum rapidity"};
228226
Configurable<float> cfgKstarMinRap{"cfgKstarMinRap", -0.5, "Kstar minimum rapidity"};
229-
} KstarCuts;
227+
} KstarCuts;
230228

231229
// Bkg estimation
232230
struct : ConfigurableGroup {
@@ -237,14 +235,13 @@ struct Chk892pp {
237235
Configurable<int> cfgNrotBkg{"cfgNrotBkg", 4, "Number of rotated copies (background) per each original candidate"};
238236
} BkgEstimationConfig;
239237

240-
241238
float lCentrality;
242239

243240
// PDG code
244241
int kPDGK0s = kK0Short;
245242
int kPDGK0 = kK0;
246243
int kKstarPlus = o2::constants::physics::Pdg::kKPlusStar892;
247-
//int kPiPlus = 211;
244+
// int kPiPlus = 211;
248245

249246
void init(o2::framework::InitContext&)
250247
{
@@ -259,7 +256,7 @@ struct Chk892pp {
259256
colCuts.setApplyPileupRejection(EventCuts.cfgEvtPileupRejection);
260257
colCuts.setApplyNoITSROBorderCut(EventCuts.cfgEvtNoITSROBorderCut);
261258
colCuts.setApplyCollInTimeRangeStandard(EventCuts.cfgEvtCollInTimeRangeStandard);
262-
colCuts.printCuts();
259+
colCuts.printCuts();
263260

264261
rctChecker.init(EventCuts.cfgEvtRCTFlagCheckerLabel, EventCuts.cfgEvtRCTFlagCheckerZDCCheck, EventCuts.cfgEvtRCTFlagCheckerLimitAcceptAsBad);
265262

@@ -276,7 +273,7 @@ struct Chk892pp {
276273
AxisSpec dcaxyAxis = {200, 0, 2, "DCA_{#it{xy}} (cm)"};
277274
AxisSpec dcazAxis = {200, 0, 2, "DCA_{#it{z}} (cm)"};
278275
AxisSpec yAxis = {100, -1, 1, "Rapidity"};
279-
AxisSpec invMassAxisK0s = {800 / AxisConfig.cNbinsDiv, 0.46, 0.54, "Invariant Mass (GeV/#it{c}^2)"}; // K0s ~497.611
276+
AxisSpec invMassAxisK0s = {800 / AxisConfig.cNbinsDiv, 0.46, 0.54, "Invariant Mass (GeV/#it{c}^2)"}; // K0s ~497.611
280277
AxisSpec invMassAxisReso = {900 / AxisConfig.cNbinsDiv, 0.5f, 1.4f, "Invariant Mass (GeV/#it{c}^2)"}; // chK(892) ~892
281278
AxisSpec pidQAAxis = {130 / AxisConfig.cNbinsDivQA, -6.5, 6.5};
282279
AxisSpec dataTypeAxis = {9, 0, 9, "Histogram types"};
@@ -510,14 +507,14 @@ struct Chk892pp {
510507
return false;
511508
if (TrackCuts.cfgPrimaryTrack && !track.isPrimaryTrack())
512509
return false;
513-
if (TrackCuts.cfgpTdepDCAxyCut) {
510+
if (TrackCuts.cfgpTdepDCAxyCut) {
514511
// Tuned on the LHC22f anchored MC LHC23d1d on primary pions. 7 Sigmas of the resolution
515512
if (std::abs(track.dcaXY()) > (0.004 + (0.013 / track.pt())))
516-
return false;
517-
} else {
518-
if (std::abs(track.dcaXY()) > TrackCuts.cfgMaxbDCArToPVcut)
519-
return false;
520-
}
513+
return false;
514+
} else {
515+
if (std::abs(track.dcaXY()) > TrackCuts.cfgMaxbDCArToPVcut)
516+
return false;
517+
}
521518
if (std::abs(track.dcaZ()) > TrackCuts.cfgMaxbDCAzToPVcut)
522519
return false;
523520
return true;
@@ -560,10 +557,10 @@ struct Chk892pp {
560557
auto lCPA = candidate.v0cosPA();
561558
auto lPropTauK0s = candidate.distovertotmom(collision.posX(), collision.posY(), collision.posZ()) * MassK0Short;
562559
auto lMk0s = candidate.mK0Short();
563-
auto lMLambda = candidate.mLambda();
560+
auto lMLambda = candidate.mLambda();
564561
auto lMALambda = candidate.mAntiLambda();
565562

566-
auto checkCommonCuts = [&]() {
563+
auto checkCommonCuts = [&]() {
567564
if (lDauDCA > SecondaryCuts.cfgSecondaryDauDCAMax)
568565
return false;
569566
if (lDauPosDCAtoPV < SecondaryCuts.cfgSecondaryDauPosDCAtoPVMin)
@@ -592,7 +589,7 @@ struct Chk892pp {
592589
return true;
593590
};
594591

595-
if (SecondaryCuts.cfgReturnFlag) { // For cut study
592+
if (SecondaryCuts.cfgReturnFlag) { // For cut study
596593
bool returnFlag = true;
597594
histos.fill(HIST("QA/K0sCutCheck"), 0);
598595
if (lDauDCA > SecondaryCuts.cfgSecondaryDauDCAMax) {
@@ -836,7 +833,7 @@ struct Chk892pp {
836833
histos.fill(HIST("QA/after/hCPASecondary"), trkkCPA);
837834
histos.fill(HIST("QA/after/hPropTauSecondary"), trkkPropTau);
838835
histos.fill(HIST("QA/after/hInvmassSecondary"), trkkMass);
839-
histos.fill(HIST("hInvmass_K0s"), lCentrality, lResoSecondary.Pt(), lResoSecondary.M());
836+
histos.fill(HIST("hInvmass_K0s"), lCentrality, lResoSecondary.Pt(), lResoSecondary.M());
840837
}
841838
k0sIndicies.push_back(k0sCand.index());
842839
}
@@ -845,7 +842,7 @@ struct Chk892pp {
845842
for (auto& k0sIndex : k0sIndicies) {
846843
auto bTrack = dTracks1.rawIteratorAt(trackIndex);
847844
auto k0sCand = dTracks2.rawIteratorAt(k0sIndex);
848-
auto trkkMass = k0sCand.mK0Short();
845+
auto trkkMass = k0sCand.mK0Short();
849846

850847
lDecayDaughter_bach.SetXYZM(bTrack.px(), bTrack.py(), bTrack.pz(), MassPionCharged);
851848
lResoSecondary.SetXYZM(k0sCand.px(), k0sCand.py(), k0sCand.pz(), trkkMass);
@@ -867,23 +864,23 @@ struct Chk892pp {
867864
histos.fill(HIST("QA/after/kstarinvmass"), lResoKstar.M());
868865
histos.fill(HIST("hInvmass_Kstar"), typeKstar, lCentrality, lResoKstar.Pt(), lResoKstar.M());
869866

870-
if (BkgEstimationConfig.cfgFillRotBkg) {
871-
for (int i = 0; i < BkgEstimationConfig.cfgNrotBkg; i++) {
872-
auto lRotAngle = BkgEstimationConfig.cfgMinRot + i * ((BkgEstimationConfig.cfgMaxRot - BkgEstimationConfig.cfgMinRot) / (BkgEstimationConfig.cfgNrotBkg - 1));
873-
histos.fill(HIST("QA/RotBkg/hRotBkg"), lRotAngle);
874-
if (BkgEstimationConfig.cfgRotPion) {
875-
lDaughterRot = lDecayDaughter_bach;
876-
lDaughterRot.RotateZ(lRotAngle);
877-
lResonanceRot = lDaughterRot + lResoSecondary;
878-
} else {
879-
lDaughterRot = lResoSecondary;
880-
lDaughterRot.RotateZ(lRotAngle);
881-
lResonanceRot = lDecayDaughter_bach + lDaughterRot;
882-
}
883-
typeKstar = bTrack.sign() > 0 ? BinType::kKstarP_Rot : BinType::kKstarN_Rot;
884-
histos.fill(HIST("hInvmass_Kstar"), typeKstar, lCentrality, lResonanceRot.Pt(), lResonanceRot.M());
885-
}
886-
}
867+
if (BkgEstimationConfig.cfgFillRotBkg) {
868+
for (int i = 0; i < BkgEstimationConfig.cfgNrotBkg; i++) {
869+
auto lRotAngle = BkgEstimationConfig.cfgMinRot + i * ((BkgEstimationConfig.cfgMaxRot - BkgEstimationConfig.cfgMinRot) / (BkgEstimationConfig.cfgNrotBkg - 1));
870+
histos.fill(HIST("QA/RotBkg/hRotBkg"), lRotAngle);
871+
if (BkgEstimationConfig.cfgRotPion) {
872+
lDaughterRot = lDecayDaughter_bach;
873+
lDaughterRot.RotateZ(lRotAngle);
874+
lResonanceRot = lDaughterRot + lResoSecondary;
875+
} else {
876+
lDaughterRot = lResoSecondary;
877+
lDaughterRot.RotateZ(lRotAngle);
878+
lResonanceRot = lDecayDaughter_bach + lDaughterRot;
879+
}
880+
typeKstar = bTrack.sign() > 0 ? BinType::kKstarP_Rot : BinType::kKstarN_Rot;
881+
histos.fill(HIST("hInvmass_Kstar"), typeKstar, lCentrality, lResonanceRot.Pt(), lResonanceRot.M());
882+
}
883+
}
887884
} // IsMix
888885
} // K0scand
889886
} // bTrack

0 commit comments

Comments
 (0)