#include <GEMDigiTrackMatch.h>

Inheritance diagram for GEMDigiTrackMatch:

Public Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

	GEMDigiTrackMatch (const edm::ParameterSet &ps)

	~GEMDigiTrackMatch () override

Public Member Functions inherited from GEMTrackMatch
void	buildLUT (const int maxChamberId)

void	FillWithTrigger (MonitorElement *me[3], Float_t eta)

void	FillWithTrigger (MonitorElement *me[3][3], Float_t eta, Float_t phi, bool odd[3], bool even[3])

void	FillWithTrigger (MonitorElement *me[4][3], bool array[3][2], Float_t value)

void	FillWithTrigger (MonitorElement *me[4][3][3], bool array[3][2], Float_t eta, Float_t phi, bool odd[3], bool even[3])

	GEMTrackMatch (const edm::ParameterSet &cfg)

std::pair< int, int >	getClosestChambers (const int maxChamberId, int region, float phi)

std::pair< double, double >	getEtaRange (int station, int chamber)

std::pair< double, double >	getEtaRangeForPhi (int station)

bool	isSimTrackGood (const SimTrack &)

void	setGeometry (const GEMGeometry &geom)

	~GEMTrackMatch () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Private Attributes
MonitorElement *	copad_eta [4][3]

MonitorElement *	copad_phi [4][3][3]

MonitorElement *	dg_eta [4][3]

MonitorElement *	dg_phi [4][3][3]

MonitorElement *	dg_sh_eta [4][3]

MonitorElement *	dg_sh_phi [4][3][3]

edm::EDGetToken	gem_copadToken_

edm::EDGetToken	gem_digiToken_

edm::EDGetToken	gem_padToken_

MonitorElement *	pad_eta [4][3]

MonitorElement *	pad_phi [4][3][3]

MonitorElement *	track_eta [3]

MonitorElement *	track_phi [3][3]

Additional Inherited Members
Protected Attributes inherited from GEMTrackMatch
edm::ParameterSet	cfg_

float	chamberHeight_

bool	detailPlot_

std::vector< double >	etaRangeForPhi

const GEMGeometry *	gem_geom_

float	maxEta_

float	minEta_

float	minPt_

std::pair< std::vector< float >, std::vector< int > >	negativeLUT_

unsigned int	nstation

std::pair< std::vector< float >, std::vector< int > >	positiveLUT_

float	radiusCenter_

edm::EDGetToken	simHitsToken_

edm::EDGetToken	simTracksToken_

edm::EDGetToken	simVerticesToken_

int	useRoll_

Detailed Description

Definition at line 14 of file GEMDigiTrackMatch.h.

Constructor & Destructor Documentation

GEMDigiTrackMatch::GEMDigiTrackMatch ( const edm::ParameterSet & ps )

explicit

Definition at line 12 of file GEMDigiTrackMatch.cc.

References GEMTrackMatch::cfg_, GEMTrackMatch::detailPlot_, gem_copadToken_, gem_digiToken_, gem_padToken_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), GEMTrackMatch::simHitsToken_, GEMTrackMatch::simTracksToken_, GEMTrackMatch::simVerticesToken_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                               : GEMTrackMatch(ps) {
   std::string simInputLabel_ = ps.getUntrackedParameter<std::string>("simInputLabel");
   simHitsToken_ = consumes<edm::PSimHitContainer>(edm::InputTag(simInputLabel_, "MuonGEMHits"));
   simTracksToken_ = consumes<edm::SimTrackContainer>(ps.getParameter<edm::InputTag>("simTrackCollection"));
   simVerticesToken_ = consumes<edm::SimVertexContainer>(ps.getParameter<edm::InputTag>("simVertexCollection"));
 
   gem_digiToken_ = consumes<GEMDigiCollection>(ps.getParameter<edm::InputTag>("gemDigiInput"));
   gem_padToken_ = consumes<GEMPadDigiCollection>(ps.getParameter<edm::InputTag>("gemPadDigiInput"));
   gem_copadToken_ = consumes<GEMCoPadDigiCollection>(ps.getParameter<edm::InputTag>("gemCoPadDigiInput"));
   detailPlot_ = ps.getParameter<bool>("detailPlot");
   cfg_ = ps;
 }

