PatBJetTrackAnalyzer Class Reference

Inheritance diagram for PatBJetTrackAnalyzer:

Classes
struct	Plots

Public Member Functions
virtual void	analyze (const edm::Event &event, const edm::EventSetup &es) override
virtual void	beginJob () override
	PatBJetTrackAnalyzer (const edm::ParameterSet &params)
	constructor and destructor More...
	~PatBJetTrackAnalyzer ()
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
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

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

Private Attributes
edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken_
TH1 *	flavours_
double	jetEtaCut_
double	jetPtCut_
edm::EDGetTokenT < pat::JetCollection >	jetsToken_
double	maxDeltaR_
unsigned int	minPixelHits_
double	minPt_
unsigned int	minTotalHits_
unsigned int	nThTrack_
struct PatBJetTrackAnalyzer::Plots	plots_ [N_JET_TYPES]
edm::EDGetTokenT < reco::VertexCollection >	primaryVerticesToken_
edm::EDGetTokenT < reco::TrackCollection >	tracksToken_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
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

Definition at line 30 of file PatBJetTrackAnalyzer.cc.

Member Enumeration Documentation

enum PatBJetTrackAnalyzer::Flavour

private

Enumerator
ALL_JETS
UDSG_JETS
C_JETS
B_JETS
NONID_JETS
N_JET_TYPES

Definition at line 59 of file PatBJetTrackAnalyzer.cc.

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

Constructor & Destructor Documentation

PatBJetTrackAnalyzer::PatBJetTrackAnalyzer

(

const edm::ParameterSet &

params

)

constructor and destructor

Definition at line 79 of file PatBJetTrackAnalyzer.cc.

                                                                        :
    jetsToken_(consumes<pat::JetCollection>(params.getParameter<edm::InputTag>("jets"))),
    tracksToken_(consumes<reco::TrackCollection>(params.getParameter<edm::InputTag>("tracks"))),
    beamSpotToken_(consumes<reco::BeamSpot>(params.getParameter<edm::InputTag>("beamSpot"))),
    primaryVerticesToken_(consumes<reco::VertexCollection>(params.getParameter<edm::InputTag>("primaryVertices"))),
    jetPtCut_(params.getParameter<double>("jetPtCut")),
    jetEtaCut_(params.getParameter<double>("jetEtaCut")),
    maxDeltaR_(params.getParameter<double>("maxDeltaR")),
    minPt_(params.getParameter<double>("minPt")),
    minPixelHits_(params.getParameter<unsigned int>("minPixelHits")),
    minTotalHits_(params.getParameter<unsigned int>("minTotalHits")),
    nThTrack_(params.getParameter<unsigned int>("nThTrack"))
{
}

PatBJetTrackAnalyzer::~PatBJetTrackAnalyzer

(

)

Definition at line 94 of file PatBJetTrackAnalyzer.cc.

{
}

Member Function Documentation

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

overridevirtual

Implements edm::EDAnalyzer.

Definition at line 179 of file PatBJetTrackAnalyzer.cc.

References funct::abs(), ALL_JETS, PatBJetTrackAnalyzer::Plots::allDeltaR, PatBJetTrackAnalyzer::Plots::allIP, PatBJetTrackAnalyzer::Plots::allIPErr, PatBJetTrackAnalyzer::Plots::allIPSig, B_JETS, SiPixelRawToDigiRegional_cfi::beamSpot, beamSpotToken_, C_JETS, deltaR(), Measurement1D::error(), reco::flavour(), flavours_, metsig::jet, jetEtaCut_, jetPtCut_, jetsToken_, maxDeltaR_, minPixelHits_, minPt_, minTotalHits_, PatBJetTrackAnalyzer::Plots::negativeIP, PatBJetTrackAnalyzer::Plots::negativeIPErr, PatBJetTrackAnalyzer::Plots::negativeIPSig, NONID_JETS, nThTrack_, PatBJetTrackAnalyzer::Plots::nTracks, plots_, primaryVerticesToken_, Measurement1D::significance(), significanceHigher(), python.multivaluedict::sort(), PatBJetTrackAnalyzer::Plots::trackIP, PatBJetTrackAnalyzer::Plots::trackIPErr, PatBJetTrackAnalyzer::Plots::trackIPSig, tracksToken_, UDSG_JETS, and Measurement1D::value().

