EGammaCutBasedEleIdAnalyzer Class Reference

#include <EGamma/EGammaCutBasedEleIdAnalyzer/src/EGammaCutBasedEleIdAnalyzer.cc>

Inheritance diagram for EGammaCutBasedEleIdAnalyzer:

Public Types
typedef std::vector < edm::Handle< edm::ValueMap < reco::IsoDeposit > > >	IsoDepositMaps
typedef std::vector < edm::Handle< edm::ValueMap < double > > >	IsoDepositVals
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Member Functions
	EGammaCutBasedEleIdAnalyzer (const edm::ParameterSet &)
	~EGammaCutBasedEleIdAnalyzer ()
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
virtual void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	beginRun (edm::Run const &, edm::EventSetup const &)
virtual void	endJob ()
virtual void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	endRun (edm::Run const &, edm::EventSetup const &)

Private Attributes
edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken_
edm::EDGetTokenT < reco::ConversionCollection >	conversionsToken_
edm::EDGetTokenT < reco::GsfElectronCollection >	electronsToken_
TH1F *	h1_pt_
TH1F *	h1_pt_fbremeopin_
TH1F *	h1_pt_loose_
TH1F *	h1_pt_medium_
TH1F *	h1_pt_tight_
TH1F *	h1_pt_trig_
TH1F *	h1_pt_veto_
std::vector< edm::EDGetTokenT < edm::ValueMap< double > > >	isoValTokens_
edm::EDGetTokenT < reco::VertexCollection >	primaryVertexToken_
bool	printDebug_
edm::EDGetTokenT< double >	rhoIsoToken_

Additional Inherited Members
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 48 of file EGammaCutBasedEleIdAnalyzer.cc.

Member Typedef Documentation

typedef std::vector< edm::Handle< edm::ValueMap<reco::IsoDeposit> > > EGammaCutBasedEleIdAnalyzer::IsoDepositMaps

Definition at line 51 of file EGammaCutBasedEleIdAnalyzer.cc.

typedef std::vector< edm::Handle< edm::ValueMap<double> > > EGammaCutBasedEleIdAnalyzer::IsoDepositVals

Definition at line 52 of file EGammaCutBasedEleIdAnalyzer.cc.

Constructor & Destructor Documentation

EGammaCutBasedEleIdAnalyzer::EGammaCutBasedEleIdAnalyzer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 102 of file EGammaCutBasedEleIdAnalyzer.cc.

References beamSpotToken_, edm::EDConsumerBase::consumes(), conversionsToken_, electronsToken_, edm::ParameterSet::getParameter(), isoValTokens_, primaryVertexToken_, rhoIsoToken_, GlobalPosition_Frontier_DevDB_cff::tag, and edm::vector_transform().

{

    // get input parameters
    electronsToken_      = consumes<reco::GsfElectronCollection>(iConfig.getParameter<edm::InputTag>("electronsInputTag"));
    conversionsToken_    = consumes<reco::ConversionCollection>(iConfig.getParameter<edm::InputTag>("conversionsInputTag"));
    beamSpotToken_       = consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("beamSpotInputTag"));
    rhoIsoToken_         = consumes<double>(iConfig.getParameter<edm::InputTag>("rhoIsoInputTag"));
    primaryVertexToken_  = consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("primaryVertexInputTag"));
    isoValTokens_        = edm::vector_transform(iConfig.getParameter<std::vector<edm::InputTag> >("isoValInputTags"), [this](edm::InputTag const & tag){return consumes<edm::ValueMap<double> >(tag);});

    // debug
    printDebug_             = iConfig.getParameter<bool>("printDebug");

    // output histograms
    edm::Service<TFileService> fs;

    h1_pt_               = fs->make<TH1F>("h1_pt",               "pt",              100, 0.0, 100.0);
    h1_pt_veto_          = fs->make<TH1F>("h1_pt_veto",          "pt (veto)",       100, 0.0, 100.0);
    h1_pt_loose_         = fs->make<TH1F>("h1_pt_loose",         "pt (loose)",      100, 0.0, 100.0);
    h1_pt_medium_        = fs->make<TH1F>("h1_pt_medium",        "pt (medium)",     100, 0.0, 100.0);
    h1_pt_tight_         = fs->make<TH1F>("h1_pt_tight",         "pt (tight)",      100, 0.0, 100.0);
    h1_pt_trig_          = fs->make<TH1F>("h1_pt_trig",          "pt (trig)",       100, 0.0, 100.0);
    h1_pt_fbremeopin_    = fs->make<TH1F>("h1_pt_fbremeopin",    "pt (fbremeopin)", 100, 0.0, 100.0);

}

