#include <SUSY_HLT_DoubleEle_Hadronic.h>

Inheritance diagram for SUSY_HLT_DoubleEle_Hadronic:

Public Member Functions
	SUSY_HLT_DoubleEle_Hadronic (const edm::ParameterSet &ps)

	~SUSY_HLT_DoubleEle_Hadronic () override

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

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

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

	~DQMEDAnalyzer () override=default

Protected Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &eSetup) override

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

void	dqmBeginRun (edm::Run const &, edm::EventSetup const &) override

void	endRun (edm::Run const &run, edm::EventSetup const &eSetup) override

Private Member Functions
void	bookHistos (DQMStore::IBooker &)

Private Attributes
double	etaThrJet_

HLTConfigProvider	fHltConfig

MonitorElement *	h_EleTurnOn_den

MonitorElement *	h_EleTurnOn_num

MonitorElement *	h_pfHTTurnOn_den

MonitorElement *	h_pfHTTurnOn_num

MonitorElement *	h_triggerDoubleEleMass

MonitorElement *	h_triggerEleEta

MonitorElement *	h_triggerElePhi

MonitorElement *	h_triggerElePt

std::string	HLTProcess_

double	ptThrJet_

edm::EDGetTokenT< reco::CaloJetCollection >	theCaloJetCollection_

edm::EDGetTokenT< reco::GsfElectronCollection >	theElectronCollection_

edm::EDGetTokenT< reco::PFJetCollection >	thePfJetCollection_

edm::EDGetTokenT< trigger::TriggerEvent >	theTrigSummary_

edm::InputTag	triggerFilter_

std::string	triggerPath_

std::string	triggerPathAuxiliaryForElectron_

std::string	triggerPathAuxiliaryForHadronic_

edm::EDGetTokenT< edm::TriggerResults >	triggerResults_

Detailed Description

Definition at line 35 of file SUSY_HLT_DoubleEle_Hadronic.h.

Constructor & Destructor Documentation

SUSY_HLT_DoubleEle_Hadronic::SUSY_HLT_DoubleEle_Hadronic ( const edm::ParameterSet & ps )

Definition at line 9 of file SUSY_HLT_DoubleEle_Hadronic.cc.

References etaThrJet_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLTProcess_, ptThrJet_, AlCaHLTBitMon_QueryRunRegistry::string, theCaloJetCollection_, theElectronCollection_, thePfJetCollection_, theTrigSummary_, triggerFilter_, triggerPath_, triggerPathAuxiliaryForElectron_, triggerPathAuxiliaryForHadronic_, and triggerResults_.

 {
     edm::LogInfo("SUSY_HLT_DoubleEle_Hadronic") << "Constructor SUSY_HLT_DoubleEle_Hadronic::SUSY_HLT_DoubleEle_Hadronic " << std::endl;
     // Get parameters from configuration file
     theTrigSummary_ = consumes<trigger::TriggerEvent>(ps.getParameter<edm::InputTag>("trigSummary"));
     theElectronCollection_ = consumes<reco::GsfElectronCollection>(ps.getParameter<edm::InputTag>("ElectronCollection"));
     thePfJetCollection_ = consumes<reco::PFJetCollection>(ps.getParameter<edm::InputTag>("pfJetCollection"));
     theCaloJetCollection_ = consumes<reco::CaloJetCollection>(ps.getParameter<edm::InputTag>("caloJetCollection"));
     triggerResults_ = consumes<edm::TriggerResults>(ps.getParameter<edm::InputTag>("TriggerResults"));
     HLTProcess_ = ps.getParameter<std::string>("HLTProcess");
     triggerPath_ = ps.getParameter<std::string>("TriggerPath");
     triggerPathAuxiliaryForElectron_ = ps.getParameter<std::string>("TriggerPathAuxiliaryForElectron");
     triggerPathAuxiliaryForHadronic_ = ps.getParameter<std::string>("TriggerPathAuxiliaryForHadronic");
     triggerFilter_ = ps.getParameter<edm::InputTag>("TriggerFilter");
     ptThrJet_ = ps.getUntrackedParameter<double>("PtThrJet");
     etaThrJet_ = ps.getUntrackedParameter<double>("EtaThrJet");
 }

SUSY_HLT_DoubleEle_Hadronic::~SUSY_HLT_DoubleEle_Hadronic ( )

override

Definition at line 27 of file SUSY_HLT_DoubleEle_Hadronic.cc.

 {
     edm::LogInfo("SUSY_HLT_DoubleEle_Hadronic") << "Destructor SUSY_HLT_DoubleEle_Hadronic::~SUSY_HLT_DoubleEle_Hadronic " << std::endl;
 }

