Skip to content

Commit 0d30ae4

Browse files
altsybeealtsybee
authored
adding a check for FT0 nearby activity
1 parent 4fb62e4 commit 0d30ae4

File tree

1 file changed

+116
-41
lines changed

1 file changed

+116
-41
lines changed

‎DPG/Tasks/AOTEvent/lightIonsEvSelQa.cxx‎

Lines changed: 116 additions & 41 deletions
Original file line numberDiff line numberDiff line change
@@ -19,11 +19,14 @@
1919
#include "Common/Core/TrackSelectionDefaults.h"
2020
// #include "Common/DataModel/Centrality.h"
2121
#include "Common/DataModel/EventSelection.h"
22+
#include "Common/DataModel/FT0Corrected.h"
2223
#include "Common/DataModel/Multiplicity.h"
2324
#include "Common/DataModel/TrackSelectionTables.h"
2425

2526
#include "CCDB/BasicCCDBManager.h"
2627
#include "CommonDataFormat/BunchFilling.h"
28+
#include "DataFormatsFT0/Digit.h"
29+
#include "DataFormatsFT0/RecPoints.h"
2730
#include "DataFormatsParameters/AggregatedRunInfo.h"
2831
#include "DataFormatsParameters/GRPECSObject.h"
2932
#include "DataFormatsParameters/GRPLHCIFData.h"
@@ -325,23 +328,41 @@ struct LightIonsEvSelQa {
325328
histos.add("noPastActivity/hColTimeResVsNcontrib", "", kTH2F, {axisNcontrib, axisColTimeRes});
326329
histos.add("noPastActivity/hColBcDiffVsNcontrib", "", kTH2F, {axisNcontrib, axisBcDiff});
327330

328-
histos.add("noPileup_cutByVzDiff_pvTOF/hBcColNoSel8", "", kTH1F, {axisBCs});
329-
histos.add("noPileup_cutByVzDiff_pvTOF/hBcTVX", "", kTH1F, {axisBCs});
330-
histos.add("noPileup_cutByVzDiff_pvTOF/hBcFT0", "", kTH1F, {axisBCs});
331-
histos.add("noPileup_cutByVzDiff_pvTOF/hBcFV0", "", kTH1F, {axisBCs});
332-
histos.add("noPileup_cutByVzDiff_pvTOF/hBcFDD", "", kTH1F, {axisBCs});
333-
histos.add("noPileup_cutByVzDiff_pvTOF/hBcZDC", "", kTH1F, {axisBCs});
334-
histos.add("noPileup_cutByVzDiff_pvTOF/hVtxFT0VsVtxCol", "", kTH2F, {axisVtxZ, axisVtxZ});
335-
histos.add("noPileup_cutByVzDiff_pvTOF/hVtxFT0MinusVtxColVsMultT0M", "", kTH2F, {axisVtxZ, axisMultT0M});
336-
histos.add("noPileup_cutByVzDiff_pvTOF/nTracksPV_vs_V0A", "", kTH2F, {axisMultV0A, axisNtracks});
337-
histos.add("noPileup_cutByVzDiff_pvTOF/nTracksPV_vs_T0A", "", kTH2F, {axisMultFT0A, axisNtracks});
338-
histos.add("noPileup_cutByVzDiff_pvTOF/nTracksPV_vs_T0C", "", kTH2F, {axisMultFT0C, axisNtracks});
339-
histos.add("noPileup_cutByVzDiff_pvTOF/nTracksGlobal_vs_V0A", "", kTH2F, {axisMultV0A, axisNtracksGlobal});
340-
histos.add("noPileup_cutByVzDiff_pvTOF/nTracksGlobal_vs_T0A", "", kTH2F, {axisMultFT0A, axisNtracksGlobal});
341-
histos.add("noPileup_cutByVzDiff_pvTOF/nTracksGlobal_vs_T0C", "", kTH2F, {axisMultFT0C, axisNtracksGlobal});
342-
histos.add("noPileup_cutByVzDiff_pvTOF/hTVXvsBcDiff", "", kTH1F, {axisBcDiff});
343-
histos.add("noPileup_cutByVzDiff_pvTOF/hColTimeResVsNcontrib", "", kTH2F, {axisNcontrib, axisColTimeRes});
344-
histos.add("noPileup_cutByVzDiff_pvTOF/hColBcDiffVsNcontrib", "", kTH2F, {axisNcontrib, axisBcDiff});
331+
histos.add("noFT0activityNearby/hBcColNoSel8", "", kTH1F, {axisBCs});
332+
histos.add("noFT0activityNearby/hBcTVX", "", kTH1F, {axisBCs});
333+
histos.add("noFT0activityNearby/hBcFT0", "", kTH1F, {axisBCs});
334+
histos.add("noFT0activityNearby/hBcFV0", "", kTH1F, {axisBCs});
335+
histos.add("noFT0activityNearby/hBcFDD", "", kTH1F, {axisBCs});
336+
histos.add("noFT0activityNearby/hBcZDC", "", kTH1F, {axisBCs});
337+
histos.add("noFT0activityNearby/hVtxFT0VsVtxCol", "", kTH2F, {axisVtxZ, axisVtxZ});
338+
histos.add("noFT0activityNearby/hVtxFT0MinusVtxColVsMultT0M", "", kTH2F, {axisVtxZ, axisMultT0M});
339+
histos.add("noFT0activityNearby/nTracksPV_vs_V0A", "", kTH2F, {axisMultV0A, axisNtracks});
340+
histos.add("noFT0activityNearby/nTracksPV_vs_T0A", "", kTH2F, {axisMultFT0A, axisNtracks});
341+
histos.add("noFT0activityNearby/nTracksPV_vs_T0C", "", kTH2F, {axisMultFT0C, axisNtracks});
342+
histos.add("noFT0activityNearby/nTracksGlobal_vs_V0A", "", kTH2F, {axisMultV0A, axisNtracksGlobal});
343+
histos.add("noFT0activityNearby/nTracksGlobal_vs_T0A", "", kTH2F, {axisMultFT0A, axisNtracksGlobal});
344+
histos.add("noFT0activityNearby/nTracksGlobal_vs_T0C", "", kTH2F, {axisMultFT0C, axisNtracksGlobal});
345+
histos.add("noFT0activityNearby/hTVXvsBcDiff", "", kTH1F, {axisBcDiff});
346+
histos.add("noFT0activityNearby/hColTimeResVsNcontrib", "", kTH2F, {axisNcontrib, axisColTimeRes});
347+
histos.add("noFT0activityNearby/hColBcDiffVsNcontrib", "", kTH2F, {axisNcontrib, axisBcDiff});
348+
349+
histos.add("noPileup_cutByVzDiff_pvTOF_noFT0act/hBcColNoSel8", "", kTH1F, {axisBCs});
350+
histos.add("noPileup_cutByVzDiff_pvTOF_noFT0act/hBcTVX", "", kTH1F, {axisBCs});
351+
histos.add("noPileup_cutByVzDiff_pvTOF_noFT0act/hBcFT0", "", kTH1F, {axisBCs});
352+
histos.add("noPileup_cutByVzDiff_pvTOF_noFT0act/hBcFV0", "", kTH1F, {axisBCs});
353+
histos.add("noPileup_cutByVzDiff_pvTOF_noFT0act/hBcFDD", "", kTH1F, {axisBCs});
354+
histos.add("noPileup_cutByVzDiff_pvTOF_noFT0act/hBcZDC", "", kTH1F, {axisBCs});
355+
histos.add("noPileup_cutByVzDiff_pvTOF_noFT0act/hVtxFT0VsVtxCol", "", kTH2F, {axisVtxZ, axisVtxZ});
356+
histos.add("noPileup_cutByVzDiff_pvTOF_noFT0act/hVtxFT0MinusVtxColVsMultT0M", "", kTH2F, {axisVtxZ, axisMultT0M});
357+
histos.add("noPileup_cutByVzDiff_pvTOF_noFT0act/nTracksPV_vs_V0A", "", kTH2F, {axisMultV0A, axisNtracks});
358+
histos.add("noPileup_cutByVzDiff_pvTOF_noFT0act/nTracksPV_vs_T0A", "", kTH2F, {axisMultFT0A, axisNtracks});
359+
histos.add("noPileup_cutByVzDiff_pvTOF_noFT0act/nTracksPV_vs_T0C", "", kTH2F, {axisMultFT0C, axisNtracks});
360+
histos.add("noPileup_cutByVzDiff_pvTOF_noFT0act/nTracksGlobal_vs_V0A", "", kTH2F, {axisMultV0A, axisNtracksGlobal});
361+
histos.add("noPileup_cutByVzDiff_pvTOF_noFT0act/nTracksGlobal_vs_T0A", "", kTH2F, {axisMultFT0A, axisNtracksGlobal});
362+
histos.add("noPileup_cutByVzDiff_pvTOF_noFT0act/nTracksGlobal_vs_T0C", "", kTH2F, {axisMultFT0C, axisNtracksGlobal});
363+
histos.add("noPileup_cutByVzDiff_pvTOF_noFT0act/hTVXvsBcDiff", "", kTH1F, {axisBcDiff});
364+
histos.add("noPileup_cutByVzDiff_pvTOF_noFT0act/hColTimeResVsNcontrib", "", kTH2F, {axisNcontrib, axisColTimeRes});
365+
histos.add("noPileup_cutByVzDiff_pvTOF_noFT0act/hColBcDiffVsNcontrib", "", kTH2F, {axisNcontrib, axisBcDiff});
345366
}
346367

