#include <MultiTrackValidator.h>

Inheritance diagram for MultiTrackValidator:

Public Types
using	Histograms = MultiTrackValidatorHistograms

Public Types inherited from DQMGlobalEDAnalyzerBase< MultiTrackValidatorHistograms, Args...>
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &, Histograms &) const override
	Method called to book the DQM histograms. More...

void	dqmAnalyze (const edm::Event &, const edm::EventSetup &, const Histograms &) const override
	Method called once per event. More...

	MultiTrackValidator (const edm::ParameterSet &pset)
	Constructor. More...

	~MultiTrackValidator () override
	Destructor. More...

Public Member Functions inherited from DQMGlobalEDAnalyzer< MultiTrackValidatorHistograms >
virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &, MultiTrackValidatorHistogramsconst &) const

void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup) const final

Public Member Functions inherited from DQMGlobalEDAnalyzerBase< MultiTrackValidatorHistograms, Args...>
void	accumulate (edm::StreamID id, edm::Event const &event, edm::EventSetup const &setup) const final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &, MultiTrackValidatorHistograms &) const

	DQMGlobalEDAnalyzerBase ()

std::shared_ptr < MultiTrackValidatorHistograms >	globalBeginRun (edm::Run const &run, edm::EventSetup const &setup) const final

void	globalEndRun (edm::Run const &, edm::EventSetup const &) const final

Public Member Functions inherited from edm::global::EDProducer< T >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

Public Member Functions inherited from edm::global::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector < edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector < edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) 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)

Protected Attributes
std::vector< edm::InputTag >	associators

edm::EDGetTokenT< reco::BeamSpot >	bsSrc

const bool	calculateDrSingleCollection_

const bool	dodEdxPlots_

const bool	doMVAPlots_

const bool	doPlotsOnlyForTruePV_

const bool	doPVAssociationPlots_

const bool	doRecoTrackPlots_

std::vector< bool >	doResolutionPlots_

const bool	doSeedPlots_

const bool	doSimPlots_

const bool	doSimTrackPlots_

const bool	doSummaryPlots_

std::unique_ptr < MTVHistoProducerAlgoForTracker >	histoProducerAlgo_

const bool	ignoremissingtkcollection_

std::vector< edm::InputTag >	label

edm::EDGetTokenT< std::vector < PileupSummaryInfo > >	label_pileupinfo

edm::EDGetTokenT < TrackingParticleCollection >	label_tp_effic

edm::EDGetTokenT < TrackingParticleRefVector >	label_tp_effic_refvector

edm::EDGetTokenT < TrackingParticleCollection >	label_tp_fake

edm::EDGetTokenT < TrackingParticleRefVector >	label_tp_fake_refvector

edm::EDGetTokenT < TrackingVertexCollection >	label_tv

std::vector< edm::EDGetTokenT < edm::View< reco::Track > > >	labelToken

std::vector< edm::EDGetTokenT < edm::View< TrajectorySeed > > >	labelTokenSeed

edm::EDGetTokenT < edm::ValueMap < reco::DeDxData > >	m_dEdx1Tag

edm::EDGetTokenT < edm::ValueMap < reco::DeDxData > >	m_dEdx2Tag

const bool	parametersDefinerIsCosmic_

std::unique_ptr < ParametersDefinerForTP >	parametersDefinerTP_

std::vector< edm::EDGetTokenT < std::vector< PSimHit > > >	simHitTokens_

const edm::ESGetToken < TrackerTopology, TrackerTopologyRcd >	tTopoEsToken

const bool	useAssociators_

Protected Attributes inherited from DQMGlobalEDAnalyzerBase< MultiTrackValidatorHistograms, Args...>
DQMStore *	dqmstore_

edm::EDPutTokenT< DQMToken >	runToken_

Private Types
using	MVACollection = std::vector< float >

using	QualityMaskCollection = std::vector< unsigned char >

Private Member Functions
const reco::Vertex::Point *	getRecoPVPosition (const edm::Event &event, const edm::Handle< TrackingVertexCollection > &htv) const

const TrackingVertex::LorentzVector *	getSimPVPosition (const edm::Handle< TrackingVertexCollection > &htv) const

size_t	tpDR (const TrackingParticleRefVector &tPCeff, const std::vector< size_t > &selected_tPCeff, DynArray< float > &dR_tPCeff, DynArray< float > &dR_tPCeff_jet, const edm::View< reco::Candidate > *cores) const

void	tpParametersAndSelection (const Histograms &histograms, const TrackingParticleRefVector &tPCeff, const edm::Event &event, const edm::EventSetup &setup, const reco::BeamSpot &bs, std::vector< std::tuple< TrackingParticle::Vector, TrackingParticle::Point >> &momVert_tPCeff, std::vector< size_t > &selected_tPCeff) const

void	trackDR (const edm::View< reco::Track > &trackCollection, const edm::View< reco::Track > &trackCollectionDr, DynArray< float > &dR_trk, DynArray< float > &dR_trk_jet, const edm::View< reco::Candidate > *cores) const

Private Attributes
edm::EDGetTokenT < SimHitTPAssociationProducer::SimHitTPAssociationList >	_simHitTpMapTag

std::vector< edm::EDGetTokenT < reco::RecoToSimCollection > >	associatormapRtSs

std::vector< edm::EDGetTokenT < reco::SimToRecoCollection > >	associatormapStRs

std::vector< edm::EDGetTokenT < reco::TrackToTrackingParticleAssociator > >	associatorTokens

edm::EDGetTokenT< edm::View < reco::Candidate > >	cores_

CosmicTrackingParticleSelector	cosmictpSelector

std::string	dirName_

TrackingParticleSelector	dRtpSelector

std::unique_ptr < RecoTrackSelectorBase >	dRTrackSelector

edm::EDGetTokenT< edm::View < reco::Track > >	labelTokenForDrCalculation

std::vector< std::vector < std::tuple< edm::EDGetTokenT < MVACollection > , edm::EDGetTokenT < QualityMaskCollection > > > >	mvaQualityCollectionTokens_

double	ptMinJet_

edm::EDGetTokenT< edm::View < reco::Vertex > >	recoVertexToken_

const double	simPVMaxZ_

edm::EDGetTokenT < edm::ValueMap< unsigned int > >	tpNLayersToken_

edm::EDGetTokenT < edm::ValueMap< unsigned int > >	tpNPixelLayersToken_

edm::EDGetTokenT < edm::ValueMap< unsigned int > >	tpNStripStereoLayersToken_

TrackingParticleSelector	tpSelector

bool	useGsf

edm::EDGetTokenT < reco::VertexToTrackingVertexAssociator >	vertexAssociatorToken_

Additional Inherited Members
Static Public Member Functions inherited from edm::global::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from DQMGlobalEDAnalyzerBase< MultiTrackValidatorHistograms, Args...>
uint64_t	meId (edm::Run const &run) const

Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

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

Class that prodecs histrograms to validate Track Reconstruction performances

Author: cerati

Definition at line 39 of file MultiTrackValidator.h.

Member Typedef Documentation

using MultiTrackValidator::Histograms = MultiTrackValidatorHistograms

Definition at line 41 of file MultiTrackValidator.h.

using MultiTrackValidator::MVACollection = std::vector<float>

private

Definition at line 129 of file MultiTrackValidator.h.

using MultiTrackValidator::QualityMaskCollection = std::vector<unsigned char>

private

Definition at line 130 of file MultiTrackValidator.h.

Constructor & Destructor Documentation

MultiTrackValidator::MultiTrackValidator ( const edm::ParameterSet & pset )

Constructor.

Definition at line 48 of file MultiTrackValidator.cc.

