CSCNoiseMatrixDBConditions.h

Go to the documentation of this file.

#ifndef _CSCNOISEMATRIXDBCONDITIONS_H
#define _CSCNOISEMATRIXDBCONDITIONS_H

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

#include <DataFormats/MuonDetId/interface/CSCDetId.h>
#include "CondFormats/CSCObjects/interface/CSCDBNoiseMatrix.h"
#include "CondFormats/DataRecord/interface/CSCDBNoiseMatrixRcd.h"

class CSCNoiseMatrixDBConditions: public edm::ESProducer, public edm::EventSetupRecordIntervalFinder  {
 public:
  CSCNoiseMatrixDBConditions(const edm::ParameterSet&);
  ~CSCNoiseMatrixDBConditions() override;
  
  inline static CSCDBNoiseMatrix * prefillDBNoiseMatrix();

  typedef const  CSCDBNoiseMatrix * ReturnType;
  
  ReturnType produceDBNoiseMatrix(const CSCDBNoiseMatrixRcd&);
  
 private:
  // ----------member data ---------------------------
  void setIntervalFor(const edm::eventsetup::EventSetupRecordKey &, const edm::IOVSyncValue&, edm::ValidityInterval & ) override;
  CSCDBNoiseMatrix *cndbMatrix ;

};


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