GEMDigiTrackMatch::~GEMDigiTrackMatch ( )

override

Definition at line 91 of file GEMDigiTrackMatch.cc.

91 {}

Member Function Documentation

void GEMDigiTrackMatch::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Implements GEMTrackMatch.

Definition at line 93 of file GEMDigiTrackMatch.cc.

References begin, GEMTrackMatch::cfg_, relativeConstraints::chamber, GEMDigiMatcher::chamberIds(), SimHitMatcher::chamberIdsGEM(), GEMDigiMatcher::chamberIdsWithPads(), copad_eta, copad_phi, gather_cfg::cout, edmIntegrityCheck::d, GEMTrackMatch::detailPlot_, dg_eta, dg_phi, dg_sh_eta, dg_sh_phi, end, MySimTrack::eta, lumiContext::fill, GEMTrackMatch::FillWithTrigger(), gem_copadToken_, MySimTrack::gem_dg, gem_digiToken_, MySimTrack::gem_pad, gem_padToken_, MySimTrack::gem_sh, relativeConstraints::geom, edm::EventSetup::get(), edm::Event::getByToken(), MySimTrack::hitEven, MySimTrack::hitOdd, mps_fire::i, triggerObjects_cff::id, GEMTrackMatch::isSimTrackGood(), edm::HandleBase::isValid(), pad_eta, pad_phi, MySimTrack::phi, edm::Handle< T >::product(), MySimTrack::pt, GEMTrackMatch::simHitsToken_, GEMTrackMatch::simTracksToken_, GEMTrackMatch::simVerticesToken_, protons_cff::t, track_eta, and track_phi.

                                                                                    {
   edm::ESHandle<GEMGeometry> hGeom;
   iSetup.get<MuonGeometryRecord>().get(hGeom);
   const GEMGeometry &geom = *hGeom;
 
   edm::Handle<edm::PSimHitContainer> simhits;
   edm::Handle<edm::SimTrackContainer> sim_tracks;
   edm::Handle<edm::SimVertexContainer> sim_vertices;
   iEvent.getByToken(simHitsToken_, simhits);
   iEvent.getByToken(simTracksToken_, sim_tracks);
   iEvent.getByToken(simVerticesToken_, sim_vertices);
   if (!simhits.isValid() || !sim_tracks.isValid() || !sim_vertices.isValid())
     return;
 
   if (detailPlot_) {
     edm::Handle<edm::PSimHitContainer> simhits;
     edm::Handle<edm::SimTrackContainer> sim_tracks;
     edm::Handle<edm::SimVertexContainer> sim_vertices;
     iEvent.getByToken(simHitsToken_, simhits);
     iEvent.getByToken(simTracksToken_, sim_tracks);
     iEvent.getByToken(simVerticesToken_, sim_vertices);
     if (!simhits.isValid() || !sim_tracks.isValid() || !sim_vertices.isValid())
       return;
 
     // const edm::SimVertexContainer & sim_vert = *sim_vertices.product();
     const edm::SimTrackContainer &sim_trks = *sim_tracks.product();
     MySimTrack track_;
     for (auto &t : sim_trks) {
       if (!isSimTrackGood(t)) {
         continue;
       }
 
       // match hits and digis to this SimTrack
 
       const SimHitMatcher &match_sh =
           SimHitMatcher(t, iEvent, geom, cfg_, simHitsToken_, simTracksToken_, simVerticesToken_);
       const GEMDigiMatcher &match_gd =
           GEMDigiMatcher(match_sh, iEvent, geom, cfg_, gem_digiToken_, gem_padToken_, gem_copadToken_);
 
       track_.pt = t.momentum().pt();
       track_.phi = t.momentum().phi();
       track_.eta = t.momentum().eta();
       std::fill(std::begin(track_.hitOdd), std::end(track_.hitOdd), false);
       std::fill(std::begin(track_.hitEven), std::end(track_.hitEven), false);
 
       for (int i = 0; i < 3; i++) {
         for (int j = 0; j < 2; j++) {
           track_.gem_sh[i][j] = false;
           track_.gem_dg[i][j] = false;
           track_.gem_pad[i][j] = false;
         }
       }
 
       // ** GEM SimHits ** //
       const auto gem_sh_ids_ch = match_sh.chamberIdsGEM();
       for (auto d : gem_sh_ids_ch) {
         const GEMDetId id(d);
         if (id.chamber() % 2 == 0)
           track_.hitEven[id.station() - 1] = true;
         else if (id.chamber() % 2 == 1)
           track_.hitOdd[id.station() - 1] = true;
         else {
           std::cout << "Error to get chamber id" << std::endl;
         }
 
         track_.gem_sh[id.station() - 1][(id.layer() - 1)] = true;
       }
       // ** GEM Digis, Pads and CoPads ** //
       const auto gem_dg_ids_ch = match_gd.chamberIds();
 
       for (auto d : gem_dg_ids_ch) {
         const GEMDetId id(d);
         track_.gem_dg[id.station() - 1][(id.layer() - 1)] = true;
       }
       const auto gem_pad_ids_ch = match_gd.chamberIdsWithPads();
       for (auto d : gem_pad_ids_ch) {
         const GEMDetId id(d);
         track_.gem_pad[id.station() - 1][(id.layer() - 1)] = true;
       }
 
       FillWithTrigger(track_eta, fabs(track_.eta));
       FillWithTrigger(track_phi, fabs(track_.eta), track_.phi, track_.hitOdd, track_.hitEven);
 
       FillWithTrigger(dg_sh_eta, track_.gem_sh, fabs(track_.eta));
       FillWithTrigger(dg_eta, track_.gem_dg, fabs(track_.eta));
       FillWithTrigger(pad_eta, track_.gem_pad, fabs(track_.eta));
       FillWithTrigger(copad_eta, track_.gem_pad, fabs(track_.eta));
 
       // Separate station.
 
       FillWithTrigger(dg_sh_phi, track_.gem_sh, fabs(track_.eta), track_.phi, track_.hitOdd, track_.hitEven);
       FillWithTrigger(dg_phi, track_.gem_dg, fabs(track_.eta), track_.phi, track_.hitOdd, track_.hitEven);
       FillWithTrigger(pad_phi, track_.gem_pad, fabs(track_.eta), track_.phi, track_.hitOdd, track_.hitEven);
       FillWithTrigger(copad_phi, track_.gem_pad, fabs(track_.eta), track_.phi, track_.hitOdd, track_.hitEven);
     }
   }
 }

