#include <SUSYDQMAnalyzer.h>

Inheritance diagram for SUSYDQMAnalyzer:

Public Member Functions
	SUSYDQMAnalyzer (const edm::ParameterSet &)

	~SUSYDQMAnalyzer () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

Private Attributes
double	_maxAbsEta

double	_maxNJets

double	_ptThreshold

MonitorElement *	hCaloAlpha_T

MonitorElement *	hCaloHT

MonitorElement *	hCaloMET

MonitorElement *	hCaloMHT

MonitorElement *	hPFAlpha_T

MonitorElement *	hPFHT

MonitorElement *	hPFMET

MonitorElement *	hPFMHT

edm::ParameterSet	iConfig

std::string	SUSYFolder

edm::EDGetTokenT< reco::CaloJetCollection >	theCaloJetCollectionToken

edm::EDGetTokenT< reco::CaloMETCollection >	theCaloMETCollectionToken

edm::EDGetTokenT< std::vector< reco::PFJet > >	thePFJetCollectionToken

edm::EDGetTokenT< reco::PFMETCollection >	thePFMETCollectionToken

Static Private Attributes
static const char *	messageLoggerCatregory = "SUSYDQM"

Detailed Description

Definition at line 26 of file SUSYDQMAnalyzer.h.

Constructor & Destructor Documentation

SUSYDQMAnalyzer::SUSYDQMAnalyzer ( const edm::ParameterSet & pSet )

explicit

Definition at line 72 of file SUSYDQMAnalyzer.cc.

References edm::ParameterSet::getParameter(), messageLoggerCatregory, and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   iConfig = pSet;
   
   SUSYFolder = iConfig.getParameter<std::string>("folderName");
   // Load parameters 
   thePFMETCollectionToken     = consumes<reco::PFMETCollection>   (iConfig.getParameter<edm::InputTag>("PFMETCollectionLabel"));
   theCaloMETCollectionToken   = consumes<reco::CaloMETCollection> (iConfig.getParameter<edm::InputTag>("CaloMETCollectionLabel"));
 
   //remove TCMET and JPT related variables due to anticipated changes in RECO
   //theTCMETCollectionToken     = consumes<reco::METCollection>     (iConfig.getParameter<edm::InputTag>("TCMETCollectionLabel"));
 
   theCaloJetCollectionToken   = consumes<reco::CaloJetCollection>   (iConfig.getParameter<edm::InputTag>("CaloJetCollectionLabel"));
   //theJPTJetCollectionToken    = consumes<reco::JPTJetCollection>    (iConfig.getParameter<edm::InputTag>("JPTJetCollectionLabel"));
   thePFJetCollectionToken     = consumes<std::vector<reco::PFJet> > (iConfig.getParameter<edm::InputTag>("PFJetCollectionLabel"));
 
   _ptThreshold = iConfig.getParameter<double>("ptThreshold");
   _maxNJets = iConfig.getParameter<double>("maxNJets");
   _maxAbsEta = iConfig.getParameter<double>("maxAbsEta");
 
 }

SUSYDQMAnalyzer::~SUSYDQMAnalyzer ( )

override

Definition at line 273 of file SUSYDQMAnalyzer.cc.

273 {

274 }

Member Function Documentation

void SUSYDQMAnalyzer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 141 of file SUSYDQMAnalyzer.cc.

