#include <HLTTauDQMPathPlotter.h>

Inheritance diagram for HLTTauDQMPathPlotter:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &, const std::map< int, LVColl > &)
	HLTTauDQMPathPlotter (const edm::ParameterSet &, bool, std::string)
const std::string	name ()
	~HLTTauDQMPathPlotter ()
Private Member Functions
void	endJob ()
LVColl	getFilterCollection (size_t, int, const trigger::TriggerEventWithRefs &)
Private Attributes
MonitorElement *	accepted_events
MonitorElement *	accepted_events_matched
bool	doRefAnalysis_
std::vector < HLTTauDQMPlotter::FilterObject >	filterObjs_
std::vector< edm::ParameterSet >	filters_
edm::ParameterSet	reference_
double	refLeptonPt_
unsigned int	refNTriggeredLeptons_
unsigned int	refNTriggeredTaus_
double	refTauPt_
edm::InputTag	triggerEventObject_

Detailed Description

Definition at line 25 of file HLTTauDQMPathPlotter.h.

Constructor & Destructor Documentation

HLTTauDQMPathPlotter::HLTTauDQMPathPlotter	(	const edm::ParameterSet &	ps,
		bool	ref,
		std::string	dqmBaseFolder
	)

Definition at line 3 of file HLTTauDQMPathPlotter.cc.

References accepted_events, accepted_events_matched, DQMStore::book1D(), doRefAnalysis_, HLTTauDQMPlotter::dqmBaseFolder_, alignCSCRings::e, filterObjs_, filters_, edm::ParameterSet::getUntrackedParameter(), HLTTauDQMPlotter::FilterObject::isValid(), gen::k, HLTTauDQMPlotter::name_, reference_, refLeptonPt_, refNTriggeredLeptons_, refNTriggeredTaus_, refTauPt_, DQMStore::removeContents(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), HLTTauDQMPlotter::store_, tmp, triggerEventObject_, HLTTauDQMPlotter::triggerTag(), HLTTauDQMPlotter::triggerTag_, HLTTauDQMPlotter::triggerTagAlias_, HLTTauDQMPlotter::validity_, and cms::Exception::what().

                                                                                                         {
    //Initialize Plotter
    name_ = "HLTTauDQMPathPlotter";
    
    //Process PSet
    try {
        triggerEventObject_   = ps.getUntrackedParameter<edm::InputTag>("TriggerEventObject");
        triggerTag_           = ps.getUntrackedParameter<std::string>("DQMFolder");
        triggerTagAlias_      = ps.getUntrackedParameter<std::string>("Alias","");
        filters_              = ps.getUntrackedParameter<std::vector<edm::ParameterSet> >("Filters");
        reference_            = ps.getUntrackedParameter<edm::ParameterSet>("Reference");
        refNTriggeredTaus_    = reference_.getUntrackedParameter<unsigned int>("NTriggeredTaus");
        refNTriggeredLeptons_ = reference_.getUntrackedParameter<unsigned int>("NTriggeredLeptons");
        refTauPt_             = reference_.getUntrackedParameter<double>("refTauPt",20);
        refLeptonPt_          = reference_.getUntrackedParameter<double>("refLeptonPt",15);
        dqmBaseFolder_        = dqmBaseFolder;
        doRefAnalysis_        = ref;
        validity_             = true;
    } catch ( cms::Exception &e ) {
        edm::LogInfo("HLTTauDQMPathPlotter::HLTTauDQMPathPlotter") << e.what() << std::endl;
        validity_ = false;
        return;
    }
    
    for ( std::vector<edm::ParameterSet>::const_iterator iter = filters_.begin(); iter != filters_.end(); ++iter ) {
        HLTTauDQMPlotter::FilterObject tmp(*iter);
        if (tmp.isValid()) filterObjs_.push_back(tmp);
    }
    
    if (store_) {
        //Create the histograms
        store_->setCurrentFolder(triggerTag());
        store_->removeContents();
        
        accepted_events = store_->book1D("TriggerBits","Accepted Events per Path;;entries",filterObjs_.size(),0,filterObjs_.size());
        for ( size_t k = 0; k < filterObjs_.size(); ++k ) {
            accepted_events->setBinLabel(k+1,filterObjs_[k].getAlias(),1);
        }
        if (doRefAnalysis_) {
            accepted_events_matched = store_->book1D("MatchedTriggerBits","Accepted+Matched Events per Path;;entries",filterObjs_.size()+1,0,filterObjs_.size()+1);
            accepted_events_matched->setBinLabel(1,"RefEvents",1);
            for ( size_t k = 0; k < filterObjs_.size(); ++k ) {
                accepted_events_matched->setBinLabel(k+2,filterObjs_[k].getAlias(),1);
            }
        }
    }
}