347368
Preslice<FullTracksIU> perCollision = aod::track::collisionId;
@@ -409,6 +430,7 @@ struct LightIonsEvSelQa {
409430

410431
std::vector<bool> vPastActivity(nBCs, 0);
411432
std::vector<bool> vFutureActivity(nBCs, 0);
433+
std::vector<bool> vNearbyFT0activity(nBCs, 0);
412434

413435
// create maps from globalBC to bc index for TVX or FT0-OR fired bcs
414436
// to be used for closest TVX (FT0-OR) searches
@@ -436,6 +458,9 @@ struct LightIonsEvSelQa {
436458
int deltaIndex = 0; // backward move counts
437459
int deltaBC = 0; // current difference wrt globalBC
438460
int maxDeltaBC = 30; // maximum difference
461+
462+
bool nearbyFT0activity = 0;
463+
439464
// past
440465
bool pastActivityFT0 = 0;
441466
bool pastActivityFDD = 0;
@@ -452,6 +477,12 @@ struct LightIonsEvSelQa {
452477
pastActivityFV0 |= bcPast.has_fv0a();
453478
pastActivityFDD |= bcPast.has_fdd();
454479
}
480+
if (deltaBC < 2) {
481+
if (bcPast.has_ft0()) {
482+
std::bitset<8> triggers = bcPast.ft0().triggerMask();
483+
nearbyFT0activity |= triggers[o2::ft0::RecPoints::ETriggerBits::kIsActiveSideA];
484+
}
485+
}
455486
}
456487
bool pastActivity = pastActivityFT0 | pastActivityFV0 | pastActivityFDD;
457488
vPastActivity[indexBc] = pastActivity;
@@ -474,9 +505,16 @@ struct LightIonsEvSelQa {
474505
futureActivityFV0 |= bcFuture.has_fv0a();
475506
futureActivityFDD |= bcFuture.has_fdd();
476507
}
508+
if (deltaBC < 2) {
509+
if (bcFuture.has_ft0()) {
510+
std::bitset<8> triggers = bcFuture.ft0().triggerMask();
511+
nearbyFT0activity |= triggers[o2::ft0::RecPoints::ETriggerBits::kIsActiveSideA];
512+
}
513+
}
477514
}
478515
bool futureActivity = futureActivityFT0 | futureActivityFV0 | futureActivityFDD;
479516
vFutureActivity[indexBc] = futureActivity;
517+
vNearbyFT0activity[indexBc] = nearbyFT0activity;
480518