References HLT_FULL_cff::beamSpotTag, SplitLinear::begin, bsSrc, edm::EDConsumerBase::consumes(), edm::EDConsumerBase::consumesCollector(), cores_, dirName_, dodEdxPlots_, doMVAPlots_, doPlotsOnlyForTruePV_, doPVAssociationPlots_, doResolutionPlots_, TrackValidation_cff::doResolutionPlotsForLabels, doSeedPlots_, relativeConstraints::empty, dataset::end, Exception, spr::find(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), histoProducerAlgo_, cmsLHEtoEOSManager::l, edm::InputTag::label(), label, label_pileupinfo, label_tp_effic, label_tp_effic_refvector, label_tp_fake, label_tp_fake_refvector, label_tv, HLT_FULL_cff::labels, edm::EDConsumerBase::labelsForToken(), labelToken, m_dEdx1Tag, m_dEdx2Tag, edm::ProductLabels::module, mvaQualityCollectionTokens_, parametersDefinerIsCosmic_, parametersDefinerTP_, recoVertexToken_, simHitTokens_, AlCaHLTBitMon_QueryRunRegistry::string, GlobalPosition_Frontier_DevDB_cff::tag, tpNLayersToken_, tpNPixelLayersToken_, tpNStripStereoLayersToken_, edm::vector_transform(), and vertexAssociatorToken_.

     : tTopoEsToken(esConsumes()),
       parametersDefinerIsCosmic_(pset.getParameter<std::string>("parametersDefiner") == "CosmicParametersDefinerForTP"),
       associators(pset.getUntrackedParameter<std::vector<edm::InputTag>>("associators")),
       label(pset.getParameter<std::vector<edm::InputTag>>("label")),
       ignoremissingtkcollection_(pset.getUntrackedParameter<bool>("ignoremissingtrackcollection", false)),
       useAssociators_(pset.getParameter<bool>("UseAssociators")),
       calculateDrSingleCollection_(pset.getUntrackedParameter<bool>("calculateDrSingleCollection")),
       doPlotsOnlyForTruePV_(pset.getUntrackedParameter<bool>("doPlotsOnlyForTruePV")),
       doSummaryPlots_(pset.getUntrackedParameter<bool>("doSummaryPlots")),
       doSimPlots_(pset.getUntrackedParameter<bool>("doSimPlots")),
       doSimTrackPlots_(pset.getUntrackedParameter<bool>("doSimTrackPlots")),
       doRecoTrackPlots_(pset.getUntrackedParameter<bool>("doRecoTrackPlots")),
       dodEdxPlots_(pset.getUntrackedParameter<bool>("dodEdxPlots")),
       doPVAssociationPlots_(pset.getUntrackedParameter<bool>("doPVAssociationPlots")),
       doSeedPlots_(pset.getUntrackedParameter<bool>("doSeedPlots")),
       doMVAPlots_(pset.getUntrackedParameter<bool>("doMVAPlots")),
       simPVMaxZ_(pset.getUntrackedParameter<double>("simPVMaxZ")) {
   if (not(pset.getParameter<edm::InputTag>("cores").label().empty())) {
     cores_ = consumes<edm::View<reco::Candidate>>(pset.getParameter<edm::InputTag>("cores"));
   }
   if (label.empty()) {
     // Disable prefetching of everything if there are no track collections
     return;
   }
 
   const edm::InputTag& label_tp_effic_tag = pset.getParameter<edm::InputTag>("label_tp_effic");
   const edm::InputTag& label_tp_fake_tag = pset.getParameter<edm::InputTag>("label_tp_fake");
 
   if (pset.getParameter<bool>("label_tp_effic_refvector")) {
     label_tp_effic_refvector = consumes<TrackingParticleRefVector>(label_tp_effic_tag);
   } else {
     label_tp_effic = consumes<TrackingParticleCollection>(label_tp_effic_tag);
   }
   if (pset.getParameter<bool>("label_tp_fake_refvector")) {
     label_tp_fake_refvector = consumes<TrackingParticleRefVector>(label_tp_fake_tag);
   } else {
     label_tp_fake = consumes<TrackingParticleCollection>(label_tp_fake_tag);
   }
   label_pileupinfo = consumes<std::vector<PileupSummaryInfo>>(pset.getParameter<edm::InputTag>("label_pileupinfo"));
   for (const auto& tag : pset.getParameter<std::vector<edm::InputTag>>("sim")) {
     simHitTokens_.push_back(consumes<std::vector<PSimHit>>(tag));
   }
 
   std::vector<edm::InputTag> doResolutionPlotsForLabels =
       pset.getParameter<std::vector<edm::InputTag>>("doResolutionPlotsForLabels");
   doResolutionPlots_.reserve(label.size());
   for (auto& itag : label) {
     labelToken.push_back(consumes<edm::View<reco::Track>>(itag));
     const bool doResol = doResolutionPlotsForLabels.empty() ||
                          (std::find(cbegin(doResolutionPlotsForLabels), cend(doResolutionPlotsForLabels), itag) !=
                           cend(doResolutionPlotsForLabels));
     doResolutionPlots_.push_back(doResol);
   }
   {  // check for duplicates
     auto labelTmp = edm::vector_transform(label, [&](const edm::InputTag& tag) { return tag.label(); });
     std::sort(begin(labelTmp), end(labelTmp));
     std::string empty;
     const std::string* prev = &empty;
     for (const std::string& l : labelTmp) {
       if (l == *prev) {
         throw cms::Exception("Configuration") << "Duplicate InputTag in labels: " << l;
       }
       prev = &l;
     }
   }
 
   edm::InputTag beamSpotTag = pset.getParameter<edm::InputTag>("beamSpot");
   bsSrc = consumes<reco::BeamSpot>(beamSpotTag);
 
   if (parametersDefinerIsCosmic_) {
     parametersDefinerTP_ = std::make_unique<CosmicParametersDefinerForTP>(consumesCollector());
   } else {
     parametersDefinerTP_ = std::make_unique<ParametersDefinerForTP>(beamSpotTag, consumesCollector());
   }
 
   ParameterSet psetForHistoProducerAlgo = pset.getParameter<ParameterSet>("histoProducerAlgoBlock");
   histoProducerAlgo_ = std::make_unique<MTVHistoProducerAlgoForTracker>(psetForHistoProducerAlgo, doSeedPlots_);
 
   dirName_ = pset.getParameter<std::string>("dirName");
 
   tpNLayersToken_ = consumes<edm::ValueMap<unsigned int>>(pset.getParameter<edm::InputTag>("label_tp_nlayers"));
   tpNPixelLayersToken_ =
       consumes<edm::ValueMap<unsigned int>>(pset.getParameter<edm::InputTag>("label_tp_npixellayers"));
   tpNStripStereoLayersToken_ =
       consumes<edm::ValueMap<unsigned int>>(pset.getParameter<edm::InputTag>("label_tp_nstripstereolayers"));
 
   if (dodEdxPlots_) {
     m_dEdx1Tag = consumes<edm::ValueMap<reco::DeDxData>>(pset.getParameter<edm::InputTag>("dEdx1Tag"));
     m_dEdx2Tag = consumes<edm::ValueMap<reco::DeDxData>>(pset.getParameter<edm::InputTag>("dEdx2Tag"));
   }
 
   label_tv = consumes<TrackingVertexCollection>(pset.getParameter<edm::InputTag>("label_tv"));
   if (doPlotsOnlyForTruePV_ || doPVAssociationPlots_) {
     recoVertexToken_ = consumes<edm::View<reco::Vertex>>(pset.getUntrackedParameter<edm::InputTag>("label_vertex"));
     vertexAssociatorToken_ =
         consumes<reco::VertexToTrackingVertexAssociator>(pset.getUntrackedParameter<edm::InputTag>("vertexAssociator"));
   }
 
   if (doMVAPlots_) {
     mvaQualityCollectionTokens_.resize(labelToken.size());
     auto mvaPSet = pset.getUntrackedParameter<edm::ParameterSet>("mvaLabels");
     for (size_t iIter = 0; iIter < labelToken.size(); ++iIter) {
       edm::EDConsumerBase::Labels labels;
       labelsForToken(labelToken[iIter], labels);
       if (mvaPSet.exists(labels.module)) {
         mvaQualityCollectionTokens_[iIter] = edm::vector_transform(
             mvaPSet.getUntrackedParameter<std::vector<std::string>>(labels.module), [&](const std::string& tag) {
               return std::make_tuple(consumes<MVACollection>(edm::InputTag(tag, "MVAValues")),
                                      consumes<QualityMaskCollection>(edm::InputTag(tag, "QualityMasks")));
             });
       }
     }
   }
 
   tpSelector = TrackingParticleSelector(pset.getParameter<double>("ptMinTP"),
                                         pset.getParameter<double>("ptMaxTP"),
                                         pset.getParameter<double>("minRapidityTP"),
                                         pset.getParameter<double>("maxRapidityTP"),
                                         pset.getParameter<double>("tipTP"),
                                         pset.getParameter<double>("lipTP"),
                                         pset.getParameter<int>("minHitTP"),
                                         pset.getParameter<bool>("signalOnlyTP"),
                                         pset.getParameter<bool>("intimeOnlyTP"),
                                         pset.getParameter<bool>("chargedOnlyTP"),
                                         pset.getParameter<bool>("stableOnlyTP"),
                                         pset.getParameter<std::vector<int>>("pdgIdTP"),
                                         pset.getParameter<bool>("invertRapidityCutTP"));
 
   cosmictpSelector = CosmicTrackingParticleSelector(pset.getParameter<double>("ptMinTP"),
                                                     pset.getParameter<double>("minRapidityTP"),
                                                     pset.getParameter<double>("maxRapidityTP"),
                                                     pset.getParameter<double>("tipTP"),
                                                     pset.getParameter<double>("lipTP"),
                                                     pset.getParameter<int>("minHitTP"),
                                                     pset.getParameter<bool>("chargedOnlyTP"),
                                                     pset.getParameter<std::vector<int>>("pdgIdTP"));
 
   ParameterSet psetVsPhi = psetForHistoProducerAlgo.getParameter<ParameterSet>("TpSelectorForEfficiencyVsPhi");
   dRtpSelector = TrackingParticleSelector(psetVsPhi.getParameter<double>("ptMin"),
                                           psetVsPhi.getParameter<double>("ptMax"),
                                           psetVsPhi.getParameter<double>("minRapidity"),
                                           psetVsPhi.getParameter<double>("maxRapidity"),
                                           psetVsPhi.getParameter<double>("tip"),
                                           psetVsPhi.getParameter<double>("lip"),
                                           psetVsPhi.getParameter<int>("minHit"),
                                           psetVsPhi.getParameter<bool>("signalOnly"),
                                           psetVsPhi.getParameter<bool>("intimeOnly"),
                                           psetVsPhi.getParameter<bool>("chargedOnly"),
                                           psetVsPhi.getParameter<bool>("stableOnly"),
                                           psetVsPhi.getParameter<std::vector<int>>("pdgId"),
                                           psetVsPhi.getParameter<bool>("invertRapidityCut"));
 
   dRTrackSelector = MTVHistoProducerAlgoForTracker::makeRecoTrackSelectorFromTPSelectorParameters(psetVsPhi);
 
   useGsf = pset.getParameter<bool>("useGsf");
 
   _simHitTpMapTag = mayConsume<SimHitTPAssociationProducer::SimHitTPAssociationList>(
       pset.getParameter<edm::InputTag>("simHitTpMapTag"));
 
   if (calculateDrSingleCollection_) {
     labelTokenForDrCalculation =
         consumes<edm::View<reco::Track>>(pset.getParameter<edm::InputTag>("trackCollectionForDrCalculation"));
   }
 
   if (useAssociators_) {
     for (auto const& src : associators) {
       associatorTokens.push_back(consumes<reco::TrackToTrackingParticleAssociator>(src));
     }
   } else {
     for (auto const& src : associators) {
       associatormapStRs.push_back(consumes<reco::SimToRecoCollection>(src));
       associatormapRtSs.push_back(consumes<reco::RecoToSimCollection>(src));
     }
   }
 }

