Hlt gen validation 1230pre4

DQMOffline/Trigger/src/FunctionDefs.cc

@@ -20,3 +20,12 @@ std::function<float(const reco::Photon&)> hltdqm::getUnaryFuncExtraFloat<reco::P
     varFunc = [](const reco::Photon& pho) -> float { return pho.hadTowOverEm(); };
   return varFunc;
 }
+
+//hack to put the function definitions in the non-plugin libary
+//so they can be used by other packages
+#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
+#include "DataFormats/EgammaCandidates/interface/Photon.h"
+#include "DataFormats/MuonReco/interface/Muon.h"
+auto testPho = hltdqm::getUnaryFuncFloat<reco::Photon>("et");
+auto testEle = hltdqm::getUnaryFuncFloat<reco::GsfElectron>("et");
+auto testMu = hltdqm::getUnaryFuncFloat<reco::Muon>("et");

Validation/HLTrigger/BuildFile.xml

@@ -0,0 +1,13 @@
+<use name="CommonTools/TriggerUtils"/>
+<use name="DQMOffline/Trigger"/>
+<use name="FWCore/Framework"/>
+<use name="FWCore/ParameterSet"/>
+<use name="FWCore/ServiceRegistry"/>
+<use name="DataFormats/HLTReco"/>
+<use name="HLTrigger/HLTcore"/>
+<use name="DQMServices/Core"/>
+<use name="CommonTools/Utils"/>
+<use name="root"/>
+<export>
+  <lib name="1"/>
+</export>

Validation/HLTrigger/interface/HLTGenValHist.h

@@ -0,0 +1,84 @@
+#ifndef Validation_HLTrigger_HLTGenValHist_h
+#define Validation_HLTrigger_HLTGenValHist_h
+
+//********************************************************************************
+//
+// Description:
+//   Histogram holder class for the GEN-level HLT validation
+//   Handling and filling of 1D and 2D histograms is done by this class.
+//
+//
+// Author : Finn Labe, UHH, Jul. 2022
+//          (Strongly inspired by Sam Harpers HLTDQMHist class)
+//
+//***********************************************************************************
+
+#include "DQMOffline/Trigger/interface/FunctionDefs.h"
+
+#include "DQMOffline/Trigger/interface/VarRangeCutColl.h"
+
+#include "FWCore/Framework/interface/Event.h"
+
+#include "Validation/HLTrigger/interface/HLTGenValObject.h"
+
+#include <TH1.h>
+#include <TH2.h>
+
+// base histogram class, with specific implementations following below
+class HLTGenValHist {
+public:
+  HLTGenValHist() = default;
+  virtual ~HLTGenValHist() = default;
+  virtual void fill(const HLTGenValObject& objType) = 0;
+};
+
+// specific implimentation of a HLTGenValHist for 1D histograms
+// it takes the histogram which it will fill
+// it takes the variable to plot (func) and its name (varName)
+// also, it takes additional cuts (rangeCuts) applied before filling
+// to fill the histogram, an object is passed in the Fill function
+class HLTGenValHist1D : public HLTGenValHist {
+public:
+  HLTGenValHist1D(TH1* hist,
+               std::string varName,
+               std::function<float(const HLTGenValObject&)> func,
+               VarRangeCutColl<HLTGenValObject> rangeCuts)
+      : var_(std::move(func)), varName_(std::move(varName)), rangeCuts_(std::move(rangeCuts)), hist_(hist) {}
+
+  void fill(const HLTGenValObject& obj) override {
+    if(rangeCuts_(obj)) hist_->Fill(var_(obj));
+  }
+
+private:
+  std::function<float(const HLTGenValObject&)> var_;
+  std::string varName_;
+  VarRangeCutColl<HLTGenValObject> rangeCuts_;
+  TH1* hist_;  //we do not own this
+};
+
+// specific implimentation of a HLTGenValHist for 2D histograms
+// it takes the histogram which it will fill
+// it takes the two variable to plot (func) and their name (varName)
+// to fill the histogram, two objects are passed in the Fill function
+class HLTGenValHist2D : public HLTGenValHist {
+public:
+  HLTGenValHist2D(TH2* hist,
+               std::string varNameX,
+               std::string varNameY,
+               std::function<float(const HLTGenValObject&)> funcX,
+               std::function<float(const HLTGenValObject&)> funcY)
+      : varX_(std::move(funcX)), varY_(std::move(funcY)), varNameX_(std::move(varNameX)), varNameY_(std::move(varNameY)), hist_(hist) {}
+
+  void fill(const HLTGenValObject& obj) override {
+    hist_->Fill(varX_(obj), varY_(obj));
+  }
+
+private:
+  std::function<float(const HLTGenValObject&)> varX_;
+  std::function<float(const HLTGenValObject&)> varY_;
+  std::string varNameX_;
+  std::string varNameY_;
+  TH2* hist_;  //we do not own this
+};
+
+#endif

Validation/HLTrigger/interface/HLTGenValHistCollFilter.h

