#include <GeometryTranslator.h>

Public Member Functions
double	calculateBendAngle (const TriggerPrimitive &) const

double	calculateGlobalEta (const TriggerPrimitive &) const

double	calculateGlobalPhi (const TriggerPrimitive &) const

void	checkAndUpdateGeometry (const edm::EventSetup &)

	GeometryTranslator ()

	~GeometryTranslator ()

Private Member Functions
double	calcCSCSpecificBend (const TriggerPrimitive &) const

double	calcCSCSpecificEta (const TriggerPrimitive &) const

double	calcCSCSpecificPhi (const TriggerPrimitive &) const

double	calcDTSpecificBend (const TriggerPrimitive &) const

double	calcDTSpecificEta (const TriggerPrimitive &) const

double	calcDTSpecificPhi (const TriggerPrimitive &) const

GlobalPoint	calcDTSpecificPoint (const TriggerPrimitive &) const

double	calcRPCSpecificBend (const TriggerPrimitive &) const

double	calcRPCSpecificEta (const TriggerPrimitive &) const

double	calcRPCSpecificPhi (const TriggerPrimitive &) const

GlobalPoint	getCSCSpecificPoint (const TriggerPrimitive &) const

GlobalPoint	getRPCSpecificPoint (const TriggerPrimitive &) const

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

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 36 of file GeometryTranslator.h.

Constructor & Destructor Documentation

GeometryTranslator::GeometryTranslator ( )

Definition at line 20 of file GeometryTranslator.cc.

                                       :  
   _geom_cache_id(0ULL) {
 }

GeometryTranslator::~GeometryTranslator ( )

Definition at line 24 of file GeometryTranslator.cc.

24 {

25 }

Member Function Documentation

double GeometryTranslator::calcCSCSpecificBend ( const TriggerPrimitive & tp ) const

private

Definition at line 201 of file GeometryTranslator.cc.

Referenced by calculateBendAngle().

                                                                         {  
   return 0.0;
 }

double GeometryTranslator::calcCSCSpecificEta ( const TriggerPrimitive & tp ) const

private

Definition at line 191 of file GeometryTranslator.cc.

References PV3DBase< T, PVType, FrameType >::eta(), and getCSCSpecificPoint().

Referenced by calculateGlobalEta().

                                                                        {  
   return getCSCSpecificPoint(tp).eta();
 }

double GeometryTranslator::calcCSCSpecificPhi ( const TriggerPrimitive & tp ) const

private

Definition at line 196 of file GeometryTranslator.cc.

References getCSCSpecificPoint(), and PV3DBase< T, PVType, FrameType >::phi().

Referenced by calculateGlobalPhi().

                                                                        {  
   return getCSCSpecificPoint(tp).phi();
 }

double GeometryTranslator::calcDTSpecificBend ( const TriggerPrimitive & tp ) const

private

Definition at line 255 of file GeometryTranslator.cc.

References L1TMuon::TriggerPrimitive::DTData::bendingAngle, and L1TMuon::TriggerPrimitive::getDTData().

Referenced by calculateBendAngle().

                                                                        {
   int bend = tp.getDTData().bendingAngle;
   double bendf = bend/512.0;
   return bendf;
 }

double GeometryTranslator::calcDTSpecificEta ( const TriggerPrimitive & tp ) const

private

Definition at line 244 of file GeometryTranslator.cc.

References calcDTSpecificPoint(), and PV3DBase< T, PVType, FrameType >::eta().

Referenced by calculateGlobalEta().

                                                                       {  
   return calcDTSpecificPoint(tp).eta();
 }

double GeometryTranslator::calcDTSpecificPhi ( const TriggerPrimitive & tp ) const

private

Definition at line 249 of file GeometryTranslator.cc.

References calcDTSpecificPoint(), and PV3DBase< T, PVType, FrameType >::phi().

Referenced by calculateGlobalPhi().

                                                                       {
   return calcDTSpecificPoint(tp).phi();
 }

