EcalDQMStatusWriter Class Reference

#include <EcalDQMStatusWriter.h>

Inheritance diagram for EcalDQMStatusWriter:

Public Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c)
void	beginRun (const edm::Run &r, const edm::EventSetup &c)
	EcalDQMStatusWriter (const edm::ParameterSet &ps)
virtual	~EcalDQMStatusWriter ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Private Member Functions
uint32_t	convert (uint32_t c)
EcalDQMChannelStatus *	readEcalDQMChannelStatusFromFile (const char *)
EcalDQMTowerStatus *	readEcalDQMTowerStatusFromFile (const char *)

Private Attributes
std::vector< std::string >	inpFileName_
std::vector< std::string >	objectName_
std::vector< unsigned long long >	since_
bool	verbose_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 29 of file EcalDQMStatusWriter.h.

Constructor & Destructor Documentation

EcalDQMStatusWriter::EcalDQMStatusWriter

(

const edm::ParameterSet &

ps

)

Definition at line 31 of file EcalDQMStatusWriter.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), inpFileName_, objectName_, since_, convertXMLtoSQLite_cfg::toPut, and verbose_.

                                                                  {

  verbose_ = ps.getUntrackedParameter<bool>("verbose", false);

  typedef std::vector<edm::ParameterSet> Parameters;
  Parameters toPut = ps.getParameter<Parameters>("toPut");

  for ( Parameters::iterator itToPut=toPut.begin(); itToPut!=toPut.end(); itToPut++ ) {
    inpFileName_.push_back(itToPut->getUntrackedParameter<std::string>("inputFile"));
    objectName_.push_back(itToPut->getUntrackedParameter<std::string>("conditionType"));
    since_.push_back(itToPut->getUntrackedParameter<unsigned int>("since"));
  }

}

virtual EcalDQMStatusWriter::~EcalDQMStatusWriter ( )

inlinevirtual

Definition at line 34 of file EcalDQMStatusWriter.h.

{};

Member Function Documentation

void EcalDQMStatusWriter::analyze ( const edm::Event &  e, const edm::EventSetup &  c  )

inlinevirtual

Implements edm::EDAnalyzer.

Definition at line 36 of file EcalDQMStatusWriter.h.

{};