MultiTrackValidator::~MultiTrackValidator ( )

override

Destructor.

Definition at line 225 of file MultiTrackValidator.cc.

225 {}

Member Function Documentation

void MultiTrackValidator::bookHistograms	(	DQMStore::IBooker &	ibook,
		edm::Run const &	,
		edm::EventSetup const &	setup,
		Histograms &	histograms
	)		const

overridevirtual

Method called to book the DQM histograms.

Implements DQMGlobalEDAnalyzerBase< MultiTrackValidatorHistograms, Args...>.

Definition at line 227 of file MultiTrackValidator.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::NavigatorBase::cd(), TrackerOfflineValidation_Dqm_cff::dirName, dirName_, dodEdxPlots_, doMVAPlots_, doPVAssociationPlots_, doRecoTrackPlots_, doResolutionPlots_, doSeedPlots_, doSimPlots_, doSimTrackPlots_, doSummaryPlots_, MultiTrackValidatorHistograms::h_assoc2_coll, MultiTrackValidatorHistograms::h_assoc_coll, MultiTrackValidatorHistograms::h_looper_coll, MultiTrackValidatorHistograms::h_pileup_coll, MultiTrackValidatorHistograms::h_reco_coll, MultiTrackValidatorHistograms::h_simul_coll, MultiTrackValidatorHistograms::histoProducerAlgo, histoProducerAlgo_, mps_fire::i, edm::InputTag::instance(), edm::InputTag::label(), label, hlt_dqm_clientPB-live_cfg::me, mvaQualityCollectionTokens_, edm::InputTag::process(), python.rootplot.root2matplotlib::replace(), and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                                        {
   if (label.empty()) {
     // Disable histogram booking if there are no track collections
     return;
   }
 
   const auto minColl = -0.5;
   const auto maxColl = label.size() - 0.5;
   const auto nintColl = label.size();
 
   auto binLabels = [&](dqm::reco::MonitorElement* me) {
     for (size_t i = 0; i < label.size(); ++i) {
       me->setBinLabel(i + 1, label[i].label());
     }
     me->disableAlphanumeric();
     return me;
   };
 
   //Booking histograms concerning with simulated tracks
   if (doSimPlots_) {
     ibook.cd();
     ibook.setCurrentFolder(dirName_ + "simulation");
 
     histoProducerAlgo_->bookSimHistos(ibook, histograms.histoProducerAlgo);
 
     ibook.cd();
     ibook.setCurrentFolder(dirName_);
   }
 
   for (unsigned int ww = 0; ww < associators.size(); ww++) {
     ibook.cd();
     // FIXME: these need to be moved to a subdirectory whose name depends on the associator
     ibook.setCurrentFolder(dirName_);
 
     if (doSummaryPlots_) {
       if (doSimTrackPlots_) {
         histograms.h_assoc_coll.push_back(
             binLabels(ibook.book1D("num_assoc(simToReco)_coll",
                                    "N of associated (simToReco) tracks vs track collection",
                                    nintColl,
                                    minColl,
                                    maxColl)));
         histograms.h_simul_coll.push_back(binLabels(
             ibook.book1D("num_simul_coll", "N of simulated tracks vs track collection", nintColl, minColl, maxColl)));
       }
       if (doRecoTrackPlots_) {
         histograms.h_reco_coll.push_back(binLabels(
             ibook.book1D("num_reco_coll", "N of reco track vs track collection", nintColl, minColl, maxColl)));
         histograms.h_assoc2_coll.push_back(
             binLabels(ibook.book1D("num_assoc(recoToSim)_coll",
                                    "N of associated (recoToSim) tracks vs track collection",
                                    nintColl,
                                    minColl,
                                    maxColl)));
         histograms.h_looper_coll.push_back(
             binLabels(ibook.book1D("num_duplicate_coll",
                                    "N of associated (recoToSim) looper tracks vs track collection",
                                    nintColl,
                                    minColl,
                                    maxColl)));
         histograms.h_pileup_coll.push_back(
             binLabels(ibook.book1D("num_pileup_coll",
                                    "N of associated (recoToSim) pileup tracks vs track collection",
                                    nintColl,
                                    minColl,
                                    maxColl)));
       }
     }
 
     for (unsigned int www = 0; www < label.size(); www++) {
       ibook.cd();
       InputTag algo = label[www];
       string dirName = dirName_;
       if (!algo.process().empty())
         dirName += algo.process() + "_";
       if (!algo.label().empty())
         dirName += algo.label() + "_";
       if (!algo.instance().empty())
         dirName += algo.instance() + "_";
       if (dirName.find("Tracks") < dirName.length()) {
         dirName.replace(dirName.find("Tracks"), 6, "");
       }
       string assoc = associators[ww].label();
       if (assoc.find("Track") < assoc.length()) {
         assoc.replace(assoc.find("Track"), 5, "");
       }
       dirName += assoc;
       std::replace(dirName.begin(), dirName.end(), ':', '_');
 
       ibook.setCurrentFolder(dirName);
 
       const bool doResolutionPlots = doResolutionPlots_[www];
 
       if (doSimTrackPlots_) {
         histoProducerAlgo_->bookSimTrackHistos(ibook, histograms.histoProducerAlgo, doResolutionPlots);
         if (doPVAssociationPlots_)
           histoProducerAlgo_->bookSimTrackPVAssociationHistos(ibook, histograms.histoProducerAlgo);
       }
 
       //Booking histograms concerning with reconstructed tracks
       if (doRecoTrackPlots_) {
         histoProducerAlgo_->bookRecoHistos(ibook, histograms.histoProducerAlgo, doResolutionPlots);
         if (dodEdxPlots_)
           histoProducerAlgo_->bookRecodEdxHistos(ibook, histograms.histoProducerAlgo);
         if (doPVAssociationPlots_)
           histoProducerAlgo_->bookRecoPVAssociationHistos(ibook, histograms.histoProducerAlgo);
         if (doMVAPlots_)
           histoProducerAlgo_->bookMVAHistos(
               ibook, histograms.histoProducerAlgo, mvaQualityCollectionTokens_[www].size());
       }
 
       if (doSeedPlots_) {
         histoProducerAlgo_->bookSeedHistos(ibook, histograms.histoProducerAlgo);
       }
     }  //end loop www
   }    // end loop ww
 }

