CMS 3D CMS Logo

List of all members | Public Types | Public Member Functions | Private Attributes
CSCHitAssociator Class Reference

#include <CSCHitAssociator.h>

Public Types

typedef edm::DetSetVector< StripDigiSimLinkDigiSimLinks
 
typedef edm::DetSet< StripDigiSimLinkLayerLinks
 
typedef std::pair< uint32_t, EncodedEventIdSimHitIdpr
 
typedef edm::DetSetVector< StripDigiSimLinkWireDigiSimLinks
 

Public Member Functions

std::vector< SimHitIdprassociateCSCHitId (const CSCRecHit2D *) const
 
std::vector< SimHitIdprassociateHitId (const TrackingRecHit &) const
 
 CSCHitAssociator (const edm::Event &, const edm::EventSetup &, const edm::ParameterSet &)
 
 CSCHitAssociator (const edm::ParameterSet &, edm::ConsumesCollector &&iC)
 
void initEvent (const edm::Event &, const edm::EventSetup &)
 

Private Attributes

const CSCGeometrycscgeom
 
edm::InputTag linksTag
 
const DigiSimLinkstheDigiSimLinks
 

Detailed Description

Definition at line 24 of file CSCHitAssociator.h.

Member Typedef Documentation

Definition at line 26 of file CSCHitAssociator.h.

Definition at line 28 of file CSCHitAssociator.h.

typedef std::pair<uint32_t, EncodedEventId> CSCHitAssociator::SimHitIdpr

Definition at line 29 of file CSCHitAssociator.h.

Definition at line 27 of file CSCHitAssociator.h.

Constructor & Destructor Documentation

CSCHitAssociator::CSCHitAssociator ( const edm::Event event,
const edm::EventSetup setup,
const edm::ParameterSet conf 
)

Definition at line 8 of file CSCHitAssociator.cc.

References initEvent().

9  : theDigiSimLinks(nullptr), linksTag(conf.getParameter<edm::InputTag>("CSClinksTag")) {
10  initEvent(event, setup);
11 }
T getParameter(std::string const &) const
edm::InputTag linksTag
const DigiSimLinks * theDigiSimLinks
void initEvent(const edm::Event &, const edm::EventSetup &)
CSCHitAssociator::CSCHitAssociator ( const edm::ParameterSet conf,
edm::ConsumesCollector &&  iC 
)

Definition at line 13 of file CSCHitAssociator.cc.

References linksTag.

14  : theDigiSimLinks(nullptr), linksTag(conf.getParameter<edm::InputTag>("CSClinksTag")) {
16 }
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
T getParameter(std::string const &) const
edm::InputTag linksTag
edm::DetSetVector< StripDigiSimLink > DigiSimLinks
const DigiSimLinks * theDigiSimLinks

Member Function Documentation

std::vector< CSCHitAssociator::SimHitIdpr > CSCHitAssociator::associateCSCHitId ( const CSCRecHit2D cscrechit) const

Definition at line 30 of file CSCHitAssociator.cc.

References CSCLayerGeometry::channel(), CSCRecHit2D::channels(), CSCRecHit2D::cscDetId(), cscgeom, edm::DetSetVector< T >::end(), spr::find(), edm::DetSetVector< T >::find(), TrackingRecHit::geographicalId(), geometry, CSCGeometry::layer(), LogTrace, CSCRecHit2D::nStrips(), DetId::rawId(), and theDigiSimLinks.

Referenced by MuonAssociatorByHitsHelper::getMatchedIds().

30  {
31  std::vector<SimHitIdpr> simtrackids;
32 
33  unsigned int detId = cscrechit->geographicalId().rawId();
34  int nchannels = cscrechit->nStrips();
35  const CSCLayerGeometry *laygeom = cscgeom->layer(cscrechit->cscDetId())->geometry();
36 
38 
39  if (layerLinks != theDigiSimLinks->end()) {
40  for (int idigi = 0; idigi < nchannels; ++idigi) {
41  // strip and readout channel numbers may differ in ME1/1A
42  int istrip = cscrechit->channels(idigi);
43  int channel = laygeom->channel(istrip);
44 
45  for (LayerLinks::const_iterator link = layerLinks->begin(); link != layerLinks->end(); ++link) {
46  int ch = static_cast<int>(link->channel());
47  if (ch == channel) {
48  SimHitIdpr currentId(link->SimTrackId(), link->eventId());
49  if (find(simtrackids.begin(), simtrackids.end(), currentId) == simtrackids.end())
50  simtrackids.push_back(currentId);
51  }
52  }
53  }
54 
55  } else
56  LogTrace("CSCHitAssociator") << "*** WARNING in CSCHitAssociator::associateCSCHitId - CSC layer " << detId
57  << " has no DigiSimLinks !" << std::endl;
58 
59  return simtrackids;
60 }
CSCDetId cscDetId() const
Definition: CSCRecHit2D.h:52
iterator find(det_id_type id)
Definition: DetSetVector.h:290
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:50
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:20
int channels(unsigned int i) const
Extracting strip channel numbers comprising the rechit - low.
Definition: CSCRecHit2D.h:55
unsigned int nStrips() const
Definition: CSCRecHit2D.h:56
iterator end()
Return the off-the-end iterator.
Definition: DetSetVector.h:361
#define LogTrace(id)
int channel(int strip) const
std::pair< uint32_t, EncodedEventId > SimHitIdpr
const DigiSimLinks * theDigiSimLinks
ESHandle< TrackerGeometry > geometry
const CSCLayer * layer(CSCDetId id) const
Return the layer corresponding to given DetId.
Definition: CSCGeometry.cc:139
const CSCGeometry * cscgeom
DetId geographicalId() const
collection_type::const_iterator const_iterator
Definition: DetSet.h:33
collection_type::const_iterator const_iterator
Definition: DetSetVector.h:104
std::vector< CSCHitAssociator::SimHitIdpr > CSCHitAssociator::associateHitId ( const TrackingRecHit hit) const

