#include <AngleConverter.h>

Public Member Functions
	AngleConverter ()

void	checkAndUpdateGeometry (const edm::EventSetup &)
	Update the Geometry with current Event Setup. More...

int	getGlobalEta (unsigned int rawid, const L1MuDTChambPhDigi &aDigi, const L1MuDTChambThContainer *dtThDigis)
	Convert local eta coordinate to global digital microGMT scale. More...

int	getGlobalEta (unsigned int rawid, const CSCCorrelatedLCTDigi &aDigi)
	Convert local eta coordinate to global digital microGMT scale. More...

int	getGlobalEta (unsigned int rawid, const RPCDigi &aDigi)
	Convert local eta coordinate to global digital microGMT scale. More...

int	getGlobalPhi (unsigned int rawid, const L1MuDTChambPhDigi &aDigi)
	Convert local phi coordinate to global digital OMTF scale. More...

int	getGlobalPhi (unsigned int rawid, const CSCCorrelatedLCTDigi &aDigi)
	Convert local phi coordinate to global digital OMTF scale. More...

float	getGlobalPhi (unsigned int rawid, const RPCDigi &aDigi)

	~AngleConverter ()

Private Member Functions
const int	findBTIgroup (const L1MuDTChambPhDigi &aDigi, const L1MuDTChambThContainer *dtThDigis)
	Find BTI group. More...

bool	isCSCCounterClockwise (const std::unique_ptr< const CSCLayer > &layer) const
	Check orientation of strips in given CSC chamber. More...

Private Attributes
edm::ESHandle< CSCGeometry >	_geocsc

edm::ESHandle< DTGeometry >	_geodt

unsigned long long	_geom_cache_id

edm::ESHandle< RPCGeometry >	_georpc

Detailed Description

Definition at line 22 of file AngleConverter.h.

Constructor & Destructor Documentation

AngleConverter::AngleConverter ( )

Definition at line 23 of file AngleConverter.cc.

23 : _geom_cache_id(0ULL) { }

AngleConverter::_geom_cache_id

unsigned long long _geom_cache_id

Definition: AngleConverter.h:62

AngleConverter::~AngleConverter ( )

Definition at line 26 of file AngleConverter.cc.

26 { }

Member Function Documentation

void AngleConverter::checkAndUpdateGeometry ( const edm::EventSetup & es )

Update the Geometry with current Event Setup.

Definition at line 29 of file AngleConverter.cc.

References _geocsc, _geodt, _geom_cache_id, _georpc, edm::eventsetup::EventSetupRecord::cacheIdentifier(), relativeConstraints::geom, edm::EventSetup::get(), and edm::eventsetup::EventSetupRecord::get().

                                                                    {
   const MuonGeometryRecord& geom = es.get<MuonGeometryRecord>();
   unsigned long long geomid = geom.cacheIdentifier();
   if( _geom_cache_id != geomid ) {
     geom.get(_georpc);  
     geom.get(_geocsc);    
     geom.get(_geodt);
     _geom_cache_id = geomid;
   }  
 }

const int AngleConverter::findBTIgroup	(	const L1MuDTChambPhDigi &	aDigi,
		const L1MuDTChambThContainer *	dtThDigis
	)

private

Find BTI group.

If there are more than one theta digi we do not take is due to unresolvet ambiguity. In this case we take eta of the middle of the chamber.

Definition at line 276 of file AngleConverter.cc.

References L1MuDTChambPhDigi::bxNum(), L1MuDTChambThContainer::chThetaSegm(), i, L1MuDTChambThDigi::position(), L1MuDTChambPhDigi::scNum(), L1MuDTChambPhDigi::stNum(), and L1MuDTChambPhDigi::whNum().

Referenced by getGlobalEta().

                                                                {
 
   int bti_group = -1;
   
   const L1MuDTChambThDigi *theta_segm = dtThDigis->chThetaSegm(aDigi.whNum(),
                                    aDigi.stNum(),
                                    aDigi.scNum(),
                                    aDigi.bxNum());
   if(!theta_segm) return  bti_group;
   
   for(unsigned int i = 0; i < 7; ++i ){
     if(theta_segm->position(i) && bti_group<0) bti_group = i;
     else if(theta_segm->position(i) && bti_group>-1) return -1;
   }
       
   return bti_group;
 }

