MuonTCMETValueMapProducer.h

Go to the documentation of this file.

#ifndef RecoMET_MuonTCMETValueMapProducer_h
#define RecoMET_MuonTCMETValueMapProducer_h

// -*- C++ -*-
//
// Package:    Test
// Class:      Test
// 
//
// Created:  Wed Aug 29 2007
//
//

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "TH2.h"
#include "TVector3.h"
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"

class TCMETAlgo;

namespace cms {
class MuonTCMETValueMapProducer : public edm::EDProducer {
   public:
      explicit MuonTCMETValueMapProducer(const edm::ParameterSet&);
      ~MuonTCMETValueMapProducer();

//      TH2D* getResponseFunction_fit ( );
//      TH2D* getResponseFunction_mode ( );

   private:
      virtual void beginJob() ;
      virtual void produce(edm::Event&, const edm::EventSetup&);
      virtual void endJob() ;
      
      // ----------member data ---------------------------
      edm::Handle<reco::MuonCollection>    muon_h;
      edm::Handle<reco::BeamSpot>          beamSpot_h;
      edm::Handle<reco::VertexCollection>  VertexHandle;

      edm::InputTag muonInputTag_;  
      edm::InputTag beamSpotInputTag_;
      edm::InputTag vertexInputTag_;

      const class MagneticField* bField;

      const reco::VertexCollection *vertexColl;

      class TH2D* response_function;

      bool muonGlobal_;
      bool muonTracker_;
      bool useCaloMuons_;
      bool hasValidVertex;

      int     rfType_;
      int     nLayers_;
      int     nLayersTight_;
      int     vertexNdof_;
      double  vertexZ_;
      double  vertexRho_;
      double  vertexMaxDZ_;
      double  maxpt_eta25_;
      double  maxpt_eta20_;
      int     maxTrackAlgo_;
      double  minpt_;
      double  maxpt_;
      double  maxeta_;
      double  maxchi2_;
      double  minhits_;
      double  maxPtErr_;
      double  maxd0cut_;
      double  maxchi2_tight_;
      double  minhits_tight_;
      double  maxPtErr_tight_;
      double  d0cuta_;
      double  d0cutb_;
      bool    usePvtxd0_;      
      std::vector<int> trkQuality_;
      std::vector<int> trkAlgos_;

      int     muonMinValidStaHits_;
      double  muonpt_;
      double  muoneta_;
      double  muonchi2_;
      double  muonhits_;
      double  muond0_;
      double  muonDeltaR_;
      double  muon_dptrel_;
      TCMETAlgo *tcmetAlgo_;

      //functions
      bool isGoodMuon( const reco::Muon* );
      bool isGoodCaloMuon( const reco::Muon*, const unsigned int );
      bool isGoodTrack( const reco::Muon* );
      class TVector3 propagateTrack( const reco::Muon* );
      int nLayers(const reco::TrackRef);
      bool isValidVertex();
  };
}
#endif