// to workaround plugin library
inline CSCDBNoiseMatrix *  CSCNoiseMatrixDBConditions::prefillDBNoiseMatrix()
{
  //const int MAX_SIZE = 273024; //for ME1a unganged
  const int MAX_SIZE = 252288;
  const int FACTOR=1000;
  const int MAX_SHORT= 32767;
 
  int new_index, db_index;
  float db_elm33,db_elm34, db_elm44, db_elm35, db_elm45, db_elm55;
  float db_elm46, db_elm56, db_elm66, db_elm57, db_elm67, db_elm77;
  std::vector<int> db_index_id;
  std::vector<float> db_elem33;
  std::vector<float> db_elem34;
  std::vector<float> db_elem44;
  std::vector<float> db_elem45;
  std::vector<float> db_elem35;
  std::vector<float> db_elem55;
  std::vector<float> db_elem46;
  std::vector<float> db_elem56;
  std::vector<float> db_elem66;
  std::vector<float> db_elem57;
  std::vector<float> db_elem67;
  std::vector<float> db_elem77;


  float new_elm33,new_elm34, new_elm44, new_elm35, new_elm45, new_elm55;
  float  new_elm46, new_elm56, new_elm66, new_elm57, new_elm67, new_elm77;
  std::vector<int> new_cham_id;
  std::vector<int> new_index_id;
  std::vector<float> new_elem33;
  std::vector<float> new_elem34;
  std::vector<float> new_elem44;
  std::vector<float> new_elem45;
  std::vector<float> new_elem35;
  std::vector<float> new_elem55;
  std::vector<float> new_elem46;
  std::vector<float> new_elem56;
  std::vector<float> new_elem66;
  std::vector<float> new_elem57;
  std::vector<float> new_elem67;
  std::vector<float> new_elem77;

  CSCDBNoiseMatrix * cndbmatrix = new CSCDBNoiseMatrix();

 int counter;
  int db_nrlines=0;
  int new_nrlines=0;
    
  std::ifstream dbdata; 
  dbdata.open("old_dbmatrix.dat",std::ios::in); 
  if(!dbdata) {
    std::cerr <<"Error: old_dbmatrix.dat -> no such file!"<< std::endl;
    exit(1);
  }
  
  while (!dbdata.eof() ) { 
    dbdata >> db_index >> db_elm33 >> db_elm34 >> db_elm44 >> db_elm35 >> db_elm45 >> db_elm55 >> db_elm46 >> db_elm56 >> db_elm66 >> db_elm57 >> db_elm67 >> db_elm77 ; 
    db_index_id.push_back(db_index);
    db_elem33.push_back(db_elm33);
    db_elem34.push_back(db_elm34);
    db_elem35.push_back(db_elm35);
    db_elem44.push_back(db_elm44);
    db_elem45.push_back(db_elm45);
    db_elem46.push_back(db_elm46);
    db_elem55.push_back(db_elm55);
    db_elem56.push_back(db_elm56);
    db_elem57.push_back(db_elm57);
    db_elem66.push_back(db_elm66);
    db_elem67.push_back(db_elm67);
    db_elem77.push_back(db_elm77);
    db_nrlines++;
  }
  dbdata.close();

  std::ifstream newdata;
  newdata.open("matrix.dat",std::ios::in); 
  if(!newdata) {
    std::cerr <<"Error: matrix.dat -> no such file!"<< std::endl;
    exit(1);
  }
  
  while (!newdata.eof() ) { 
    newdata >> new_index>> new_elm33 >> new_elm34 >> new_elm44 >> new_elm35 >> new_elm45 >> new_elm55 >> new_elm46 >> new_elm56 >> new_elm66 >> new_elm57 >> new_elm67 >> new_elm77 ; 
    //new_cham_id.push_back(new_chamber_id);
    new_index_id.push_back(new_index);
    new_elem33.push_back(new_elm33);
    new_elem34.push_back(new_elm34);
    new_elem35.push_back(new_elm35);
    new_elem44.push_back(new_elm44);
    new_elem45.push_back(new_elm45);
    new_elem46.push_back(new_elm46);
    new_elem55.push_back(new_elm55);
    new_elem56.push_back(new_elm56);
    new_elem57.push_back(new_elm57);
    new_elem66.push_back(new_elm66);
    new_elem67.push_back(new_elm67);
    new_elem77.push_back(new_elm77); 
    new_nrlines++;
  }
  newdata.close();
  
  CSCDBNoiseMatrix::NoiseMatrixContainer & itemvector = cndbmatrix->matrix;
  itemvector.resize(MAX_SIZE);
  cndbmatrix->factor_noise= int (FACTOR);
  
  for(int i=0; i<MAX_SIZE;++i){
    itemvector[i].elem33 = (short int) (db_elem33[i]*FACTOR+0.5);
    itemvector[i].elem34 = (short int) (db_elem34[i]*FACTOR+0.5); 
    itemvector[i].elem35 = (short int) (db_elem35[i]*FACTOR+0.5);
    itemvector[i].elem44 = (short int) (db_elem44[i]*FACTOR+0.5);
    itemvector[i].elem45 = (short int) (db_elem45[i]*FACTOR+0.5);
    itemvector[i].elem46 = (short int) (db_elem46[i]*FACTOR+0.5);
    itemvector[i].elem55 = (short int) (db_elem55[i]*FACTOR+0.5);
    itemvector[i].elem56 = (short int) (db_elem56[i]*FACTOR+0.5);
    itemvector[i].elem57 = (short int) (db_elem57[i]*FACTOR+0.5);
    itemvector[i].elem66 = (short int) (db_elem66[i]*FACTOR+0.5);
    itemvector[i].elem67 = (short int) (db_elem67[i]*FACTOR+0.5);
    itemvector[i].elem77 = (short int) (db_elem77[i]*FACTOR+0.5);
  }

  
  for(int i=0; i<MAX_SIZE;++i){
    counter=db_index_id[i];  
    itemvector[i] = itemvector[counter];
    itemvector[i].elem33 = int (db_elem33[i]);
    itemvector[i].elem34 = int (db_elem34[i]);
    itemvector[i].elem35 = int (db_elem35[i]);
    itemvector[i].elem44 = int (db_elem44[i]);
    itemvector[i].elem45 = int (db_elem45[i]);
    itemvector[i].elem46 = int (db_elem46[i]);
    itemvector[i].elem55 = int (db_elem55[i]);
    itemvector[i].elem56 = int (db_elem56[i]);
    itemvector[i].elem57 = int (db_elem57[i]);
    itemvector[i].elem66 = int (db_elem66[i]);
    itemvector[i].elem67 = int (db_elem67[i]);
    itemvector[i].elem77 = int (db_elem77[i]);
    
    for (unsigned int k=0;k<new_index_id.size()-1;k++){
      if(counter==new_index_id[k]){
    if((short int) (fabs(new_elem33[k]*FACTOR+0.5))<MAX_SHORT) itemvector[counter].elem33 = int (new_elem33[k]*FACTOR+0.5);
    if((short int) (fabs(new_elem34[k]*FACTOR+0.5))<MAX_SHORT) itemvector[counter].elem34 = int (new_elem34[k]*FACTOR+0.5); 
    if((short int) (fabs(new_elem35[k]*FACTOR+0.5))<MAX_SHORT) itemvector[counter].elem35 = int (new_elem35[k]*FACTOR+0.5);
    if((short int) (fabs(new_elem44[k]*FACTOR+0.5))<MAX_SHORT) itemvector[counter].elem44 = int (new_elem44[k]*FACTOR+0.5);
    if((short int) (fabs(new_elem45[k]*FACTOR+0.5))<MAX_SHORT) itemvector[counter].elem45 = int (new_elem45[k]*FACTOR+0.5);
    if((short int) (fabs(new_elem46[k]*FACTOR+0.5))<MAX_SHORT) itemvector[counter].elem46 = int (new_elem46[k]*FACTOR+0.5);
    if((short int) (fabs(new_elem55[k]*FACTOR+0.5))<MAX_SHORT) itemvector[counter].elem55 = int (new_elem55[k]*FACTOR+0.5);
    if((short int) (fabs(new_elem56[k]*FACTOR+0.5))<MAX_SHORT) itemvector[counter].elem56 = int (new_elem56[k]*FACTOR+0.5);
    if((short int) (fabs(new_elem57[k]*FACTOR+0.5))<MAX_SHORT) itemvector[counter].elem57 = int (new_elem57[k]*FACTOR+0.5);
    if((short int) (fabs(new_elem66[k]*FACTOR+0.5))<MAX_SHORT) itemvector[counter].elem66 = int (new_elem66[k]*FACTOR+0.5);
    if((short int) (fabs(new_elem67[k]*FACTOR+0.5))<MAX_SHORT) itemvector[counter].elem67 = int (new_elem67[k]*FACTOR+0.5);
    if((short int) (fabs(new_elem77[k]*FACTOR+0.5))<MAX_SHORT) itemvector[counter].elem77 = int (new_elem77[k]*FACTOR+0.5);
    itemvector[i] = itemvector[counter];
      }  
    }

    if (counter>223968){
      itemvector[counter].elem33 = int (db_elem33[i]);
      itemvector[counter].elem34 = int (db_elem34[i]);
      itemvector[counter].elem35 = int (db_elem35[i]);
      itemvector[counter].elem44 = int (db_elem44[i]);
      itemvector[counter].elem45 = int (db_elem45[i]);
      itemvector[counter].elem46 = int (db_elem46[i]);
      itemvector[counter].elem55 = int (db_elem55[i]);
      itemvector[counter].elem56 = int (db_elem56[i]);
      itemvector[counter].elem57 = int (db_elem57[i]);
      itemvector[counter].elem66 = int (db_elem66[i]);
      itemvector[counter].elem67 = int (db_elem67[i]);
      itemvector[counter].elem77 = int (db_elem77[i]);
      itemvector[i] = itemvector[counter];
    }
  }
  
  return cndbmatrix;
}

#endif