References funct::abs(), edm::Event::getByToken(), edm::HandleBase::isValid(), metsig::jet, edm::Handle< T >::product(), reco::LeafCandidate::pt(), HT< T >::ScalarSum, and HT< T >::v.

 {
   //###########################################################
   // HTand MHT
   
   //===========================================================
   // Calo HT, MHT and alpha_T
 
   edm::Handle<reco::CaloJetCollection> CaloJetcoll;
 
   iEvent.getByToken(theCaloJetCollectionToken, CaloJetcoll);
 
   if(!CaloJetcoll.isValid()) return;
   
   std::vector<math::XYZTLorentzVector> Ps;
   for (reco::CaloJetCollection::const_iterator jet = CaloJetcoll->begin(); jet!=CaloJetcoll->end(); ++jet){
     if ((jet->pt()>_ptThreshold) && (abs(jet->eta()) < _maxAbsEta)){
       if(Ps.size()>_maxNJets) {
     edm::LogInfo(messageLoggerCatregory)<<"NMax Jets exceded..";
         break;
       }
       Ps.push_back(jet->p4());
     }
   }
 
   hCaloAlpha_T->Fill( alpha_T()(Ps));
 
   HT< reco::CaloJetCollection > CaloHT(CaloJetcoll, _ptThreshold, _maxAbsEta);
 
   hCaloHT->Fill(CaloHT.ScalarSum);
   hCaloMHT->Fill(CaloHT.v.Mod());
 
   //===========================================================
   // PF HT, MHT and alpha_T
 
   edm::Handle<reco::PFJetCollection> PFjetcoll;
 
   iEvent.getByToken(thePFJetCollectionToken, PFjetcoll);
 
   if(!PFjetcoll.isValid()) return;
 
   Ps.clear();
   for (reco::PFJetCollection::const_iterator jet = PFjetcoll->begin(); jet!=PFjetcoll->end(); ++jet){
     if ((jet->pt()>_ptThreshold) && (abs(jet->eta()) < _maxAbsEta)){
       if(Ps.size()>_maxNJets) {
     edm::LogInfo(messageLoggerCatregory)<<"NMax Jets exceded..";
     break;
       }
       Ps.push_back(jet->p4());
     }
   }
   hPFAlpha_T->Fill( alpha_T()(Ps));
 
   HT<reco::PFJetCollection> PFHT(PFjetcoll, _ptThreshold, _maxAbsEta);
 
   hPFHT->Fill(PFHT.ScalarSum);
   hPFMHT->Fill(PFHT.v.Mod());
   
   //===========================================================
   // JPT HT, MHT and alpha_T
 
   //edm::Handle<reco::JPTJetCollection> JPTjetcoll;
 
   //iEvent.getByToken(theJPTJetCollectionToken, JPTjetcoll);
 
   //if(!JPTjetcoll.isValid()) return;
 
   //Ps.clear();
   //for (reco::JPTJetCollection::const_iterator jet = JPTjetcoll->begin(); jet!=JPTjetcoll->end(); ++jet){
   //if ((jet->pt()>_ptThreshold) && (abs(jet->eta())<_maxAbsEta)){
   //  if(Ps.size()>_maxNJets) {
   //    edm::LogInfo(messageLoggerCatregory)<<"NMax Jets exceded..";
   //    break;
   //  }
   //  Ps.push_back(jet->p4());
   //}
   //}
   //hJPTAlpha_T->Fill( alpha_T()(Ps));
 
   //HT<reco::JPTJetCollection> JPTHT(JPTjetcoll, _ptThreshold, _maxAbsEta);
 
   //hJPTHT->Fill(JPTHT.ScalarSum);
   //hJPTMHT->Fill(JPTHT.v.Mod());
 
   //###########################################################
   // MET
 
   //===========================================================  
   // Calo MET
 
   edm::Handle<reco::CaloMETCollection> calometcoll;
   iEvent.getByToken(theCaloMETCollectionToken, calometcoll);
 
   if(!calometcoll.isValid()) return;
 
   const CaloMETCollection *calometcol = calometcoll.product();
   const CaloMET *calomet;
   calomet = &(calometcol->front());
   
   hCaloMET->Fill(calomet->pt());
 
   //===========================================================
   // PF MET
 
   edm::Handle<reco::PFMETCollection> pfmetcoll;
   iEvent.getByToken(thePFMETCollectionToken, pfmetcoll);
   
   if(!pfmetcoll.isValid()) return;
 
   const PFMETCollection *pfmetcol = pfmetcoll.product();
   const PFMET *pfmet;
   pfmet = &(pfmetcol->front());
 
   hPFMET->Fill(pfmet->pt());
 
   //===========================================================
   // TC MET
 
   //edm::Handle<reco::METCollection> tcmetcoll;
   //iEvent.getByToken(theTCMETCollectionToken, tcmetcoll);
   
   //if(!tcmetcoll.isValid()) return;
 
   //const METCollection *tcmetcol = tcmetcoll.product();
   //const MET *tcmet;
   //tcmet = &(tcmetcol->front());
 
   //hTCMET->Fill(tcmet->pt());
 
 }

void SUSYDQMAnalyzer::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	iRun,
		edm::EventSetup const &
	)

overrideprivate

Definition at line 97 of file SUSYDQMAnalyzer.cc.