GlobalPoint GeometryTranslator::calcDTSpecificPoint ( const TriggerPrimitive & tp ) const

private

Definition at line 206 of file GeometryTranslator.cc.

References _geodt, L1TMuon::TriggerPrimitive::detId(), L1TMuon::TriggerPrimitive::getDTData(), M_PI, PV3DBase< T, PVType, FrameType >::mag(), phi(), L1TMuon::TriggerPrimitive::DTData::radialAngle, L1TMuon::TriggerPrimitive::DTData::sector, DTChamberId::station(), PV3DBase< T, PVType, FrameType >::theta(), and L1TMuon::TriggerPrimitive::DTData::theta_bti_group.

Referenced by calcDTSpecificEta(), and calcDTSpecificPhi().

                                                                         {  
   const DTChamberId baseid(tp.detId<DTChamberId>());
   // 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 = tp.getDTData().theta_bti_group;
   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);
   
   // local phi in sector -> global phi
   double phi = ((double)tp.getDTData().radialAngle)/4096.0; 
   phi += tp.getDTData().sector*M_PI/6.0; // add sector offset  
 
   return GlobalPoint( GlobalPoint::Polar( theta_gp.theta(),
                       phi,
                       theta_gp.mag() ) );
 }

double GeometryTranslator::calcRPCSpecificBend ( const TriggerPrimitive & tp ) const

private

Definition at line 118 of file GeometryTranslator.cc.

Referenced by calculateBendAngle().

                                                                         {
   return 0.0;
 }

double GeometryTranslator::calcRPCSpecificEta ( const TriggerPrimitive & tp ) const

private

Definition at line 106 of file GeometryTranslator.cc.

References PV3DBase< T, PVType, FrameType >::eta(), and getRPCSpecificPoint().

Referenced by calculateGlobalEta().

                                                                        {  
   return getRPCSpecificPoint(tp).eta();
 }

double GeometryTranslator::calcRPCSpecificPhi ( const TriggerPrimitive & tp ) const

private

Definition at line 111 of file GeometryTranslator.cc.

References getRPCSpecificPoint(), and PV3DBase< T, PVType, FrameType >::phi().

Referenced by calculateGlobalPhi().

                                                                        {  
   return getRPCSpecificPoint(tp).phi();
 }

double GeometryTranslator::calculateBendAngle ( const TriggerPrimitive & tp ) const

Definition at line 64 of file GeometryTranslator.cc.

References calcCSCSpecificBend(), calcDTSpecificBend(), calcRPCSpecificBend(), L1TMuon::TriggerPrimitive::kCSC, L1TMuon::TriggerPrimitive::kDT, L1TMuon::TriggerPrimitive::kRPC, and L1TMuon::TriggerPrimitive::subsystem().

                                                                        {
   switch(tp.subsystem()) {
   case TriggerPrimitive::kDT:
     return calcDTSpecificBend(tp);
     break;
   case TriggerPrimitive::kCSC:
     return calcCSCSpecificBend(tp);
     break;
   case TriggerPrimitive::kRPC:
     return calcRPCSpecificBend(tp);
     break;
   default:
     return std::nan("Invalid TP type!");
     break;
   }
 }

double GeometryTranslator::calculateGlobalEta ( const TriggerPrimitive & tp ) const

Definition at line 28 of file GeometryTranslator.cc.

References calcCSCSpecificEta(), calcDTSpecificEta(), calcRPCSpecificEta(), L1TMuon::TriggerPrimitive::kCSC, L1TMuon::TriggerPrimitive::kDT, L1TMuon::TriggerPrimitive::kRPC, and L1TMuon::TriggerPrimitive::subsystem().

                                                                        {
   switch(tp.subsystem()) {
   case TriggerPrimitive::kDT:
     return calcDTSpecificEta(tp);
     break;
   case TriggerPrimitive::kCSC:
     return calcCSCSpecificEta(tp);
     break;
   case TriggerPrimitive::kRPC:
     return calcRPCSpecificEta(tp);
     break;
   default:
     return std::nan("Invalid TP type!"); 
     break;
   }
 }

