From 51dc87f4488f7e90bbac9e9f2e303ba08ac3ae3e Mon Sep 17 00:00:00 2001 From: jaimenorman Date: Tue, 14 Jan 2025 16:10:00 +0100 Subject: [PATCH 1/3] add MB gap event rejection --- PWGJE/Tasks/jetFinderQA.cxx | 16 ++++++++++++++++ 1 file changed, 16 insertions(+) diff --git a/PWGJE/Tasks/jetFinderQA.cxx b/PWGJE/Tasks/jetFinderQA.cxx index 67806db1115..45687620736 100644 --- a/PWGJE/Tasks/jetFinderQA.cxx +++ b/PWGJE/Tasks/jetFinderQA.cxx @@ -72,6 +72,7 @@ struct JetFinderQATask { Configurable checkMcCollisionIsMatched{"checkMcCollisionIsMatched", false, "0: count whole MCcollisions, 1: select MCcollisions which only have their correspond collisions"}; Configurable trackOccupancyInTimeRangeMax{"trackOccupancyInTimeRangeMax", 999999, "maximum occupancy of tracks in neighbouring collisions in a given time range; only applied to reconstructed collisions (data and mcd jets), not mc collisions (mcp jets)"}; Configurable trackOccupancyInTimeRangeMin{"trackOccupancyInTimeRangeMin", -999999, "minimum occupancy of tracks in neighbouring collisions in a given time range; only applied to reconstructed collisions (data and mcd jets), not mc collisions (mcp jets)"}; + Configurable rejectMbEvents{"rejectMbEvents", false, "flag to choose to reject min. bias events based on weight (MB = 1, weighted = anything else)"}; std::vector filledJetR_Both; std::vector filledJetR_Low; @@ -827,6 +828,9 @@ struct JetFinderQATask { if (!isAcceptedJet(jet)) { continue; } + if(rejectMbEvents && jet.eventWeight()==1) { + continue; + } double pTHat = 10. / (std::pow(jet.eventWeight(), 1.0 / pTHatExponent)); for (int N = 1; N < 21; N++) { if (jet.pt() < N * 0.25 * pTHat && jet.r() == round(selectedJetsRadius * 100.0f)) { @@ -865,6 +869,9 @@ struct JetFinderQATask { if (!isAcceptedJet(jet)) { return; } + if(rejectMbEvents && jet.eventWeight()==1) { + return; + } double pTHat = 10. / (std::pow(jet.eventWeight(), 1.0 / pTHatExponent)); for (int N = 1; N < 21; N++) { if (jet.pt() < N * 0.25 * pTHat && jet.r() == round(selectedJetsRadius * 100.0f)) { @@ -917,6 +924,9 @@ struct JetFinderQATask { if (!isAcceptedJet(mcdjet)) { continue; } + if(rejectMbEvents && mcdjet.eventWeight()==1) { + continue; + } fillMatchedHistograms::iterator, soa::Join>(mcdjet, mcdjet.eventWeight()); } } @@ -924,6 +934,9 @@ struct JetFinderQATask { void processMCCollisionsWeighted(aod::JetMcCollision const& collision) { + if(rejectMbEvents && collision.weight()==1) { + return; + } registry.fill(HIST("h_collision_eventweight_part"), collision.weight()); } PROCESS_SWITCH(JetFinderQATask, processMCCollisionsWeighted, "collision QA for weighted events", false); @@ -1078,6 +1091,9 @@ struct JetFinderQATask { soa::Filtered> const& tracks) { float eventWeight = collision.mcCollision().weight(); + if(rejectMbEvents && eventWeight==1) { + return; + } registry.fill(HIST("h_collisions"), 0.5); registry.fill(HIST("h_collisions_weighted"), 0.5, eventWeight); if (!jetderiveddatautilities::selectCollision(collision, eventSelection)) { From 40a39f006e16fa0e97342b8f06b8c5476fa19b0d Mon Sep 17 00:00:00 2001 From: jaimenorman Date: Tue, 14 Jan 2025 16:24:00 +0100 Subject: [PATCH 2/3] fix event matching for particle level weighted process function --- PWGJE/Tasks/jetFinderQA.cxx | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/PWGJE/Tasks/jetFinderQA.cxx b/PWGJE/Tasks/jetFinderQA.cxx index 45687620736..b0809330f41 100644 --- a/PWGJE/Tasks/jetFinderQA.cxx +++ b/PWGJE/Tasks/jetFinderQA.cxx @@ -880,7 +880,7 @@ struct JetFinderQATask { } if (checkMcCollisionIsMatched) { auto collisionspermcpjet = collisions.sliceBy(CollisionsPerMCPCollision, jet.mcCollisionId()); - if (collisionspermcpjet.size() >= 1) { + if (collisionspermcpjet.size() >= 1 && jetderiveddatautilities::selectCollision(collisionspermcpjet.begin(), eventSelection)) { fillMCPHistograms(jet, jet.eventWeight()); } } else { From b25f1faa72ea46d3822d6c8d4e9279724930269c Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Wed, 15 Jan 2025 10:25:29 +0000 Subject: [PATCH 3/3] Please consider the following formatting changes --- PWGJE/Tasks/jetFinderQA.cxx | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/PWGJE/Tasks/jetFinderQA.cxx b/PWGJE/Tasks/jetFinderQA.cxx index b0809330f41..50d57adb870 100644 --- a/PWGJE/Tasks/jetFinderQA.cxx +++ b/PWGJE/Tasks/jetFinderQA.cxx @@ -828,7 +828,7 @@ struct JetFinderQATask { if (!isAcceptedJet(jet)) { continue; } - if(rejectMbEvents && jet.eventWeight()==1) { + if (rejectMbEvents && jet.eventWeight() == 1) { continue; } double pTHat = 10. / (std::pow(jet.eventWeight(), 1.0 / pTHatExponent)); @@ -869,7 +869,7 @@ struct JetFinderQATask { if (!isAcceptedJet(jet)) { return; } - if(rejectMbEvents && jet.eventWeight()==1) { + if (rejectMbEvents && jet.eventWeight() == 1) { return; } double pTHat = 10. / (std::pow(jet.eventWeight(), 1.0 / pTHatExponent)); @@ -880,7 +880,7 @@ struct JetFinderQATask { } if (checkMcCollisionIsMatched) { auto collisionspermcpjet = collisions.sliceBy(CollisionsPerMCPCollision, jet.mcCollisionId()); - if (collisionspermcpjet.size() >= 1 && jetderiveddatautilities::selectCollision(collisionspermcpjet.begin(), eventSelection)) { + if (collisionspermcpjet.size() >= 1 && jetderiveddatautilities::selectCollision(collisionspermcpjet.begin(), eventSelection)) { fillMCPHistograms(jet, jet.eventWeight()); } } else { @@ -924,7 +924,7 @@ struct JetFinderQATask { if (!isAcceptedJet(mcdjet)) { continue; } - if(rejectMbEvents && mcdjet.eventWeight()==1) { + if (rejectMbEvents && mcdjet.eventWeight() == 1) { continue; } fillMatchedHistograms::iterator, soa::Join>(mcdjet, mcdjet.eventWeight()); @@ -934,7 +934,7 @@ struct JetFinderQATask { void processMCCollisionsWeighted(aod::JetMcCollision const& collision) { - if(rejectMbEvents && collision.weight()==1) { + if (rejectMbEvents && collision.weight() == 1) { return; } registry.fill(HIST("h_collision_eventweight_part"), collision.weight()); @@ -1091,7 +1091,7 @@ struct JetFinderQATask { soa::Filtered> const& tracks) { float eventWeight = collision.mcCollision().weight(); - if(rejectMbEvents && eventWeight==1) { + if (rejectMbEvents && eventWeight == 1) { return; } registry.fill(HIST("h_collisions"), 0.5);