Member Function Documentation

void SUSY_HLT_DoubleEle_Hadronic::analyze	(	edm::Event const &	e,
		edm::EventSetup const &	eSetup
	)

overrideprotected

Definition at line 71 of file SUSY_HLT_DoubleEle_Hadronic.cc.

References edm::HLTGlobalStatus::accept(), nanoDQM_cfi::Electron, trigger::TriggerObject::eta(), etaThrJet_, MonitorElement::Fill(), HcalObjRepresent::Fill(), trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), edm::Event::getByToken(), trigger::TriggerEvent::getObjects(), h_EleTurnOn_den, h_EleTurnOn_num, h_pfHTTurnOn_den, h_pfHTTurnOn_num, h_triggerDoubleEleMass, h_triggerEleEta, h_triggerElePhi, h_triggerElePt, trigger::TriggerObject::id(), edm::HandleBase::isValid(), relativeConstraints::keys, trigger::TriggerObject::phi(), trigger::TriggerObject::pt(), ptThrJet_, edm::TriggerNames::size(), trigger::TriggerEvent::sizeFilters(), mathSSE::sqrt(), theCaloJetCollection_, theElectronCollection_, thePfJetCollection_, theTrigSummary_, triggerFilter_, edm::TriggerNames::triggerName(), edm::Event::triggerNames(), TriggerAnalyzer::triggerObjects, triggerPath_, triggerPathAuxiliaryForElectron_, triggerPathAuxiliaryForHadronic_, triggerResults_, trigNames, edm::HLTGlobalStatus::wasrun(), and x.

                                                                                        {
     edm::LogInfo("SUSY_HLT_DoubleEle_Hadronic") << "SUSY_HLT_DoubleEle_Hadronic::analyze" << std::endl;
     
     //-------------------------------
     //--- Jets
     //-------------------------------
     edm::Handle<reco::PFJetCollection> pfJetCollection;
     e.getByToken (thePfJetCollection_,pfJetCollection);
     if ( !pfJetCollection.isValid() ){
         edm::LogError ("SUSY_HLT_DoubleEle_Hadronic") << "invalid collection: PFJets" << "\n";
         return;
     }
     edm::Handle<reco::CaloJetCollection> caloJetCollection;
     e.getByToken (theCaloJetCollection_,caloJetCollection);
     if ( !caloJetCollection.isValid() ){
         edm::LogError ("SUSY_HLT_DoubleEle_Hadronic") << "invalid collection: CaloJets" << "\n";
         return;
     }
     
     //-------------------------------
     //--- Electron
     //-------------------------------
     edm::Handle<reco::GsfElectronCollection> ElectronCollection;
     e.getByToken (theElectronCollection_, ElectronCollection);
     if ( !ElectronCollection.isValid() ){
         edm::LogError ("SUSY_HLT_DoubleEle_Hadronic") << "invalid collection: Electrons " << "\n";
         return;
     }
     
     
     //-------------------------------
     //--- Trigger
     //-------------------------------
     edm::Handle<edm::TriggerResults> hltresults;
     e.getByToken(triggerResults_,hltresults);
     if(!hltresults.isValid()){
         edm::LogError ("SUSY_HLT_DoubleEle_Hadronic") << "invalid collection: TriggerResults" << "\n";
         return;
     }
     edm::Handle<trigger::TriggerEvent> triggerSummary;
     e.getByToken(theTrigSummary_, triggerSummary);
     if(!triggerSummary.isValid()) {
         edm::LogError ("SUSY_HLT_DoubleEle_Hadronic") << "invalid collection: TriggerSummary" << "\n";
         return;
     }
     
     
     //get online objects
     std::vector<float> ptElectron, etaElectron, phiElectron;
     size_t filterIndex = triggerSummary->filterIndex( triggerFilter_ );
     trigger::TriggerObjectCollection triggerObjects = triggerSummary->getObjects();
     if( !(filterIndex >= triggerSummary->sizeFilters()) ){
         const trigger::Keys& keys = triggerSummary->filterKeys( filterIndex );
         for( size_t j = 0; j < keys.size(); ++j ){
             trigger::TriggerObject foundObject = triggerObjects[keys[j]];
             if(fabs(foundObject.id()) == 11){ //It's an electron
             
                 bool same= false;
                 for(unsigned int x=0;x<ptElectron.size();x++){
                     if(fabs(ptElectron[x] - foundObject.pt()) < 0.01 || fabs(etaElectron[x] - foundObject.eta()) < 0.001 || fabs(phiElectron[x] - foundObject.phi()) < 0.001)
                         same = true;
                 }
                 
                 if(!same){
                     h_triggerElePt->Fill(foundObject.pt());
                     h_triggerEleEta->Fill(foundObject.eta());
                     h_triggerElePhi->Fill(foundObject.phi());
                     ptElectron.push_back(foundObject.pt());
                     etaElectron.push_back(foundObject.eta());
                     phiElectron.push_back(foundObject.phi());
                 }
             }
         }
         if (ptElectron.size()>=2) {
             math::PtEtaPhiMLorentzVectorD* ele1 = new math::PtEtaPhiMLorentzVectorD(ptElectron[0],etaElectron[0],phiElectron[0],0.0005);
             math::PtEtaPhiMLorentzVectorD* ele2 = new math::PtEtaPhiMLorentzVectorD(ptElectron[1],etaElectron[1],phiElectron[1],0.0005);
             (*ele1)+=(*ele2);
             h_triggerDoubleEleMass->Fill(ele1->M());
             delete ele1;
             delete ele2;
         } else {
             h_triggerDoubleEleMass->Fill(-1.);
         }
     }
     
     
     bool hasFired = false;
     bool hasFiredAuxiliaryForElectronLeg = false;
     bool hasFiredAuxiliaryForHadronicLeg = false;
     const edm::TriggerNames& trigNames = e.triggerNames(*hltresults);
     unsigned int numTriggers = trigNames.size();
     for( unsigned int hltIndex=0; hltIndex<numTriggers; ++hltIndex ){
         if (trigNames.triggerName(hltIndex).find(triggerPath_) != std::string::npos && hltresults->wasrun(hltIndex) && hltresults->accept(hltIndex)) hasFired = true;
         if (trigNames.triggerName(hltIndex).find(triggerPathAuxiliaryForElectron_) != std::string::npos && hltresults->wasrun(hltIndex) && hltresults->accept(hltIndex)) hasFiredAuxiliaryForElectronLeg = true;
         if (trigNames.triggerName(hltIndex).find(triggerPathAuxiliaryForHadronic_) != std::string::npos && hltresults->wasrun(hltIndex) && hltresults->accept(hltIndex)) hasFiredAuxiliaryForHadronicLeg = true;
     }
     
     
     
     if(hasFiredAuxiliaryForElectronLeg || hasFiredAuxiliaryForHadronicLeg) {
         
         //Matching the Electron
         int indexOfMatchedElectron[2] = {-1};
         int matchedCounter = 0;
         int offlineCounter = 0;
         for(reco::GsfElectronCollection::const_iterator Electron = ElectronCollection->begin(); (Electron != ElectronCollection->end() && matchedCounter < 2) ; ++Electron) {
             for(size_t off_i = 0; off_i < ptElectron.size(); ++off_i) {
                 if(sqrt((Electron->phi()-phiElectron[off_i])*(Electron->phi()-phiElectron[off_i]) + (Electron->eta()-etaElectron[off_i])*(Electron->eta()-etaElectron[off_i])) < 0.5) {
                     indexOfMatchedElectron[matchedCounter] = offlineCounter;
                     matchedCounter++;
                     break;
                 }
             }
             offlineCounter++;
         }
         
         float caloHT = 0.0;
         float pfHT = 0.0;
         for (reco::PFJetCollection::const_iterator i_pfjet = pfJetCollection->begin(); i_pfjet != pfJetCollection->end(); ++i_pfjet){
             if (i_pfjet->pt() < ptThrJet_) continue;
             if (fabs(i_pfjet->eta()) > etaThrJet_) continue;
             pfHT += i_pfjet->pt();
         }
         for (reco::CaloJetCollection::const_iterator i_calojet = caloJetCollection->begin(); i_calojet != caloJetCollection->end(); ++i_calojet){
             if (i_calojet->pt() < ptThrJet_) continue;
             if (fabs(i_calojet->eta()) > etaThrJet_) continue;
             caloHT += i_calojet->pt();
         }
         
         if(hasFiredAuxiliaryForElectronLeg && ElectronCollection->size()>1) {
             if(hasFired && indexOfMatchedElectron[1] >= 0) { //fill trailing leg
                 h_EleTurnOn_num-> Fill(ElectronCollection->at(indexOfMatchedElectron[1]).pt());
             }
             h_EleTurnOn_den-> Fill(ElectronCollection->at(1).pt());
         }
         if(hasFiredAuxiliaryForHadronicLeg) {
             if(hasFired) {
                 h_pfHTTurnOn_num-> Fill(pfHT);
             }
             h_pfHTTurnOn_den-> Fill(pfHT);
         }
     }
 }

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

