#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 39 of file GEMRecHitMatcher.h.

Member Typedef Documentation

typedef matching::Digi GEMRecHitMatcher::RecHit

Definition at line 41 of file GEMRecHitMatcher.h.

typedef matching::DigiContainer GEMRecHitMatcher::RecHitContainer

Definition at line 42 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 30 of file GEMRecHitMatcher.cc.

30 {}

Member Function Documentation

std::set< unsigned int > GEMRecHitMatcher::chamberIds ( ) const

Definition at line 87 of file GEMRecHitMatcher.cc.

References chamber_to_recHits_, AlCaHLTBitMon_ParallelJobs::p, and mps_fire::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 80 of file GEMRecHitMatcher.cc.

References detid_to_recHits_, AlCaHLTBitMon_ParallelJobs::p, and mps_fire::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(), and verbose().

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

GEMRecHitMatcher::matchRecHitsToSimTrack

void matchRecHitsToSimTrack(const GEMRecHitCollection &recHits)

Definition: GEMRecHitMatcher.cc:34

edm::Handle::product

T const * product() const

Definition: Handle.h:74

GEMRecHitMatcher::gem_rechits_

edm::Handle< GEMRecHitCollection > gem_rechits_

Definition: GEMRecHitMatcher.h:88

void GEMRecHitMatcher::matchRecHitsToSimTrack ( const GEMRecHitCollection & recHits )

private

Definition at line 34 of file GEMRecHitMatcher.cc.

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

Referenced by init(), and verbose().

                                                                                 {
   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 119 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 144 of file GEMRecHitMatcher.cc.

References detIds(), mps_fire::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 170 of file GEMRecHitMatcher.cc.

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

Referenced by verbose().

                                                                                     {
   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 155 of file GEMRecHitMatcher.cc.

References GEMEtaPartition::centreOfStrip(), matching::digi_channel(), matching::digi_id(), matching::digi_type(), GEMGeometry::etaPartition(), gem_geo_, matching::GEM_STRIP, runTauDisplay::gp, GEMGeometry::idToDet(), digitizers_cfi::strip, GeomDet::surface(), protons_cff::t, and Surface::toGlobal().

Referenced by recHitMeanPosition(), and verbose().

                                                                        {
   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 107 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 101 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 113 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 132 of file GEMRecHitMatcher.cc.

References mps_fire::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 94 of file GEMRecHitMatcher.cc.

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

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