Inheritance diagram for PatBJetTrackAnalyzer:

Classes
struct	Plots

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &es) override

void	beginJob () override

	PatBJetTrackAnalyzer (const edm::ParameterSet &params)
	constructor and destructor More...

	~PatBJetTrackAnalyzer () 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

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)

virtual	~EDConsumerBase () noexcept(false)

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

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 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 ( )

override

Definition at line 94 of file PatBJetTrackAnalyzer.cc.

95 {

96 }

Member Function Documentation

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

override

Definition at line 179 of file PatBJetTrackAnalyzer.cc.

 {
     // 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.

 {
     // 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().

Classes

Public Member Functions

Private Types

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation