SiStripDetVOffFakeBuilder.cc

Go to the documentation of this file.

// system include files
#include <memory>

// user include files

#include "CondCore/DBOutputService/interface/PoolDBOutputService.h"


#include "CondFormats/SiStripObjects/interface/SiStripDetVOff.h"

#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h" 
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "Geometry/CommonDetUnit/interface/GeomDetUnit.h"
#include "Geometry/CommonTopologies/interface/StripTopology.h"
#include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetUnit.h"


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

using namespace std;
using namespace cms;

SiStripDetVOffFakeBuilder::SiStripDetVOffFakeBuilder( const edm::ParameterSet& iConfig ):
  printdebug_(iConfig.getUntrackedParameter<bool>("printDebug",false)){}

SiStripDetVOffFakeBuilder::~SiStripDetVOffFakeBuilder(){}


void SiStripDetVOffFakeBuilder::initialize( const edm::EventSetup& iSetup ) {

  edm::ESHandle<TrackerGeometry> pDD;
  iSetup.get<TrackerDigiGeometryRecord>().get( pDD );
  edm::LogInfo("SiStripDetVOffFakeBuilder") <<" There are "<<pDD->detUnits().size() <<" detectors"<<std::endl;
  
  for(TrackerGeometry::DetUnitContainer::const_iterator it = pDD->detUnits().begin(); it != pDD->detUnits().end(); it++){
  
    if( dynamic_cast<StripGeomDetUnit const*>((*it))!=0){
      uint32_t detid=((*it)->geographicalId()).rawId();            
      const StripTopology& p = dynamic_cast<StripGeomDetUnit const*>((*it))->specificTopology();
      unsigned short Nstrips = p.nstrips();
      if(Nstrips<1 || Nstrips>768 ) {
    edm::LogError("SiStripDetVOffFakeBuilder")<<" Problem with Number of strips in detector.. "<< p.nstrips() <<" Exiting program"<<endl;
    exit(1);
      }
      detids.push_back(detid);
      if (printdebug_)
    edm::LogInfo("SiStripDetVOffFakeBuilder")<< "detid " << detid;
    }
  }
}

void SiStripDetVOffFakeBuilder::analyze(const edm::Event& evt, const edm::EventSetup& iSetup)
{
  initialize(iSetup);

  unsigned int run=evt.id().run();

  edm::LogInfo("SiStripDetVOffFakeBuilder") << "... creating dummy SiStripDetVOff Data for Run " << run << "\n " << std::endl;



  SiStripDetVOff* SiStripDetVOff_ = new SiStripDetVOff();

   // std::vector<uint32_t> TheDetIdHVVector;

    for(std::vector<uint32_t>::const_iterator it = detids.begin(); it != detids.end(); it++){

    //Generate HV and LV for each channel, if at least one of the two is off fill the value
    int hv=rand() % 20;
    int lv=rand() % 20;
    if( hv<=2 ) {
      edm::LogInfo("SiStripDetVOffFakeBuilder") << "detid: " <<  *it << " HV\t OFF" << std::endl;
      SiStripDetVOff_->put( *it, 1, -1 );
      // TheDetIdHVVector.push_back(*it);
    }
    if( lv<=2 ) {
      edm::LogInfo("SiStripDetVOffFakeBuilder") << "detid: " <<  *it << " LV\t OFF" << std::endl;
      SiStripDetVOff_->put( *it, -1, 1 );
      // TheDetIdHVVector.push_back(*it);
    }
    if( lv<=2 || hv<=2 ) edm::LogInfo("SiStripDetVOffFakeBuilder") << "detid: " <<  *it << " V\t OFF" << std::endl;
  }

  // SiStripDetVOff_->put(TheDetIdHVVector);



  //End now write DetVOff data in DB
  edm::Service<cond::service::PoolDBOutputService> mydbservice;
  
  if( mydbservice.isAvailable() ){
    try{
      if( mydbservice->isNewTagRequest("SiStripDetVOffRcd") ){
    mydbservice->createNewIOV<SiStripDetVOff>(SiStripDetVOff_,mydbservice->beginOfTime(),mydbservice->endOfTime(),"SiStripDetVOffRcd");      
      } else {
    mydbservice->appendSinceTime<SiStripDetVOff>(SiStripDetVOff_,mydbservice->currentTime(),"SiStripDetVOffRcd");      
      }
    }catch(const cond::Exception& er){
      edm::LogError("SiStripDetVOffFakeBuilder")<<er.what()<<std::endl;
    }catch(const std::exception& er){
      edm::LogError("SiStripDetVOffFakeBuilder")<<"caught std::exception "<<er.what()<<std::endl;
    }
  }else{
    edm::LogError("SiStripDetVOffFakeBuilder")<<"Service is unavailable"<<std::endl;
  }
}