L1TMuon::GeometryTranslator Class Reference

#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 ()
const CSCGeometry &	getCSCGeometry () const
const DTGeometry &	getDTGeometry () const
const GEMGeometry &	getGEMGeometry () const
GlobalPoint	getGlobalPoint (const TriggerPrimitive &) const
const MagneticField &	getMagneticField () const
const RPCGeometry &	getRPCGeometry () const
	~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	calcGEMSpecificBend (const TriggerPrimitive &) const
double	calcGEMSpecificEta (const TriggerPrimitive &) const
double	calcGEMSpecificPhi (const TriggerPrimitive &) const
double	calcRPCSpecificBend (const TriggerPrimitive &) const
double	calcRPCSpecificEta (const TriggerPrimitive &) const
double	calcRPCSpecificPhi (const TriggerPrimitive &) const
GlobalPoint	getCSCSpecificPoint (const TriggerPrimitive &) const
GlobalPoint	getGEMSpecificPoint (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
edm::ESHandle< GEMGeometry >	_geogem
unsigned long long	_geom_cache_id
edm::ESHandle< RPCGeometry >	_georpc
edm::ESHandle< MagneticField >	_magfield
unsigned long long	_magfield_cache_id

Detailed Description

Definition at line 38 of file GeometryTranslator.h.

Constructor & Destructor Documentation

GeometryTranslator::GeometryTranslator

(

)

Definition at line 24 of file GeometryTranslator.cc.

                                      :
  _geom_cache_id(0ULL), _magfield_cache_id(0ULL) {
}

GeometryTranslator::~GeometryTranslator

(

)

Definition at line 28 of file GeometryTranslator.cc.

                                        {  
}

Member Function Documentation

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

private

Definition at line 277 of file GeometryTranslator.cc.

References L1TMuon::TriggerPrimitive::CSCData::bend, and L1TMuon::TriggerPrimitive::getCSCData().

Referenced by calculateBendAngle().

                                                                        {
  return tp.getCSCData().bend;
}

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

private

Definition at line 267 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 272 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 331 of file GeometryTranslator.cc.

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

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 320 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 325 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 282 of file GeometryTranslator.cc.

References _geodt, DTGeometry::chamber(), L1TMuon::TriggerPrimitive::detId(), L1TMuon::TriggerPrimitive::getDTData(), M_PI, PV3DBase< T, PVType, FrameType >::mag(), 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(), calcDTSpecificPhi(), and getGlobalPoint().

                                                                        {  
  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::calcGEMSpecificBend ( const TriggerPrimitive &  tp ) const

private

Definition at line 161 of file GeometryTranslator.cc.

Referenced by calculateBendAngle().

                                                                        {
  return 0.0;
}

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

private

Definition at line 149 of file GeometryTranslator.cc.

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

Referenced by calculateGlobalEta().

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

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

private

Definition at line 154 of file GeometryTranslator.cc.

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

Referenced by calculateGlobalPhi().

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

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

private

Definition at line 194 of file GeometryTranslator.cc.

Referenced by calculateBendAngle().

                                                                        {
  return 0.0;
}

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

private

Definition at line 182 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 187 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 74 of file GeometryTranslator.cc.

References calcCSCSpecificBend(), calcDTSpecificBend(), calcGEMSpecificBend(), calcRPCSpecificBend(), L1TMuon::TriggerPrimitive::kCSC, L1TMuon::TriggerPrimitive::kDT, L1TMuon::TriggerPrimitive::kGEM, 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;
  case TriggerPrimitive::kGEM:
    return calcGEMSpecificBend(tp);
    break;
  default:
    return std::nan("Invalid TP type!");
    break;
  }
}

double GeometryTranslator::calculateGlobalEta

(

const TriggerPrimitive &

tp

)

const

Definition at line 32 of file GeometryTranslator.cc.

