EcalTPGFineGrainStripfromFile.cc

Go to the documentation of this file.

#include "CondTools/Ecal/interface/EcalTPGFineGrainStripfromFile.h"

#include "FWCore/Framework/interface/ESHandle.h"

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

#include "Geometry/EcalMapping/interface/EcalElectronicsMapping.h"
#include "Geometry/EcalMapping/interface/EcalMappingRcd.h"

#include<iostream>
#include<fstream>
#include <sstream>

popcon::EcalTPGFineGrainStripfromFile::EcalTPGFineGrainStripfromFile(const edm::ParameterSet & ps)
  :    m_name(ps.getUntrackedParameter<std::string>("name","EcalTPGFineGrainStripfromFile")) {
  fname = ps.getParameter<std::string>("FileName");

  std::cout << "EcalTPGFineGrainStripfromFile constructor\n" << std::endl;
}

popcon::EcalTPGFineGrainStripfromFile::~EcalTPGFineGrainStripfromFile(){
  // do nothing
}

void popcon::EcalTPGFineGrainStripfromFile::getNewObjects() {
  std::cout << "------- Ecal -> getNewObjects\n";
    edm::LogInfo("EcalTPGFineGrainStripfromFile") << "Started GetNewObjects!!!";
  
  EcalTPGFineGrainStripEE * fgrStripEE = new EcalTPGFineGrainStripEE;
  int fileIOV;
  std::cout << "LinPed which input IOV do you want " << std::endl;
  std::cin >> fileIOV; 
  std::ifstream fLin;
  std::ostringstream oss;
  oss << fileIOV;
  //  std::string fname = "/afs/cern.ch/cms/ECAL/triggerTransp/TPG_beamv6_trans_" + oss.str() + "_spikekill.txt";
  fLin.open(fname.c_str());
  if(!fLin.is_open()) {
    std::cout << "ERROR : can't open file '" << fname << std::endl;
    return;
  }
  std::cout << " file " << fname << " opened" << std::endl;
  std::string line;
  for(int i = 0; i < 304517; i++) getline (fLin, line);
  char strip[8];
  unsigned int stripId, Threshold, LUTFgr;
  for (int istrip = 0; istrip < 15176; istrip++) {
    getline (fLin, line);
    sscanf(line.c_str(), "%s %u", strip, &stripId);
    if(istrip < 10 || (istrip > 12239 && istrip < 12250)) std::cout << " string " << strip << " Id " << stripId;
    getline (fLin, line);    // sliding_window
    getline (fLin, line);    // weightGroupId
    getline (fLin, line);    // threshold_sfg lut_sfg
    if(istrip < 10 || (istrip > 12239 && istrip < 12250)) std::cout << " line " << line;
    sscanf(line.c_str(), "%x %x", &Threshold, &LUTFgr);
    if(istrip < 10 || (istrip > 12239 && istrip < 12250)) std::cout  << " Threshold " << Threshold << std::endl;
    EcalTPGFineGrainStripEE::Item item;
    item.threshold = Threshold;
    item.lut = LUTFgr;

    fgrStripEE->setValue(stripId, item);
    if(istrip == 12239) getline (fLin, line);    // 1 empty line between EB and EE
  }   // end loop over EB + EE strips
  fLin.close();

  m_to_transfer.push_back(std::make_pair((EcalTPGFineGrainStripEE*)fgrStripEE, fileIOV));

  std::cout << "Ecal -> end of getNewObjects -----------\n";    
}