Skip to content

EXOnanoAODv1 producers and setup #48502

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

Open
kerstinlovisa wants to merge 3 commits into
base: master
Choose a base branch
from
+4,732 −0
Open

EXOnanoAODv1 producers and setup #48502

Show file tree
Hide file tree
Changes from 1 commit
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
c51bf4e
Including EXOnanoAODv1 custom constent
Dec 8, 2025
3913626
Muon PR comments
Jan 30, 2026
81669af
Code format
Jan 30, 2026
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
12 changes: 12 additions & 0 deletions Configuration/PyReleaseValidation/python/relval_nano.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -306,6 +306,9 @@ def subnext(self):
steps['BTVNANO_mc15.0'] = merge([{'-s': 'NANO:@BTV,DQM:@nanoAODDQM', '-n': '1000'},
steps['NANO_mc15.0']])

steps['EXONANO_mc15.0'] = merge([{'-s': 'PAT,NANO:@EXO', '-n': '1000'},
steps['NANO_mc15.0']])

steps['lepTrackInfoNANO_mc15.0'] = merge([{'-s': 'NANO:@LepTrackInfo,DQM:@nanoAODDQM', '-n': '1000'},
steps['NANO_mc15.0']])

Expand All @@ -323,10 +326,14 @@ def subnext(self):

# ====== DATA ======
lumis_Run2025C = {392175: [[95, 542]]}
lumis_Run2025C_Muon0 = {392997: [[146, 202]]}

steps['JetMET1_Run2025C_MINIAOD_150X'] = {'INPUT': InputInfo(
location='STD', ls=lumis_Run2025C, dataSet='/JetMET1/Run2025C-PromptReco-v1/MINIAOD')}

steps['Muon0_Run2025C_AOD_150X'] = {'INPUT': InputInfo(
location='STD',ls=lumis_Run2025C_Muon0, dataSet='/Muon0/Run2025C-PromptReco-v1/AOD')}

steps['ScoutingPFRun3_Run2025C_HLTSCOUT_150X'] = {'INPUT': InputInfo(location='STD', ls=lumis_Run2025C,
dataSet='/ScoutingPFRun3/Run2025C-v1/HLTSCOUT')}

Expand All @@ -350,6 +357,9 @@ def subnext(self):
steps['BTVNANO_data15.0'] = merge([{'-s': 'NANO:@BTV,DQM:@nanoAODDQM', '-n': '1000'},
steps['NANO_data15.0']])

steps['EXONANO_data15.0'] = merge([{'-s': 'PAT,NANO:@EXO', '-n': '1000'},
steps['NANO_data15.0']])

steps['lepTrackInfoNANO_data15.0'] = merge([{'-s': 'NANO:@LepTrackInfo,DQM:@nanoAODDQM', '-n': '1000'},
steps['NANO_data15.0']])

Expand Down Expand Up @@ -495,6 +505,7 @@ def subnext(self):
workflows[_wfn()] = ['ScoutingNANOmc150X', ['TTbar_13p6_Summer24_MINIAOD_150X', 'scoutingNANO_mc15.0']]
workflows[_wfn()] = ['ScoutingNANOwithPromptmc150X', ['TTbar_13p6_Summer24_MINIAOD_150X', 'scoutingNANO_withPrompt_mc15.0']]
workflows[_wfn()] = ['BPHNANOmc150X', ['TTbar_13p6_Summer24_MINIAOD_150X', 'BPHNANO_mc15.0']]
workflows[_wfn()] = ['EXONANOmc150X', ['TTbar_13p6_Summer24_AOD_140X', 'EXONANO_mc15.0']]

# POG/PAG custom NANOs, data
_wfn.subnext()
Expand All @@ -507,6 +518,7 @@ def subnext(self):
workflows[_wfn()] = ['ScoutingNANOdata150Xrun3', ['ScoutingPFRun3_Run2025C_HLTSCOUT_150X', 'scoutingNANO_data15.0']]
workflows[_wfn()] = ['ScoutingNANOwithPromptdata150Xrun3', ['ScoutingPFMonitor_Run2025C_MINIAOD_150X', 'scoutingNANO_withPrompt_data15.0']] # noqa
workflows[_wfn()] = ['BPHNANOdata150Xrun3', ['JetMET1_Run2025C_MINIAOD_150X', 'BPHNANO_data15.0']]
workflows[_wfn()] = ['EXONANOdata150Xrun3', ['Muon0_Run2025C_AOD_150X', 'EXONANO_data15.0']]