void EcalDQMStatusWriter::beginRun ( const edm::Run &  r, const edm::EventSetup &  c  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 46 of file EcalDQMStatusWriter.cc.

References cond::service::PoolDBOutputService::appendSinceTime(), cond::service::PoolDBOutputService::beginOfTime(), gather_cfg::cout, cond::service::PoolDBOutputService::createNewIOV(), cond::service::PoolDBOutputService::endOfTime(), i, Numbers::initGeometry(), inpFileName_, edm::Service< T >::isAvailable(), cond::service::PoolDBOutputService::isNewTagRequest(), objectName_, readEcalDQMChannelStatusFromFile(), readEcalDQMTowerStatusFromFile(), since_, ntuplemaker::status, and verbose_.

                                                                          {

  Numbers::initGeometry(c, verbose_);

  edm::Service<cond::service::PoolDBOutputService> dbservice;
  if ( !dbservice.isAvailable() ){
    std::cout << "PoolDBOutputService is unavailable" << std::endl;
    return;
  }

  for ( unsigned int i=0; i<objectName_.size(); i++ ) {

    bool toAppend;
    cond::Time_t newTime;

    if ( dbservice->isNewTagRequest( objectName_[i]+std::string("Rcd") ) ) {
      // This is the first object for this tag.
      // Append mode should be off.
      // newTime is the end of this new objects IOV.
      toAppend = false;
      newTime = dbservice->beginOfTime();
    } else {
      // There should already be an object in the DB for this tag.
      // Append IOV mode should be on.
      // newTime is the beginning of this new objects IOV.
      toAppend = true;
      newTime = (cond::Time_t) since_[i];
    }

    std::cout << "Reading " << objectName_[i] << " from file and writing to DB with newTime " << newTime << std::endl;

    if (objectName_[i]  ==  "EcalDQMChannelStatus") {

      EcalDQMChannelStatus* status = readEcalDQMChannelStatusFromFile(inpFileName_[i].c_str());

      if ( !toAppend ) {
        dbservice->createNewIOV<EcalDQMChannelStatus>(status, newTime, dbservice->endOfTime(), "EcalDQMChannelStatusRcd");
      } else {
        dbservice->appendSinceTime<EcalDQMChannelStatus>(status, newTime, "EcalDQMChannelStatusRcd");
      }

    } else if (objectName_[i]  ==  "EcalDQMTowerStatus") {

      EcalDQMTowerStatus* status = readEcalDQMTowerStatusFromFile(inpFileName_[i].c_str());

      if ( !toAppend ) {
        dbservice->createNewIOV<EcalDQMTowerStatus>(status, newTime, dbservice->endOfTime(), "EcalDQMTowerStatusRcd");
      } else {
        dbservice->appendSinceTime<EcalDQMTowerStatus>(status, newTime, "EcalDQMTowerStatusRcd");
      }

    } else {

      std::cout << "Object " << objectName_[i]  << " is not supported by this program." << std::endl;

    }

  }

}

uint32_t EcalDQMStatusWriter::convert ( uint32_t  c )

private

Definition at line 518 of file EcalDQMStatusWriter.cc.

References EcalDQMStatusHelper::LASER_MEAN_ERROR, EcalDQMStatusHelper::LASER_RMS_ERROR, EcalDQMStatusHelper::LASER_TIMING_MEAN_ERROR, EcalDQMStatusHelper::LASER_TIMING_RMS_ERROR, EcalDQMStatusHelper::LED_MEAN_ERROR, EcalDQMStatusHelper::LED_RMS_ERROR, EcalDQMStatusHelper::LED_TIMING_MEAN_ERROR, EcalDQMStatusHelper::LED_TIMING_RMS_ERROR, EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_MEAN_ERROR, EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_RMS_ERROR, EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_MEAN_ERROR, EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_RMS_ERROR, EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_MEAN_ERROR, EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_RMS_ERROR, EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_MEAN_ERROR, EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR, EcalDQMStatusHelper::PHYSICS_BAD_CHANNEL_WARNING, EcalDQMStatusHelper::STATUS_FLAG_ERROR, EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_MEAN_ERROR, EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_RMS_ERROR, EcalDQMStatusHelper::TESTPULSE_LOW_GAIN_MEAN_ERROR, EcalDQMStatusHelper::TESTPULSE_LOW_GAIN_RMS_ERROR, EcalDQMStatusHelper::TESTPULSE_MIDDLE_GAIN_MEAN_ERROR, and EcalDQMStatusHelper::TESTPULSE_MIDDLE_GAIN_RMS_ERROR.

Referenced by readEcalDQMChannelStatusFromFile().

                                                       {

  if ( chStatus == 1 || (chStatus >= 8 && chStatus <= 12 )) {

    chStatus = 0;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_LOW_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_MIDDLE_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_LOW_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_MIDDLE_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LASER_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LASER_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LASER_TIMING_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LASER_TIMING_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LED_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LED_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LED_TIMING_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LED_TIMING_RMS_ERROR;

  } else if ( chStatus == 2 ) {

    chStatus = 0;
    chStatus |= 1 << EcalDQMStatusHelper::LASER_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LASER_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LASER_TIMING_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LASER_TIMING_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LED_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LED_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LED_TIMING_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LED_TIMING_RMS_ERROR;

  } else if ( chStatus == 3 ) {

    chStatus = 0;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PHYSICS_BAD_CHANNEL_WARNING;

  } else if ( chStatus == 4 ) {

    chStatus = 0;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR;

  } else if ( chStatus == 13 || chStatus == 14 ) {

    chStatus = 0;
    chStatus |= 1 << EcalDQMStatusHelper::STATUS_FLAG_ERROR;

  } else if ( (chStatus&0x3f) == 13 || (chStatus&0x7f) == 13 ||
              (chStatus&0x3f) == 14 || (chStatus&0x7f) == 14 ) {

    chStatus = 0;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_LOW_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_MIDDLE_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_LOW_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_MIDDLE_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LASER_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LASER_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LASER_TIMING_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LASER_TIMING_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LED_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LED_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LED_TIMING_MEAN_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::LED_TIMING_RMS_ERROR;
    chStatus |= 1 << EcalDQMStatusHelper::STATUS_FLAG_ERROR;

  }

  return( chStatus );

}

EcalDQMChannelStatus * EcalDQMStatusWriter::readEcalDQMChannelStatusFromFile ( const char *  inputFile )

private

Definition at line 107 of file EcalDQMStatusWriter.cc.

References convert(), gather_cfg::cout, EcalBarrel, EcalCondObjectContainer< T >::end(), edm::hlt::Exception, EcalCondObjectContainer< T >::find(), EcalDQMStatusDictionary::getDictionary(), i, errorMatrix2Lands_multiChannel::id, compare_using_db::ifile, getHLTprescales::index, ecalpyutils::ism(), Numbers::iSM(), combine::key, geometryCSVtoXML::line, python.rootplot.argparse::module, EcalCondObjectContainer< T >::setValue(), EBDetId::SMCRYSTALMODE, ntuplemaker::status, EBDetId::unhashIndex(), EEDetId::validDetId(), and verbose_.

Referenced by beginRun().

                                                                                                 {

  EcalDQMChannelStatus* status = new EcalDQMChannelStatus();

  std::vector<EcalDQMStatusDictionary::codeDef> dictionary;
  EcalDQMStatusDictionary::getDictionary( dictionary );

  // barrel
  for ( int ism=1; ism<=36; ism++ ) {
    int jsm = Numbers::iSM(ism, EcalBarrel);
    for ( int ic=1; ic<=1700; ic++ ) {
      EBDetId id(jsm, ic, EBDetId::SMCRYSTALMODE);
      status->setValue(id, 0);
    }
  }

  // endcap
  for ( int ix=1; ix<=100; ix++ ) {
    for ( int iy=1; iy<=100; iy++ ) {
      if ( EEDetId::validDetId(ix, iy, +1) ) {
        EEDetId id(ix, iy, +1);
        status->setValue(id, 0);
      }
      if ( EEDetId::validDetId(ix, iy, -1) ) {
        EEDetId id(ix, iy, -1);
        status->setValue(id, 0);
      }
    }
  }

  std::cout << "Reading channel status from file " << inputFile << std::endl;
  FILE *ifile = fopen( inputFile ,"r" );

  if ( !ifile ) throw cms::Exception ("Cannot open file") ;

  char line[256];

  int ii = 0;
  while ( fgets(line, 255, ifile) ) {

    std::stringstream stream;

    ii++;
    std::string key;
    stream << line;
    if ( verbose_ ) std::cout << line;
    stream >> key;

    if ( key.size() == 0 || strcmp(key.c_str(), " ") == 0 || strncmp(key.c_str(), "#", 1) == 0 ) {

      // skip comments

    } else if ( strcmp(key.c_str(), "EB") == 0 ) {

      int index;
      uint32_t code;
      stream >> index >> code;

      EBDetId id = EBDetId::unhashIndex(index);
      code = convert(code);

      EcalDQMChannelStatus::const_iterator it = status->find(id);
      if ( it != status->end() ) code |= it->getStatusCode();
      status->setValue(id, code);

    } else if ( strcmp(key.c_str(), "EE") == 0 ) {

      int ix, iy, iz;
      uint32_t code;
      stream >> ix >> iy >> iz >> code;

      EEDetId id(ix, iy, iz);
      code = convert(code);

      EcalDQMChannelStatus::const_iterator it = status->find(id);
      if ( it != status->end() ) code |= it->getStatusCode();
      status->setValue(id, code);

    } else if ( strcmp(key.c_str(), "Crystal") == 0 ) {

      std::string module;
      stream >> module;

      if ( strncmp(module.c_str(), "EB+", 3) == 0 || strncmp(module.c_str(), "EB-", 3) == 0 ) {

#if 0
        int ie, ip;
        std::string token;
        stream >> ie >> ip >> token;

        int sm = atoi( module.substr(2, module.size()-2).c_str() );

        int ism = (sm>=1&&sm<=18) ? sm+18 : -sm;
        int jsm = Numbers::iSM(ism, EcalBarrel);
        int ic = 20*(ie-1)+(ip-1)+1;
#else
        int ic;
        std::string token;
        stream >> ic >> token;

        int sm = atoi( module.substr(2, module.size()-2).c_str() );

        int ism = (sm>=1&&sm<=18) ? sm+18 : -sm;
        int jsm = Numbers::iSM(ism, EcalBarrel);
#endif

        EBDetId id(jsm, ic, EBDetId::SMCRYSTALMODE);
        uint32_t code = 0;
        for ( unsigned int i=0; i<dictionary.size(); i++ ) {
          if ( strcmp(token.c_str(), dictionary[i].desc) == 0 ) {
            code = dictionary[i].code;
          }
        }
        if ( code == 0 ) {
          std::cout << " --> not found in the dictionary: " << token << std::endl;
          continue;
        }

        EcalDQMChannelStatus::const_iterator it = status->find(id);
        if ( it != status->end() ) code |= it->getStatusCode();
        status->setValue(id, code);

      } else if ( strncmp(module.c_str(), "EE+", 3) == 0 || strncmp(module.c_str(), "EE-", 3) == 0 ) {

#if 0
        int jx, jy;
        std::string token;
        stream >> jx >> jy >> token;
#else
        int index;
        std::string token;
        stream >> index >> token;
        int jx = index/1000;
        int jy = index%1000;
#endif

        int sm = atoi( module.substr(2, module.size()-2).c_str() );

        int ism = 0;
        if( sm == -99 ) sm = -1;
        if( sm == +99 ) sm = +1;
        switch ( sm ) {
          case -7: ism =  1; break;
          case -8: ism =  2; break;
          case -9: ism =  3; break;
          case -1: ism =  4; break;
          case -2: ism =  5; break;
          case -3: ism =  6; break;
          case -4: ism =  7; break;
          case -5: ism =  8; break;
          case -6: ism =  9; break;
          case +7: ism = 10; break;
          case +8: ism = 11; break;
          case +9: ism = 12; break;
          case +1: ism = 13; break;
          case +2: ism = 14; break;
          case +3: ism = 15; break;
          case +4: ism = 16; break;
          case +5: ism = 17; break;
          case +6: ism = 18; break;
        }

        EEDetId id(jx, jy, (ism>=1&&ism<=9)?-1:+1);
        uint32_t code = 0;
        for ( unsigned int i=0; i<dictionary.size(); i++ ) {
          if ( strcmp(token.c_str(), dictionary[i].desc) == 0 ) {
            code = dictionary[i].code;
          }
        }
        if ( code == 0 ) {
          std::cout << " --> not found in the dictionary: " << token << std::endl;
          continue;
        }

        EcalDQMChannelStatus::const_iterator it = status->find(id);
        if ( it != status->end() ) code |= it->getStatusCode();
        status->setValue(id, code);

      } else {

        std:: cout << "--> unknown token at line #" << ii << " : " << line;

      }

    } else if ( strcmp(key.c_str(), "TT") == 0 ) {

      // see readEcalDQMTowerStatusFromFile

    } else if ( strcmp(key.c_str(), "PN") == 0 || strcmp(key.c_str(), "MemCh") == 0 || strcmp(key.c_str(), "MemTT") == 0 ) {

      std::cout << "--> unsupported key at line #" << ii << " : " << line;

    } else {

      std:: cout << "--> skipped line #" << ii << " : " << line;

    }

  }

  fclose(ifile);

  return status;

}

EcalDQMTowerStatus * EcalDQMStatusWriter::readEcalDQMTowerStatusFromFile ( const char *  inputFile )

private

Definition at line 313 of file EcalDQMStatusWriter.cc.

References abs, gather_cfg::cout, Numbers::crystals(), EcalBarrel, EcalCondTowerObjectContainer< T >::end(), edm::hlt::Exception, EcalCondTowerObjectContainer< T >::find(), EcalDQMStatusDictionary::getDictionary(), i, errorMatrix2Lands_multiChannel::id, compare_using_db::ifile, ecalpyutils::ism(), combine::key, geometryCSVtoXML::line, python.rootplot.argparse::module, EcalCondTowerObjectContainer< T >::setValue(), ntuplemaker::status, EcalTrigTowerDetId::validDetId(), EcalScDetId::validDetId(), and verbose_.

Referenced by beginRun().

                                                                                             {

  EcalDQMTowerStatus* status = new EcalDQMTowerStatus();

  std::vector<EcalDQMStatusDictionary::codeDef> dictionary;
  EcalDQMStatusDictionary::getDictionary( dictionary );

  // barrel
  for ( int ix=1; ix<=17; ix++ ) {
    for ( int iy=1; iy<=72; iy++ ) {
      if ( EcalTrigTowerDetId::validDetId(+1, EcalBarrel, ix, iy) ) {
        EcalTrigTowerDetId id(+1, EcalBarrel, ix, iy);
        status->setValue(id, 0);
      }
      if ( EcalTrigTowerDetId::validDetId(-1, EcalBarrel, ix, iy) ) {
        EcalTrigTowerDetId id(-1, EcalBarrel, ix, iy);
        status->setValue(id, 0);
      }
    }
  }
  
  // endcap
  for ( int ix=1; ix<=20; ix++ ) {
    for ( int iy=1; iy<=20; iy++ ) {
      if ( EcalScDetId::validDetId(ix, iy, +1) ) {
        EcalScDetId id(ix, iy, +1);
        status->setValue(id, 0);
      }
      if ( EcalScDetId::validDetId(ix, iy, -1) ) {
        EcalScDetId id(ix, iy, -1);
        status->setValue(id, 0);
      }
    }
  }

  std::cout << "Reading tower status from file " << inputFile << std::endl;
  FILE *ifile = fopen( inputFile ,"r" );

  if ( !ifile ) throw cms::Exception ("Cannot open file") ;

  char line[256];

  int ii = 0;
  while ( fgets(line, 255, ifile) ) {

    std::stringstream stream;

    ii++;
    std::string key;
    stream << line;
    if ( verbose_ ) std::cout << line;
    stream >> key;

    if ( key.size() == 0 || strcmp(key.c_str(), " ") == 0 || strncmp(key.c_str(), "#", 1) == 0 ) {

      // skip comments

    } else if ( strcmp(key.c_str(), "EB") == 0 ) {

      // see readEcalDQMChannelStatusFromFile

    } else if ( strcmp(key.c_str(), "EE") == 0 ) {

      // see readEcalDQMChannelStatusFromFile

    } else if ( strcmp(key.c_str(), "Crystal") == 0 ) {

      // see readEcalDQMChannelStatusFromFile

    } else if ( strcmp(key.c_str(), "TT") == 0 ) {

      std::string module;
      stream >> module;

      if ( strncmp(module.c_str(), "EB+", 3) == 0 || strncmp(module.c_str(), "EB-", 3) == 0 ) {

        int itt;
        std::string token;
        stream >> itt >> token;

        if ( itt >= 1 && itt <= 68 ) {

          int sm = atoi( module.substr(2, module.size()-2).c_str() );

          int iet = (itt-1)/4+1;
          int ipt = (itt-1)%4+1;

          if ( sm<0 ) {
            ipt = ipt+(std::abs(sm)-1)*4-2;
            if ( ipt < 1 ) ipt = ipt+72;
            if ( ipt > 72 ) ipt = ipt-72;
          } else {
            ipt = (5-ipt)+(std::abs(sm)-1)*4-2;
            if ( ipt < 1 ) ipt = ipt+72;
            if ( ipt > 72 ) ipt = ipt-72;
          }

          EcalTrigTowerDetId id((sm<0)?-1:+1, EcalBarrel, iet, ipt);
          uint32_t code = 0;
          for ( unsigned int i=0; i<dictionary.size(); i++ ) {
            if ( strcmp(token.c_str(), dictionary[i].desc) == 0 ) {
              code = dictionary[i].code;
            }
          }
          if ( code == 0 ) {
            std::cout << " --> not found in the dictionary: " << token << std::endl;
            continue;
          }

          EcalDQMTowerStatus::const_iterator it = status->find(id);
          if ( it != status->end() ) code |= it->getStatusCode();
          status->setValue(id, code);

        } else {

          std::cout << "--> unsupported configuration at line #" << ii << " : " << line;

        }

      } else if ( strncmp(module.c_str(), "EE+", 3) == 0 || strncmp(module.c_str(), "EE-", 3) == 0 ) {

        int isc;
        std::string token;
        stream >> isc >> token;

        if ( isc >= 1 && isc <= 68 ) {

          int sm = atoi( module.substr(2, module.size()-2).c_str() );

          int ism = 0;
          switch ( sm ) {
            case -7: ism =  1; break;
            case -8: ism =  2; break;
            case -9: ism =  3; break;
            case -1: ism =  4; break;
            case -2: ism =  5; break;
            case -3: ism =  6; break;
            case -4: ism =  7; break;
            case -5: ism =  8; break;
            case -6: ism =  9; break;
            case +7: ism = 10; break;
            case +8: ism = 11; break;
            case +9: ism = 12; break;
            case +1: ism = 13; break;
            case +2: ism = 14; break;
            case +3: ism = 15; break;
            case +4: ism = 16; break;
            case +5: ism = 17; break;
            case +6: ism = 18; break;
          }

          int idcc = (ism>=1&&ism<=9) ? ism : ism-9+45;

          std::vector<DetId>* crystals = Numbers::crystals(idcc, isc);

          for ( unsigned int i=0; i<crystals->size(); i++ ) {

            EcalScDetId id = ((EEDetId) (*crystals)[i]).sc();
            uint32_t code = 0;
            for ( unsigned int i=0; i<dictionary.size(); i++ ) {
              if ( strcmp(token.c_str(), dictionary[i].desc) == 0 ) {
                code = dictionary[i].code;
              }
            }
            if ( code == 0 ) {
              std::cout << " --> not found in the dictionary: " << token << std::endl;
              continue;
            }

            EcalDQMTowerStatus::const_iterator it = status->find(id);
            if ( it != status->end() ) code |= it->getStatusCode();
            status->setValue(id, code);

          }

        } else {

          std::cout << "--> unsupported configuration at line #" << ii << " : " << line;

        }

      } else {

        std:: cout << "--> unknown token at line #" << ii << " : " << line;

      }

    } else if ( strcmp(key.c_str(), "PN") == 0 || strcmp(key.c_str(), "MemCh") == 0 || strcmp(key.c_str(), "MemTT") == 0 ) {

      std::cout << "--> unsupported key at line #" << ii << " : " << line;

    } else {

      std:: cout << "--> skipped line #" << ii << " : " << line;

    }

  }

  fclose(ifile);

  return status;

}

Member Data Documentation

std::vector<std::string> EcalDQMStatusWriter::inpFileName_

private

Definition at line 49 of file EcalDQMStatusWriter.h.

Referenced by beginRun(), and EcalDQMStatusWriter().

std::vector<std::string> EcalDQMStatusWriter::objectName_

private

Definition at line 48 of file EcalDQMStatusWriter.h.

Referenced by beginRun(), and EcalDQMStatusWriter().

std::vector<unsigned long long> EcalDQMStatusWriter::since_

private

Definition at line 50 of file EcalDQMStatusWriter.h.

Referenced by beginRun(), and EcalDQMStatusWriter().

bool EcalDQMStatusWriter::verbose_

private

Definition at line 47 of file EcalDQMStatusWriter.h.

Referenced by beginRun(), EcalDQMStatusWriter(), readEcalDQMChannelStatusFromFile(), and readEcalDQMTowerStatusFromFile().