void MultiTrackValidator::dqmAnalyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup,
		const Histograms &	histograms
	)		const

overridevirtual

Method called once per event.

Implements DQMGlobalEDAnalyzerBase< MultiTrackValidatorHistograms, Args...>.

Reimplemented in MultiTrackValidatorGenPs.

Definition at line 583 of file MultiTrackValidator.cc.

References _simHitTpMapTag, funct::abs(), associationMapFilterValues(), associatormapRtSs, associatormapStRs, associatorTokens, cms::cuda::bs, bsSrc, EncodedEventId::bunchCrossing(), calculateDrSingleCollection_, HLT_FULL_cff::cores, cores_, cosmictpSelector, declareDynArray, dodEdxPlots_, doMVAPlots_, doPlotsOnlyForTruePV_, doPVAssociationPlots_, doRecoTrackPlots_, doResolutionPlots_, doSeedPlots_, doSimTrackPlots_, doSummaryPlots_, dRtpSelector, dRTrackSelector, TrackingParticleIP::dxy(), TrackingParticleIP::dz(), relativeConstraints::empty, edm::AssociationMap< Tag >::end(), EncodedEventId::event(), TrackingParticle::eventId(), Exception, edm::AssociationMap< Tag >::find(), first, edm::Event::getByToken(), edm::EventSetup::getData(), getRecoPVPosition(), getSimPVPosition(), MultiTrackValidatorHistograms::h_assoc2_coll, MultiTrackValidatorHistograms::h_assoc_coll, MultiTrackValidatorHistograms::h_looper_coll, MultiTrackValidatorHistograms::h_pileup_coll, MultiTrackValidatorHistograms::h_reco_coll, MultiTrackValidatorHistograms::h_simul_coll, reco::TrackBase::highPurity, MultiTrackValidatorHistograms::histoProducerAlgo, histoProducerAlgo_, mps_fire::i, ignoremissingtkcollection_, CosmicTrackingParticleSelector::initEvent(), edm::EDGetTokenT< T >::isUninitialized(), edm::HandleBase::isValid(), label, label_pileupinfo, label_tp_effic, label_tp_effic_refvector, label_tp_fake, label_tp_fake_refvector, label_tv, labelToken, labelTokenForDrCalculation, LogDebug, LogTrace, reco::TrackBase::loose, m_dEdx1Tag, m_dEdx2Tag, SiStripPI::max, TrackingParticle::momentum(), beam_dqm_sourceclient-live_cfg::mva, mvaQualityCollectionTokens_, TrackingParticle::numberOfTrackerHits(), parametersDefinerIsCosmic_, parametersDefinerTP_, edm::Handle< T >::product(), edm::RefToBaseVector< T >::push_back(), dt_dqm_sourceclient_common_cff::reco, edm::View< T >::refAt(), hcal_runs::rt, simPVMaxZ_, edm::View< T >::size(), findQualityFiles::size, mathSSE::sqrt(), cmsswSequenceInfo::tp, tpDR(), tpNLayersToken_, tpNPixelLayersToken_, tpNStripStereoLayersToken_, tpParametersAndSelection(), HLT_FULL_cff::track, HLT_FULL_cff::trackCollection, trackDR(), trackFromSeedFitFailed(), tTopoEsToken, useAssociators_, TrackingParticle::vertex(), and w.

                                                                          {
   if (label.empty()) {
     // Disable if there are no track collections
     return;
   }
 
   using namespace reco;
 
   LogDebug("TrackValidator") << "\n===================================================="
                              << "\n"
                              << "Analyzing new event"
                              << "\n"
                              << "====================================================\n"
                              << "\n";
 
   const TrackerTopology& ttopo = setup.getData(tTopoEsToken);
 
   // FIXME: we really need to move to edm::View for reading the
   // TrackingParticles... Unfortunately it has non-trivial
   // consequences on the associator/association interfaces etc.
   TrackingParticleRefVector tmpTPeff;
   TrackingParticleRefVector tmpTPfake;
   const TrackingParticleRefVector* tmpTPeffPtr = nullptr;
   const TrackingParticleRefVector* tmpTPfakePtr = nullptr;
 
   edm::Handle<TrackingParticleCollection> TPCollectionHeff;
   edm::Handle<TrackingParticleRefVector> TPCollectionHeffRefVector;
 
   const bool tp_effic_refvector = label_tp_effic.isUninitialized();
   if (!tp_effic_refvector) {
     event.getByToken(label_tp_effic, TPCollectionHeff);
     tmpTPeff.reserve(TPCollectionHeff->size());
     for (size_t i = 0, size = TPCollectionHeff->size(); i < size; ++i) {
       tmpTPeff.push_back(TrackingParticleRef(TPCollectionHeff, i));
     }
     tmpTPeffPtr = &tmpTPeff;
   } else {
     event.getByToken(label_tp_effic_refvector, TPCollectionHeffRefVector);
     tmpTPeffPtr = TPCollectionHeffRefVector.product();
   }
   if (!label_tp_fake.isUninitialized()) {
     edm::Handle<TrackingParticleCollection> TPCollectionHfake;
     event.getByToken(label_tp_fake, TPCollectionHfake);
     tmpTPfake.reserve(TPCollectionHfake->size());
     for (size_t i = 0, size = TPCollectionHfake->size(); i < size; ++i) {
       tmpTPfake.push_back(TrackingParticleRef(TPCollectionHfake, i));
     }
     tmpTPfakePtr = &tmpTPfake;
   } else {
     edm::Handle<TrackingParticleRefVector> TPCollectionHfakeRefVector;
     event.getByToken(label_tp_fake_refvector, TPCollectionHfakeRefVector);
     tmpTPfakePtr = TPCollectionHfakeRefVector.product();
   }
 
   TrackingParticleRefVector const& tPCeff = *tmpTPeffPtr;
   TrackingParticleRefVector const& tPCfake = *tmpTPfakePtr;
 
 #ifdef EDM_ML_DEBUG
   ensureEffIsSubsetOfFake(tPCeff, tPCfake);
 #endif
 
   if (parametersDefinerIsCosmic_) {
     edm::Handle<SimHitTPAssociationProducer::SimHitTPAssociationList> simHitsTPAssoc;
     //warning: make sure the TP collection used in the map is the same used in the MTV!
     event.getByToken(_simHitTpMapTag, simHitsTPAssoc);
     parametersDefinerTP_->initEvent(simHitsTPAssoc);
     cosmictpSelector.initEvent(simHitsTPAssoc);
   }
 
   // Find the sim PV and tak its position
   edm::Handle<TrackingVertexCollection> htv;
   event.getByToken(label_tv, htv);
   const TrackingVertex::LorentzVector* theSimPVPosition = getSimPVPosition(htv);
   if (simPVMaxZ_ >= 0) {
     if (!theSimPVPosition)
       return;
     if (std::abs(theSimPVPosition->z()) > simPVMaxZ_)
       return;
   }
 
   // Check, when necessary, if reco PV matches to sim PV
   const reco::Vertex::Point* thePVposition = nullptr;
   if (doPlotsOnlyForTruePV_ || doPVAssociationPlots_) {
     thePVposition = getRecoPVPosition(event, htv);
     if (doPlotsOnlyForTruePV_ && !thePVposition)
       return;
 
     // Rest of the code assumes that if thePVposition is non-null, the
     // PV-association histograms get filled. In above, the "nullness"
     // is used to deliver the information if the reco PV is matched to
     // the sim PV.
     if (!doPVAssociationPlots_)
       thePVposition = nullptr;
   }
 
   edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
   event.getByToken(bsSrc, recoBeamSpotHandle);
   reco::BeamSpot const& bs = *recoBeamSpotHandle;
 
   edm::Handle<std::vector<PileupSummaryInfo>> puinfoH;
   event.getByToken(label_pileupinfo, puinfoH);
   PileupSummaryInfo puinfo;
 
   for (unsigned int puinfo_ite = 0; puinfo_ite < (*puinfoH).size(); ++puinfo_ite) {
     if ((*puinfoH)[puinfo_ite].getBunchCrossing() == 0) {
       puinfo = (*puinfoH)[puinfo_ite];
       break;
     }
   }
 
   // Number of 3D layers for TPs
   edm::Handle<edm::ValueMap<unsigned int>> tpNLayersH;
   event.getByToken(tpNLayersToken_, tpNLayersH);
   const auto& nLayers_tPCeff = *tpNLayersH;
 
   event.getByToken(tpNPixelLayersToken_, tpNLayersH);
   const auto& nPixelLayers_tPCeff = *tpNLayersH;
 
   event.getByToken(tpNStripStereoLayersToken_, tpNLayersH);
   const auto& nStripMonoAndStereoLayers_tPCeff = *tpNLayersH;
 
   // Precalculate TP selection (for efficiency), and momentum and vertex wrt PCA
   //
   // TODO: ParametersDefinerForTP ESProduct needs to be changed to
   // EDProduct because of consumes.
   //
   // In principle, we could just precalculate the momentum and vertex
   // wrt PCA for all TPs for once and put that to the event. To avoid
   // repetitive calculations those should be calculated only once for
   // each TP. That would imply that we should access TPs via Refs
   // (i.e. View) in here, since, in general, the eff and fake TP
   // collections can be different (and at least HI seems to use that
   // feature). This would further imply that the
   // RecoToSimCollection/SimToRecoCollection should be changed to use
   // View<TP> instead of vector<TP>, and migrate everything.
   //
   // Or we could take only one input TP collection, and do another
   // TP-selection to obtain the "fake" collection like we already do
   // for "efficiency" TPs.
   std::vector<size_t> selected_tPCeff;
   std::vector<std::tuple<TrackingParticle::Vector, TrackingParticle::Point>> momVert_tPCeff;
   tpParametersAndSelection(histograms, tPCeff, event, setup, bs, momVert_tPCeff, selected_tPCeff);
 
   //calculate dR for TPs
   declareDynArray(float, tPCeff.size(), dR_tPCeff);
 
   //calculate dR_jet for TPs
   const edm::View<reco::Candidate>* coresVector = nullptr;
   if (not cores_.isUninitialized()) {
     Handle<edm::View<reco::Candidate>> cores;
     event.getByToken(cores_, cores);
     if (cores.isValid()) {
       coresVector = cores.product();
     }
   }
   declareDynArray(float, tPCeff.size(), dR_tPCeff_jet);
 
   size_t n_selTP_dr = tpDR(tPCeff, selected_tPCeff, dR_tPCeff, dR_tPCeff_jet, coresVector);
 
   edm::Handle<View<Track>> trackCollectionForDrCalculation;
   if (calculateDrSingleCollection_) {
     event.getByToken(labelTokenForDrCalculation, trackCollectionForDrCalculation);
   }
 
   // dE/dx
   // at some point this could be generalized, with a vector of tags and a corresponding vector of Handles
   // I'm writing the interface such to take vectors of ValueMaps
   std::vector<const edm::ValueMap<reco::DeDxData>*> v_dEdx;
   if (dodEdxPlots_) {
     edm::Handle<edm::ValueMap<reco::DeDxData>> dEdx1Handle;
     edm::Handle<edm::ValueMap<reco::DeDxData>> dEdx2Handle;
     event.getByToken(m_dEdx1Tag, dEdx1Handle);
     event.getByToken(m_dEdx2Tag, dEdx2Handle);
     v_dEdx.push_back(dEdx1Handle.product());
     v_dEdx.push_back(dEdx2Handle.product());
   }
 
   std::vector<const MVACollection*> mvaCollections;
   std::vector<const QualityMaskCollection*> qualityMaskCollections;
   std::vector<float> mvaValues;
 
   int w = 0;  //counter counting the number of sets of histograms
   for (unsigned int ww = 0; ww < associators.size(); ww++) {
     // run value filtering of recoToSim map already here as it depends only on the association, not track collection
     reco::SimToRecoCollection const* simRecCollPFull = nullptr;
     reco::RecoToSimCollection const* recSimCollP = nullptr;
     reco::RecoToSimCollection recSimCollL;
     if (!useAssociators_) {
       Handle<reco::SimToRecoCollection> simtorecoCollectionH;
       event.getByToken(associatormapStRs[ww], simtorecoCollectionH);
       simRecCollPFull = simtorecoCollectionH.product();
 
       Handle<reco::RecoToSimCollection> recotosimCollectionH;
       event.getByToken(associatormapRtSs[ww], recotosimCollectionH);
       recSimCollP = recotosimCollectionH.product();
 
       // We need to filter the associations of the fake-TrackingParticle
       // collection only from RecoToSim collection, otherwise the
       // RecoToSim histograms get false entries
       recSimCollL = associationMapFilterValues(*recSimCollP, tPCfake);
       recSimCollP = &recSimCollL;
     }
 
     for (unsigned int www = 0; www < label.size();
          www++, w++) {  // need to increment w here, since there are many continues in the loop body
       //
       //get collections from the event
       //
       edm::Handle<View<Track>> trackCollectionHandle;
       if (!event.getByToken(labelToken[www], trackCollectionHandle) && ignoremissingtkcollection_)
         continue;
       const edm::View<Track>& trackCollection = *trackCollectionHandle;
 
       reco::SimToRecoCollection const* simRecCollP = nullptr;
       reco::SimToRecoCollection simRecCollL;
 
       //associate tracks
       LogTrace("TrackValidator") << "Analyzing " << label[www] << " with " << associators[ww] << "\n";
       if (useAssociators_) {
         edm::Handle<reco::TrackToTrackingParticleAssociator> theAssociator;
         event.getByToken(associatorTokens[ww], theAssociator);
 
         // The associator interfaces really need to be fixed...
         edm::RefToBaseVector<reco::Track> trackRefs;
         // trackRefs.vectorHolder()->reserve(trackCollection.size());  NOT a good idea
         for (edm::View<Track>::size_type i = 0; i < trackCollection.size(); ++i) {
           trackRefs.push_back(trackCollection.refAt(i));
         }
 
         LogTrace("TrackValidator") << "Calling associateRecoToSim method"
                                    << "\n";
         recSimCollL = theAssociator->associateRecoToSim(trackRefs, tPCfake);
         recSimCollP = &recSimCollL;
         LogTrace("TrackValidator") << "Calling associateSimToReco method"
                                    << "\n";
         // It is necessary to do the association wrt. fake TPs,
         // because this SimToReco association is used also for
         // duplicates. Since the set of efficiency TPs are required to
         // be a subset of the set of fake TPs, for efficiency
         // histograms it doesn't matter if the association contains
         // associations of TPs not in the set of efficiency TPs.
         simRecCollL = theAssociator->associateSimToReco(trackRefs, tPCfake);
         simRecCollP = &simRecCollL;
       } else {
         // We need to filter the associations of the current track
         // collection only from SimToReco collection, otherwise the
         // SimToReco histograms get false entries. The filtering must
         // be done separately for each track collection.
         simRecCollL = associationMapFilterValues(*simRecCollPFull, trackCollection);
         simRecCollP = &simRecCollL;
       }
 
       reco::RecoToSimCollection const& recSimColl = *recSimCollP;
       reco::SimToRecoCollection const& simRecColl = *simRecCollP;
 
       // read MVA collections
       if (doMVAPlots_ && !mvaQualityCollectionTokens_[www].empty()) {
         edm::Handle<MVACollection> hmva;
         edm::Handle<QualityMaskCollection> hqual;
         for (const auto& tokenTpl : mvaQualityCollectionTokens_[www]) {
           event.getByToken(std::get<0>(tokenTpl), hmva);
           event.getByToken(std::get<1>(tokenTpl), hqual);
 
           mvaCollections.push_back(hmva.product());
           qualityMaskCollections.push_back(hqual.product());
           if (mvaCollections.back()->size() != trackCollection.size()) {
             throw cms::Exception("Configuration")
                 << "Inconsistency in track collection and MVA sizes. Track collection " << www << " has "
                 << trackCollection.size() << " tracks, whereas the MVA " << (mvaCollections.size() - 1)
                 << " for it has " << mvaCollections.back()->size() << " entries. Double-check your configuration.";
           }
           if (qualityMaskCollections.back()->size() != trackCollection.size()) {
             throw cms::Exception("Configuration")
                 << "Inconsistency in track collection and quality mask sizes. Track collection " << www << " has "
                 << trackCollection.size() << " tracks, whereas the quality mask " << (qualityMaskCollections.size() - 1)
                 << " for it has " << qualityMaskCollections.back()->size()
                 << " entries. Double-check your configuration.";
           }
         }
       }
 
       // ########################################################
       // fill simulation histograms (LOOP OVER TRACKINGPARTICLES)
       // ########################################################
 
       //compute number of tracks per eta interval
       //
       LogTrace("TrackValidator") << "\n# of TrackingParticles: " << tPCeff.size() << "\n";
       int ats(0);  //This counter counts the number of simTracks that are "associated" to recoTracks
       int st(0);   //This counter counts the number of simulated tracks passing the MTV selection (i.e. tpSelector(tp) )
 
       //loop over already-selected TPs for tracking efficiency
       for (size_t i = 0; i < selected_tPCeff.size(); ++i) {
         size_t iTP = selected_tPCeff[i];
         const TrackingParticleRef& tpr = tPCeff[iTP];
         const TrackingParticle& tp = *tpr;
 
         auto const& momVert = momVert_tPCeff[i];
         TrackingParticle::Vector momentumTP;
         TrackingParticle::Point vertexTP;
 
         double dxySim(0);
         double dzSim(0);
         double dxyPVSim = 0;
         double dzPVSim = 0;
         double dR = dR_tPCeff[iTP];
         double dR_jet = dR_tPCeff_jet[iTP];
 
         //---------- THIS PART HAS TO BE CLEANED UP. THE PARAMETER DEFINER WAS NOT MEANT TO BE USED IN THIS WAY ----------
         //If the TrackingParticle is collison like, get the momentum and vertex at production state
         if (!parametersDefinerIsCosmic_) {
           momentumTP = tp.momentum();
           vertexTP = tp.vertex();
           //Calcualte the impact parameters w.r.t. PCA
           const TrackingParticle::Vector& momentum = std::get<TrackingParticle::Vector>(momVert);
           const TrackingParticle::Point& vertex = std::get<TrackingParticle::Point>(momVert);
           dxySim = TrackingParticleIP::dxy(vertex, momentum, bs.position());
           dzSim = TrackingParticleIP::dz(vertex, momentum, bs.position());
 
           if (theSimPVPosition) {
             dxyPVSim = TrackingParticleIP::dxy(vertex, momentum, *theSimPVPosition);
             dzPVSim = TrackingParticleIP::dz(vertex, momentum, *theSimPVPosition);
           }
         }
         //If the TrackingParticle is comics, get the momentum and vertex at PCA
         else {
           momentumTP = std::get<TrackingParticle::Vector>(momVert);
           vertexTP = std::get<TrackingParticle::Point>(momVert);
           dxySim = TrackingParticleIP::dxy(vertexTP, momentumTP, bs.position());
           dzSim = TrackingParticleIP::dz(vertexTP, momentumTP, bs.position());
 
           // Do dxy and dz vs. PV make any sense for cosmics? I guess not
         }
         //---------- THE PART ABOVE HAS TO BE CLEANED UP. THE PARAMETER DEFINER WAS NOT MEANT TO BE USED IN THIS WAY ----------
 
         // in the coming lines, histos are filled using as input
         // - momentumTP
         // - vertexTP
         // - dxySim
         // - dzSim
         if (!doSimTrackPlots_)
           continue;
 
         // ##############################################
         // fill RecoAssociated SimTracks' histograms
         // ##############################################
         const reco::Track* matchedTrackPointer = nullptr;
         const reco::Track* matchedSecondTrackPointer = nullptr;
         unsigned int selectsLoose = mvaCollections.size();
         unsigned int selectsHP = mvaCollections.size();
         if (simRecColl.find(tpr) != simRecColl.end()) {
           auto const& rt = simRecColl[tpr];
           if (!rt.empty()) {
             ats++;  //This counter counts the number of simTracks that have a recoTrack associated
             // isRecoMatched = true; // UNUSED
             matchedTrackPointer = rt.begin()->first.get();
             if (rt.size() >= 2) {
               matchedSecondTrackPointer = (rt.begin() + 1)->first.get();
             }
             LogTrace("TrackValidator") << "TrackingParticle #" << st << " with pt=" << sqrt(momentumTP.perp2())
                                        << " associated with quality:" << rt.begin()->second << "\n";
 
             if (doMVAPlots_) {
               // for each MVA we need to take the value of the track
               // with largest MVA value (for the cumulative histograms)
               //
               // also identify the first MVA that possibly selects any
               // track matched to this TrackingParticle, separately
               // for loose and highPurity qualities
               for (size_t imva = 0; imva < mvaCollections.size(); ++imva) {
                 const auto& mva = *(mvaCollections[imva]);
                 const auto& qual = *(qualityMaskCollections[imva]);
 
                 auto iMatch = rt.begin();
                 float maxMva = mva[iMatch->first.key()];
                 for (; iMatch != rt.end(); ++iMatch) {
                   auto itrk = iMatch->first.key();
                   maxMva = std::max(maxMva, mva[itrk]);
 
                   if (selectsLoose >= imva && trackSelected(qual[itrk], reco::TrackBase::loose))
                     selectsLoose = imva;
                   if (selectsHP >= imva && trackSelected(qual[itrk], reco::TrackBase::highPurity))
                     selectsHP = imva;
                 }
                 mvaValues.push_back(maxMva);
               }
             }
           }
         } else {
           LogTrace("TrackValidator") << "TrackingParticle #" << st << " with pt,eta,phi: " << sqrt(momentumTP.perp2())
                                      << " , " << momentumTP.eta() << " , " << momentumTP.phi() << " , "
                                      << " NOT associated to any reco::Track"
                                      << "\n";
         }
 
         int nSimHits = tp.numberOfTrackerHits();
         int nSimLayers = nLayers_tPCeff[tpr];
         int nSimPixelLayers = nPixelLayers_tPCeff[tpr];
         int nSimStripMonoAndStereoLayers = nStripMonoAndStereoLayers_tPCeff[tpr];
         histoProducerAlgo_->fill_recoAssociated_simTrack_histos(histograms.histoProducerAlgo,
                                                                 w,
                                                                 tp,
                                                                 momentumTP,
                                                                 vertexTP,
                                                                 dxySim,
                                                                 dzSim,
                                                                 dxyPVSim,
                                                                 dzPVSim,
                                                                 nSimHits,
                                                                 nSimLayers,
                                                                 nSimPixelLayers,
                                                                 nSimStripMonoAndStereoLayers,
                                                                 matchedTrackPointer,
                                                                 puinfo.getPU_NumInteractions(),
                                                                 dR,
                                                                 dR_jet,
                                                                 thePVposition,
                                                                 theSimPVPosition,
                                                                 bs.position(),
                                                                 mvaValues,
                                                                 selectsLoose,
                                                                 selectsHP);
         mvaValues.clear();
 
         if (matchedTrackPointer && matchedSecondTrackPointer) {
           histoProducerAlgo_->fill_duplicate_histos(
               histograms.histoProducerAlgo, w, *matchedTrackPointer, *matchedSecondTrackPointer);
         }
 
         if (doSummaryPlots_) {
           if (dRtpSelector(tp)) {
             histograms.h_simul_coll[ww]->Fill(www);
             if (matchedTrackPointer) {
               histograms.h_assoc_coll[ww]->Fill(www);
             }
           }
         }
 
       }  // End  for (TrackingParticleCollection::size_type i=0; i<tPCeff.size(); i++){
 
       // ##############################################
       // fill recoTracks histograms (LOOP OVER TRACKS)
       // ##############################################
       if (!doRecoTrackPlots_)
         continue;
       LogTrace("TrackValidator") << "\n# of reco::Tracks with " << label[www].process() << ":" << label[www].label()
                                  << ":" << label[www].instance() << ": " << trackCollection.size() << "\n";
 
       int sat(0);  //This counter counts the number of recoTracks that are associated to SimTracks from Signal only
       int at(0);   //This counter counts the number of recoTracks that are associated to SimTracks
       int rT(0);   //This counter counts the number of recoTracks in general
       int seed_fit_failed = 0;
       size_t n_selTrack_dr = 0;
 
       //calculate dR for tracks
       declareDynArray(float, trackCollection.size(), dR_trk);
       declareDynArray(float, trackCollection.size(), dR_trk_jet);
 #ifndef NO_TRACK_DR
       // this accounts for most of the time spent in MTV and it is used to fill just one histo that is of doubtful usefulness but (maybe) for the whole collection
       const edm::View<Track>* trackCollectionDr = &trackCollection;
       if (calculateDrSingleCollection_) {
         trackCollectionDr = trackCollectionForDrCalculation.product();
       }
       trackDR(trackCollection, *trackCollectionDr, dR_trk, dR_trk_jet, coresVector);
 #endif
 
       for (View<Track>::size_type i = 0; i < trackCollection.size(); ++i) {
         auto track = trackCollection.refAt(i);
         rT++;
         if (trackFromSeedFitFailed(*track))
           ++seed_fit_failed;
         if ((*dRTrackSelector)(*track, bs.position()))
           ++n_selTrack_dr;
 
         bool isSigSimMatched(false);
         bool isSimMatched(false);
         bool isChargeMatched(true);
         int numAssocRecoTracks = 0;
         int nSimHits = 0;
         double sharedFraction = 0.;
 
         auto tpFound = recSimColl.find(track);
         isSimMatched = tpFound != recSimColl.end();
         if (isSimMatched) {
           const auto& tp = tpFound->val;
           nSimHits = tp[0].first->numberOfTrackerHits();
           sharedFraction = tp[0].second;
           if (tp[0].first->charge() != track->charge())
             isChargeMatched = false;
           if (simRecColl.find(tp[0].first) != simRecColl.end())
             numAssocRecoTracks = simRecColl[tp[0].first].size();
           at++;
           for (unsigned int tp_ite = 0; tp_ite < tp.size(); ++tp_ite) {
             TrackingParticle trackpart = *(tp[tp_ite].first);
             if ((trackpart.eventId().event() == 0) && (trackpart.eventId().bunchCrossing() == 0)) {
               isSigSimMatched = true;
               sat++;
               break;
             }
           }
           LogTrace("TrackValidator") << "reco::Track #" << rT << " with pt=" << track->pt()
                                      << " associated with quality:" << tp.begin()->second << "\n";
         } else {
           LogTrace("TrackValidator") << "reco::Track #" << rT << " with pt=" << track->pt()
                                      << " NOT associated to any TrackingParticle"
                                      << "\n";
         }
 
         // set MVA values for this track
         // take also the indices of first MVAs to select by loose and
         // HP quality
         unsigned int selectsLoose = mvaCollections.size();
         unsigned int selectsHP = mvaCollections.size();
         if (doMVAPlots_) {
           for (size_t imva = 0; imva < mvaCollections.size(); ++imva) {
             const auto& mva = *(mvaCollections[imva]);
             const auto& qual = *(qualityMaskCollections[imva]);
             mvaValues.push_back(mva[i]);
 
             if (selectsLoose >= imva && trackSelected(qual[i], reco::TrackBase::loose))
               selectsLoose = imva;
             if (selectsHP >= imva && trackSelected(qual[i], reco::TrackBase::highPurity))
               selectsHP = imva;
           }
         }
 
         double dR = dR_trk[i];
         double dR_jet = dR_trk_jet[i];
         histoProducerAlgo_->fill_generic_recoTrack_histos(histograms.histoProducerAlgo,
                                                           w,
                                                           *track,
                                                           ttopo,
                                                           bs.position(),
                                                           thePVposition,
                                                           theSimPVPosition,
                                                           isSimMatched,
                                                           isSigSimMatched,
                                                           isChargeMatched,
                                                           numAssocRecoTracks,
                                                           puinfo.getPU_NumInteractions(),
                                                           nSimHits,
                                                           sharedFraction,
                                                           dR,
                                                           dR_jet,
                                                           mvaValues,
                                                           selectsLoose,
                                                           selectsHP);
         mvaValues.clear();
 
         if (doSummaryPlots_) {
           histograms.h_reco_coll[ww]->Fill(www);
           if (isSimMatched) {
             histograms.h_assoc2_coll[ww]->Fill(www);
             if (numAssocRecoTracks > 1) {
               histograms.h_looper_coll[ww]->Fill(www);
             }
             if (!isSigSimMatched) {
               histograms.h_pileup_coll[ww]->Fill(www);
             }
           }
         }
 
         // dE/dx
         if (dodEdxPlots_)
           histoProducerAlgo_->fill_dedx_recoTrack_histos(histograms.histoProducerAlgo, w, track, v_dEdx);
 
         //Fill other histos
         if (!isSimMatched)
           continue;
 
         histoProducerAlgo_->fill_simAssociated_recoTrack_histos(histograms.histoProducerAlgo, w, *track);
 
         /* TO BE FIXED LATER
     if (associators[ww]=="trackAssociatorByChi2"){
       //association chi2
       double assocChi2 = -tp.begin()->second;//in association map is stored -chi2
       h_assochi2[www]->Fill(assocChi2);
       h_assochi2_prob[www]->Fill(TMath::Prob((assocChi2)*5,5));
     }
     else if (associators[ww]=="quickTrackAssociatorByHits"){
       double fraction = tp.begin()->second;
       h_assocFraction[www]->Fill(fraction);
       h_assocSharedHit[www]->Fill(fraction*track->numberOfValidHits());
     }
     */
 
         if (doResolutionPlots_[www]) {
           //Get tracking particle parameters at point of closest approach to the beamline
           TrackingParticleRef tpr = tpFound->val.begin()->first;
           TrackingParticle::Vector momentumTP = parametersDefinerTP_->momentum(event, setup, tpr);
           TrackingParticle::Point vertexTP = parametersDefinerTP_->vertex(event, setup, tpr);
           int chargeTP = tpr->charge();
 
           histoProducerAlgo_->fill_ResoAndPull_recoTrack_histos(
               histograms.histoProducerAlgo, w, momentumTP, vertexTP, chargeTP, *track, bs.position());
         }
 
         //TO BE FIXED
         //std::vector<PSimHit> simhits=tpr.get()->trackPSimHit(DetId::Tracker);
         //nrecHit_vs_nsimHit_rec2sim[w]->Fill(track->numberOfValidHits(), (int)(simhits.end()-simhits.begin() ));
 
       }  // End of for(View<Track>::size_type i=0; i<trackCollection.size(); ++i){
       mvaCollections.clear();
       qualityMaskCollections.clear();
 
       histoProducerAlgo_->fill_trackBased_histos(histograms.histoProducerAlgo, w, at, rT, n_selTrack_dr, n_selTP_dr);
       // Fill seed-specific histograms
       if (doSeedPlots_) {
         histoProducerAlgo_->fill_seed_histos(
             histograms.histoProducerAlgo, www, seed_fit_failed, trackCollection.size());
       }
 
       LogTrace("TrackValidator") << "Collection " << www << "\n"
                                  << "Total Simulated (selected): " << n_selTP_dr << "\n"
                                  << "Total Reconstructed (selected): " << n_selTrack_dr << "\n"
                                  << "Total Reconstructed: " << rT << "\n"
                                  << "Total Associated (recoToSim): " << at << "\n"
                                  << "Total Fakes: " << rT - at << "\n";
     }  // End of  for (unsigned int www=0;www<label.size();www++){
   }    //END of for (unsigned int ww=0;ww<associators.size();ww++){
 }

