#include <GEMDigiMatcher.h>

Public Member Functions
std::set< unsigned int >	chamberIds () const

std::set< unsigned int >	chamberIdsWithPads () const

std::set< int >	coPadNumbersInDetId (unsigned int) const

const DigiContainer &	coPadsInDetId (unsigned int) const

const DigiContainer &	coPadsInSuperChamber (unsigned int) const

std::set< unsigned int >	detIds () const

std::set< unsigned int >	detIdsWithCoPads () const

const DigiContainer &	digisInChamber (unsigned int) const

const DigiContainer &	digisInDetId (unsigned int) const

const DigiContainer &	digisInSuperChamber (unsigned int) const

	GEMDigiMatcher (const SimHitMatcher &sh, const edm::Event &, const GEMGeometry &geom, const edm::ParameterSet &cfg, edm::EDGetToken &, edm::EDGetToken &, edm::EDGetToken &)

int	nCoPads () const
	How many coincidence pads in GEM did this simtrack get in total? More...

int	nLayersWithDigisInSuperChamber (unsigned int) const

int	nLayersWithPadsInSuperChamber (unsigned int) const

int	nPads () const
	How many pads in GEM did this simtrack get in total? More...

std::set< int >	padNumbersInDetId (unsigned int) const

const DigiContainer &	padsInChamber (unsigned int) const

const DigiContainer &	padsInDetId (unsigned int) const

const DigiContainer &	padsInSuperChamber (unsigned int) const

std::set< int >	partitionNumbers () const

std::set< int >	partitionNumbersWithCoPads () const

std::set< int >	stripNumbersInDetId (unsigned int) const

std::set< unsigned int >	superChamberIds () const

std::set< unsigned int >	superChamberIdsWithCoPads () const

std::set< unsigned int >	superChamberIdsWithPads () const

	~GEMDigiMatcher ()

Private Types
using	DigiContainer = matching::DigiContainer

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

void	matchCoPadsToSimTrack (const GEMCoPadDigiCollection &co_pads)

void	matchDigisToSimTrack (const GEMDigiCollection &digis)

void	matchPadsToSimTrack (const GEMPadDigiCollection &pads)

Private Attributes
std::map< unsigned int, DigiContainer >	chamber_to_copads_

std::map< unsigned int, DigiContainer >	chamber_to_digis_

std::map< unsigned int, DigiContainer >	chamber_to_pads_

std::map< unsigned int, DigiContainer >	detid_to_copads_

std::map< unsigned int, DigiContainer >	detid_to_digis_

std::map< unsigned int, DigiContainer >	detid_to_pads_

edm::Handle< GEMCoPadDigiCollection >	gem_co_pads_

edm::Handle< GEMDigiCollection >	gem_digis_

const GEMGeometry &	gem_geo_

edm::Handle< GEMPadDigiCollection >	gem_pads_

int	matchDeltaStrip_

int	maxBXGEM_

int	minBXGEM_

const DigiContainer	no_digis_

const SimHitMatcher &	simhit_matcher_

std::map< unsigned int, DigiContainer >	superchamber_to_copads_

std::map< unsigned int, DigiContainer >	superchamber_to_digis_

std::map< unsigned int, DigiContainer >	superchamber_to_pads_

bool	verbose_

Detailed Description

Definition at line 34 of file GEMDigiMatcher.h.

Member Typedef Documentation

using GEMDigiMatcher::DigiContainer = matching::DigiContainer

private

Definition at line 35 of file GEMDigiMatcher.h.

Constructor & Destructor Documentation

GEMDigiMatcher::GEMDigiMatcher	(	const SimHitMatcher &	sh,
		const edm::Event &	e,
		const GEMGeometry &	geom,
		const edm::ParameterSet &	cfg,
		edm::EDGetToken &	gem_digiToken,
		edm::EDGetToken &	gem_padToken,
		edm::EDGetToken &	gem_copadToken
	)

