PatBJetTagAnalyzer Class Reference

Inheritance diagram for PatBJetTagAnalyzer:

Classes
struct	Plots

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &es) override
void	beginJob () override
	PatBJetTagAnalyzer (const edm::ParameterSet &params)
	constructor and destructor More...
	~PatBJetTagAnalyzer () override
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
	~EDAnalyzer () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Types
enum	Flavour {   ALL_JETS = 0, UDSG_JETS, C_JETS, B_JETS,   NONID_JETS, N_JET_TYPES }

Private Attributes
TH1 *	flavours_
double	jetEtaCut_
double	jetPtCut_
edm::EDGetTokenT< pat::JetCollection >	jetsToken_
struct PatBJetTagAnalyzer::Plots	plots_ [N_JET_TYPES]

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
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

Definition at line 20 of file PatBJetTagAnalyzer.cc.

Member Enumeration Documentation

enum PatBJetTagAnalyzer::Flavour

private

Enumerator
ALL_JETS
UDSG_JETS
C_JETS
B_JETS
NONID_JETS
N_JET_TYPES

Definition at line 39 of file PatBJetTagAnalyzer.cc.

                 {
        ALL_JETS = 0,
        UDSG_JETS,
        C_JETS,
        B_JETS,
        NONID_JETS,
        N_JET_TYPES
    };

Constructor & Destructor Documentation

PatBJetTagAnalyzer::PatBJetTagAnalyzer

(

const edm::ParameterSet &

params

)

constructor and destructor

Definition at line 56 of file PatBJetTagAnalyzer.cc.

                                                                    :
    jetsToken_(consumes<pat::JetCollection>(params.getParameter<edm::InputTag>("jets"))),
    jetPtCut_(params.getParameter<double>("jetPtCut")),
    jetEtaCut_(params.getParameter<double>("jetEtaCut"))
{
}

PatBJetTagAnalyzer::~PatBJetTagAnalyzer ( )

override

Definition at line 63 of file PatBJetTagAnalyzer.cc.

{
}

Member Function Documentation

void PatBJetTagAnalyzer::analyze ( const edm::Event &  event, const edm::EventSetup &  es  )

override

Definition at line 115 of file PatBJetTagAnalyzer.cc.

References funct::abs(), ALL_JETS, B_JETS, C_JETS, DEFINE_FWK_MODULE, PatBJetTagAnalyzer::Plots::discrCSV, PatBJetTagAnalyzer::Plots::discrSSV, PatBJetTagAnalyzer::Plots::discrTC, GenHFHadronMatcher_cff::flavour, flavours_, metsig::jet, jetEtaCut_, jetPtCut_, jetsToken_, NONID_JETS, plots_, and UDSG_JETS.

{
    // handle to the jets collection
    edm::Handle<pat::JetCollection> jetsHandle;
    event.getByToken(jetsToken_, jetsHandle);

    // now go through all jets
    for(pat::JetCollection::const_iterator jet = jetsHandle->begin();
        jet != jetsHandle->end(); ++jet) {

        // only look at jets that pass the pt and eta cut
        if (jet->pt() < jetPtCut_ ||
            std::abs(jet->eta()) > jetEtaCut_)
            continue;

        Flavour flavour;
        // find out the jet flavour (differs between quark and anti-quark)
        switch(std::abs(jet->partonFlavour())) {
            case 1:
            case 2:
            case 3:
            case 21:
            flavour = UDSG_JETS;
            break;
            case 4:
            flavour = C_JETS;
            break;
            case 5:
            flavour = B_JETS;
            break;
            default:
            flavour = NONID_JETS;
        }

        // simply count the number of accepted jets
        flavours_->Fill(ALL_JETS);
        flavours_->Fill(flavour);

        double discrTC = jet->bDiscriminator("trackCountingHighEffBJetTags");
        double discrSSV = jet->bDiscriminator("simpleSecondaryVertexBJetTags");
        double discrCSV = jet->bDiscriminator("combinedSecondaryVertexBJetTags");

        plots_[ALL_JETS].discrTC->Fill(discrTC);
        plots_[flavour].discrTC->Fill(discrTC);

        plots_[ALL_JETS].discrSSV->Fill(discrSSV);
        plots_[flavour].discrSSV->Fill(discrSSV);

        plots_[ALL_JETS].discrCSV->Fill(discrCSV);
        plots_[flavour].discrCSV->Fill(discrCSV);
    }
}

void PatBJetTagAnalyzer::beginJob ( void  )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 67 of file PatBJetTagAnalyzer.cc.

References ALL_JETS, B_JETS, C_JETS, PatBJetTagAnalyzer::Plots::discrCSV, PatBJetTagAnalyzer::Plots::discrSSV, PatBJetTagAnalyzer::Plots::discrTC, GenHFHadronMatcher_cff::flavour, flavours_, mps_fire::i, TFileService::make(), N_JET_TYPES, dataset::name, nanoDQM_cff::plots, plots_, and UDSG_JETS.

{
    // retrieve handle to auxiliary service
    //  used for storing histograms into ROOT file
    edm::Service<TFileService> fs;

    flavours_ = fs->make<TH1F>("flavours", "jet flavours", 5, 0, 5);

    // book histograms for all jet flavours
    for(unsigned int i = 0; i < N_JET_TYPES; i++) {
        Plots &plots = plots_[i];
        const char *flavour, *name;

        switch((Flavour)i) {
            case ALL_JETS:
            flavour = "all jets";
            name = "all";
            break;
            case UDSG_JETS:
            flavour = "light flavour jets";
            name = "udsg";
            break;
            case C_JETS:
            flavour = "charm jets";
            name = "c";
            break;
            case B_JETS:
            flavour = "bottom jets";
            name = "b";
            break;
            default:
            flavour = "unidentified jets";
            name = "ni";
            break;
        }

        plots.discrTC = fs->make<TH1F>(Form("discrTC_%s", name),
                                       Form("track counting (\"high efficiency\") in %s", flavour),
                                       100, 0, 20);
        plots.discrSSV = fs->make<TH1F>(Form("discrSSV_%s", name),
                                        Form("simple secondary vertex in %s", flavour),
                                        100, 0, 10);
        plots.discrCSV = fs->make<TH1F>(Form("discrCSV_%s", name),
                                        Form("combined secondary vertex in %s", flavour),
                                        100, 0, 1);
    }
}

Member Data Documentation

TH1* PatBJetTagAnalyzer::flavours_

private

Definition at line 48 of file PatBJetTagAnalyzer.cc.

Referenced by analyze(), and beginJob().

double PatBJetTagAnalyzer::jetEtaCut_

private

Definition at line 35 of file PatBJetTagAnalyzer.cc.

Referenced by analyze().

double PatBJetTagAnalyzer::jetPtCut_

private

Definition at line 34 of file PatBJetTagAnalyzer.cc.

Referenced by analyze().

edm::EDGetTokenT<pat::JetCollection> PatBJetTagAnalyzer::jetsToken_

private

Definition at line 32 of file PatBJetTagAnalyzer.cc.

Referenced by analyze().

struct PatBJetTagAnalyzer::Plots PatBJetTagAnalyzer::plots_[N_JET_TYPES]

private

Referenced by analyze(), and beginJob().