TriggerValidator.h

Go to the documentation of this file.

#ifndef SusyBsmTriggerPerformance_TriggerValidator_TriggerValidator_h
#define SusyBsmTriggerPerformance_TriggerValidator_TriggerValidator_h

// -*- C++ -*-
//
// Package:    TriggerValidator
// Class:      TriggerValidator
// 
//
// Original Author:  Massimiliano Chiorboli
//                   Maurizio Pierini
//                   Maria Spiropulu
//         Created:  Wed Aug 29 15:10:56 CEST 2007
//
//

// system include files
#include <memory>
#include <fstream>
#include <iostream>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "HLTrigger/HLTcore/interface/HLTConfigProvider.h"

#include "HLTriggerOffline/SUSYBSM/interface/RecoSelector.h"
#include "HLTriggerOffline/SUSYBSM/interface/McSelector.h"
#include "HLTriggerOffline/SUSYBSM/interface/PlotMakerL1.h"
#include "HLTriggerOffline/SUSYBSM/interface/PlotMakerReco.h"
#include "HLTriggerOffline/SUSYBSM/interface/MuonAnalyzer.h"
//#include "HLTriggerOffline/SUSYBSM/interface/TurnOnMaker.h"

//data formats
//Added by Max for the Trigger
#include "DataFormats/HLTReco/interface/TriggerEventWithRefs.h"
#include "DataFormats/Common/interface/TriggerResults.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutRecord.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutRecord.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerObjectMapRecord.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerObjectMap.h"


//included for DQM
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"




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

      TFile* theHistoFile;

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


      void beginRun(const edm::Run& run, const edm::EventSetup& c);
      void endRun(const edm::Run& run, const edm::EventSetup& c);

      
      std::map<int,std::string> l1NameMap;
      // ----------member data ---------------------------
      DQMStore * dbe_;
      std::string dirname_;
      HLTConfigProvider hltConfig_;
      // --- names of the folders in the dbe ---
      std::string triggerBitsDir ;
      std::string recoSelBitsDir ; 
      std::string mcSelBitsDir   ;      



      unsigned int nHltPaths;
      int nL1Bits;


      //RecoSelector
      std::vector<RecoSelector*> myRecoSelector;
      std::vector<McSelector*> myMcSelector;
      
      //For the moment I switch off the more complex plots
       PlotMakerL1* myPlotMakerL1; 
       PlotMakerReco* myPlotMakerReco; 
/*       TurnOnMaker* myTurnOnMaker; */

      //Histo
      std::string HistoFileName;
      std::string StatFileName;
      edm::EDGetTokenT<L1GlobalTriggerObjectMapRecord> l1Label;
      edm::InputTag hltLabel;
      edm::EDGetTokenT<edm::TriggerResults> hlt_token_;

      //McFlag
      bool mcFlag;
      bool l1Flag;

      //Cut parameters decided by the user
      std::vector<edm::ParameterSet> reco_parametersets;
      std::vector<edm::ParameterSet> mc_parametersets;
      edm::ParameterSet turnOn_params;
      edm::ParameterSet plotMakerL1Input;
      edm::ParameterSet plotMakerRecoInput;

      edm::InputTag muonTag_;
      edm::InputTag triggerTag_;
      std::string processName_;
      std::string triggerName_;
      edm::EDGetTokenT<L1GlobalTriggerReadoutRecord> gtDigis_token_;

      // name of each L1 algorithm
      std::vector<std::string> l1Names_;    
      // name of each hlt algorithm
      std::vector<std::string>  hlNames_;  


      //Counters for L1 and HLT
      std::vector<int> numTotL1BitsBeforeCuts;
      std::vector<int> numTotHltBitsBeforeCuts;
      std::vector< std::vector<int> > numTotL1BitsAfterRecoCuts;
      std::vector< std::vector<int> > numTotHltBitsAfterRecoCuts;
      std::vector< std::vector<int> > numTotL1BitsAfterMcCuts;
      std::vector< std::vector<int> > numTotHltBitsAfterMcCuts;

      std::vector<double> effL1BeforeCuts;
      std::vector<double> effHltBeforeCuts;
      std::vector<double> effL1AfterRecoCuts;
      std::vector<double> effHltAfterRecoCuts;
      std::vector<double> effL1AfterMcCuts;
      std::vector<double> effHltAfterMcCuts;
      
      

      std::vector< std::vector<int> > vCorrL1;
      std::vector< std::vector<int> > vCorrHlt;
      std::vector< std::vector<double> > vCorrNormL1;
      std::vector< std::vector<double> > vCorrNormHlt;

      int nEvTot;
      std::vector<int> nEvRecoSelected;
      std::vector<int> nEvMcSelected;



      
      //Histos for L1 and HLT bits
      MonitorElement* hL1BitsBeforeCuts;    
      MonitorElement* hHltBitsBeforeCuts;   
      std::vector<MonitorElement*> hL1BitsAfterRecoCuts;    
      std::vector<MonitorElement*> hHltBitsAfterRecoCuts;  
      std::vector<MonitorElement*> hL1BitsAfterMcCuts;  
      std::vector<MonitorElement*> hHltBitsAfterMcCuts;  

      MonitorElement* hL1PathsBeforeCuts;
      MonitorElement* hHltPathsBeforeCuts;
      std::vector<MonitorElement*> hL1PathsAfterRecoCuts;
      std::vector<MonitorElement*> hHltPathsAfterRecoCuts;
      std::vector<MonitorElement*> hL1PathsAfterMcCuts;
      std::vector<MonitorElement*> hHltPathsAfterMcCuts;
        

      //if we want to keep these, probably thay have to move to the client      
/*       TH2D* hL1OverlapNormToTotal; */
/*       TH2D* hHltOverlapNormToTotal; */
/*       TH2D* hL1OverlapNormToLargestPath; */
/*       TH2D* hHltOverlapNormToLargestPath; */

      std::vector<int> l1bits;
      std::vector<int> hltbits;

      bool firstEvent;

      MuonAnalyzerSBSM* myMuonAnalyzer;
};



#endif