CSCGainsDBConditions.h

Go to the documentation of this file.

#ifndef _CSCGAINSDBCONDITIONS_H
#define _CSCGAINSDBCONDITIONS_H

#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/ESProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/EventSetupRecordIntervalFinder.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/SourceFactory.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include <cmath>
#include <memory>

#include "CondFormats/CSCObjects/interface/CSCDBGains.h"
#include "CondFormats/DataRecord/interface/CSCDBGainsRcd.h"
#include <DataFormats/MuonDetId/interface/CSCDetId.h>

class CSCGainsDBConditions : public edm::ESProducer, public edm::EventSetupRecordIntervalFinder {
public:
  CSCGainsDBConditions(const edm::ParameterSet &);
  ~CSCGainsDBConditions() override;

  inline static CSCDBGains *prefillDBGains();

  typedef std::unique_ptr<CSCDBGains> ReturnType;

  ReturnType produceDBGains(const CSCDBGainsRcd &);

private:
  // ----------member data ---------------------------
  void setIntervalFor(const edm::eventsetup::EventSetupRecordKey &,
                      const edm::IOVSyncValue &,
                      edm::ValidityInterval &) override;
  CSCDBGains *cndbGains;
};

#include <fstream>
#include <iostream>
#include <vector>

// to workaround plugin library
inline CSCDBGains *CSCGainsDBConditions::prefillDBGains() {
  // const int MAX_SIZE = 273024;//for ME1a unganged
  const int MAX_SIZE = 252288;
  const int FACTOR = 1000;
  const int MAX_SHORT = 32767;
  CSCDBGains *cndbgains = new CSCDBGains();

  int db_index;
  float db_gainslope;  // db_intercpt, db_chisq;
  std::vector<int> db_index_id;
  std::vector<float> db_slope;
  std::vector<float> db_intercept;
  std::vector<float> db_chi2;
  int new_index;
  float new_gainslope, new_intercpt, new_chisq;
  std::vector<int> new_cham_id;
  std::vector<int> new_index_id;
  std::vector<int> new_strips;
  std::vector<float> new_slope;
  std::vector<float> new_intercept;
  std::vector<float> new_chi2;

  int counter;

  std::ifstream dbdata;
  dbdata.open("old_dbgains.dat", std::ios::in);
  if (!dbdata) {
    std::cerr << "Error: old_dbgains.dat -> no such file!" << std::endl;
    exit(1);
  }

  while (!dbdata.eof()) {
    dbdata >> db_index >> db_gainslope;
    db_index_id.push_back(db_index);
    db_slope.push_back(db_gainslope);
    // db_intercept.push_back(db_intercpt);
    // db_chi2.push_back(db_chisq);
  }
  dbdata.close();

  std::ifstream newdata;
  newdata.open("gains.dat", std::ios::in);
  if (!newdata) {
    std::cerr << "Error: gains.dat -> no such file!" << std::endl;
    exit(1);
  }

  while (!newdata.eof()) {
    newdata >> new_index >> new_gainslope >> new_intercpt >> new_chisq;
    new_index_id.push_back(new_index);
    new_slope.push_back(new_gainslope);
    new_intercept.push_back(new_intercpt);
    new_chi2.push_back(new_chisq);
  }
  newdata.close();

  CSCDBGains::GainContainer &itemvector = cndbgains->gains;
  itemvector.resize(MAX_SIZE);
  cndbgains->factor_gain = (short int)(FACTOR);
  std::cout << " myfactor " << cndbgains->factor_gain << std::endl;

  for (int i = 0; i < MAX_SIZE; ++i) {
    itemvector[i].gain_slope = int(db_slope[i] * FACTOR + 0.5);
  }

  for (int i = 0; i < MAX_SIZE; ++i) {
    counter = db_index_id[i];
    itemvector[i] = itemvector[counter];
    itemvector[i].gain_slope = int(db_slope[i]);

    for (unsigned int k = 0; k < new_index_id.size() - 1; k++) {
      if (counter == new_index_id[k]) {
        if ((short int)(fabs(new_slope[k] * FACTOR + 0.5)) < MAX_SHORT)
          itemvector[counter].gain_slope = int(new_slope[k] * FACTOR + 0.5);
        itemvector[i] = itemvector[counter];
      }
    }
    if (counter > 223968) {
      itemvector[counter].gain_slope = int(db_slope[i]);
      itemvector[i] = itemvector[counter];
    }
  }
  return cndbgains;
}

#endif