HLTTauDQMPathPlotter::~HLTTauDQMPathPlotter ( )

Definition at line 51 of file HLTTauDQMPathPlotter.cc.

{
}

Member Function Documentation

void HLTTauDQMPathPlotter::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		const std::map< int, LVColl > &	refC
	)

Definition at line 58 of file HLTTauDQMPathPlotter.cc.

References accepted_events, accepted_events_matched, doRefAnalysis_, MonitorElement::Fill(), filterObjs_, edm::Event::getByLabel(), getFilterCollection(), i, j, EgammaValidation_Wenu_cff::leptons, HLTTauDQMPlotter::match(), dt_dqm_sourceclient_common_cff::reco, refLeptonPt_, refNTriggeredLeptons_, refNTriggeredTaus_, refTauPt_, and triggerEventObject_.

                                                                                                                        {
    using namespace std;
    using namespace edm;
    using namespace reco;
    using namespace l1extra;
    using namespace trigger;
    
    bool isGoodReferenceEvent = false;
    LVColl refTaus, refLeptons;
    
    if (doRefAnalysis_) {
        unsigned int highPtTaus = 0;
        unsigned int highPtElectrons = 0;
        unsigned int highPtMuons = 0;
        
        bool tau_ok = true;
        bool leptons_ok = true;
        
        std::map<int,LVColl>::const_iterator iref;
        
        //Tau reference
        iref = refC.find(15);
        if ( iref != refC.end() ) {
            for ( LVColl::const_iterator lvi = iref->second.begin(); lvi != iref->second.end(); ++lvi ) {
                if ( lvi->Et() > refTauPt_ ) {
                    highPtTaus++;
                }
                refTaus.push_back(*lvi);
            }
        }
        if ( highPtTaus < refNTriggeredTaus_ ) tau_ok = false;
        
        //Electron reference
        iref = refC.find(11);
        if ( iref != refC.end() ) {
            for ( LVColl::const_iterator lvi = iref->second.begin(); lvi != iref->second.end(); ++lvi ) {
                if ( lvi->Et() > refLeptonPt_ ) {
                    highPtElectrons++;
                }
                refLeptons.push_back(*lvi);
            }
        }
        if ( filterObjs_.size() > 0 && filterObjs_.back().leptonId() == 11 && highPtElectrons < refNTriggeredLeptons_ ) leptons_ok = false;
        
        //Muon reference
        iref = refC.find(13);
        if ( iref != refC.end() ) {
            for ( LVColl::const_iterator lvi = iref->second.begin(); lvi != iref->second.end(); ++lvi ) {
                if ( lvi->Et() > refLeptonPt_ ) {
                    highPtMuons++;
                }
                refLeptons.push_back(*lvi);
            }
        }
        if ( filterObjs_.size() > 0 && filterObjs_.back().leptonId() == 13 && highPtMuons < refNTriggeredLeptons_ ) leptons_ok = false;
        
        if ( tau_ok && leptons_ok ) {
            accepted_events_matched->Fill(0.5);
            isGoodReferenceEvent = true;
        }
    }
    
    Handle<TriggerEventWithRefs> trigEv;
    bool gotTEV = iEvent.getByLabel(triggerEventObject_,trigEv) && trigEv.isValid();
    
    if (gotTEV) {
        //Loop through the filters
        for ( size_t i = 0; i < filterObjs_.size(); ++i ) {             
            size_t ID = trigEv->filterIndex(filterObjs_[i].getFilterName());
            
            if ( ID != trigEv->size() ) {
                LVColl leptons = getFilterCollection(ID,filterObjs_[i].getLeptonType(),*trigEv);
                LVColl taus = getFilterCollection(ID,filterObjs_[i].getTauType(),*trigEv);
                //Fired
                if ( leptons.size() >= filterObjs_[i].getNTriggeredLeptons() && taus.size() >= filterObjs_[i].getNTriggeredTaus() ) {
                    accepted_events->Fill(i+0.5);
                    //Now do the matching only though if we have a good reference event
                    if ( doRefAnalysis_ && isGoodReferenceEvent ) {
                        size_t nT = 0;
                        for ( size_t j = 0; j < taus.size(); ++j ) {
                            if( match(taus[j],refTaus,filterObjs_[i].getMatchDeltaR()).first ) nT++;
                        }
                        
                        size_t nL = 0;
                        for ( size_t j = 0; j < leptons.size(); ++j ) {
                            if ( match(leptons[j],refLeptons,filterObjs_[i].getMatchDeltaR()).first ) nL++;
                        }
                        
                        if ( nT >= filterObjs_[i].getNTriggeredTaus() && nL >= filterObjs_[i].getNTriggeredLeptons() ) {
                            accepted_events_matched->Fill(i+1.5);
                        }
                    }
                }
            }
        }
    }
}

