#include <MultiTrackValidator.h>

Inheritance diagram for MultiTrackValidator:

Public Types
using	Histograms = MultiTrackValidatorHistograms

Public Types inherited from DQMGlobalEDAnalyzer< MultiTrackValidatorHistograms >
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

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

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 edm::global::EDAnalyzer< edm::RunCache< MultiTrackValidatorHistograms >, Args... >
	EDAnalyzer ()=default

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

	~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

ESProxyIndex const *	esGetTokenIndices (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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

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

std::string	parametersDefiner

const bool	parametersDefinerIsCosmic_

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

const bool	useAssociators_

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 ParametersDefinerForTP &parametersDefinerTP, 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::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)

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<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

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 127 of file MultiTrackValidator.h.

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

private

Definition at line 128 of file MultiTrackValidator.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 49 of file MultiTrackValidator.cc.

References _simHitTpMapTag, associatormapRtSs, associatormapStRs, associators, associatorTokens, AlignmentProducer_cff::beamSpotTag, begin, bsSrc, calculateDrSingleCollection_, edm::EDConsumerBase::consumes(), cores_, cosmictpSelector, dirName_, dodEdxPlots_, doMVAPlots_, doPlotsOnlyForTruePV_, doPVAssociationPlots_, doResolutionPlots_, MultiTrackValidator_cfi::doResolutionPlotsForLabels, doSeedPlots_, dRtpSelector, dRTrackSelector, relativeConstraints::empty, 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, SummaryClient_cfi::labels, edm::EDConsumerBase::labelsForToken(), labelToken, labelTokenForDrCalculation, m_dEdx1Tag, m_dEdx2Tag, MTVHistoProducerAlgoForTracker::makeRecoTrackSelectorFromTPSelectorParameters(), edm::ProductLabels::module, mvaQualityCollectionTokens_, recoVertexToken_, simHitTokens_, TrackRefitter_38T_cff::src, AlCaHLTBitMon_QueryRunRegistry::string, GlobalPosition_Frontier_DevDB_cff::tag, tpNLayersToken_, tpNPixelLayersToken_, tpNStripStereoLayersToken_, tpSelector, useAssociators_, useGsf, edm::vector_transform(), and vertexAssociatorToken_.

     : associators(pset.getUntrackedParameter<std::vector<edm::InputTag>>("associators")),
       label(pset.getParameter<std::vector<edm::InputTag>>("label")),
       parametersDefiner(pset.getParameter<std::string>("parametersDefiner")),
       parametersDefinerIsCosmic_(parametersDefiner == "CosmicParametersDefinerForTP"),
       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);
 
   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 220 of file MultiTrackValidator.cc.

220 {}

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 DQMGlobalEDAnalyzer< MultiTrackValidatorHistograms >.

Definition at line 222 of file MultiTrackValidator.cc.