# DPG custom NANOs, data
_wfn.subnext()
Expand Down
191 changes: 191 additions & 0 deletions PhysicsTools/NanoAOD/interface/DisTauTagScaling.h
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,191 @@
namespace Scaling {
Copy link
Contributor

@makortel makortel Dec 17, 2025

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Namespaces should start with lower case letter (see 2.7 in https://cms-sw.github.io/cms_coding_rules.html)

Also the include guard is missing (4.1 in https://cms-sw.github.io/cms_coding_rules.html)

struct PfCand {
inline static const std::vector<std::vector<float>> mean = {{0, 0},
{50.0, 50.0},
{0.0, 0.0},
{0.0, 0.0},
{0.09469, 0.09469},
{0, 0},
{0, 0},
{0, 0},
{5.0, 5.0},
{5.0, 5.0},
{6.442, 6.442},
{0, 0},
{-0.01635, -0.01635},
{0.02704, 0.02704},
{-0.01443, -0.01443},
{0.05551, 0.05551},
{11.16, 11.16},
{13.53, 13.53},
{0.5, 0.5},
{0.5, 0.5},
{1.3, 1.3},
{0.5, 0.5},
{0, 0},
{0, 0}};
inline static const std::vector<std::vector<float>> std = {{1, 1}, {50.0, 50.0},
{3.0, 3.0}, {3.141592653589793, 3.141592653589793},
{0.0651, 0.0651}, {1, 1},
{1, 1}, {1, 1},
{5.0, 5.0}, {5.0, 5.0},
{8.344, 8.344}, {1, 1},
{2.4, 2.4}, {0.08913, 0.08913},
{7.444, 7.444}, {0.5998, 0.5998},
{60.39, 60.39}, {6.44, 6.44},
{0.5, 0.5}, {0.5, 0.5},
{1.3, 1.3}, {0.5, 0.5},
{1, 1}, {1, 1}};
inline static const std::vector<std::vector<float>> lim_min = {
{-std::numeric_limits<float>::infinity(), -std::numeric_limits<float>::infinity()},
{-1.0, -1.0},
{-1.0, -1.0},
{-1.0, -1.0},
{-5, -5},
{-std::numeric_limits<float>::infinity(), -std::numeric_limits<float>::infinity()},
{-std::numeric_limits<float>::infinity(), -std::numeric_limits<float>::infinity()},
{-std::numeric_limits<float>::infinity(), -std::numeric_limits<float>::infinity()},
{-1.0, -1.0},
{-1.0, -1.0},
{-5, -5},
{-std::numeric_limits<float>::infinity(), -std::numeric_limits<float>::infinity()},
{-5, -5},
{-5, -5},
{-5, -5},
{-5, -5},
{-5, -5},
{-5, -5},
{-1.0, -1.0},
{-1.0, -1.0},
{-1.0, -1.0},
{-1.0, -1.0},
{-std::numeric_limits<float>::infinity(), -std::numeric_limits<float>::infinity()},
{-std::numeric_limits<float>::infinity(), -std::numeric_limits<float>::infinity()}};
inline static const std::vector<std::vector<float>> lim_max = {
{std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity()},
{1.0, 1.0},
{1.0, 1.0},
{1.0, 1.0},
{5, 5},
{std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity()},
{std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity()},
{std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity()},
{1.0, 1.0},
{1.0, 1.0},
{5, 5},
{std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity()},
{5, 5},
{5, 5},
{5, 5},
{5, 5},
{5, 5},
{5, 5},
{1.0, 1.0},
{1.0, 1.0},
{1.0, 1.0},
{1.0, 1.0},
{std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity()},
{std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity()}};
};
struct PfCandCategorical {
inline static const std::vector<std::vector<float>> mean = {{0, 0}, {0, 0}, {0, 0}};
inline static const std::vector<std::vector<float>> std = {{1, 1}, {1, 1}, {1, 1}};
inline static const std::vector<std::vector<float>> lim_min = {
{-std::numeric_limits<float>::infinity(), -std::numeric_limits<float>::infinity()},
{-std::numeric_limits<float>::infinity(), -std::numeric_limits<float>::infinity()},
{-std::numeric_limits<float>::infinity(), -std::numeric_limits<float>::infinity()}};
inline static const std::vector<std::vector<float>> lim_max = {
{std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity()},
{std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity()},
{std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity()}};
};
}; // namespace Scaling

namespace Setup {
Copy link
Contributor

@makortel makortel Dec 17, 2025

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Setup is a very common name (that should start with lower case), so there is a risk for name collisions.

const inline long int output_classes = 2;
const inline size_t n_Global = 5;
const inline size_t n_PfCand = 24;
const inline size_t nSeq_PfCand = 50;
const inline size_t n_PfCandCategorical = 3;
const inline size_t nSeq_PfCandCategorical = 50;
const inline std::vector<std::string> CellObjectTypes{"PfCand", "PfCandCategorical"};
}; // namespace Setup

enum class Global_Features { jet_pt = 0, jet_eta = 1, Lxy = 2, Lz = 3, Lrel = 4 };

enum class PfCand_Features {
pfCand_valid = 0,
pfCand_pt = 1,
pfCand_eta = 2,
pfCand_phi = 3,
pfCand_mass = 4,
pfCand_charge = 5,
pfCand_puppiWeight = 6,
pfCand_puppiWeightNoLep = 7,
pfCand_lostInnerHits = 8,
pfCand_nPixelHits = 9,
pfCand_nHits = 10,
pfCand_hasTrackDetails = 11,
pfCand_dxy = 12,
pfCand_dxy_error = 13,
pfCand_dz = 14,
pfCand_dz_error = 15,
pfCand_track_chi2 = 16,
pfCand_track_ndof = 17,
pfCand_caloFraction = 18,
pfCand_hcalFraction = 19,
pfCand_rawCaloFraction = 20,
pfCand_rawHcalFraction = 21,
pfCand_deta = 22,
pfCand_dphi = 23
};

enum class PfCandCategorical_Features { pfCand_particleType = 0, pfCand_pvAssociationQuality = 1, pfCand_fromPV = 2 };

enum class CellObjectType { PfCand, PfCandCategorical };

template <typename T>
struct FeaturesHelper;

template <>
struct FeaturesHelper<PfCand_Features> {
static constexpr CellObjectType object_type = CellObjectType::PfCand;
static constexpr size_t size = 24;
static constexpr size_t length = 50;
using scaler_type = Scaling::PfCand;
};

template <>
struct FeaturesHelper<PfCandCategorical_Features> {
static constexpr CellObjectType object_type = CellObjectType::PfCandCategorical;
static constexpr size_t size = 3;
static constexpr size_t length = 50;
using scaler_type = Scaling::PfCandCategorical;
};

using FeatureTuple = std::tuple<PfCand_Features, PfCandCategorical_Features>;

enum class PFParticleType {
Undefined = 0, // undefined
h = 1, // charged hadron
e = 2, // electron
mu = 3, // muon
gamma = 4, // photon
h0 = 5, // neutral hadron
h_HF = 6, // HF tower identified as a hadron
egamma_HF = 7, // HF tower identified as an EM particle
};

inline PFParticleType TranslatePdgIdToPFParticleType(int pdgId) {
static const std::map<int, PFParticleType> type_map = {
{11, PFParticleType::e},
{13, PFParticleType::mu},
{22, PFParticleType::gamma},
{211, PFParticleType::h},
{130, PFParticleType::h0},
{1, PFParticleType::h_HF},
{2, PFParticleType::egamma_HF},
};
auto iter = type_map.find(std::abs(pdgId));
return iter == type_map.end() ? PFParticleType::Undefined : iter->second;
}
10 changes: 10 additions & 0 deletions PhysicsTools/NanoAOD/plugins/BuildFile.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -22,8 +22,17 @@
<use name="FWCore/ParameterSet"/>
<use name="FWCore/ServiceRegistry"/>
<use name="FWCore/Utilities"/>
<use name="FWCore/PluginManager"/>
<use name="Geometry/Records"/>
<use name="Geometry/TrackerGeometryBuilder"/>
<use name="Geometry/RPCGeometry"/>
<use name="Geometry/DTGeometry"/>
<use name="Geometry/CSCGeometry"/>
<use name="IOPool/Provenance"/>
<use name="MagneticField/Records"/>
<use name="PhysicsTools/PatAlgos"/>
<use name="PhysicsTools/RecoUtils"/>
<use name="RecoVertex/KalmanVertexFit"/>
<use name="RecoVertex/VertexPrimitives"/>
<use name="RecoVertex/VertexTools"/>
<library file="*.cc" name="PhysicsToolsNanoAODPlugins">
Expand All @@ -34,6 +43,7 @@
<use name="PhysicsTools/JetMCUtils"/>
<use name="RecoBTag/FeatureTools"/>
<use name="RecoJets/JetAlgorithms"/>
<use name="TrackingTools/IPTools"/>
<use name="TrackingTools/Records"/>
<use name="roottmva"/>
<flags EDM_PLUGIN="1"/>
Expand Down
Loading