L1ValidatorHists.h

Go to the documentation of this file.

#ifndef L1ValidatorHists_h
#define L1ValidatorHists_h

#include "TH1F.h"
//#include <TH2F.h>

//#include <DataFormats/HepMCCandidate/interface/GenParticle.h>
#include <DataFormats/Candidate/interface/LeafCandidate.h>

#include <string>

#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"

/*#define DECLARE_HISTS(TYPE) \
TH2F * ## TYPE ## _N_Pt; \
TH2F * ## TYPE ## _N_Eta; \
TH2F * ## TYPE ## _Eff_Pt; \
TH2F * ## TYPE ## _Eff_Eta; \
TH2F * ## TYPE ## _dR; \
TH2F * ## TYPE ## _dPt;
*/

class L1ValidatorHists {
  public:
    L1ValidatorHists();
    ~L1ValidatorHists();

    void Book(DQMStore::IBooker &, std::string dirname);
    void Write();

    int NEvents;

    //enum Type{ IsoEG, NonIsoEG, CenJet, ForJet, TauJet, Muon, Number };
    enum Type{ Egamma, Jet, Tau, Muon, Number };
    std::string Name[4];

    // Each object has gen (pt, eta, phi, pdg), reco (pt, eta, phi, pdg)
    // split by pdg (CenJet, ETM, ETT, ForJet, HTM, HTT, IsoEG, Mu, NoIsoEG, TauJet): 10
    // needs efficiency (by pt, eta), fake rate (by pt, eta), dR (by pt, (eta?))
    /*DECLARE_HISTS(IsoEG)
    DECLARE_HISTS(NonIsoEG)
    DECLARE_HISTS(CenJet)
    DECLARE_HISTS(ForJet)
    DECLARE_HISTS(TauJet)
    DECLARE_HISTS(Muon)*/

    MonitorElement *N[Type::Number];

    MonitorElement *Eff_Pt[Type::Number];
    MonitorElement *Eff_Pt_Denom[Type::Number];
    MonitorElement *Eff_Pt_Nomin[Type::Number];
    MonitorElement *Eff_Eta[Type::Number];
    MonitorElement *Eff_Eta_Denom[Type::Number];
    MonitorElement *Eff_Eta_Nomin[Type::Number];
    MonitorElement *TurnOn_15[Type::Number];
    MonitorElement *TurnOn_15_Denom[Type::Number];
    MonitorElement *TurnOn_15_Nomin[Type::Number];
    MonitorElement *TurnOn_30[Type::Number];
    MonitorElement *TurnOn_30_Denom[Type::Number];
    MonitorElement *TurnOn_30_Nomin[Type::Number];
    MonitorElement *dR[Type::Number];
    MonitorElement *dPt[Type::Number];
    MonitorElement *dR_vs_Pt[Type::Number];
    MonitorElement *dPt_vs_Pt[Type::Number];


    // add the rest...
    //TH2F *ETM_Delta, *ETT_Delta, *HTM_Delta, *HTT_Delta;

    void Fill(int, const reco::LeafCandidate *, const reco::LeafCandidate *);
    void FillNumber(int, int);
//  private:
//    void NormalizeSlices(TH2F *Hist);
};


#endif