#include <GEMTrackMatch.h>

Inheritance diagram for GEMTrackMatch:

Public Member Functions
virtual void	analyze (const edm::Event &e, const edm::EventSetup &)=0

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	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &)=0

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 ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Protected Attributes
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_

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)

Detailed Description

Definition at line 42 of file GEMTrackMatch.h.

Constructor & Destructor Documentation

GEMTrackMatch::GEMTrackMatch ( const edm::ParameterSet & cfg )

explicit

Definition at line 6 of file GEMTrackMatch.cc.

References edm::ParameterSet::getUntrackedParameter(), maxEta_, minEta_, and minPt_.

 {
   minPt_  = ps.getUntrackedParameter<double>("gemMinPt",5.0);
   minEta_ = ps.getUntrackedParameter<double>("gemMinEta",1.55);
   maxEta_ = ps.getUntrackedParameter<double>("gemMaxEta",2.45);
 }

GEMTrackMatch::~GEMTrackMatch ( )

virtual

Definition at line 14 of file GEMTrackMatch.cc.

14 {

15 }

Member Function Documentation

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

pure virtual

Implements edm::stream::EDAnalyzerBase.

Implemented in GEMRecHitTrackMatch, GEMDigiTrackMatch, and GEMSimTrackMatch.

void GEMTrackMatch::buildLUT ( const int maxChamberId )

Definition at line 140 of file GEMTrackMatch.cc.

References i, negativeLUT_, positiveLUT_, and useRoll_.

Referenced by setGeometry().

 {
 
   edm::LogInfo("GEMTrackMatch")<<"max chamber "<<maxChamberId<<"\n";
 
   std::vector<int> pos_ids;
   pos_ids.push_back(GEMDetId(1,1,1,1,maxChamberId,useRoll_).rawId());
 
   std::vector<int> neg_ids;
   neg_ids.push_back(GEMDetId(-1,1,1,1,maxChamberId,useRoll_).rawId());
 
   // VK: I would really suggest getting phis from GEMGeometry
 
   std::vector<float> phis;
   phis.push_back(0.);
   for(int i=1; i<maxChamberId+1; ++i)
   {
     pos_ids.push_back(GEMDetId(1,1,1,1,i,useRoll_).rawId());
     neg_ids.push_back(GEMDetId(-1,1,1,1,i,useRoll_).rawId());
     phis.push_back(i*10.);
   }
   positiveLUT_ = std::make_pair(phis,pos_ids);
   negativeLUT_ = std::make_pair(phis,neg_ids);
 }

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

Definition at line 17 of file GEMTrackMatch.cc.

References MonitorElement::Fill(), i, and nstation.

Referenced by GEMSimTrackMatch::analyze(), GEMDigiTrackMatch::analyze(), and GEMRecHitTrackMatch::analyze().

 {
   for( unsigned int i=0 ; i<nstation ; i++) {
     hist[i]->Fill(eta);
   }
   return;
 }

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

Definition at line 25 of file GEMTrackMatch.cc.

References MonitorElement::Fill(), plotBeamSpotDB::first, getEtaRange(), i, nstation, edm::second(), and relativeConstraints::station.

 {
   for( unsigned int i=0 ; i<nstation ; i++) {
     int station = i+1;
     if ( odd[i] && eta > getEtaRange(station,1).first&& eta < getEtaRange(station,1).second ) {
       hist[i][0]->Fill(phi);
       if( hist[i][1] != nullptr) hist[i][1]->Fill(phi);
     }
     if ( even[i] && eta > getEtaRange(station,2).first&& eta < getEtaRange(station,2).second ) {
       hist[i][0]->Fill(phi);
       if( hist[i][2] != nullptr) hist[i][2]->Fill(phi);
     }
   }
   return;
 }

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

Definition at line 41 of file GEMTrackMatch.cc.

References MonitorElement::Fill(), i, and nstation.

 {
   for( unsigned int i=0 ; i<nstation ; i++) {
     if ( array[i][0] ) hist[0][i]->Fill(value);
     if ( array[i][1] ) hist[1][i]->Fill(value);
     if ( array[i][0] || array[i][1] ) hist[2][i]->Fill(value);
     if ( array[i][0] && array[i][1] ) hist[3][i]->Fill(value);
   } 
   return;
 }

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