const reco::Vertex::Point * MultiTrackValidator::getRecoPVPosition	(	const edm::Event &	event,
		const edm::Handle< TrackingVertexCollection > &	htv
	)		const

private

Definition at line 396 of file MultiTrackValidator.cc.

References EncodedEventId::bunchCrossing(), EncodedEventId::event(), TrackingVertex::eventId(), recoVertexToken_, and vertexAssociatorToken_.

Referenced by dqmAnalyze().

                                                                                  {
   edm::Handle<edm::View<reco::Vertex>> hvertex;
   event.getByToken(recoVertexToken_, hvertex);
 
   edm::Handle<reco::VertexToTrackingVertexAssociator> hvassociator;
   event.getByToken(vertexAssociatorToken_, hvassociator);
 
   auto v_r2s = hvassociator->associateRecoToSim(hvertex, htv);
   auto pvPtr = hvertex->refAt(0);
   if (pvPtr->isFake() || pvPtr->ndof() < 0)  // skip junk vertices
     return nullptr;
 
   auto pvFound = v_r2s.find(pvPtr);
   if (pvFound == v_r2s.end())
     return nullptr;
 
   for (const auto& vertexRefQuality : pvFound->val) {
     const TrackingVertex& tv = *(vertexRefQuality.first);
     if (tv.eventId().event() == 0 && tv.eventId().bunchCrossing() == 0) {
       return &(pvPtr->position());
     }
   }
 
   return nullptr;
 }