int AngleConverter::getGlobalEta	(	unsigned int	rawid,
		const L1MuDTChambPhDigi &	aDigi,
		const L1MuDTChambThContainer *	dtThDigis
	)

Convert local eta coordinate to global digital microGMT scale.

Definition at line 145 of file AngleConverter.cc.

References _geodt, PV3DBase< T, PVType, FrameType >::eta(), findBTIgroup(), L1MuDTChambPhDigi::scNum(), L1MuDTChambPhDigi::stNum(), and L1MuDTChambPhDigi::whNum().

                                                          {
 
   
   const DTChamberId baseid(aDigi.whNum(),aDigi.stNum(),aDigi.scNum()+1);
   
   // do not use this pointer for anything other than creating a trig geom
   std::unique_ptr<DTChamber> chamb(const_cast<DTChamber*>(_geodt->chamber(baseid)));
   
   std::unique_ptr<DTTrigGeom> trig_geom( new DTTrigGeom(chamb.get(),false) );
   chamb.release(); // release it here so no one gets funny ideas  
   // super layer one is the theta superlayer in a DT chamber
   // station 4 does not have a theta super layer
   // the BTI index from the theta trigger is an OR of some BTI outputs
   // so, we choose the BTI that's in the middle of the group
   // as the BTI that we get theta from
   // TODO:::::>>> need to make sure this ordering doesn't flip under wheel sign
   const int NBTI_theta = ( (baseid.station() != 4) ? 
                trig_geom->nCell(2) : trig_geom->nCell(3) );
   const int bti_group = findBTIgroup(aDigi,dtThDigis);
   const unsigned bti_actual = bti_group*NBTI_theta/7 + NBTI_theta/14 + 1;  
   DTBtiId thetaBTI;  
   if ( baseid.station() != 4 && bti_group != -1) {
     thetaBTI = DTBtiId(baseid,2,bti_actual);
   } else {
     // since this is phi oriented it'll give us theta in the middle
     // of the chamber
     thetaBTI = DTBtiId(baseid,3,1); 
   }
   const GlobalPoint theta_gp = trig_geom->CMSPosition(thetaBTI);
   int iEta = theta_gp.eta()/2.61*240;
   return iEta;
 }

int AngleConverter::getGlobalEta	(	unsigned int	rawid,
		const CSCCorrelatedLCTDigi &	aDigi
	)

Convert local eta coordinate to global digital microGMT scale.

Code taken from GeometryTranslator.

Will be replaced by direct CSC phi local to global scale transformation as used in FPGA implementation

Definition at line 181 of file AngleConverter.cc.

References _geocsc, PV3DBase< T, PVType, FrameType >::eta(), CSCPatternLUT::get2007Position(), CSCCorrelatedLCTDigi::getKeyWG(), CSCCorrelatedLCTDigi::getPattern(), CSCCorrelatedLCTDigi::getStrip(), isCSCCounterClockwise(), CSCConstants::KEY_ALCT_LAYER, PV3DBase< T, PVType, FrameType >::mag(), hltrates_dqm_sourceclient-live_cfg::offset, listBenchmarks::pattern, PV3DBase< T, PVType, FrameType >::phi(), PV3DBase< T, PVType, FrameType >::theta(), and Geom::Phi< T >::value().

                                                                                      {
 
 
   
   // alot of this is transcription and consolidation of the CSC
   // global phi calculation code
   // this works directly with the geometry 
   // rather than using the old phi luts
   const CSCDetId id(rawid); 
   // we should change this to weak_ptrs at some point
   // requires introducing std::shared_ptrs to geometry
   std::unique_ptr<const CSCChamber> chamb(_geocsc->chamber(id));
   std::unique_ptr<const CSCLayerGeometry> layer_geom(
                              chamb->layer(CSCConstants::KEY_ALCT_LAYER)->geometry()
                              );
   std::unique_ptr<const CSCLayer> layer(
                     chamb->layer(CSCConstants::KEY_ALCT_LAYER)
                     );
   
   const uint16_t halfstrip = aDigi.getStrip();
   const uint16_t pattern = aDigi.getPattern();
   const uint16_t keyWG = aDigi.getKeyWG();
   //const unsigned maxStrips = layer_geom->numberOfStrips();  
 
   // so we can extend this later 
   // assume TMB2007 half-strips only as baseline
   double offset = 0.0;
   switch(1) {
   case 1:
     offset = CSCPatternLUT::get2007Position(pattern);
   }
   const unsigned halfstrip_offs = unsigned(0.5 + halfstrip + offset);
   const unsigned strip = halfstrip_offs/2 + 1; // geom starts from 1
 
   // the rough location of the hit at the ALCT key layer
   // we will refine this using the half strip information
   const LocalPoint coarse_lp = 
     layer_geom->stripWireGroupIntersection(strip,keyWG);  
   const GlobalPoint coarse_gp = layer->surface().toGlobal(coarse_lp);  
   
   // the strip width/4.0 gives the offset of the half-strip
   // center with respect to the strip center
   const double hs_offset = layer_geom->stripPhiPitch()/4.0;
   
   // determine handedness of the chamber
   const bool ccw = isCSCCounterClockwise(layer);
   // we need to subtract the offset of even half strips and add the odd ones
   const double phi_offset = ( ( halfstrip_offs%2 ? 1 : -1)*
                   ( ccw ? -hs_offset : hs_offset ) );
   
   // the global eta calculation uses the middle of the strip
   // so no need to increment it
   const GlobalPoint final_gp( GlobalPoint::Polar( coarse_gp.theta(),
                           (coarse_gp.phi().value() + 
                            phi_offset),
                           coarse_gp.mag() ) );  
   // release ownership of the pointers
   chamb.release();
   layer_geom.release();
   layer.release();
 
   int iEta =  final_gp.eta()/2.61*240;
   return iEta;  
 }