References calcCSCSpecificEta(), calcDTSpecificEta(), calcGEMSpecificEta(), calcRPCSpecificEta(), L1TMuon::TriggerPrimitive::kCSC, L1TMuon::TriggerPrimitive::kDT, L1TMuon::TriggerPrimitive::kGEM, 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;
  case TriggerPrimitive::kGEM:
    return calcGEMSpecificEta(tp);
    break;
  default:
    return std::nan("Invalid TP type!"); 
    break;
  }
}

double GeometryTranslator::calculateGlobalPhi

(

const TriggerPrimitive &

tp

)

const

Definition at line 53 of file GeometryTranslator.cc.

References calcCSCSpecificPhi(), calcDTSpecificPhi(), calcGEMSpecificPhi(), calcRPCSpecificPhi(), L1TMuon::TriggerPrimitive::kCSC, L1TMuon::TriggerPrimitive::kDT, L1TMuon::TriggerPrimitive::kGEM, 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;
  case TriggerPrimitive::kGEM:
    return calcGEMSpecificPhi(tp);
    break;
  default:
    return std::nan("Invalid TP type!");
    break;
  }
}

void GeometryTranslator::checkAndUpdateGeometry

(

const edm::EventSetup &

es

)

Definition at line 116 of file GeometryTranslator.cc.

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

Referenced by TrackFinder::process().

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

  const IdealMagneticFieldRecord& magfield = es.get<IdealMagneticFieldRecord>();
  unsigned long long magfieldid = magfield.cacheIdentifier();
  if( _magfield_cache_id != magfieldid ) {
    magfield.get(_magfield);
    _magfield_cache_id = magfieldid;
  }
}

const CSCGeometry& L1TMuon::GeometryTranslator::getCSCGeometry ( ) const

inline

Definition at line 53 of file GeometryTranslator.h.

{ return *_geocsc; }

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

private

Definition at line 204 of file GeometryTranslator.cc.

References _geocsc, CSCGeometry::chamber(), L1TMuon::TriggerPrimitive::detId(), CSCPatternLUT::get2007Position(), L1TMuon::TriggerPrimitive::getCSCData(), triggerObjects_cff::id, isCSCCounterClockwise(), CSCConstants::KEY_ALCT_LAYER, L1TMuon::TriggerPrimitive::CSCData::keywire, PV3DBase< T, PVType, FrameType >::mag(), PFRecoTauDiscriminationByIsolation_cfi::offset, listBenchmarks::pattern, L1TMuon::TriggerPrimitive::CSCData::pattern, PV3DBase< T, PVType, FrameType >::phi(), digi_MixPreMix_cfi::strip, L1TMuon::TriggerPrimitive::CSCData::strip, and PV3DBase< T, PVType, FrameType >::theta().

Referenced by calcCSCSpecificEta(), calcCSCSpecificPhi(), and getGlobalPoint().

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

const DTGeometry& L1TMuon::GeometryTranslator::getDTGeometry ( ) const

inline

Definition at line 54 of file GeometryTranslator.h.

{ return *_geodt;  }

const GEMGeometry& L1TMuon::GeometryTranslator::getGEMGeometry ( ) const

inline

Definition at line 51 of file GeometryTranslator.h.

{ return *_geogem; }

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

private

Definition at line 136 of file GeometryTranslator.cc.

References _geogem, GEMEtaPartition::centreOfPad(), L1TMuon::TriggerPrimitive::detId(), GEMGeometry::etaPartition(), L1TMuon::TriggerPrimitive::getGEMData(), runTauDisplay::gp, triggerObjects_cff::id, L1TMuon::TriggerPrimitive::GEMData::pad, and GeomDet::toGlobal().

Referenced by calcGEMSpecificEta(), calcGEMSpecificPhi(), and getGlobalPoint().

                                                                        {
  const GEMDetId id(tp.detId<GEMDetId>());
  const GEMEtaPartition * roll = _geogem->etaPartition(id);
  const uint16_t pad = tp.getGEMData().pad;
  const LocalPoint lp = roll->centreOfPad(pad);
  const GlobalPoint gp = roll->toGlobal(lp);

  //roll.release();

  return gp;
}

GlobalPoint GeometryTranslator::getGlobalPoint

