GEMRecHitMatcher Class Reference

#include <GEMRecHitMatcher.h>

Public Types
typedef matching::Digi	RecHit
typedef matching::DigiContainer	RecHitContainer

Public Member Functions
std::set< unsigned int >	chamberIds () const
std::set< unsigned int >	detIds () const
	GEMRecHitMatcher (const SimHitMatcher &sh, const edm::Event &, const GEMGeometry &geom, const edm::ParameterSet &cfg, edm::EDGetToken &)
int	nLayersWithRecHitsInSuperChamber (unsigned int) const
int	nRecHits () const
	How many recHits in GEM did this simtrack get in total? More...
std::set< int >	partitionNumbers () const
GlobalPoint	recHitMeanPosition (const RecHitContainer &rechits) const
GlobalPoint	recHitPosition (const RecHit &rechit) const
const RecHitContainer &	recHitsInChamber (unsigned int) const
const RecHitContainer &	recHitsInDetId (unsigned int) const
const RecHitContainer &	recHitsInSuperChamber (unsigned int) const
std::set< int >	stripNumbersInDetId (unsigned int) const
std::set< unsigned int >	superChamberIds () const
bool	verbose () const
	~GEMRecHitMatcher ()

Private Member Functions
void	init (const edm::Event &)
void	matchRecHitsToSimTrack (const GEMRecHitCollection &recHits)

Private Attributes
std::map< unsigned int, RecHitContainer >	chamber_to_recHits_
std::map< unsigned int, RecHitContainer >	detid_to_recHits_
const GEMGeometry &	gem_geo_
edm::Handle< GEMRecHitCollection >	gem_rechits_
int	matchDeltaStrip_
int	maxBXGEM_
int	minBXGEM_
const RecHitContainer	no_recHits_
const SimHitMatcher &	simhit_matcher_
std::map< unsigned int, RecHitContainer >	superchamber_to_recHits_
bool	verbose_

Detailed Description

Definition at line 42 of file GEMRecHitMatcher.h.

Member Typedef Documentation

typedef matching::Digi GEMRecHitMatcher::RecHit

Definition at line 46 of file GEMRecHitMatcher.h.

typedef matching::DigiContainer GEMRecHitMatcher::RecHitContainer

Definition at line 47 of file GEMRecHitMatcher.h.

Constructor & Destructor Documentation

GEMRecHitMatcher::GEMRecHitMatcher	(	const SimHitMatcher &	sh,
		const edm::Event &	e,
		const GEMGeometry &	geom,
		const edm::ParameterSet &	cfg,
		edm::EDGetToken &	gem_recHitToken
	)

Definition at line 11 of file GEMRecHitMatcher.cc.

References gem_rechits_, edm::Event::getByToken(), edm::ParameterSet::getUntrackedParameter(), init(), matchDeltaStrip_, maxBXGEM_, minBXGEM_, and verbose_.

                                   :simhit_matcher_(sh),gem_geo_(geom)
{
  minBXGEM_ = cfg.getUntrackedParameter<int>("minBXGEM", -1);
  maxBXGEM_ = cfg.getUntrackedParameter<int>("maxBXGEM", 1);

  matchDeltaStrip_ = cfg.getUntrackedParameter<int>("matchDeltaStripGEM", 1);

  //setVerbose(cfg.getUntrackedParameter<int>("verboseGEMDigi", 0));
  verbose_=false;

  e.getByToken(gem_recHitToken, gem_rechits_);

  init(e);
}

GEMRecHitMatcher::~GEMRecHitMatcher

(

)

Definition at line 28 of file GEMRecHitMatcher.cc.

{}

Member Function Documentation

std::set< unsigned int > GEMRecHitMatcher::chamberIds

(

)

const

Definition at line 100 of file GEMRecHitMatcher.cc.

References chamber_to_recHits_, AlCaHLTBitMon_ParallelJobs::p, and query::result.

Referenced by GEMRecHitTrackMatch::analyze().

{
  std::set<unsigned int> result;
  for (auto& p: chamber_to_recHits_) result.insert(p.first);
  return result;
}

std::set< unsigned int > GEMRecHitMatcher::detIds

(

void

)

const

Definition at line 91 of file GEMRecHitMatcher.cc.

References detid_to_recHits_, AlCaHLTBitMon_ParallelJobs::p, and query::result.

Referenced by partitionNumbers().

{
  std::set<unsigned int> result;
  for (auto& p: detid_to_recHits_) result.insert(p.first);
  return result;
}

void GEMRecHitMatcher::init ( const edm::Event &  e )

private

Definition at line 32 of file GEMRecHitMatcher.cc.

References gem_rechits_, matchRecHitsToSimTrack(), and edm::Handle< T >::product().

Referenced by GEMRecHitMatcher().

{
  matchRecHitsToSimTrack(*gem_rechits_.product());
}

void GEMRecHitMatcher::matchRecHitsToSimTrack ( const GEMRecHitCollection &  recHits )

private

Definition at line 39 of file GEMRecHitMatcher.cc.