References ecalcalib_dqm_sourceclient-live_cfg::algo, trackingPlots::assoc, associators, dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), dqm::dqmstoreimpl::DQMStore::IBooker::cd(), TrackerOfflineValidation_Dqm_cff::dirName, dirName_, dodEdxPlots_, doMVAPlots_, doPVAssociationPlots_, doRecoTrackPlots_, doResolutionPlots_, doSeedPlots_, doSimPlots_, doSimTrackPlots_, doSummaryPlots_, edm::RefVector< C, T, F >::empty(), Exception, MultiTrackValidatorHistograms::h_assoc2_coll, MultiTrackValidatorHistograms::h_assoc_coll, MultiTrackValidatorHistograms::h_looper_coll, MultiTrackValidatorHistograms::h_pileup_coll, MultiTrackValidatorHistograms::h_reco_coll, MultiTrackValidatorHistograms::h_simul_coll, edm::IndexSet::has(), MultiTrackValidatorHistograms::histoProducerAlgo, histoProducerAlgo_, mps_fire::i, edm::RefVector< C, T, F >::id(), edm::IndexSet::insert(), edm::InputTag::instance(), edm::InputTag::label(), label, hlt_dqm_clientPB-live_cfg::me, mvaQualityCollectionTokens_, edm::InputTag::process(), python.rootplot.root2matplotlib::replace(), edm::IndexSet::reserve(), dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), and edm::RefVector< C, T, F >::size().

                                                                        {
   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

override

Method called once per event.

Definition at line 572 of file MultiTrackValidator.cc.

References _simHitTpMapTag, funct::abs(), reco::TrackToTrackingParticleAssociator::associateRecoToSim(), reco::TrackToTrackingParticleAssociator::associateSimToReco(), associationMapFilterValues(), associatormapRtSs, associatormapStRs, associators, associatorTokens, bsSrc, EncodedEventId::bunchCrossing(), calculateDrSingleCollection_, ParametersDefinerForTP::clone(), HLT_FULL_cff::cores, cores_, cosmictpSelector, declareDynArray, dodEdxPlots_, doMVAPlots_, doPlotsOnlyForTruePV_, doPVAssociationPlots_, doRecoTrackPlots_, doResolutionPlots_, doSeedPlots_, doSimTrackPlots_, doSummaryPlots_, HGC3DClusterGenMatchSelector_cfi::dR, dRtpSelector, dRTrackSelector, TrackingParticleIP::dxy(), TrackingParticleIP::dz(), relativeConstraints::empty, edm::AssociationMap< Tag >::end(), EncodedEventId::event(), TrackingParticle::eventId(), Exception, edm::AssociationMap< Tag >::find(), dqmdumpme::first, edm::EventSetup::get(), edm::Event::getByToken(), PileupSummaryInfo::getPU_NumInteractions(), 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(), parametersDefiner, parametersDefinerIsCosmic_, reco::BeamSpot::position(), edm::Handle< T >::product(), edm::RefToBaseVector< T >::push_back(), edm::RefVector< C, T, F >::push_back(), edm::View< T >::refAt(), simPVMaxZ_, edm::RefVector< C, T, F >::size(), edm::View< T >::size(), findQualityFiles::size, mathSSE::sqrt(), tpDR(), tpNLayersToken_, tpNPixelLayersToken_, tpNStripStereoLayersToken_, tpParametersAndSelection(), HLT_2018_cff::track, duplicaterechits_cfi::trackCollection, TrackValidationHeavyIons_cff::trackCollectionForDrCalculation, trackDR(), trackFromSeedFitFailed(), useAssociators_, bphysicsOniaDQM_cfi::vertex, 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";
 
   edm::ESHandle<ParametersDefinerForTP> parametersDefinerTPHandle;
   setup.get<TrackAssociatorRecord>().get(parametersDefiner, parametersDefinerTPHandle);
   //Since we modify the object, we must clone it
   auto parametersDefinerTP = parametersDefinerTPHandle->clone();
 
   edm::ESHandle<TrackerTopology> httopo;
   setup.get<TrackerTopologyRcd>().get(httopo);
   const TrackerTopology& ttopo = *httopo;
 
   // 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);
     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);
     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;
 
   ensureEffIsSubsetOfFake(tPCeff, tPCfake);
 
   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, *parametersDefinerTP, 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;
         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
       const edm::View<Track>* trackCollectionDr = &trackCollection;
       if (calculateDrSingleCollection_) {
         trackCollectionDr = trackCollectionForDrCalculation.product();
       }
       declareDynArray(float, trackCollection.size(), dR_trk);
       declareDynArray(float, trackCollection.size(), dR_trk_jet);
       trackDR(trackCollection, *trackCollectionDr, dR_trk, dR_trk_jet, coresVector);
 
       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 389 of file MultiTrackValidator.cc.

References reco::VertexToTrackingVertexAssociator::associateRecoToSim(), 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 377 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 476 of file MultiTrackValidator.cc.

References reco::deltaR2(), HGC3DClusterGenMatchSelector_cfi::dR, dRtpSelector, 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(), and mathSSE::sqrt().

Referenced by dqmAnalyze().

                                                                             {
   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 ParametersDefinerForTP &	parametersDefinerTP,
		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 416 of file MultiTrackValidator.cc.

References EncodedEventId::bunchCrossing(), cosmictpSelector, doSimPlots_, TrackingParticle::eventId(), MultiTrackValidatorHistograms::histoProducerAlgo, histoProducerAlgo_, dqmiolumiharvest::j, ParametersDefinerForTP::momentum(), TrackingParticle::momentum(), parametersDefinerIsCosmic_, edm::RefVector< C, T, F >::size(), tpSelector, bphysicsOniaDQM_cfi::vertex, ParametersDefinerForTP::vertex(), 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];
 
       TrackingParticle::Vector momentum = parametersDefinerTP.momentum(event, setup, tpr);
       TrackingParticle::Point vertex = parametersDefinerTP.vertex(event, setup, tpr);
       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);
         TrackingParticle::Vector momentum = parametersDefinerTP.momentum(event, setup, tpr);
         TrackingParticle::Point vertex = parametersDefinerTP.vertex(event, setup, tpr);
         momVert_tPCeff.emplace_back(momentum, vertex);
       }
       ++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 525 of file MultiTrackValidator.cc.

References reco::deltaR2(), HGC3DClusterGenMatchSelector_cfi::dR, 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_2018_cff::track, and trackFromSeedFitFailed().

Referenced by dqmAnalyze().

                                                                              {
   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);
   }
 }