@@ -705,12 +705,12 @@ struct StrangeCascTrack {
705705 histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAxyPrimaryXi" ), cascade.dcaXYCascToPV (), pt, mult);
706706 histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAzPrimaryXi" ), cascade.dcaZCascToPV (), pt, mult);
707707 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- }
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+ }
714714 }
715715 if (passedAllSelsXi) {
716716 histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec/GenRecRapidityXi" ), genYXi, cascade.yXi ());
@@ -720,11 +720,11 @@ struct StrangeCascTrack {
720720 histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAzPrimaryXi" ), cascade.dcaZCascToPV (), pt, mult);
721721 histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/Radius/PrimaryXi" ), cascade.cascradius (), pt);
722722 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);
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);
725725 } else {
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);
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);
728728 }
729729 }
730730 }
@@ -733,12 +733,12 @@ struct StrangeCascTrack {
733733 histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAxyPrimaryOmega" ), cascade.dcaXYCascToPV (), pt, mult);
734734 histos.fill (HIST (TypeNames[Type]) + HIST (" /NoSel-Truth/DCA/DCAzPrimaryOmega" ), cascade.dcaZCascToPV (), pt, mult);
735735 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- }
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+ }
742742 }
743743 if (passedAllSelsOmega) {
744744 histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec/GenRecRapidityOmega" ), genYOmega, cascade.yOmega ());
@@ -748,12 +748,12 @@ struct StrangeCascTrack {
748748 histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/DCA/DCAzPrimaryOmega" ), cascade.dcaZCascToPV (), pt, mult);
749749 histos.fill (HIST (TypeNames[Type]) + HIST (" /Rec-Truth/Radius/PrimaryOmega" ), cascade.cascradius (), pt);
750750 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- }
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+ }
757757 }
758758 }
759759 }
@@ -1194,11 +1194,11 @@ struct StrangeCascTrack {
11941194 void processDerivedMCRec (DerMCRecCollisions const & collisions, DerMCRecCascDatas const & allCascs, DerMCRecTraCascDatas const & traCascs, DauTracks const &, DerMCGenCascades const &)
11951195 {
11961196 // custom group tracked cascades - a temporary fix for wrong ordering in OO MC
1197- traCascsGrouped.clear ();
1198- traCascsGrouped.resize (collisions.size ());
1199- for (const auto & casc : traCascs) {
1200- traCascsGrouped[casc.straCollisionId ()].push_back (casc.globalIndex ());
1201- }
1197+ traCascsGrouped.clear ();
1198+ traCascsGrouped.resize (collisions.size ());
1199+ for (const auto & casc : traCascs) {
1200+ traCascsGrouped[casc.straCollisionId ()].push_back (casc.globalIndex ());
1201+ }
12021202 for (const auto & collision : collisions) {
12031203 fillEvents (collision); // save info about all processed events
12041204 auto slicedAllCascs = allCascs.sliceBy (cascsPerCollision, collision.globalIndex ());
@@ -1210,8 +1210,8 @@ struct StrangeCascTrack {
12101210 histos.fill (HIST (" Rec-Events/Mult" ), mult);
12111211 analyseCascs (collision, slicedAllCascs); // process all cascades
12121212 for (int const & idx : traCascsGrouped[collision.globalIndex ()]) {
1213- auto casc = traCascs.rawIteratorAt (idx);
1214- analyseCascade (collision, casc);
1213+ auto casc = traCascs.rawIteratorAt (idx);
1214+ analyseCascade (collision, casc);
12151215 }
12161216 }
12171217 }
0 commit comments