References GEMDetId::chamber(), chamber_to_recHits_, GEMDetId::chamberId(), filterCSVwithJSON::copy, gather_cfg::cout, ztail::d, detid_to_recHits_, SimHitMatcher::detIdsGEM(), matching::GEM_STRIP, SimHitMatcher::hitStripsInDetId(), i, matching::make_digi(), matchDeltaStrip_, maxBXGEM_, minBXGEM_, DetId::rawId(), GEMDetId::region(), GEMDetId::ring(), simhit_matcher_, GEMDetId::station(), superchamber_to_recHits_, and verbose().

Referenced by init().

{
  
  auto det_ids = simhit_matcher_.detIdsGEM();
  for (auto id: det_ids)
  {
    GEMDetId p_id(id);
    GEMDetId superch_id(p_id.region(), p_id.ring(), p_id.station(), 1, p_id.chamber(), 0);

    auto hit_strips = simhit_matcher_.hitStripsInDetId(id, matchDeltaStrip_);
    if (verbose())
    {
      cout<<"hit_strips_fat ";
      copy(hit_strips.begin(), hit_strips.end(), ostream_iterator<int>(cout, " "));
      cout<<endl;
    }

    auto rechits_in_det = rechits.get(GEMDetId(id));

    for (auto d = rechits_in_det.first; d != rechits_in_det.second; ++d)
    {
      if (verbose()) cout<<"recHit "<<p_id<<" "<<*d<<endl;
      // check that the rechit is within BX range
      if (d->BunchX() < minBXGEM_ || d->BunchX() > maxBXGEM_) continue;
      // check that it matches a strip that was hit by SimHits from our track

      int firstStrip = d->firstClusterStrip();
      int cls = d->clusterSize();
      bool stripFound = false;

      for(int i = firstStrip; i < (firstStrip + cls); i++){

    if (hit_strips.find(i) != hit_strips.end()) stripFound = true;
        //std::cout<<i<<" "<<firstStrip<<" "<<cls<<" "<<stripFound<<std::endl;
    
      }

      if (!stripFound) continue;
      if (verbose()) cout<<"oki"<<endl;

      auto myrechit = make_digi(id, d->firstClusterStrip(), d->BunchX(), GEM_STRIP, d->clusterSize());
      detid_to_recHits_[id].push_back(myrechit);
      chamber_to_recHits_[ p_id.chamberId().rawId() ].push_back(myrechit);
      superchamber_to_recHits_[ superch_id() ].push_back(myrechit);

    }
  }
  
}

int GEMRecHitMatcher::nLayersWithRecHitsInSuperChamber

(

unsigned int

detid

)

const

Definition at line 138 of file GEMRecHitMatcher.cc.

References LayerTriplets::layers().

{
  set<int> layers;
  /*
  auto recHits = recHitsInSuperChamber(detid);
  for (auto& d: recHits)
  {
    GEMDetId idd(digi_id(d));
    layers.insert(idd.layer());
  }
  */
  return layers.size();
}

int GEMRecHitMatcher::nRecHits

(

)

const

How many recHits in GEM did this simtrack get in total?

std::set< int > GEMRecHitMatcher::partitionNumbers

(

)

const

Definition at line 168 of file GEMRecHitMatcher.cc.

References detIds(), query::result, and GEMDetId::roll().

{
  std::set<int> result;

  auto detids = detIds();
  for (auto id: detids)
  {
    GEMDetId idd(id);
    result.insert( idd.roll() );
  }
  return result;
}

GlobalPoint GEMRecHitMatcher::recHitMeanPosition

(

const RecHitContainer &

rechits

)

const

Definition at line 201 of file GEMRecHitMatcher.cc.

