Skip to content

[PWGLF] Added zorro for processing of skimmed data and removed magic numbers #10735

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
alcaliva merged 4 commits into from
Apr 3, 2025
Merged

[PWGLF] Added zorro for processing of skimmed data and removed magic numbers #10735

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
1ca651c
[PWGLF] Added zorro for processing of skimmed data
alcaliva Apr 2, 2025
e53a291
added EventFilteringUtils to public link libraries
alcaliva Apr 2, 2025
8318fb3
moved zorroSummary to init func
alcaliva Apr 2, 2025
c4129c7
removed whitespace
alcaliva Apr 2, 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
2 changes: 1 addition & 1 deletion PWGLF/Tasks/Strangeness/CMakeLists.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,3 +1,3 @@
# Copyright 2019-2020 CERN and copyright holders of ALICE O2.
# See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
# All rights not expressly granted are reserved.
Expand All @@ -9,7 +9,7 @@
# granted to it by virtue of its status as an Intergovernmental Organization
# or submit itself to any jurisdiction.

o2physics_add_dpl_workflow(hyperon-reco-test

Check warning on line 12 in PWGLF/Tasks/Strangeness/CMakeLists.txt

View workflow job for this annotation

GitHub Actions / O2 linter

[name/o2-workflow] 

Workflow name hyperon-reco-test does not match its file name hyperon-reco-test.cxx. (Matches hyperonRecoTest.cxx.)
SOURCES hyperon-reco-test.cxx
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
COMPONENT_NAME Analysis)
Expand All @@ -19,7 +19,7 @@
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore O2Physics::MLCore O2Physics::AnalysisCCDB
COMPONENT_NAME Analysis)

o2physics_add_dpl_workflow(lambdakzeroanalysis-mc

Check warning on line 22 in PWGLF/Tasks/Strangeness/CMakeLists.txt

View workflow job for this annotation

GitHub Actions / O2 linter

[name/o2-workflow] 

Workflow name lambdakzeroanalysis-mc does not match its file name lambdakzeroanalysisMC.cxx. (Matches lambdakzeroanalysisMc.cxx.)
SOURCES lambdakzeroanalysisMC.cxx
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
COMPONENT_NAME Analysis)
Expand All @@ -29,12 +29,12 @@
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
COMPONENT_NAME Analysis)

o2physics_add_dpl_workflow(sigmaminus-task

Check warning on line 32 in PWGLF/Tasks/Strangeness/CMakeLists.txt

View workflow job for this annotation

GitHub Actions / O2 linter

[name/o2-workflow] 

Workflow name sigmaminus-task does not match its file name sigmaminustask.cxx. (Matches sigmaminusTask.cxx.)
SOURCES sigmaminustask.cxx
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
COMPONENT_NAME Analysis)

o2physics_add_dpl_workflow(cascadeanalysismc

Check warning on line 37 in PWGLF/Tasks/Strangeness/CMakeLists.txt

View workflow job for this annotation

GitHub Actions / O2 linter

[name/o2-workflow] 

Workflow name cascadeanalysismc does not match its file name cascadeanalysisMC.cxx. (Matches cascadeanalysismc.cxx.)
SOURCES cascadeanalysisMC.cxx
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
COMPONENT_NAME Analysis)
Expand All @@ -54,12 +54,12 @@
PUBLIC_LINK_LIBRARIES O2::Framework O2::ReconstructionDataFormats O2Physics::AnalysisCore O2::DetectorsBase O2Physics::EventFilteringUtils
COMPONENT_NAME Analysis)

o2physics_add_dpl_workflow(k0mixedevents

Check warning on line 57 in PWGLF/Tasks/Strangeness/CMakeLists.txt

View workflow job for this annotation

GitHub Actions / O2 linter

[name/o2-workflow] 

Workflow name k0mixedevents does not match its file name k0_mixed_events.cxx. (Matches k0mixedevents.cxx.)
SOURCES k0_mixed_events.cxx
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
COMPONENT_NAME Analysis)

o2physics_add_dpl_workflow(vzero-cascade-absorption

Check warning on line 62 in PWGLF/Tasks/Strangeness/CMakeLists.txt

View workflow job for this annotation

GitHub Actions / O2 linter

[name/o2-workflow] 

Workflow name vzero-cascade-absorption does not match its file name vzero_cascade_absorption.cxx. (Matches vzeroCascadeAbsorption.cxx.)
SOURCES vzero_cascade_absorption.cxx
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
COMPONENT_NAME Analysis)
Expand All @@ -74,7 +74,7 @@
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
COMPONENT_NAME Analysis)

