/data/refman/pasoursint/CMSSW_4_4_5_patch3/src/CondTools/SiStrip/plugins/SiStripBadFiberBuilder.cc

Go to the documentation of this file.

#include "CondTools/SiStrip/plugins/SiStripBadFiberBuilder.h"

#include "CalibTracker/SiStripCommon/interface/SiStripDetInfoFileReader.h"

#include <iostream>
#include <fstream>
#include <sstream>

SiStripBadFiberBuilder::SiStripBadFiberBuilder(const edm::ParameterSet& iConfig) : ConditionDBWriter<SiStripBadStrip>(iConfig){

  fp_ = iConfig.getUntrackedParameter<edm::FileInPath>("file",edm::FileInPath("CalibTracker/SiStripCommon/data/SiStripDetInfo.dat"));
  printdebug_ = iConfig.getUntrackedParameter<bool>("printDebug",false);
  BadComponentList_ =  iConfig.getUntrackedParameter<Parameters>("BadComponentList");
}


SiStripBadFiberBuilder::~SiStripBadFiberBuilder(){
}

void SiStripBadFiberBuilder::algoAnalyze(const edm::Event & event, const edm::EventSetup& iSetup){
  
  unsigned int run=event.id().run();

  edm::LogInfo("SiStripBadFiberBuilder") << "... creating dummy SiStripBadStrip Data for Run " << run << "\n " << std::endl;
  
  obj = new SiStripBadStrip();

  SiStripDetInfoFileReader reader(fp_.fullPath());
  
  const std::vector<uint32_t> DetIds = reader.getAllDetIds();

  std::stringstream ss;
  for(Parameters::iterator iBadComponent = BadComponentList_.begin(); iBadComponent != BadComponentList_.end(); ++iBadComponent ) {
    
    uint32_t BadModule_ = iBadComponent->getParameter<uint32_t>("BadModule");
    std::vector<uint32_t> BadApvList_ = iBadComponent->getParameter<std::vector<uint32_t> >("BadApvList");

    std::vector<unsigned int> theSiStripVector;
    
    unsigned short firstBadStrip=0, NconsecutiveBadStrips=0;
    unsigned int theBadStripRange;
    
    for(std::vector<uint32_t>::const_iterator is=BadApvList_.begin(); is!=BadApvList_.end(); ++is){

      firstBadStrip=(*is)*128;
      NconsecutiveBadStrips=128;
      
      theBadStripRange = obj->encode(firstBadStrip,NconsecutiveBadStrips);
      
      if (printdebug_)
        ss << "detid " << BadModule_ << " \t"
           << " firstBadStrip " << firstBadStrip << "\t "
           << " NconsecutiveBadStrips " << NconsecutiveBadStrips << "\t "
           << " packed integer " << std::hex << theBadStripRange  << std::dec
           << std::endl;            
      
      theSiStripVector.push_back(theBadStripRange);
    }      
     
    SiStripBadStrip::Range range(theSiStripVector.begin(),theSiStripVector.end());
    if ( ! obj->put(BadModule_,range) )
      edm::LogError("SiStripBadFiberBuilder")<<"[SiStripBadFiberBuilder::analyze] detid already exists"<<std::endl;
  }
  if (printdebug_)
    edm::LogInfo("SiStripBadFiberBuilder") << ss.str();
}