(

const TriggerPrimitive &

tp

)

const

Definition at line 95 of file GeometryTranslator.cc.

References calcDTSpecificPoint(), getCSCSpecificPoint(), getGEMSpecificPoint(), getRPCSpecificPoint(), L1TMuon::TriggerPrimitive::kCSC, L1TMuon::TriggerPrimitive::kDT, L1TMuon::TriggerPrimitive::kGEM, L1TMuon::TriggerPrimitive::kRPC, and L1TMuon::TriggerPrimitive::subsystem().

Referenced by PrimitiveConversion::convert_csc(), PrimitiveConversion::convert_gem(), and PrimitiveConversion::convert_rpc().

                                                                   {
  switch(tp.subsystem()) {
  case TriggerPrimitive::kDT:
    return calcDTSpecificPoint(tp);
    break;
  case TriggerPrimitive::kCSC:
    return getCSCSpecificPoint(tp);
    break;
  case TriggerPrimitive::kRPC:
    return getRPCSpecificPoint(tp);
    break;
  case TriggerPrimitive::kGEM:
    return getGEMSpecificPoint(tp);
    break;
  default:
    GlobalPoint ret(GlobalPoint::Polar(std::nan("Invalid TP type!"), std::nan("Invalid TP type!"), std::nan("Invalid TP type!")));
    return ret;
    break;
  }
}

const MagneticField& L1TMuon::GeometryTranslator::getMagneticField ( ) const

inline

Definition at line 56 of file GeometryTranslator.h.

{ return *_magfield; }

const RPCGeometry& L1TMuon::GeometryTranslator::getRPCGeometry ( ) const

inline

Definition at line 52 of file GeometryTranslator.h.

{ return *_georpc; }

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

private

Definition at line 167 of file GeometryTranslator.cc.

References _georpc, RPCRoll::centreOfStrip(), L1TMuon::TriggerPrimitive::detId(), L1TMuon::TriggerPrimitive::getRPCData(), runTauDisplay::gp, triggerObjects_cff::id, RPCGeometry::roll(), digi_MixPreMix_cfi::strip, L1TMuon::TriggerPrimitive::RPCData::strip_hi, L1TMuon::TriggerPrimitive::RPCData::strip_low, and GeomDet::toGlobal().

Referenced by calcRPCSpecificEta(), calcRPCSpecificPhi(), and getGlobalPoint().

                                                                        {
  const RPCDetId id(tp.detId<RPCDetId>());
  const RPCRoll * roll = _georpc->roll(id);
  //const int strip = static_cast<int>(tp.getRPCData().strip);
  // Use half-strip precision, - 0.5 at the end to get the center of the strip
  const float strip = (0.5 * static_cast<float>(tp.getRPCData().strip_low + tp.getRPCData().strip_hi)) - 0.5;
  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 338 of file GeometryTranslator.cc.

References funct::abs(), and M_PI.

Referenced by calcDTSpecificBend(), and 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)     );  
}

Member Data Documentation

edm::ESHandle<CSCGeometry> L1TMuon::GeometryTranslator::_geocsc

private

Definition at line 63 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry(), and getCSCSpecificPoint().

edm::ESHandle<DTGeometry> L1TMuon::GeometryTranslator::_geodt

private

Definition at line 64 of file GeometryTranslator.h.

Referenced by calcDTSpecificPoint(), and checkAndUpdateGeometry().

edm::ESHandle<GEMGeometry> L1TMuon::GeometryTranslator::_geogem

private

Definition at line 61 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry(), and getGEMSpecificPoint().

unsigned long long L1TMuon::GeometryTranslator::_geom_cache_id

private

Definition at line 60 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().

edm::ESHandle<RPCGeometry> L1TMuon::GeometryTranslator::_georpc

private

Definition at line 62 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry(), and getRPCSpecificPoint().

edm::ESHandle<MagneticField> L1TMuon::GeometryTranslator::_magfield

private

Definition at line 67 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().

unsigned long long L1TMuon::GeometryTranslator::_magfield_cache_id

private

Definition at line 66 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().