const TrackingVertex::LorentzVector * MultiTrackValidator::getSimPVPosition ( const edm::Handle< TrackingVertexCollection > & htv ) const

private

Definition at line 384 of file MultiTrackValidator.cc.

Referenced by dqmAnalyze().

                                                           {
   for (const auto& simV : *htv) {
     if (simV.eventId().bunchCrossing() != 0)
       continue;  // remove OOTPU
     if (simV.eventId().event() != 0)
       continue;  // pick the PV of hard scatter
     return &(simV.position());
   }
   return nullptr;
 }

size_t MultiTrackValidator::tpDR	(	const TrackingParticleRefVector &	tPCeff,
		const std::vector< size_t > &	selected_tPCeff,
		DynArray< float > &	dR_tPCeff,
		DynArray< float > &	dR_tPCeff_jet,
		const edm::View< reco::Candidate > *	cores
	)		const

private

Definition at line 481 of file MultiTrackValidator.cc.

References reco::deltaR2(), dRtpSelector, edm::RefVector< C, T, F >::empty(), PVValHelper::eta, reco::Candidate::eta(), metsig::jet, SiStripPI::max, AlCaHLTBitMon_ParallelJobs::p, phi, reco::Candidate::phi(), edm::RefVector< C, T, F >::size(), edm::View< T >::size(), mathSSE::sqrt(), and cmsswSequenceInfo::tp.