Definition at line 10 of file GEMDigiMatcher.cc.

References gem_co_pads_, gem_digis_, gem_pads_, edm::Event::getByToken(), edm::ParameterSet::getUntrackedParameter(), init(), edm::HandleBase::isValid(), 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_digiToken, gem_digis_);
   e.getByToken(gem_padToken, gem_pads_);
   e.getByToken(gem_copadToken, gem_co_pads_);
   if (!gem_digis_.isValid())
     return;
 
   init(e);
 }

GEMDigiMatcher::~GEMDigiMatcher ( )

Definition at line 35 of file GEMDigiMatcher.cc.

35 {}

Member Function Documentation

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

Definition at line 160 of file GEMDigiMatcher.cc.

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

Referenced by GEMDigiTrackMatch::analyze().

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

std::set< unsigned int > GEMDigiMatcher::chamberIdsWithPads ( ) const

Definition at line 166 of file GEMDigiMatcher.cc.

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

Referenced by GEMDigiTrackMatch::analyze().

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

std::set< int > GEMDigiMatcher::coPadNumbersInDetId ( unsigned int detid ) const

Definition at line 301 of file GEMDigiMatcher.cc.

References coPadsInDetId(), edmIntegrityCheck::d, matching::digi_channel(), and mps_fire::result.

                                                                         {
   set<int> result;
   auto digis = coPadsInDetId(detid);
   for (auto &d : digis) {
     result.insert(digi_channel(d));
   }
   return result;
 }

const matching::DigiContainer & GEMDigiMatcher::coPadsInDetId ( unsigned int detid ) const

Definition at line 233 of file GEMDigiMatcher.cc.

References detid_to_copads_, and no_digis_.

Referenced by coPadNumbersInDetId().

                                                                                    {
   if (detid_to_copads_.find(detid) == detid_to_copads_.end())
     return no_digis_;
   return detid_to_copads_.at(detid);
 }

const matching::DigiContainer & GEMDigiMatcher::coPadsInSuperChamber ( unsigned int detid ) const

Definition at line 239 of file GEMDigiMatcher.cc.

References no_digis_, and superchamber_to_copads_.

Referenced by nCoPads().

                                                                                           {
   if (superchamber_to_copads_.find(detid) == superchamber_to_copads_.end())
     return no_digis_;
   return superchamber_to_copads_.at(detid);
 }

std::set< unsigned int > GEMDigiMatcher::detIds ( void ) const

Definition at line 153 of file GEMDigiMatcher.cc.

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

Referenced by partitionNumbers().

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

std::set< unsigned int > GEMDigiMatcher::detIdsWithCoPads ( ) const

Definition at line 190 of file GEMDigiMatcher.cc.

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

Referenced by partitionNumbersWithCoPads().

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

const matching::DigiContainer & GEMDigiMatcher::digisInChamber ( unsigned int detid ) const

Definition at line 203 of file GEMDigiMatcher.cc.

References chamber_to_digis_, and no_digis_.

                                                                                     {
   if (chamber_to_digis_.find(detid) == chamber_to_digis_.end())
     return no_digis_;
   return chamber_to_digis_.at(detid);
 }

const matching::DigiContainer & GEMDigiMatcher::digisInDetId ( unsigned int detid ) const

Definition at line 197 of file GEMDigiMatcher.cc.

References detid_to_digis_, and no_digis_.

Referenced by stripNumbersInDetId().

                                                                                   {
   if (detid_to_digis_.find(detid) == detid_to_digis_.end())
     return no_digis_;
   return detid_to_digis_.at(detid);
 }

const matching::DigiContainer & GEMDigiMatcher::digisInSuperChamber ( unsigned int detid ) const

Definition at line 209 of file GEMDigiMatcher.cc.

References no_digis_, and superchamber_to_digis_.