Definition at line 62 of file CSCHitAssociator.cc.

References CSCLayerGeometry::channel(), CSCRecHit2D::channels(), CSCRecHit2D::cscDetId(), cscgeom, edm::DetSetVector< T >::end(), spr::find(), edm::DetSetVector< T >::find(), TrackingRecHit::geographicalId(), geometry, CSCGeometry::layer(), LogTrace, CSCRecHit2D::nStrips(), DetId::rawId(), and theDigiSimLinks.

62  {
63  std::vector<SimHitIdpr> simtrackids;
64 
65  const TrackingRecHit *hitp = &hit;
66  const CSCRecHit2D *cscrechit = dynamic_cast<const CSCRecHit2D *>(hitp);
67 
68  if (cscrechit) {
69  unsigned int detId = cscrechit->geographicalId().rawId();
70  int nchannels = cscrechit->nStrips();
71  const CSCLayerGeometry *laygeom = cscgeom->layer(cscrechit->cscDetId())->geometry();
72 
74 
75  if (layerLinks != theDigiSimLinks->end()) {
76  for (int idigi = 0; idigi < nchannels; ++idigi) {
77  // strip and readout channel numbers may differ in ME1/1A
78  int istrip = cscrechit->channels(idigi);
79  int channel = laygeom->channel(istrip);
80 
81  for (LayerLinks::const_iterator link = layerLinks->begin(); link != layerLinks->end(); ++link) {
82  int ch = static_cast<int>(link->channel());
83  if (ch == channel) {
84  SimHitIdpr currentId(link->SimTrackId(), link->eventId());
85  if (find(simtrackids.begin(), simtrackids.end(), currentId) == simtrackids.end())
86  simtrackids.push_back(currentId);
87  }
88  }
89  }
90 
91  } else
92  LogTrace("CSCHitAssociator") << "*** WARNING in CSCHitAssociator::associateHitId - CSC layer " << detId
93  << " has no DigiSimLinks !" << std::endl;
94 
95  } else
96  LogTrace("CSCHitAssociator") << "*** WARNING in CSCHitAssociator::associateHitId, null dynamic_cast "
97  "!";
98 
99  return simtrackids;
100 }
CSCDetId cscDetId() const
Definition: CSCRecHit2D.h:52
iterator find(det_id_type id)
Definition: DetSetVector.h:290
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:50
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:20
int channels(unsigned int i) const
Extracting strip channel numbers comprising the rechit - low.
Definition: CSCRecHit2D.h:55
unsigned int nStrips() const
Definition: CSCRecHit2D.h:56
iterator end()
Return the off-the-end iterator.
Definition: DetSetVector.h:361
#define LogTrace(id)
int channel(int strip) const
std::pair< uint32_t, EncodedEventId > SimHitIdpr
const DigiSimLinks * theDigiSimLinks
ESHandle< TrackerGeometry > geometry
const CSCLayer * layer(CSCDetId id) const
Return the layer corresponding to given DetId.
Definition: CSCGeometry.cc:139
const CSCGeometry * cscgeom
DetId geographicalId() const
collection_type::const_iterator const_iterator
Definition: DetSet.h:33
collection_type::const_iterator const_iterator
Definition: DetSetVector.h:104
void CSCHitAssociator::initEvent ( const edm::Event event,
const edm::EventSetup setup 
)

Definition at line 18 of file CSCHitAssociator.cc.

References cscgeom, TTClusterAssociation_cfi::digiSimLinks, edm::EventSetup::get(), linksTag, LogTrace, edm::Handle< T >::product(), and theDigiSimLinks.

Referenced by CSCHitAssociator().

18  {
20  LogTrace("CSCHitAssociator") << "getting CSC Strip DigiSimLink collection - " << linksTag;
21  event.getByLabel(linksTag, digiSimLinks);
22  theDigiSimLinks = digiSimLinks.product();
23 
24  // get CSC Geometry to use CSCLayer methods
26  setup.get<MuonGeometryRecord>().get(mugeom);
27  cscgeom = &*mugeom;
28 }
edm::InputTag linksTag
#define LogTrace(id)
T const * product() const
Definition: Handle.h:74
const DigiSimLinks * theDigiSimLinks
T get() const
Definition: EventSetup.h:71
const CSCGeometry * cscgeom

Member Data Documentation

const CSCGeometry* CSCHitAssociator::cscgeom
private

Definition at line 44 of file CSCHitAssociator.h.

Referenced by associateCSCHitId(), associateHitId(), and initEvent().

edm::InputTag CSCHitAssociator::linksTag
private

Definition at line 42 of file CSCHitAssociator.h.

Referenced by CSCHitAssociator(), and initEvent().

const DigiSimLinks* CSCHitAssociator::theDigiSimLinks
private

Definition at line 40 of file CSCHitAssociator.h.

Referenced by associateCSCHitId(), associateHitId(), and initEvent().