MuonDetLayerMeasurements.h

Go to the documentation of this file.

#ifndef MeasurementDet_MuonDetLayerMeasurements_H
#define MeasurementDet_MuonDetLayerMeasurements_H

#include "FWCore/Framework/interface/Event.h"
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
#include "TrackingTools/GeomPropagators/interface/Propagator.h"
#include "TrackingTools/PatternTools/interface/MeasurementEstimator.h"
#include "RecoMuon/TransientTrackingRecHit/interface/MuonTransientTrackingRecHit.h"
//#include "TrackingTools/ementDet/interface/TrajectoryMeasurement.h"
#include "TrackingTools/MeasurementDet/interface/TrajectoryMeasurementGroup.h"
#include "DataFormats/DTRecHit/interface/DTRecSegment4DCollection.h"
#include "DataFormats/CSCRecHit/interface/CSCSegmentCollection.h"
#include "DataFormats/RPCRecHit/interface/RPCRecHitCollection.h"
#include "DataFormats/RPCRecHit/interface/RPCRecHitCollection.h"
#include "FWCore/Utilities/interface/InputTag.h"

#include <vector>

class DetLayer;
class GeomDet;
class TrajectoryMeasurement;


//FIXME: these typedefs MUST GO inside the scope of MuonDetLayerMeasurements
typedef std::vector<TrajectoryMeasurement>          MeasurementContainer;
typedef std::pair<const GeomDet*,TrajectoryStateOnSurface> DetWithState;


class MuonDetLayerMeasurements {
 public:
  typedef MuonTransientTrackingRecHit::MuonRecHitContainer MuonRecHitContainer;

  MuonDetLayerMeasurements(edm::InputTag dtlabel,
               edm::InputTag csclabel,
               edm::InputTag rpclabel,
               bool enableDT = true,
               bool enableCSC = true,
               bool enableRPC = true);
  
  virtual ~MuonDetLayerMeasurements();
  
  // for a given det and state.  Not clear when the fastMeasurements below
  //  should be used, since it isn't passed a GeomDet
  MeasurementContainer
    measurements( const DetLayer* layer,
                  const GeomDet * det,
                  const TrajectoryStateOnSurface& stateOnDet,
                  const MeasurementEstimator& est,
                  const edm::Event& iEvent);

  MeasurementContainer
    measurements( const DetLayer* layer,
          const TrajectoryStateOnSurface& startingState,
          const Propagator& prop,
          const MeasurementEstimator& est,
          const edm::Event& iEvent);

  MeasurementContainer
    fastMeasurements( const DetLayer* layer,
              const TrajectoryStateOnSurface& theStateOnDet,
              const TrajectoryStateOnSurface& startingState,
              const Propagator& prop,
              const MeasurementEstimator& est,
              const edm::Event& iEvent);

  MeasurementContainer
    measurements( const DetLayer* layer,
          const TrajectoryStateOnSurface& startingState,
          const Propagator& prop,
          const MeasurementEstimator& est);

  MeasurementContainer
    fastMeasurements( const DetLayer* layer,
              const TrajectoryStateOnSurface& theStateOnDet,
              const TrajectoryStateOnSurface& startingState,
              const Propagator& prop,
              const MeasurementEstimator& est);

  std::vector<TrajectoryMeasurementGroup>
    groupedMeasurements( const DetLayer* layer,
                  const TrajectoryStateOnSurface& startingState,
                  const Propagator& prop,
                  const MeasurementEstimator& est,
                  const edm::Event& iEvent);

  std::vector<TrajectoryMeasurementGroup>
    groupedMeasurements( const DetLayer* layer,
                  const TrajectoryStateOnSurface& startingState,
                  const Propagator& prop,
                  const MeasurementEstimator& est);
 
  void setEvent(const edm::Event &);  

  MuonRecHitContainer recHits(const DetLayer* layer, const edm::Event& iEvent);

  MuonRecHitContainer recHits(const DetLayer* layer);


 private:

  MuonRecHitContainer recHits(const GeomDet*, const edm::Event& iEvent);

  void checkEvent() const;

  edm::InputTag theDTRecHitLabel;
  edm::InputTag theCSCRecHitLabel;
  edm::InputTag theRPCRecHitLabel;

  bool enableDTMeasurement;
  bool enableCSCMeasurement;
  bool enableRPCMeasurement;
  
  // caches that should get filled once per event
  edm::Handle<DTRecSegment4DCollection> theDTRecHits;
  edm::Handle<CSCSegmentCollection>     theCSCRecHits;
  edm::Handle<RPCRecHitCollection>      theRPCRecHits;

  void checkDTRecHits();
  void checkCSCRecHits();
  void checkRPCRecHits();

  // keeps track of which event the cache holds
  edm::EventID theDTEventID;
  edm::EventID theCSCEventID;
  edm::EventID theRPCEventID;

  const edm::Event* theEvent;   

  // strings to uniquely identify current process
  std::string theDTCheckName;
  std::string theRPCCheckName;
  std::string theCSCCheckName;
};
#endif