{
    // handle to the primary vertex collection
    edm::Handle<reco::VertexCollection> pvHandle;
    event.getByToken(primaryVerticesToken_, pvHandle);

    // handle to the tracks collection
    edm::Handle<reco::TrackCollection> tracksHandle;
    event.getByToken(tracksToken_, tracksHandle);

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

    // handle to the beam spot
    edm::Handle<reco::BeamSpot> beamSpot;
    event.getByToken(beamSpotToken_, beamSpot);

    // rare case of no reconstructed primary vertex
    if (pvHandle->empty())
        return;

    // extract the position of the (most probable) reconstructed vertex
    math::XYZPoint pv = (*pvHandle)[0].position();

    // 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);

        // this vector will contain IP value / error pairs
        std::vector<Measurement1D> ipValErr;

        // Note: PAT is also able to store associated tracks
        //       within the jet object, so we don't have to do the
        //       matching ourselves
        // (see ->associatedTracks() method)
        // However, using this we can't play with the DeltaR cone
        // withour rerunning the PAT producer

        // now loop through all tracks
        for(reco::TrackCollection::const_iterator track = tracksHandle->begin();
            track != tracksHandle->end(); ++track) {

            // check the quality criteria
            if (track->pt() < minPt_ ||
                track->hitPattern().numberOfValidHits() < (int)minTotalHits_ ||
                track->hitPattern().numberOfValidPixelHits() < (int)minPixelHits_)
                continue;

            // check the Delta R between jet axis and track
            // (Delta_R^2 = Delta_Eta^2 + Delta_Phi^2)
            double deltaR = ROOT::Math::VectorUtil::DeltaR(
                    jet->momentum(), track->momentum());

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

            // only look at tracks in jet cone
            if (deltaR > maxDeltaR_)
                continue;

            // What follows here is an approximation!
            //
            // The dxy() method of the tracks does a linear
            // extrapolation from the track reference position
            // given as the closest point to the beam spot
            // with respect to the given vertex.
            // Since we are using primary vertices, this
            // approximation works well
            //
            // In order to get better results, the
            // "TransientTrack" and specialised methods have
            // to be used.
            // Or look at the "impactParameterTagInfos",
            // which contains the precomputed information
            // from the official b-tagging algorithms
            //
            // see ->tagInfoTrackIP() method

            double ipError = track->dxyError();
            double ipValue = std::abs(track->dxy(pv));

            // in order to compute the sign, we check if
            // the point of closest approach to the vertex
            // is in front or behind the vertex.
            // Again, we a linear approximation
            //
            // dot product between reference point and jet axis

            math::XYZVector closestPoint = track->referencePoint() - beamSpot->position();
            // only interested in transverse component, z -> 0
            closestPoint.SetZ(0.);
            double sign = closestPoint.Dot(jet->momentum());

            if (sign < 0)
                ipValue = -ipValue;

            ipValErr.push_back(Measurement1D(ipValue, ipError));
        }

        // now order all tracks by significance (highest first)
        std::sort(ipValErr.begin(), ipValErr.end(), significanceHigher);

        plots_[ALL_JETS].nTracks->Fill(ipValErr.size());
        plots_[flavour].nTracks->Fill(ipValErr.size());

        // plot all tracks

        for(std::vector<Measurement1D>::const_iterator iter = ipValErr.begin();
            iter != ipValErr.end(); ++iter) {
            plots_[ALL_JETS].allIP->Fill(iter->value());
            plots_[flavour].allIP->Fill(iter->value());

            plots_[ALL_JETS].allIPErr->Fill(iter->error());
            plots_[flavour].allIPErr->Fill(iter->error());

            // significance (is really just value / error)
            plots_[ALL_JETS].allIPSig->Fill(iter->significance());
            plots_[flavour].allIPSig->Fill(iter->significance());
        }

        // check if we have enough tracks to fulfill the
        // n-th track requirement
        if (ipValErr.size() < nThTrack_)
            continue;

        // pick the n-th highest track
        const Measurement1D *trk = &ipValErr[nThTrack_ - 1];

        plots_[ALL_JETS].trackIP->Fill(trk->value());
        plots_[flavour].trackIP->Fill(trk->value());

        plots_[ALL_JETS].trackIPErr->Fill(trk->error());
        plots_[flavour].trackIPErr->Fill(trk->error());

        plots_[ALL_JETS].trackIPSig->Fill(trk->significance());
        plots_[flavour].trackIPSig->Fill(trk->significance());

        // here we define a "negative tagger", i.e. we take
        // the track with the n-lowest signed IP
        // (i.e. preferrably select tracks that appear to become
        //  from "behind" the primary vertex, supposedly mismeasured
        //  tracks really coming from the primary vertex, and
        //  the contamination of displaced tracks should be small)
        trk = &ipValErr[ipValErr.size() - nThTrack_];

        plots_[ALL_JETS].negativeIP->Fill(trk->value());
        plots_[flavour].negativeIP->Fill(trk->value());

        plots_[ALL_JETS].negativeIPErr->Fill(trk->error());
        plots_[flavour].negativeIPErr->Fill(trk->error());

        plots_[ALL_JETS].negativeIPSig->Fill(trk->significance());
        plots_[flavour].negativeIPSig->Fill(trk->significance());
    }
}

