AlphaTVarAnalyzer.cc

Go to the documentation of this file.

#include "DQM/DataScouting/plugins/AlphaTVarAnalyzer.h"

#include "DataFormats/JetReco/interface/CaloJet.h"
#include "DataFormats/JetReco/interface/CaloJetCollection.h"

#include "DataFormats/EgammaCandidates/interface/Electron.h"
#include "DataFormats/EgammaCandidates/interface/ElectronFwd.h"

#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/RecoCandidate/interface/RecoChargedCandidate.h"

#include <cmath>

// A simple constructor which takes as input only the name of the PF jet collection
AlphaTVarAnalyzer::AlphaTVarAnalyzer( const edm::ParameterSet & conf ):
  ScoutingAnalyzerBase(conf),
  m_jetCollectionTag(conf.getUntrackedParameter<edm::InputTag>("jetCollectionName",edm::InputTag("hltCaloJetIDPassed"))),
  m_alphaTVarCollectionTag(conf.getUntrackedParameter<edm::InputTag>("alphaTVarCollectionName")) {
  //set Token(-s)
  m_alphaTVarCollectionTagToken_ = consumes<std::vector<double> >(conf.getUntrackedParameter<edm::InputTag>("alphaTVarCollectionName"));
}

AlphaTVarAnalyzer::~AlphaTVarAnalyzer() {}

// Function to book the Monitoring Elements.
void AlphaTVarAnalyzer::bookHistograms(DQMStore::IBooker & iBooker, edm::Run const &, edm::EventSetup const &) {
  ScoutingAnalyzerBase::prepareBooking(iBooker);
  //the full inclusive histograms
  m_HTAlphaT = bookH2withSumw2(
      iBooker,
      "HTvsAlphaT",
      "H_{T} vs #alpha_{T} (All Events)",
      400,0.,4000.,
      50,0.,1.,
      "H_{T} [GeV]",
      "#alpha_{T}");
  m_HTAlphaTg0p55 = bookH1withSumw2(
      iBooker,
      "HTvsAlphaTg0p55",
      "H_{T} (#alpha_{T} > 0.55)",
      400,0.,4000.,
      "H_{T} [GeV]");
  m_HTAlphaTg0p60 = bookH1withSumw2(
      iBooker,
      "HTvsAlphaTg0p60",
      "H_{T} (#alpha_{T} > 0.60)",
      400,0.,4000.,
      "H_{T} [GeV]");
} 

// Usual analyze method
void AlphaTVarAnalyzer::analyze( const edm::Event & iEvent, const edm::EventSetup & c ) {
  edm::Handle<std::vector<double> > alphaTvar_handle;
  iEvent.getByToken(m_alphaTVarCollectionTagToken_, alphaTvar_handle);

  if(alphaTvar_handle->size() > 1){
    const double AlphaT = alphaTvar_handle->at(0);
    const double HT = alphaTvar_handle->at(1);
    m_HTAlphaT->Fill(HT,AlphaT);
    if(AlphaT > 0.55) m_HTAlphaTg0p55->Fill(HT);
    if(AlphaT > 0.60) m_HTAlphaTg0p60->Fill(HT);
  }
}