@@ -185,6 +185,7 @@ struct Chargedkstaranalysis {
185185 Configurable<bool > confincludeCentralityMC{" confincludeCentralityMC" , false , " Include centrality in MC" };
186186 Configurable<bool > confEvtCollInTimeRangeStandard{" confEvtCollInTimeRangeStandard" , true , " Evt sel: apply NoCollInTimeRangeStandard" };
187187 Configurable<bool > confEvtCollInTimeRangeNarrow{" confEvtCollInTimeRangeNarrow" , false , " Evt sel: apply NoCollInTimeRangeNarrow" };
188+
188189 Configurable<bool > confEvtCollInTimeRangeStrict{" confEvtCollInTimeRangeStrict" , false , " Evt sel: apply NoCollInTimeRangeStrict" };
189190 Configurable<bool > confEvtNoCollInRofStandard{" confEvtNoCollInRofStandard" , false , " Evt sel: apply NoCollInRofStandard" };
190191 Configurable<bool > confEvtNoCollInRofStrict{" confEvtNoCollInRofStrict" , false , " Evt sel: apply NoCollInRofStrict" };
@@ -295,9 +296,7 @@ struct Chargedkstaranalysis {
295296 {
296297 centrality = -999 ;
297298
298- rctCut.rctChecker .init (rctCut.cfgEvtRCTFlagCheckerLabel , rctCut.cfgEvtRCTFlagCheckerZDCCheck , rctCut.cfgEvtRCTFlagCheckerLimitAcceptAsBad );
299-
300- colCuts.setCuts (eventCutCfgs.confEvtZvtx , eventCutCfgs.confEvtTriggerCheck , eventCutCfgs.confEvtOfflineCheck , /* checkRun3*/ true , /* triggerTVXsel*/ false , eventCutCfgs.confEvtOccupancyInTimeRangeMax , eventCutCfgs.confEvtOccupancyInTimeRangeMin );
299+ colCuts.setCuts (eventCutCfgs.confEvtZvtx , eventCutCfgs.confEvtTriggerCheck , eventCutCfgs.confEvtOfflineCheck , /* checkRun3*/ true , eventCutCfgs.confEvtTriggerTVXSel , eventCutCfgs.confEvtOccupancyInTimeRangeMax , eventCutCfgs.confEvtOccupancyInTimeRangeMin );
301300 colCuts.init (&histos);
302301 colCuts.setTriggerTVX (eventCutCfgs.confEvtTriggerTVXSel );
303302 colCuts.setApplyTFBorderCut (eventCutCfgs.confEvtTFBorderCut );
@@ -318,6 +317,8 @@ struct Chargedkstaranalysis {
318317 colCuts.setApplyBBT0A (eventCutCfgs.confEvtBBT0A );
319318 colCuts.setApplyBBT0C (eventCutCfgs.confEvtBBT0C );
320319
320+ rctCut.rctChecker .init (rctCut.cfgEvtRCTFlagCheckerLabel , rctCut.cfgEvtRCTFlagCheckerZDCCheck , rctCut.cfgEvtRCTFlagCheckerLimitAcceptAsBad );
321+
321322 AxisSpec centAxis = {axisCfgs.cfgBinsCent , " T0M (%)" };
322323 AxisSpec vtxzAxis = {axisCfgs.cfgBinsVtxZ , " Z Vertex (cm)" };
323324 AxisSpec ptAxis = {axisCfgs.cfgBinsPt , " #it{p}_{T} (GeV/#it{c})" };
@@ -329,7 +330,7 @@ struct Chargedkstaranalysis {
329330 AxisSpec dcaxyAxis = {200 , 0 , 2 , " DCA_{#it{xy}} (cm)" };
330331 AxisSpec dcazAxis = {200 , 0 , 2 , " DCA_{#it{z}} (cm)" };
331332 AxisSpec yAxis = {100 , -1 , 1 , " Rapidity" };
332- AxisSpec invMassAxisK0s = {400 / cNbinsDiv, 0.3 , 0.7 , " Invariant Mass (GeV/#it{c}^2)" }; // K0s ~497.611
333+ AxisSpec invMassAxisK0s = {800 / cNbinsDiv, 0.46 , 0.54 , " Invariant Mass (GeV/#it{c}^2)" }; // K0s ~497.611
333334 AxisSpec invMassAxisReso = {900 / cNbinsDiv, 0 .5f , 1 .4f , " Invariant Mass (GeV/#it{c}^2)" }; // chK(892) ~892
334335 AxisSpec invMassAxisScan = {150 , 0 , 1.5 , " Invariant Mass (GeV/#it{c}^2)" }; // For selection
335336 AxisSpec pidQAAxis = {130 , -6.5 , 6.5 };
@@ -702,13 +703,13 @@ struct Chargedkstaranalysis {
702703 bool selectionK0s (CollisionType const & collision, K0sType const & candidate)
703704 {
704705 int ibin = 1 ;
705- auto dauDCA = candidate.dcaV0daughters ();
706- auto dauPosDCAtoPV = candidate.dcapostopv ();
707- auto dauNegDCAtoPV = candidate.dcanegtopv ();
706+ auto dauDCA = std::fabs ( candidate.dcaV0daughters () );
707+ auto dauPosDCAtoPV = std::fabs ( candidate.dcapostopv () );
708+ auto dauNegDCAtoPV = std::fabs ( candidate.dcanegtopv () );
708709 auto pT = candidate.pt ();
709710 auto rapidity = candidate.yK0Short ();
710711 auto v0Radius = candidate.v0radius ();
711- auto dcaToPV = candidate.dcav0topv ();
712+ auto dcaToPV = std::fabs ( candidate.dcav0topv () );
712713 auto cosPA = candidate.v0cosPA ();
713714 auto propTauK0s = candidate.distovertotmom (collision.posX (), collision.posY (), collision.posZ ()) * massK0s;
714715 auto mK0s = candidate.mK0Short ();
@@ -717,20 +718,20 @@ struct Chargedkstaranalysis {
717718
718719 bool returnFlag = true ;
719720 histos.fill (HIST (" QA/K0sCutCheck" ), ibin);
720- if (dauDCA > secondaryCutsCfgs.cSecondaryDauDCAMax ) {
721+ if (std::fabs ( dauDCA) > secondaryCutsCfgs.cSecondaryDauDCAMax ) {
721722 returnFlag = false ;
722723 }
723724 ibin++;
724725 if (returnFlag == true )
725726 histos.fill (HIST (" QA/K0sCutCheck" ), ibin);
726- if (dauPosDCAtoPV < secondaryCutsCfgs.cSecondaryDauPosDCAtoPVMin ) {
727+ if (std::fabs ( dauPosDCAtoPV) < secondaryCutsCfgs.cSecondaryDauPosDCAtoPVMin ) {
727728 returnFlag = false ;
728729 }
729730 ibin++;
730731 if (returnFlag == true )
731732 histos.fill (HIST (" QA/K0sCutCheck" ), ibin);
732733
733- if (dauNegDCAtoPV < secondaryCutsCfgs.cSecondaryDauNegDCAtoPVMin ) {
734+ if (std::fabs ( dauNegDCAtoPV) < secondaryCutsCfgs.cSecondaryDauNegDCAtoPVMin ) {
734735 returnFlag = false ;
735736 }
736737 ibin++;
0 commit comments