overrideprotected

Definition at line 61 of file SUSY_HLT_DoubleEle_Hadronic.cc.

References bookHistos().

 {
     edm::LogInfo("SUSY_HLT_DoubleEle_Hadronic") << "SUSY_HLT_DoubleEle_Hadronic::bookHistograms" << std::endl;
     //book at beginRun
     bookHistos(ibooker_);
 }

void SUSY_HLT_DoubleEle_Hadronic::bookHistos ( DQMStore::IBooker & ibooker_ )

private

Definition at line 223 of file SUSY_HLT_DoubleEle_Hadronic.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::cd(), DEFINE_FWK_MODULE, h_EleTurnOn_den, h_EleTurnOn_num, h_pfHTTurnOn_den, h_pfHTTurnOn_num, h_triggerDoubleEleMass, h_triggerEleEta, h_triggerElePhi, h_triggerElePt, DQMStore::IBooker::setCurrentFolder(), and triggerPath_.

Referenced by bookHistograms().

 {
     ibooker_.cd();
     ibooker_.setCurrentFolder("HLT/SUSYBSM/" + triggerPath_);
     
     //offline quantities
     
     //online quantities
     h_triggerElePt = ibooker_.book1D("triggerElePt", "Trigger Electron Pt; GeV", 50, 0.0, 500.0);
     h_triggerEleEta = ibooker_.book1D("triggerEleEta", "Trigger Electron Eta", 20, -3.0, 3.0);
     h_triggerElePhi = ibooker_.book1D("triggerElePhi", "Trigger Electron Phi", 20, -3.5, 3.5);
 
     h_triggerDoubleEleMass = ibooker_.book1D("triggerDoubleEleMass", "Trigger DoubleElectron Mass", 202, -2, 200);
 
     //num and den hists to be divided in harvesting step to make turn on curves
     h_pfHTTurnOn_num = ibooker_.book1D("pfHTTurnOn_num", "PF HT Turn On Numerator", 30, 0.0, 1500.0 );
     h_pfHTTurnOn_den = ibooker_.book1D("pfHTTurnOn_den", "PF HT Turn On Denominator", 30, 0.0, 1500.0 );
     h_EleTurnOn_num = ibooker_.book1D("EleTurnOn_num", "Electron Turn On Numerator", 30, 0.0, 150 );
     h_EleTurnOn_den = ibooker_.book1D("EleTurnOn_den", "Electron Turn On Denominator", 30, 0.0, 150.0 );
     
     ibooker_.cd();
 }