Referenced by dqmAnalyze().

                                                                             {
   if (tPCeff.empty()) {
     return 0;
   }
   float etaL[tPCeff.size()], phiL[tPCeff.size()];
   size_t n_selTP_dr = 0;
   for (size_t iTP : selected_tPCeff) {
     //calculare dR wrt inclusive collection (also with PU, low pT, displaced)
     auto const& tp2 = *(tPCeff[iTP]);
     auto&& p = tp2.momentum();
     etaL[iTP] = etaFromXYZ(p.x(), p.y(), p.z());
     phiL[iTP] = atan2f(p.y(), p.x());
   }
   for (size_t iTP1 : selected_tPCeff) {
     auto const& tp = *(tPCeff[iTP1]);
     double dR = std::numeric_limits<double>::max();
     double dR_jet = std::numeric_limits<double>::max();
     if (dRtpSelector(tp)) {  //only for those needed for efficiency!
       ++n_selTP_dr;
       float eta = etaL[iTP1];
       float phi = phiL[iTP1];
       for (size_t iTP2 : selected_tPCeff) {
         //calculare dR wrt inclusive collection (also with PU, low pT, displaced)
         if (iTP1 == iTP2) {
           continue;
         }
         auto dR_tmp = reco::deltaR2(eta, phi, etaL[iTP2], phiL[iTP2]);
         if (dR_tmp < dR)
           dR = dR_tmp;
       }  // ttp2 (iTP)
       if (cores != nullptr) {
         for (unsigned int ji = 0; ji < cores->size(); ji++) {  //jet loop
           const reco::Candidate& jet = (*cores)[ji];
           double jet_eta = jet.eta();
           double jet_phi = jet.phi();
           auto dR_jet_tmp = reco::deltaR2(eta, phi, jet_eta, jet_phi);
           if (dR_jet_tmp < dR_jet)
             dR_jet = dR_jet_tmp;
         }
       }
     }
     dR_tPCeff[iTP1] = std::sqrt(dR);
     dR_tPCeff_jet[iTP1] = std::sqrt(dR_jet);
 
   }  // tp
   return n_selTP_dr;
 }

