/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/CalibTracker/SiStripQuality/plugins/SiStripBadStripFromASCIIFile.cc

Go to the documentation of this file.

// system include files
#include <memory>
#include <cstdio>
#include <string>


// user include files


#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h" 
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "Geometry/CommonDetUnit/interface/GeomDetUnit.h"
#include "Geometry/CommonDetUnit/interface/GeomDetType.h"
#include "Geometry/CommonTopologies/interface/StripTopology.h"
#include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetUnit.h"
#include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetType.h"
#include "DataFormats/FEDRawData/interface/FEDNumbering.h"


#include "CLHEP/Random/RandFlat.h"
#include "CLHEP/Random/RandGauss.h"



#include "CalibTracker/SiStripQuality/plugins/SiStripBadStripFromASCIIFile.h"



using namespace std;

SiStripBadStripFromASCIIFile::SiStripBadStripFromASCIIFile( const edm::ParameterSet& iConfig ):
  ConditionDBWriter<SiStripBadStrip>( iConfig),
  printdebug_(iConfig.getUntrackedParameter<bool>("printDebug",false))
{
  fp_ = iConfig.getUntrackedParameter<edm::FileInPath>("file",edm::FileInPath("CalibTracker/SiStripQuality/data/DefectsFromConstructionDB.dat"));
 }


SiStripBadStrip *SiStripBadStripFromASCIIFile::getNewObject() {

  SiStripBadStrip* SiStripBadStrip_ = new SiStripBadStrip();

   // open file and fill DB
  ifstream infile((fp_.fullPath()).c_str());
  if(!infile){std::cout << "Problem while trying to open File: " << (fp_.fullPath()).c_str() << std::endl;}
  

  //variables needed for reading file and filling of SiStripBadStripObject
  uint32_t detid;
  short flag;
  short channel;
 
  bool firstrun=true;
  short tempchannel=0;
  int count =0;
  std::vector<unsigned int> theSiStripVector;
  short tempflag=0;
  uint32_t tempdetid=0;


   while(!infile.eof()){

    // get data from file: 
    //infile >> detid >> channel >> flag;

    //if no flag is available, use the following:
    infile >> detid >> channel;
    flag = 1;

    unsigned int theBadStripRange=0;

    // first loop ?
    if(firstrun) {
      tempdetid=detid;
      tempchannel=channel;
      tempflag=flag;
      count =0;
      firstrun=false;
    }


    if(detid==tempdetid){
      if (channel!=tempchannel+count || flag != tempflag){
                                          // 1.badstrip, nconsectrips, flag
        theBadStripRange = SiStripBadStrip_->encode(tempchannel-1, count, tempflag); // In the quality object, strips are counted from 0 to 767!!! Therefore "tempchannel-1"!
                                                                                     // In the input txt-file, they have to be from 1 to 768 instead!!!
          edm::LogInfo("SiStripBadStripFromASCIIFile")<< "detid " << tempdetid << " \t"
                                                           << " firstBadStrip " << tempchannel << "\t "
                                                           << " NconsecutiveBadStrips " << count  << "\t "
                                                           <<"flag " << tempflag << "\t"
                                                           << " packed integer " << std::hex << theBadStripRange  << std::dec
                                                           << std::endl;            

          theSiStripVector.push_back(theBadStripRange);

          if (infile.eof()){ // Don't forget to save the last strip before eof!!!
            SiStripBadStrip::Range range(theSiStripVector.begin(),theSiStripVector.end());
            if ( ! SiStripBadStrip_->put(tempdetid,range) )
              edm::LogError("SiStripBadStripFromASCIIFile")<<"[SiStripBadStripFromASCIIFile::GetNewObject] detid already exists"<<std::endl;
            theSiStripVector.clear();
          }

          count=1;
          tempchannel=channel;
          tempflag=flag;
                
      }else{count++;}
    }
    
    if(detid!=tempdetid){
                                                    // 1.badstrip, nconsectrips, flag
          theBadStripRange = SiStripBadStrip_->encode(tempchannel-1, count, tempflag); // In the quality object, strips are counted from 0 to 767!!! Therefore "tempchannel-1"!
                                                                                       // In the input txt-file, they have to be from 1 to 768 instead!!!
          edm::LogInfo("SiStripBadStripFromASCIIFile")<< "detid " << tempdetid << " \t"
                                                           << " firstBadStrip " << tempchannel << "\t "
                                                           << " NconsecutiveBadStrips " << count  << "\t "
                                                           <<"flag " << tempflag << "\t"
                                                           << " packed integer " << std::hex << theBadStripRange  << std::dec
                                                           << std::endl;            
          
          theSiStripVector.push_back(theBadStripRange);
          
          // populate db  object
          SiStripBadStrip::Range range(theSiStripVector.begin(),theSiStripVector.end());
    if ( ! SiStripBadStrip_->put(tempdetid,range) )
      edm::LogError("SiStripBadStripFromASCIIFile")<<"[SiStripBadStripFromASCIIFile::GetNewObject] detid already exists"<<std::endl;
    theSiStripVector.clear();
    
    count=1;
    tempdetid=detid;
    tempchannel=channel;
    tempflag=flag;
    }
   }
 
  return SiStripBadStrip_;
}