void SUSY_HLT_DoubleEle_Hadronic::dqmBeginRun	(	edm::Run const &	run,
		edm::EventSetup const &	e
	)

overrideprotected

Definition at line 32 of file SUSY_HLT_DoubleEle_Hadronic.cc.

References fHltConfig, spr::find(), HLTProcess_, HLTConfigProvider::init(), LogDebug, HLTConfigProvider::triggerNames(), and triggerPath_.

 {
     
     bool changed;
     
     if (!fHltConfig.init(run, e, HLTProcess_, changed)) {
         edm::LogError("SUSY_HLT_DoubleEle_Hadronic") << "Initialization of HLTConfigProvider failed!!";
         return;
     }
     
     bool pathFound = false;
     const std::vector<std::string> allTrigNames = fHltConfig.triggerNames();
     for(size_t j = 0; j <allTrigNames.size(); ++j) {
         if(allTrigNames[j].find(triggerPath_) != std::string::npos) {
             pathFound = true;
         }
     }
     
     if(!pathFound) {
         LogDebug("SUSY_HLT_DoubleEle_Hadronic") << "Path not found" << "\n";
         return;
     }
     //std::vector<std::string> filtertags = fHltConfig.moduleLabels( triggerPath_ );
     //triggerFilter_ = edm::InputTag(filtertags[filtertags.size()-1],"",fHltConfig.processName());
     //triggerFilter_ = edm::InputTag("hltPFMET120Mu5L3PreFiltered", "", fHltConfig.processName());
     
     edm::LogInfo("SUSY_HLT_DoubleEle_Hadronic") << "SUSY_HLT_DoubleEle_Hadronic::beginRun" << std::endl;
 }

void SUSY_HLT_DoubleEle_Hadronic::endRun	(	edm::Run const &	run,
		edm::EventSetup const &	eSetup
	)