void PatBJetTrackAnalyzer::beginJob ( void  )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 98 of file PatBJetTrackAnalyzer.cc.

References ALL_JETS, PatBJetTrackAnalyzer::Plots::allDeltaR, PatBJetTrackAnalyzer::Plots::allIP, PatBJetTrackAnalyzer::Plots::allIPErr, PatBJetTrackAnalyzer::Plots::allIPSig, B_JETS, C_JETS, reco::flavour(), flavours_, i, TFileService::make(), N_JET_TYPES, mergeVDriftHistosByStation::name, PatBJetTrackAnalyzer::Plots::negativeIP, PatBJetTrackAnalyzer::Plots::negativeIPErr, PatBJetTrackAnalyzer::Plots::negativeIPSig, PatBJetTrackAnalyzer::Plots::nTracks, RecoTauValidation_cfi::plots, plots_, PatBJetTrackAnalyzer::Plots::trackIP, PatBJetTrackAnalyzer::Plots::trackIPErr, PatBJetTrackAnalyzer::Plots::trackIPSig, 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.allIP = fs->make<TH1F>(Form("allIP_%s", name),
                                     Form("signed IP for all tracks in %s", flavour),
                                     100, -0.1, 0.2);
        plots.allIPErr = fs->make<TH1F>(Form("allIPErr_%s", name),
                                        Form("error of signed IP for all tracks in %s", flavour),
                                        100, 0, 0.05);
        plots.allIPSig = fs->make<TH1F>(Form("allIPSig_%s", name),
                                        Form("signed IP significance for all tracks in %s", flavour),
                                        100, -10, 20);

        plots.trackIP = fs->make<TH1F>(Form("trackIP_%s", name),
                                       Form("signed IP for selected positive track in %s", flavour),
                                       100, -0.1, 0.2);
        plots.trackIPErr = fs->make<TH1F>(Form("trackIPErr_%s", name),
                                          Form("error of signed IP for selected positive track in %s", flavour),
                                          100, 0, 0.05);
        plots.trackIPSig = fs->make<TH1F>(Form("trackIPSig_%s", name),
                                          Form("signed IP significance for selected positive track in %s", flavour),
                                          100, -10, 20);

        plots.negativeIP = fs->make<TH1F>(Form("negativeIP_%s", name),
                                          Form("signed IP for selected negative track in %s", flavour),
                                          100, -0.2, 0.1);
        plots.negativeIPErr = fs->make<TH1F>(Form("negativeIPErr_%s", name),
                                             Form("error of signed IP for selected negative track in %s", flavour),
                                             100, 0, 0.05);
        plots.negativeIPSig = fs->make<TH1F>(Form("negativeIPSig_%s", name),
                                             Form("signed IP significance for selected negative track in %s", flavour),
                                             100, -20, 10);

        plots.nTracks = fs->make<TH1F>(Form("nTracks_%s", name),
                                       Form("number of usable tracks in %s", flavour),
                                       30, 0, 30);
        plots.allDeltaR = fs->make<TH1F>(Form("allDeltaR_%s", name),
                                         Form("\\DeltaR between track and %s", flavour),
                                         100, 0, 1);
    }
}

Member Data Documentation

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

private

Definition at line 44 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

TH1* PatBJetTrackAnalyzer::flavours_

private

Definition at line 68 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze(), and beginJob().

double PatBJetTrackAnalyzer::jetEtaCut_

private

Definition at line 48 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

double PatBJetTrackAnalyzer::jetPtCut_

private

Definition at line 47 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

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

private

Definition at line 42 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

double PatBJetTrackAnalyzer::maxDeltaR_

private

Definition at line 49 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

unsigned int PatBJetTrackAnalyzer::minPixelHits_

private

Definition at line 52 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

double PatBJetTrackAnalyzer::minPt_

private

Definition at line 51 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

unsigned int PatBJetTrackAnalyzer::minTotalHits_

private

Definition at line 53 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

unsigned int PatBJetTrackAnalyzer::nThTrack_

private

Definition at line 55 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

struct PatBJetTrackAnalyzer::Plots PatBJetTrackAnalyzer::plots_[N_JET_TYPES]

private

Referenced by analyze(), and beginJob().

edm::EDGetTokenT<reco::VertexCollection> PatBJetTrackAnalyzer::primaryVerticesToken_

private

Definition at line 45 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

edm::EDGetTokenT<reco::TrackCollection> PatBJetTrackAnalyzer::tracksToken_

private

Definition at line 43 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().