MuonFixedId.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/
 
/***************************************************************************
 * Muon Calibration Utilities
 * -----------------------------------------
 *
 * Authors       : Martin Woudstra, Zdenko van Kesteren, Peter Kluit
 * Creation Date: 21 March 2005
 * Updated : 03 December 2007
 * Updated : 16 January  2015 E. Diehl Add BME/BOE/BMG chambers 
 *   Note: BOE=BOL in offline, so do not need to explicitly add BOE, it will be handled as BOL
 ***************************************************************************/
 
 
#ifndef MUONCALIBIDENTIFIER_MUONFIXEDID_H
#define MUONCALIBIDENTIFIER_MUONFIXEDID_H
// std
#include <iostream>
#include <iomanip>
#include <string>
 
#include <cstdlib>
 
// other packages
 
// this package
//#include "MuonCalibIdentifier/MuonStationName.h"
 
 
namespace MuonCalib {
 
class MuonFixedId {
public:
   MuonFixedId(); 
   explicit MuonFixedId( unsigned int id );
 
   unsigned int getIdInt() const;
   void clear();
   bool isValid() const;
 
   bool is_mdt() const;
   bool is_tgc() const;
   bool is_csc() const;
   bool is_rpc() const;
 
   void set_mdt();
   void set_tgc();
   void set_csc();
   void set_rpc();
 
 
   bool operator==( const MuonFixedId& rhs ) const;
   bool operator!=( const MuonFixedId& rhs ) const;
   bool operator<( const MuonFixedId& rhs ) const;
   
   std::ostream& dump(std::ostream& os) const;
 
//   friend std::istream& operator>>( std::istream& is, MuonCalib::MuonFixedId& id );
 
// Muon generic code
 
   bool setTechnologyIndex( unsigned int idx );
   bool setStationNameIndex( unsigned int idx );
   bool setStationEtaIndex( unsigned int idx );
   bool setStationPhiIndex( unsigned int idx );
 
   bool setTechnology( int num );
   bool setStationName( int num );
   bool setStationEta( int num );
   bool setStationPhi( int num );
 
   int technologyIndex() const;
   int stationNameIndex() const;
   int etaIndex() const;
   int phiIndex() const;
 
   int technology() const;
   std::string technologyString() const;
   int stationName() const;
   std::string stationNameString() const;
   int eta() const;
   int phi() const;
   
   static int technologyMin();
   static int stationNameMin();
   static int etaMin();
   static int phiMin();
 
   static int technologyMax();
   static int stationNameMax();
   static int etaMax();
   static int phiMax();
 
   static int stationStringToFixedStationNumber(const std::string& station);
   static std::string stationNumberToFixedStationString(const int station);
 
 
//  Mdt specific code
 
   bool setMdtTubeIndex( unsigned int idx ) ;
   bool setMdtTubeLayerIndex( unsigned int idx ) ;
   bool setMdtMultilayerIndex( unsigned int idx ) ;
 
   bool setMdtTube( int num );
   bool setMdtTubeLayer( int num );
   bool setMdtMultilayer( int num );
 
   int mdtTubeIndex() const;
   int mdtTubeLayerIndex() const;
   int mdtMultilayerIndex() const;
 
   int mdtTube() const;
   int mdtTubeLayer() const;
   int mdtMultilayer() const;
   int mdtMezzanine() const;
 
   MuonFixedId mdtChamberId() const;
 
   MuonFixedId mdtMultilayerId() const;
 
   static int mdtTubeMin();
   static int mdtTubeLayerMin();
   static int mdtMultilayerMin();
 
   static int mdtTubeMax();
   static int mdtTubeLayerMax();
   static int mdtMultilayerMax();
 
   
//  Csc specific code
 
   bool setCscChamberLayerIndex( unsigned int idx );
   bool setCscWireLayerIndex( unsigned int idx );
   bool setCscMeasuresPhiIndex( unsigned int idx );
   bool setCscStripIndex( unsigned int idx );
 
   bool setCscChamberLayer( int num );
   bool setCscWireLayer( int num );
   bool setCscMeasuresPhi( int num );
   bool setCscStrip( int num );  
 
   int cscChamberLayerIndex() const;
   int cscWireLayerIndex() const;
   int cscMeasuresPhiIndex() const;
   int cscStripIndex() const;
 
