EcalTPGSpikeThresholdfromFile.cc

Go to the documentation of this file.

#include "CondTools/Ecal/interface/EcalTPGSpikeThresholdfromFile.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::EcalTPGSpikeThresholdfromFile::EcalTPGSpikeThresholdfromFile(const edm::ParameterSet & ps)
  :    m_name(ps.getUntrackedParameter<std::string>("name","EcalTPGSpikeThresholdfromFile")) {

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

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

void popcon::EcalTPGSpikeThresholdfromFile::getNewObjects() {
  std::cout << "------- Ecal -> getNewObjects\n";
    edm::LogInfo("EcalTPGSpikeThresholdfromFile") << "Started GetNewObjects!!!";
  
  Ref payload= lastPayload();
  
  // here popcon tells us which is the last since of the last object in the offline DB

  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 < 365223; i++) getline (fLin, line);
  char tow[8];
  unsigned int towerId, Threshold;
  EcalTPGSpike *lut=new EcalTPGSpike();
  for (int istrip = 0; istrip < 2448; istrip++) {
    getline (fLin, line);
    sscanf(line.c_str(), "%s %u", tow, &towerId);
    if(istrip < 10) std::cout << " string " << tow << " Id " << towerId;
    getline (fLin, line);    // LUTGroupId
    getline (fLin, line);    // FgGroupId
    getline (fLin, line);    // spike_killing_threshold
    if(istrip < 10) std::cout << " line " << line;
    sscanf(line.c_str(), "%u", &Threshold);
    if(istrip < 10) std::cout  << " Threshold " << Threshold << std::endl;

    lut->setValue(towerId, Threshold);
  }   // end loop over EB towers
  fLin.close();

  m_to_transfer.push_back(std::make_pair(lut, fileIOV));

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