---

HLTTauL1DQMOfflineSource.h

Go to the documentation of this file.

// Original Author:  Chi Nhan Nguyen
//         Created:  Fri Feb 22 09:20:55 CST 2008

// system include files
#include <memory>
#include <vector>
#include <string>
#include <iostream>
#include <iomanip>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

// L1 Trigger data formats
#include "DataFormats/L1Trigger/interface/L1EmParticle.h"
#include "DataFormats/L1Trigger/interface/L1EmParticleFwd.h"
#include "DataFormats/L1Trigger/interface/L1JetParticle.h"
#include "DataFormats/L1Trigger/interface/L1JetParticleFwd.h"
#include "DataFormats/L1Trigger/interface/L1MuonParticle.h"
#include "DataFormats/L1Trigger/interface/L1MuonParticleFwd.h"
#include "DataFormats/L1Trigger/interface/L1EtMissParticle.h"
#include "DataFormats/L1Trigger/interface/L1EtMissParticleFwd.h"
#include "DataFormats/L1Trigger/interface/L1ParticleMap.h"
#include "DataFormats/L1Trigger/interface/L1ParticleMapFwd.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutRecord.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutSetupFwd.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerObjectMapRecord.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerObjectMapFwd.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerObjectMap.h"
#include "DataFormats/L1GlobalMuonTrigger/interface/L1MuGMTExtendedCand.h"

#include "FWCore/ServiceRegistry/interface/Service.h" // Framework services
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"

