L1ScalesTester.cc

Go to the documentation of this file.

#include "L1TriggerConfig/L1ScalesProducers/interface/L1ScalesTester.h"

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

#include "CondFormats/L1TObjects/interface/L1CaloEtScale.h"
#include "CondFormats/DataRecord/interface/L1EmEtScaleRcd.h"
#include "CondFormats/DataRecord/interface/L1JetEtScaleRcd.h"
#include "CondFormats/L1TObjects/interface/L1CaloEcalScale.h"
#include "CondFormats/L1TObjects/interface/L1CaloHcalScale.h"
#include "CondFormats/DataRecord/interface/L1CaloEcalScaleRcd.h"
#include "CondFormats/DataRecord/interface/L1CaloHcalScaleRcd.h"

#include <iostream>

using std::cout;
using std::endl;

L1ScalesTester::L1ScalesTester(const edm::ParameterSet& ps) { cout << "Constructing a L1ScalesTester" << endl; }

L1ScalesTester::~L1ScalesTester() {}

void L1ScalesTester::analyze(const edm::Event& e, const edm::EventSetup& es) {
  using namespace edm;

  ESHandle<L1CaloEtScale> emScale;
  es.get<L1EmEtScaleRcd>().get(emScale);

  cout << "L1EmEtScaleRcd :" << endl;
  emScale->print(cout);
  cout << endl;

  ESHandle<L1CaloEcalScale> ecalScale;
  es.get<L1CaloEcalScaleRcd>().get(ecalScale);

  ESHandle<L1CaloHcalScale> hcalScale;
  es.get<L1CaloHcalScaleRcd>().get(hcalScale);

  cout << " L1ColoEcalScale  :" << endl;
  ecalScale->print(cout);
  cout << endl;

  cout << " L1ColoHcalScale  :" << endl;
  hcalScale->print(cout);
  cout << endl;

  ESHandle<L1CaloEtScale> jetScale;
  es.get<L1JetEtScaleRcd>().get(jetScale);

  cout << "L1JetEtScaleRcd :" << endl;
  jetScale->print(cout);
  cout << endl;

  // test EM lin-rank conversion
  cout << "Testing EM linear-to-rank conversion" << endl;
  for (unsigned short i = 0; i < 32; i++) {
    unsigned rank = emScale->rank(i);
    cout << "EM linear : " << i << ", Et : " << i * emScale->linearLsb() << " GeV, rank : " << rank << endl;
  }
  cout << endl;

  // test jet lin-rank conversion
  cout << "Testing jet linear-to-rank conversion" << endl;
  for (unsigned short i = 0; i < 32; i++) {
    unsigned rank = jetScale->rank(i);
    cout << "jet linear : " << i << ", Et : " << i * jetScale->linearLsb() << " GeV, rank : " << rank << endl;
  }
  cout << endl;

  // test EM rank-et conversion
  cout << "Testing EM rank-to-Et conversion" << endl;
  for (unsigned i = 0; i < 32; i++) {
    double et = emScale->et(i);
    cout << "EM rank : " << i << " Et : " << et << " GeV" << endl;
  }
  cout << endl;

  // test jet rank-et conversion
  cout << "Testing jet rank-to-Et conversion" << endl;
  for (unsigned i = 0; i < 32; i++) {
    double et = jetScale->et(i);
    cout << "jet rank : " << i << " Et : " << et << " GeV" << endl;
  }
  cout << endl;
}