int AngleConverter::getGlobalEta	(	unsigned int	rawid,
		const RPCDigi &	aDigi
	)

Convert local eta coordinate to global digital microGMT scale.

Definition at line 250 of file AngleConverter.cc.

References _georpc, PV3DBase< T, PVType, FrameType >::eta(), and RPCDigi::strip().

                                                                         {
 
   const RPCDetId id(rawid);
   std::unique_ptr<const RPCRoll>  roll(_georpc->roll(id));
   const uint16_t strip = aDigi.strip();
   const LocalPoint lp = roll->centreOfStrip(strip);
   const GlobalPoint gp = roll->toGlobal(lp);
   roll.release();
 
   float iEta = gp.eta()/2.61*240;
 
   return iEta;
 
 } 

int AngleConverter::getGlobalPhi	(	unsigned int	rawid,
		const L1MuDTChambPhDigi &	aDigi
	)

Convert local phi coordinate to global digital OMTF scale.

Definition at line 60 of file AngleConverter.cc.

References M_PI, OMTFConfiguration::nPhiBins, phi, L1MuDTChambPhDigi::phi(), and L1MuDTChambPhDigi::scNum().

                                                                                    {
 
   // local phi in sector -> global phi
   float phi = aDigi.phi()/4096.0; 
   phi += aDigi.scNum()*M_PI/6.0; // add sector offset
   if(phi>M_PI) phi-=2*M_PI;
   phi/=2.0*M_PI;
       
   int iPhi = phi*OMTFConfiguration::nPhiBins;
    
   return iPhi;
 }

int AngleConverter::getGlobalPhi	(	unsigned int	rawid,
		const CSCCorrelatedLCTDigi &	aDigi
	)

Convert local phi coordinate to global digital OMTF scale.

Code taken from GeometryTranslator. Will be replaced by direct CSC phi local to global scale transformation as used in FPGA implementation

Definition at line 74 of file AngleConverter.cc.