References ztail::d, gen::n, recHitPosition(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

{
  GlobalPoint point_zero;
  if (rechit.empty()) return point_zero; // point "zero"

  float sumx, sumy, sumz;
  sumx = sumy = sumz = 0.f;
  size_t n = 0;
  for (auto& d: rechit)
  {
    GlobalPoint gp = recHitPosition(d);
    if (gp == point_zero) continue;

    sumx += gp.x();
    sumy += gp.y();
    sumz += gp.z();
    ++n;
  }
  if (n == 0) return GlobalPoint();
  return GlobalPoint(sumx/n, sumy/n, sumz/n);
}

GlobalPoint GEMRecHitMatcher::recHitPosition

(

const RecHit &

rechit

)

const

Definition at line 182 of file GEMRecHitMatcher.cc.

References GEMEtaPartition::centreOfStrip(), matching::digi_channel(), matching::digi_id(), matching::digi_type(), GEMGeometry::etaPartition(), gem_geo_, matching::GEM_STRIP, GEMGeometry::idToDet(), GeomDet::surface(), lumiQTWidget::t, and Surface::toGlobal().

Referenced by recHitMeanPosition().

{
  unsigned int id = digi_id(rechit);
  int strip = digi_channel(rechit);
  DigiType t = digi_type(rechit);

  GlobalPoint gp;
  if ( t == GEM_STRIP )
  {
    GEMDetId idd(id);
    LocalPoint lp = gem_geo_.etaPartition(idd)->centreOfStrip(strip);
    gp = gem_geo_.idToDet(id)->surface().toGlobal(lp);
  }

  return gp;
}

const GEMRecHitMatcher::RecHitContainer & GEMRecHitMatcher::recHitsInChamber

(

unsigned int

detid

)

const

Definition at line 124 of file GEMRecHitMatcher.cc.

References chamber_to_recHits_, and no_recHits_.

{
  if (chamber_to_recHits_.find(detid) == chamber_to_recHits_.end()) return no_recHits_;
  return chamber_to_recHits_.at(detid);
}

const GEMRecHitMatcher::RecHitContainer & GEMRecHitMatcher::recHitsInDetId

(

unsigned int

detid

)

const

Definition at line 117 of file GEMRecHitMatcher.cc.

References detid_to_recHits_, and no_recHits_.

{
  if (detid_to_recHits_.find(detid) == detid_to_recHits_.end()) return no_recHits_;
  return detid_to_recHits_.at(detid);
}

const GEMRecHitMatcher::RecHitContainer & GEMRecHitMatcher::recHitsInSuperChamber

(

unsigned int

detid

)

const

Definition at line 131 of file GEMRecHitMatcher.cc.

References no_recHits_, and superchamber_to_recHits_.

{
  if (superchamber_to_recHits_.find(detid) == superchamber_to_recHits_.end()) return no_recHits_;
  return superchamber_to_recHits_.at(detid);
}

std::set< int > GEMRecHitMatcher::stripNumbersInDetId

(

unsigned int

detid

)

const

Definition at line 154 of file GEMRecHitMatcher.cc.

References query::result.

{
  set<int> result;
  /*
  auto recHits = recHitsInDetId(detid);
  for (auto& d: recHits)
  {
    result.insert( digi_channel(d) );
  }
  */
  return result;
}

std::set< unsigned int > GEMRecHitMatcher::superChamberIds

(

)

const

Definition at line 108 of file GEMRecHitMatcher.cc.

References AlCaHLTBitMon_ParallelJobs::p, query::result, and superchamber_to_recHits_.

{
  std::set<unsigned int> result;
  for (auto& p: superchamber_to_recHits_) result.insert(p.first);
  return result;
}

bool GEMRecHitMatcher::verbose ( ) const

inline

Definition at line 79 of file GEMRecHitMatcher.h.

References verbose_.

Referenced by matchRecHitsToSimTrack().

{ return verbose_; }

Member Data Documentation

std::map<unsigned int, RecHitContainer> GEMRecHitMatcher::chamber_to_recHits_

private

Definition at line 101 of file GEMRecHitMatcher.h.

Referenced by chamberIds(), matchRecHitsToSimTrack(), and recHitsInChamber().

std::map<unsigned int, RecHitContainer> GEMRecHitMatcher::detid_to_recHits_

private

Definition at line 100 of file GEMRecHitMatcher.h.

Referenced by detIds(), matchRecHitsToSimTrack(), and recHitsInDetId().

const GEMGeometry& GEMRecHitMatcher::gem_geo_

private

Definition at line 93 of file GEMRecHitMatcher.h.

Referenced by recHitPosition().

edm::Handle<GEMRecHitCollection> GEMRecHitMatcher::gem_rechits_

private

Definition at line 90 of file GEMRecHitMatcher.h.

Referenced by GEMRecHitMatcher(), and init().

int GEMRecHitMatcher::matchDeltaStrip_

private

Definition at line 98 of file GEMRecHitMatcher.h.

Referenced by GEMRecHitMatcher(), and matchRecHitsToSimTrack().

int GEMRecHitMatcher::maxBXGEM_

private

Definition at line 95 of file GEMRecHitMatcher.h.

Referenced by GEMRecHitMatcher(), and matchRecHitsToSimTrack().

int GEMRecHitMatcher::minBXGEM_

private

Definition at line 95 of file GEMRecHitMatcher.h.

Referenced by GEMRecHitMatcher(), and matchRecHitsToSimTrack().

const RecHitContainer GEMRecHitMatcher::no_recHits_

private

Definition at line 104 of file GEMRecHitMatcher.h.

Referenced by recHitsInChamber(), recHitsInDetId(), and recHitsInSuperChamber().

const SimHitMatcher& GEMRecHitMatcher::simhit_matcher_

private

Definition at line 92 of file GEMRecHitMatcher.h.

Referenced by matchRecHitsToSimTrack().

std::map<unsigned int, RecHitContainer> GEMRecHitMatcher::superchamber_to_recHits_

private

Definition at line 102 of file GEMRecHitMatcher.h.

Referenced by matchRecHitsToSimTrack(), recHitsInSuperChamber(), and superChamberIds().

bool GEMRecHitMatcher::verbose_

private

Definition at line 96 of file GEMRecHitMatcher.h.

Referenced by GEMRecHitMatcher(), and verbose().