Egamma/plugins/HLTEgammaGenericQuadraticFilter.h

0001 #ifndef HLTEgammaGenericQuadraticFilter_h
0002 #define HLTEgammaGenericQuadraticFilter_h
0003
0004 /** \class HLTEgammaGenericQuadraticFilter
0005  *
0006  *  \author Roberto Covarelli (CERN)
0007  *  modified by Chris Tully (Princeton)
0008  */
0009
0010 #include "HLTrigger/HLTcore/interface/HLTFilter.h"
0011
0012 #include "DataFormats/HLTReco/interface/TriggerFilterObjectWithRefs.h"
0013 #include "DataFormats/RecoCandidate/interface/RecoEcalCandidateIsolation.h"
0014
0015 namespace edm {
0016   class ConfigurationDescriptions;
0017 }
0018
0019 //
0020 // class declaration
0021 //
0022
0023 class HLTEgammaGenericQuadraticFilter : public HLTFilter {
0024 public:
0025   explicit HLTEgammaGenericQuadraticFilter(const edm::ParameterSet&);
0026   ~HLTEgammaGenericQuadraticFilter() override;
0027   bool hltFilter(edm::Event&,
0028                  const edm::EventSetup&,
0029                  trigger::TriggerFilterObjectWithRefs& filterproduct) const override;
0030   static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
0031
0032 private:
0033   edm::InputTag candTag_;  // input tag identifying product that contains filtered candidates
0034   edm::InputTag varTag_;   // input tag identifying product that contains the variable map
0035   edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> candToken_;
0036   edm::EDGetTokenT<reco::RecoEcalCandidateIsolationMap> varToken_;
0037
0038   bool lessThan_;  // the cut is "<" or ">" ?
0039   bool useEt_;     // use E or Et in relative isolation cuts
0040   /*  Barrel quadratic threshold function:
0041       vali (<= or >=) thrRegularEB_ + (E or Et)*thrOverEEB_ + (E or Et)*(E or Et)*thrOverE2EB_
0042     Endcap quadratic threshold function:
0043       vali (<= or >=) thrRegularEE_ + (E or Et)*thrOverEEE_ + (E or Et)*(E or Et)*thrOverE2EE_
0044       */
0045   std::vector<double> energyLowEdges_;  // lower bin edges for energy-dependent cuts
0046   std::vector<double> thrRegularEB_;    // threshold value for zeroth order term - ECAL barrel
0047   std::vector<double> thrRegularEE_;    // threshold value for zeroth order term - ECAL endcap
0048   std::vector<double> thrOverEEB_;      // coefficient for first order term - ECAL barrel
0049   std::vector<double> thrOverEEE_;      // coefficient for first order term - ECAL endcap
0050   std::vector<double> thrOverE2EB_;     // coefficient for second order term - ECAL barrel
0051   std::vector<double> thrOverE2EE_;     // coefficient for second order term - ECAL endcap
0052   int ncandcut_;                        // number of candidates required
0053
0054   edm::InputTag l1EGTag_;
0055
0056   edm::InputTag rhoTag_;  // input tag identifying rho producer
0057   edm::EDGetTokenT<double> rhoToken_;
0058   bool doRhoCorrection_;
0059   double rhoMax_;
0060   double rhoScale_;
0061   std::vector<double> effectiveAreas_;
0062   std::vector<double> absEtaLowEdges_;
0063 };
0064
0065 #endif  //HLTEgammaGenericQuadraticFilter_h