@@ -285,6 +285,7 @@ struct SGCandProducer {
285285 {
286286 if (verboseInfo)
287287 LOGF (debug, " <SGCandProducer> collision %d" , collision.globalIndex ());
288+
288289 getHist (TH1, histdir + " /Stat" )->Fill (0 ., 1 .);
289290 // reject collisions at TF boundaries
290291 if (rejectAtTFBoundary && !collision.selection_bit (aod::evsel::kNoTimeFrameBorder )) {
@@ -389,6 +390,8 @@ struct SGCandProducer {
389390 fitInfo.BBFV0Apf , fitInfo.BGFV0Apf ,
390391 fitInfo.BBFDDApf , fitInfo.BBFDDCpf , fitInfo.BGFDDApf , fitInfo.BGFDDCpf );
391392 outputCollisionSelExtras (chFT0A, chFT0C, chFDDA, chFDDC, chFV0A, occ, ir, trs, trofs, hmpr, tfb, itsROFb, sbp, zVtxFT0vPv, vtxITSTPC, collision.rct_raw ());
393+ if (verboseInfo)
394+ LOGF (info, " %s Coll GID %d" , histdir, collision.globalIndex ());
392395 outputCollsLabels (collision.globalIndex ());
393396 if (newbc.has_zdc ()) {
394397 auto zdc = newbc.zdc ();
@@ -481,7 +484,7 @@ struct SGCandProducer {
481484
482485 if (std::find (generatorIds->begin (), generatorIds->end (), mccol.getGeneratorId ()) != generatorIds->end ()) {
483486 if (verboseInfo)
484- LOGF (info, " Event with good generatorId" );
487+ LOGF (info, " Event with good generatorId %d " , mccol. getGeneratorId () );
485488 processReco (std::string (" MCreco" ), collision, bcs, tracks, fwdtracks, fv0as, ft0s, fdds);
486489 }
487490 }
@@ -688,6 +691,8 @@ struct McSGCandProducer {
688691 // loop over McCollisions and UDCCs simultaneously
689692 auto mccol = mccols.iteratorAt (0 );
690693 auto mcOfInterest = std::find (generatorIds->begin (), generatorIds->end (), mccol.getGeneratorId ()) != generatorIds->end ();
694+ if (verboseInfoMC)
695+ LOGF (info, " Is Generator ID OK %d, MCcoll GenId %d, Set in json ID %d" , mcOfInterest, mccol.getGeneratorId (), *(generatorIds->begin ()));
691696 auto lastmccol = mccols.iteratorAt (mccols.size () - 1 );
692697 auto mccolAtEnd = false ;
693698
@@ -701,6 +706,7 @@ struct McSGCandProducer {
701706 bool goon = true ;
702707 while (goon) {
703708 auto globBC = mccol.bc_as <BCs>().globalBC ();
709+
704710 // check if dgcand has an associated McCollision
705711 if (sgcand.has_collision ()) {
706712 auto sgcandCol = sgcand.collision_as <CCs>();
@@ -811,7 +817,7 @@ struct McSGCandProducer {
811817 // advance mccol
812818 if (mccol != lastmccol) {
813819 mccol++;
814- mcOfInterest = std::find (generatorIds->begin (), generatorIds->end (), mccol.getGeneratorId ()) != generatorIds->end ();
820+ // mcOfInterest = std::find(generatorIds->begin(), generatorIds->end(), mccol.getGeneratorId()) != generatorIds->end();
815821 mccolId = mccol.globalIndex ();
816822 } else {
817823 mccolAtEnd = true ;
@@ -875,6 +881,7 @@ struct McSGCandProducer {
875881 LOGF (info, " Number of McCollisions %d" , mccols.size ());
876882 LOGF (info, " Number of SG candidates %d" , sgcands.size ());
877883 LOGF (info, " Number of UD tracks %d" , udtracks.size ());
884+ LOGF (info, " Number of McParticles %d" , mcparts.size ());
878885 }
879886 if (mccols.size () > 0 ) {
880887 if (sgcands.size () > 0 ) {
0 commit comments