DTChamberEfficiency.h

Go to the documentation of this file.

#ifndef DTCHAMBEREFFICIENCY_H
#define DTCHAMBEREFFICIENCY_H

#include "DataFormats/Common/interface/Handle.h"


#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "DataFormats/TrackReco/interface/Track.h"

#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Utilities/interface/InputTag.h"

#include "FWCore/Framework/interface/Frameworkfwd.h"

#include "Geometry/DTGeometry/interface/DTGeometry.h"
#include "Geometry/CommonDetUnit/interface/GlobalTrackingGeometry.h"
#include "MagneticField/Engine/interface/MagneticField.h"

#include "RecoMuon/MeasurementDet/interface/MuonDetLayerMeasurements.h"
#include <vector>

namespace reco {
  class TransientTrack;
}



class Chi2MeasurementEstimator;
class MuonServiceProxy;

class DQMStore;
class MonitorElement;
class FreeTrajectoryState;
class DetLayer;
class DetId;

class DTChamberEfficiency : public edm::EDAnalyzer
{

 public:
  //Constructor 
  DTChamberEfficiency(const edm::ParameterSet& pset) ;

  //Destructor
  ~DTChamberEfficiency() ;

  //Operations
  void analyze(const edm::Event & event, const edm::EventSetup& eventSetup);
  void beginJob();
  void beginRun(const edm::Run& , const edm::EventSetup&);
  void endJob();

 private:

  //functions
  std::vector<const DetLayer*> compatibleLayers(const DetLayer *initialLayer,
                        const FreeTrajectoryState& fts, PropagationDirection propDir);


  void bookHistos();
  MeasurementContainer segQualityCut(const MeasurementContainer seg_list) const;
  bool chamberSelection(const DetId& idDetLay, reco::TransientTrack& trans_track) const;
  inline edm::ESHandle<Propagator> propagator() const;

  //data members
  bool debug;

  edm::InputTag theTracksLabel;

  edm::InputTag labelRPCRecHits;
  edm::InputTag thedt4DSegments;
  edm::InputTag thecscSegments;

  double theMaxChi2;
  double theNSigma;
  int theMinNrec;

  std::string theNavigationType;

  edm::ESHandle<DTGeometry> dtGeom;

  DQMStore* theDbe;

  MuonServiceProxy* theService;
  MuonDetLayerMeasurements* theMeasurementExtractor;
  Chi2MeasurementEstimator* theEstimator;

  edm::ESHandle<MagneticField> magfield;
  edm::ESHandle<GlobalTrackingGeometry> theTrackingGeometry;

  //std::map<DTChamberId, std::vector<MonitorElement*> > histosPerW;

  std::vector<std::vector<MonitorElement*> > histosPerW;

 protected:

};

#endif // DTANALYZER_H