//
typedef math::XYZTLorentzVectorD   LV;
typedef std::vector<LV>            LVColl;


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


   private:
      virtual void beginJob(const edm::EventSetup&) ;
      virtual void analyze(const edm::Event&, const edm::EventSetup&);
      virtual void endJob() ;


  void getL1extraObjects(const edm::Event&);
  void evalL1extraDecisions();
  void fillL1Histograms();
  void fillL1MCTauMatchedHists(const edm::Event& iEvent);

  //  void convertToIntegratedEff(MonitorElement*,double);
  
  // ----------member data ---------------------------

  edm::InputTag     _refTauColl; 
  edm::InputTag     _refElecColl; 
  edm::InputTag     _refMuonColl; 
  
  edm::InputTag _L1extraTauJetSource;
  edm::InputTag _L1extraCenJetSource;
  edm::InputTag _L1extraForJetSource;
  edm::InputTag _L1extraMuonSource;
  edm::InputTag _L1extraMETSource;
  edm::InputTag _L1extraNonIsoEgammaSource;
  edm::InputTag _L1extraIsoEgammaSource;

  // Thresholds of L1 menu
  double _SingleTauThreshold;
  double _DoubleTauThreshold;
  std::vector<double> _SingleTauMETThresholds;
  std::vector<double> _MuTauThresholds;
  std::vector<double> _IsoEgTauThresholds;

  // Cuts
  double _L1MCTauMinDeltaR;
  double _MCTauHadMinEt;
  double _MCTauHadMaxAbsEta;

  double _L1MCElecMinDeltaR;
  double _MCElecMinEt;
  double _MCElecMaxAbsEta;

  double _L1MCMuonMinDeltaR;
  double _MCMuonMinEt;
  double _MCMuonMaxAbsEta;

  //Output file
  std::string _triggerTag;//tag for dqm flder
  std::string _outFile;

  // L1extra Objects
  LVColl _L1Taus;
  LVColl _L1CenJets;
  LVColl _L1ForJets;
  LVColl _L1NonIsoEgammas;
  LVColl _L1IsoEgammas;
  LVColl _L1METs;
  LVColl _L1Muons;
  std::vector<int> _L1MuQuals;
  
  // histograms
  MonitorElement* h_L1TauEt;  
  MonitorElement* h_L1TauEta;
  MonitorElement* h_L1TauPhi;

  MonitorElement* h_L1Tau1Et;
  MonitorElement* h_L1Tau1Eta;
  MonitorElement* h_L1Tau1Phi;

  MonitorElement* h_L1Tau2Et;
  MonitorElement* h_L1Tau2Eta;
  MonitorElement* h_L1Tau2Phi;

  MonitorElement* h_L1IsoEg1Et;  
  MonitorElement* h_L1IsoEg1Eta;
  MonitorElement* h_L1IsoEg1Phi;

  MonitorElement* h_L1Muon1Et;  
  MonitorElement* h_L1Muon1Eta;
  MonitorElement* h_L1Muon1Phi;

  MonitorElement* h_L1Met;  
  MonitorElement* h_L1MetEta;  
  MonitorElement* h_L1MetPhi;  

  // L1 response
  MonitorElement* h_L1MCTauDeltaR;
  MonitorElement* h_L1minusMCTauEt;
  MonitorElement* h_L1minusMCoverMCTauEt;

  // MC w/o cuts
  MonitorElement* h_GenTauHadEt;
  MonitorElement* h_GenTauHadEta;
  MonitorElement* h_GenTauHadPhi;

  MonitorElement* h_GenTauElecEt;
  MonitorElement* h_GenTauElecEta;
  MonitorElement* h_GenTauElecPhi;

  MonitorElement* h_GenTauMuonEt;
  MonitorElement* h_GenTauMuonEta;
  MonitorElement* h_GenTauMuonPhi;

  // Tau -> Electron
  // MC matching efficiencies
  //  MonitorElement* h_EffMCTauElecEt;
  //MonitorElement* h_EffMCTauElecEta;
  //MonitorElement* h_EffMCTauElecPhi;
  // Numerators
  MonitorElement* h_L1MCMatchedTauElecEt;
  MonitorElement* h_L1MCMatchedTauElecEta;
  MonitorElement* h_L1MCMatchedTauElecPhi;
  // Denominators
  MonitorElement* h_MCTauElecEt;
  MonitorElement* h_MCTauElecEta;
  MonitorElement* h_MCTauElecPhi;

  // Tau -> Muon
  // MC matching efficiencies
  //MonitorElement* h_EffMCTauMuonEt;
  //MonitorElement* h_EffMCTauMuonEta;
  //MonitorElement* h_EffMCTauMuonPhi;
  // Numerators
  MonitorElement* h_L1MCMatchedTauMuonEt;
  MonitorElement* h_L1MCMatchedTauMuonEta;
  MonitorElement* h_L1MCMatchedTauMuonPhi;
  // Denominators
  MonitorElement* h_MCTauMuonEt;
  MonitorElement* h_MCTauMuonEta;
  MonitorElement* h_MCTauMuonPhi;

  // Tau -> Hadr
  // MC matching efficiencies
  //MonitorElement* h_EffMCTauEt;
  //MonitorElement* h_EffMCTauEta;
  //MonitorElement* h_EffMCTauPhi;
  // Numerators
  MonitorElement* h_L1MCMatchedTauEt;
  MonitorElement* h_L1MCMatchedTauEta;
  MonitorElement* h_L1MCMatchedTauPhi;
  // Denominators
  MonitorElement* h_MCTauHadEt;
  MonitorElement* h_MCTauHadEta;
  MonitorElement* h_MCTauHadPhi;
  
  // Event based efficiencies as a function of thresholds
  MonitorElement* h_L1SingleTauEffEt;
  MonitorElement* h_L1DoubleTauEffEt;
  MonitorElement* h_L1SingleTauEffMCMatchEt;
  MonitorElement* h_L1DoubleTauEffMCMatchEt;

  MonitorElement* h_L1TauMETfixEffEt;
  MonitorElement* h_L1TauMETfixEffMCMatchEt;
  MonitorElement* h_L1METTaufixEffEt;
  MonitorElement* h_L1METTaufixEffMCMatchEt;

  MonitorElement* h_L1TauIsoEgfixEffEt;
  MonitorElement* h_L1TauIsoEgfixEffMCMatchEt;
  MonitorElement* h_L1IsoEgTaufixEffEt;
  MonitorElement* h_L1IsoEgTaufixEffMCMatchEt;

  MonitorElement* h_L1TauMuonfixEffEt;
  MonitorElement* h_L1TauMuonfixEffMCMatchEt;
  MonitorElement* h_L1MuonTaufixEffEt;
  MonitorElement* h_L1MuonTaufixEffMCMatchEt;


  //Counter Histogram
  MonitorElement* h_nfidCounter;

  // Counters for event based efficiencies
  int _nEvents; // all events processed

  int _nfidEventsGenTauHad; 
  int _nfidEventsDoubleGenTauHads; 
  int _nfidEventsGenTauMuonTauHad; 
  int _nfidEventsGenTauElecTauHad;   

};

//

---