overrideprotected

Definition at line 218 of file SUSY_HLT_DoubleEle_Hadronic.cc.

 {
     edm::LogInfo("SUSY_HLT_DoubleEle_Hadronic") << "SUSY_HLT_DoubleEle_Hadronic::endRun" << std::endl;
 }

Member Data Documentation

double SUSY_HLT_DoubleEle_Hadronic::etaThrJet_

private

Definition at line 67 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by analyze(), and SUSY_HLT_DoubleEle_Hadronic().

HLTConfigProvider SUSY_HLT_DoubleEle_Hadronic::fHltConfig

private

Definition at line 59 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by dqmBeginRun().

MonitorElement* SUSY_HLT_DoubleEle_Hadronic::h_EleTurnOn_den

private

Definition at line 77 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by analyze(), and bookHistos().

MonitorElement* SUSY_HLT_DoubleEle_Hadronic::h_EleTurnOn_num

private

Definition at line 76 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by analyze(), and bookHistos().

MonitorElement* SUSY_HLT_DoubleEle_Hadronic::h_pfHTTurnOn_den

private

Definition at line 75 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by analyze(), and bookHistos().

MonitorElement* SUSY_HLT_DoubleEle_Hadronic::h_pfHTTurnOn_num

private

Definition at line 74 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by analyze(), and bookHistos().

MonitorElement* SUSY_HLT_DoubleEle_Hadronic::h_triggerDoubleEleMass

private

Definition at line 73 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by analyze(), and bookHistos().

MonitorElement* SUSY_HLT_DoubleEle_Hadronic::h_triggerEleEta

private

Definition at line 71 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by analyze(), and bookHistos().

MonitorElement* SUSY_HLT_DoubleEle_Hadronic::h_triggerElePhi

private

Definition at line 72 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by analyze(), and bookHistos().

MonitorElement* SUSY_HLT_DoubleEle_Hadronic::h_triggerElePt

private

Definition at line 70 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by analyze(), and bookHistos().

std::string SUSY_HLT_DoubleEle_Hadronic::HLTProcess_

private

Definition at line 61 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by dqmBeginRun(), and SUSY_HLT_DoubleEle_Hadronic().

double SUSY_HLT_DoubleEle_Hadronic::ptThrJet_

private

Definition at line 66 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by analyze(), and SUSY_HLT_DoubleEle_Hadronic().

edm::EDGetTokenT<reco::CaloJetCollection> SUSY_HLT_DoubleEle_Hadronic::theCaloJetCollection_

private

Definition at line 54 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by analyze(), and SUSY_HLT_DoubleEle_Hadronic().

edm::EDGetTokenT<reco::GsfElectronCollection> SUSY_HLT_DoubleEle_Hadronic::theElectronCollection_

private

Definition at line 52 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by analyze(), and SUSY_HLT_DoubleEle_Hadronic().

edm::EDGetTokenT<reco::PFJetCollection> SUSY_HLT_DoubleEle_Hadronic::thePfJetCollection_

private

Definition at line 53 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by analyze(), and SUSY_HLT_DoubleEle_Hadronic().

edm::EDGetTokenT<trigger::TriggerEvent> SUSY_HLT_DoubleEle_Hadronic::theTrigSummary_

private

Definition at line 56 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by analyze(), and SUSY_HLT_DoubleEle_Hadronic().

edm::InputTag SUSY_HLT_DoubleEle_Hadronic::triggerFilter_

private

Definition at line 65 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by analyze(), and SUSY_HLT_DoubleEle_Hadronic().

std::string SUSY_HLT_DoubleEle_Hadronic::triggerPath_

private

Definition at line 62 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by analyze(), bookHistos(), dqmBeginRun(), and SUSY_HLT_DoubleEle_Hadronic().

std::string SUSY_HLT_DoubleEle_Hadronic::triggerPathAuxiliaryForElectron_

private

Definition at line 63 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by analyze(), and SUSY_HLT_DoubleEle_Hadronic().

std::string SUSY_HLT_DoubleEle_Hadronic::triggerPathAuxiliaryForHadronic_

private

Definition at line 64 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by analyze(), and SUSY_HLT_DoubleEle_Hadronic().

edm::EDGetTokenT<edm::TriggerResults> SUSY_HLT_DoubleEle_Hadronic::triggerResults_

private

Definition at line 55 of file SUSY_HLT_DoubleEle_Hadronic.h.

Referenced by analyze(), and SUSY_HLT_DoubleEle_Hadronic().

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation