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::one::EDAnalyzer< edm::one::SharedResources >
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const noexcept final
bool	wantsGlobalRuns () const noexcept final
bool	wantsInputProcessBlocks () const noexcept final
bool	wantsProcessBlocks () const noexcept final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const noexcept
bool	wantsStreamRuns () const noexcept
	~EDAnalyzerBase () 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
ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESResolverIndex > 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
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProductResolverIndices 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::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
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 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 ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
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

Flavour

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::PatBJetTrackAnalyzer

(

const edm::ParameterSet &

params

)

constructor and destructor

Definition at line 72 of file PatBJetTrackAnalyzer.cc.

References TFileService::kSharedResource.

    : 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")) {
  usesResource(TFileService::kSharedResource);
}

~PatBJetTrackAnalyzer()

PatBJetTrackAnalyzer::~PatBJetTrackAnalyzer ( )

override

Definition at line 87 of file PatBJetTrackAnalyzer.cc.

{}

Member Function Documentation

analyze()

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

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 169 of file PatBJetTrackAnalyzer.cc.

References funct::abs(), ALL_JETS, PatBJetTrackAnalyzer::Plots::allDeltaR, PatBJetTrackAnalyzer::Plots::allIP, PatBJetTrackAnalyzer::Plots::allIPErr, PatBJetTrackAnalyzer::Plots::allIPSig, B_JETS, pwdgSkimBPark_cfi::beamSpot, beamSpotToken_, C_JETS, electronAnalyzer_cfi::DeltaR, PbPb_ZMuSkimMuonDPG_cff::deltaR, Measurement1D::error(), flavours_, createfilelist::int, 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_, Validation_hcalonly_cfi::sign, Measurement1D::significance(), significanceHigher(), jetUpdater_cfi::sort, HLT_2024v13_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());
  }
}

beginJob()

void PatBJetTrackAnalyzer::beginJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 89 of file PatBJetTrackAnalyzer.cc.

References ALL_JETS, B_JETS, C_JETS, flavours_, compareTotals::fs, mps_fire::i, N_JET_TYPES, Skims_PA_cff::name, HLTObjectsMonitor_cfi::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.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

beamSpotToken_

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

private

Definition at line 44 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

flavours_

TH1* PatBJetTrackAnalyzer::flavours_

private

Definition at line 61 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze(), and beginJob().

jetEtaCut_

double PatBJetTrackAnalyzer::jetEtaCut_

private

Definition at line 48 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

jetPtCut_

double PatBJetTrackAnalyzer::jetPtCut_

private

Definition at line 47 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

jetsToken_

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

private

Definition at line 42 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

maxDeltaR_

double PatBJetTrackAnalyzer::maxDeltaR_

private

Definition at line 49 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

minPixelHits_

unsigned int PatBJetTrackAnalyzer::minPixelHits_

private

Definition at line 52 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

minPt_

double PatBJetTrackAnalyzer::minPt_

private

Definition at line 51 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

minTotalHits_

unsigned int PatBJetTrackAnalyzer::minTotalHits_

private

Definition at line 53 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

nThTrack_

unsigned int PatBJetTrackAnalyzer::nThTrack_

private

Definition at line 55 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

plots_

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

private

Referenced by analyze(), and beginJob().

primaryVerticesToken_

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

private

Definition at line 45 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().

tracksToken_

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

private

Definition at line 43 of file PatBJetTrackAnalyzer.cc.

Referenced by analyze().