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() {}