btau/plugins/HLTmumutkVtxProducer.h

0001 #ifndef HLTmumutkVtxProducer_h
0002 #define HLTmumutkVtxProducer_h
0003 //
0004 // Package:    HLTstaging
0005 // Class:      HLTmumutkVtxProducer
0006 //
0007 /**\class HLTmumutkVtxProducer
0008
0009  HLT Filter for b to (mumu) + X
0010
0011  Implementation:
0012      <Notes on implementation>
0013 */
0014
0015 // system include files
0016 #include <memory>
0017
0018 #include "FWCore/Framework/interface/stream/EDProducer.h"
0019 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0020 #include "FWCore/Framework/interface/Event.h"
0021 #include "FWCore/Framework/interface/EventSetup.h"
0022 #include "DataFormats/RecoCandidate/interface/RecoChargedCandidateFwd.h"
0023 #include "DataFormats/TrackReco/interface/TrackFwd.h"
0024 #include "DataFormats/HLTReco/interface/TriggerFilterObjectWithRefs.h"
0025 #include "DataFormats/HLTReco/interface/TriggerRefsCollections.h"
0026 #include "DataFormats/BeamSpot/interface/BeamSpot.h"
0027 #include "TrackingTools/Records/interface/TransientTrackRecord.h"
0028 #include "TrackingTools/TransientTrack/interface/TransientTrackBuilder.h"
0029 #include "MagneticField/Records/interface/IdealMagneticFieldRecord.h"
0030 #include "MagneticField/Engine/interface/MagneticField.h"
0031 #include <vector>
0032
0033 namespace edm {
0034   class ConfigurationDescriptions;
0035 }
0036
0037 // ----------------------------------------------------------------------
0038
0039 namespace reco {
0040   class Candidate;
0041   class Track;
0042 }  // namespace reco
0043
0044 class FreeTrajectoryState;
0045 class MagneticField;
0046
0047 class HLTmumutkVtxProducer : public edm::stream::EDProducer<> {
0048 public:
0049   explicit HLTmumutkVtxProducer(const edm::ParameterSet&);
0050   ~HLTmumutkVtxProducer() override;
0051   static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
0052   void produce(edm::Event&, const edm::EventSetup&) override;
0053
0054 private:
0055   bool overlap(const reco::TrackRef& trackref1, const reco::TrackRef& trackref2);
0056   static FreeTrajectoryState initialFreeState(const reco::Track&, const MagneticField*);
0057   bool checkPreviousCand(const reco::TrackRef& trackref, const std::vector<reco::RecoChargedCandidateRef>& ref2) const;
0058
0059   const edm::ESGetToken<TransientTrackBuilder, TransientTrackRecord> transientTrackRecordToken_;
0060
0061   const edm::InputTag muCandTag_;
0062   const edm::EDGetTokenT<reco::RecoChargedCandidateCollection> muCandToken_;
0063   const edm::InputTag trkCandTag_;
0064   const edm::EDGetTokenT<reco::RecoChargedCandidateCollection> trkCandToken_;
0065   const edm::InputTag previousCandTag_;
0066   const edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> previousCandToken_;
0067
0068   const std::string mfName_;
0069   const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> idealMagneticFieldRecordToken_;
0070
0071   const double thirdTrackMass_;
0072   const double maxEta_;
0073   const double minPt_;
0074   const double minInvMass_;
0075   const double maxInvMass_;
0076   const double minD0Significance_;
0077   const double overlapDR_;
0078
0079   const edm::InputTag beamSpotTag_;
0080   const edm::EDGetTokenT<reco::BeamSpot> beamSpotToken_;
0081 };
0082 #endif