   int cscChamberLayer() const;
   int cscWireLayer() const;
   int cscMeasuresPhi() const;
   int cscStrip() const;
 
   static int cscChamberLayerMin(); 
   static int cscWireLayerMin(); 
   static int cscMeasuresPhiMin(); 
   static int cscStripMin();
   
   static int cscChamberLayerMax(); 
   static int cscWireLayerMax(); 
   static int cscMeasuresPhiMax(); 
   static int cscStripMax();
 
//  Rpc specific code 
 
   bool setRpcDoubletRIndex( unsigned int idx );
   bool setRpcDoubletZIndex( unsigned int idx );
   bool setRpcDoubletPhiIndex( unsigned int idx );
   bool setRpcGasGapIndex( unsigned int idx );
   bool setRpcMeasuresPhiIndex( unsigned int idx );
   bool setRpcStripIndex( unsigned int idx );
 
   bool setRpcDoubletR( int num );
   bool setRpcDoubletZ( int num );
   bool setRpcDoubletPhi( int num );
   bool setRpcGasGap( int num );
   bool setRpcMeasuresPhi( int num );
   bool setRpcStrip( int num );
   
   int rpcDoubletRIndex() const;
   int rpcDoubletZIndex() const;
   int rpcDoubletPhiIndex() const;
   int rpcGasGapIndex() const;
   int rpcMeasuresPhiIndex() const;
   int rpcStripIndex() const;
   
   int rpcDoubletR() const;
   int rpcDoubletZ() const;
   int rpcDoubletPhi() const;
   int rpcGasGap() const;
   int rpcMeasuresPhi() const;
   int rpcStrip() const;
 
//  Tgc specific code
 
   bool setTgcGasGapIndex( unsigned int idx );  
   bool setTgcIsStripIndex( unsigned int idx );
   bool setTgcChannelIndex( unsigned int idx );  
 
   bool setTgcGasGap( int num );
   bool setTgcIsStrip( int num );
   bool setTgcChannel( int num );
 
   int tgcGasGapIndex() const;   
   int tgcIsStripIndex() const;   
   int tgcChannelIndex() const;   
 
   int tgcGasGap() const;
   int tgcIsStrip() const;
   int tgcChannel() const;
 
   static const int technologyMDT = 0;
   static const int technologyCSC = 1;
   static const int technologyTGC = 2;
   static const int technologyRPC = 3;
   
private:
   // the member holding the packed fields
   unsigned int m_data;
 
   static const unsigned int kInvalidData = 0xFFFFFFFF;
 
//   Muon generic
 
   static const unsigned int kTechnologyMask   = 3;
   static const unsigned int kTechnologyShift  = 30;
   static const unsigned int kTechnologyMin    = 0;
   
   static const int kNumberOfTechnologies = 4;
   static const char kTechnologyStrings[kNumberOfTechnologies][4];
 
   static const unsigned int kStationNameMask      = 63;
   static const unsigned int kStationNameShift     = 24;
   static const unsigned int kStationNameMin       = 1;
 
   static const unsigned int kStationEtaMask          = 31;
   static const unsigned int kStationEtaShift         = 19;
   static const int kStationEtaMin                    = -8;
 
   static const unsigned int kStationPhiMask          = 63;
   static const unsigned int kStationPhiShift         = 13;
   static const unsigned int kStationPhiMin           = 1;
 
   static const int kNumberOfStationNames = 35;
   static const char kStationNameStrings[kNumberOfStationNames][4];
 
   // the full station information for making a station identifier
   static const unsigned int kStationShift = kStationPhiShift; // the smallest shift
   static const unsigned int kStationMask = 
     ( (kTechnologyMask  << kTechnologyShift)  |
       (kStationNameMask << kStationNameShift) |
       (kStationEtaMask  << kStationEtaShift)  |
       (kStationPhiMask  << kStationPhiShift)    ) >> kStationShift; // kStationShift for consistent meaning
 
//   Mdt specific code 
 
   static const unsigned int kMdtMultilayerMask    = 1;
   static const unsigned int kMdtMultilayerShift   = 9;
   static const unsigned int kMdtMultilayerMin     = 1;
 
   static const unsigned int kMdtTubeLayerMask     = 3;
   static const unsigned int kMdtTubeLayerShift    = 7;
   static const unsigned int kMdtTubeLayerMin      = 1;
 