void MultiTrackValidator::tpParametersAndSelection	(	const Histograms &	histograms,
		const TrackingParticleRefVector &	tPCeff,
		const edm::Event &	event,
		const edm::EventSetup &	setup,
		const reco::BeamSpot &	bs,
		std::vector< std::tuple< TrackingParticle::Vector, TrackingParticle::Point >> &	momVert_tPCeff,
		std::vector< size_t > &	selected_tPCeff
	)		const

private

Definition at line 423 of file MultiTrackValidator.cc.

References EncodedEventId::bunchCrossing(), cosmictpSelector, doSimPlots_, TrackingParticle::eventId(), MultiTrackValidatorHistograms::histoProducerAlgo, histoProducerAlgo_, dqmiolumiharvest::j, TrackingParticle::momentum(), parametersDefinerIsCosmic_, parametersDefinerTP_, edm::RefVector< C, T, F >::size(), cmsswSequenceInfo::tp, tpSelector, and TrackingParticle::vertex().

Referenced by dqmAnalyze().

                                               {
   selected_tPCeff.reserve(tPCeff.size());
   momVert_tPCeff.reserve(tPCeff.size());
   int nIntimeTPs = 0;
   if (parametersDefinerIsCosmic_) {
     for (size_t j = 0; j < tPCeff.size(); ++j) {
       const TrackingParticleRef& tpr = tPCeff[j];
 
       auto const& rec = parametersDefinerTP_->momentumAndVertex(event, setup, tpr);
       TrackingParticle::Vector const& momentum = std::get<0>(rec);
       TrackingParticle::Point const& vertex = std::get<1>(rec);
       if (doSimPlots_) {
         histoProducerAlgo_->fill_generic_simTrack_histos(
             histograms.histoProducerAlgo, momentum, vertex, tpr->eventId().bunchCrossing());
       }
       if (tpr->eventId().bunchCrossing() == 0)
         ++nIntimeTPs;
 
       if (cosmictpSelector(tpr, &bs, event, setup)) {
         selected_tPCeff.push_back(j);
         momVert_tPCeff.emplace_back(momentum, vertex);
       }
     }
   } else {
     size_t j = 0;
     for (auto const& tpr : tPCeff) {
       const TrackingParticle& tp = *tpr;
 
       // TODO: do we want to fill these from all TPs that include IT
       // and OOT (as below), or limit to IT+OOT TPs passing tpSelector
       // (as it was before)? The latter would require another instance
       // of tpSelector with intimeOnly=False.
       if (doSimPlots_) {
         histoProducerAlgo_->fill_generic_simTrack_histos(
             histograms.histoProducerAlgo, tp.momentum(), tp.vertex(), tp.eventId().bunchCrossing());
       }
       if (tp.eventId().bunchCrossing() == 0)
         ++nIntimeTPs;
 
       if (tpSelector(tp)) {
         selected_tPCeff.push_back(j);
         momVert_tPCeff.emplace_back(parametersDefinerTP_->momentumAndVertex(event, setup, tpr));
       }
       ++j;
     }
   }
   if (doSimPlots_) {
     histoProducerAlgo_->fill_simTrackBased_histos(histograms.histoProducerAlgo, nIntimeTPs);
   }
 }

void MultiTrackValidator::trackDR	(	const edm::View< reco::Track > &	trackCollection,
		const edm::View< reco::Track > &	trackCollectionDr,
		DynArray< float > &	dR_trk,
		DynArray< float > &	dR_trk_jet,
		const edm::View< reco::Candidate > *	cores
	)		const

private

Definition at line 533 of file MultiTrackValidator.cc.

References reco::deltaR2(), edm::View< T >::empty(), PVValHelper::eta, reco::Candidate::eta(), mps_fire::i, dqmiolumiharvest::j, metsig::jet, SiStripPI::max, min(), AlCaHLTBitMon_ParallelJobs::p, phi, reco::Candidate::phi(), edm::View< T >::size(), mathSSE::sqrt(), HLT_FULL_cff::track, and trackFromSeedFitFailed().

Referenced by dqmAnalyze().

                                                                              {
   if (trackCollectionDr.empty()) {
     return;
   }
   int i = 0;
   float etaL[trackCollectionDr.size()];
   float phiL[trackCollectionDr.size()];
   bool validL[trackCollectionDr.size()];
   for (auto const& track2 : trackCollectionDr) {
     auto&& p = track2.momentum();
     etaL[i] = etaFromXYZ(p.x(), p.y(), p.z());
     phiL[i] = atan2f(p.y(), p.x());
     validL[i] = !trackFromSeedFitFailed(track2);
     ++i;
   }
   for (View<reco::Track>::size_type i = 0; i < trackCollection.size(); ++i) {
     auto const& track = trackCollection[i];
     auto dR = std::numeric_limits<float>::max();
     auto dR_jet = std::numeric_limits<float>::max();
     if (!trackFromSeedFitFailed(track)) {
       auto&& p = track.momentum();
       float eta = etaFromXYZ(p.x(), p.y(), p.z());
       float phi = atan2f(p.y(), p.x());
       for (View<reco::Track>::size_type j = 0; j < trackCollectionDr.size(); ++j) {
         if (!validL[j])
           continue;
         auto dR_tmp = reco::deltaR2(eta, phi, etaL[j], phiL[j]);
         if ((dR_tmp < dR) & (dR_tmp > std::numeric_limits<float>::min()))
           dR = dR_tmp;
       }
       if (cores != nullptr) {
         for (unsigned int ji = 0; ji < cores->size(); ji++) {  //jet loop
           const reco::Candidate& jet = (*cores)[ji];
           double jet_eta = jet.eta();
           double jet_phi = jet.phi();
           auto dR_jet_tmp = reco::deltaR2(eta, phi, jet_eta, jet_phi);
           if (dR_jet_tmp < dR_jet)
             dR_jet = dR_jet_tmp;
         }
       }
     }
     dR_trk[i] = std::sqrt(dR);
     dR_trk_jet[i] = std::sqrt(dR_jet);
   }
 }