/data/refman/pasoursint/CMSSW_5_3_3/src/CaloOnlineTools/HcalOnlineDb/src/HcalChannelIterator.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:     HcalOnlineDb
// Class  :     HcalChannelIterator
// 
// Implementation:
//     <Notes on implementation>
//
// Original Author:  Gena Kukartsev
//         Created:  Mon Jul 13 12:15:33 CEST 2009
// $Id: HcalChannelIterator.cc,v 1.5 2010/08/06 20:24:11 wmtan Exp $
//

#include <fstream>
#include "CaloOnlineTools/HcalOnlineDb/interface/HcalChannelIterator.h"
#include "CaloOnlineTools/HcalOnlineDb/interface/RooGKCounter.h"
#include "DataFormats/HcalDetId/interface/HcalDetId.h"

HcalChannelIterator::HcalChannelIterator()
{
}


HcalChannelIterator::~HcalChannelIterator()
{
}



int HcalChannelIterator::size(void){
  return channel_list.size();
}


int HcalChannelIterator::clearChannelList(void){
  channel_list.clear();
  return 0;
}


int HcalChannelIterator::addListFromLmapAscii(std::string filename){
  RooGKCounter lines;
  int _current_size = size();
  std::string _row;
  std::ifstream inFile( filename . c_str(), std::ios::in );
  if (!inFile){
    std::cout << "Unable to open file with the logical map: " << filename << std::endl;
  }
  else{
    std::cout << "File with the logical map opened successfully: " << filename << std::endl;
  }
  while ( getline( inFile, _row ) > 0 ){
    //#   side    eta    phi   dphi  depth    det
    int _num, _side, _eta, _phi, _dphi, _depth;
    char subdetbuf[32];
    
    int _read;
    const char * _format = "%d %d %d %d %d %d %s";
    _read = sscanf( _row . c_str(), _format,
                    &_num, &_side, &_eta, &_phi, &_dphi, &_depth,
                    subdetbuf );
    if ( _read == 7 ){
      lines . count();
      
      std::string subdet(subdetbuf);

      HcalSubdetector _det;
      if ( subdet.find("HB")!=std::string::npos ) _det = HcalBarrel;
      else if ( subdet.find("HE")!=std::string::npos ) _det = HcalEndcap;
      else if ( subdet.find("HO")!=std::string::npos ) _det = HcalOuter;
      else if ( subdet.find("HF")!=std::string::npos ) _det = HcalForward;
      else _det = HcalOther;

      HcalDetId _detid(_det, _side*_eta, _phi, _depth);
      
      if (_det == HcalBarrel || _det == HcalEndcap || _det == HcalOuter || _det == HcalForward){
        channel_list . push_back( _detid );
      }
    }  
  }
  inFile.close();
  std::cout << "Logical map file: " << lines . getCount() << " lines read" << std::endl;
  std::cout << "Logical map file: " << size() - _current_size << " lines added to the list" << std::endl;
  //
  return 0;
}


int HcalChannelIterator::begin(void){
  const_iterator = channel_list.begin();
  return 0;
}


int HcalChannelIterator::next(void){
  const_iterator++;
  return 0;
}


bool HcalChannelIterator::end(void){
  if (const_iterator==channel_list.end()){
    return true;
  }
  else{
    return false;
  }
}


HcalGenericDetId HcalChannelIterator::getHcalGenericDetId(void){
  if (const_iterator!=channel_list.end()){
    return *const_iterator;
  }
  else{
    return 0;
  }
}


HcalSubdetector HcalChannelIterator::getHcalSubdetector(void){
  if (const_iterator!=channel_list.end()){
    HcalDetId _id(*const_iterator);
    return _id.subdet();
  }
  else return HcalOther;
}
 

int HcalChannelIterator::getIeta(void){
  if (const_iterator!=channel_list.end()){
    HcalDetId _id(*const_iterator);
    return _id.ieta();
  }
  else return -1000;
}


int HcalChannelIterator::getIphi(void){
  if (const_iterator!=channel_list.end()){
    HcalDetId _id(*const_iterator);
    return _id.iphi();
  }
  else return -1000;
}


int HcalChannelIterator::getDepth(void){
  if (const_iterator!=channel_list.end()){
    HcalDetId _id(*const_iterator);
    return _id.depth();
  }
  else return -1000;
}


int HcalChannelIterator::initHBEFListFromLmapAscii(void){
  clearChannelList();
  addListFromLmapAscii("HCALmapHBEF_Jan.27.2009.txt");
  addListFromLmapAscii("HCALmapHO_Jan.27.2009.txt");
  return channel_list.size();
}


int HcalChannelIterator::init(std::vector<HcalGenericDetId> map){
  channel_list.clear();
  channel_list = map;
  return channel_list.size();
}