GEMDigiTrackMatch Class Reference

#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 ()
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)
virtual	~GEMTrackMatch ()
Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final
virtual void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer (void)
virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
uint32_t	streamId () const
Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default
Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

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
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
static std::shared_ptr< dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
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 15 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

(

)

Definition at line 95 of file GEMDigiTrackMatch.cc.

{  }

Member Function Documentation

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

overridevirtual

Implements GEMTrackMatch.

Definition at line 97 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, hcalTTPDigis_cfi::id, GEMTrackMatch::isSimTrackGood(), edm::HandleBase::isValid(), pad_eta, pad_phi, MySimTrack::phi, edm::Handle< T >::product(), MySimTrack::pt, GEMTrackMatch::simHitsToken_, GEMTrackMatch::simTracksToken_, GEMTrackMatch::simVerticesToken_, lumiQTWidget::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  )

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 26 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) ;

        }
      }
    }
  }
}

Member Data Documentation

MonitorElement* GEMDigiTrackMatch::copad_eta[4][3]

private

Definition at line 30 of file GEMDigiTrackMatch.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* GEMDigiTrackMatch::copad_phi[4][3][3]

private

Definition at line 35 of file GEMDigiTrackMatch.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* GEMDigiTrackMatch::dg_eta[4][3]

private

Definition at line 27 of file GEMDigiTrackMatch.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* GEMDigiTrackMatch::dg_phi[4][3][3]

private

Definition at line 32 of file GEMDigiTrackMatch.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* GEMDigiTrackMatch::dg_sh_eta[4][3]

private

Definition at line 28 of file GEMDigiTrackMatch.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* GEMDigiTrackMatch::dg_sh_phi[4][3][3]

private

Definition at line 33 of file GEMDigiTrackMatch.h.

Referenced by analyze(), and bookHistograms().

edm::EDGetToken GEMDigiTrackMatch::gem_copadToken_

private

Definition at line 39 of file GEMDigiTrackMatch.h.

Referenced by analyze(), and GEMDigiTrackMatch().

edm::EDGetToken GEMDigiTrackMatch::gem_digiToken_

private

Definition at line 37 of file GEMDigiTrackMatch.h.

Referenced by analyze(), and GEMDigiTrackMatch().

edm::EDGetToken GEMDigiTrackMatch::gem_padToken_

private

Definition at line 38 of file GEMDigiTrackMatch.h.

Referenced by analyze(), and GEMDigiTrackMatch().

MonitorElement* GEMDigiTrackMatch::pad_eta[4][3]

private

Definition at line 29 of file GEMDigiTrackMatch.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* GEMDigiTrackMatch::pad_phi[4][3][3]

private

Definition at line 34 of file GEMDigiTrackMatch.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* GEMDigiTrackMatch::track_eta[3]

private

Definition at line 24 of file GEMDigiTrackMatch.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* GEMDigiTrackMatch::track_phi[3][3]

private

Definition at line 25 of file GEMDigiTrackMatch.h.

Referenced by analyze(), and bookHistograms().