#include <MultiTrackValidator.h>

Inheritance diagram for MultiTrackValidator:

Public Types
using	Histograms = MultiTrackValidatorHistograms

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::ConcurrentBooker &, 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

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

virtual	~EDConsumerBase () noexcept(false)

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

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_

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 ESProduct , Transition Tr = Transition::Event>
auto	esConsumes (eventsetup::EventSetupRecordKey const &, 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 37 of file MultiTrackValidator.h.

Member Typedef Documentation

using MultiTrackValidator::Histograms = MultiTrackValidatorHistograms

Definition at line 39 of file MultiTrackValidator.h.

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

private

Definition at line 121 of file MultiTrackValidator.h.

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

private

Definition at line 122 of file MultiTrackValidator.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 52 of file MultiTrackValidator.cc.

References _simHitTpMapTag, associatormapRtSs, associatormapStRs, associators, associatorTokens, AlignmentProducer_cff::beamSpotTag, begin, bsSrc, calculateDrSingleCollection_, edm::EDConsumerBase::consumes(), 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_, checklumidiff::l, edm::InputTag::label(), label, label_pileupinfo, label_tp_effic, label_tp_effic_refvector, label_tp_fake, label_tp_fake_refvector, label_tv, tablePrinter::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(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"));
 
   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"));
 
   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 215 of file MultiTrackValidator.cc.

215 {}

Member Function Documentation

void MultiTrackValidator::bookHistograms	(	DQMStore::ConcurrentBooker &	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 218 of file MultiTrackValidator.cc.

References patPFMETCorrections_cff::algo, trackingPlots::assoc, associators, DQMStore::ConcurrentBooker::book1D(), 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, mvaQualityCollectionTokens_, edm::InputTag::process(), python.rootplot.root2matplotlib::replace(), edm::IndexSet::reserve(), 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 = [&](ConcurrentMonitorElement 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 492 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(), cosmictpSelector, declareDynArray, dodEdxPlots_, doMVAPlots_, doPlotsOnlyForTruePV_, doPVAssociationPlots_, doRecoTrackPlots_, doResolutionPlots_, doSeedPlots_, doSimTrackPlots_, doSummaryPlots_, PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, dRtpSelector, dRTrackSelector, TrackingParticleIP::dxy(), TrackingParticleIP::dz(), relativeConstraints::empty, edm::AssociationMap< Tag >::end(), EncodedEventId::event(), TrackingParticle::eventId(), Exception, edm::AssociationMap< Tag >::find(), plotBeamSpotDB::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(), 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(), DetachedQuadStep_cff::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(), HiIsolationCommonParameters_cff::track, findElectronsInSiStrips_cfi::trackCollection, TrackValidationHeavyIons_cff::trackCollectionForDrCalculation, trackDR(), trackFromSeedFitFailed(), 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";
 
 
   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);
   size_t n_selTP_dr = tpDR(tPCeff, selected_tPCeff, dR_tPCeff);
 
   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];
 
     //---------- 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, 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);
       trackDR(trackCollection, *trackCollectionDr, dR_trk);
 
       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];
     histoProducerAlgo_->fill_generic_recoTrack_histos(histograms.histoProducerAlgo,w,*track, ttopo, bs.position(), thePVposition, theSimPVPosition, isSimMatched,isSigSimMatched, isChargeMatched, numAssocRecoTracks, puinfo.getPU_NumInteractions(), nSimHits, sharedFraction, dR, 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 347 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 338 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
	)		const

private

Definition at line 431 of file MultiTrackValidator.cc.

References reco::deltaR2(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, dRtpSelector, PVValHelper::eta, SiStripPI::max, AlCaHLTBitMon_ParallelJobs::p, phi, edm::RefVector< C, T, F >::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();
     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)
     }
     dR_tPCeff[iTP1] = std::sqrt(dR);
   }  // 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 373 of file MultiTrackValidator.cc.

References EncodedEventId::bunchCrossing(), cosmictpSelector, doSimPlots_, TrackingParticle::eventId(), MultiTrackValidatorHistograms::histoProducerAlgo, histoProducerAlgo_, ParametersDefinerForTP::momentum(), TrackingParticle::momentum(), parametersDefinerIsCosmic_, edm::RefVector< C, T, F >::size(), tpSelector, 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
	)		const

private

Definition at line 462 of file MultiTrackValidator.cc.

References reco::deltaR2(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, PVValHelper::eta, mps_fire::i, SiStripPI::max, min(), AlCaHLTBitMon_ParallelJobs::p, phi, edm::View< T >::size(), mathSSE::sqrt(), HiIsolationCommonParameters_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();
     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;
       }
     }
     dR_trk[i] = std::sqrt(dR);
   }
 }