Referenced by nLayersWithDigisInSuperChamber().

                                                                                          {
   if (superchamber_to_digis_.find(detid) == superchamber_to_digis_.end())
     return no_digis_;
   return superchamber_to_digis_.at(detid);
 }

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

private

Definition at line 37 of file GEMDigiMatcher.cc.

References gem_co_pads_, gem_digis_, gem_pads_, edm::HandleBase::isValid(), matchCoPadsToSimTrack(), matchDigisToSimTrack(), matchPadsToSimTrack(), and edm::Handle< T >::product().

Referenced by GEMDigiMatcher().

                                            {
   matchDigisToSimTrack(*gem_digis_.product());
   if (gem_pads_.isValid())
     matchPadsToSimTrack(*gem_pads_.product());
   if (gem_co_pads_.isValid())
     matchCoPadsToSimTrack(*gem_co_pads_.product());
 }

void GEMDigiMatcher::matchCoPadsToSimTrack ( const GEMCoPadDigiCollection & co_pads )

private

Definition at line 122 of file GEMDigiMatcher.cc.

References GEMDetId::chamber(), detid_to_copads_, SimHitMatcher::detIdsGEMCoincidences(), matching::GEM_COPAD, SimHitMatcher::hitCoPadsInDetId(), triggerObjects_cff::id, matching::make_digi(), maxBXGEM_, minBXGEM_, GEMDetId::region(), GEMDetId::ring(), simhit_matcher_, GEMDetId::station(), and superchamber_to_copads_.

Referenced by init().

                                                                                 {
   auto det_ids = simhit_matcher_.detIdsGEMCoincidences();
   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_co_pads = simhit_matcher_.hitCoPadsInDetId(id);
     auto co_pads_in_det = co_pads.get(p_id);
 
     for (auto pad = co_pads_in_det.first; pad != co_pads_in_det.second; ++pad) {
       // check that the pad BX is within the range
       if (pad->bx(1) < minBXGEM_ || pad->bx(1) > maxBXGEM_)
         continue;
       if (pad->bx(2) < minBXGEM_ || pad->bx(2) > maxBXGEM_)
         continue;
       // check that it matches a coincidence pad that was hit by SimHits from
       // our track
       if (hit_co_pads.find(pad->pad(1)) == hit_co_pads.end())
         continue;
       if (hit_co_pads.find(pad->pad(2)) == hit_co_pads.end())
         continue;
 
       auto mydigi = make_digi(id, pad->pad(1), pad->bx(1), GEM_COPAD);
       // auto mydigi = make_digi(id, pad->pad(2), pad->bx(2), GEM_COPAD);  //
       // FIXIT : Solve duplicate problem.
       detid_to_copads_[id].push_back(mydigi);
       superchamber_to_copads_[superch_id()].push_back(mydigi);
     }
   }
 }

void GEMDigiMatcher::matchDigisToSimTrack ( const GEMDigiCollection & digis )

private

Definition at line 45 of file GEMDigiMatcher.cc.

References GEMDetId::chamber(), chamber_to_digis_, GEMDetId::chamberId(), popcon2dropbox::copy(), gather_cfg::cout, edmIntegrityCheck::d, detid_to_digis_, SimHitMatcher::detIdsGEM(), matching::GEM_STRIP, SimHitMatcher::hitStripsInDetId(), triggerObjects_cff::id, matching::make_digi(), matchDeltaStrip_, maxBXGEM_, minBXGEM_, DetId::rawId(), GEMDetId::region(), GEMDetId::ring(), simhit_matcher_, GEMDetId::station(), superchamber_to_digis_, 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 digis_in_det = digis.get(GEMDetId(id));
 
     for (auto d = digis_in_det.first; d != digis_in_det.second; ++d) {
       if (verbose_)
         cout << "gdigi " << p_id << " " << *d << endl;
       // check that the digi is within BX range
       if (d->bx() < minBXGEM_ || d->bx() > maxBXGEM_)
         continue;
       // check that it matches a strip that was hit by SimHits from our track
       if (hit_strips.find(d->strip()) == hit_strips.end())
         continue;
       if (verbose_)
         cout << "oki" << endl;
 
       auto mydigi = make_digi(id, d->strip(), d->bx(), GEM_STRIP);
       detid_to_digis_[id].push_back(mydigi);
       chamber_to_digis_[p_id.chamberId().rawId()].push_back(mydigi);
       superchamber_to_digis_[superch_id()].push_back(mydigi);
 
       // int pad_num = 1 + static_cast<int>( roll->padOfStrip(d->strip()) ); //
       // d->strip() is int digi_map[ make_pair(pad_num, d->bx()) ].push_back(
       // d->strip() );
     }
   }
 }

