#include <GEMSimTrackMatch.h>

Inheritance diagram for GEMSimTrackMatch:

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

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

	GEMSimTrackMatch (const edm::ParameterSet &ps)

	~GEMSimTrackMatch () 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 DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
bool	detailPlot_

MonitorElement *	sh_eta [4][3]

MonitorElement *	sh_phi [4][3][3]

MonitorElement *	track_eta [3]

MonitorElement *	track_phi [3][3]

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
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_

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 7 of file GEMSimTrackMatch.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 12 of file GEMSimTrackMatch.cc.

References GEMTrackMatch::cfg_, detailPlot_, 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"));
   detailPlot_ = ps.getParameter<bool>("detailPlot");
 
   cfg_ = ps;
 }

GEMSimTrackMatch::~GEMSimTrackMatch ( )

override

Definition at line 74 of file GEMSimTrackMatch.cc.

74 {

75 }

Member Function Documentation

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

overridevirtual

Implements GEMTrackMatch.

Definition at line 77 of file GEMSimTrackMatch.cc.

References begin, GEMTrackMatch::cfg_, relativeConstraints::chamber, SimHitMatcher::chamberIdsGEM(), gather_cfg::cout, edmIntegrityCheck::d, detailPlot_, end, MySimTrack::eta, lumiContext::fill, GEMTrackMatch::FillWithTrigger(), 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(), MySimTrack::phi, edm::Handle< T >::product(), MySimTrack::pt, sh_eta, sh_phi, 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;
 
   MySimTrack track_;
   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();
 
 
   if ( detailPlot_) {
 
     for (auto& t: sim_trks)
     {
       if (!isSimTrackGood(t)) 
       { continue; } 
 
       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;
         }
       }
 
       // match hits to this SimTrack
       const SimHitMatcher match_sh = SimHitMatcher( t, iEvent, geom, cfg_, simHitsToken_, simTracksToken_, simVerticesToken_ );
 
       // check for hit chambers
       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;
       }
       FillWithTrigger( track_eta, fabs(track_.eta)) ;
       FillWithTrigger( track_phi, fabs(track_.eta), track_.phi, track_.hitOdd, track_.hitEven);
       FillWithTrigger( sh_eta, track_.gem_sh  , fabs( track_.eta) );
       FillWithTrigger( sh_phi, track_.gem_sh ,fabs(track_.eta), track_.phi , track_.hitOdd, track_.hitEven);
     }
   }
 }

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 24 of file GEMSimTrackMatch.cc.

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

 {
   // Mandantory
   edm::ESHandle<GEMGeometry> hGeom;
   iSetup.get<MuonGeometryRecord>().get(hGeom);
   const GEMGeometry& geom = *hGeom;
   setGeometry(geom);
     
   ibooker.setCurrentFolder("MuonGEMHitsV/GEMHitsTask");
   LogDebug("GEMSimTrackMatch")<<"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])+l_suffix[i];
 
         string sh_eta_name = string("sh_eta")+suffix;
         string sh_eta_title = sh_eta_name+"; tracks |#eta|; # of tracks";
         sh_eta[i][j] = ibooker.book1D( sh_eta_name.c_str(), sh_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 sh_phi_name = string("sh_phi")+suffix;
           string sh_phi_title = sh_phi_name+"; tracks #phi; # of tracks";
           sh_phi[i][j][k] = ibooker.book1D( (sh_phi_name).c_str(), sh_phi_title.c_str(), 200,-PI,PI) ;
 
         }
       }
     }
   }
 }