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/GeomDet.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"));
}

std::unique_ptr<SiStripBadStrip> SiStripBadStripFromASCIIFile::getNewObject() {
  auto SiStripBadStrip_ = std::make_unique<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_;
}