L1CaloInputScalesGenerator.cc

Go to the documentation of this file.

#include "L1TriggerConfig/L1ScalesProducers/interface/L1CaloInputScalesGenerator.h"
 
// system include files
#include <memory>
#include <iostream>
using std::cout;
using std::endl;
#include <iomanip>
using std::setprecision;
#include <fstream>
using std::ofstream;
 
// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
 
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
 
#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
#include "CalibFormats/CaloTPG/interface/CaloTPGTranscoder.h"
#include "CalibFormats/CaloTPG/interface/CaloTPGRecord.h"
#include "CalibCalorimetry/EcalTPGTools/interface/EcalTPGScale.h"
#include "DataFormats/EcalDetId/interface/EcalTrigTowerDetId.h"
 
//
// constants, enums and typedefs
//
 
//
// static data member definitions
//
 
//
// constructors and destructor
//
L1CaloInputScalesGenerator::L1CaloInputScalesGenerator(const edm::ParameterSet& iConfig)
 
{
  //now do what ever initialization is needed
}
 
L1CaloInputScalesGenerator::~L1CaloInputScalesGenerator() {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
}
 
//
// member functions
//
 
// ------------ method called to for each event  ------------
void L1CaloInputScalesGenerator::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
  using namespace edm;
 
  ESHandle<CaloTPGTranscoder> caloTPGTranscoder;
  iSetup.get<CaloTPGRecord>().get(caloTPGTranscoder);
 
  EcalTPGScale* ecalTPGScale = new EcalTPGScale();
  ecalTPGScale->setEventSetup(iSetup);
 
  double output;
  ofstream scalesFile("L1CaloInputScales_cfi.py");
 
  // Write the ecal scales, positive eta
 
  scalesFile << "import FWCore.ParameterSet.Config as cms\n" << endl;
 
  scalesFile << "L1CaloInputScalesProducer =cms.ESProducer(\"L1CaloInputScalesProducer\"," << endl;
  scalesFile << "L1EcalEtThresholdsPositiveEta = cms.vdouble(" << endl;
 
  //Python does not support arrays over 255 entries so we neeed ton accomodate it by creating new array after 255 entries
  int nEntries = 0;
 
  // loop over ietas, barrel
  for (unsigned short absIeta = 1; absIeta <= 28; absIeta++) {
    EcalSubdetector subdet = (absIeta <= 17) ? EcalBarrel : EcalEndcap;
    // 8 bits of input energy
    for (unsigned short input = 0; input <= 0xFF; input++) {
      output = ecalTPGScale->getTPGInGeV((unsigned int)input, EcalTrigTowerDetId(1, subdet, absIeta, 1));
      scalesFile << setprecision(8) << output;
      nEntries++;
 
      if (absIeta == 28 && input == 0xFF) {
        scalesFile << "),";
      } else if (nEntries > 254) {
        scalesFile << ")+cms.vdouble(";
        nEntries = 0;
      } else {
        scalesFile << ", ";
      }
    }
    scalesFile << endl;
  }
 
  // Write the ecal scales, negative eta
 
  scalesFile << endl << "\tL1EcalEtThresholdsNegativeEta = cms.vdouble(" << endl;
 
  nEntries = 0;
  // loop over ietas, barrel first
  for (unsigned short absIeta = 1; absIeta <= 28; absIeta++) {
    EcalSubdetector subdet = (absIeta <= 17) ? EcalBarrel : EcalEndcap;
    // 8 bits of input energy
    for (unsigned short input = 0; input <= 0xFF; input++) {
      // negative eta
      output = ecalTPGScale->getTPGInGeV((unsigned int)input, EcalTrigTowerDetId(-1, subdet, absIeta, 2));
      scalesFile << setprecision(8) << output;
      nEntries++;
 
      if (absIeta == 28 && input == 0xFF) {
        scalesFile << "),";
      } else if (nEntries > 254) {
        scalesFile << ")+cms.vdouble(";
        nEntries = 0;
      } else {
        scalesFile << ", ";
      }
    }
    scalesFile << endl;
  }
 
  // Write the hcal scales (Positive Eta)
 
  scalesFile << endl << "\tL1HcalEtThresholdsPositiveEta = cms.vdouble(" << endl;
 
  // loop over ietas
 
  nEntries = 0;
  for (unsigned short absIeta = 1; absIeta <= 32; absIeta++) {
    for (unsigned short input = 0; input <= 0xFF; input++) {
      output = caloTPGTranscoder->hcaletValue(absIeta, input);
      scalesFile << setprecision(8) << output;
      nEntries++;
 
      if (absIeta == 32 && input == 0xFF) {
        scalesFile << "),";
      } else if (nEntries > 254) {
        scalesFile << ")+cms.vdouble(";
        nEntries = 0;
      } else {
        scalesFile << ", ";
      }
    }
    scalesFile << endl;
  }
 
  // Write the hcal scales (Negative Eta)
 
  scalesFile << endl << "\tL1HcalEtThresholdsNegativeEta = cms.vdouble(" << endl;
 
  nEntries = 0;
  // loop over ietas
  for (unsigned short absIeta = 1; absIeta <= 32; absIeta++) {
    for (unsigned short input = 0; input <= 0xFF; input++) {
      output = caloTPGTranscoder->hcaletValue(-absIeta, input);
      scalesFile << setprecision(8) << output;
      nEntries++;
 
      if (absIeta == 32 && input == 0xFF) {
        scalesFile << ")";
      } else if (nEntries > 254) {
        scalesFile << ")+cms.vdouble(";
        nEntries = 0;
      } else {
        scalesFile << ", ";
      }
    }
    scalesFile << endl;
  }
 
  scalesFile << ")" << endl;
 
  scalesFile.close();
}
 
// ------------ method called once each job just before starting event loop  ------------
void L1CaloInputScalesGenerator::beginJob() {}
 
// ------------ method called once each job just after ending the event loop  ------------
void L1CaloInputScalesGenerator::endJob() {}