MuonSimHitMatcher.h

Go to the documentation of this file.

#ifndef Validation_MuonHits_MuonSimHitMatcher_h
#define Validation_MuonHits_MuonSimHitMatcher_h

#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/Math/interface/deltaPhi.h"
#include "FWCore/Framework/interface/ConsumesCollector.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "Geometry/CommonDetUnit/interface/TrackingGeometry.h"
#include "Geometry/CommonTopologies/interface/TrapezoidalStripTopology.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"
#include "SimDataFormats/Track/interface/SimTrackContainer.h"
#include "SimDataFormats/TrackingHit/interface/PSimHitContainer.h"
#include "SimDataFormats/Vertex/interface/SimVertexContainer.h"
#include "Validation/MuonHits/interface/MuonHitHelper.h"

#include <map>
#include <set>
#include <vector>

class MuonSimHitMatcher {
public:
  // constructor
  MuonSimHitMatcher(const edm::ParameterSet& iPS, edm::ConsumesCollector&& iC);

  // destructor
  ~MuonSimHitMatcher() {}

  // initialize the event
  void init(const edm::Event& e, const edm::EventSetup& eventSetup);

  // do the matching
  void match(const SimTrack& t, const SimVertex& v);

  // access to all the Muon SimHits (use MuonSubdetId::SubSystem)
  const edm::PSimHitContainer& simHits(int) const;

  // partitions' detIds with SimHits
  std::set<unsigned int> detIds(int type = 0) const;

  // chamber detIds with SimHits
  std::set<unsigned int> chamberIds(int type = 0) const;

  // simhits from a particular partition, chamber
  const edm::PSimHitContainer& hitsInDetId(unsigned int) const;
  const edm::PSimHitContainer& hitsInChamber(unsigned int) const;

  // calculate Global average position for a provided collection of simhits
  GlobalPoint simHitsMeanPosition(const edm::PSimHitContainer& sim_hits) const;

  // calculate Global average momentum for a provided collection of simhits in
  // CSC
  GlobalVector simHitsMeanMomentum(const edm::PSimHitContainer& sim_hits) const;

  // calculate the average position at the second station
  GlobalPoint simHitsMeanPositionStation(int n) const;

  const TrackingGeometry* geometry() { return geometry_; }

protected:
  std::vector<unsigned int> getIdsOfSimTrackShower(unsigned trk_id,
                                                   const edm::SimTrackContainer& simTracks,
                                                   const edm::SimVertexContainer& simVertices);

  void clear();

  bool verboseSimTrack_;
  bool simMuOnly_;
  bool discardEleHits_;
  bool verbose_;
  bool hasGeometry_;

  const TrackingGeometry* geometry_;

  edm::EDGetTokenT<edm::SimVertexContainer> simVertexInput_;
  edm::EDGetTokenT<edm::SimTrackContainer> simTrackInput_;
  edm::EDGetTokenT<edm::PSimHitContainer> simHitInput_;

  edm::Handle<edm::SimTrackContainer> simTracksH_;
  edm::Handle<edm::SimVertexContainer> simVerticesH_;
  edm::Handle<edm::PSimHitContainer> simHitsH_;

  edm::SimTrackContainer simTracks_;
  edm::SimVertexContainer simVertices_;
  // input collection
  edm::PSimHitContainer simHits_;

  std::vector<unsigned> track_ids_;
  std::map<unsigned int, unsigned int> trkid_to_index_;

  edm::PSimHitContainer no_hits_;

  // selected hits
  edm::PSimHitContainer hits_;
  std::map<unsigned int, edm::PSimHitContainer> detid_to_hits_;
  std::map<unsigned int, edm::PSimHitContainer> chamber_to_hits_;

  edm::ParameterSet simHitPSet_;
};

#endif