void GEMDigiMatcher::matchPadsToSimTrack ( const GEMPadDigiCollection & pads )

private

Definition at line 84 of file GEMDigiMatcher.cc.

References GEMDetId::chamber(), chamber_to_pads_, GEMDetId::chamberId(), popcon2dropbox::copy(), gather_cfg::cout, detid_to_pads_, SimHitMatcher::detIdsGEM(), SoftLeptonByDistance_cfi::distance, matching::GEM_PAD, SimHitMatcher::hitPadsInDetId(), triggerObjects_cff::id, matching::make_digi(), maxBXGEM_, minBXGEM_, DetId::rawId(), GEMDetId::region(), GEMDetId::ring(), simhit_matcher_, GEMDetId::station(), superchamber_to_pads_, 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_pads = simhit_matcher_.hitPadsInDetId(id);
     auto pads_in_det = pads.get(p_id);
 
     if (verbose_) {
       cout << "checkpads " << hit_pads.size() << " " << std::distance(pads_in_det.first, pads_in_det.second)
            << " hit_pads: ";
       copy(hit_pads.begin(), hit_pads.end(), ostream_iterator<int>(cout, " "));
       cout << endl;
     }
 
     for (auto pad = pads_in_det.first; pad != pads_in_det.second; ++pad) {
       if (verbose_)
         cout << "chp " << *pad << endl;
       // check that the pad BX is within the range
       if (pad->bx() < minBXGEM_ || pad->bx() > maxBXGEM_)
         continue;
       if (verbose_)
         cout << "chp1" << endl;
       // check that it matches a pad that was hit by SimHits from our track
       if (hit_pads.find(pad->pad()) == hit_pads.end())
         continue;
       if (verbose_)
         cout << "chp2" << endl;
 
       auto mydigi = make_digi(id, pad->pad(), pad->bx(), GEM_PAD);
       detid_to_pads_[id].push_back(mydigi);
       chamber_to_pads_[p_id.chamberId().rawId()].push_back(mydigi);
       superchamber_to_pads_[superch_id()].push_back(mydigi);
     }
   }
 }

int GEMDigiMatcher::nCoPads ( ) const

How many coincidence pads in GEM did this simtrack get in total?

Definition at line 265 of file GEMDigiMatcher.cc.

References coPadsInSuperChamber(), photons_cff::ids, gen::n, and superChamberIdsWithCoPads().

                                   {
   int n = 0;
   auto ids = superChamberIdsWithCoPads();
   for (auto id : ids) {
     n += coPadsInSuperChamber(id).size();
   }
   return n;
 }

int GEMDigiMatcher::nLayersWithDigisInSuperChamber ( unsigned int detid ) const

Definition at line 245 of file GEMDigiMatcher.cc.

References edmIntegrityCheck::d, matching::digi_id(), digisInSuperChamber(), GEMDetId::layer(), and LayerTriplets::layers().

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

int GEMDigiMatcher::nLayersWithPadsInSuperChamber ( unsigned int detid ) const

Definition at line 255 of file GEMDigiMatcher.cc.

