@@ -649,14 +649,14 @@ struct TaskPolarisationCharmHadrons {
649649 hEPaxes.insert (hEPaxes.end (), {thnAxisAbsEtaTrackMin, thnAxisNumItsClsMin, thnAxisNumTpcClsMin});
650650 }
651651 if (doprocessDstarMc || doprocessDstarMcWithMl || doprocessDstarMcInPbPb || doprocessDstarMcWithMlInPbPb) {
652- std::vector<AxisSpec> RecoPromptEPAxes (hEPaxes);
653- RecoPromptEPAxes .insert (RecoPromptEPAxes .end (), {thnAxisDauToMuons});
652+ std::vector<AxisSpec> hRecoPromptEPAxes (hEPaxes);
653+ hRecoPromptEPAxes .insert (hRecoPromptEPAxes .end (), {thnAxisDauToMuons});
654654 std::vector<AxisSpec> hRecoNonPromptEPAxes (hEPaxes);
655655 hRecoNonPromptEPAxes.insert (hRecoNonPromptEPAxes.end (), {thnAxisDauToMuons, thnAxisPtB});
656- registry.add (" hRecoPromptEP" , " THn for polarisation studies with cosThStar w.r.t. event plane axis and BDT scores for reconstructed prompt D*+ candidates" , HistType::kTHnSparseF , RecoPromptEPAxes );
656+ registry.add (" hRecoPromptEP" , " THn for polarisation studies with cosThStar w.r.t. event plane axis and BDT scores for reconstructed prompt D*+ candidates" , HistType::kTHnSparseF , hRecoPromptEPAxes );
657657 registry.add (" hRecoNonPromptEP" , " THn for polarisation studies with cosThStar w.r.t. event plane axis and BDT scores for reconstructed non-prompt D*+ candidates" , HistType::kTHnSparseF , hRecoNonPromptEPAxes);
658658 if (activatePartRecoDstar) {
659- registry.add (" hPartRecoPromptEP" , " THn for polarisation studies with cosThStar w.r.t. event plane axis and BDT scores for partially reconstructed prompt D*+ candidates" , HistType::kTHnSparseF , RecoPromptEPAxes );
659+ registry.add (" hPartRecoPromptEP" , " THn for polarisation studies with cosThStar w.r.t. event plane axis and BDT scores for partially reconstructed prompt D*+ candidates" , HistType::kTHnSparseF , hRecoPromptEPAxes );
660660 registry.add (" hPartRecoNonPromptEP" , " THn for polarisation studies with cosThStar w.r.t. event plane axis and BDT scores for partially reconstructed non-prompt D*+ candidates" , HistType::kTHnSparseF , hRecoNonPromptEPAxes);
661661 }
662662 } else {
@@ -1406,12 +1406,18 @@ struct TaskPolarisationCharmHadrons {
14061406 template <charm_polarisation::DecayChannel channel>
14071407 bool isInSignalRegion (float invMass)
14081408 {
1409+ float minInvMass = 0 .f ;
1410+ float maxInvMass = 100 .f ;
14091411 if constexpr (channel == charm_polarisation::DecayChannel::DstarToDzeroPi) { // D*+
1410- if (0 .142f < invMass && invMass < 0 .15f ) {
1412+ minInvMass = 0 .142f ;
1413+ maxInvMass = 0 .15f ;
1414+ if (minInvMass < invMass && invMass < maxInvMass) {
14111415 return true ;
14121416 }
14131417 } else if constexpr (channel == charm_polarisation::DecayChannel::LcToPKPi) { // Lc+ (to be tuned!)
1414- if (2 .25f < invMass && invMass < 2 .35f ) {
1418+ minInvMass = 2 .25f ;
1419+ maxInvMass = 2 .35f ;
1420+ if (minInvMass < invMass && invMass < maxInvMass) {
14151421 return true ;
14161422 }
14171423 }
0 commit comments