DTHVStatus.cc

Go to the documentation of this file.

/*
 *  See header file for a description of this class.
 *
 *  \author Paolo Ronchese INFN Padova
 *
 */

//----------------------
// This Class' Header --
//----------------------
#include "CondFormats/DTObjects/interface/DTHVStatus.h"

//-------------------------------
// Collaborating Class Headers --
//-------------------------------
#include "CondFormats/DTObjects/interface/DTBufferTree.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "DataFormats/MuonDetId/interface/DTLayerId.h"
#include "DataFormats/MuonDetId/interface/DTChamberId.h"

//---------------
// C++ Headers --
//---------------
#include <iostream>
#include <sstream>

//-------------------
// Initializations --
//-------------------


//----------------
// Constructors --
//----------------
DTHVStatus::DTHVStatus():
  dataVersion( " " ),
  dBuf(new DTBufferTree<int,int>) {
  dataList.reserve( 10 );
}


DTHVStatus::DTHVStatus( const std::string& version ):
  dataVersion( version ),
  dBuf(new DTBufferTree<int,int>) {
  dataList.reserve( 10 );
}


DTHVStatusId::DTHVStatusId() :
    wheelId( 0 ),
  stationId( 0 ),
   sectorId( 0 ),
       slId( 0 ),
    layerId( 0 ),
     partId( 0 ) {
}


DTHVStatusData::DTHVStatusData() :
  flagA( 0 ),
  flagC( 0 ),
  flagS( 0 ) {
}


//--------------
// Destructor --
//--------------
DTHVStatus::~DTHVStatus() {
}


DTHVStatusId::~DTHVStatusId() {
}


DTHVStatusData::~DTHVStatusData() {
}


//--------------
// Operations --
//--------------
int DTHVStatus::get( int   wheelId,
                     int stationId,
                     int  sectorId,
                     int      slId,
                     int   layerId,
                     int    partId,
                     int&    fCell,
                     int&    lCell,
                     int&    flagA,
                     int&    flagC,
                     int&    flagS ) const {
  fCell =
  lCell =
  flagA =
  flagC = 
  flagS = 0;

  std::vector<int> chanKey;
  chanKey.reserve(6);
  chanKey.push_back(   wheelId );
  chanKey.push_back( stationId );
  chanKey.push_back(  sectorId );
  chanKey.push_back(      slId );
  chanKey.push_back(   layerId );
  chanKey.push_back(    partId );
  int ientry;
  int searchStatus = dBuf->find( chanKey.begin(), chanKey.end(), ientry );
  if ( !searchStatus ) {
    const DTHVStatusData& data( dataList[ientry].second );
    fCell = data.fCell;
    lCell = data.lCell;
    flagA = data.flagA;
    flagC = data.flagC;
    flagS = data.flagS;
  }

  return searchStatus;

}


int DTHVStatus::get( const DTLayerId& id,
                     int    partId,
                     int&    fCell,
                     int&    lCell,
                     int&    flagA,
                     int&    flagC,
                     int&    flagS ) const {
  return get( id.wheel(),
              id.station(),
              id.sector(),
              id.superLayer(),
              id.layer(),
              partId,
              fCell, lCell, flagA, flagC, flagS );
}


int DTHVStatus::get( const DTWireId& id,
                     int&         flagA,
                     int&         flagC,
                     int&         flagS ) const {
  flagA = flagC = flagS = 0;
  int iCell = id.wire();
  int fCell;
  int lCell;
  int
  fCheck = get( id.wheel(),
                id.station(),
                id.sector(),
                id.superLayer(),
                id.layer(),
                0,     fCell, lCell,
                flagA, flagC, flagS );
  if ( ( fCheck == 0 ) &&
       ( fCell <= iCell ) &&
       ( lCell >= iCell ) ) return 0;
  fCheck = get( id.wheel(),
                id.station(),
                id.sector(),
                id.superLayer(),
                id.layer(),
                1,     fCell, lCell,
                flagA, flagC, flagS );
  if ( ( fCheck == 0 ) &&
       ( fCell <= iCell ) &&
       ( lCell >= iCell ) ) return 0;
  flagA = flagC = flagS = 0;
  return 1;
}


