L1ABCDebugger.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    SiStripTools
// Class:      L1ABCDebugger
//
//
// Original Author:  Andrea Venturi
//         Created:  Tue Jul 19 11:56:00 CEST 2009
//
//

// system include files
#include <memory>

// user include files
#include "TH2F.h"
#include "TProfile.h"
#include "DPGAnalysis/SiStripTools/interface/RunHistogramManager.h"

#include <vector>
#include <string>

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/Run.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "FWCore/Utilities/interface/InputTag.h"

#include "DataFormats/Scalers/interface/L1AcceptBunchCrossing.h"
//
// class decleration
//

class L1ABCDebugger : public edm::EDAnalyzer {
public:
  explicit L1ABCDebugger(const edm::ParameterSet&);
  ~L1ABCDebugger() override;

private:
  void beginJob() override;
  void analyze(const edm::Event&, const edm::EventSetup&) override;
  void beginRun(const edm::Run&, const edm::EventSetup&) override;
  void endJob() override;

  // ----------member data ---------------------------

  edm::EDGetTokenT<L1AcceptBunchCrossingCollection> m_l1abccollectionToken;
  const unsigned int m_maxLS;
  const unsigned int m_LSfrac;

  RunHistogramManager m_rhm;

  TH2F** m_hoffsets;
  TProfile** m_horboffvsorb;
  TProfile** m_hbxoffvsorb;
};

//
// constants, enums and typedefs
//

//
// static data member definitions
//

//
// constructors and destructor
//
L1ABCDebugger::L1ABCDebugger(const edm::ParameterSet& iConfig)
    : m_l1abccollectionToken(
          consumes<L1AcceptBunchCrossingCollection>(iConfig.getParameter<edm::InputTag>("l1ABCCollection"))),
      m_maxLS(iConfig.getUntrackedParameter<unsigned int>("maxLSBeforeRebin", 250)),
      m_LSfrac(iConfig.getUntrackedParameter<unsigned int>("startingLSFraction", 16)),
      m_rhm(consumesCollector()) {
  //now do what ever initialization is needed

  m_hoffsets = m_rhm.makeTH2F(
      "offsets", "Orbit vs BX offsets between SCAL and Event", 2 * 3564 + 1, -3564.5, 3564.5, 201, -100.5, 100.5);
  m_horboffvsorb =
      m_rhm.makeTProfile("orboffvsorb", "SCAL Orbit offset vs orbit number", m_LSfrac * m_maxLS, 0, m_maxLS * 262144);
  m_hbxoffvsorb =
      m_rhm.makeTProfile("bxoffvsorb", "SCAL BX offset vs orbit number", m_LSfrac * m_maxLS, 0, m_maxLS * 262144);
}

L1ABCDebugger::~L1ABCDebugger() {
  // 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 L1ABCDebugger::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
  using namespace edm;

  Handle<L1AcceptBunchCrossingCollection> pIn;
  iEvent.getByToken(m_l1abccollectionToken, pIn);

  // offset computation
  for (L1AcceptBunchCrossingCollection::const_iterator l1abc = pIn->begin(); l1abc != pIn->end(); ++l1abc) {
    if (l1abc->l1AcceptOffset() == 0) {
      if (m_hoffsets && *m_hoffsets)
        (*m_hoffsets)
            ->Fill((int)l1abc->bunchCrossing() - (int)iEvent.bunchCrossing(),
                   (long long)l1abc->orbitNumber() - (long long)iEvent.orbitNumber());
      if (m_horboffvsorb && *m_horboffvsorb)
        (*m_horboffvsorb)->Fill(iEvent.orbitNumber(), (long long)l1abc->orbitNumber() - (long long)iEvent.orbitNumber());
      if (m_hbxoffvsorb && *m_hbxoffvsorb)
        (*m_hbxoffvsorb)->Fill(iEvent.orbitNumber(), (int)l1abc->bunchCrossing() - (int)iEvent.bunchCrossing());
    }
  }

  // dump of L1ABC collection

  edm::LogInfo("L1ABCDebug") << "Dump of L1AcceptBunchCrossing Collection for event in orbit " << iEvent.orbitNumber()
                             << " and BX " << iEvent.bunchCrossing();

  for (L1AcceptBunchCrossingCollection::const_iterator l1abc = pIn->begin(); l1abc != pIn->end(); ++l1abc) {
    edm::LogVerbatim("L1ABCDebug") << *l1abc;
  }
}

void L1ABCDebugger::beginRun(const edm::Run& iRun, const edm::EventSetup& iSetup) {
  m_rhm.beginRun(iRun);

  if (m_hoffsets && *m_hoffsets) {
    (*m_hoffsets)->GetXaxis()->SetTitle("#Delta BX (SCAL-Event)");
    (*m_hoffsets)->GetYaxis()->SetTitle("#Delta orbit (SCAL-Event)");
  }
  if (m_horboffvsorb && *m_horboffvsorb) {
    (*m_horboffvsorb)->GetXaxis()->SetTitle("Orbit");
    (*m_horboffvsorb)->GetYaxis()->SetTitle("#Delta orbit (SCAL-Event)");
    (*m_horboffvsorb)->SetCanExtend(TH1::kXaxis);
  }
  if (m_hbxoffvsorb && *m_hbxoffvsorb) {
    (*m_hbxoffvsorb)->GetXaxis()->SetTitle("Orbit");
    (*m_hbxoffvsorb)->GetYaxis()->SetTitle("#Delta BX (SCAL-Event)");
    (*m_hbxoffvsorb)->SetCanExtend(TH1::kXaxis);
  }
}
// ------------ method called once each job just before starting event loop  ------------
void L1ABCDebugger::beginJob() {}

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

//define this as a plug-in
DEFINE_FWK_MODULE(L1ABCDebugger);