   static const unsigned int kMdtTubeMask          = 127;
   static const unsigned int kMdtTubeShift         = 0;
   static const unsigned int kMdtTubeMin           = 1;
 
//  Csc specific code
 
   static const unsigned int kCscChamberLayerMask  = 1;
   static const unsigned int kCscChamberLayerShift = 11;
   static const unsigned int kCscChamberLayerMin   = 1;
 
   static const unsigned int kCscWireLayerMask     = 3;
   static const unsigned int kCscWireLayerShift    = 9;
   static const unsigned int kCscWireLayerMin      = 1;
 
   static const unsigned int kCscMeasuresPhiMask   = 1;
   static const unsigned int kCscMeasuresPhiShift  = 8;
   static const unsigned int kCscMeasuresPhiMin    = 0;
 
   static const unsigned int kCscStripMask         = 255;
   static const unsigned int kCscStripShift        = 0;
   static const unsigned int kCscStripMin          = 1;
 
//  Rpc specific code
 
   static const unsigned int kRpcDoubletRMask      = 1;
   static const unsigned int kRpcDoubletRShift     = 12; 
   static const unsigned int kRpcDoubletRMin       = 1;
 
   static const unsigned int kRpcDoubletZMask      = 3;
   static const unsigned int kRpcDoubletZShift     = 10; 
   static const unsigned int kRpcDoubletZMin       = 1;
 
   static const unsigned int kRpcDoubletPhiMask    = 1;
   static const unsigned int kRpcDoubletPhiShift   = 9; 
   static const unsigned int kRpcDoubletPhiMin     = 1;
 
   static const unsigned int kRpcGasGapMask        = 1;
   static const unsigned int kRpcGasGapShift       = 8;  
   static const unsigned int kRpcGasGapMin         = 1;
 
   static const unsigned int kRpcMeasuresPhiMask   = 1;
   static const unsigned int kRpcMeasuresPhiShift  = 7;  
   static const unsigned int kRpcMeasuresPhiMin    = 0;
 
   static const unsigned int kRpcStripMask         = 127;
   static const unsigned int kRpcStripShift        = 0; 
   static const unsigned int kRpcStripMin          = 1;
 
// Tgc specific code 
   static const unsigned int kTgcGasGapMask        = 3;
   static const unsigned int kTgcGasGapShift       = 9;
   static const unsigned int kTgcGasGapMin         = 1;
 
   static const unsigned int kTgcIsStripMask       = 1;
   static const unsigned int kTgcIsStripShift      = 8;
   static const unsigned int kTgcIsStripMin        = 0;
 
