special/plugins/HLTHFAsymmetryFilter.h

0001 #ifndef _HLTHFAsymmetryFilter_H
0002 #define _HLTHFAsymmetryFilter_H
0003
0004 ///////////////////////////////////////////////////////
0005 //
0006 // HLTHFAsymetryFilter
0007 //
0008 // Filter definition
0009 //
0010 // We perform a selection on HF energy repartition
0011 //
0012 // This filter is primarily used to select Beamgas (aka PKAM) events
0013 //
0014 // An asymmetry parameter, based on the pixel clusters, is computed as follows
0015 //
0016 //  asym1 = E_HF-/(E_HF- + E_HF+) for beam1
0017 //  asym2 = E_HF+/(E_HF- + E_HF+) for beam2
0018 //
0019 // where E_HF is the total energy of clusters passing a certain threshold (given by eCut_HF_)
0020 //
0021 //  Usually for PKAM events, asym1 is close to 1. for B1 BGas events, and close to 0 for B2 BGAS events
0022 //
0023 //
0024 // More details:
0025 // http://sviret.web.cern.ch/sviret/Welcome.php?n=CMS.MIB
0026 //
0027 // S.Viret: 12/01/2011 (viret@in2p3.fr)
0028 //
0029 ///////////////////////////////////////////////////////
0030
0031 // system include files
0032 #include <memory>
0033
0034 // user include files
0035 #include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
0036 #include "FWCore/Framework/interface/Frameworkfwd.h"
0037 #include "FWCore/Framework/interface/Event.h"
0038 #include "FWCore/Framework/interface/global/EDFilter.h"
0039 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0040 #include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
0041 #include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
0042 #include "FWCore/Utilities/interface/InputTag.h"
0043
0044 //
0045 // class decleration
0046 //
0047
0048 class HLTHFAsymmetryFilter : public edm::global::EDFilter<> {
0049 public:
0050   explicit HLTHFAsymmetryFilter(const edm::ParameterSet &);
0051   ~HLTHFAsymmetryFilter() override;
0052
0053   static void fillDescriptions(edm::ConfigurationDescriptions &descriptions);
0054   bool filter(edm::StreamID, edm::Event &, const edm::EventSetup &) const override;
0055
0056 private:
0057   // ----------member data ---------------------------
0058
0059   const edm::EDGetTokenT<HFRecHitCollection> HFHitsToken_;
0060   const double eCut_HF_;
0061   const double os_asym_;
0062   const double ss_asym_;
0063 };
0064
0065 #endif  // _HLTHFAsymmetryFilter_H