Definition at line 52 of file GEMTrackMatch.cc.

References MonitorElement::Fill(), plotBeamSpotDB::first, getEtaRange(), i, nstation, edm::second(), and relativeConstraints::station.

 {
   for( unsigned int i=0 ; i<nstation ; i++) {
     int station = i+1;
     if ( odd[i] && eta > getEtaRange(station,1).first&& eta < getEtaRange(station,1).second ) {
       if ( array[i][0] ) { 
         hist[0][i][0]->Fill(phi);
         if ( hist[0][i][1] != nullptr) hist[0][i][1]->Fill(phi); 
       }
       if ( array[i][1] ) { 
         hist[1][i][0]->Fill(phi); 
         if( hist[1][i][1] != nullptr ) hist[1][i][1]->Fill(phi); 
       }
       if ( array[i][0] || array[i][1] ) { 
         hist[2][i][0]->Fill(phi); 
         if ( hist[2][i][1] != nullptr) hist[2][i][1]->Fill(phi); 
       }
       if ( array[i][0] && array[i][1] ) { 
         hist[3][i][0]->Fill(phi); 
         if ( hist[3][i][1] != nullptr) hist[3][i][1]->Fill(phi); 
       }
     }
     if ( even[i] && eta > getEtaRange(station,2).first&& eta < getEtaRange(station,2).second ) {
       if ( array[i][0] ) { 
         hist[0][i][0]->Fill(phi); 
         if ( hist[0][i][2]!=nullptr) hist[0][i][2]->Fill(phi); 
       }
       if ( array[i][1] ) { 
         hist[1][i][0]->Fill(phi); 
         if( hist[1][i][2]!=nullptr ) hist[1][i][2]->Fill(phi); 
       }
       if ( array[i][0] || array[i][1] ) { 
         hist[2][i][0]->Fill(phi); 
         if( hist[2][i][2]!=nullptr) hist[2][i][2]->Fill(phi); 
       }
       if ( array[i][0] && array[i][1] ) { 
         hist[3][i][0]->Fill(phi); 
         if( hist[3][i][2]!=nullptr) hist[3][i][2]->Fill(phi);
       }
     }
   }
   return;
 }

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

Definition at line 194 of file GEMTrackMatch.cc.

References negativeLUT_, phi, positiveLUT_, and pileupCalc::upper.

 {
   auto& phis(positiveLUT_.first);
   auto upper = std::upper_bound(phis.begin(), phis.end(), phi);
   auto& LUT = (region == 1 ? positiveLUT_.second : negativeLUT_.second);
   return std::make_pair(LUT.at(upper - phis.begin()), (LUT.at((upper - phis.begin() + 1)%maxChamberId)));
 }

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

Definition at line 96 of file GEMTrackMatch.cc.

References gather_cfg::cout, gem_geom_, and GEMGeometry::regions().

Referenced by FillWithTrigger().

 {
   if( gem_geom_ != nullptr) {
     auto& ch = gem_geom_->regions()[0]->stations()[station-1]->rings()[0]->superChambers()[chamber-1]->chambers()[0];
     auto& roll1 = ch->etaPartitions()[0]; //.begin();
     auto& roll2 = ch->etaPartitions()[ch->nEtaPartitions()-1];
     const BoundPlane& bSurface1(roll1->surface());
     const BoundPlane& bSurface2(roll2->surface());
     auto& parameters1( roll1->specs()->parameters());
     float height1(parameters1[2]);
     auto& parameters2( roll2->specs()->parameters());
     float height2(parameters2[2]);
     LocalPoint lTop1( 0., height1, 0.);
     GlobalPoint gTop1(bSurface1.toGlobal(lTop1));
     //LocalPoint lBottom1( 0., -height1, 0.);
     //GlobalPoint gBottom1(bSurface1.toGlobal(lBottom1));
     //LocalPoint lTop2( 0., height2, 0.);
     //GlobalPoint gTop2(bSurface2.toGlobal(lTop2));
     LocalPoint lBottom2( 0., -height2, 0.);
     GlobalPoint gBottom2(bSurface2.toGlobal(lBottom2));
     double eta1 = fabs(gTop1.eta()) - 0.01;
     double eta2 = fabs(gBottom2.eta()) + 0.01;
     return std::make_pair(eta1,eta2);
   }
   else { std::cout<<"Failed to get geometry information"<<std::endl;
     return std::make_pair(0,0);
   }
 }

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