EGammaCutBasedEleIdAnalyzer::~EGammaCutBasedEleIdAnalyzer

(

)

Definition at line 130 of file EGammaCutBasedEleIdAnalyzer.cc.

{

    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)

}

Member Function Documentation

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 145 of file EGammaCutBasedEleIdAnalyzer.cc.

References SiPixelRawToDigiRegional_cfi::beamSpot, beamSpotToken_, conversionsToken_, electronsToken_, edm::EventID::event(), edm::Event::getByToken(), h1_pt_, h1_pt_fbremeopin_, h1_pt_loose_, h1_pt_medium_, h1_pt_tight_, h1_pt_trig_, h1_pt_veto_, i, edm::EventBase::id(), isoValTokens_, j, EgammaCutBasedEleId::LOOSE, edm::EventBase::luminosityBlock(), EgammaCutBasedEleId::MEDIUM, n, EgammaCutBasedEleId::PassEoverPCuts(), EgammaCutBasedEleId::PassTriggerCuts(), EgammaCutBasedEleId::PassWP(), primaryVertexToken_, printDebug_, edm::Handle< T >::product(), rhoIsoToken_, edm::EventID::run(), EgammaCutBasedEleId::TIGHT, EgammaCutBasedEleId::TRIGGERTIGHT, EgammaCutBasedEleId::TRIGGERWP70, and EgammaCutBasedEleId::VETO.

{
    // electrons
    edm::Handle<reco::GsfElectronCollection> els_h;
    iEvent.getByToken(electronsToken_, els_h);

    // conversions
    edm::Handle<reco::ConversionCollection> conversions_h;
    iEvent.getByToken(conversionsToken_, conversions_h);

    // iso deposits
    IsoDepositVals isoVals(isoValTokens_.size());
    for (size_t j = 0; j < isoValTokens_.size(); ++j) {
        iEvent.getByToken(isoValTokens_[j], isoVals[j]);
    }

    // beam spot
    edm::Handle<reco::BeamSpot> beamspot_h;
    iEvent.getByToken(beamSpotToken_, beamspot_h);
    const reco::BeamSpot &beamSpot = *(beamspot_h.product());

    // vertices
    edm::Handle<reco::VertexCollection> vtx_h;
    iEvent.getByToken(primaryVertexToken_, vtx_h);

    // rho for isolation
    edm::Handle<double> rhoIso_h;
    iEvent.getByToken(rhoIsoToken_, rhoIso_h);
    double rhoIso = *(rhoIso_h.product());

    // loop on electrons
    unsigned int n = els_h->size();
    for(unsigned int i = 0; i < n; ++i) {

        // get reference to electron
        reco::GsfElectronRef ele(els_h, i);

        //
        // get particle flow isolation
        //

        double iso_ch =  (*(isoVals)[0])[ele];
        double iso_em = (*(isoVals)[1])[ele];
        double iso_nh = (*(isoVals)[2])[ele];

        //
        // test ID
        //

        // working points
        bool veto       = EgammaCutBasedEleId::PassWP(EgammaCutBasedEleId::VETO, ele, conversions_h, beamSpot, vtx_h, iso_ch, iso_em, iso_nh, rhoIso);
        bool loose      = EgammaCutBasedEleId::PassWP(EgammaCutBasedEleId::LOOSE, ele, conversions_h, beamSpot, vtx_h, iso_ch, iso_em, iso_nh, rhoIso);
        bool medium     = EgammaCutBasedEleId::PassWP(EgammaCutBasedEleId::MEDIUM, ele, conversions_h, beamSpot, vtx_h, iso_ch, iso_em, iso_nh, rhoIso);
        bool tight      = EgammaCutBasedEleId::PassWP(EgammaCutBasedEleId::TIGHT, ele, conversions_h, beamSpot, vtx_h, iso_ch, iso_em, iso_nh, rhoIso);

        // eop/fbrem cuts for extra tight ID
        bool fbremeopin = EgammaCutBasedEleId::PassEoverPCuts(ele);

        // cuts to match tight trigger requirements
        bool trigtight = EgammaCutBasedEleId::PassTriggerCuts(EgammaCutBasedEleId::TRIGGERTIGHT, ele);

        // for 2011 WP70 trigger
        bool trigwp70 = EgammaCutBasedEleId::PassTriggerCuts(EgammaCutBasedEleId::TRIGGERWP70, ele);

        //
        // fill histograms
        //

        h1_pt_->Fill(ele->pt());
        if (veto)       h1_pt_veto_         ->Fill(ele->pt());
        if (loose)      h1_pt_loose_        ->Fill(ele->pt());
        if (medium)     h1_pt_medium_       ->Fill(ele->pt());
        if (tight)      h1_pt_tight_        ->Fill(ele->pt());
        if (trigtight)  h1_pt_trig_         ->Fill(ele->pt());
        if (fbremeopin) h1_pt_fbremeopin_   ->Fill(ele->pt());

        //
        // print decisions
        //

        if (printDebug_) {
            printf("%u %u %u : ",       iEvent.id().run(), iEvent.luminosityBlock(), iEvent.id().event());
            printf("veto(%i), ",        veto);
            printf("loose(%i), ",       loose);
            printf("medium(%i), ",      medium);
            printf("tight(%i), ",       tight);
            printf("trigtight(%i), ",   trigtight);
            printf("trigwp70(%i), ",    trigwp70);
            printf("fbremeopin(%i)\n",  fbremeopin);
        }

    }

}