   static const unsigned int kTgcChannelMask       = 255;
   static const unsigned int kTgcChannelShift      = 0;
   static const unsigned int kTgcChannelMin        = 1;
 
};
 
// Muon generic methods
 
inline MuonFixedId::MuonFixedId()
  : m_data(kInvalidData){
  }
 
inline MuonFixedId::MuonFixedId( unsigned int id ) : m_data(id){
}
 
inline unsigned int MuonFixedId::getIdInt() const{
  return m_data;
}
/* void MuonFixedId::setMuonIdentifier( int id ){ */
/*   if( id>-1 || id<kInvalidData ) m_data = id; */
/*   else  m_data = kInvalidData; */
/* } */
 
inline void MuonFixedId::clear() {
   m_data = kInvalidData;
}
 
inline bool MuonFixedId::isValid() const {
  return m_data != kInvalidData;
}
 
inline bool MuonFixedId::is_mdt() const{
  return isValid() && technology() == technologyMDT;
}
 
inline bool MuonFixedId::is_tgc() const{
  return isValid() && technology() == technologyTGC;
}
 
inline bool MuonFixedId::is_csc() const{
  return isValid() && technology() == technologyCSC;
}
 
inline bool MuonFixedId::is_rpc() const{
  return isValid() && technology() == technologyRPC;
}
 
inline void MuonFixedId::set_mdt() {
  setTechnology( technologyMDT );
}
 
inline void MuonFixedId::set_csc() {
  setTechnology( technologyCSC );
}
 
inline void MuonFixedId::set_tgc() {
  setTechnology( technologyTGC );
}
 
inline void MuonFixedId::set_rpc() {
  setTechnology( technologyRPC );
}
 
inline bool MuonFixedId::operator==( const MuonFixedId& rhs ) const {
   return m_data == rhs.m_data;
}
 
inline bool MuonFixedId::operator!=( const MuonFixedId& rhs ) const {
   return m_data != rhs.m_data;
}
 
inline bool MuonFixedId::operator<( const MuonFixedId& rhs ) const {
   return m_data < rhs.m_data;
}
 
inline bool MuonFixedId::setTechnologyIndex( unsigned int idx ) { 
  if ( idx > kTechnologyMask ) {
      clear();
      return false;
   }
   m_data &= ~(kTechnologyMask << kTechnologyShift);
   m_data |= (idx & kTechnologyMask) << kTechnologyShift;
   return true;
}
 
inline bool MuonFixedId::setTechnology( int num ) { 
   return setTechnologyIndex( num - kTechnologyMin );
}
 
inline int MuonFixedId::technologyIndex() const { 
   return (m_data >> kTechnologyShift) & kTechnologyMask;
}
 
inline int MuonFixedId::technology() const {
   return technologyIndex() + kTechnologyMin;
}
 
inline std::string MuonFixedId::technologyString() const {
  int index = technologyIndex();
  if ( index >= 0 && index < kNumberOfTechnologies ) {
    return kTechnologyStrings[index];
  }
  return "???";
}
 
inline bool MuonFixedId::setStationNameIndex( unsigned int idx ) {
  if ( idx > kStationNameMask ) {
      clear();
      return false;
   }
   m_data &= ~(kStationNameMask << kStationNameShift);
   m_data |= (idx & kStationNameMask) << kStationNameShift;
   return true;
}
 
inline bool MuonFixedId::setStationName( int num ) {
   return setStationNameIndex( num - kStationNameMin );
}
 
inline int MuonFixedId::stationNameIndex() const {
   return (m_data >> kStationNameShift) & kStationNameMask;
}
 
inline int MuonFixedId::stationName() const {
   return stationNameIndex() + kStationNameMin;
}
 
inline std::string MuonFixedId::stationNameString() const {
  int index = stationNameIndex();
  if ( index >= 0 && index < kNumberOfStationNames ) {
    return kStationNameStrings[index];
  } 
  return "XXX";    
}
 
inline bool MuonFixedId::setStationEtaIndex( unsigned int idx ) {
  if ( idx > kStationEtaMask ) {
      clear();
      return false;
   }
   m_data &= ~(kStationEtaMask << kStationEtaShift);
   m_data |= (idx & kStationEtaMask) << kStationEtaShift;
   return true;
}
 
inline bool MuonFixedId::setStationEta( int num ) {
   return setStationEtaIndex( num - kStationEtaMin );
}
 
inline int MuonFixedId::etaIndex() const {
   return (m_data >> kStationEtaShift) & kStationEtaMask;
}
 
inline int MuonFixedId::eta() const {
   return etaIndex() + kStationEtaMin;
}
 
 
inline bool MuonFixedId::setStationPhiIndex( unsigned int idx ) {
  if ( idx > kStationPhiMask ) {
      clear();
      return false;
   }
   m_data &= ~(kStationPhiMask << kStationPhiShift);
   m_data |= (idx & kStationPhiMask) << kStationPhiShift;
   return true;
}
 
inline bool MuonFixedId::setStationPhi( int num ) {
   return setStationPhiIndex( num - kStationPhiMin );
}
 
inline int MuonFixedId::phiIndex() const {
   return (m_data >> kStationPhiShift) & kStationPhiMask;
}
 
inline int MuonFixedId::phi() const {
   return phiIndex() + kStationPhiMin;
}
  
inline int MuonFixedId::technologyMin(){
  return kTechnologyMin;
}
 
inline int MuonFixedId::stationNameMin(){
  return kStationNameMin;
}
 
inline int MuonFixedId::etaMin(){ 
  return kStationEtaMin;
}
 
inline int MuonFixedId::phiMin(){
  return kStationPhiMin;
}
 
inline int MuonFixedId::technologyMax(){
  return kTechnologyMin + kTechnologyMask;
}
 
inline int MuonFixedId::stationNameMax(){ 
  return kStationNameMin + kStationNameMask;
}
 
inline int MuonFixedId::etaMax(){ 
  return kStationEtaMin + kStationEtaMask;
}
 
inline int MuonFixedId::phiMax(){
   return kStationPhiMin + kStationPhiMask; 
}
 
inline int MuonFixedId::stationStringToFixedStationNumber(const std::string& station ) {
  for ( int i = 0; i < kNumberOfStationNames; ++i ) {
    if ( station == kStationNameStrings[i] ) return i + kStationNameMin;
  }
  return -1;  // signal error if not found
}
 
inline std::string MuonFixedId::stationNumberToFixedStationString(const int station) {
  int index = station - kStationNameMin;
  if ( index >= 0 && index < kNumberOfStationNames ) {
    return kStationNameStrings[index];
  }
  return "XXX";   // signal error if not found
}
 
// Mdt specific methods
 
inline bool MuonFixedId::setMdtTubeIndex( unsigned int idx ) {
  if ( idx > kMdtTubeMask ) {
      clear();
      return false;
   }
   m_data &= ~(kMdtTubeMask << kMdtTubeShift);
   m_data |= (idx & kMdtTubeMask) << kMdtTubeShift;
   return true;
}
 
inline bool MuonFixedId::setMdtTube( int num ) {
   return setMdtTubeIndex( num - kMdtTubeMin );
}
 
inline int MuonFixedId::mdtTubeIndex() const {
   return (m_data >> kMdtTubeShift) & kMdtTubeMask;
}
 
inline int MuonFixedId::mdtTube() const {
   return mdtTubeIndex() + kMdtTubeMin;
}
 
inline int MuonFixedId::mdtMezzanine() const {
   int Ichamber =  stationName() ;
   int Ieta = eta() ;
   int Iphi = phi() ;
   int Iml  = mdtMultilayer() ;
   int Itube = mdtTube() ;
   if (Ieta < 0 ) Ieta = -Ieta*10 ;
   int ImezzTubes = 8 ;
   if (Ichamber==1 || Ichamber==2 || Ichamber==8 || Ichamber==13 || Ichamber==7 || Ichamber==15 || Ichamber==24) ImezzTubes=6;
//exception BIS eta=8 is a 3-layer chamber
   if(Ichamber==2 && std::abs(Ieta)==8) ImezzTubes=8;
   int Imezz = ((Itube-1)/ImezzTubes)*2+(Iml-1) ;
   Imezz = 9*100000000 + Ichamber*1000000 + Ieta*10000 + Iphi*100 + Imezz ;
   return Imezz ;
}
 
inline bool MuonFixedId::setMdtTubeLayerIndex( unsigned int idx ) {
  if ( idx > kMdtTubeLayerMask ) {
      clear();
      return false;
   }
   m_data &= ~(kMdtTubeLayerMask << kMdtTubeLayerShift);
   m_data |= (idx & kMdtTubeLayerMask) << kMdtTubeLayerShift;
   return true;
}
 
inline bool MuonFixedId::setMdtTubeLayer( int num ) {
   return setMdtTubeLayerIndex( num - kMdtTubeLayerMin );
}
 
inline int MuonFixedId::mdtTubeLayerIndex() const {
   return (m_data >> kMdtTubeLayerShift) & kMdtTubeLayerMask;
}
 
inline int MuonFixedId::mdtTubeLayer() const {
   return mdtTubeLayerIndex() + kMdtTubeLayerMin;
}
 
inline bool MuonFixedId::setMdtMultilayerIndex( unsigned int idx ) {
  if ( idx > kMdtMultilayerMask ) {
      clear();
      return false;
   }
   m_data &= ~(kMdtMultilayerMask << kMdtMultilayerShift);
   m_data |= (idx & kMdtMultilayerMask) << kMdtMultilayerShift;
   return true;
}
 
inline bool MuonFixedId::setMdtMultilayer( int num ) {
   return setMdtMultilayerIndex( num - kMdtMultilayerMin );
}
 
inline int MuonFixedId::mdtMultilayerIndex() const {
   return (m_data >> kMdtMultilayerShift) & kMdtMultilayerMask;
}
 
inline int MuonFixedId::mdtMultilayer() const {
   return mdtMultilayerIndex() + kMdtMultilayerMin;
}
 
inline int MuonFixedId::mdtTubeMin(){
  return kMdtTubeMin;
}
 
inline int MuonFixedId::mdtTubeMax(){
  return kMdtTubeMin + kMdtTubeMask;
}
 
inline int MuonFixedId::mdtTubeLayerMin(){
  return kMdtTubeLayerMin;
}
 
inline int MuonFixedId::mdtTubeLayerMax(){
  return kMdtTubeLayerMin + kMdtTubeLayerMask;
}
 
inline int MuonFixedId::mdtMultilayerMin(){
  return kMdtMultilayerMin;
}
 
inline int MuonFixedId::mdtMultilayerMax(){
  return kMdtMultilayerMin + kMdtMultilayerMask;
}
 
inline MuonFixedId MuonFixedId::mdtChamberId() const {
  // mdt chamber id = muon station id
  return MuonFixedId( m_data & (kStationMask << kStationShift) );
}
 
inline MuonFixedId MuonFixedId::mdtMultilayerId() const {
  // mdt multilayer id = muon station id + multilayer field
  return MuonFixedId( m_data & ( (kStationMask << kStationShift) | (kMdtMultilayerMask << kMdtMultilayerShift) ) );
}
 
 
// Csc specific methods
 
inline bool MuonFixedId::setCscChamberLayerIndex( unsigned int idx ) {
  if ( idx > kCscChamberLayerMask ) {
      clear();
      return false;
   }
   m_data &= ~( kCscChamberLayerMask <<  kCscChamberLayerShift );
   m_data |= (idx &  kCscChamberLayerMask) << kCscChamberLayerShift;
   return true;
}
 
inline bool MuonFixedId::setCscChamberLayer( int num ) {
   return setCscChamberLayerIndex( num -  kCscChamberLayerMin );
}
 
inline int MuonFixedId::cscChamberLayerIndex() const {
   return (m_data >> kCscChamberLayerShift) & kCscChamberLayerMask;
}
 
inline int MuonFixedId::cscChamberLayer() const {
   return cscChamberLayerIndex() + kCscChamberLayerMin;
}
 
 
inline bool MuonFixedId::setCscWireLayerIndex( unsigned int idx ) {
  if ( idx > kCscWireLayerMask ) {
      clear();
      return false;
   }
   m_data &= ~( kCscWireLayerMask <<  kCscWireLayerShift );
   m_data |= (idx &  kCscWireLayerMask) << kCscWireLayerShift;
   return true;
}
 
inline bool MuonFixedId::setCscWireLayer( int num ) {
  return setCscWireLayerIndex( num -  kCscWireLayerMin );
}
 
inline int MuonFixedId::cscWireLayerIndex() const {
   return (m_data >> kCscWireLayerShift) & kCscWireLayerMask;
}
 
inline int MuonFixedId::cscWireLayer() const {
   return cscWireLayerIndex() + kCscWireLayerMin;
}
 
inline bool MuonFixedId::setCscMeasuresPhiIndex( unsigned int idx ) {
  if ( idx > kCscMeasuresPhiMask ) {
      clear();
      return false;
   }
   m_data &= ~( kCscMeasuresPhiMask <<  kCscMeasuresPhiShift );
   m_data |= (idx &  kCscMeasuresPhiMask) << kCscMeasuresPhiShift;
   return true;
}
 
inline bool MuonFixedId::setCscMeasuresPhi( int num ) {
  return setCscMeasuresPhiIndex( num -  kCscMeasuresPhiMin );
}
 
inline int MuonFixedId::cscMeasuresPhiIndex() const {
   return (m_data >> kCscMeasuresPhiShift) & kCscMeasuresPhiMask;
}
 
inline int MuonFixedId::cscMeasuresPhi() const {
   return cscMeasuresPhiIndex() + kCscMeasuresPhiMin;
}
 
inline bool MuonFixedId::setCscStripIndex( unsigned int idx ) {
  if ( idx > kCscStripMask ) {
      clear();
      return false;
   }
   m_data &= ~( kCscStripMask <<  kCscStripShift );
   m_data |= (idx &  kCscStripMask) << kCscStripShift;
   return true;
}
 
inline bool MuonFixedId::setCscStrip( int num ) {
  return setCscStripIndex( num -  kCscStripMin );
}
 
inline int MuonFixedId::cscStripIndex() const {
   return (m_data >> kCscStripShift) & kCscStripMask;
}
 
inline int MuonFixedId::cscStrip() const {
   return cscStripIndex() + kCscStripMin;
}
 
inline int MuonFixedId::cscChamberLayerMin(){
  return kCscChamberLayerMin;
}
 
inline int MuonFixedId::cscChamberLayerMax(){
  return kCscChamberLayerMin + kCscChamberLayerMask;
}
 
inline int MuonFixedId::cscWireLayerMin(){
  return kCscWireLayerMin;
}
 
inline int MuonFixedId::cscWireLayerMax(){
  return kCscWireLayerMin + kCscWireLayerMask;
}
 
inline int MuonFixedId::cscMeasuresPhiMin(){
  return kCscMeasuresPhiMin;
}
 
inline int MuonFixedId::cscMeasuresPhiMax(){
  return kCscMeasuresPhiMin + kCscMeasuresPhiMask;
}
 
inline int MuonFixedId::cscStripMin(){
  return kCscStripMin;
}
 
inline int MuonFixedId::cscStripMax(){
  return kCscStripMin + kCscStripMask;
}
 
// Rpc specific methods
inline bool MuonFixedId::setRpcDoubletRIndex( unsigned int idx ) {
  if ( idx > kRpcDoubletRMask ) {
      clear();
      return false;
   }
   m_data &= ~( kRpcDoubletRMask <<  kRpcDoubletRShift );
   m_data |= (idx &  kRpcDoubletRMask) << kRpcDoubletRShift;
   return true;
}
 
inline bool MuonFixedId::setRpcDoubletR( int num ) {
  return setRpcDoubletRIndex( num -  kRpcDoubletRMin );
}
 
inline int MuonFixedId::rpcDoubletRIndex() const {
   return (m_data >> kRpcDoubletRShift) & kRpcDoubletRMask;
}
 
inline int MuonFixedId::rpcDoubletR() const {
   return rpcDoubletRIndex() + kRpcDoubletRMin;
}
 
inline bool MuonFixedId::setRpcDoubletZIndex( unsigned int idx ) {
  if ( idx > kRpcDoubletZMask ) {
      clear();
      return false;
   }
   m_data &= ~( kRpcDoubletZMask <<  kRpcDoubletZShift );
   m_data |= (idx &  kRpcDoubletZMask) << kRpcDoubletZShift;
   return true;
}
 
inline bool MuonFixedId::setRpcDoubletZ( int num ) {
  return setRpcDoubletZIndex( num -  kRpcDoubletZMin );
}
 
inline int MuonFixedId::rpcDoubletZIndex() const {
   return (m_data >> kRpcDoubletZShift) & kRpcDoubletZMask;
}
 
inline int MuonFixedId::rpcDoubletZ() const {
   return rpcDoubletZIndex() + kRpcDoubletZMin;
}
 
inline bool MuonFixedId::setRpcDoubletPhiIndex( unsigned int idx ) {
  if ( idx > kRpcDoubletPhiMask ) {
      clear();
      return false;
   }
   m_data &= ~( kRpcDoubletPhiMask <<  kRpcDoubletPhiShift );
   m_data |= (idx &  kRpcDoubletPhiMask) << kRpcDoubletPhiShift;
   return true;
}
 
inline bool MuonFixedId::setRpcDoubletPhi( int num ) {
  return setRpcDoubletPhiIndex( num -  kRpcDoubletPhiMin );
}
 
inline int MuonFixedId::rpcDoubletPhiIndex() const {
   return (m_data >> kRpcDoubletPhiShift) & kRpcDoubletPhiMask;
}
 
inline int MuonFixedId::rpcDoubletPhi() const {
   return rpcDoubletPhiIndex() + kRpcDoubletPhiMin;
}
 
inline bool MuonFixedId::setRpcGasGapIndex( unsigned int idx ) {
  if ( idx > kRpcGasGapMask ) { 
      clear();
      return false;
   }
   m_data &= ~( kRpcGasGapMask <<  kRpcGasGapShift );
   m_data |= (idx &  kRpcGasGapMask) << kRpcGasGapShift;
   return true;
}
 
inline bool MuonFixedId::setRpcGasGap( int num ) {
  return setRpcGasGapIndex( num -  kRpcGasGapMin );
}
 
inline int MuonFixedId::rpcGasGapIndex() const {
   return (m_data >> kRpcGasGapShift) & kRpcGasGapMask;
}
 
inline int MuonFixedId::rpcGasGap() const {
   return rpcGasGapIndex() + kRpcGasGapMin;
}
 
inline bool MuonFixedId::setRpcMeasuresPhiIndex( unsigned int idx ) {
  if ( idx > kRpcMeasuresPhiMask ) {
      clear();
      return false;
   }
   m_data &= ~( kRpcMeasuresPhiMask <<  kRpcMeasuresPhiShift );
   m_data |= (idx &  kRpcMeasuresPhiMask) << kRpcMeasuresPhiShift;
   return true;
}
 
inline bool MuonFixedId::setRpcMeasuresPhi( int num ) {
  return setRpcMeasuresPhiIndex( num -  kRpcMeasuresPhiMin );
}
 
inline int MuonFixedId::rpcMeasuresPhiIndex() const {
   return (m_data >> kRpcMeasuresPhiShift) & kRpcMeasuresPhiMask;
}
 
inline int MuonFixedId::rpcMeasuresPhi() const {
   return rpcMeasuresPhiIndex() + kRpcMeasuresPhiMin;
}
 
inline bool MuonFixedId::setRpcStripIndex( unsigned int idx ) {
  if ( idx > kRpcStripMask ) {
      clear();
      return false;
   }
   m_data &= ~( kRpcStripMask <<  kRpcStripShift );
   m_data |= (idx &  kRpcStripMask) << kRpcStripShift;
   return true;
}
 
inline bool MuonFixedId::setRpcStrip( int num ) {
  return setRpcStripIndex( num -  kRpcStripMin );
}
 
inline int MuonFixedId::rpcStripIndex() const {
   return (m_data >> kRpcStripShift) & kRpcStripMask;
}
 
inline int MuonFixedId::rpcStrip() const {
   return rpcStripIndex() + kRpcStripMin;
}
 
 
//  Tgc specific methods
inline bool MuonFixedId::setTgcGasGapIndex( unsigned int idx ) {
  if ( idx > kTgcGasGapMask ) {
      clear();
      return false;
   }
   m_data &= ~( kTgcGasGapMask <<  kTgcGasGapShift );
   m_data |= (idx &  kTgcGasGapMask) << kTgcGasGapShift;
   return true;
}
 
inline bool MuonFixedId::setTgcGasGap( int num ) {
  return setTgcGasGapIndex( num -  kTgcGasGapMin );
}
 
inline int MuonFixedId::tgcGasGapIndex() const {
   return (m_data >> kTgcGasGapShift) & kTgcGasGapMask;
}
 
inline int MuonFixedId::tgcGasGap() const {
   return tgcGasGapIndex() + kTgcGasGapMin;
}
 
inline bool MuonFixedId::setTgcIsStripIndex( unsigned int idx ) {
  if ( idx > kTgcIsStripMask ) {
      clear();
      return false;
   }
   m_data &= ~( kTgcIsStripMask <<  kTgcIsStripShift );
   m_data |= (idx &  kTgcIsStripMask) << kTgcIsStripShift;
   return true;
}
 
inline bool MuonFixedId::setTgcIsStrip( int num ) {
  return setTgcIsStripIndex( num -  kTgcIsStripMin );
}
 
inline int MuonFixedId::tgcIsStripIndex() const {
   return (m_data >> kTgcIsStripShift) & kTgcIsStripMask;
}
 
inline int MuonFixedId::tgcIsStrip() const {
   return tgcIsStripIndex() + kTgcIsStripMin;
}
 
 
inline bool MuonFixedId::setTgcChannelIndex( unsigned int idx ) {
  if ( idx > kTgcChannelMask ) {
      clear();
      return false;
   }
   m_data &= ~( kTgcChannelMask <<  kTgcChannelShift );
   m_data |= (idx &  kTgcChannelMask) << kTgcChannelShift;
   return true;
}
 
inline bool MuonFixedId::setTgcChannel( int num ) {
  return setTgcChannelIndex( num -  kTgcChannelMin );
}
 
inline int MuonFixedId::tgcChannelIndex() const {
   return (m_data >> kTgcChannelShift) & kTgcChannelMask;
}
 
inline int MuonFixedId::tgcChannel() const {
   return tgcChannelIndex() + kTgcChannelMin;
}
 
 
inline std::ostream& operator<<( std::ostream& os, const MuonCalib::MuonFixedId& id ) {
  return id.dump( os );
}
} //MuonCalib namespace 
 
/*
***** must solve conflicts between MuonCalib namespace and friendship
inline std::istream& operator>>( std::istream& is, MuonCalib::MuonFixedId& id){
  is >> id.m_data ;
  return is;
}
*/
#endif // MUONCALIBIDENTIFIER_MUONFIXEDID_H