References edmIntegrityCheck::d, matching::digi_id(), GEMDetId::layer(), LayerTriplets::layers(), and padsInSuperChamber().

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

int GEMDigiMatcher::nPads ( ) const

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

Definition at line 274 of file GEMDigiMatcher.cc.

References photons_cff::ids, gen::n, padsInSuperChamber(), and superChamberIdsWithPads().

                                 {
   int n = 0;
   auto ids = superChamberIdsWithPads();
   for (auto id : ids) {
     n += padsInSuperChamber(id).size();
   }
   return n;
 }

std::set< int > GEMDigiMatcher::padNumbersInDetId ( unsigned int detid ) const

Definition at line 292 of file GEMDigiMatcher.cc.

References edmIntegrityCheck::d, matching::digi_channel(), padsInDetId(), and mps_fire::result.

                                                                       {
   set<int> result;
   auto digis = padsInDetId(detid);
   for (auto &d : digis) {
     result.insert(digi_channel(d));
   }
   return result;
 }

const matching::DigiContainer & GEMDigiMatcher::padsInChamber ( unsigned int detid ) const

Definition at line 221 of file GEMDigiMatcher.cc.

References chamber_to_pads_, and no_digis_.

                                                                                    {
   if (chamber_to_pads_.find(detid) == chamber_to_pads_.end())
     return no_digis_;
   return chamber_to_pads_.at(detid);
 }

const matching::DigiContainer & GEMDigiMatcher::padsInDetId ( unsigned int detid ) const

Definition at line 215 of file GEMDigiMatcher.cc.

References detid_to_pads_, and no_digis_.

Referenced by padNumbersInDetId().

                                                                                  {
   if (detid_to_pads_.find(detid) == detid_to_pads_.end())
     return no_digis_;
   return detid_to_pads_.at(detid);
 }

const matching::DigiContainer & GEMDigiMatcher::padsInSuperChamber ( unsigned int detid ) const

Definition at line 227 of file GEMDigiMatcher.cc.

References no_digis_, and superchamber_to_pads_.

Referenced by nLayersWithPadsInSuperChamber(), and nPads().

                                                                                         {
   if (superchamber_to_pads_.find(detid) == superchamber_to_pads_.end())
     return no_digis_;
   return superchamber_to_pads_.at(detid);
 }

std::set< int > GEMDigiMatcher::partitionNumbers ( ) const

Definition at line 310 of file GEMDigiMatcher.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;
 }

std::set< int > GEMDigiMatcher::partitionNumbersWithCoPads ( ) const

Definition at line 321 of file GEMDigiMatcher.cc.

References detIdsWithCoPads(), mps_fire::result, and GEMDetId::roll().

                                                              {
   std::set<int> result;
 
   auto detids = detIdsWithCoPads();
   for (auto id : detids) {
     GEMDetId idd(id);
     result.insert(idd.roll());
   }
   return result;
 }

std::set< int > GEMDigiMatcher::stripNumbersInDetId ( unsigned int detid ) const

Definition at line 283 of file GEMDigiMatcher.cc.

References edmIntegrityCheck::d, matching::digi_channel(), digisInDetId(), and mps_fire::result.

                                                                         {
   set<int> result;
   auto digis = digisInDetId(detid);
   for (auto &d : digis) {
     result.insert(digi_channel(d));
   }
   return result;
 }

std::set< unsigned int > GEMDigiMatcher::superChamberIds ( ) const

Definition at line 172 of file GEMDigiMatcher.cc.

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

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

std::set< unsigned int > GEMDigiMatcher::superChamberIdsWithCoPads ( ) const

Definition at line 184 of file GEMDigiMatcher.cc.

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

Referenced by nCoPads().

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

std::set< unsigned int > GEMDigiMatcher::superChamberIdsWithPads ( ) const

Definition at line 178 of file GEMDigiMatcher.cc.

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

Referenced by nPads().

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