void HLTTauDQMPathPlotter::endJob ( ) [private]

LVColl HLTTauDQMPathPlotter::getFilterCollection	(	size_t	filterID,
		int	id,
		const trigger::TriggerEventWithRefs &	trigEv
	)		`[private]`

Definition at line 156 of file HLTTauDQMPathPlotter.cc.

References trigger::TriggerEventWithRefs::getObjects(), i, VarParsing::obj, dbtoconf::out, p4, trigger::TriggerElectron, trigger::TriggerL1CenJet, trigger::TriggerL1IsoEG, trigger::TriggerL1Mu, trigger::TriggerL1NoIsoEG, trigger::TriggerL1TauJet, trigger::TriggerMuon, and trigger::TriggerTau.

Referenced by analyze().

                                                                                                                     {
    using namespace trigger;
    LVColl out;
    
    if ( id == trigger::TriggerL1IsoEG || id == trigger::TriggerL1NoIsoEG ) {
        VRl1em obj;
        trigEv.getObjects(filterID,id,obj);
        for (size_t i=0;i<obj.size();++i)
            if (obj.at(i).isAvailable())
                out.push_back(obj[i]->p4());
    }
    
    if ( id == trigger::TriggerL1Mu ) {
        VRl1muon obj;
        trigEv.getObjects(filterID,id,obj);
        for (size_t i=0;i<obj.size();++i)
            if (obj.at(i).isAvailable())
                out.push_back(obj[i]->p4());
    }
    
    if ( id == trigger::TriggerMuon ) {
        VRmuon obj;
        trigEv.getObjects(filterID,id,obj);
        for (size_t i=0;i<obj.size();++i)
            if (obj.at(i).isAvailable())
                out.push_back(obj[i]->p4());
    }
    
    if ( id == trigger::TriggerElectron ) {
        VRelectron obj;
        trigEv.getObjects(filterID,id,obj);
        for (size_t i=0;i<obj.size();++i)
            if (obj.at(i).isAvailable())
                out.push_back(obj[i]->p4());
    }
    
    if ( id == trigger::TriggerL1TauJet ) {
        VRl1jet obj;
        trigEv.getObjects(filterID,id,obj);
        for (size_t i=0;i<obj.size();++i)
            if (obj.at(i).isAvailable())
                out.push_back(obj[i]->p4());
        trigEv.getObjects(filterID,trigger::TriggerL1CenJet,obj);
        for (size_t i=0;i<obj.size();++i)
            if (obj.at(i).isAvailable())
                out.push_back(obj[i]->p4());
        
    }
    
    if ( id == trigger::TriggerTau ) {
        VRjet obj;
        trigEv.getObjects(filterID,id,obj);
        for (size_t i = 0; i < obj.size(); ++i) {
            if (obj.at(i).isAvailable()) {
                out.push_back(obj[i]->p4());
            }
        }
        VRpfjet pfjetobj;
        trigEv.getObjects(filterID,id,pfjetobj);
        for (size_t i = 0; i < pfjetobj.size(); ++i) {
            if (pfjetobj.at(i).isAvailable()) {
                out.push_back(pfjetobj[i]->p4());
            }
        }
        VRpftau pftauobj;
        trigEv.getObjects(filterID,id,pftauobj);
        for (size_t i = 0; i < pftauobj.size(); ++i) {
            if (pftauobj.at(i).isAvailable()) {
                out.push_back(pftauobj[i]->p4());
            }
        }
    }

    return out;
}