@@ -0,0 +1,61 @@
+#ifndef Validation_HLTrigger_HLTGenValHistCollFilter_h
+#define Validation_HLTrigger_HLTGenValHistCollFilter_h
+
+//********************************************************************************
+//
+// Description:
+//   This class contains a collection of HLTGenvalHists used to measure the efficiency of a
+//   specified filter. It is resonsible for booking and filling the histograms of all vsVars
+//   histograms that are created for a specific filter.
+//
+// Author : Finn Labe, UHH, Jul. 2022
+//          (Strongly inspired by Sam Harpers HLTDQMFilterEffHists class)
+//
+//***********************************************************************************
+
+#include "FWCore/ParameterSet/interface/ParameterSet.h"
+#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
+
+#include "DQMServices/Core/interface/DQMStore.h"
+
+#include "DataFormats/HLTReco/interface/TriggerEvent.h"
+
+#include "Validation/HLTrigger/interface/HLTGenValHist.h"
+#include "DQMOffline/Trigger/interface/FunctionDefs.h"
+#include "DQMOffline/Trigger/interface/UtilFuncs.h"
+
+#include "Validation/HLTrigger/interface/HLTGenValObject.h"
+#include "Validation/HLTrigger/interface/HLTGenValPathSpecificSettingParser.h"
+
+#include "DataFormats/Math/interface/deltaR.h"
+
+#include <utility>
+
+// class containing a collection of HLTGenValHist for a specific filter
+// functions for initial booking of hists, and filling of hists for a single object are available
+class HLTGenValHistCollFilter {
+public:
+  typedef dqm::legacy::MonitorElement MonitorElement;
+  typedef dqm::legacy::DQMStore DQMStore;
+
+  explicit HLTGenValHistCollFilter(edm::ParameterSet filterCollConfig);
+
+  static edm::ParameterSetDescription makePSetDescription();
+
+  void bookHists(DQMStore::IBooker& iBooker, const std::vector<edm::ParameterSet>& histConfigs, const std::vector<edm::ParameterSet>& histConfigs2D);
+  void fillHists(const HLTGenValObject& obj, edm::Handle<trigger::TriggerEvent>& triggerEvent);
+
+private:
+  void book1D(DQMStore::IBooker& iBooker, const edm::ParameterSet& histConfig);
+  void book2D(DQMStore::IBooker& iBooker, const edm::ParameterSet& histConfig2D);
+
+  std::vector<std::unique_ptr<HLTGenValHist>> hists_; // the collection of histograms
+  std::string objType_;
+  std::string tag_;
+  std::string filter_;
+  std::string path_;
+  std::string hltProcessName_;
+  double dR2limit_;
+};
+
+#endif

Validation/HLTrigger/interface/HLTGenValHistCollPath.h

@@ -0,0 +1,59 @@
+#ifndef Validation_HLTrigger_HLTGenValHistCollPath_h
+#define Validation_HLTrigger_HLTGenValHistCollPath_h
+
+//********************************************************************************
+//
+// Description:
+//   This contains a collection of HLTGenValHistCollFilter used to measure the efficiencies of all
+//   filters in a specified path. The actual booking and filling happens in the respective HLTGenValHistCollFilters.
+//
+// Author : Finn Labe, UHH, Oct. 2021
+// (Heavily borrowed from Sam Harpers HLTDQMFilterEffHists)
+//
+//***********************************************************************************
+
+#include "FWCore/ParameterSet/interface/ParameterSet.h"
+#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
+
+#include "DQMServices/Core/interface/DQMStore.h"
+
+#include "DataFormats/HLTReco/interface/TriggerEvent.h"
+
+#include "DQMOffline/Trigger/interface/FunctionDefs.h"
+#include "DQMOffline/Trigger/interface/UtilFuncs.h"
+
+#include "Validation/HLTrigger/interface/HLTGenValHistCollFilter.h"
+#include "HLTrigger/HLTcore/interface/HLTConfigProvider.h"
+
+#include "Validation/HLTrigger/interface/HLTGenValObject.h"
+
+// class containing a collection of HLTGenValHistCollFilters for a specific path
+// at object creation time, the object type (used for systematically naming the histogram),
+// triggerPath, hltConfig and dR2limit (for deltaR matching) need to be specified
+// functions for initial booking of hists, and filling of hists for a single object, are available
+class HLTGenValHistCollPath {
+public:
+  typedef dqm::legacy::MonitorElement MonitorElement;
+  typedef dqm::legacy::DQMStore DQMStore;
+
+  explicit HLTGenValHistCollPath(edm::ParameterSet pathCollConfig, HLTConfigProvider& hltConfig);
+
+  static edm::ParameterSetDescription makePSetDescription();
+
+  void bookHists(DQMStore::IBooker& iBooker, std::vector<edm::ParameterSet>& histConfigs, std::vector<edm::ParameterSet>& histConfigs2D);
+  void fillHists(const HLTGenValObject& obj, edm::Handle<trigger::TriggerEvent>& triggerEvent);
+
+private:
+  std::string triggerPath_;
+  std::vector<HLTGenValHistCollFilter> collectionFilter_;
+  std::vector<std::string> filters_;
+  HLTConfigProvider hltConfig_;
+  bool doOnlyLastFilter_;
+
+  // we will add a string to the root file that is named after the path
+  // it will contain the filters of that path divided by semicola
+  std::string pathStringName_;
+  std::string pathString_;
+};
+
+#endif

