@@ -419,6 +419,7 @@ struct FullJetSpectra {
419419 }
420420 if (doprocessJetsMCP || doprocessJetsMCPWeighted) {
421421 registry.add (" hPartcollisionCounter" , " event status;event status;entries" , {HistType::kTH1F , {{10 , 0.0 , 10.0 }}});
422+ registry.add (" hRecoMatchesPerMcCollision" , " split vertices QA;;entries" , {HistType::kTH1F , {{5 , 0.0 , 5.0 }}});
422423
423424 registry.add (" h_full_mcpjet_tablesize" , " " , {HistType::kTH1F , {{4 , 0 ., 5 .}}});
424425 registry.add (" h_full_mcpjet_ntracks" , " " , {HistType::kTH1F , {{200 , -0.5 , 200 .}}});
@@ -524,7 +525,7 @@ struct FullJetSpectra {
524525
525526 if (doprocessMBMCPCollisionsWithMultiplicity || doprocessMBMCPCollisionsWeightedWithMultiplicity) {
526527 registry.add (" hPartEventmultiplicityCounter" , " event status;event status;entries" , {HistType::kTH1F , {{11 , 0.0 , 11.0 }}});
527- registry.add (" hRecoMatchesPerMcCollision " , " split vertices QA;;entries" , {HistType::kTH1F , {{5 , 0.0 , 5.0 }}});
528+ registry.add (" hRecoMatchesPerMcCollisionMult " , " split vertices QA;;entries" , {HistType::kTH1F , {{5 , 0.0 , 5.0 }}});
528529 registry.add (" hMCCollMatchedFT0Mult" , " " , {HistType::kTH1F , {{3500 , 0 ., 3500 .}}});
529530 registry.add (" hMCCollMatchedFT0Cent" , " " , {HistType::kTH1F , {{105 , 0 ., 105 .}}});
530531
@@ -900,19 +901,19 @@ struct FullJetSpectra {
900901 registry.fill (HIST (" hDetcollisionCounter" ), 1.5 ); // DetCollWithVertexZ
901902
902903 if (!jetderiveddatautilities::selectCollision (collision, eventSelectionBits, doMBGapTrigger)) {
903- registry.fill (HIST (" hDetcollisionCounter" ), 2 .5 ); // EventsNotSatisfyingEventSelection
904+ registry.fill (HIST (" hDetcollisionCounter" ), 4 .5 ); // EventsNotSatisfyingEventSelection
904905 return ;
905906 }
906907 if (doEMCALEventWorkaround) {
907908 if (collision.isEmcalReadout () && !collision.isAmbiguous ()) { // i.e. EMCAL has a cell content
908909 eventAccepted = true ;
909910 if (collision.alias_bit (kTVXinEMC )) {
910- registry.fill (HIST (" hDetcollisionCounter" ), 3 .5 ); // EMCreadoutDetEventsWithkTVXinEMC
911+ registry.fill (HIST (" hDetcollisionCounter" ), 5 .5 ); // EMCreadoutDetEventsWithkTVXinEMC
911912 }
912913 }
913914 } else {
914915 if (!collision.isAmbiguous () && jetderiveddatautilities::eventEMCAL (collision) && collision.alias_bit (kTVXinEMC )) {
915- registry.fill (HIST (" hDetcollisionCounter" ), 3 .5 ); // EMCreadoutDetEventsWithkTVXinEMC
916+ registry.fill (HIST (" hDetcollisionCounter" ), 5 .5 ); // EMCreadoutDetEventsWithkTVXinEMC
916917 eventAccepted = true ;
917918 }
918919 }
@@ -923,10 +924,10 @@ struct FullJetSpectra {
923924 fillRejectedJetHistograms (jet, 1.0 );
924925 }
925926 }
926- registry.fill (HIST (" hDetcollisionCounter" ), 4 .5 ); // AllRejectedDetEventsAfterEMCEventSelection
927+ registry.fill (HIST (" hDetcollisionCounter" ), 6 .5 ); // AllRejectedDetEventsAfterEMCEventSelection
927928 return ;
928929 }
929- registry.fill (HIST (" hDetcollisionCounter" ), 5 .5 ); // EMCAcceptedDetColl
930+ registry.fill (HIST (" hDetcollisionCounter" ), 7 .5 ); // EMCAcceptedDetColl
930931
931932 for (auto const & jet : jets) {
932933 if (!jetfindingutilities::isInEtaAcceptance (jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
@@ -2390,7 +2391,7 @@ struct FullJetSpectra {
23902391 // Perform MC Collision matching, i.e. match the current MC collision to its associated reco (MCD) collision
23912392 // to get the corresponding FT0M component at the particle level
23922393 auto collisionspermcpjet = collisions.sliceBy (CollisionsPerMCPCollision, mccollision.globalIndex ());
2393- registry.fill (HIST (" hRecoMatchesPerMcCollision " ), collisionspermcpjet.size ()); // for split vertices QA
2394+ registry.fill (HIST (" hRecoMatchesPerMcCollisionMult " ), collisionspermcpjet.size ()); // for split vertices QA
23942395
23952396 if (collisionspermcpjet.size () == 0 || collisionspermcpjet.size () < 1 ) {
23962397 registry.fill (HIST (" hPartEventmultiplicityCounter" ), 4.5 ); // RejectedPartCollForDetCollWithSize0or<1
@@ -2555,7 +2556,7 @@ struct FullJetSpectra {
25552556 // Perform MC Collision matching, i.e. match the current MC collision to its associated reco (MCD) collision
25562557 // to get the corresponding FT0M component at the particle level
25572558 auto collisionspermcpjet = collisions.sliceBy (CollisionsPerMCPCollision, mccollision.globalIndex ());
2558- registry.fill (HIST (" hRecoMatchesPerMcCollision " ), collisionspermcpjet.size (), mccollision.weight ()); // for split vertices QA
2559+ registry.fill (HIST (" hRecoMatchesPerMcCollisionMult " ), collisionspermcpjet.size (), mccollision.weight ()); // for split vertices QA
25592560
25602561 if (collisionspermcpjet.size () == 0 || collisionspermcpjet.size () < 1 ) {
25612562 registry.fill (HIST (" hPartEventmultiplicityCounter" ), 4.5 , mccollision.weight ()); // RejectedWeightedPartCollForDetCollWithSize0or<1
0 commit comments