int DTHVStatus::offChannelsNumber() const {
  int offNum = 0;
  DTHVStatus::const_iterator iter = begin();
  DTHVStatus::const_iterator iend = end();
  while ( iter != iend ) {
    const std::pair<DTHVStatusId,DTHVStatusData>& entry = *iter++;
    DTHVStatusId   hvId = entry.first;
    DTHVStatusData data = entry.second;
    int offA = data.flagA & 1;
    int offC = data.flagC & 1;
    int offS = data.flagS & 1;
    if ( offA || offC || offS )
         offNum += ( 1 + data.lCell - data.fCell );
  }
  return offNum;
}


int DTHVStatus::offChannelsNumber( const DTChamberId& id ) const {
  int offNum = 0;
  DTHVStatus::const_iterator iter = begin();
  DTHVStatus::const_iterator iend = end();
  while ( iter != iend ) {
    const std::pair<DTHVStatusId,DTHVStatusData>& entry = *iter++;
    DTHVStatusId   hvId = entry.first;
    DTHVStatusData data = entry.second;
    if ( hvId.  wheelId != id.  wheel() ) continue;
    if ( hvId.stationId != id.station() ) continue;
    if ( hvId. sectorId != id. sector() ) continue;
    int offA = data.flagA & 1;
    int offC = data.flagC & 1;
    int offS = data.flagS & 1;
    if ( offA || offC || offS )
         offNum += ( 1 + data.lCell - data.fCell );
  }
  return offNum;
}


int DTHVStatus::badChannelsNumber() const {
  int offNum = 0;
  DTHVStatus::const_iterator iter = begin();
  DTHVStatus::const_iterator iend = end();
  while ( iter != iend ) {
    const std::pair<DTHVStatusId,DTHVStatusData>& entry = *iter++;
    DTHVStatusId   hvId = entry.first;
    DTHVStatusData data = entry.second;
    if ( data.flagA || data.flagC || data.flagS )
         offNum += ( 1 + data.lCell - data.fCell );
  }
  return offNum;
}


int DTHVStatus::badChannelsNumber( const DTChamberId& id ) const {
  int offNum = 0;
  DTHVStatus::const_iterator iter = begin();
  DTHVStatus::const_iterator iend = end();
  while ( iter != iend ) {
    const std::pair<DTHVStatusId,DTHVStatusData>& entry = *iter++;
    DTHVStatusId   hvId = entry.first;
    DTHVStatusData data = entry.second;
    if ( hvId.  wheelId != id.  wheel() ) continue;
    if ( hvId.stationId != id.station() ) continue;
    if ( hvId. sectorId != id. sector() ) continue;
    if ( data.flagA || data.flagC || data.flagS )
         offNum += ( 1 + data.lCell - data.fCell );
  }
  return offNum;
}


const
std::string& DTHVStatus::version() const {
  return dataVersion;
}


std::string& DTHVStatus::version() {
  return dataVersion;
}


void DTHVStatus::clear() {
  dataList.clear();
  dataList.reserve( 10 );
  initialize();
  return;
}


int DTHVStatus::set( int   wheelId,
                     int stationId,
                     int  sectorId,
                     int      slId,
                     int   layerId,
                     int    partId,
                     int     fCell,
                     int     lCell,
                     int     flagA,
                     int     flagC,
                     int     flagS ) {

  std::vector<int> chanKey;
  chanKey.reserve(6);
  chanKey.push_back(   wheelId );
  chanKey.push_back( stationId );
  chanKey.push_back(  sectorId );
  chanKey.push_back(      slId );
  chanKey.push_back(   layerId );
  chanKey.push_back(    partId );
  int ientry;
  int searchStatus = dBuf->find( chanKey.begin(), chanKey.end(), ientry );

  if ( !searchStatus ) {
    DTHVStatusData& data( dataList[ientry].second );
    data.fCell = fCell;
    data.lCell = lCell;
    data.flagA = flagA;
    data.flagC = flagC;
    data.flagS = flagS;
    return -1;
  }
  else {
    DTHVStatusId key;
    key.  wheelId =   wheelId;
    key.stationId = stationId;
    key. sectorId =  sectorId;
    key.     slId =      slId;
    key.  layerId =   layerId;
    key.   partId =    partId;
    DTHVStatusData data;
    data.fCell = fCell;
    data.lCell = lCell;
    data.flagA = flagA;
    data.flagC = flagC;
    data.flagS = flagS;
    ientry = dataList.size();
    dataList.push_back( std::pair<DTHVStatusId,
                                  DTHVStatusData>( key, data ) );
    dBuf->insert( chanKey.begin(), chanKey.end(), ientry );
    return 0;
  }

  return 99;

}