Definition at line 202 of file GEMTrackMatch.cc.

References etaRangeForPhi.

 {
   std::pair<double, double> range;
   if( station== 0 )      range = std::make_pair( etaRangeForPhi[0],etaRangeForPhi[1]) ; 
   else if( station== 1 ) range = std::make_pair( etaRangeForPhi[2],etaRangeForPhi[3]) ; 
   else if( station== 2 ) range = std::make_pair( etaRangeForPhi[4],etaRangeForPhi[5]) ; 
 
   return range;
 }

bool GEMTrackMatch::isSimTrackGood ( const SimTrack & t )

Definition at line 126 of file GEMTrackMatch.cc.

References funct::abs(), eta, maxEta_, minEta_, minPt_, CoreSimTrack::momentum(), SimTrack::noGenpart(), SimTrack::noVertex(), and CoreSimTrack::type().

Referenced by GEMSimTrackMatch::analyze(), GEMDigiTrackMatch::analyze(), and GEMRecHitTrackMatch::analyze().

 {
 
   // SimTrack selection
   if (t.noVertex())   return false; 
   if (t.noGenpart()) return false;
   if (std::abs(t.type()) != 13) return false; // only interested in direct muon simtracks
   if (t.momentum().pt() < minPt_ ) return false;
   const float eta(std::abs(t.momentum().eta()));
   if (eta > maxEta_ || eta < minEta_ ) return false; // no GEMs could be in such eta
   return true;
 }

void GEMTrackMatch::setGeometry ( const GEMGeometry & geom )

Definition at line 166 of file GEMTrackMatch.cc.

References buildLUT(), GEMGeometry::chamber(), chamberHeight_, GEMGeometry::etaPartition(), gem_geom_, relativeConstraints::geom, GEMGeometry::idToDetUnit(), GEMChamber::nEtaPartitions(), PV3DBase< T, PVType, FrameType >::perp(), radiusCenter_, GEMGeometry::regions(), GEMEtaPartitionSpecs::specificTopology(), GEMEtaPartition::specs(), StripTopology::stripLength(), and useRoll_.

Referenced by GEMSimTrackMatch::bookHistograms(), GEMDigiTrackMatch::bookHistograms(), and GEMRecHitTrackMatch::bookHistograms().

 {
   gem_geom_ = &geom;
   bool isEvenOK = true;
   bool isOddOK  = true;
   useRoll_=1;
   if ( geom.etaPartition( GEMDetId(1,1,1,1,1,1) ) == nullptr) isOddOK = false;
   if ( geom.etaPartition( GEMDetId(1,1,1,1,2,1) ) == nullptr) isEvenOK = false;
   if ( !isEvenOK || !isOddOK) useRoll_=2;
 
   const auto top_chamber = static_cast<const GEMEtaPartition*>(geom.idToDetUnit(GEMDetId(1,1,1,1,1,useRoll_))); 
   const int nEtaPartitions(geom.chamber(GEMDetId(1,1,1,1,1,useRoll_))->nEtaPartitions());
   const auto bottom_chamber = static_cast<const GEMEtaPartition*>(geom.idToDetUnit(GEMDetId(1,1,1,1,1,nEtaPartitions)));
   const float top_half_striplength = top_chamber->specs()->specificTopology().stripLength()/2.;
   const float bottom_half_striplength = bottom_chamber->specs()->specificTopology().stripLength()/2.;
   const LocalPoint lp_top(0., top_half_striplength, 0.);
   const LocalPoint lp_bottom(0., -bottom_half_striplength, 0.);
   const GlobalPoint gp_top = top_chamber->toGlobal(lp_top);
   const GlobalPoint gp_bottom = bottom_chamber->toGlobal(lp_bottom);
 
   radiusCenter_ = (gp_bottom.perp() + gp_top.perp())/2.;
   chamberHeight_ = gp_top.perp() - gp_bottom.perp();
 
   const int maxChamberId = geom.regions()[0]->stations()[0]->superChambers().size();
   buildLUT(maxChamberId);
 }  