void EGammaCutBasedEleIdAnalyzer::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 243 of file EGammaCutBasedEleIdAnalyzer.cc.

{
}

void EGammaCutBasedEleIdAnalyzer::beginLuminosityBlock ( edm::LuminosityBlock const &  , edm::EventSetup const &    )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 267 of file EGammaCutBasedEleIdAnalyzer.cc.

{
}

void EGammaCutBasedEleIdAnalyzer::beginRun ( edm::Run const &  , edm::EventSetup const &    )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 255 of file EGammaCutBasedEleIdAnalyzer.cc.

{
}

void EGammaCutBasedEleIdAnalyzer::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 249 of file EGammaCutBasedEleIdAnalyzer.cc.

{
}

void EGammaCutBasedEleIdAnalyzer::endLuminosityBlock ( edm::LuminosityBlock const &  , edm::EventSetup const &    )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 273 of file EGammaCutBasedEleIdAnalyzer.cc.

{
}

void EGammaCutBasedEleIdAnalyzer::endRun ( edm::Run const &  , edm::EventSetup const &    )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 261 of file EGammaCutBasedEleIdAnalyzer.cc.

{
}

void EGammaCutBasedEleIdAnalyzer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 279 of file EGammaCutBasedEleIdAnalyzer.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

                                                                                        {
    //The following says we do not know what parameters are allowed so do no validation
    // Please change this to state exactly what you do use, even if it is no parameters
    edm::ParameterSetDescription desc;
    desc.setUnknown();
    descriptions.addDefault(desc);
}

Member Data Documentation

edm::EDGetTokenT<reco::BeamSpot> EGammaCutBasedEleIdAnalyzer::beamSpotToken_

private

Definition at line 76 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().

edm::EDGetTokenT<reco::ConversionCollection> EGammaCutBasedEleIdAnalyzer::conversionsToken_

private

Definition at line 75 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().

edm::EDGetTokenT<reco::GsfElectronCollection> EGammaCutBasedEleIdAnalyzer::electronsToken_

private

Definition at line 74 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().

TH1F* EGammaCutBasedEleIdAnalyzer::h1_pt_

private

Definition at line 85 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze().

TH1F* EGammaCutBasedEleIdAnalyzer::h1_pt_fbremeopin_

private

Definition at line 91 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze().

TH1F* EGammaCutBasedEleIdAnalyzer::h1_pt_loose_

private

Definition at line 87 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze().

TH1F* EGammaCutBasedEleIdAnalyzer::h1_pt_medium_

private

Definition at line 88 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze().

TH1F* EGammaCutBasedEleIdAnalyzer::h1_pt_tight_

private

Definition at line 89 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze().

TH1F* EGammaCutBasedEleIdAnalyzer::h1_pt_trig_

private

Definition at line 90 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze().

TH1F* EGammaCutBasedEleIdAnalyzer::h1_pt_veto_

private

Definition at line 86 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze().

std::vector<edm::EDGetTokenT<edm::ValueMap<double> > > EGammaCutBasedEleIdAnalyzer::isoValTokens_

private

Definition at line 79 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().

edm::EDGetTokenT<reco::VertexCollection> EGammaCutBasedEleIdAnalyzer::primaryVertexToken_

private

Definition at line 78 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().

bool EGammaCutBasedEleIdAnalyzer::printDebug_

private

Definition at line 82 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze().

edm::EDGetTokenT<double> EGammaCutBasedEleIdAnalyzer::rhoIsoToken_

private

Definition at line 77 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().