int DTHVStatus::set( const DTLayerId& id,
                     int    partId,
                     int     fCell,
                     int     lCell,
                     int     flagA,
                     int     flagC,
                     int     flagS ) {
  return set( id.wheel(),
              id.station(),
              id.sector(),
              id.superLayer(),
              id.layer(),
              partId,
              fCell, lCell, flagA, flagC, flagS );
}


int DTHVStatus::setFlagA( int   wheelId,
                          int stationId,
                          int  sectorId,
                          int      slId,
                          int   layerId,
                          int    partId,
                          int      flag ) {
  int fCell;
  int lCell;
  int flagA;
  int flagC;
  int flagS;
  get(   wheelId,
       stationId,
        sectorId,
            slId,
         layerId,
          partId,
           fCell,
           lCell,
           flagA,
           flagC,
           flagS );
  return set(   wheelId,
              stationId,
               sectorId,
                   slId,
                layerId,
                 partId,
                  fCell,
                  lCell,
                   flag,
                  flagC,
                  flagS );
}


int DTHVStatus::setFlagA( const DTLayerId& id,
                          int    partId,
                          int      flag ) {
  return setFlagA( id.wheel(),
                   id.station(),
                   id.sector(),
                   id.superLayer(),
                   id.layer(),
                   partId,
                   flag );
}


int DTHVStatus::setFlagC( int   wheelId,
                          int stationId,
                          int  sectorId,
                          int      slId,
                          int   layerId,
                          int    partId,
                          int      flag ) {
  int fCell;
  int lCell;
  int flagA;
  int flagC;
  int flagS;
  get(   wheelId,
       stationId,
        sectorId,
            slId,
         layerId,
          partId,
           fCell,
           lCell,
           flagA,
           flagC,
           flagS );
  return set(   wheelId,
              stationId,
               sectorId,
               slId,
                layerId,
                 partId,
                  fCell,
                  lCell,
                  flagA,
                   flag,
                  flagS );
}


int DTHVStatus::setFlagC( const DTLayerId& id,
                          int    partId,
                          int      flag ) {
  return setFlagC( id.wheel(),
                   id.station(),
                   id.sector(),
                   id.superLayer(),
                   id.layer(),
                   partId,
                   flag );
}


int DTHVStatus::setFlagS( int   wheelId,
                          int stationId,
                          int  sectorId,
                          int      slId,
                          int   layerId,
                          int    partId,
                          int      flag ) {
  int fCell;
  int lCell;
  int flagA;
  int flagC;
  int flagS;
  get(   wheelId,
       stationId,
        sectorId,
            slId,
         layerId,
          partId,
           fCell,
           lCell,
           flagA,
           flagC,
           flagS );
  return set(   wheelId,
              stationId,
               sectorId,
                   slId,
                layerId,
                 partId,
                  fCell,
                  lCell,
                  flagA,
                  flagC,
                   flag );
}


int DTHVStatus::setFlagS( const DTLayerId& id,
                          int    partId,
                          int      flag ) {
  return setFlagS( id.wheel(),
                   id.station(),
                   id.sector(),
                   id.superLayer(),
                   id.layer(),
                   partId,
                   flag );
}


DTHVStatus::const_iterator DTHVStatus::begin() const {
  return dataList.begin();
}


DTHVStatus::const_iterator DTHVStatus::end() const {
  return dataList.end();
}


std::string DTHVStatus::mapName() const {
  std::stringstream name;
  name << dataVersion << "_map_HV" << this;
  return name.str();
}


void DTHVStatus::initialize() {

  dBuf->clear();

  int entryNum = 0;
  int entryMax = dataList.size();
  std::vector<int> chanKey;
  chanKey.reserve(6);
  while ( entryNum < entryMax ) {

    const DTHVStatusId& chan = dataList[entryNum].first;

    chanKey.clear();
    chanKey.push_back( chan.  wheelId );
    chanKey.push_back( chan.stationId );
    chanKey.push_back( chan. sectorId );
    chanKey.push_back( chan.     slId );
    chanKey.push_back( chan.  layerId );
    chanKey.push_back( chan.   partId );
    dBuf->insert( chanKey.begin(), chanKey.end(), entryNum++ );

  }

  return;

}