SiStripApvGainReader.cc

Go to the documentation of this file.

#include "CondFormats/SiStripObjects/interface/SiStripApvGain.h"
#include "CondFormats/DataRecord/interface/SiStripApvGainRcd.h"
#include "CalibFormats/SiStripObjects/interface/SiStripGain.h"
#include "CalibTracker/Records/interface/SiStripGainRcd.h"

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

#include <iostream>
#include <cstdio>
#include <sys/time.h>

using namespace cms;

SiStripApvGainReader::SiStripApvGainReader(const edm::ParameterSet& iConfig)
    : printdebug_(iConfig.getUntrackedParameter<bool>("printDebug", true)),
      formatedOutput_(iConfig.getUntrackedParameter<std::string>("outputFile", "")),
      gainType_(iConfig.getUntrackedParameter<uint32_t>("gainType", 1)),
      tree_(nullptr) {
  if (fs_.isAvailable()) {
    tree_ = fs_->make<TTree>("Gains", "Gains");

    tree_->Branch("Index", &id_, "Index/I");
    tree_->Branch("DetId", &detId_, "DetId/I");
    tree_->Branch("APVId", &apvId_, "APVId/I");
    tree_->Branch("Gain", &gain_, "Gain/D");
  }
}

SiStripApvGainReader::~SiStripApvGainReader() {}

void SiStripApvGainReader::analyze(const edm::Event& e, const edm::EventSetup& iSetup) {
  edm::ESHandle<SiStripGain> SiStripApvGain_;
  iSetup.get<SiStripGainRcd>().get(SiStripApvGain_);
  edm::LogInfo("SiStripApvGainReader") << "[SiStripApvGainReader::analyze] End Reading SiStripApvGain" << std::endl;
  std::vector<uint32_t> detid;
  SiStripApvGain_->getDetIds(detid);
  edm::LogInfo("Number of detids ") << detid.size() << std::endl;

  FILE* pFile = nullptr;
  if (!formatedOutput_.empty())
    pFile = fopen(formatedOutput_.c_str(), "w");
  for (size_t id = 0; id < detid.size(); id++) {
    SiStripApvGain::Range range = SiStripApvGain_->getRange(detid[id], gainType_);
    if (printdebug_) {
      int apv = 0;
      for (int it = 0; it < range.second - range.first; it++) {
        edm::LogInfo("SiStripApvGainReader") << "detid " << detid[id] << " \t " << apv++ << " \t "
                                             << SiStripApvGain_->getApvGain(it, range) << std::endl;
        id_++;

        if (tree_) {
          detId_ = detid[id];
          apvId_ = apv;
          gain_ = SiStripApvGain_->getApvGain(it, range);
          tree_->Fill();
        }
      }
    }

    if (pFile) {
      fprintf(pFile, "%i ", detid[id]);
      for (int it = 0; it < range.second - range.first; it++) {
        fprintf(pFile, "%f ", SiStripApvGain_->getApvGain(it, range));
      }
      fprintf(pFile, "\n");
    }
  }

  if (pFile)
    fclose(pFile);
}