L1CaloInputScaleTester.cc

Go to the documentation of this file.

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

// system include files
#include <memory>
#include <iostream>
using std::cout;
using std::endl;
#include <iomanip>
using std::setprecision;

// 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 "CondFormats/L1TObjects/interface/L1CaloEcalScale.h"
#include "CondFormats/DataRecord/interface/L1CaloEcalScaleRcd.h"
#include "CondFormats/L1TObjects/interface/L1CaloHcalScale.h"
#include "CondFormats/DataRecord/interface/L1CaloHcalScaleRcd.h"
#include "DataFormats/EcalDetId/interface/EcalTrigTowerDetId.h"

//
// constants, enums and typedefs
//

//
// static data member definitions
//

//
// constructors and destructor
//
L1CaloInputScaleTester::L1CaloInputScaleTester(const edm::ParameterSet& iConfig)
    : hcalScaleToken_(esConsumes<L1CaloHcalScale, L1CaloHcalScaleRcd>()),
      ecalScaleToken_(esConsumes<L1CaloEcalScale, L1CaloEcalScaleRcd>()),
      transcoderToken_(esConsumes<CaloTPGTranscoder, CaloTPGRecord>()),
      tokens_(consumesCollector()) {
  //now do what ever initialization is needed
}

L1CaloInputScaleTester::~L1CaloInputScaleTester() {
  // do anything here that needs to be done at destruction time
  // (e.g. close files, deallocate resources etc.)
}

//
// member functions
//

// ------------ method called to for each event  ------------
void L1CaloInputScaleTester::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
  using namespace edm;

#ifdef THIS_IS_AN_EVENT_EXAMPLE
  Handle<ExampleData> pIn;
  iEvent.getByLabel("example", pIn);
#endif

#ifdef THIS_IS_AN_EVENTSETUP_EXAMPLE
  ESHandle<SetupData> pSetup;
  iSetup.get<SetupRecord>().get(pSetup);
#endif

  ESHandle<L1CaloEcalScale> caloEcalScale = iSetup.getHandle(ecalScaleToken_);
  ESHandle<L1CaloHcalScale> caloHcalScale = iSetup.getHandle(hcalScaleToken_);
  ESHandle<CaloTPGTranscoder> caloTPGTranscoder = iSetup.getHandle(transcoderToken_);

  EcalTPGScale ecalTPGScale(tokens_, iSetup);

  bool ecalIsConsistent = true;
  bool hcalIsConsistent = true;

  double ecal1;
  double ecal2;
  double hcal1;
  double hcal2;
  double hcal3;
  double hcal4;

  // compare the ecal scales

  // 8 bits of input energy
  for (unsigned short input = 0; input <= 0xFF; input++) {
    // loop over ietas, barrel first
    for (unsigned short absIeta = 1; absIeta <= 17; absIeta++) {
      // positive eta
      ecal1 = ecalTPGScale.getTPGInGeV((unsigned int)input, EcalTrigTowerDetId(1, EcalBarrel, absIeta, 1));
      ecal2 = caloEcalScale->et(input, absIeta, 1);

      if (!(ecal1 == ecal2)) {
        ecalIsConsistent = false;
        /*LogWarning("InconsistentData") 
         << "ECAL scales not consistent! (pos eta barrel)"
         << "old ECAL is " << setprecision (8) << ecal1
         << " new ECAL is " << setprecision (8) << ecal2 
         << " difference is " << (ecal1 - ecal2) ; */
      }

      // negative eta
      ecal1 = ecalTPGScale.getTPGInGeV((unsigned int)input, EcalTrigTowerDetId(-1, EcalBarrel, absIeta, 2));
      ecal2 = caloEcalScale->et(input, absIeta, -1);

      if (!(ecal1 == ecal2)) {
        ecalIsConsistent = false;
        /*LogWarning("InconsistentData") 
         << "ECAL scales not consistent! (neg eta barrel)"
         << "old ECAL is " << setprecision (8) << ecal1
         << " new ECAL is " << setprecision (8) << ecal2           
         << " difference is " << (ecal1 - ecal2) ; */
      }
    }
    // now loop over endcap ietas
    for (unsigned short absIeta = 18; absIeta <= 28; absIeta++) {
      // positive eta
      ecal1 = ecalTPGScale.getTPGInGeV((unsigned int)input, EcalTrigTowerDetId(1, EcalEndcap, absIeta, 1));
      ecal2 = caloEcalScale->et(input, absIeta, 1);

      if (!(ecal1 == ecal2)) {
        ecalIsConsistent = false;
        /*LogWarning("InconsistentData") 
         << "ECAL scales not consistent! (pos eta endcap)"
         << "old ECAL is " << setprecision (8) << ecal1
         << " new ECAL is " << setprecision (8) << ecal2 
         << " difference is " << (ecal1 - ecal2) ; */
      }

      // negative eta
      ecal1 = ecalTPGScale.getTPGInGeV((unsigned int)input, EcalTrigTowerDetId(-1, EcalEndcap, absIeta, 2));
      ecal2 = caloEcalScale->et(input, absIeta, -1);

      if (!(ecal1 == ecal2)) {
        ecalIsConsistent = false;
        /*LogWarning("InconsistentData") 
         << "ECAL scales not consistent! (neg eta endcap)"
         << "old ECAL is " << setprecision (8) << ecal1
         << " new ECAL is " << setprecision (8) << ecal2 
         << " difference is " << (ecal1 - ecal2) ; */
      }
    }
  }

  if (!ecalIsConsistent) {
    // do something
    //cout << "WARNING: ECAL scales not consistent!" << endl;
    LogWarning("InconsistentData") << "ECAL scales not consistent!";
  } else {
    // do something else
    //cout << "ECAL scales okay" << endl;
  }

  // compare the hcal scales

  for (unsigned short input = 0; input <= 0xFF; input++) {
    // loop over ietas
    for (unsigned short absIeta = 1; absIeta <= 32; absIeta++) {
      hcal1 = caloTPGTranscoder->hcaletValue(absIeta, input);   // no eta-
      hcal2 = caloHcalScale->et(input, absIeta, 1);             // sign in transcoder
      hcal3 = caloTPGTranscoder->hcaletValue(-absIeta, input);  // no eta-
      hcal4 = caloHcalScale->et(input, absIeta, -1);            // sign in transcoder

      if ((!(hcal1 == hcal2)) || (!(hcal3 == hcal4))) {
        hcalIsConsistent = false;
        /*LogWarning("InconsistentData") 
         << "HCAL scales not consistent!"
         << "old HCAL is " << hcal1
         << " new HCAL is " << hcal2 ; */
      }
    }
  }
  if (!hcalIsConsistent) {
    // do something
    //cout << "WARNING: HCAL scales not consistent!" << endl;
    LogWarning("InconsistentData") << "HCAL scales not consistent!";
  } else {
    // do something else
    //cout << "HCAL scales okay" << endl;
  }
}

// ------------ method called once each job just before starting event loop  ------------
void L1CaloInputScaleTester::beginJob() {}

// ------------ method called once each job just after ending the event loop  ------------
void L1CaloInputScaleTester::endJob() {}