L2MuonSeedGeneratorFromL1T.h

Go to the documentation of this file.

#ifndef RecoMuon_L2MuonSeedGenerator_L2MuonSeedGeneratorFromL1T_H
#define RecoMuon_L2MuonSeedGenerator_L2MuonSeedGeneratorFromL1T_H

//-------------------------------------------------
//
//L2MuonSeedGeneratorFromL1T
//--------------------------------------------------

#include "FWCore/Framework/interface/stream/EDProducer.h"
#include "FWCore/Utilities/interface/InputTag.h"

// Data Formats 
#include "DataFormats/MuonSeed/interface/L2MuonTrajectorySeed.h"
#include "DataFormats/MuonSeed/interface/L2MuonTrajectorySeedCollection.h"
#include "DataFormats/TrajectoryState/interface/PTrajectoryStateOnDet.h"
#include "DataFormats/MuonDetId/interface/DTChamberId.h"
#include "DataFormats/MuonDetId/interface/CSCDetId.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/GeometrySurface/interface/BoundCylinder.h"
#include "DataFormats/Math/interface/deltaR.h"
#include "DataFormats/L1Trigger/interface/Muon.h"

#include "CLHEP/Vector/ThreeVector.h"

#include "Geometry/CommonDetUnit/interface/GeomDetEnumerators.h"

class MuonServiceProxy;
class MeasurementEstimator;
class TrajectorySeed;
class TrajectoryStateOnSurface;

namespace edm {class ParameterSet; class Event; class EventSetup;}

class L2MuonSeedGeneratorFromL1T : public edm::stream::EDProducer<> {
 
 public:
  
  explicit L2MuonSeedGeneratorFromL1T(const edm::ParameterSet&);

  ~L2MuonSeedGeneratorFromL1T() override;

  static void fillDescriptions(edm::ConfigurationDescriptions & descriptions);
  void produce(edm::Event&, const edm::EventSetup&) override;
  
 private:

  edm::InputTag theSource;
  edm::InputTag theL1GMTReadoutCollection;
  edm::InputTag theOfflineSeedLabel;
  std::string   thePropagatorName;

  edm::EDGetTokenT<l1t::MuonBxCollection> muCollToken_;
  edm::EDGetTokenT<edm::View<TrajectorySeed> > offlineSeedToken_;

  const double theL1MinPt;
  const double theL1MaxEta;
  const unsigned theL1MinQuality;
  const double theMinPtBarrel;
  const double theMinPtEndcap;
  const bool useOfflineSeed;
  const bool useUnassociatedL1;
  std::vector<double> matchingDR;
  std::vector<double> etaBins;

  bool centralBxOnly_;
  const unsigned matchType;
  const unsigned sortType;

  MuonServiceProxy *theService;  

  MeasurementEstimator *theEstimator;

  const TrajectorySeed* associateOfflineSeedToL1( edm::Handle<edm::View<TrajectorySeed> > &, 
                          std::vector<int> &, 
                          TrajectoryStateOnSurface &,
                          double );

  bool isAssociateOfflineSeedToL1( edm::Handle<edm::View<TrajectorySeed> > &,
              std::vector< std::vector<double> > &,
              TrajectoryStateOnSurface &,
              unsigned int,
              std::vector< std::vector<const TrajectorySeed *> > &,
              double );

};

#endif