References DQMStore::IBooker::book1D(), dir, DQMStore::IBooker::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                              {
   //  if( dqm ) {
     //===========================================================                                                                                 
     // book HT histos.                                                                                                                           
     std::string dir=SUSYFolder;
     dir+="HT"; 
     ibooker.setCurrentFolder(dir);
     hCaloHT = ibooker.book1D("Calo_HT", "", 500, 0., 2000);
     hPFHT   = ibooker.book1D("PF_HT"  , "", 500, 0., 2000);
     //hJPTHT  = ibooker.book1D("JPT_HT" , "", 500, 0., 2000);
     //===========================================================                                                                                 
     // book MET histos.                                                                                                                           
 
     dir=SUSYFolder;
     dir+="MET";
     ibooker.setCurrentFolder(dir);
     hCaloMET = ibooker.book1D("Calo_MET", "", 500, 0., 1000);
     hPFMET   = ibooker.book1D("PF_MET"  , "", 500, 0., 1000);
     //hTCMET   = ibooker.book1D("TC_MET"  , "", 500, 0., 1000);
 
     //===========================================================                                                                                 
     // book MHT histos.                                                                                                                           
 
     dir=SUSYFolder;
     dir+="MHT";
     ibooker.setCurrentFolder(dir);
     hCaloMHT = ibooker.book1D("Calo_MHT", "", 500, 0., 1000);
     hPFMHT   = ibooker.book1D("PF_MHT"  , "", 500, 0., 1000);
     //hJPTMHT  = ibooker.book1D("JPT_MHT" , "", 500, 0., 1000);
 
     //===========================================================                                                                                 
     // book alpha_T histos.                                                                                                                       
 
     dir=SUSYFolder;
     dir+="Alpha_T";
     ibooker.setCurrentFolder(dir);
     hCaloAlpha_T = ibooker.book1D("Calo_AlphaT", "", 100, 0., 1.);
     //hJPTAlpha_T  = ibooker.book1D("PF_AlphaT"  , "", 100, 0., 1.);
     hPFAlpha_T   = ibooker.book1D("PF_AlphaT"  , "", 100, 0., 1.);
     //  }
 }

Member Data Documentation

double SUSYDQMAnalyzer::_maxAbsEta

private

Definition at line 47 of file SUSYDQMAnalyzer.h.

double SUSYDQMAnalyzer::_maxNJets

private

Definition at line 46 of file SUSYDQMAnalyzer.h.

double SUSYDQMAnalyzer::_ptThreshold

private

Definition at line 45 of file SUSYDQMAnalyzer.h.

MonitorElement* SUSYDQMAnalyzer::hCaloAlpha_T

private

Definition at line 67 of file SUSYDQMAnalyzer.h.

MonitorElement* SUSYDQMAnalyzer::hCaloHT

private

Definition at line 55 of file SUSYDQMAnalyzer.h.

MonitorElement* SUSYDQMAnalyzer::hCaloMET

private

Definition at line 59 of file SUSYDQMAnalyzer.h.

MonitorElement* SUSYDQMAnalyzer::hCaloMHT

private

Definition at line 63 of file SUSYDQMAnalyzer.h.

MonitorElement* SUSYDQMAnalyzer::hPFAlpha_T

private

Definition at line 69 of file SUSYDQMAnalyzer.h.

MonitorElement* SUSYDQMAnalyzer::hPFHT

private

Definition at line 57 of file SUSYDQMAnalyzer.h.

MonitorElement* SUSYDQMAnalyzer::hPFMET

private

Definition at line 60 of file SUSYDQMAnalyzer.h.

MonitorElement* SUSYDQMAnalyzer::hPFMHT

private

Definition at line 65 of file SUSYDQMAnalyzer.h.

edm::ParameterSet SUSYDQMAnalyzer::iConfig

private

Definition at line 32 of file SUSYDQMAnalyzer.h.

const char * SUSYDQMAnalyzer::messageLoggerCatregory = "SUSYDQM"

staticprivate

Definition at line 50 of file SUSYDQMAnalyzer.h.

Referenced by SUSYDQMAnalyzer().

std::string SUSYDQMAnalyzer::SUSYFolder

private

Definition at line 49 of file SUSYDQMAnalyzer.h.

edm::EDGetTokenT<reco::CaloJetCollection> SUSYDQMAnalyzer::theCaloJetCollectionToken

private

Definition at line 39 of file SUSYDQMAnalyzer.h.

edm::EDGetTokenT<reco::CaloMETCollection> SUSYDQMAnalyzer::theCaloMETCollectionToken

private

Definition at line 41 of file SUSYDQMAnalyzer.h.

edm::EDGetTokenT<std::vector<reco::PFJet> > SUSYDQMAnalyzer::thePFJetCollectionToken

private

Definition at line 38 of file SUSYDQMAnalyzer.h.

edm::EDGetTokenT<reco::PFMETCollection> SUSYDQMAnalyzer::thePFMETCollectionToken

private

Definition at line 37 of file SUSYDQMAnalyzer.h.

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation