@@ -680,19 +680,31 @@ class O2HitMerger : public FairMQDevice
680680 int eta1Point2Counter = 0 ;
681681 int eta1Point0Counter = 0 ;
682682 int eta0Point8Counter = 0 ;
683+ int eta1Point2CounterPi = 0 ;
684+ int eta1Point0CounterPi = 0 ;
685+ int eta0Point8CounterPi = 0 ;
683686 int prims = 0 ;
684687 for (auto & tr : tracks) {
685688 if (tr.isPrimary ()) {
686689 prims++;
687690 const auto eta = tr.GetEta ();
688691 if (eta < 1.2 ) {
689692 eta1Point2Counter++;
693+ if (std::abs (tr.GetPdgCode ()) == 211 ) {
694+ eta1Point2CounterPi++;
695+ }
690696 }
691697 if (eta < 1.0 ) {
692698 eta1Point0Counter++;
699+ if (std::abs (tr.GetPdgCode ()) == 211 ) {
700+ eta1Point0CounterPi++;
701+ }
693702 }
694703 if (eta < 0.8 ) {
695704 eta0Point8Counter++;
705+ if (std::abs (tr.GetPdgCode ()) == 211 ) {
706+ eta0Point8CounterPi++;
707+ }
696708 }
697709 } else {
698710 break ; // track layout is such that all prims are first anyway
@@ -703,6 +715,9 @@ class O2HitMerger : public FairMQDevice
703715 eventheader->putInfo (" prims_eta_1.2" , eta1Point2Counter);
704716 eventheader->putInfo (" prims_eta_1.0" , eta1Point0Counter);
705717 eventheader->putInfo (" prims_eta_0.8" , eta0Point8Counter);
718+ eventheader->putInfo (" prims_eta_1.2_pi" , eta1Point2CounterPi);
719+ eventheader->putInfo (" prims_eta_1.0_pi" , eta1Point0CounterPi);
720+ eventheader->putInfo (" prims_eta_0.8_pi" , eta0Point8CounterPi);
706721 eventheader->putInfo (" prims_total" , prims);
707722 };
708723
0 commit comments