double GeometryTranslator::calculateGlobalPhi ( const TriggerPrimitive & tp ) const

Definition at line 46 of file GeometryTranslator.cc.

References calcCSCSpecificPhi(), calcDTSpecificPhi(), calcRPCSpecificPhi(), L1TMuon::TriggerPrimitive::kCSC, L1TMuon::TriggerPrimitive::kDT, L1TMuon::TriggerPrimitive::kRPC, and L1TMuon::TriggerPrimitive::subsystem().

                                                                        {
   switch(tp.subsystem()) {
   case TriggerPrimitive::kDT:
     return calcDTSpecificPhi(tp);
     break;
   case TriggerPrimitive::kCSC:
     return calcCSCSpecificPhi(tp);
     break;
   case TriggerPrimitive::kRPC:
     return calcRPCSpecificPhi(tp);
     break;
   default:
     return std::nan("Invalid TP type!");
     break;
   }
 }

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

Definition at line 81 of file GeometryTranslator.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;
   }  
 }

GlobalPoint GeometryTranslator::getCSCSpecificPoint ( const TriggerPrimitive & tp ) const

private

Definition at line 128 of file GeometryTranslator.cc.

References _geocsc, L1TMuon::TriggerPrimitive::detId(), CSCPatternLUT::get2007Position(), L1TMuon::TriggerPrimitive::getCSCData(), isCSCCounterClockwise(), CSCConstants::KEY_ALCT_LAYER, L1TMuon::TriggerPrimitive::CSCData::keywire, PV3DBase< T, PVType, FrameType >::mag(), hltrates_dqm_sourceclient-live_cfg::offset, listBenchmarks::pattern, L1TMuon::TriggerPrimitive::CSCData::pattern, PV3DBase< T, PVType, FrameType >::phi(), L1TMuon::TriggerPrimitive::CSCData::strip, PV3DBase< T, PVType, FrameType >::theta(), and Geom::Phi< T >::value().

Referenced by calcCSCSpecificEta(), and calcCSCSpecificPhi().

                                                                         {
   const CSCDetId id(tp.detId<CSCDetId>()); 
   // 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 = tp.getCSCData().strip;
   const uint16_t pattern = tp.getCSCData().pattern;
   const uint16_t keyWG = tp.getCSCData().keywire; 
   //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() ) );
     
   // We need to add in some notion of the 'error' on trigger primitives
   // like the width of the wire group by the width of the strip
   // or something similar      
 
   // release ownership of the pointers
   chamb.release();
   layer_geom.release();
   layer.release();
   
   return final_gp;
 }

GlobalPoint GeometryTranslator::getRPCSpecificPoint ( const TriggerPrimitive & tp ) const

private

Definition at line 93 of file GeometryTranslator.cc.

References _georpc, RPCRoll::centreOfStrip(), L1TMuon::TriggerPrimitive::detId(), L1TMuon::TriggerPrimitive::getRPCData(), L1TMuon::TriggerPrimitive::RPCData::strip, and GeomDet::toGlobal().

Referenced by calcRPCSpecificEta(), and calcRPCSpecificPhi().

                                                                         {
   const RPCDetId id(tp.detId<RPCDetId>());
   const RPCRoll * roll = _georpc->roll(id);
   const uint16_t strip = tp.getRPCData().strip;
   const LocalPoint lp = roll->centreOfStrip(strip);
   const GlobalPoint gp = roll->toGlobal(lp);
   
   //roll.release();
  
   return gp;
 }

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

private

Definition at line 262 of file GeometryTranslator.cc.

References funct::abs(), and M_PI.

Referenced by getCSCSpecificPoint().

                                                                         {
   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)     );  
 }