481519
// monitor BCs with nearby activity:
482520

@@ -556,7 +594,7 @@ struct LightIonsEvSelQa {
556594
// ### count tracks of different types
557595
int nPVtracks = 0;
558596
int nGlobalTracks = 0;
559-
int nTOFtracks = 0;
597+
// int nTOFtracks = 0;
560598
auto tracksGrouped = tracks.sliceBy(perCollision, col.globalIndex());
561599
for (const auto& track : tracksGrouped) {
562600
if (!track.isPVContributor()) {
@@ -570,7 +608,7 @@ struct LightIonsEvSelQa {
570608
continue;
571609

572610
nPVtracks++;
573-
nTOFtracks += track.hasTOF();
611+
// nTOFtracks += track.hasTOF();
574612

575613
if (track.hasITS() && track.hasTPC() && track.tpcNClsFound() > 50 && track.tpcNClsCrossedRows() > 50 && track.tpcChi2NCl() < 4)
576614
nGlobalTracks++;
@@ -579,6 +617,15 @@ struct LightIonsEvSelQa {
579617
bool hasFT0 = foundBC.has_ft0();
580618
bool hasFV0A = foundBC.has_fv0a();
581619

620+
// bool noFT0activityNearby = false;
621+
bool noFT0activityNearby = (vNearbyFT0activity[bcIndex] == 0);
622+
// check kIsFlangeEvent
623+
if (hasFT0) {
624+
std::bitset<8> triggers = foundBC.ft0().triggerMask();
625+
if (triggers[o2::ft0::RecPoints::ETriggerBits::kIsFlangeEvent])
626+
noFT0activityNearby = false;
627+
}
628+
582629
float vZ = col.posZ();
583630
float vZft0 = hasFT0 ? foundBC.ft0().posZ() : -1000;
584631
float diffVz = vZft0 - vZ;
@@ -656,6 +703,9 @@ struct LightIonsEvSelQa {
656703
if (noPastActivity) {
657704
histos.fill(HIST("noPastActivity/hBcColNoSel8"), localBC);
658705
}
706+
if (noFT0activityNearby) {
707+
histos.fill(HIST("noFT0activityNearby/hBcColNoSel8"), localBC);
708+
}
659709
if (badVzDiff) {
660710
histos.fill(HIST("badVzDiff/hBcColNoSel8"), localBC);
661711
}
@@ -674,8 +724,8 @@ struct LightIonsEvSelQa {
674724
if (noPU && grassOnTheRight) {
675725
histos.fill(HIST("noPileup_HighMultCloudCut/hBcColNoSel8"), localBC);
676726
}
677-
if (noPU && pvTOFmatched && !badVzDiff) { // noPileup_cutByVzDiff_pvTOF
678-
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF/hBcColNoSel8"), localBC);
727+
if (noPU && pvTOFmatched && !badVzDiff && noFT0activityNearby) { // noPileup_cutByVzDiff_pvTOF_noFT0act
728+
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF_noFT0act/hBcColNoSel8"), localBC);
679729
}
680730

681731
// only here cut on sel8:
@@ -706,6 +756,11 @@ struct LightIonsEvSelQa {
706756
histos.fill(HIST("noPastActivity/hColBcDiffVsNcontrib"), nContributors, bcToClosestTVXdiff);
707757
histos.fill(HIST("noPastActivity/hColTimeResVsNcontrib"), nContributors, timeRes);
708758
}
759+
if (noFT0activityNearby) {
760+
histos.fill(HIST("noFT0activityNearby/hBcTVX"), localBC);
761+
histos.fill(HIST("noFT0activityNearby/hColBcDiffVsNcontrib"), nContributors, bcToClosestTVXdiff);
762+
histos.fill(HIST("noFT0activityNearby/hColTimeResVsNcontrib"), nContributors, timeRes);
763+
}
709764
if (badVzDiff) {
710765
histos.fill(HIST("badVzDiff/hBcTVX"), localBC);
711766
histos.fill(HIST("badVzDiff/hColBcDiffVsNcontrib"), nContributors, bcToClosestTVXdiff);
@@ -736,10 +791,10 @@ struct LightIonsEvSelQa {
736791
histos.fill(HIST("noPileup_HighMultCloudCut/hColBcDiffVsNcontrib"), nContributors, bcToClosestTVXdiff);
737792
histos.fill(HIST("noPileup_HighMultCloudCut/hColTimeResVsNcontrib"), nContributors, timeRes);
738793
}
739-
if (noPU && pvTOFmatched && !badVzDiff) { // noPileup_cutByVzDiff_pvTOF
740-
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF/hBcTVX"), localBC);
741-
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF/hColBcDiffVsNcontrib"), nContributors, bcToClosestTVXdiff);
742-
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF/hColTimeResVsNcontrib"), nContributors, timeRes);
794+
if (noPU && pvTOFmatched && !badVzDiff && noFT0activityNearby) { // noPileup_cutByVzDiff_pvTOF_noFT0act
795+
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF_noFT0act/hBcTVX"), localBC);
796+
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF_noFT0act/hColBcDiffVsNcontrib"), nContributors, bcToClosestTVXdiff);
797+
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF_noFT0act/hColTimeResVsNcontrib"), nContributors, timeRes);
743798
}
744799

745800
if (foundBC.has_ft0()) {
@@ -808,6 +863,15 @@ struct LightIonsEvSelQa {
808863
histos.fill(HIST("noPastActivity/nTracksPV_vs_T0C"), multT0C, nPVtracks);
809864
histos.fill(HIST("noPastActivity/nTracksGlobal_vs_T0C"), multT0C, nGlobalTracks);
810865
}
866+
if (noFT0activityNearby) {
867+
histos.fill(HIST("noFT0activityNearby/hBcFT0"), localBC);
868+
histos.fill(HIST("noFT0activityNearby/hVtxFT0VsVtxCol"), vZft0, vZ);
869+
histos.fill(HIST("noFT0activityNearby/hVtxFT0MinusVtxColVsMultT0M"), diffVz, multT0A + multT0C);
870+
histos.fill(HIST("noFT0activityNearby/nTracksPV_vs_T0A"), multT0A, nPVtracks);
871+
histos.fill(HIST("noFT0activityNearby/nTracksGlobal_vs_T0A"), multT0A, nGlobalTracks);
872+
histos.fill(HIST("noFT0activityNearby/nTracksPV_vs_T0C"), multT0C, nPVtracks);
873+
histos.fill(HIST("noFT0activityNearby/nTracksGlobal_vs_T0C"), multT0C, nGlobalTracks);
874+
}
811875
if (narrowDeltaTimeVeto) {
812876
histos.fill(HIST("narrowTimeVeto/hBcFT0"), localBC);
813877
histos.fill(HIST("narrowTimeVeto/hVtxFT0VsVtxCol"), vZft0, vZ);
@@ -844,14 +908,14 @@ struct LightIonsEvSelQa {
844908
histos.fill(HIST("noPileup_HighMultCloudCut/nTracksPV_vs_T0C"), multT0C, nPVtracks);
845909
histos.fill(HIST("noPileup_HighMultCloudCut/nTracksGlobal_vs_T0C"), multT0C, nGlobalTracks);
846910
}
847-
if (noPU && pvTOFmatched && !badVzDiff) { // noPileup_cutByVzDiff_pvTOF
848-
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF/hBcFT0"), localBC);
849-
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF/hVtxFT0VsVtxCol"), vZft0, vZ);
850-
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF/hVtxFT0MinusVtxColVsMultT0M"), diffVz, multT0A + multT0C);
851-
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF/nTracksPV_vs_T0A"), multT0A, nPVtracks);
852-
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF/nTracksGlobal_vs_T0A"), multT0A, nGlobalTracks);
853-
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF/nTracksPV_vs_T0C"), multT0C, nPVtracks);
854-
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF/nTracksGlobal_vs_T0C"), multT0C, nGlobalTracks);
911+
if (noPU && pvTOFmatched && !badVzDiff && noFT0activityNearby) { // noPileup_cutByVzDiff_pvTOF_noFT0act
912+
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF_noFT0act/hBcFT0"), localBC);
913+
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF_noFT0act/hVtxFT0VsVtxCol"), vZft0, vZ);
914+
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF_noFT0act/hVtxFT0MinusVtxColVsMultT0M"), diffVz, multT0A + multT0C);
915+
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF_noFT0act/nTracksPV_vs_T0A"), multT0A, nPVtracks);
916+
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF_noFT0act/nTracksGlobal_vs_T0A"), multT0A, nGlobalTracks);
917+
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF_noFT0act/nTracksPV_vs_T0C"), multT0C, nPVtracks);
918+
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF_noFT0act/nTracksGlobal_vs_T0C"), multT0C, nGlobalTracks);
855919
}
856920
}
857921

@@ -889,6 +953,11 @@ struct LightIonsEvSelQa {
889953
histos.fill(HIST("noPastActivity/nTracksPV_vs_V0A"), multV0A, nPVtracks);
890954
histos.fill(HIST("noPastActivity/nTracksGlobal_vs_V0A"), multV0A, nGlobalTracks);
891955
}
956+
if (noFT0activityNearby) {
957+
histos.fill(HIST("noFT0activityNearby/hBcFV0"), localBC);
958+
histos.fill(HIST("noFT0activityNearby/nTracksPV_vs_V0A"), multV0A, nPVtracks);
959+
histos.fill(HIST("noFT0activityNearby/nTracksGlobal_vs_V0A"), multV0A, nGlobalTracks);
960+
}
892961
if (narrowDeltaTimeVeto) {
893962
histos.fill(HIST("narrowTimeVeto/hBcFV0"), localBC);
894963
histos.fill(HIST("narrowTimeVeto/nTracksPV_vs_V0A"), multV0A, nPVtracks);
@@ -909,10 +978,10 @@ struct LightIonsEvSelQa {
909978
histos.fill(HIST("noPileup_HighMultCloudCut/nTracksPV_vs_V0A"), multV0A, nPVtracks);
910979
histos.fill(HIST("noPileup_HighMultCloudCut/nTracksGlobal_vs_V0A"), multV0A, nGlobalTracks);
911980
}
912-
if (noPU && pvTOFmatched && !badVzDiff) { // noPileup_cutByVzDiff_pvTOF
913-
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF/hBcFV0"), localBC);
914-
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF/nTracksPV_vs_V0A"), multV0A, nPVtracks);
915-
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF/nTracksGlobal_vs_V0A"), multV0A, nGlobalTracks);
981+
if (noPU && pvTOFmatched && !badVzDiff && noFT0activityNearby) { // noPileup_cutByVzDiff_pvTOF_noFT0act
982+
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF_noFT0act/hBcFV0"), localBC);
983+
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF_noFT0act/nTracksPV_vs_V0A"), multV0A, nPVtracks);
984+
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF_noFT0act/nTracksGlobal_vs_V0A"), multV0A, nGlobalTracks);
916985
}
917986
}
918987
if (foundBC.has_zdc()) {
@@ -935,6 +1004,9 @@ struct LightIonsEvSelQa {
9351004
if (noPastActivity) {
9361005
histos.fill(HIST("noPastActivity/hBcZDC"), localBC);
9371006
}
1007+
if (noFT0activityNearby) {
1008+
histos.fill(HIST("noFT0activityNearby/hBcZDC"), localBC);
1009+
}
9381010
if (narrowDeltaTimeVeto) {
9391011
histos.fill(HIST("narrowTimeVeto/hBcZDC"), localBC);
9401012
}
@@ -947,8 +1019,8 @@ struct LightIonsEvSelQa {
9471019
if (noPU && grassOnTheRight) {
9481020
histos.fill(HIST("noPileup_HighMultCloudCut/hBcZDC"), localBC);
9491021
}
950-
if (noPU && pvTOFmatched && !badVzDiff) { // noPileup_cutByVzDiff_pvTOF
951-
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF/hBcZDC"), localBC);
1022+
if (noPU && pvTOFmatched && !badVzDiff && noFT0activityNearby) { // noPileup_cutByVzDiff_pvTOF_noFT0act
1023+
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF_noFT0act/hBcZDC"), localBC);
9521024
}
9531025
}
9541026

@@ -976,6 +1048,9 @@ struct LightIonsEvSelQa {
9761048
if (noPastActivity) {
9771049
histos.fill(HIST("noPastActivity/hTVXvsBcDiff"), bcDiff);
9781050
}
1051+
if (noFT0activityNearby) {
1052+
histos.fill(HIST("noFT0activityNearby/hTVXvsBcDiff"), bcDiff);
1053+
}
9791054
if (narrowDeltaTimeVeto) {
9801055
histos.fill(HIST("narrowTimeVeto/hTVXvsBcDiff"), bcDiff);
9811056
}
@@ -988,8 +1063,8 @@ struct LightIonsEvSelQa {
9881063
if (noPU && grassOnTheRight) {
9891064
histos.fill(HIST("noPileup_HighMultCloudCut/hTVXvsBcDiff"), bcDiff);
9901065
}
991-
if (noPU && pvTOFmatched && !badVzDiff) { // noPileup_cutByVzDiff_pvTOF
992-
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF/hTVXvsBcDiff"), bcDiff);
1066+
if (noPU && pvTOFmatched && !badVzDiff && noFT0activityNearby) { // noPileup_cutByVzDiff_pvTOF_noFT0act
1067+
histos.fill(HIST("noPileup_cutByVzDiff_pvTOF_noFT0act/hTVXvsBcDiff"), bcDiff);
9931068
}
9941069

9951070
} // end of collisions loop

0 commit comments

Comments
 (0)