PatBJetTrackAnalyzer Class Reference

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
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
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
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	wantsInputProcessBlocks ()
static bool	wantsProcessBlocks ()
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)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
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)
void	resetItemsToGetFrom (BranchType iType)

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 72 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 85 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 167 of file PatBJetTrackAnalyzer.cc.

References funct::abs(), ALL_JETS, PatBJetTrackAnalyzer::Plots::allDeltaR, PatBJetTrackAnalyzer::Plots::allIP, PatBJetTrackAnalyzer::Plots::allIPErr, PatBJetTrackAnalyzer::Plots::allIPSig, B_JETS, beam_dqm_sourceclient-live_cfg::beamSpot, beamSpotToken_, C_JETS, HLT_FULL_cff::DeltaR, HLT_FULL_cff::deltaR, Measurement1D::error(), 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_, MetAnalyzer::pv(), jetcorrextractor::sign(), Measurement1D::significance(), significanceHigher(), HLT_FULL_cff::track, 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 87 of file PatBJetTrackAnalyzer.cc.

References ALL_JETS, PatBJetTrackAnalyzer::Plots::allDeltaR, PatBJetTrackAnalyzer::Plots::allIP, PatBJetTrackAnalyzer::Plots::allIPErr, PatBJetTrackAnalyzer::Plots::allIPSig, B_JETS, C_JETS, flavours_, mps_fire::i, TFileService::make(), N_JET_TYPES, mergeVDriftHistosByStation::name, PatBJetTrackAnalyzer::Plots::negativeIP, PatBJetTrackAnalyzer::Plots::negativeIPErr, PatBJetTrackAnalyzer::Plots::negativeIPSig, PatBJetTrackAnalyzer::Plots::nTracks, L1TEGammaDiff_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 61 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().