References _geocsc, CSCPatternLUT::get2007Position(), CSCCorrelatedLCTDigi::getKeyWG(), CSCCorrelatedLCTDigi::getPattern(), CSCCorrelatedLCTDigi::getStrip(), isCSCCounterClockwise(), CSCConstants::KEY_ALCT_LAYER, M_PI, PV3DBase< T, PVType, FrameType >::mag(), OMTFConfiguration::nPhiBins, hltrates_dqm_sourceclient-live_cfg::offset, listBenchmarks::pattern, phi, PV3DBase< T, PVType, FrameType >::phi(), PV3DBase< T, PVType, FrameType >::theta(), and Geom::Phi< T >::value().

                                                                                       {
 
 
   
   // alot of this is transcription and consolidation of the CSC
   // global phi calculation code
   // this works directly with the geometry 
   // rather than using the old phi luts
   const CSCDetId id(rawid); 
   // we should change this to weak_ptrs at some point
   // requires introducing std::shared_ptrs to geometry
   std::unique_ptr<const CSCChamber> chamb(_geocsc->chamber(id));
   std::unique_ptr<const CSCLayerGeometry> layer_geom(
                              chamb->layer(CSCConstants::KEY_ALCT_LAYER)->geometry()
                              );
   std::unique_ptr<const CSCLayer> layer(
                     chamb->layer(CSCConstants::KEY_ALCT_LAYER)
                     );
   
   const uint16_t halfstrip = aDigi.getStrip();
   const uint16_t pattern = aDigi.getPattern();
   const uint16_t keyWG = aDigi.getKeyWG();
   //const unsigned maxStrips = layer_geom->numberOfStrips();  
 
   // so we can extend this later 
   // assume TMB2007 half-strips only as baseline
   double offset = 0.0;
   switch(1) {
   case 1:
     offset = CSCPatternLUT::get2007Position(pattern);
   }
   const unsigned halfstrip_offs = unsigned(0.5 + halfstrip + offset);
   const unsigned strip = halfstrip_offs/2 + 1; // geom starts from 1
 
   // the rough location of the hit at the ALCT key layer
   // we will refine this using the half strip information
   const LocalPoint coarse_lp = 
     layer_geom->stripWireGroupIntersection(strip,keyWG);  
   const GlobalPoint coarse_gp = layer->surface().toGlobal(coarse_lp);  
   
   // the strip width/4.0 gives the offset of the half-strip
   // center with respect to the strip center
   const double hs_offset = layer_geom->stripPhiPitch()/4.0;
   
   // determine handedness of the chamber
   const bool ccw = isCSCCounterClockwise(layer);
   // we need to subtract the offset of even half strips and add the odd ones
   const double phi_offset = ( ( halfstrip_offs%2 ? 1 : -1)*
                   ( ccw ? -hs_offset : hs_offset ) );
   
   // the global eta calculation uses the middle of the strip
   // so no need to increment it
   const GlobalPoint final_gp( GlobalPoint::Polar( coarse_gp.theta(),
                           (coarse_gp.phi().value() + 
                            phi_offset),
                           coarse_gp.mag() ) );  
   // release ownership of the pointers
   chamb.release();
   layer_geom.release();
   layer.release();
 
   float phi = final_gp.phi()/(2.0*M_PI);
   int iPhi =  phi*OMTFConfiguration::nPhiBins;
   
   return iPhi;
 }

float AngleConverter::getGlobalPhi	(	unsigned int	rawid,
		const RPCDigi &	aDigi
	)

Convert local phi coordinate to global digital OMTF scale. To maintain backward comtability return float value fo global phi. Later whewn LUT will be used will return int as other methods.

Will be replaced by LUT based transformation as used in FPGA implementation

Definition at line 41 of file AngleConverter.cc.

References _georpc, PV3DBase< T, PVType, FrameType >::phi(), and RPCDigi::strip().

                                                                            {
 
 
   const RPCDetId id(rawid);
   std::unique_ptr<const RPCRoll>  roll(_georpc->roll(id));
   const uint16_t strip = aDigi.strip();
   const LocalPoint lp = roll->centreOfStrip(strip);
   const GlobalPoint gp = roll->toGlobal(lp);
   roll.release();
 
   //float phi = gp.phi()/(2.0*M_PI);
   //int iPhi = phi*OMTFConfiguration::nPhiBins;
 
   return gp.phi();
 }

bool AngleConverter::isCSCCounterClockwise ( const std::unique_ptr< const CSCLayer > & layer ) const

private

Check orientation of strips in given CSC chamber.

Definition at line 267 of file AngleConverter.cc.

References funct::abs(), and M_PI.

Referenced by getGlobalEta(), and getGlobalPhi().

                                                                         {
   const int nStrips = layer->geometry()->numberOfStrips();
   const double phi1 = layer->centerOfStrip(1).phi();
   const double phiN = layer->centerOfStrip(nStrips).phi();
   return ( (std::abs(phi1 - phiN) < M_PI  && phi1 >= phiN) || 
        (std::abs(phi1 - phiN) >= M_PI && phi1 < phiN)     );  
 }