CSCTFConfigProducer.cc

Go to the documentation of this file.

#include <L1TriggerConfig/CSCTFConfigProducers/interface/CSCTFConfigProducer.h>
#include <FWCore/MessageLogger/interface/MessageLogger.h>

#include <stdio.h>
#include <errno.h>
#include <iostream>
#include <fstream>

CSCTFConfigProducer::CSCTFConfigProducer(const edm::ParameterSet& pset) {
  const char *name[12] = {"registersSP1", "registersSP2", "registersSP3", "registersSP4",
              "registersSP5", "registersSP6", "registersSP7", "registersSP8",
              "registersSP9", "registersSP10","registersSP11","registersSP12"};
    
  for(int sp=0; sp<12; sp++){
    std::vector<std::string> regs = pset.getParameter<std::vector<std::string> >(name[sp]);
    for(std::vector<std::string>::const_iterator line=regs.begin(); line!=regs.end(); line++)
      registers[sp] += *line + "\n";
  }
    
  alignment  = pset.getParameter< std::vector<double> >("alignment");
  ptLUT_path = pset.getParameter< std::string >        ("ptLUT_path");
  setWhatProduced(this, &CSCTFConfigProducer::produceL1MuCSCTFConfigurationRcd);
  setWhatProduced(this, &CSCTFConfigProducer::produceL1MuCSCTFAlignmentRcd);
  setWhatProduced(this, &CSCTFConfigProducer::produceL1MuCSCPtLutRcd);

}

std::auto_ptr<L1MuCSCTFConfiguration> CSCTFConfigProducer::produceL1MuCSCTFConfigurationRcd(const L1MuCSCTFConfigurationRcd& iRecord){

  edm::LogInfo( "L1-O2O: CSCTFConfigProducer" ) << "Producing "
                        << " L1MuCSCTFConfiguration from PSET";
                        
  std::auto_ptr<L1MuCSCTFConfiguration> config = std::auto_ptr<L1MuCSCTFConfiguration>( new L1MuCSCTFConfiguration(registers) );
  return config;
}

std::auto_ptr<L1MuCSCTFAlignment> CSCTFConfigProducer::produceL1MuCSCTFAlignmentRcd(const L1MuCSCTFAlignmentRcd& iRecord){
  edm::LogInfo( "L1-O2O: CSCTFConfigProducer" ) << "Producing "
                        << " L1MuCSCTFAlignment from PSET";
  

  std::auto_ptr<L1MuCSCTFAlignment> al = std::auto_ptr<L1MuCSCTFAlignment>( new L1MuCSCTFAlignment(alignment) );
  return al;
}

std::auto_ptr<L1MuCSCPtLut> CSCTFConfigProducer::produceL1MuCSCPtLutRcd(const L1MuCSCPtLutRcd& iRecord){
  edm::LogInfo( "L1-O2O: CSCTFConfigProducer" ) << "Producing "
                        << " L1MuCSCPtLut from PSET";
  


  std::auto_ptr<L1MuCSCPtLut> pt_lut = std::auto_ptr<L1MuCSCPtLut>( new L1MuCSCPtLut() );

  if( ptLUT_path.length() ){
    readLUT(ptLUT_path, (unsigned short *)pt_lut->pt_lut, 1<<21); //CSCBitWidths::kPtAddressWidth
  } else {
    throw cms::Exception("Undefined pT LUT")<<"CSCTFConfigProducer is unable to generate LUTs on the fly.\n"
      "Specify full LUT file names or just avoid using CSCTFConfigProducer by uncommenting PTLUT "
      "parameter sets in L1Trigger/CSCTrackFinder configuration."<<std::endl;
  }
  return pt_lut;
}

void CSCTFConfigProducer::readLUT(std::string path, unsigned short* lut, unsigned long length){
  // Reading
  if( path.find(".bin") != std::string::npos ) { // Binary format
    std::ifstream file(path.c_str(), std::ios::binary);
    file.read((char*)lut,length*sizeof(unsigned short));
    if( file.fail() )
      throw cms::Exception("Reading error")<<"CSCTFConfigProducer cannot read "<<length<<" words from "<<path<<" (errno="<<errno<<")"<<std::endl;
    if( (unsigned int)file.gcount() != length*sizeof(unsigned short) )
      throw cms::Exception("Incorrect LUT size")<<"CSCTFConfigProducer read "<<(file.gcount()/sizeof(unsigned short))<<" words from "<<path<<" instead of "<<length<<" (errno="<<errno<<")"<<std::endl;
    file.close();
  } else {
    std::ifstream file(path.c_str());
    if( file.fail() )
      throw cms::Exception("Cannot open file")<<"CSCTFConfigProducer cannot open "<<path<<" (errno="<<errno<<")"<<std::endl;
    unsigned int address=0;
    for(address=0; !file.eof() && address<length; address++)
      file >> lut[address]; // Warning: this may throw non-cms like exception
    if( address!=length ) throw cms::Exception("Incorrect LUT size")<<"CSCTFConfigProducer read "<<address<<" words from "<<path<<" instead of "<<length<<std::endl;
    file.close();
  }
  LogDebug("CSCTFConfigProducer::readLUT")<<" read from "<<path<<" "<<length<<" words"<<std::endl;
}