Skip to content

Please consider the following formatting changes to #9484 #7

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

Closed
alibuild wants to merge 2 commits into from
Closed

Please consider the following formatting changes to #9484 #7

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
bdc015f
add MB gap event rejection for tracks and collisions
jaimenorman Jan 14, 2025
47a44f4
Please consider the following formatting changes
alibuild Jan 23, 2025
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
13 changes: 11 additions & 2 deletions PWGJE/Tasks/jetFinderFullQA.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -74,6 +74,7 @@ struct JetFinderFullQATask {
Configurable<float> pTHatMaxMCD{"pTHatMaxMCD", 999.0, "maximum fraction of hard scattering for jet acceptance in detector MC"};
Configurable<float> pTHatMaxMCP{"pTHatMaxMCP", 999.0, "maximum fraction of hard scattering for jet acceptance in particle MC"};
Configurable<float> pTHatExponent{"pTHatExponent", 6.0, "exponent of the event weight for the calculation of pTHat"};
Configurable<bool> skipMBGapEvents{"skipMBGapEvents", false, "flag to choose to reject min. bias gap events; jet-level rejection applied at the jet finder level, here rejection is applied for collision and track process functions"};

std::vector<bool> filledJetR;
std::vector<double> jetRadiiValues;
Expand Down Expand Up @@ -566,7 +567,9 @@ struct JetFinderFullQATask {

void processMCCollisionsWeighted(aod::JetMcCollision const& collision)
{

if (skipMBGapEvents && collision.subGeneratorId() == jetderiveddatautilities::JCollisionSubGeneratorId::mbGap) {
return;
}
registry.fill(HIST("h_collision_eventweight_part"), collision.weight());
}
PROCESS_SWITCH(JetFinderFullQATask, processMCCollisionsWeighted, "collision QA for weighted events", false);
Expand All @@ -575,6 +578,9 @@ struct JetFinderFullQATask {
soa::Filtered<aod::JetTracks> const& tracks,
soa::Filtered<aod::JetClusters> const& clusters)
{
if (skipMBGapEvents && collision.subGeneratorId() == jetderiveddatautilities::JCollisionSubGeneratorId::mbGap) {
return;
}
registry.fill(HIST("h_collisions"), 0.5);
registry.fill(HIST("h_centrality_collisions"), collision.centrality(), 0.5);
if (!jetderiveddatautilities::eventEMCAL(collision)) {
Expand All @@ -586,12 +592,15 @@ struct JetFinderFullQATask {
}
PROCESS_SWITCH(JetFinderFullQATask, processTracks, "QA for charged tracks", false);

void processTracksWeighted(soa::Join<aod::JCollisions, aod::JMcCollisionLbs>::iterator const& collision,
void processTracksWeighted(soa::Join<aod::JetCollisions, aod::JMcCollisionLbs>::iterator const& collision,
aod::JMcCollisions const&,
soa::Filtered<aod::JetTracks> const& tracks,
soa::Filtered<aod::JetClusters> const& clusters)
{
float eventWeight = collision.mcCollision().weight();
if (skipMBGapEvents && collision.subGeneratorId() == jetderiveddatautilities::JCollisionSubGeneratorId::mbGap) {
return;
}
registry.fill(HIST("h_collisions"), 0.5);
registry.fill(HIST("h_collisions_weighted"), 0.5, eventWeight);
if (!jetderiveddatautilities::eventEMCAL(collision)) {
Expand Down
22 changes: 22 additions & 0 deletions PWGJE/Tasks/jetFinderHFQA.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -69,6 +69,7 @@ struct JetFinderHFQATask {
Configurable<float> pTHatMaxMCP{"pTHatMaxMCP", 999.0, "maximum fraction of hard scattering for jet acceptance in particle MC"};
Configurable<float> pTHatExponent{"pTHatExponent", 6.0, "exponent of the event weight for the calculation of pTHat"};
Configurable<float> randomConeR{"randomConeR", 0.4, "size of random Cone for estimating background fluctuations"};
Configurable<bool> skipMBGapEvents{"skipMBGapEvents", false, "flag to choose to reject min. bias gap events; jet-level rejection applied at the jet finder level, here rejection is applied for collision and track process functions"};

HfHelper hfHelper;
std::vector<bool> filledJetR_Both;
Expand Down Expand Up @@ -1173,6 +1174,9 @@ struct JetFinderHFQATask {

void processMCCollisionsWeighted(aod::JetMcCollision const& collision)
{
if (skipMBGapEvents && collision.subGeneratorId() == jetderiveddatautilities::JCollisionSubGeneratorId::mbGap) {
return;
}
registry.fill(HIST("h_collision_eventweight_part"), collision.weight());
}
PROCESS_SWITCH(JetFinderHFQATask, processMCCollisionsWeighted, "collision QA for weighted events", false);
Expand Down Expand Up @@ -1463,6 +1467,9 @@ struct JetFinderHFQATask {
void processTracks(soa::Filtered<aod::JetCollisions>::iterator const& collision,
soa::Filtered<aod::JetTracks> const& tracks)
{
if (skipMBGapEvents && collision.subGeneratorId() == jetderiveddatautilities::JCollisionSubGeneratorId::mbGap) {
return;
}
registry.fill(HIST("h_collisions"), 0.5);
registry.fill(HIST("h2_centrality_collisions"), collision.centrality(), 0.5);
if (!jetderiveddatautilities::selectCollision(collision, eventSelection)) {
Expand All @@ -1479,6 +1486,9 @@ struct JetFinderHFQATask {
soa::Filtered<aod::JetTracks> const& tracks)
{
float eventWeight = collision.mcCollision().weight();
if (skipMBGapEvents && collision.subGeneratorId() == jetderiveddatautilities::JCollisionSubGeneratorId::mbGap) {
return;
}
registry.fill(HIST("h_collisions"), 0.5);
registry.fill(HIST("h_collisions_weighted"), 0.5, eventWeight);
if (!jetderiveddatautilities::selectCollision(collision, eventSelection)) {
Expand All @@ -1494,6 +1504,9 @@ struct JetFinderHFQATask {
CandidateTableData const& candidates,
soa::Filtered<JetTracksDataSub> const& tracks)
{
if (skipMBGapEvents && collision.subGeneratorId() == jetderiveddatautilities::JCollisionSubGeneratorId::mbGap) {
return;
}
if (!jetderiveddatautilities::selectCollision(collision, eventSelection)) {
return;
}
Expand All @@ -1511,6 +1524,9 @@ struct JetFinderHFQATask {

void processRho(aod::JetCollision const& collision, soa::Join<CandidateTableData, BkgRhoTable> const& candidates, soa::Filtered<aod::JetTracks> const& tracks)
{
if (skipMBGapEvents && collision.subGeneratorId() == jetderiveddatautilities::JCollisionSubGeneratorId::mbGap) {
return;
}
if (!jetderiveddatautilities::selectCollision(collision, eventSelection)) {
return;
}
Expand Down Expand Up @@ -1539,13 +1555,19 @@ struct JetFinderHFQATask {

void processRandomConeMCD(soa::Filtered<aod::JetCollisions>::iterator const& collision, JetTableMCDJoined const& jets, soa::Join<CandidateTableMCD, BkgRhoTable> const& candidates, soa::Filtered<aod::JetTracks> const& tracks)
{
if (skipMBGapEvents && collision.subGeneratorId() == jetderiveddatautilities::JCollisionSubGeneratorId::mbGap) {
return;
}
randomCone(collision, jets, candidates, tracks);
}
PROCESS_SWITCH(JetFinderHFQATask, processRandomConeMCD, "QA for random cone estimation of background fluctuations in mcd", false);

void processCandidates(soa::Filtered<aod::JetCollisions>::iterator const& collision, CandidateTableData const& candidates)
{

if (skipMBGapEvents && collision.subGeneratorId() == jetderiveddatautilities::JCollisionSubGeneratorId::mbGap) {
return;
}
for (auto const& candidate : candidates) {
registry.fill(HIST("h_candidate_invmass"), jetcandidateutilities::getCandidateInvariantMass(candidate));
registry.fill(HIST("h_candidate_pt"), candidate.pt());
Expand Down
19 changes: 19 additions & 0 deletions PWGJE/Tasks/jetFinderQA.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -72,6 +72,7 @@ struct JetFinderQATask {
Configurable<bool> checkMcCollisionIsMatched{"checkMcCollisionIsMatched", false, "0: count whole MCcollisions, 1: select MCcollisions which only have their correspond collisions"};
Configurable<int> 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<int> 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<bool> skipMBGapEvents{"skipMBGapEvents", false, "flag to choose to reject min. bias gap events; jet-level rejection applied at the jet finder level, here rejection is applied for collision and track process functions"};

std::vector<bool> filledJetR_Both;
std::vector<bool> filledJetR_Low;
Expand Down Expand Up @@ -924,6 +925,9 @@ struct JetFinderQATask {

void processMCCollisionsWeighted(aod::JetMcCollision const& collision)
{
if (skipMBGapEvents && collision.subGeneratorId() == jetderiveddatautilities::JCollisionSubGeneratorId::mbGap) {
return;
}
registry.fill(HIST("h_collision_eventweight_part"), collision.weight());
}
PROCESS_SWITCH(JetFinderQATask, processMCCollisionsWeighted, "collision QA for weighted events", false);
Expand Down Expand Up @@ -1057,6 +1061,9 @@ struct JetFinderQATask {
void processTracks(soa::Filtered<aod::JetCollisions>::iterator const& collision,
soa::Filtered<soa::Join<aod::JetTracks, aod::JTrackExtras>> const& tracks)
{
if (skipMBGapEvents && collision.subGeneratorId() == jetderiveddatautilities::JCollisionSubGeneratorId::mbGap) {
return;
}
registry.fill(HIST("h_collisions"), 0.5);
registry.fill(HIST("h2_centrality_collisions"), collision.centrality(), 0.5);
if (!jetderiveddatautilities::selectCollision(collision, eventSelection)) {
Expand All @@ -1078,6 +1085,9 @@ struct JetFinderQATask {
soa::Filtered<soa::Join<aod::JetTracks, aod::JTrackExtras>> const& tracks)
{
float eventWeight = collision.mcCollision().weight();
if (skipMBGapEvents && collision.subGeneratorId() == jetderiveddatautilities::JCollisionSubGeneratorId::mbGap) {
return;
}
registry.fill(HIST("h_collisions"), 0.5);
registry.fill(HIST("h_collisions_weighted"), 0.5, eventWeight);
if (!jetderiveddatautilities::selectCollision(collision, eventSelection)) {
Expand All @@ -1097,6 +1107,9 @@ struct JetFinderQATask {
void processTracksSub(soa::Filtered<aod::JetCollisions>::iterator const& collision,
soa::Filtered<aod::JetTracksSub> const& tracks)
{
if (skipMBGapEvents && collision.subGeneratorId() == jetderiveddatautilities::JCollisionSubGeneratorId::mbGap) {
return;
}
if (!jetderiveddatautilities::selectCollision(collision, eventSelection)) {
return;
}
Expand All @@ -1113,6 +1126,9 @@ struct JetFinderQATask {

void processRho(soa::Filtered<soa::Join<aod::JetCollisions, aod::BkgChargedRhos>>::iterator const& collision, soa::Filtered<aod::JetTracks> const& tracks)
{
if (skipMBGapEvents && collision.subGeneratorId() == jetderiveddatautilities::JCollisionSubGeneratorId::mbGap) {
return;
}
if (!jetderiveddatautilities::selectCollision(collision, eventSelection)) {
return;
}
Expand Down Expand Up @@ -1141,6 +1157,9 @@ struct JetFinderQATask {

void processRandomConeMCD(soa::Filtered<soa::Join<aod::JetCollisions, aod::BkgChargedRhos>>::iterator const& collision, soa::Join<aod::ChargedMCDetectorLevelJets, aod::ChargedMCDetectorLevelJetConstituents> const& jets, soa::Filtered<aod::JetTracks> const& tracks)
{
if (skipMBGapEvents && collision.subGeneratorId() == jetderiveddatautilities::JCollisionSubGeneratorId::mbGap) {
return;
}
randomCone(collision, jets, tracks);
}
PROCESS_SWITCH(JetFinderQATask, processRandomConeMCD, "QA for random cone estimation of background fluctuations in mcd", false);
Expand Down
Loading