o2physics_add_dpl_workflow(hstrangecorrelation

Check warning on line 77 in PWGLF/Tasks/Strangeness/CMakeLists.txt

View workflow job for this annotation

GitHub Actions / O2 linter

[name/o2-workflow] 

Workflow name hstrangecorrelation does not match its file name hStrangeCorrelation.cxx. (Matches hstrangecorrelation.cxx.)
SOURCES hStrangeCorrelation.cxx
PUBLIC_LINK_LIBRARIES O2::Framework O2::DetectorsBase O2Physics::AnalysisCore O2Physics::EventFilteringUtils
COMPONENT_NAME Analysis)
Expand Down Expand Up @@ -107,7 +107,7 @@
if(FastJet_FOUND)
o2physics_add_dpl_workflow(strangeness-in-jets
SOURCES strangenessInJets.cxx
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore O2Physics::PWGJECore FastJet::FastJet FastJet::Contrib
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore O2Physics::PWGJECore FastJet::FastJet FastJet::Contrib O2Physics::EventFilteringUtils
COMPONENT_NAME Analysis)
endif()

Expand Down Expand Up @@ -136,7 +136,7 @@
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore
COMPONENT_NAME Analysis)

o2physics_add_dpl_workflow(lambdajetpolarization

Check warning on line 139 in PWGLF/Tasks/Strangeness/CMakeLists.txt

View workflow job for this annotation

GitHub Actions / O2 linter

[name/o2-workflow] 

Workflow name lambdajetpolarization does not match its file name lambdaJetpolarization.cxx. (Matches lambdajetpolarization.cxx.)
SOURCES lambdaJetpolarization.cxx
PUBLIC_LINK_LIBRARIES O2::Framework O2Physics::AnalysisCore
COMPONENT_NAME Analysis)
78 changes: 59 additions & 19 deletions PWGLF/Tasks/Strangeness/strangenessInJets.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -52,6 +52,9 @@
#include "PWGJE/DataModel/JetReducedData.h"
#include "PWGJE/DataModel/Jet.h"

#include "EventFiltering/Zorro.h"
#include "EventFiltering/ZorroSummary.h"

using namespace std;
using namespace o2;
using namespace o2::soa;
Expand All @@ -71,6 +74,12 @@ using MCTracks = soa::Join<StrHadronDaughterTracks, aod::McTrackLabels>;