void GEMDigiTrackMatch::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	run,
		edm::EventSetup const &	iSetup
	)

override

Definition at line 25 of file GEMDigiTrackMatch.cc.

References DQMStore::IBooker::book1D(), GEMDetLabel::c_suffix, copad_eta, copad_phi, GEMTrackMatch::detailPlot_, dg_eta, dg_phi, dg_sh_eta, dg_sh_phi, relativeConstraints::geom, edm::EventSetup::get(), mps_fire::i, gen::k, GEMDetLabel::l_suffix, LogDebug, GEMTrackMatch::maxEta_, GEMTrackMatch::minEta_, GEMTrackMatch::nstation, pad_eta, pad_phi, Pi, PI, GEMGeometry::regions(), GEMDetLabel::s_suffix, DQMStore::IBooker::setCurrentFolder(), GEMTrackMatch::setGeometry(), AlCaHLTBitMon_QueryRunRegistry::string, createPayload::suffix, track_eta, and track_phi.

                                                                                                                {
   edm::ESHandle<GEMGeometry> hGeom;
   iSetup.get<MuonGeometryRecord>().get(hGeom);
   const GEMGeometry &geom = *hGeom;
   setGeometry(geom);
 
   ibooker.setCurrentFolder("MuonGEMDigisV/GEMDigisTask");
   LogDebug("GEMDigiTrackMatch") << "ibooker set current folder\n";
 
   const float PI = TMath::Pi();
 
   nstation = geom.regions()[0]->stations().size();
   if (detailPlot_) {
     for (unsigned int j = 0; j < nstation; j++) {
       string track_eta_name = string("track_eta") + s_suffix[j];
       string track_eta_title = string("track_eta") + ";SimTrack |#eta|;# of tracks";
       track_eta[j] = ibooker.book1D(track_eta_name.c_str(), track_eta_title.c_str(), 140, minEta_, maxEta_);
 
       for (unsigned int k = 0; k < 3; k++) {
         string suffix = string(s_suffix[j]) + string(c_suffix[k]);
         string track_phi_name = string("track_phi") + suffix;
         string track_phi_title = string("track_phi") + suffix + ";SimTrack #phi;# of tracks";
         track_phi[j][k] = ibooker.book1D(track_phi_name.c_str(), track_phi_title.c_str(), 200, -PI, PI);
       }
 
       for (unsigned int i = 0; i < 4; i++) {
         string suffix = string(s_suffix[j]) + string(l_suffix[i]);
         string dg_eta_name = string("dg_eta") + suffix;
         string dg_eta_title = dg_eta_name + "; tracks |#eta|; # of tracks";
         dg_eta[i][j] = ibooker.book1D(dg_eta_name.c_str(), dg_eta_title.c_str(), 140, minEta_, maxEta_);
 
         string dg_sh_eta_name = string("dg_sh_eta") + suffix;
         string dg_sh_eta_title = dg_sh_eta_name + "; tracks |#eta|; # of tracks";
         dg_sh_eta[i][j] = ibooker.book1D(dg_sh_eta_name.c_str(), dg_sh_eta_title.c_str(), 140, minEta_, maxEta_);
 
         string pad_eta_name = string("pad_eta") + suffix;
         string pad_eta_title = pad_eta_name + "; tracks |#eta|; # of tracks";
         pad_eta[i][j] = ibooker.book1D(pad_eta_name.c_str(), pad_eta_title.c_str(), 140, minEta_, maxEta_);
 
         string copad_eta_name = string("copad_eta") + suffix;
         string copad_eta_title = copad_eta_name + "; tracks |#eta|; # of tracks";
         copad_eta[i][j] = ibooker.book1D(copad_eta_name.c_str(), copad_eta_title.c_str(), 140, minEta_, maxEta_);
 
         for (unsigned int k = 0; k < 3; k++) {
           suffix = string(s_suffix[j]) + string(l_suffix[i]) + string(c_suffix[k]);
           string dg_phi_name = string("dg_phi") + suffix;
           string dg_phi_title = dg_phi_name + "; tracks #phi; # of tracks";
           dg_phi[i][j][k] = ibooker.book1D((dg_phi_name).c_str(), dg_phi_title.c_str(), 200, -PI, PI);
 
           string dg_sh_phi_name = string("dg_sh_phi") + suffix;
           string dg_sh_phi_title = dg_sh_phi_name + "; tracks #phi; # of tracks";
           dg_sh_phi[i][j][k] = ibooker.book1D((dg_sh_phi_name).c_str(), dg_sh_phi_title.c_str(), 200, -PI, PI);
 
           string pad_phi_name = string("pad_phi") + suffix;
           string pad_phi_title = pad_phi_name + "; tracks #phi; # of tracks";
           pad_phi[i][j][k] = ibooker.book1D((pad_phi_name).c_str(), pad_phi_title.c_str(), 200, -PI, PI);
 
           string copad_phi_name = string("copad_phi") + suffix;
           string copad_phi_title = copad_phi_name + "; tracks #phi; # of tracks";
           copad_phi[i][j][k] = ibooker.book1D((copad_phi_name).c_str(), copad_phi_title.c_str(), 200, -PI, PI);
         }
       }
     }
   }
 }