Validation/HLTrigger/interface/HLTGenValObject.h

@@ -0,0 +1,51 @@
+#ifndef Validation_HLTrigger_HLTGenValObject_h
+#define Validation_HLTrigger_HLTGenValObject_h
+
+//********************************************************************************
+//
+// Description:
+//   This class is an object wrapper for the Generator level validation code.
+//   It handles the different type of objects the code needs to run on: GenParticles, GenJets and event-level energy sums
+//
+// Author : Finn Labe, UHH, Nov. 2021
+//
+//********************************************************************************
+
+
+#include "DataFormats/HepMCCandidate/interface/GenParticle.h"
+#include "DataFormats/JetReco/interface/GenJet.h"
+
+#include "DataFormats/Math/interface/LorentzVector.h"
+#include <vector>
+
+class HLTGenValObject {
+public:
+
+  // empty constructor
+  HLTGenValObject() {}
+
+  // constructor from GenParticle
+  HLTGenValObject(const reco::GenParticle &p)
+      : p4Polar_(p.p4()), p4Cartesian_(p.p4()) {}
+
+  // constructor from GenJet
+  HLTGenValObject(const reco::GenJet &p)
+      : p4Polar_(p.p4()), p4Cartesian_(p.p4()) {}
+
+  // constructor from LorentzVector (for energy sums)
+  HLTGenValObject(const reco::Candidate::PolarLorentzVector& p)
+      : p4Polar_(p), p4Cartesian_(p) {}
+
+  // object functions, for usage of HLTGenValObjects by other modules
+  double pt() const { return p4Polar_.pt(); }
+  double eta() const { return p4Polar_.eta(); }
+  double phi() const { return p4Polar_.phi(); }
+  double et() const { return (pt() <= 0) ? 0 : p4Cartesian_.Et(); }
+
+private:
+  // containing information in two "shapes"
+  math::PtEtaPhiMLorentzVector p4Polar_;
+  math::XYZTLorentzVector p4Cartesian_;
+};
+
+#endif

Validation/HLTrigger/interface/HLTGenValPathSpecificSettingParser.h

@@ -0,0 +1,40 @@
+#ifndef Validation_HLTrigger_HLTGenValPathSpecificSettingParser_h
+#define Validation_HLTrigger_HLTGenValSpecificCutParser_h
+
+//********************************************************************************
+//
+// Description:
+//   This class handles parsing of additional settings that can be set for each path in the generator-level validation module
+//   Mainly, these are cuts in addition to the ones specified in the module. Passing a pre-defined region is also possible
+//   The binning of a certain variable can be changed through this class, as well as setting a tag for all histograms of a path.
+//
+// Author : Finn Labe, UHH, Jul. 2022
+//
+//********************************************************************************
+
+#include "FWCore/ParameterSet/interface/ParameterSet.h"
+#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
+
+#include "DQMServices/Core/interface/DQMStore.h"
+
+#include <vector>
+
+class HLTGenValPathSpecificSettingParser {
+public:
+
+  // constructor
+  HLTGenValPathSpecificSettingParser(std::string pathSpecificSettings, std::vector<edm::ParameterSet> binnings, std::string vsVar);
+
+  std::vector<edm::ParameterSet> getPathSpecificCuts() {return pathSpecificCutsVector_;}
+  std::vector<double> getPathSpecificBins() {return pathSpecificBins_;}
+  bool havePathSpecificBins() { return (!pathSpecificBins_.empty()); }
+  std::string getTag() {return tag_;}
+
+private:
+  std::vector<edm::ParameterSet> pathSpecificCutsVector_;
+  std::vector<double> pathSpecificBins_;
+  std::string tag_ = "";
+
+};
+
+#endif

Validation/HLTrigger/plugins/BuildFile.xml

@@ -0,0 +1,15 @@
+<use name="DataFormats/TrackReco"/>
+<use name="FWCore/Framework"/>
+<use name="FWCore/ParameterSet"/>
+<use name="FWCore/PluginManager"/>
+<use name="CommonTools/UtilAlgos"/>
+<use name="PhysicsTools/UtilAlgos"/>
+<use name="DQMServices/Core"/>
+<use name="CommonTools/TriggerUtils"/>
+<use name="DQMOffline/Trigger"/>
+<use name="FWCore/ServiceRegistry"/>
+<use name="DataFormats/HLTReco"/>
+<use name="HLTrigger/HLTcore"/>
+<use name="CommonTools/Utils"/>
+<use name="Validation/HLTrigger"/>
+<flags EDM_PLUGIN="1"/>