struct StrangenessInJets {

Service<o2::ccdb::BasicCCDBManager> ccdb;
o2::ccdb::CcdbApi ccdbApi;

Zorro zorro;
OutputObj<ZorroSummary> zorroSummary{"zorroSummary"};

HistogramRegistry registryData{"registryData", {}, OutputObjHandlingPolicy::AnalysisObject, true, true};
HistogramRegistry registryMC{"registryMC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true};
HistogramRegistry registryQC{"registryQC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true};
Expand All @@ -81,6 +90,7 @@ struct StrangenessInJets {
Configurable<double> rJet{"rJet", 0.3, "Jet resolution parameter R"};
Configurable<double> zVtx{"zVtx", 10.0, "Maximum zVertex"};
Configurable<double> deltaEtaEdge{"deltaEtaEdge", 0.05, "eta gap from the edge"};
Configurable<bool> cfgSkimmedProcessing{"cfgSkimmedProcessing", false, "Skimmed dataset processing"};

// Axis parameters
struct : ConfigurableGroup {
Expand Down Expand Up @@ -167,9 +177,6 @@ struct StrangenessInJets {
TH1F* weightsAntiXiInJet;
TH1F* weightsAntiXiInUe;

Service<o2::ccdb::BasicCCDBManager> ccdb;
o2::ccdb::CcdbApi ccdbApi;

// List of Particles
enum Option { KZeroLambda,
CascadePart,
Expand All @@ -180,8 +187,19 @@ struct StrangenessInJets {
// Jet background subtraction
JetBkgSubUtils backgroundSub;

void initCCDB(aod::BCsWithTimestamps::iterator const& bc)
{
if (cfgSkimmedProcessing) {
zorro.initCCDB(ccdb.service, bc.runNumber(), bc.timestamp(), "fOmega");
}
}

void init(InitContext const&)
{
if (cfgSkimmedProcessing) {
zorroSummary.setObject(zorro.getZorroSummary());
}

ccdb->setURL(urlToCcdb.value);
ccdb->setCaching(true);
ccdb->setLocalObjectValidityChecking();
Expand Down Expand Up @@ -501,27 +519,38 @@ struct StrangenessInJets {
template <typename JetTrack>
bool passedTrackSelectionForJetReconstruction(const JetTrack& track)
{
const int minTpcCr = 70;
const double minCrFindable = 0.8;
const double maxChi2Tpc = 4.0;
const double maxChi2Its = 36.0;
const double maxPseudorapidity = 0.8;
const double minPtTrack = 0.1;
const double dcaxyMaxTrackPar0 = 0.0105;
const double dcaxyMaxTrackPar1 = 0.035;
const double dcaxyMaxTrackPar2 = 1.1;
const double dcazMaxTrack = 2.0;

if (!track.hasITS())
return false;
if ((!hasITSHit(track, 1)) && (!hasITSHit(track, 2)) && (!hasITSHit(track, 3)))
return false;
if (!track.hasTPC())
return false;
if (track.tpcNClsCrossedRows() < 70)
if (track.tpcNClsCrossedRows() < minTpcCr)
return false;
if ((static_cast<double>(track.tpcNClsCrossedRows()) / static_cast<double>(track.tpcNClsFindable())) < 0.8)
if ((static_cast<double>(track.tpcNClsCrossedRows()) / static_cast<double>(track.tpcNClsFindable())) < minCrFindable)
return false;
if (track.tpcChi2NCl() > 4)
if (track.tpcChi2NCl() > maxChi2Tpc)
return false;
if (track.itsChi2NCl() > 36)
if (track.itsChi2NCl() > maxChi2Its)
return false;
if (track.eta() < -0.8 || track.eta() > 0.8)
if (track.eta() < -maxPseudorapidity || track.eta() > maxPseudorapidity)
return false;
if (track.pt() < 0.1)
if (track.pt() < minPtTrack)
return false;
if (std::fabs(track.dcaXY()) > (0.0105 + 0.035 / std::pow(track.pt(), 1.1)))
if (std::fabs(track.dcaXY()) > (dcaxyMaxTrackPar0 + dcaxyMaxTrackPar1 / std::pow(track.pt(), dcaxyMaxTrackPar2)))
return false;
if (std::fabs(track.dcaZ()) > 2.0)
if (std::fabs(track.dcaZ()) > dcazMaxTrack)
return false;
return true;
}
Expand Down Expand Up @@ -932,9 +961,11 @@ struct StrangenessInJets {
template <typename pionTrack>
bool isHighPurityPion(const pionTrack& track, const float nsigmaTPC, const float nsigmaTOF)
{
if (track.p() < 0.6 && std::fabs(nsigmaTPC) < 2.0)
double nsigmaPi = 2.0;
double pThreshold = 0.6;
if (track.p() < pThreshold && std::fabs(nsigmaTPC) < nsigmaPi)
return true;
if (track.p() > 0.6 && std::fabs(nsigmaTPC) < 2.0 && std::fabs(nsigmaTOF) < 2.0)
if (track.p() > pThreshold && std::fabs(nsigmaTPC) < nsigmaPi && std::fabs(nsigmaTOF) < nsigmaPi)
return true;
return false;
}
Expand Down Expand Up @@ -1016,6 +1047,13 @@ struct StrangenessInJets {
// event counter: before event selection
registryData.fill(HIST("number_of_events_data"), 0.5);

auto bc = collision.template bc_as<aod::BCsWithTimestamps>();
initCCDB(bc);

if (cfgSkimmedProcessing && !zorro.isSelected(collision.template bc_as<aod::BCsWithTimestamps>().globalBC())) {
return;
}

// event selection
if (!collision.sel8() || std::fabs(collision.posZ()) > zVtx)
return;
Expand Down Expand Up @@ -1421,7 +1459,7 @@ struct StrangenessInJets {
continue;

registryMC.fill(HIST("number_of_events_mc"), 1.5);
if (std::fabs(collision.posZ()) > 10.0)
if (std::fabs(collision.posZ()) > zVtx)
continue;

registryMC.fill(HIST("number_of_events_mc"), 2.5);
Expand Down Expand Up @@ -1886,7 +1924,7 @@ struct StrangenessInJets {
for (const auto& mccollision : mcCollisions) {

// Selection on z_{vertex}
if (std::fabs(mccollision.posZ()) > 10.0)
if (std::fabs(mccollision.posZ()) > zVtx)
continue;

// MC Particles per Collision
Expand All @@ -1898,9 +1936,10 @@ struct StrangenessInJets {

if (!particle.isPhysicalPrimary())
continue;
if (std::fabs(particle.eta()) > 0.8)
if (std::fabs(particle.eta()) > etaMax)
continue;
if (particle.pt() < 0.1)
double ptMinPart = 0.1;
if (particle.pt() < ptMinPart)
continue;

// 4-momentum representation of a particle
Expand Down Expand Up @@ -1956,9 +1995,10 @@ struct StrangenessInJets {

if (!particle.isPhysicalPrimary())
continue;
if (std::fabs(particle.eta()) > 0.8)
if (std::fabs(particle.eta()) > etaMax)
continue;
if (particle.pt() < 0.1)
double ptMinPart = 0.1;
if (particle.pt() < ptMinPart)
continue;

TVector3 particleDir(particle.px(), particle.py(), particle.pz());
Expand Down
Loading