Skip to content

Updated EMCAL event selection for Full Jet Finders #7345

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
mfasDa merged 5 commits into from
Aug 19, 2024
Merged

Updated EMCAL event selection for Full Jet Finders #7345

Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
2c65d83
updated EMCAL event selection for all the Jets process functions
Archita-Dash Aug 18, 2024
555df28
Merge branch 'master' of https://github.com/AliceO2Group/O2Physics
Archita-Dash Aug 18, 2024
fc9474a
updated EMCAL event selection for all the Jets process functions
Archita-Dash Aug 18, 2024
0cb953b
Please consider the following formatting changes
alibuild Aug 18, 2024
5e67b1f
Merge pull request #7 from alibuild/alibot-cleanup-7345
Archita-Dash Aug 18, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
169 changes: 144 additions & 25 deletions PWGJE/Tasks/fulljetspectrapp.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -117,14 +117,24 @@ struct FullJetSpectrapp {
{
auto h_collisions_unweighted = registry.get<TH1>(HIST("h_collisions_unweighted"));
h_collisions_unweighted->GetXaxis()->SetBinLabel(2, "total events");
h_collisions_unweighted->GetXaxis()->SetBinLabel(3, "EMC events with kTVXinEMC");
h_collisions_unweighted->GetXaxis()->SetBinLabel(4, "EMC events w/o kTVXinEMC");
h_collisions_unweighted->GetXaxis()->SetBinLabel(3, "JetsData with kTVXinEMC");
h_collisions_unweighted->GetXaxis()->SetBinLabel(4, "JetsMCD with kTVXinEMC");
h_collisions_unweighted->GetXaxis()->SetBinLabel(5, "Tracks with kTVXinEMC");
h_collisions_unweighted->GetXaxis()->SetBinLabel(6, "JetsMCPMCDMatched with kTVXinEMC");
h_collisions_unweighted->GetXaxis()->SetBinLabel(7, "JetsData w/o kTVXinEMC");
h_collisions_unweighted->GetXaxis()->SetBinLabel(8, "JetsMCD w/o kTVXinEMC");
h_collisions_unweighted->GetXaxis()->SetBinLabel(9, "Tracks w/o kTVXinEMC");
h_collisions_unweighted->GetXaxis()->SetBinLabel(10, "JetsMCPMCDMatched w/o kTVXinEMC");

if (doprocessTracksWeighted) {
auto h_collisions_weighted = registry.get<TH1>(HIST("h_collisions_weighted"));
h_collisions_weighted->GetXaxis()->SetBinLabel(2, "total events");
h_collisions_weighted->GetXaxis()->SetBinLabel(3, "EMC events with kTVXinEMC");
h_collisions_weighted->GetXaxis()->SetBinLabel(4, "EMC events w/o kTVXinEMC");
h_collisions_weighted->GetXaxis()->SetBinLabel(3, "JetsMCDWeighted with kTVXinEMC");
h_collisions_weighted->GetXaxis()->SetBinLabel(4, "JetsMCPMCDMatchedWeighted with kTVXinEMC");
h_collisions_weighted->GetXaxis()->SetBinLabel(5, "TracksWeighted with kTVXinEMC");
h_collisions_weighted->GetXaxis()->SetBinLabel(6, "JetsMCDWeighted w/o kTVXinEMC");
h_collisions_weighted->GetXaxis()->SetBinLabel(7, "JetsMCPMCDMatchedWeighted w/o kTVXinEMC");
h_collisions_weighted->GetXaxis()->SetBinLabel(8, "TracksWeighted w/o kTVXinEMC");
}
}

Expand All @@ -147,7 +157,7 @@ struct FullJetSpectrapp {

// JetTrack QA histograms
if (doprocessTracks || doprocessTracksWeighted) {
registry.add("h_collisions_unweighted", "event status; event status;entries", {HistType::kTH1F, {{4, 0., 4.0}}});
registry.add("h_collisions_unweighted", "event status; event status;entries", {HistType::kTH1F, {{11, 0., 11.0}}});

registry.add("h_track_pt", "track pT;#it{p}_{T,track} (GeV/#it{c});entries", {HistType::kTH1F, {{350, 0., 350.}}});
registry.add("h_track_eta", "track #eta;#eta_{track};entries", {HistType::kTH1F, {{100, -1., 1.}}});
Expand All @@ -172,7 +182,7 @@ struct FullJetSpectrapp {
// registry.add("h_gaptrig_cluster_energy", "gap triggered cluster #varphi;#varphi_{cluster};entries", {HistType::kTH1F, {{160, 0., 7.}}});

if (doprocessTracksWeighted) {
registry.add("h_collisions_weighted", "event status;event status;entries", {HistType::kTH1F, {{5, 0.0, 5.0}}});
registry.add("h_collisions_weighted", "event status;event status;entries", {HistType::kTH1F, {{10, 0.0, 10.0}}});
registry.add("h_gaptrig_collisions", "event status; event status; entries", {HistType::kTH1F, {{4, 0.0, 4.0}}});

// registry.add("h_gaptrig_track_pt", "gap triggered track pT;#it{p}_{T,track} (GeV/#it{c});entries", {HistType::kTH1F, {{350, 0., 350.}}});
Expand Down Expand Up @@ -394,8 +404,30 @@ struct FullJetSpectrapp {
}
PROCESS_SWITCH(FullJetSpectrapp, processDummy, "dummy task", true);

void processJetsData(soa::Filtered<EMCCollisions>::iterator const&, FullJetTableDataJoined const& jets, JetTracks const&, JetClusters const&)
void processJetsData(soa::Filtered<EMCCollisions>::iterator const& collision, FullJetTableDataJoined const& jets, JetTracks const&, JetClusters const&)
{
registry.fill(HIST("h_collisions_unweighted"), 1.0);
bool eventAccepted = false;

if (doEMCALEventWorkaround) {
if (collision.isEmcalReadout() && !collision.isAmbiguous()) { // i.e. EMCAL has a cell content
eventAccepted = true;
if (collision.alias_bit(kTVXinEMC)) {
registry.fill(HIST("h_collisions_unweighted"), 2.0);
}
}
} else {
if (!collision.isAmbiguous() && jetderiveddatautilities::eventEMCAL(collision) && collision.alias_bit(kTVXinEMC)) {
eventAccepted = true;
registry.fill(HIST("h_collisions_unweighted"), 2.0);
}
}

if (!eventAccepted) {
registry.fill(HIST("h_collisions_unweighted"), 6.0);
return;
}

for (auto const& jet : jets) {
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
continue;
Expand All @@ -408,8 +440,30 @@ struct FullJetSpectrapp {
}
PROCESS_SWITCH(FullJetSpectrapp, processJetsData, "Full Jets Data", false);

void processJetsMCD(soa::Filtered<EMCCollisions>::iterator const&, JetTableMCDJoined const& jets, JetTracks const&, JetClusters const&)
void processJetsMCD(soa::Filtered<EMCCollisions>::iterator const& collision, JetTableMCDJoined const& jets, JetTracks const&, JetClusters const&)
{
registry.fill(HIST("h_collisions_unweighted"), 1.0);
bool eventAccepted = false;

if (doEMCALEventWorkaround) {
if (collision.isEmcalReadout() && !collision.isAmbiguous()) { // i.e. EMCAL has a cell content
eventAccepted = true;
if (collision.alias_bit(kTVXinEMC)) {
registry.fill(HIST("h_collisions_unweighted"), 3.0);
}
}
} else {
if (!collision.isAmbiguous() && jetderiveddatautilities::eventEMCAL(collision) && collision.alias_bit(kTVXinEMC)) {
eventAccepted = true;
registry.fill(HIST("h_collisions_unweighted"), 3.0);
}
}

if (!eventAccepted) {
registry.fill(HIST("h_collisions_unweighted"), 7.0);
return;
}

for (auto const& jet : jets) {
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
continue;
Expand All @@ -422,8 +476,30 @@ struct FullJetSpectrapp {
}
PROCESS_SWITCH(FullJetSpectrapp, processJetsMCD, "Full Jets at Detector Level", false);

void processJetsMCDWeighted(soa::Filtered<EMCCollisions>::iterator const&, JetTableMCDWeightedJoined const& jets, JetTracks const&, JetClusters const&)
void processJetsMCDWeighted(soa::Filtered<EMCCollisions>::iterator const& collision, JetTableMCDWeightedJoined const& jets, JetTracks const&, JetClusters const&)
{
registry.fill(HIST("h_collisions_weighted"), 1.0);
bool eventAccepted = false;

if (doEMCALEventWorkaround) {
if (collision.isEmcalReadout() && !collision.isAmbiguous()) { // i.e. EMCAL has a cell content
eventAccepted = true;
if (collision.alias_bit(kTVXinEMC)) {
registry.fill(HIST("h_collisions_weighted"), 2.0);
}
}
} else {
if (!collision.isAmbiguous() && jetderiveddatautilities::eventEMCAL(collision) && collision.alias_bit(kTVXinEMC)) {
eventAccepted = true;
registry.fill(HIST("h_collisions_weighted"), 2.0);
}
}

if (!eventAccepted) {
registry.fill(HIST("h_collisions_weighted"), 5.0);
return;
}

for (auto const& jet : jets) {
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
continue;
Expand All @@ -434,9 +510,9 @@ struct FullJetSpectrapp {
fillJetHistograms(jet, jet.eventWeight());
}
}
PROCESS_SWITCH(FullJetSpectrapp, processJetsMCDWeighted, "jet finder HF QA mcd on weighted events", false);
PROCESS_SWITCH(FullJetSpectrapp, processJetsMCDWeighted, "Full Jets at Detector Level on weighted events", false);

void processJetsMCP(typename JetTableMCPJoined::iterator const& jet, JetParticles const&)
void processJetsMCP(typename JetTableMCPJoined::iterator const& jet, JetParticles const&, JetMcCollisions const&)
{
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
return;
Expand All @@ -458,7 +534,7 @@ struct FullJetSpectrapp {
}
fillMCPHistograms(jet, jet.eventWeight());
}
PROCESS_SWITCH(FullJetSpectrapp, processJetsMCPWeighted, "jet finder HF QA mcp on weighted events", false);
PROCESS_SWITCH(FullJetSpectrapp, processJetsMCPWeighted, "Full Jets at Particle Level on weighted events", false);

void processTracks(soa::Filtered<EMCCollisions>::iterator const& collision, soa::Filtered<JetTracks> const& tracks, soa::Filtered<JetClusters> const& clusters)
{
Expand All @@ -475,29 +551,51 @@ struct FullJetSpectrapp {
if (collision.isEmcalReadout() && !collision.isAmbiguous()) { // i.e. EMCAL has a cell content
eventAccepted = true;
if (collision.alias_bit(kTVXinEMC)) {
registry.fill(HIST("h_collisions_unweighted"), 2.0);
registry.fill(HIST("h_collisions_unweighted"), 4.0);
}
}
} else {
// Check if EMCAL was readout with the MB trigger(kTVXinEMC) fired. If not then reject the event and exit the function.
// This is the default check for the simulations with proper trigger flags not requiring the above workaround.
if (!collision.isAmbiguous() && jetderiveddatautilities::eventEMCAL(collision) && collision.alias_bit(kTVXinEMC)) {
eventAccepted = true;
registry.fill(HIST("h_collisions_unweighted"), 2.0);
registry.fill(HIST("h_collisions_unweighted"), 4.0);
}
}

if (!eventAccepted) {
registry.fill(HIST("h_collisions_unweighted"), 3.0);
registry.fill(HIST("h_collisions_unweighted"), 8.0);
return;
}
// Fill Accepted events histos
fillTrackHistograms(tracks, clusters, 1.0);
}
PROCESS_SWITCH(FullJetSpectrapp, processTracks, "QA for fulljet tracks", false);
PROCESS_SWITCH(FullJetSpectrapp, processTracks, "Full Jet tracks", false);

void processJetsMCPMCDMatched(soa::Filtered<EMCCollisions>::iterator const&, JetTableMCDMatchedJoined const& mcdjets, JetTableMCPMatchedJoined const&, JetTracks const&, JetClusters const&, JetParticles const&)
void processJetsMCPMCDMatched(soa::Filtered<EMCCollisions>::iterator const& collision, JetTableMCDMatchedJoined const& mcdjets, JetTableMCPMatchedJoined const&, JetTracks const&, JetClusters const&, JetParticles const&)
{
registry.fill(HIST("h_collisions_unweighted"), 1.0);
bool eventAccepted = false;

if (doEMCALEventWorkaround) {
if (collision.isEmcalReadout() && !collision.isAmbiguous()) { // i.e. EMCAL has a cell content
eventAccepted = true;
if (collision.alias_bit(kTVXinEMC)) {
registry.fill(HIST("h_collisions_unweighted"), 5.0);
}
}
} else {
if (!collision.isAmbiguous() && jetderiveddatautilities::eventEMCAL(collision) && collision.alias_bit(kTVXinEMC)) {
eventAccepted = true;
registry.fill(HIST("h_collisions_unweighted"), 5.0);
}
}

if (!eventAccepted) {
registry.fill(HIST("h_collisions_unweighted"), 9.0);
return;
}

for (const auto& mcdjet : mcdjets) {
if (!jetfindingutilities::isInEtaAcceptance(mcdjet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
continue;
Expand All @@ -508,10 +606,32 @@ struct FullJetSpectrapp {
fillMatchedHistograms<typename JetTableMCDMatchedJoined::iterator, JetTableMCPMatchedJoined>(mcdjet);
}
}
PROCESS_SWITCH(FullJetSpectrapp, processJetsMCPMCDMatched, "full jet finder MCP matched to MCD", false);
PROCESS_SWITCH(FullJetSpectrapp, processJetsMCPMCDMatched, "Full Jet finder MCP matched to MCD", false);

void processJetsMCPMCDMatchedWeighted(soa::Filtered<EMCCollisions>::iterator const&, JetTableMCDMatchedWeightedJoined const& mcdjets, JetTableMCPMatchedWeightedJoined const&, JetTracks const&, JetClusters const&, JetParticles const&)
void processJetsMCPMCDMatchedWeighted(soa::Filtered<EMCCollisions>::iterator const& collision, JetTableMCDMatchedWeightedJoined const& mcdjets, JetTableMCPMatchedWeightedJoined const&, JetTracks const&, JetClusters const&, JetParticles const&)
{
registry.fill(HIST("h_collisions_weighted"), 1.0);
bool eventAccepted = false;

if (doEMCALEventWorkaround) {
if (collision.isEmcalReadout() && !collision.isAmbiguous()) { // i.e. EMCAL has a cell content
eventAccepted = true;
if (collision.alias_bit(kTVXinEMC)) {
registry.fill(HIST("h_collisions_weighted"), 3.0);
}
}
} else {
if (!collision.isAmbiguous() && jetderiveddatautilities::eventEMCAL(collision) && collision.alias_bit(kTVXinEMC)) {
eventAccepted = true;
registry.fill(HIST("h_collisions_weighted"), 3.0);
}
}

if (!eventAccepted) {
registry.fill(HIST("h_collisions_weighted"), 6.0);
return;
}

for (const auto& mcdjet : mcdjets) {
if (!jetfindingutilities::isInEtaAcceptance(mcdjet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
continue;
Expand All @@ -522,7 +642,7 @@ struct FullJetSpectrapp {
fillMatchedHistograms<typename JetTableMCDMatchedWeightedJoined::iterator, JetTableMCPMatchedWeightedJoined>(mcdjet, mcdjet.eventWeight());
}
}
PROCESS_SWITCH(FullJetSpectrapp, processJetsMCPMCDMatchedWeighted, "full jet finder MCP matched to MCD on weighted events", false);
PROCESS_SWITCH(FullJetSpectrapp, processJetsMCPMCDMatchedWeighted, "Full Jet finder MCP matched to MCD on weighted events", false);

void processTracksWeighted(soa::Filtered<soa::Join<EMCCollisions, aod::JMcCollisionLbs>>::iterator const& collision,
aod::JMcCollisions const&,
Expand All @@ -531,7 +651,6 @@ struct FullJetSpectrapp {
{
bool eventAccepted = false;
float eventWeight = collision.mcCollision().weight();
// registry.fill(HIST("h_collisions_unweighted"), 1.0);
registry.fill(HIST("h_collisions_weighted"), 1.0, eventWeight);

// set "doMBGapTrigger" to true only if you are testing with MB Gap-triggers
Expand All @@ -544,26 +663,26 @@ struct FullJetSpectrapp {
eventAccepted = true;
fillTrackHistograms(tracks, clusters, eventWeight);
if (collision.alias_bit(kTVXinEMC)) {
registry.fill(HIST("h_collisions_weighted"), 2.0, eventWeight);
registry.fill(HIST("h_collisions_weighted"), 4.0, eventWeight);
}
}
} else {
// Check if EMCAL was readout with the MB trigger(kTVXinEMC) fired. If not then reject the event and exit the function.
// This is the default check for the simulations with proper trigger flags not requiring the above workaround.
if (!collision.isAmbiguous() && jetderiveddatautilities::eventEMCAL(collision) && collision.alias_bit(kTVXinEMC)) {
eventAccepted = true;
registry.fill(HIST("h_collisions_weighted"), 2.0, eventWeight);
registry.fill(HIST("h_collisions_weighted"), 4.0, eventWeight);
}
}

if (!eventAccepted) {
registry.fill(HIST("h_collisions_weighted"), 3.0, eventWeight);
registry.fill(HIST("h_collisions_weighted"), 7.0, eventWeight);
return;
}
// registry.fill(HIST("h_gaptrig_collisions"), 1.0, eventWeight);
fillTrackHistograms(tracks, clusters, eventWeight);
}
PROCESS_SWITCH(FullJetSpectrapp, processTracksWeighted, "QA for fulljet tracks weighted", false);
PROCESS_SWITCH(FullJetSpectrapp, processTracksWeighted, "Full Jet tracks weighted", false);

}; // struct

Expand Down