@@ -683,6 +683,7 @@ struct StrangeCascTrack {
683683 }
684684
685685 // statistics - compare gen and reco pt and rapidity
686+ int pdgCharmBeautyCodes = 4000 ;
686687 if constexpr (requires { collision.straMCCollisionId (); }) {
687688 if constexpr (requires { stdCasc.has_cascMCCore (); }) {
688689 auto cascmccore = stdCasc.template cascMCCore_as <DerMCGenCascades>();
@@ -703,13 +704,13 @@ struct StrangeCascTrack {
703704 histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/GenRecRapidityXi" ), genYXi, cascade.yXi ());
704705 histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAxyPrimaryXi" ), cascade.dcaXYCascToPV (), pt, mult);
705706 histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAzPrimaryXi" ), cascade.dcaZCascToPV (), pt, mult);
706- if (std::abs (cascmccore.pdgCodeMother ()) > 4000 ) {
707- histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAxyDecayXi" ), cascade.dcaXYCascToPV (), pt, mult);
708- histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAzDecayXi" ), cascade.dcaZCascToPV (), pt, mult);
709- } else {
710- histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAxyDirectXi" ), cascade.dcaXYCascToPV (), pt, mult);
711- histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAzDirectXi" ), cascade.dcaZCascToPV (), pt, mult);
712- }
707+ if (std::abs (cascmccore.pdgCodeMother ()) > pdgCharmBeautyCodes ) {
708+ histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAxyDecayXi" ), cascade.dcaXYCascToPV (), pt, mult);
709+ histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAzDecayXi" ), cascade.dcaZCascToPV (), pt, mult);
710+ } else {
711+ histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAxyDirectXi" ), cascade.dcaXYCascToPV (), pt, mult);
712+ histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAzDirectXi" ), cascade.dcaZCascToPV (), pt, mult);
713+ }
713714 }
714715 if (passedAllSelsXi) {
715716 histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec/GenRecRapidityXi" ), genYXi, cascade.yXi ());
@@ -718,26 +719,26 @@ struct StrangeCascTrack {
718719 histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAxyPrimaryXi" ), cascade.dcaXYCascToPV (), pt, mult);
719720 histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAzPrimaryXi" ), cascade.dcaZCascToPV (), pt, mult);
720721 histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/Radius/PrimaryXi" ), cascade.cascradius (), pt);
721- if (std::abs (cascmccore.pdgCodeMother ()) > 4000 ) {
722- histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAxyDecayXi" ), cascade.dcaXYCascToPV (), pt, mult);
723- histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAzDecayXi" ), cascade.dcaZCascToPV (), pt, mult);
722+ if (std::abs (cascmccore.pdgCodeMother ()) > pdgCharmBeautyCodes ) {
723+ histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAxyDecayXi" ), cascade.dcaXYCascToPV (), pt, mult);
724+ histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAzDecayXi" ), cascade.dcaZCascToPV (), pt, mult);
724725 } else {
725- histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAxyDirectXi" ), cascade.dcaXYCascToPV (), pt, mult);
726- histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAzDirectXi" ), cascade.dcaZCascToPV (), pt, mult);
726+ histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAxyDirectXi" ), cascade.dcaXYCascToPV (), pt, mult);
727+ histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAzDirectXi" ), cascade.dcaZCascToPV (), pt, mult);
727728 }
728729 }
729730 }
730731 if (fillTruthOmega) {
731732 histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/GenRecRapidityOmega" ), genYOmega, cascade.yOmega ());
732733 histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAxyPrimaryOmega" ), cascade.dcaXYCascToPV (), pt, mult);
733734 histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAzPrimaryOmega" ), cascade.dcaZCascToPV (), pt, mult);
734- if (std::abs (cascmccore.pdgCodeMother ()) > 4000 ) {
735- histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAxyDecayOmega" ), cascade.dcaXYCascToPV (), pt, mult);
736- histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAzDecayOmega" ), cascade.dcaZCascToPV (), pt, mult);
737- } else {
738- histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAxyDirectOmega" ), cascade.dcaXYCascToPV (), pt, mult);
739- histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAzDirectOmega" ), cascade.dcaZCascToPV (), pt, mult);
740- }
735+ if (std::abs (cascmccore.pdgCodeMother ()) > pdgCharmBeautyCodes ) {
736+ histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAxyDecayOmega" ), cascade.dcaXYCascToPV (), pt, mult);
737+ histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAzDecayOmega" ), cascade.dcaZCascToPV (), pt, mult);
738+ } else {
739+ histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAxyDirectOmega" ), cascade.dcaXYCascToPV (), pt, mult);
740+ histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAzDirectOmega" ), cascade.dcaZCascToPV (), pt, mult);
741+ }
741742 }
742743 if (passedAllSelsOmega) {
743744 histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec/GenRecRapidityOmega" ), genYOmega, cascade.yOmega ());
@@ -746,13 +747,13 @@ struct StrangeCascTrack {
746747 histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAxyPrimaryOmega" ), cascade.dcaXYCascToPV (), pt, mult);
747748 histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAzPrimaryOmega" ), cascade.dcaZCascToPV (), pt, mult);
748749 histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/Radius/PrimaryOmega" ), cascade.cascradius (), pt);
749- if (std::abs (cascmccore.pdgCodeMother ()) > 4000 ) {
750- histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAxyDecayOmega" ), cascade.dcaXYCascToPV (), pt, mult);
751- histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAzDecayOmega" ), cascade.dcaZCascToPV (), pt, mult);
752- } else {
753- histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAxyDirectOmega" ), cascade.dcaXYCascToPV (), pt, mult);
754- histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAzDirectOmega" ), cascade.dcaZCascToPV (), pt, mult);
755- }
750+ if (std::abs (cascmccore.pdgCodeMother ()) > pdgCharmBeautyCodes ) {
751+ histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAxyDecayOmega" ), cascade.dcaXYCascToPV (), pt, mult);
752+ histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAzDecayOmega" ), cascade.dcaZCascToPV (), pt, mult);
753+ } else {
754+ histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAxyDirectOmega" ), cascade.dcaXYCascToPV (), pt, mult);
755+ histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAzDirectOmega" ), cascade.dcaZCascToPV (), pt, mult);
756+ }
756757 }
757758 }
758759 }
@@ -1193,11 +1194,11 @@ struct StrangeCascTrack {
11931194 void processDerivedMCRec (DerMCRecCollisions const & collisions, DerMCRecCascDatas const & allCascs, DerMCRecTraCascDatas const & traCascs, DauTracks const &, DerMCGenCascades const &)
11941195 {
11951196 // custom group tracked cascades - a temporary fix for wrong ordering in OO MC
1196- traCascsGrouped.clear ();
1197- traCascsGrouped.resize (collisions.size ());
1198- for (const auto & casc : traCascs) {
1199- traCascsGrouped[casc.straCollisionId ()].push_back (casc.globalIndex ());
1200- }
1197+ traCascsGrouped.clear ();
1198+ traCascsGrouped.resize (collisions.size ());
1199+ for (const auto & casc : traCascs) {
1200+ traCascsGrouped[casc.straCollisionId ()].push_back (casc.globalIndex ());
1201+ }
12011202 for (const auto & collision : collisions) {
12021203 fillEvents (collision); // save info about all processed events
12031204 auto slicedAllCascs = allCascs.sliceBy (cascsPerCollision, collision.globalIndex ());
@@ -1209,8 +1210,8 @@ struct StrangeCascTrack {
12091210 histos.fill (HIST (" Rec-Events/Mult" ), mult);
12101211 analyseCascs (collision, slicedAllCascs); // process all cascades
12111212 for (int const & idx : traCascsGrouped[collision.globalIndex ()]) {
1212- auto casc = traCascs.rawIteratorAt (idx);
1213- analyseCascade (collision, casc);
1213+ auto casc = traCascs.rawIteratorAt (idx);
1214+ analyseCascade (collision, casc);
12141215 }
12151216 }
12161217 }
0 commit comments