Primary4DVertexValidation Class Reference

Inheritance diagram for Primary4DVertexValidation:

Classes
struct	recoPrimaryVertex
struct	simPrimaryVertex

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &) override
	Primary4DVertexValidation (const edm::ParameterSet &)
	~Primary4DVertexValidation () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Private Types
typedef math::XYZTLorentzVector	LorentzVector

Private Member Functions
std::pair< const edm::Ref< std::vector< TrackingParticle > > *, int >	getMatchedTP (const reco::TrackBaseRef &, const TrackingVertexRef &)
std::vector< Primary4DVertexValidation::recoPrimaryVertex >	getRecoPVs (const edm::Handle< edm::View< reco::Vertex >> &)
std::vector< Primary4DVertexValidation::simPrimaryVertex >	getSimPVs (const edm::Handle< TrackingVertexCollection > &)
void	getWosWnt (const reco::Vertex &, const reco::TrackBaseRef &, const edm::ValueMap< float > &, const edm::Handle< reco::BeamSpot > &, double &, double &)
void	isParticle (const reco::TrackBaseRef &, const edm::ValueMap< float > &, const edm::ValueMap< float > &, const edm::ValueMap< float > &, unsigned int &, bool &, bool &, bool &, bool &)
void	matchReco2Sim (std::vector< recoPrimaryVertex > &, std::vector< simPrimaryVertex > &, const edm::ValueMap< float > &, const edm::ValueMap< float > &, const edm::Handle< reco::BeamSpot > &)
bool	matchRecoTrack2SimSignal (const reco::TrackBaseRef &)
const bool	mvaRecSel (const reco::TrackBase &, const reco::Vertex &, const double &, const double &)
const bool	mvaTPSel (const TrackingParticle &)
void	observablesFromJets (const std::vector< reco::Track > &, const std::vector< double > &, const std::vector< int > &, const std::string &, unsigned int &, double &, double &, double &, double &)
void	printMatchedRecoTrackInfo (const reco::Vertex &, const reco::TrackBaseRef &, const TrackingParticleRef &, const unsigned int &)
void	printSimVtxRecoVtxInfo (const struct Primary4DVertexValidation::simPrimaryVertex &, const struct Primary4DVertexValidation::recoPrimaryVertex &)
bool	select (const reco::Vertex &, int level=0)
double	timeFromTrueMass (double, double, double, double)

Private Attributes
bool	debug_
const std::string	folder_
const std::vector< double >	lineDensityPar_
MonitorElement *	meBarrelNoPIDtype_
MonitorElement *	meBarrelPIDp_
MonitorElement *	meBarrelTrueKAsK_
MonitorElement *	meBarrelTrueKAsP_
MonitorElement *	meBarrelTrueKAsPi_
MonitorElement *	meBarrelTrueKNoPID_
MonitorElement *	meBarrelTruePAsK_
MonitorElement *	meBarrelTruePAsP_
MonitorElement *	meBarrelTruePAsPi_
MonitorElement *	meBarrelTruePiAsK_
MonitorElement *	meBarrelTruePiAsP_
MonitorElement *	meBarrelTruePiAsPi_
MonitorElement *	meBarrelTruePiNoPID_
MonitorElement *	meBarrelTruePNoPID_
MonitorElement *	meDeltaTfakefake_
MonitorElement *	meDeltaTfakereal_
MonitorElement *	meDeltaTrealreal_
MonitorElement *	meDeltaZfakefake_
MonitorElement *	meDeltaZfakereal_
MonitorElement *	meDeltaZrealreal_
MonitorElement *	meEndcapNoPIDtype_
MonitorElement *	meEndcapPIDp_
MonitorElement *	meEndcapTrueKAsK_
MonitorElement *	meEndcapTrueKAsP_
MonitorElement *	meEndcapTrueKAsPi_
MonitorElement *	meEndcapTrueKNoPID_
MonitorElement *	meEndcapTruePAsK_
MonitorElement *	meEndcapTruePAsP_
MonitorElement *	meEndcapTruePAsPi_
MonitorElement *	meEndcapTruePiAsK_
MonitorElement *	meEndcapTruePiAsP_
MonitorElement *	meEndcapTruePiAsPi_
MonitorElement *	meEndcapTruePiNoPID_
MonitorElement *	meEndcapTruePNoPID_
MonitorElement *	meFakeTrackRecLVRelMult_
MonitorElement *	meFakeTrackRecLVRelMultvsMult_
MonitorElement *	meFakeTrackRecLVRelSumWos_
MonitorElement *	meFakeTrackRecLVRelSumWosvsSumWos_
MonitorElement *	meFakeTrackRelMult_
MonitorElement *	meFakeTrackRelMultvsMult_
MonitorElement *	meFakeTrackRelSumWos_
MonitorElement *	meFakeTrackRelSumWosvsSumWos_
MonitorElement *	meJetsFakeRelSumPt2_
MonitorElement *	meJetsFakeRelSumPt2vsSumPt2_
MonitorElement *	meJetsPUHt_
MonitorElement *	meJetsPUMetPt_
MonitorElement *	meJetsPUMult_
MonitorElement *	meJetsPURelHt_
MonitorElement *	meJetsPURelHtvsHt_
MonitorElement *	meJetsPURelMetPt_
MonitorElement *	meJetsPURelMetPtvsMetPt_
MonitorElement *	meJetsPURelMult_
MonitorElement *	meJetsPURelMultvsMult_
MonitorElement *	meJetsPURelSumPt2_
MonitorElement *	meJetsPURelSumPt2vsSumPt2_
MonitorElement *	meJetsPURelSumPz_
MonitorElement *	meJetsPURelSumPzvsSumPz_
MonitorElement *	meJetsPUSumPt2_
MonitorElement *	meJetsPUSumPz_
MonitorElement *	meJetsRecLVFakeRelSumPt2_
MonitorElement *	meJetsRecLVFakeRelSumPt2vsSumPt2_
MonitorElement *	meJetsRecLVPUHt_
MonitorElement *	meJetsRecLVPUMetPt_
MonitorElement *	meJetsRecLVPUMult_
MonitorElement *	meJetsRecLVPURelHt_
MonitorElement *	meJetsRecLVPURelHtvsHt_
MonitorElement *	meJetsRecLVPURelMetPt_
MonitorElement *	meJetsRecLVPURelMetPtvsMetPt_
MonitorElement *	meJetsRecLVPURelMult_
MonitorElement *	meJetsRecLVPURelMultvsMult_
MonitorElement *	meJetsRecLVPURelSumPt2_
MonitorElement *	meJetsRecLVPURelSumPt2vsSumPt2_
MonitorElement *	meJetsRecLVPURelSumPz_
MonitorElement *	meJetsRecLVPURelSumPzvsSumPz_
MonitorElement *	meJetsRecLVPUSumPt2_
MonitorElement *	meJetsRecLVPUSumPz_
MonitorElement *	meMatchQual_
MonitorElement *	meMVATrackEffEtaTot_
MonitorElement *	meMVATrackEffPtTot_
MonitorElement *	meMVATrackMatchedEffEtaMtd_
MonitorElement *	meMVATrackMatchedEffEtaTot_
MonitorElement *	meMVATrackMatchedEffPtMtd_
MonitorElement *	meMVATrackMatchedEffPtTot_
MonitorElement *	meMVATrackPullTot_
MonitorElement *	meMVATrackResTot_
MonitorElement *	meMVATrackZposResTot_
MonitorElement *	mePUTrackMult_
MonitorElement *	mePUTrackRecLVMult_
MonitorElement *	mePUTrackRecLVRelMult_
MonitorElement *	mePUTrackRecLVRelMultvsMult_
MonitorElement *	mePUTrackRecLVRelSumPt2_
MonitorElement *	mePUTrackRecLVRelSumPt2vsSumPt2_
MonitorElement *	mePUTrackRecLVRelSumPt_
MonitorElement *	mePUTrackRecLVRelSumPtvsSumPt_
MonitorElement *	mePUTrackRecLVRelSumWnt_
MonitorElement *	mePUTrackRecLVRelSumWntvsSumWnt_
MonitorElement *	mePUTrackRecLVRelSumWos_
MonitorElement *	mePUTrackRecLVRelSumWosvsSumWos_
MonitorElement *	mePUTrackRecLVSumPt2_
MonitorElement *	mePUTrackRecLVSumPt_
MonitorElement *	mePUTrackRecLVSumWnt_
MonitorElement *	mePUTrackRecLVSumWos_
MonitorElement *	mePUTrackRecLVWnt_
MonitorElement *	mePUTrackRelMult_
MonitorElement *	mePUTrackRelMultvsMult_
MonitorElement *	mePUTrackRelSumPt2_
MonitorElement *	mePUTrackRelSumPt2vsSumPt2_
MonitorElement *	mePUTrackRelSumPt_
MonitorElement *	mePUTrackRelSumPtvsSumPt_
MonitorElement *	mePUTrackRelSumWnt_
MonitorElement *	mePUTrackRelSumWntvsSumWnt_
MonitorElement *	mePUTrackRelSumWos_
MonitorElement *	mePUTrackRelSumWosvsSumWos_
MonitorElement *	mePUTrackSumPt2_
MonitorElement *	mePUTrackSumPt_
MonitorElement *	mePUTrackSumWnt_
MonitorElement *	mePUTrackSumWos_
MonitorElement *	mePUTrackWnt_
MonitorElement *	mePUvsOtherFakeV_
MonitorElement *	mePUvsRealV_
MonitorElement *	mePUvsSplitV_
MonitorElement *	meRecoPosInRecoOrigCollection_
MonitorElement *	meRecoPosInSimCollection_
MonitorElement *	meRecoPVPosSignal_
MonitorElement *	meRecoPVPosSignalNotHighestPt_
MonitorElement *	meRecoVtxVsLineDensity_
MonitorElement *	meRecPVT_
MonitorElement *	meRecPVZ_
MonitorElement *	meRecVerNumber_
MonitorElement *	meSecTrackRecLVRelSumWos_
MonitorElement *	meSecTrackRecLVRelSumWosvsSumWos_
MonitorElement *	meSecTrackRecLVSumWos_
MonitorElement *	meSecTrackRelSumWos_
MonitorElement *	meSecTrackRelSumWosvsSumWos_
MonitorElement *	meSecTrackSumWos_
MonitorElement *	meSimPosInSimOrigCollection_
MonitorElement *	meSimPVZ_
MonitorElement *	meTimePull_
MonitorElement *	meTimeRes_
MonitorElement *	meTimeSignalPull_
MonitorElement *	meTimeSignalRes_
MonitorElement *	meTrack3DposRes_ [3]
MonitorElement *	meTrackPull_ [3]
MonitorElement *	meTrackPullHighP_ [3]
MonitorElement *	meTrackPullHighPTot_
MonitorElement *	meTrackPullLowP_ [3]
MonitorElement *	meTrackPullLowPTot_
MonitorElement *	meTrackPullTot_
MonitorElement *	meTrackRes_ [3]
MonitorElement *	meTrackResHighP_ [3]
MonitorElement *	meTrackResHighPTot_
MonitorElement *	meTrackResLowP_ [3]
MonitorElement *	meTrackResLowPTot_
MonitorElement *	meTrackResMass_ [3]
MonitorElement *	meTrackResMassPions_ [3]
MonitorElement *	meTrackResMassProtons_ [3]
MonitorElement *	meTrackResMassTrue_ [3]
MonitorElement *	meTrackResMassTruePions_ [3]
MonitorElement *	meTrackResMassTrueProtons_ [3]
MonitorElement *	meTrackResTot_
MonitorElement *	meTrackZposRes_ [3]
MonitorElement *	meTrackZposResTot_
const double	minProbHeavy_
edm::EDGetTokenT< edm::ValueMap< float > >	momentumToken_
const double	mvaTh_
bool	optionalPlots_
edm::EDGetTokenT< edm::ValueMap< float > >	pathLengthToken_
edm::ESGetToken< HepPDT::ParticleDataTable, edm::DefaultRecord >	pdtToken_
edm::EDGetTokenT< edm::ValueMap< float > >	probKToken_
edm::EDGetTokenT< edm::ValueMap< float > >	probPiToken_
edm::EDGetTokenT< edm::ValueMap< float > >	probPToken_
const reco::RecoToSimCollection *	r2s_
edm::EDGetTokenT< edm::View< reco::Vertex > >	Rec4DVerToken_
edm::EDGetTokenT< reco::BeamSpot >	RecBeamSpotToken_
edm::EDGetTokenT< reco::RecoToSimCollection >	recoToSimAssociationToken_
edm::EDGetTokenT< reco::TrackCollection >	RecTrackToken_
const reco::SimToRecoCollection *	s2r_
edm::EDGetTokenT< edm::ValueMap< float > >	sigmat0SafePidToken_
edm::EDGetTokenT< reco::SimToRecoCollection >	simToRecoAssociationToken_
edm::EDGetTokenT< edm::ValueMap< float > >	t0PidToken_
edm::EDGetTokenT< edm::ValueMap< float > >	t0SafePidToken_
edm::EDGetTokenT< edm::ValueMap< float > >	timeToken_
edm::EDGetTokenT< edm::ValueMap< float > >	tmtdToken_
edm::EDGetTokenT< edm::ValueMap< float > >	tofKToken_
edm::EDGetTokenT< edm::ValueMap< float > >	tofPiToken_
edm::EDGetTokenT< edm::ValueMap< float > >	tofPToken_
edm::EDGetTokenT< edm::ValueMap< int > >	trackAssocToken_
edm::EDGetTokenT< TrackingParticleCollection >	trackingParticleCollectionToken_
edm::EDGetTokenT< TrackingVertexCollection >	trackingVertexCollectionToken_
edm::EDGetTokenT< edm::ValueMap< float > >	trackMVAQualToken_
const double	trackweightTh_
bool	use3dNoTime_
bool	use_only_charged_tracks_
edm::EDGetTokenT< std::vector< PileupSummaryInfo > >	vecPileupSummaryInfoToken_

Static Private Attributes
static constexpr double	c_ = 2.99792458e1
static constexpr float	c_cm_ns = geant_units::operators::convertMmToCm(CLHEP::c_light)
static constexpr double	deltaZcut_ = 0.1
static constexpr double	etacutGEN_ = 4.
static constexpr double	etacutREC_ = 3.
static constexpr double	maxRank_ = 8.
static constexpr double	maxTry_ = 10.
static constexpr double	minThrMetPt_ = 1.e-3
static constexpr double	minThrSumPt2_ = 1.e-3
static constexpr double	minThrSumPt_ = 0.01
static constexpr double	minThrSumPz_ = 1.e-4
static constexpr double	minThrSumWnt_ = 0.01
static constexpr double	minThrSumWos_ = 0.1
static constexpr double	mvaH_ = 0.8
static constexpr double	mvaL_ = 0.5
static constexpr unsigned int	NOT_MATCHED = 66666
static constexpr double	pTcut_ = 0.7
static constexpr double	selNdof_ = 4.
static constexpr double	simUnit_ = 1e9
static constexpr double	tol_ = 1.e-4
static constexpr double	trackMaxBtlEta_ = 1.5
static constexpr double	trackMaxEtlEta_ = 3.
static constexpr double	trackMinEtlEta_ = 1.6
static constexpr double	zWosMatchMax_ = 1.

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 66 of file Primary4DVertexValidation.cc.

Member Typedef Documentation

LorentzVector

typedef math::XYZTLorentzVector Primary4DVertexValidation::LorentzVector

private

Definition at line 67 of file Primary4DVertexValidation.cc.

Constructor & Destructor Documentation

Primary4DVertexValidation()

Primary4DVertexValidation::Primary4DVertexValidation ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 510 of file Primary4DVertexValidation.cc.

References edm::ParameterSet::getParameter(), ProducerED_cfi::InputTag, momentumToken_, pathLengthToken_, probKToken_, probPiToken_, probPToken_, Rec4DVerToken_, RecBeamSpotToken_, recoToSimAssociationToken_, RecTrackToken_, sigmat0SafePidToken_, simToRecoAssociationToken_, AlCaHLTBitMon_QueryRunRegistry::string, t0PidToken_, t0SafePidToken_, timeToken_, tmtdToken_, tofKToken_, tofPiToken_, tofPToken_, trackAssocToken_, trackingParticleCollectionToken_, trackingVertexCollectionToken_, trackMVAQualToken_, and vecPileupSummaryInfoToken_.

    : folder_(iConfig.getParameter<std::string>("folder")),
      use_only_charged_tracks_(iConfig.getParameter<bool>("useOnlyChargedTracks")),
      debug_(iConfig.getUntrackedParameter<bool>("debug")),
      optionalPlots_(iConfig.getUntrackedParameter<bool>("optionalPlots")),
      use3dNoTime_(iConfig.getParameter<bool>("use3dNoTime")),
      minProbHeavy_(iConfig.getParameter<double>("minProbHeavy")),
      trackweightTh_(iConfig.getParameter<double>("trackweightTh")),
      mvaTh_(iConfig.getParameter<double>("mvaTh")),
      lineDensityPar_(iConfig.getParameter<std::vector<double>>("lineDensityPar")),
      pdtToken_(esConsumes<HepPDT::ParticleDataTable, edm::DefaultRecord>()) {
  vecPileupSummaryInfoToken_ = consumes<std::vector<PileupSummaryInfo>>(edm::InputTag(std::string("addPileupInfo")));
  trackingParticleCollectionToken_ =
      consumes<TrackingParticleCollection>(iConfig.getParameter<edm::InputTag>("SimTag"));
  trackingVertexCollectionToken_ = consumes<TrackingVertexCollection>(iConfig.getParameter<edm::InputTag>("SimTag"));
  simToRecoAssociationToken_ =
      consumes<reco::SimToRecoCollection>(iConfig.getParameter<edm::InputTag>("TPtoRecoTrackAssoc"));
  recoToSimAssociationToken_ =
      consumes<reco::RecoToSimCollection>(iConfig.getParameter<edm::InputTag>("TPtoRecoTrackAssoc"));
  RecTrackToken_ = consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("mtdTracks"));
  RecBeamSpotToken_ = consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("offlineBS"));
  Rec4DVerToken_ = consumes<edm::View<reco::Vertex>>(iConfig.getParameter<edm::InputTag>("offline4DPV"));
  trackAssocToken_ = consumes<edm::ValueMap<int>>(iConfig.getParameter<edm::InputTag>("trackAssocSrc"));
  pathLengthToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("pathLengthSrc"));
  momentumToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("momentumSrc"));
  timeToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("timeSrc"));
  t0PidToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("t0PID"));
  t0SafePidToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("t0SafePID"));
  sigmat0SafePidToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("sigmat0SafePID"));
  trackMVAQualToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("trackMVAQual"));
  tmtdToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("tmtd"));
  tofPiToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("tofPi"));
  tofKToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("tofK"));
  tofPToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("tofP"));
  probPiToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("probPi"));
  probKToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("probK"));
  probPToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("probP"));
}

~Primary4DVertexValidation()

Primary4DVertexValidation::~Primary4DVertexValidation ( )

override

Definition at line 549 of file Primary4DVertexValidation.cc.

{}

Member Function Documentation

analyze()

void Primary4DVertexValidation::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 1885 of file Primary4DVertexValidation.cc.

References funct::abs(), gather_cfg::cout, PVValHelper::d3D, debug_, deltaZcut_, dt, PVValHelper::dz, l1ctLayer1_cff::dZ, JetChargeProducer_cfi::exp, dqm::impl::MonitorElement::Fill(), edm::EventSetup::getHandle(), getMatchedTP(), getRecoPVs(), getSimPVs(), getWosWnt(), caHitNtupletGeneratorKernels::iev, iEvent, edm::isNotFinite(), isParticle(), edm::HandleBase::isValid(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::iv, lineDensityPar_, LogTrace, EgHLTOffHistBins_cfi::mass, matchReco2Sim(), SiStripPI::max, meBarrelNoPIDtype_, meBarrelPIDp_, meBarrelTrueKAsK_, meBarrelTrueKAsP_, meBarrelTrueKAsPi_, meBarrelTrueKNoPID_, meBarrelTruePAsK_, meBarrelTruePAsP_, meBarrelTruePAsPi_, meBarrelTruePiAsK_, meBarrelTruePiAsP_, meBarrelTruePiAsPi_, meBarrelTruePiNoPID_, meBarrelTruePNoPID_, meDeltaTfakefake_, meDeltaTfakereal_, meDeltaTrealreal_, meDeltaZfakefake_, meDeltaZfakereal_, meDeltaZrealreal_, meEndcapNoPIDtype_, meEndcapPIDp_, meEndcapTrueKAsK_, meEndcapTrueKAsP_, meEndcapTrueKAsPi_, meEndcapTrueKNoPID_, meEndcapTruePAsK_, meEndcapTruePAsP_, meEndcapTruePAsPi_, meEndcapTruePiAsK_, meEndcapTruePiAsP_, meEndcapTruePiAsPi_, meEndcapTruePiNoPID_, meEndcapTruePNoPID_, meFakeTrackRecLVRelMult_, meFakeTrackRecLVRelMultvsMult_, meFakeTrackRecLVRelSumWos_, meFakeTrackRecLVRelSumWosvsSumWos_, meFakeTrackRelMult_, meFakeTrackRelMultvsMult_, meFakeTrackRelSumWos_, meFakeTrackRelSumWosvsSumWos_, meJetsFakeRelSumPt2_, meJetsFakeRelSumPt2vsSumPt2_, meJetsPUHt_, meJetsPUMetPt_, meJetsPUMult_, meJetsPURelHt_, meJetsPURelHtvsHt_, meJetsPURelMetPt_, meJetsPURelMetPtvsMetPt_, meJetsPURelMult_, meJetsPURelMultvsMult_, meJetsPURelSumPt2_, meJetsPURelSumPt2vsSumPt2_, meJetsPURelSumPz_, meJetsPURelSumPzvsSumPz_, meJetsPUSumPt2_, meJetsPUSumPz_, meJetsRecLVFakeRelSumPt2_, meJetsRecLVFakeRelSumPt2vsSumPt2_, meJetsRecLVPUHt_, meJetsRecLVPUMetPt_, meJetsRecLVPUMult_, meJetsRecLVPURelHt_, meJetsRecLVPURelHtvsHt_, meJetsRecLVPURelMetPt_, meJetsRecLVPURelMetPtvsMetPt_, meJetsRecLVPURelMult_, meJetsRecLVPURelMultvsMult_, meJetsRecLVPURelSumPt2_, meJetsRecLVPURelSumPt2vsSumPt2_, meJetsRecLVPURelSumPz_, meJetsRecLVPURelSumPzvsSumPz_, meJetsRecLVPUSumPt2_, meJetsRecLVPUSumPz_, meMatchQual_, meMVATrackEffEtaTot_, meMVATrackEffPtTot_, meMVATrackMatchedEffEtaMtd_, meMVATrackMatchedEffEtaTot_, meMVATrackMatchedEffPtMtd_, meMVATrackMatchedEffPtTot_, meMVATrackPullTot_, meMVATrackResTot_, meMVATrackZposResTot_, mePUTrackMult_, mePUTrackRecLVMult_, mePUTrackRecLVRelMult_, mePUTrackRecLVRelMultvsMult_, mePUTrackRecLVRelSumPt2_, mePUTrackRecLVRelSumPt2vsSumPt2_, mePUTrackRecLVRelSumPt_, mePUTrackRecLVRelSumPtvsSumPt_, mePUTrackRecLVRelSumWnt_, mePUTrackRecLVRelSumWntvsSumWnt_, mePUTrackRecLVRelSumWos_, mePUTrackRecLVRelSumWosvsSumWos_, mePUTrackRecLVSumPt2_, mePUTrackRecLVSumPt_, mePUTrackRecLVSumWnt_, mePUTrackRecLVSumWos_, mePUTrackRecLVWnt_, mePUTrackRelMult_, mePUTrackRelMultvsMult_, mePUTrackRelSumPt2_, mePUTrackRelSumPt2vsSumPt2_, mePUTrackRelSumPt_, mePUTrackRelSumPtvsSumPt_, mePUTrackRelSumWnt_, mePUTrackRelSumWntvsSumWnt_, mePUTrackRelSumWos_, mePUTrackRelSumWosvsSumWos_, mePUTrackSumPt2_, mePUTrackSumPt_, mePUTrackSumWnt_, mePUTrackSumWos_, mePUTrackWnt_, mePUvsOtherFakeV_, mePUvsRealV_, mePUvsSplitV_, meRecoPosInRecoOrigCollection_, meRecoPosInSimCollection_, meRecoPVPosSignal_, meRecoPVPosSignalNotHighestPt_, meRecoVtxVsLineDensity_, meRecPVT_, meRecPVZ_, meRecVerNumber_, meSecTrackRecLVRelSumWos_, meSecTrackRecLVRelSumWosvsSumWos_, meSecTrackRecLVSumWos_, meSecTrackRelSumWos_, meSecTrackRelSumWosvsSumWos_, meSecTrackSumWos_, meSimPosInSimOrigCollection_, meSimPVZ_, meTimePull_, meTimeRes_, meTimeSignalPull_, meTimeSignalRes_, meTrack3DposRes_, meTrackPull_, meTrackPullHighP_, meTrackPullHighPTot_, meTrackPullLowP_, meTrackPullLowPTot_, meTrackPullTot_, meTrackRes_, meTrackResHighP_, meTrackResHighPTot_, meTrackResLowP_, meTrackResLowPTot_, meTrackResMass_, meTrackResMassPions_, meTrackResMassProtons_, meTrackResMassTrue_, meTrackResMassTruePions_, meTrackResMassTrueProtons_, meTrackResTot_, meTrackZposRes_, meTrackZposResTot_, minThrMetPt_, minThrSumPt2_, minThrSumPt_, minThrSumPz_, minThrSumWnt_, minThrSumWos_, momentumToken_, mvaH_, mvaL_, mvaRecSel(), mvaTPSel(), nJets, NOT_MATCHED, nt, observablesFromJets(), optionalPlots_, LHEGenericFilter_cfi::ParticleID, genfragment_ptgun_cfg::PartID, pathLengthToken_, pdtToken_, printMatchedRecoTrackInfo(), printSimVtxRecoVtxInfo(), probKToken_, probPiToken_, probPToken_, edm::Handle< T >::product(), Primary4DVertexValidation::simPrimaryVertex::ptsq, r2s_, Rec4DVerToken_, RecBeamSpotToken_, recoToSimAssociationToken_, s2r_, select(), selNdof_, sigmat0SafePidToken_, simToRecoAssociationToken_, simUnit_, submitPVValidationJobs::split(), mathSSE::sqrt(), TtFullHadEvtBuilder_cfi::sumPt, t0PidToken_, t0SafePidToken_, hcalRecHitTable_cff::time, timeFromTrueMass(), timeToken_, tmtdToken_, tofKToken_, tofPiToken_, tofPToken_, tol_, trackAssocToken_, trackingParticleCollectionToken_, trackingVertexCollectionToken_, trackMaxBtlEta_, trackMaxEtlEta_, trackMinEtlEta_, trackMVAQualToken_, trackweightTh_, use3dNoTime_, vecPileupSummaryInfoToken_, trackerHitRTTI::vector, bphysicsOniaDQM_cfi::vertex, and z.

                                                                                           {
  using edm::Handle;
  using edm::View;
  using std::cout;
  using std::endl;
  using std::vector;
  using namespace reco;

  std::vector<float> pileUpInfo_z;

  // get the pileup information
  edm::Handle<std::vector<PileupSummaryInfo>> puinfoH;
  if (iEvent.getByToken(vecPileupSummaryInfoToken_, puinfoH)) {
    for (auto const& pu_info : *puinfoH.product()) {
      if (pu_info.getBunchCrossing() == 0) {
        pileUpInfo_z = pu_info.getPU_zpositions();
        break;
      }
    }
  }

  edm::Handle<TrackingParticleCollection> TPCollectionH;
  iEvent.getByToken(trackingParticleCollectionToken_, TPCollectionH);
  if (!TPCollectionH.isValid())
    edm::LogWarning("Primary4DVertexValidation") << "TPCollectionH is not valid";

  edm::Handle<TrackingVertexCollection> TVCollectionH;
  iEvent.getByToken(trackingVertexCollectionToken_, TVCollectionH);
  if (!TVCollectionH.isValid())
    edm::LogWarning("Primary4DVertexValidation") << "TVCollectionH is not valid";

  edm::Handle<reco::SimToRecoCollection> simToRecoH;
  iEvent.getByToken(simToRecoAssociationToken_, simToRecoH);
  if (simToRecoH.isValid())
    s2r_ = simToRecoH.product();
  else
    edm::LogWarning("Primary4DVertexValidation") << "simToRecoH is not valid";

  edm::Handle<reco::RecoToSimCollection> recoToSimH;
  iEvent.getByToken(recoToSimAssociationToken_, recoToSimH);
  if (recoToSimH.isValid())
    r2s_ = recoToSimH.product();
  else
    edm::LogWarning("Primary4DVertexValidation") << "recoToSimH is not valid";

  edm::Handle<reco::BeamSpot> BeamSpotH;
  iEvent.getByToken(RecBeamSpotToken_, BeamSpotH);
  if (!BeamSpotH.isValid())
    edm::LogWarning("Primary4DVertexValidation") << "BeamSpotH is not valid";

  std::vector<simPrimaryVertex> simpv;  // a list of simulated primary MC vertices
  simpv = getSimPVs(TVCollectionH);
  // this bool check if first vertex in that with highest pT
  bool signal_is_highest_pt =
      std::max_element(simpv.begin(), simpv.end(), [](const simPrimaryVertex& lhs, const simPrimaryVertex& rhs) {
        return lhs.ptsq < rhs.ptsq;
      }) == simpv.begin();

  std::vector<recoPrimaryVertex> recopv;  // a list of reconstructed primary MC vertices
  edm::Handle<edm::View<reco::Vertex>> recVtxs;
  iEvent.getByToken(Rec4DVerToken_, recVtxs);
  if (!recVtxs.isValid())
    edm::LogWarning("Primary4DVertexValidation") << "recVtxs is not valid";
  recopv = getRecoPVs(recVtxs);

  const auto& trackAssoc = iEvent.get(trackAssocToken_);
  const auto& pathLength = iEvent.get(pathLengthToken_);
  const auto& momentum = iEvent.get(momentumToken_);
  const auto& time = iEvent.get(timeToken_);
  const auto& t0Pid = iEvent.get(t0PidToken_);
  const auto& t0Safe = iEvent.get(t0SafePidToken_);
  const auto& sigmat0Safe = iEvent.get(sigmat0SafePidToken_);
  const auto& mtdQualMVA = iEvent.get(trackMVAQualToken_);
  const auto& tMtd = iEvent.get(tmtdToken_);
  const auto& tofPi = iEvent.get(tofPiToken_);
  const auto& tofK = iEvent.get(tofKToken_);
  const auto& tofP = iEvent.get(tofPToken_);
  const auto& probPi = iEvent.get(probPiToken_);
  const auto& probK = iEvent.get(probKToken_);
  const auto& probP = iEvent.get(probPToken_);
  const auto& fPDGTable = iSetup.getHandle(pdtToken_);

  // I have simPV and recoPV collections
  matchReco2Sim(recopv, simpv, sigmat0Safe, mtdQualMVA, BeamSpotH);

  // Loop on tracks
  for (unsigned int iv = 0; iv < recopv.size(); iv++) {
    if (recopv.at(iv).ndof > selNdof_) {
      const reco::Vertex* vertex = recopv.at(iv).recVtx;

      for (unsigned int iev = 0; iev < simpv.size(); iev++) {
        auto vsim = simpv.at(iev).sim_vertex;

        bool selectedVtxMatching = recopv.at(iv).sim == iev && simpv.at(iev).rec == iv;
        bool selectedLV = simpv.at(iev).eventId.bunchCrossing() == 0 && simpv.at(iev).eventId.event() == 0 &&
                          recopv.at(iv).OriginalIndex == 0;
        bool selectedLVMatching = selectedVtxMatching && selectedLV;  // bool for reco vtx leading match
        if (selectedLVMatching && !recopv.at(iv).is_signal()) {
          edm::LogWarning("Primary4DVertexValidation")
              << "Reco vtx leading match inconsistent: BX/ID " << simpv.at(iev).eventId.bunchCrossing() << " "
              << simpv.at(iev).eventId.event();
        }
        if (selectedLVMatching && debug_) {
          printSimVtxRecoVtxInfo(simpv.at(iev), recopv.at(iv));
        }
        double vzsim = simpv.at(iev).z;
        double vtsim = simpv.at(iev).t * simUnit_;

        double wnt = 0, wos = 0;
        double PUsumWnt = 0, PUsumWos = 0, SecsumWos = 0, FakesumWos = 0, PUsumPt = 0, PUsumPt2 = 0;
        double sumWnt = 0, sumWos = 0, sumPt = 0, sumPt2 = 0;
        unsigned int nt = 0, PUnt = 0, Fakent = 0;

        std::vector<double> massVector;
        std::vector<reco::Track> recotracks;
        std::vector<int> categoryVector;
        double sumEtJets = 0, sumPt2Jets = 0, metPt = 0, sumPzJets = 0;
        double sumEtJetsnoPU = 0, sumPt2JetsnoPU = 0, metPtnoPU = 0, sumPzJetsnoPU = 0;
        double sumEtJetsnoFake = 0, sumPt2JetsnoFake = 0, metPtnoFake = 0, sumPzJetsnoFake = 0;
        unsigned int nJets = 0, nJetsnoPU = 0, nJetsnoFake = 0;
        for (auto iTrack = vertex->tracks_begin(); iTrack != vertex->tracks_end(); ++iTrack) {
          if (trackAssoc[*iTrack] == -1) {
            LogTrace("mtdTracks") << "Extended track not associated";
            continue;
          }

          if (vertex->trackWeight(*iTrack) < trackweightTh_)
            continue;
          bool noCrack = std::abs((*iTrack)->eta()) < trackMaxBtlEta_ || std::abs((*iTrack)->eta()) > trackMinEtlEta_;

          bool selectRecoTrk = mvaRecSel(**iTrack, *vertex, t0Safe[*iTrack], sigmat0Safe[*iTrack]);
          if (selectedLVMatching && selectRecoTrk) {
            if (noCrack) {
              meMVATrackEffPtTot_->Fill((*iTrack)->pt());
            }
            meMVATrackEffEtaTot_->Fill(std::abs((*iTrack)->eta()));
          }

          auto tp_info = getMatchedTP(*iTrack, vsim).first;
          int matchCategory = getMatchedTP(*iTrack, vsim).second;

          // PU, fake and secondary tracks
          if (selectedVtxMatching) {
            unsigned int no_PIDtype = 0;
            bool no_PID, is_Pi, is_K, is_P;
            int PartID = 211;  // pion
            isParticle(*iTrack, probPi, probK, probP, no_PIDtype, no_PID, is_Pi, is_K, is_P);
            if (!use3dNoTime_) {
              if (no_PID || is_Pi) {
                PartID = 211;
              } else if (is_K) {
                PartID = 321;
              } else if (is_P) {
                PartID = 2212;
              }
            }
            const HepPDT::ParticleData* PData = fPDGTable->particle(HepPDT::ParticleID(abs(PartID)));
            double mass = PData->mass().value();
            massVector.push_back(mass);
            recotracks.push_back(**iTrack);
            getWosWnt(*vertex, *iTrack, sigmat0Safe, BeamSpotH, wos, wnt);
            // reco track matched to any TP
            if (tp_info != nullptr) {
              if (debug_ && selectedLV) {
                printMatchedRecoTrackInfo(*vertex, *iTrack, *tp_info, matchCategory);
              }
              // matched TP not associated to any daughter track of a given sim vertex but having the same eventID (track from secondary vtx)
              if (matchCategory == 1) {
                categoryVector.push_back(matchCategory);
                SecsumWos += wos;
              }
              // matched TP not associated to any sim vertex of a given simulated event (PU track)
              if (matchCategory == 2) {
                if (optionalPlots_) {
                  mePUTrackWnt_->Fill(wnt);
                  if (selectedLV) {
                    mePUTrackRecLVWnt_->Fill(wnt);
                  }
                }
                PUsumWnt += wnt;
                PUsumWos += wos;
                PUsumPt += (*iTrack)->pt();
                PUsumPt2 += ((*iTrack)->pt() * (*iTrack)->pt());
                PUnt++;
                categoryVector.push_back(2);
              }
            }
            // reco track not matched to any TP (fake tracks)
            else {
              categoryVector.push_back(matchCategory);
              FakesumWos += wos;
              Fakent++;
            }
            nt++;
            sumWnt += wnt;
            sumWos += wos;
            sumPt += (*iTrack)->pt();
            sumPt2 += ((*iTrack)->pt() * (*iTrack)->pt());
          }

          // matched TP equal to any TP of a given sim vertex
          if (tp_info != nullptr && matchCategory == 0) {
            categoryVector.push_back(matchCategory);
            double mass = (*tp_info)->mass();
            double tsim = (*tp_info)->parentVertex()->position().t() * simUnit_;
            double tEst = timeFromTrueMass(mass, pathLength[*iTrack], momentum[*iTrack], time[*iTrack]);

            double xsim = (*tp_info)->parentVertex()->position().x();
            double ysim = (*tp_info)->parentVertex()->position().y();
            double zsim = (*tp_info)->parentVertex()->position().z();
            double xPCA = (*iTrack)->vx();
            double yPCA = (*iTrack)->vy();
            double zPCA = (*iTrack)->vz();

            double dZ = zPCA - zsim;
            double d3D = std::sqrt((xPCA - xsim) * (xPCA - xsim) + (yPCA - ysim) * (yPCA - ysim) + dZ * dZ);
            // orient d3D according to the projection of RECO - SIM onto simulated momentum
            if ((xPCA - xsim) * ((*tp_info)->px()) + (yPCA - ysim) * ((*tp_info)->py()) + dZ * ((*tp_info)->pz()) <
                0.) {
              d3D = -d3D;
            }

            bool selectTP = mvaTPSel(**tp_info);

            if (selectedLVMatching && selectRecoTrk && selectTP) {
              meMVATrackZposResTot_->Fill((*iTrack)->vz() - vzsim);
              if (noCrack) {
                meMVATrackMatchedEffPtTot_->Fill((*iTrack)->pt());
              }
              meMVATrackMatchedEffEtaTot_->Fill(std::abs((*iTrack)->eta()));
            }

            if (sigmat0Safe[*iTrack] == -1)
              continue;

            if (selectedLVMatching && selectRecoTrk && selectTP) {
              meMVATrackResTot_->Fill(t0Safe[*iTrack] - vtsim);
              meMVATrackPullTot_->Fill((t0Safe[*iTrack] - vtsim) / sigmat0Safe[*iTrack]);
              if (noCrack) {
                meMVATrackMatchedEffPtMtd_->Fill((*iTrack)->pt());
              }
              meMVATrackMatchedEffEtaMtd_->Fill(std::abs((*iTrack)->eta()));

              unsigned int noPIDtype = 0;
              bool noPID = false, isPi = false, isK = false, isP = false;
              isParticle(*iTrack, probPi, probK, probP, noPIDtype, noPID, isPi, isK, isP);

              if ((isPi && std::abs(tMtd[*iTrack] - tofPi[*iTrack] - t0Pid[*iTrack]) > tol_) ||
                  (isK && std::abs(tMtd[*iTrack] - tofK[*iTrack] - t0Pid[*iTrack]) > tol_) ||
                  (isP && std::abs(tMtd[*iTrack] - tofP[*iTrack] - t0Pid[*iTrack]) > tol_)) {
                edm::LogWarning("Primary4DVertexValidation")
                    << "No match between mass hyp. and time: " << std::abs((*tp_info)->pdgId()) << " mass hyp pi/k/p "
                    << isPi << " " << isK << " " << isP << " t0/t0safe " << t0Pid[*iTrack] << " " << t0Safe[*iTrack]
                    << " tMtd - tof pi/K/p " << tMtd[*iTrack] - tofPi[*iTrack] << " " << tMtd[*iTrack] - tofK[*iTrack]
                    << " " << tMtd[*iTrack] - tofP[*iTrack] << " Prob pi/K/p " << probPi[*iTrack] << " "
                    << probK[*iTrack] << " " << probP[*iTrack];
              }

              if (std::abs((*iTrack)->eta()) < trackMaxBtlEta_) {
                meBarrelPIDp_->Fill((*iTrack)->p());
                meBarrelNoPIDtype_->Fill(noPIDtype + 0.5);
                if (std::abs((*tp_info)->pdgId()) == 211) {
                  if (noPID) {
                    meBarrelTruePiNoPID_->Fill((*iTrack)->p());
                  } else if (isPi) {
                    meBarrelTruePiAsPi_->Fill((*iTrack)->p());
                  } else if (isK) {
                    meBarrelTruePiAsK_->Fill((*iTrack)->p());
                  } else if (isP) {
                    meBarrelTruePiAsP_->Fill((*iTrack)->p());
                  } else {
                    edm::LogWarning("Primary4DVertexValidation")
                        << "No PID class: " << std::abs((*tp_info)->pdgId()) << " t0/t0safe " << t0Pid[*iTrack] << " "
                        << t0Safe[*iTrack] << " Prob pi/K/p " << probPi[*iTrack] << " " << probK[*iTrack] << " "
                        << probP[*iTrack];
                  }
                } else if (std::abs((*tp_info)->pdgId()) == 321) {
                  if (noPID) {
                    meBarrelTrueKNoPID_->Fill((*iTrack)->p());
                  } else if (isPi) {
                    meBarrelTrueKAsPi_->Fill((*iTrack)->p());
                  } else if (isK) {
                    meBarrelTrueKAsK_->Fill((*iTrack)->p());
                  } else if (isP) {
                    meBarrelTrueKAsP_->Fill((*iTrack)->p());
                  } else {
                    edm::LogWarning("Primary4DVertexValidation")
                        << "No PID class: " << std::abs((*tp_info)->pdgId()) << " t0/t0safe " << t0Pid[*iTrack] << " "
                        << t0Safe[*iTrack] << " Prob pi/K/p " << probPi[*iTrack] << " " << probK[*iTrack] << " "
                        << probP[*iTrack];
                  }
                } else if (std::abs((*tp_info)->pdgId()) == 2212) {
                  if (noPID) {
                    meBarrelTruePNoPID_->Fill((*iTrack)->p());
                  } else if (isPi) {
                    meBarrelTruePAsPi_->Fill((*iTrack)->p());
                  } else if (isK) {
                    meBarrelTruePAsK_->Fill((*iTrack)->p());
                  } else if (isP) {
                    meBarrelTruePAsP_->Fill((*iTrack)->p());
                  } else {
                    edm::LogWarning("Primary4DVertexValidation")
                        << "No PID class: " << std::abs((*tp_info)->pdgId()) << " t0/t0safe " << t0Pid[*iTrack] << " "
                        << t0Safe[*iTrack] << " Prob pi/K/p " << probPi[*iTrack] << " " << probK[*iTrack] << " "
                        << probP[*iTrack];
                  }
                }
              } else if (std::abs((*iTrack)->eta()) > trackMinEtlEta_ && std::abs((*iTrack)->eta()) < trackMaxEtlEta_) {
                meEndcapPIDp_->Fill((*iTrack)->p());
                meEndcapNoPIDtype_->Fill(noPIDtype + 0.5);
                if (std::abs((*tp_info)->pdgId()) == 211) {
                  if (noPID) {
                    meEndcapTruePiNoPID_->Fill((*iTrack)->p());
                  } else if (isPi) {
                    meEndcapTruePiAsPi_->Fill((*iTrack)->p());
                  } else if (isK) {
                    meEndcapTruePiAsK_->Fill((*iTrack)->p());
                  } else if (isP) {
                    meEndcapTruePiAsP_->Fill((*iTrack)->p());
                  } else {
                    edm::LogWarning("Primary4DVertexValidation")
                        << "No PID class: " << std::abs((*tp_info)->pdgId()) << " t0/t0safe " << t0Pid[*iTrack] << " "
                        << t0Safe[*iTrack] << " Prob pi/K/p " << probPi[*iTrack] << " " << probK[*iTrack] << " "
                        << probP[*iTrack];
                  }
                } else if (std::abs((*tp_info)->pdgId()) == 321) {
                  if (noPID) {
                    meEndcapTrueKNoPID_->Fill((*iTrack)->p());
                  } else if (isPi) {
                    meEndcapTrueKAsPi_->Fill((*iTrack)->p());
                  } else if (isK) {
                    meEndcapTrueKAsK_->Fill((*iTrack)->p());
                  } else if (isP) {
                    meEndcapTrueKAsP_->Fill((*iTrack)->p());
                  } else {
                    edm::LogWarning("Primary4DVertexValidation")
                        << "No PID class: " << std::abs((*tp_info)->pdgId()) << " t0/t0safe " << t0Pid[*iTrack] << " "
                        << t0Safe[*iTrack] << " Prob pi/K/p " << probPi[*iTrack] << " " << probK[*iTrack] << " "
                        << probP[*iTrack];
                  }
                } else if (std::abs((*tp_info)->pdgId()) == 2212) {
                  if (noPID) {
                    meEndcapTruePNoPID_->Fill((*iTrack)->p());
                  } else if (isPi) {
                    meEndcapTruePAsPi_->Fill((*iTrack)->p());
                  } else if (isK) {
                    meEndcapTruePAsK_->Fill((*iTrack)->p());
                  } else if (isP) {
                    meEndcapTruePAsP_->Fill((*iTrack)->p());
                  } else {
                    edm::LogWarning("Primary4DVertexValidation")
                        << "No PID class: " << std::abs((*tp_info)->pdgId()) << " t0/t0safe " << t0Pid[*iTrack] << " "
                        << t0Safe[*iTrack] << " Prob pi/K/p " << probPi[*iTrack] << " " << probK[*iTrack] << " "
                        << probP[*iTrack];
                  }
                }
              }
            }
            meTrackResTot_->Fill(t0Safe[*iTrack] - tsim);
            meTrackPullTot_->Fill((t0Safe[*iTrack] - tsim) / sigmat0Safe[*iTrack]);
            meTrackZposResTot_->Fill(dZ);
            if ((*iTrack)->p() <= 2) {
              meTrackResLowPTot_->Fill(t0Safe[*iTrack] - tsim);
              meTrackPullLowPTot_->Fill((t0Safe[*iTrack] - tsim) / sigmat0Safe[*iTrack]);
            } else {
              meTrackResHighPTot_->Fill(t0Safe[*iTrack] - tsim);
              meTrackPullHighPTot_->Fill((t0Safe[*iTrack] - tsim) / sigmat0Safe[*iTrack]);
            }

            if (mtdQualMVA[(*iTrack)] < mvaL_) {
              meTrackZposRes_[0]->Fill(dZ);
              meTrack3DposRes_[0]->Fill(d3D);
              meTrackRes_[0]->Fill(t0Safe[*iTrack] - tsim);
              meTrackPull_[0]->Fill((t0Safe[*iTrack] - tsim) / sigmat0Safe[*iTrack]);

              if (optionalPlots_) {
                meTrackResMass_[0]->Fill(t0Safe[*iTrack] - tEst);
                meTrackResMassTrue_[0]->Fill(tEst - tsim);
              }

              if ((*iTrack)->p() <= 2) {
                meTrackResLowP_[0]->Fill(t0Safe[*iTrack] - tsim);
                meTrackPullLowP_[0]->Fill((t0Safe[*iTrack] - tsim) / sigmat0Safe[*iTrack]);
              } else if ((*iTrack)->p() > 2) {
                meTrackResHighP_[0]->Fill(t0Safe[*iTrack] - tsim);
                meTrackPullHighP_[0]->Fill((t0Safe[*iTrack] - tsim) / sigmat0Safe[*iTrack]);
              }

              if (optionalPlots_) {
                if (std::abs((*tp_info)->pdgId()) == 2212) {
                  meTrackResMassProtons_[0]->Fill(t0Safe[*iTrack] - tEst);
                  meTrackResMassTrueProtons_[0]->Fill(tEst - tsim);
                } else if (std::abs((*tp_info)->pdgId()) == 211) {
                  meTrackResMassPions_[0]->Fill(t0Safe[*iTrack] - tEst);
                  meTrackResMassTruePions_[0]->Fill(tEst - tsim);
                }
              }

            } else if (mtdQualMVA[(*iTrack)] > mvaL_ && mtdQualMVA[(*iTrack)] < mvaH_) {
              meTrackZposRes_[1]->Fill(dZ);
              meTrack3DposRes_[1]->Fill(d3D);
              meTrackRes_[1]->Fill(t0Safe[*iTrack] - tsim);
              meTrackPull_[1]->Fill((t0Safe[*iTrack] - tsim) / sigmat0Safe[*iTrack]);

              if (optionalPlots_) {
                meTrackResMass_[1]->Fill(t0Safe[*iTrack] - tEst);
                meTrackResMassTrue_[1]->Fill(tEst - tsim);
              }

              if ((*iTrack)->p() <= 2) {
                meTrackResLowP_[1]->Fill(t0Safe[*iTrack] - tsim);
                meTrackPullLowP_[1]->Fill((t0Safe[*iTrack] - tsim) / sigmat0Safe[*iTrack]);
              } else if ((*iTrack)->p() > 2) {
                meTrackResHighP_[1]->Fill(t0Safe[*iTrack] - tsim);
                meTrackPullHighP_[1]->Fill((t0Safe[*iTrack] - tsim) / sigmat0Safe[*iTrack]);
              }

              if (optionalPlots_) {
                if (std::abs((*tp_info)->pdgId()) == 2212) {
                  meTrackResMassProtons_[1]->Fill(t0Safe[*iTrack] - tEst);
                  meTrackResMassTrueProtons_[1]->Fill(tEst - tsim);
                } else if (std::abs((*tp_info)->pdgId()) == 211) {
                  meTrackResMassPions_[1]->Fill(t0Safe[*iTrack] - tEst);
                  meTrackResMassTruePions_[1]->Fill(tEst - tsim);
                }
              }

            } else if (mtdQualMVA[(*iTrack)] > mvaH_) {
              meTrackZposRes_[2]->Fill(dZ);
              meTrack3DposRes_[2]->Fill(d3D);
              meTrackRes_[2]->Fill(t0Safe[*iTrack] - tsim);
              meTrackPull_[2]->Fill((t0Safe[*iTrack] - tsim) / sigmat0Safe[*iTrack]);

              if (optionalPlots_) {
                meTrackResMass_[2]->Fill(t0Safe[*iTrack] - tEst);
                meTrackResMassTrue_[2]->Fill(tEst - tsim);
              }

              if ((*iTrack)->p() <= 2) {
                meTrackResLowP_[2]->Fill(t0Safe[*iTrack] - tsim);
                meTrackPullLowP_[2]->Fill((t0Safe[*iTrack] - tsim) / sigmat0Safe[*iTrack]);
              } else if ((*iTrack)->p() > 2) {
                meTrackResHighP_[2]->Fill(t0Safe[*iTrack] - tsim);
                meTrackPullHighP_[2]->Fill((t0Safe[*iTrack] - tsim) / sigmat0Safe[*iTrack]);
              }

              if (optionalPlots_) {
                if (std::abs((*tp_info)->pdgId()) == 2212) {
                  meTrackResMassProtons_[2]->Fill(t0Safe[*iTrack] - tEst);
                  meTrackResMassTrueProtons_[2]->Fill(tEst - tsim);
                } else if (std::abs((*tp_info)->pdgId()) == 211) {
                  meTrackResMassPions_[2]->Fill(t0Safe[*iTrack] - tEst);
                  meTrackResMassTruePions_[2]->Fill(tEst - tsim);
                }
              }
            }
          }  // if tp_info != nullptr && MatchCategory == 0
        }    // loop on reco tracks
        if (selectedVtxMatching) {
          mePUTrackRelMult_->Fill(static_cast<double>(PUnt) / nt);
          meFakeTrackRelMult_->Fill(static_cast<double>(Fakent) / nt);
          mePUTrackRelSumWnt_->Fill(PUsumWnt / sumWnt);
          mePUTrackRelSumWos_->Fill(PUsumWos / sumWos);
          meSecTrackRelSumWos_->Fill(SecsumWos / sumWos);
          meFakeTrackRelSumWos_->Fill(FakesumWos / sumWos);
          mePUTrackRelSumPt_->Fill(PUsumPt / sumPt);
          mePUTrackRelSumPt2_->Fill(PUsumPt2 / sumPt2);

          observablesFromJets(
              recotracks, massVector, categoryVector, "use_allTracks", nJets, sumEtJets, sumPt2Jets, metPt, sumPzJets);
          observablesFromJets(recotracks,
                              massVector,
                              categoryVector,
                              "skip_Fake",
                              nJetsnoFake,
                              sumEtJetsnoFake,
                              sumPt2JetsnoFake,
                              metPtnoFake,
                              sumPzJetsnoFake);
          observablesFromJets(recotracks,
                              massVector,
                              categoryVector,
                              "skip_PU",
                              nJetsnoPU,
                              sumEtJetsnoPU,
                              sumPt2JetsnoPU,
                              metPtnoPU,
                              sumPzJetsnoPU);

          meJetsPURelMult_->Fill(static_cast<double>(nJets - nJetsnoPU) / nJets);
          meJetsPURelHt_->Fill((sumEtJets - sumEtJetsnoPU) / sumEtJets);
          meJetsPURelSumPt2_->Fill((sumPt2Jets - sumPt2JetsnoPU) / sumPt2Jets);
          meJetsFakeRelSumPt2_->Fill((sumPt2Jets - sumPt2JetsnoFake) / sumPt2Jets);
          meJetsPURelMetPt_->Fill((metPt - metPtnoPU) / metPt);
          meJetsPURelSumPz_->Fill((sumPzJets - sumPzJetsnoPU) / sumPzJets);

          if (optionalPlots_) {
            mePUTrackMult_->Fill(PUnt);
            mePUTrackSumWnt_->Fill(log10(std::max(minThrSumWnt_, PUsumWnt)));
            mePUTrackSumWos_->Fill(log10(std::max(minThrSumWos_, PUsumWos)));
            meSecTrackSumWos_->Fill(log10(std::max(minThrSumWos_, SecsumWos)));
            mePUTrackSumPt_->Fill(log10(std::max(minThrSumPt_, PUsumPt)));
            mePUTrackSumPt2_->Fill(log10(std::max(minThrSumPt2_, PUsumPt2)));

            mePUTrackRelMultvsMult_->Fill(nt, static_cast<double>(PUnt) / nt);
            meFakeTrackRelMultvsMult_->Fill(nt, static_cast<double>(Fakent) / nt);
            mePUTrackRelSumWntvsSumWnt_->Fill(log10(std::max(minThrSumWnt_, sumWnt)), PUsumWnt / sumWnt);
            mePUTrackRelSumWosvsSumWos_->Fill(log10(std::max(minThrSumWos_, sumWos)), PUsumWos / sumWos);
            meSecTrackRelSumWosvsSumWos_->Fill(log10(std::max(minThrSumWos_, sumWos)), SecsumWos / sumWos);
            meFakeTrackRelSumWosvsSumWos_->Fill(log10(std::max(minThrSumWos_, sumWos)), FakesumWos / sumWos);
            mePUTrackRelSumPtvsSumPt_->Fill(log10(std::max(minThrSumPt_, sumPt)), PUsumPt / sumPt);
            mePUTrackRelSumPt2vsSumPt2_->Fill(log10(std::max(minThrSumPt2_, sumPt2)), PUsumPt2 / sumPt2);

            meJetsPUMult_->Fill(nJets - nJetsnoPU);
            meJetsPUHt_->Fill(log10(std::max(minThrSumPt_, sumEtJets - sumEtJetsnoPU)));
            meJetsPUSumPt2_->Fill(log10(std::max(minThrSumPt2_, sumPt2Jets - sumPt2JetsnoPU)));
            meJetsPUMetPt_->Fill(log10(std::max(minThrMetPt_, metPt - metPtnoPU)));
            meJetsPUSumPz_->Fill(log10(std::max(minThrSumPz_, std::abs(sumPzJets - sumPzJetsnoPU))));

            meJetsPURelMultvsMult_->Fill(nJets, static_cast<double>(nJets - nJetsnoPU) / nJets);
            meJetsPURelHtvsHt_->Fill(log10(std::max(minThrSumPt_, sumEtJets)), (sumEtJets - sumEtJetsnoPU) / sumEtJets);
            meJetsPURelSumPt2vsSumPt2_->Fill(log10(std::max(minThrSumPt2_, sumPt2Jets)),
                                             (sumPt2Jets - sumPt2JetsnoPU) / sumPt2Jets);
            meJetsFakeRelSumPt2vsSumPt2_->Fill(log10(std::max(minThrSumPt2_, sumPt2Jets)),
                                               (sumPt2Jets - sumPt2JetsnoFake) / sumPt2Jets);
            meJetsPURelMetPtvsMetPt_->Fill(log10(std::max(minThrMetPt_, metPt)), (metPt - metPtnoPU) / metPt);
            meJetsPURelSumPzvsSumPz_->Fill(log10(std::max(minThrSumPz_, std::abs(sumPzJets))),
                                           (sumPzJets - sumPzJetsnoPU) / sumPzJets);
          }
          if (selectedLV) {
            mePUTrackRecLVRelMult_->Fill(static_cast<double>(PUnt) / nt);
            meFakeTrackRecLVRelMult_->Fill(static_cast<double>(Fakent) / nt);
            mePUTrackRecLVRelSumWnt_->Fill(PUsumWnt / sumWnt);
            mePUTrackRecLVRelSumWos_->Fill(PUsumWos / sumWos);
            meSecTrackRecLVRelSumWos_->Fill(SecsumWos / sumWos);
            meFakeTrackRecLVRelSumWos_->Fill(FakesumWos / sumWos);
            mePUTrackRecLVRelSumPt_->Fill(PUsumPt / sumPt);
            mePUTrackRecLVRelSumPt2_->Fill(PUsumPt2 / sumPt2);

            meJetsRecLVPURelMult_->Fill(static_cast<double>(nJets - nJetsnoPU) / nJets);
            meJetsRecLVPURelHt_->Fill((sumEtJets - sumEtJetsnoPU) / sumEtJets);
            meJetsRecLVPURelSumPt2_->Fill((sumPt2Jets - sumPt2JetsnoPU) / sumPt2Jets);
            meJetsRecLVFakeRelSumPt2_->Fill((sumPt2Jets - sumPt2JetsnoFake) / sumPt2Jets);
            meJetsRecLVPURelMetPt_->Fill((metPt - metPtnoPU) / metPt);
            meJetsRecLVPURelSumPz_->Fill((sumPzJets - sumPzJetsnoPU) / sumPzJets);

            if (debug_) {
              edm::LogPrint("Primary4DVertexValidation")
                  << "#PUTrks = " << PUnt << " #Trks = " << nt << " PURelMult = " << std::setprecision(3)
                  << static_cast<double>(PUnt) / nt;
              edm::LogPrint("Primary4DVertexValidation")
                  << "PUsumWnt = " << std::setprecision(3) << PUsumWnt << " sumWnt = " << std::setprecision(3) << sumWnt
                  << " PURelsumWnt = " << std::setprecision(3) << PUsumWnt / sumWnt;
              edm::LogPrint("Primary4DVertexValidation")
                  << "PUsumWos = " << std::setprecision(3) << PUsumWos << " sumWos = " << std::setprecision(3) << sumWos
                  << " PURelsumWos = " << std::setprecision(3) << PUsumWos / sumWos;
              edm::LogPrint("Primary4DVertexValidation")
                  << "PuSumPt = " << std::setprecision(3) << PUsumPt << " SumPt = " << std::setprecision(4) << sumPt
                  << " PURelSumPt = " << std::setprecision(3) << PUsumPt / sumPt;
              edm::LogPrint("Primary4DVertexValidation")
                  << "PuSumPt2 = " << std::setprecision(3) << PUsumPt2 << " SumPt2 = " << std::setprecision(4) << sumPt2
                  << " PURelSumPt2 = " << std::setprecision(3) << PUsumPt2 / sumPt2;
            }
            if (optionalPlots_) {
              mePUTrackRecLVMult_->Fill(PUnt);
              mePUTrackRecLVSumWnt_->Fill(log10(std::max(minThrSumWnt_, PUsumWnt)));
              mePUTrackRecLVSumWos_->Fill(log10(std::max(minThrSumWos_, PUsumWos)));
              meSecTrackRecLVSumWos_->Fill(log10(std::max(minThrSumWos_, PUsumWos)));
              mePUTrackRecLVSumPt_->Fill(log10(std::max(minThrSumPt_, PUsumPt)));
              mePUTrackRecLVSumPt2_->Fill(log10(std::max(minThrSumPt2_, PUsumPt2)));

              mePUTrackRecLVRelMultvsMult_->Fill(nt, static_cast<double>(PUnt) / nt);
              meFakeTrackRecLVRelMultvsMult_->Fill(nt, static_cast<double>(Fakent) / nt);
              mePUTrackRecLVRelSumWntvsSumWnt_->Fill(log10(std::max(minThrSumWnt_, sumWnt)), PUsumWnt / sumWnt);
              mePUTrackRecLVRelSumWosvsSumWos_->Fill(log10(std::max(minThrSumWos_, sumWos)), PUsumWos / sumWos);
              meSecTrackRecLVRelSumWosvsSumWos_->Fill(log10(std::max(minThrSumWos_, sumWos)), SecsumWos / sumWos);
              meFakeTrackRecLVRelSumWosvsSumWos_->Fill(log10(std::max(minThrSumWos_, sumWos)), FakesumWos / sumWos);
              mePUTrackRecLVRelSumPtvsSumPt_->Fill(log10(std::max(minThrSumPt_, sumPt)), PUsumPt / sumPt);
              mePUTrackRecLVRelSumPt2vsSumPt2_->Fill(log10(std::max(minThrSumPt2_, sumPt2)), PUsumPt2 / sumPt2);

              meJetsRecLVPUMult_->Fill(nJets - nJetsnoPU);
              meJetsRecLVPUHt_->Fill(log10(std::max(minThrSumPt_, sumEtJets - sumEtJetsnoPU)));
              meJetsRecLVPUSumPt2_->Fill(log10(std::max(minThrSumPt2_, sumPt2Jets - sumPt2JetsnoPU)));
              meJetsRecLVPUMetPt_->Fill(log10(std::max(minThrMetPt_, metPt - metPtnoPU)));
              meJetsRecLVPUSumPz_->Fill(log10(std::max(minThrSumPz_, std::abs(sumPzJets - sumPzJetsnoPU))));

              meJetsRecLVPURelMultvsMult_->Fill(nJets, static_cast<double>(nJets - nJetsnoPU) / nJets);
              meJetsRecLVPURelHtvsHt_->Fill(log10(std::max(minThrSumPt_, sumEtJets)),
                                            (sumEtJets - sumEtJetsnoPU) / sumEtJets);
              meJetsRecLVPURelSumPt2vsSumPt2_->Fill(log10(std::max(minThrSumPt2_, sumPt2Jets)),
                                                    (sumPt2Jets - sumPt2JetsnoPU) / sumPt2Jets);
              meJetsRecLVFakeRelSumPt2vsSumPt2_->Fill(log10(std::max(minThrSumPt2_, sumPt2Jets)),
                                                      (sumPt2Jets - sumPt2JetsnoFake) / sumPt2Jets);
              meJetsRecLVPURelMetPtvsMetPt_->Fill(log10(std::max(minThrMetPt_, metPt)), (metPt - metPtnoPU) / metPt);
              meJetsRecLVPURelSumPzvsSumPz_->Fill(log10(std::max(minThrSumPz_, std::abs(sumPzJets))),
                                                  (sumPzJets - sumPzJetsnoPU) / sumPzJets);
            }
          }
        }
      }  // loop on simpv
    }    // ndof
  }      // loop on recopv

  int real = 0;
  int fake = 0;
  int other_fake = 0;
  int split = 0;

  auto puLineDensity = [&](double z) {
    // gaussian parameterization of line density vs z, z in cm, parameters in mm
    double argl = (z * 10. - lineDensityPar_[1]) / lineDensityPar_[2];
    return lineDensityPar_[0] * exp(-0.5 * argl * argl);
  };

  meRecVerNumber_->Fill(recopv.size());
  for (unsigned int ir = 0; ir < recopv.size(); ir++) {
    if (recopv.at(ir).ndof > selNdof_) {
      meRecoVtxVsLineDensity_->Fill(puLineDensity(recopv.at(ir).z));
      meRecPVZ_->Fill(recopv.at(ir).z, 1. / puLineDensity(recopv.at(ir).z));
      if (recopv.at(ir).recVtx->tError() > 0.) {
        meRecPVT_->Fill(recopv.at(ir).recVtx->t());
      }
      if (debug_) {
        edm::LogPrint("Primary4DVertexValidation") << "************* IR: " << ir;
        edm::LogPrint("Primary4DVertexValidation")
            << "z: " << recopv.at(ir).z << " corresponding to line density: " << puLineDensity(recopv.at(ir).z);
        edm::LogPrint("Primary4DVertexValidation") << "is_real: " << recopv.at(ir).is_real();
        edm::LogPrint("Primary4DVertexValidation") << "is_fake: " << recopv.at(ir).is_fake();
        edm::LogPrint("Primary4DVertexValidation") << "is_signal: " << recopv.at(ir).is_signal();
        edm::LogPrint("Primary4DVertexValidation") << "split_from: " << recopv.at(ir).split_from();
        edm::LogPrint("Primary4DVertexValidation") << "other fake: " << recopv.at(ir).other_fake();
      }
      if (recopv.at(ir).is_real())
        real++;
      if (recopv.at(ir).is_fake())
        fake++;
      if (recopv.at(ir).other_fake())
        other_fake++;
      if (recopv.at(ir).split_from() != -1) {
        split++;
      }
    }  // ndof
  }

  if (debug_) {
    edm::LogPrint("Primary4DVertexValidation") << "is_real: " << real;
    edm::LogPrint("Primary4DVertexValidation") << "is_fake: " << fake;
    edm::LogPrint("Primary4DVertexValidation") << "split_from: " << split;
    edm::LogPrint("Primary4DVertexValidation") << "other fake: " << other_fake;
  }
  mePUvsRealV_->Fill(simpv.size(), real);
  mePUvsOtherFakeV_->Fill(simpv.size(), other_fake);
  mePUvsSplitV_->Fill(simpv.size(), split);

  // fill vertices histograms here in a new loop
  for (unsigned int is = 0; is < simpv.size(); is++) {
    // protect against particle guns with very displaced vertices
    if (edm::isNotFinite(1. / puLineDensity(simpv.at(is).z))) {
      continue;
    }
    meSimPVZ_->Fill(simpv.at(is).z, 1. / puLineDensity(simpv.at(is).z));
    if (is == 0 && optionalPlots_) {
      meSimPosInSimOrigCollection_->Fill(simpv.at(is).OriginalIndex);
    }

    if (simpv.at(is).rec == NOT_MATCHED) {
      if (debug_) {
        edm::LogPrint("Primary4DVertexValidation") << "sim vertex: " << is << " is not matched with any reco";
      }
      continue;
    }

    for (unsigned int ir = 0; ir < recopv.size(); ir++) {
      if (recopv.at(ir).ndof > selNdof_) {
        if (recopv.at(ir).sim == is && simpv.at(is).rec == ir) {
          meTimeRes_->Fill(recopv.at(ir).recVtx->t() - simpv.at(is).t * simUnit_);
          meTimePull_->Fill((recopv.at(ir).recVtx->t() - simpv.at(is).t * simUnit_) / recopv.at(ir).recVtx->tError());
          meMatchQual_->Fill(recopv.at(ir).matchQuality - 0.5);
          if (ir == 0) {  // signal vertex plots
            meTimeSignalRes_->Fill(recopv.at(ir).recVtx->t() - simpv.at(is).t * simUnit_);
            meTimeSignalPull_->Fill((recopv.at(ir).recVtx->t() - simpv.at(is).t * simUnit_) /
                                    recopv.at(ir).recVtx->tError());
            if (optionalPlots_) {
              meRecoPosInSimCollection_->Fill(recopv.at(ir).sim);
              meRecoPosInRecoOrigCollection_->Fill(recopv.at(ir).OriginalIndex);
            }
          }
          if (simpv.at(is).eventId.bunchCrossing() == 0 && simpv.at(is).eventId.event() == 0) {
            if (!recopv.at(ir).is_signal()) {
              edm::LogWarning("Primary4DVertexValidation")
                  << "Reco vtx leading match inconsistent: BX/ID " << simpv.at(is).eventId.bunchCrossing() << " "
                  << simpv.at(is).eventId.event();
            }
            meRecoPVPosSignal_->Fill(
                recopv.at(ir).OriginalIndex);  // position in reco vtx correction associated to sim signal
            if (!signal_is_highest_pt) {
              meRecoPVPosSignalNotHighestPt_->Fill(
                  recopv.at(ir).OriginalIndex);  // position in reco vtx correction associated to sim signal
            }
          }
          if (debug_) {
            edm::LogPrint("Primary4DVertexValidation") << "*** Matching RECO: " << ir << "with SIM: " << is << " ***";
            edm::LogPrint("Primary4DVertexValidation") << "Match Quality is " << recopv.at(ir).matchQuality;
            edm::LogPrint("Primary4DVertexValidation") << "****";
          }
        }
      }  // ndof
    }
  }

  // dz histos
  for (unsigned int iv = 0; iv < recVtxs->size() - 1; iv++) {
    if (recVtxs->at(iv).ndof() > selNdof_) {
      double mindistance_realreal = 1e10;

      for (unsigned int jv = iv; jv < recVtxs->size(); jv++) {
        if ((!(jv == iv)) && select(recVtxs->at(jv))) {
          double dz = recVtxs->at(iv).z() - recVtxs->at(jv).z();
          double dtsigma = std::sqrt(recVtxs->at(iv).covariance(3, 3) + recVtxs->at(jv).covariance(3, 3));
          double dt = (std::abs(dz) <= deltaZcut_ && dtsigma > 0.)
                          ? (recVtxs->at(iv).t() - recVtxs->at(jv).t()) / dtsigma
                          : -9999.;
          if (recopv.at(iv).is_real() && recopv.at(jv).is_real()) {
            meDeltaZrealreal_->Fill(std::abs(dz));
            if (dt != -9999.) {
              meDeltaTrealreal_->Fill(std::abs(dt));
            }
            if (std::abs(dz) < std::abs(mindistance_realreal)) {
              mindistance_realreal = dz;
            }
          } else if (recopv.at(iv).is_fake() && recopv.at(jv).is_fake()) {
            meDeltaZfakefake_->Fill(std::abs(dz));
            if (dt != -9999.) {
              meDeltaTfakefake_->Fill(std::abs(dt));
            }
          }
        }
      }

      double mindistance_fakereal = 1e10;
      double mindistance_realfake = 1e10;
      for (unsigned int jv = 0; jv < recVtxs->size(); jv++) {
        if ((!(jv == iv)) && select(recVtxs->at(jv))) {
          double dz = recVtxs->at(iv).z() - recVtxs->at(jv).z();
          double dtsigma = std::sqrt(recVtxs->at(iv).covariance(3, 3) + recVtxs->at(jv).covariance(3, 3));
          double dt = (std::abs(dz) <= deltaZcut_ && dtsigma > 0.)
                          ? (recVtxs->at(iv).t() - recVtxs->at(jv).t()) / dtsigma
                          : -9999.;
          if (recopv.at(iv).is_fake() && recopv.at(jv).is_real()) {
            meDeltaZfakereal_->Fill(std::abs(dz));
            if (dt != -9999.) {
              meDeltaTfakereal_->Fill(std::abs(dt));
            }
            if (std::abs(dz) < std::abs(mindistance_fakereal)) {
              mindistance_fakereal = dz;
            }
          }

          if (recopv.at(iv).is_real() && recopv.at(jv).is_fake() && (std::abs(dz) < std::abs(mindistance_realfake))) {
            mindistance_realfake = dz;
          }
        }
      }
    }  // ndof
  }

}  // end of analyze

bookHistograms()

void Primary4DVertexValidation::bookHistograms ( DQMStore::IBooker &  i, edm::Run const &  iRun, edm::EventSetup const &  iSetup  )

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 554 of file Primary4DVertexValidation.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::bookProfile(), folder_, meBarrelNoPIDtype_, meBarrelPIDp_, meBarrelTrueKAsK_, meBarrelTrueKAsP_, meBarrelTrueKAsPi_, meBarrelTrueKNoPID_, meBarrelTruePAsK_, meBarrelTruePAsP_, meBarrelTruePAsPi_, meBarrelTruePiAsK_, meBarrelTruePiAsP_, meBarrelTruePiAsPi_, meBarrelTruePiNoPID_, meBarrelTruePNoPID_, meDeltaTfakefake_, meDeltaTfakereal_, meDeltaTrealreal_, meDeltaZfakefake_, meDeltaZfakereal_, meDeltaZrealreal_, meEndcapNoPIDtype_, meEndcapPIDp_, meEndcapTrueKAsK_, meEndcapTrueKAsP_, meEndcapTrueKAsPi_, meEndcapTrueKNoPID_, meEndcapTruePAsK_, meEndcapTruePAsP_, meEndcapTruePAsPi_, meEndcapTruePiAsK_, meEndcapTruePiAsP_, meEndcapTruePiAsPi_, meEndcapTruePiNoPID_, meEndcapTruePNoPID_, meFakeTrackRecLVRelMult_, meFakeTrackRecLVRelMultvsMult_, meFakeTrackRecLVRelSumWos_, meFakeTrackRecLVRelSumWosvsSumWos_, meFakeTrackRelMult_, meFakeTrackRelMultvsMult_, meFakeTrackRelSumWos_, meFakeTrackRelSumWosvsSumWos_, meJetsFakeRelSumPt2_, meJetsFakeRelSumPt2vsSumPt2_, meJetsPUHt_, meJetsPUMetPt_, meJetsPUMult_, meJetsPURelHt_, meJetsPURelHtvsHt_, meJetsPURelMetPt_, meJetsPURelMetPtvsMetPt_, meJetsPURelMult_, meJetsPURelMultvsMult_, meJetsPURelSumPt2_, meJetsPURelSumPt2vsSumPt2_, meJetsPURelSumPz_, meJetsPURelSumPzvsSumPz_, meJetsPUSumPt2_, meJetsPUSumPz_, meJetsRecLVFakeRelSumPt2_, meJetsRecLVFakeRelSumPt2vsSumPt2_, meJetsRecLVPUHt_, meJetsRecLVPUMetPt_, meJetsRecLVPUMult_, meJetsRecLVPURelHt_, meJetsRecLVPURelHtvsHt_, meJetsRecLVPURelMetPt_, meJetsRecLVPURelMetPtvsMetPt_, meJetsRecLVPURelMult_, meJetsRecLVPURelMultvsMult_, meJetsRecLVPURelSumPt2_, meJetsRecLVPURelSumPt2vsSumPt2_, meJetsRecLVPURelSumPz_, meJetsRecLVPURelSumPzvsSumPz_, meJetsRecLVPUSumPt2_, meJetsRecLVPUSumPz_, meMatchQual_, meMVATrackEffEtaTot_, meMVATrackEffPtTot_, meMVATrackMatchedEffEtaMtd_, meMVATrackMatchedEffEtaTot_, meMVATrackMatchedEffPtMtd_, meMVATrackMatchedEffPtTot_, meMVATrackPullTot_, meMVATrackResTot_, meMVATrackZposResTot_, mePUTrackMult_, mePUTrackRecLVMult_, mePUTrackRecLVRelMult_, mePUTrackRecLVRelMultvsMult_, mePUTrackRecLVRelSumPt2_, mePUTrackRecLVRelSumPt2vsSumPt2_, mePUTrackRecLVRelSumPt_, mePUTrackRecLVRelSumPtvsSumPt_, mePUTrackRecLVRelSumWnt_, mePUTrackRecLVRelSumWntvsSumWnt_, mePUTrackRecLVRelSumWos_, mePUTrackRecLVRelSumWosvsSumWos_, mePUTrackRecLVSumPt2_, mePUTrackRecLVSumPt_, mePUTrackRecLVSumWnt_, mePUTrackRecLVSumWos_, mePUTrackRecLVWnt_, mePUTrackRelMult_, mePUTrackRelMultvsMult_, mePUTrackRelSumPt2_, mePUTrackRelSumPt2vsSumPt2_, mePUTrackRelSumPt_, mePUTrackRelSumPtvsSumPt_, mePUTrackRelSumWnt_, mePUTrackRelSumWntvsSumWnt_, mePUTrackRelSumWos_, mePUTrackRelSumWosvsSumWos_, mePUTrackSumPt2_, mePUTrackSumPt_, mePUTrackSumWnt_, mePUTrackSumWos_, mePUTrackWnt_, mePUvsOtherFakeV_, mePUvsRealV_, mePUvsSplitV_, meRecoPosInRecoOrigCollection_, meRecoPosInSimCollection_, meRecoPVPosSignal_, meRecoPVPosSignalNotHighestPt_, meRecoVtxVsLineDensity_, meRecPVT_, meRecPVZ_, meRecVerNumber_, meSecTrackRecLVRelSumWos_, meSecTrackRecLVRelSumWosvsSumWos_, meSecTrackRecLVSumWos_, meSecTrackRelSumWos_, meSecTrackRelSumWosvsSumWos_, meSecTrackSumWos_, meSimPosInSimOrigCollection_, meSimPVZ_, meTimePull_, meTimeRes_, meTimeSignalPull_, meTimeSignalRes_, meTrack3DposRes_, meTrackPull_, meTrackPullHighP_, meTrackPullHighPTot_, meTrackPullLowP_, meTrackPullLowPTot_, meTrackPullTot_, meTrackRes_, meTrackResHighP_, meTrackResHighPTot_, meTrackResLowP_, meTrackResLowPTot_, meTrackResMass_, meTrackResMassPions_, meTrackResMassProtons_, meTrackResMassTrue_, meTrackResMassTruePions_, meTrackResMassTrueProtons_, meTrackResTot_, meTrackZposRes_, meTrackZposResTot_, optionalPlots_, and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                                            {
  ibook.setCurrentFolder(folder_);
  // --- histograms booking
  meMVATrackEffPtTot_ = ibook.book1D("MVAEffPtTot", "Pt of tracks associated to LV; track pt [GeV] ", 110, 0., 11.);
  meMVATrackEffEtaTot_ = ibook.book1D("MVAEffEtaTot", "Eta of tracks associated to LV; track eta ", 66, 0., 3.3);
  meMVATrackMatchedEffPtTot_ =
      ibook.book1D("MVAMatchedEffPtTot", "Pt of tracks associated to LV matched to TP; track pt [GeV] ", 110, 0., 11.);
  meMVATrackMatchedEffPtMtd_ = ibook.book1D(
      "MVAMatchedEffPtMtd", "Pt of tracks associated to LV matched to TP with time; track pt [GeV] ", 110, 0., 11.);
  meMVATrackMatchedEffEtaTot_ =
      ibook.book1D("MVAMatchedEffEtaTot", "Eta of tracks associated to LV matched to TP; track eta ", 66, 0., 3.3);
  meMVATrackMatchedEffEtaMtd_ = ibook.book1D(
      "MVAMatchedEffEtaMtd", "Eta of tracks associated to LV matched to TP with time; track eta ", 66, 0., 3.3);
  meMVATrackResTot_ = ibook.book1D(
      "MVATrackRes", "t_{rec} - t_{sim} for tracks from LV MVA sel.; t_{rec} - t_{sim} [ns] ", 120, -0.15, 0.15);
  meTrackResTot_ = ibook.book1D("TrackRes", "t_{rec} - t_{sim} for tracks; t_{rec} - t_{sim} [ns] ", 120, -0.15, 0.15);
  meTrackRes_[0] = ibook.book1D(
      "TrackRes-LowMVA", "t_{rec} - t_{sim} for tracks with MVA < 0.5; t_{rec} - t_{sim} [ns] ", 100, -1., 1.);
  meTrackRes_[1] = ibook.book1D(
      "TrackRes-MediumMVA", "t_{rec} - t_{sim} for tracks with 0.5 < MVA < 0.8; t_{rec} - t_{sim} [ns] ", 100, -1., 1.);
  meTrackRes_[2] = ibook.book1D(
      "TrackRes-HighMVA", "t_{rec} - t_{sim} for tracks with MVA > 0.8; t_{rec} - t_{sim} [ns] ", 100, -1., 1.);
  if (optionalPlots_) {
    meTrackResMass_[0] = ibook.book1D(
        "TrackResMass-LowMVA", "t_{rec} - t_{est} for tracks with MVA < 0.5; t_{rec} - t_{est} [ns] ", 100, -1., 1.);
    meTrackResMass_[1] = ibook.book1D("TrackResMass-MediumMVA",
                                      "t_{rec} - t_{est} for tracks with 0.5 < MVA < 0.8; t_{rec} - t_{est} [ns] ",
                                      100,
                                      -1.,
                                      1.);
    meTrackResMass_[2] = ibook.book1D(
        "TrackResMass-HighMVA", "t_{rec} - t_{est} for tracks with MVA > 0.8; t_{rec} - t_{est} [ns] ", 100, -1., 1.);
    meTrackResMassTrue_[0] = ibook.book1D(
        "TrackResMassTrue-LowMVA", "t_{est} - t_{sim} for tracks with MVA < 0.5; t_{est} - t_{sim} [ns] ", 100, -1., 1.);
    meTrackResMassTrue_[1] = ibook.book1D("TrackResMassTrue-MediumMVA",
                                          "t_{est} - t_{sim} for tracks with 0.5 < MVA < 0.8; t_{est} - t_{sim} [ns] ",
                                          100,
                                          -1.,
                                          1.);
    meTrackResMassTrue_[2] = ibook.book1D("TrackResMassTrue-HighMVA",
                                          "t_{est} - t_{sim} for tracks with MVA > 0.8; t_{est} - t_{sim} [ns] ",
                                          100,
                                          -1.,
                                          1.);
  }
  meMVATrackPullTot_ =
      ibook.book1D("MVATrackPull", "Pull for tracks from LV MAV sel.; (t_{rec}-t_{sim})/#sigma_{t}", 50, -5., 5.);
  meTrackPullTot_ = ibook.book1D("TrackPull", "Pull for tracks; (t_{rec}-t_{sim})/#sigma_{t}", 100, -10., 10.);
  meTrackPull_[0] =
      ibook.book1D("TrackPull-LowMVA", "Pull for tracks with MVA < 0.5; (t_{rec}-t_{sim})/#sigma_{t}", 100, -10., 10.);
  meTrackPull_[1] = ibook.book1D(
      "TrackPull-MediumMVA", "Pull for tracks with 0.5 < MVA < 0.8; (t_{rec}-t_{sim})/#sigma_{t}", 100, -10., 10.);
  meTrackPull_[2] =
      ibook.book1D("TrackPull-HighMVA", "Pull for tracks with MVA > 0.8; (t_{rec}-t_{sim})/#sigma_{t}", 100, -10., 10.);
  meMVATrackZposResTot_ = ibook.book1D(
      "MVATrackZposResTot", "Z_{PCA} - Z_{sim} for tracks from LV MVA sel.;Z_{PCA} - Z_{sim} [cm] ", 50, -0.1, 0.1);
  meTrackZposResTot_ =
      ibook.book1D("TrackZposResTot", "Z_{PCA} - Z_{sim} for tracks;Z_{PCA} - Z_{sim} [cm] ", 50, -0.5, 0.5);
  meTrackZposRes_[0] = ibook.book1D(
      "TrackZposRes-LowMVA", "Z_{PCA} - Z_{sim} for tracks with MVA < 0.5;Z_{PCA} - Z_{sim} [cm] ", 50, -0.5, 0.5);
  meTrackZposRes_[1] = ibook.book1D("TrackZposRes-MediumMVA",
                                    "Z_{PCA} - Z_{sim} for tracks with 0.5 < MVA < 0.8 ;Z_{PCA} - Z_{sim} [cm] ",
                                    50,
                                    -0.5,
                                    0.5);
  meTrackZposRes_[2] = ibook.book1D(
      "TrackZposRes-HighMVA", "Z_{PCA} - Z_{sim} for tracks with MVA > 0.8 ;Z_{PCA} - Z_{sim} [cm] ", 50, -0.5, 0.5);
  meTrack3DposRes_[0] =
      ibook.book1D("Track3DposRes-LowMVA",
                   "3dPos_{PCA} - 3dPos_{sim} for tracks with MVA < 0.5 ;3dPos_{PCA} - 3dPos_{sim} [cm] ",
                   50,
                   -0.5,
                   0.5);
  meTrack3DposRes_[1] =
      ibook.book1D("Track3DposRes-MediumMVA",
                   "3dPos_{PCA} - 3dPos_{sim} for tracks with 0.5 < MVA < 0.8 ;3dPos_{PCA} - 3dPos_{sim} [cm] ",
                   50,
                   -0.5,
                   0.5);
  meTrack3DposRes_[2] =
      ibook.book1D("Track3DposRes-HighMVA",
                   "3dPos_{PCA} - 3dPos_{sim} for tracks with MVA > 0.8;3dPos_{PCA} - 3dPos_{sim} [cm] ",
                   50,
                   -0.5,
                   0.5);
  meTimeRes_ = ibook.book1D("TimeRes", "t_{rec} - t_{sim} ;t_{rec} - t_{sim} [ns] ", 40, -0.2, 0.2);
  meTimePull_ = ibook.book1D("TimePull", "Pull; t_{rec} - t_{sim}/#sigma_{t rec}", 100, -10., 10.);
  meTimeSignalRes_ =
      ibook.book1D("TimeSignalRes", "t_{rec} - t_{sim} for signal ;t_{rec} - t_{sim} [ns] ", 40, -0.2, 0.2);
  meTimeSignalPull_ =
      ibook.book1D("TimeSignalPull", "Pull for signal; t_{rec} - t_{sim}/#sigma_{t rec}", 100, -10., 10.);
  mePUvsRealV_ =
      ibook.bookProfile("PUvsReal", "#PU vertices vs #real matched vertices;#PU;#real ", 100, 0, 300, 100, 0, 200);
  mePUvsOtherFakeV_ = ibook.bookProfile(
      "PUvsOtherFake", "#PU vertices vs #other fake matched vertices;#PU;#other fake ", 100, 0, 300, 100, 0, 20);
  mePUvsSplitV_ =
      ibook.bookProfile("PUvsSplit", "#PU vertices vs #split matched vertices;#PU;#split ", 100, 0, 300, 100, 0, 20);
  meMatchQual_ = ibook.book1D("MatchQuality", "RECO-SIM vertex match quality; ", 8, 0, 8.);
  meDeltaZrealreal_ = ibook.book1D("DeltaZrealreal", "#Delta Z real-real; |#Delta Z (r-r)| [cm]", 100, 0, 0.5);
  meDeltaZfakefake_ = ibook.book1D("DeltaZfakefake", "#Delta Z fake-fake; |#Delta Z (f-f)| [cm]", 100, 0, 0.5);
  meDeltaZfakereal_ = ibook.book1D("DeltaZfakereal", "#Delta Z fake-real; |#Delta Z (f-r)| [cm]", 100, 0, 0.5);
  meDeltaTrealreal_ = ibook.book1D("DeltaTrealreal", "#Delta T real-real; |#Delta T (r-r)| [sigma]", 60, 0., 30.);
  meDeltaTfakefake_ = ibook.book1D("DeltaTfakefake", "#Delta T fake-fake; |#Delta T (f-f)| [sigma]", 60, 0., 30.);
  meDeltaTfakereal_ = ibook.book1D("DeltaTfakereal", "#Delta T fake-real; |#Delta T (f-r)| [sigma]", 60, 0., 30.);
  if (optionalPlots_) {
    meRecoPosInSimCollection_ = ibook.book1D(
        "RecoPosInSimCollection", "Sim signal vertex index associated to Reco signal vertex; Sim PV index", 200, 0, 200);
    meRecoPosInRecoOrigCollection_ =
        ibook.book1D("RecoPosInRecoOrigCollection", "Reco signal index in OrigCollection; Reco index", 200, 0, 200);
    meSimPosInSimOrigCollection_ =
        ibook.book1D("SimPosInSimOrigCollection", "Sim signal index in OrigCollection; Sim index", 200, 0, 200);
  }
  meRecoPVPosSignal_ =
      ibook.book1D("RecoPVPosSignal", "Position in reco collection of PV associated to sim signal", 20, 0, 20);
  meRecoPVPosSignalNotHighestPt_ =
      ibook.book1D("RecoPVPosSignalNotHighestPt",
                   "Position in reco collection of PV associated to sim signal not highest Pt",
                   20,
                   0,
                   20);
  meRecoVtxVsLineDensity_ =
      ibook.book1D("RecoVtxVsLineDensity", "#Reco vertices/mm/event; line density [#vtx/mm/event]", 160, 0., 4.);
  meRecVerNumber_ = ibook.book1D("RecVerNumber", "RECO Vertex Number: Number of vertices", 50, 0, 250);
  meRecPVZ_ = ibook.book1D("recPVZ", "Weighted #Rec vertices/mm", 400, -20., 20.);
  meRecPVT_ = ibook.book1D("recPVT", "#Rec vertices/10 ps", 200, -1., 1.);
  meSimPVZ_ = ibook.book1D("simPVZ", "Weighted #Sim vertices/mm", 400, -20., 20.);

  mePUTrackRelMult_ = ibook.book1D(
      "PUTrackRelMult", "Relative multiplicity of PU tracks for matched vertices; #PUTrks/#Trks", 50, 0., 1.);
  meFakeTrackRelMult_ = ibook.book1D(
      "FakeTrackRelMult", "Relative multiplicity of fake tracks for matched vertices; #fakeTrks/#Trks", 50, 0., 1.);
  mePUTrackRelSumWnt_ =
      ibook.book1D("PUTrackRelSumWnt",
                   "Relative Sum of wnt of PU tracks for matched vertices; PUSumW*min(Pt, 1.)/SumW*min(Pt, 1.)",
                   50,
                   0.,
                   1.);
  mePUTrackRelSumWos_ = ibook.book1D(
      "PUTrackRelSumWos", "Relative Sum of wos of PU tracks for matched vertices; PUSumWos/SumWos", 50, 0., 1.);
  meSecTrackRelSumWos_ =
      ibook.book1D("SecTrackRelSumWos",
                   "Relative Sum of wos of tracks from secondary vtx for matched vertices; SecSumWos/SumWos",
                   50,
                   0.,
                   1.);
  meFakeTrackRelSumWos_ = ibook.book1D(
      "FakeTrackRelSumWos", "Relative Sum of wos of fake tracks for matched vertices; FakeSumWos/SumWos", 50, 0., 1.);
  mePUTrackRelSumPt_ = ibook.book1D(
      "PUTrackRelSumPt", "Relative Sum of Pt of PU tracks for matched vertices; PUSumPt/SumPt", 50, 0., 1.);
  mePUTrackRelSumPt2_ = ibook.book1D(
      "PUTrackRelSumPt2", "Relative Sum of Pt2 for PU tracks for matched vertices; PUSumPt2/SumPt2", 50, 0., 1.);

  mePUTrackRecLVRelMult_ = ibook.book1D(
      "PUTrackRecLVRelMult", "Relative multiplicity of PU tracks for matched LV; #PUTrks/#Trks", 50, 0., 1.);
  meFakeTrackRecLVRelMult_ = ibook.book1D(
      "FakeTrackRecLVRelMult", "Relative multiplicity of fake tracks for matched LV; #FakeTrks/#Trks", 50, 0., 1.);
  mePUTrackRecLVRelSumWnt_ =
      ibook.book1D("PUTrackRecLVRelSumWnt",
                   "Relative Sum of Wnt of PU tracks for matched LV; PUSumW*min(Pt, 1.)/SumW*min(Pt, 1.)",
                   50,
                   0.,
                   1.);
  mePUTrackRecLVRelSumWos_ = ibook.book1D(
      "PUTrackRecLVRelSumWos", "Relative Sum of Wos of PU tracks for matched LV; PUSumWos/SumWos", 50, 0., 1.);
  meSecTrackRecLVRelSumWos_ =
      ibook.book1D("SecTrackRecLVRelSumWos",
                   "Relative Sum of wos of tracks from secondary vtx for matched LV; SecSumWos/SumWos",
                   50,
                   0.,
                   1.);
  meFakeTrackRecLVRelSumWos_ = ibook.book1D(
      "FakeTrackRecLVRelSumWos", "Relative Sum of wos of fake tracks for matched LV; FakeSumWos/SumWos", 50, 0., 1.);
  mePUTrackRecLVRelSumPt_ =
      ibook.book1D("PUTrackRecLVRelSumPt", "Relative Sum of Pt of PU tracks for matched LV; PUSumPt/SumPt", 50, 0., 1.);
  mePUTrackRecLVRelSumPt2_ = ibook.book1D(
      "PUTrackRecLVRelSumPt2", "Relative Sum of Pt2 of PU tracks for matched LV; PUSumPt2/SumPt2", 50, 0., 1.);

  if (optionalPlots_) {
    mePUTrackMult_ = ibook.book1D("PUTrackMult", "Number of PU tracks for matched vertices; #PUTrks", 50, 0., 100.);
    mePUTrackWnt_ = ibook.book1D("PUTrackWnt", "Wnt of PU tracks for matched vertices; PUTrkW*min(Pt, 1.)", 50, 0., 1.);
    mePUTrackSumWnt_ = ibook.book1D(
        "PUTrackSumWnt", "Sum of wnt of PU tracks for matched vertices; log10(PUSumW*min(Pt, 1.))", 50, -2., 3.);
    mePUTrackSumWos_ =
        ibook.book1D("PUTrackSumWos", "Sum of wos of PU tracks for matched vertices; log10(PUSumWos)", 50, -1., 7.);
    meSecTrackSumWos_ = ibook.book1D(
        "SecTrackSumWos", "Sum of wos of tracks from secondary vtx for matched vertices; log10(SecSumWos)", 50, -1., 7.);
    mePUTrackSumPt_ =
        ibook.book1D("PUTrackSumPt", "Sum of Pt of PU tracks for matched vertices; log10(PUSumPt)", 50, -2., 3.);
    mePUTrackSumPt2_ =
        ibook.book1D("PUTrackSumPt2", "Sum of Pt2 of PU tracks for matched vertices; log10(PUSumPt2)", 50, -3., 3.);

    mePUTrackRelMultvsMult_ =
        ibook.bookProfile("PUTrackRelMultvsMult",
                          "Relative PU multiplicity vs Number of tracks for matched vertices; #Trks; #PUTrks/#Trks",
                          50,
                          0.,
                          200.,
                          0.,
                          1.,
                          "s");
    meFakeTrackRelMultvsMult_ = ibook.bookProfile(
        "FakeTrackRelMultvsMult",
        "Relative multiplicity of fake tracks vs Number of tracks for matched vertices; #Trks; #FakeTrks/#Trks",
        50,
        0.,
        200.,
        0.,
        1.,
        "s");
    mePUTrackRelSumWntvsSumWnt_ =
        ibook.bookProfile("PUTrackRelSumWntvsSumWnt",
                          "Relative PU Sum of Wnt vs Sum of Wnt of tracks for matched vertices; log10(SumW*min(Pt, "
                          "1.)); PUSumW*min(Pt, 1.)/SumW*min(Pt, 1.)",
                          50,
                          0.,
                          2.5,
                          0.,
                          1.,
                          "s");
    mePUTrackRelSumWosvsSumWos_ = ibook.bookProfile(
        "PUTrackRelSumWosvsSumWos",
        "Relative PU Sum of Wos vs Sum of Wos of tracks for matched vertices; log10(SumWos); PUSumWos/SumWos",
        50,
        2.5,
        7.,
        0.,
        1.,
        "s");
    meSecTrackRelSumWosvsSumWos_ = ibook.bookProfile("SecTrackRelSumWosvsSumWos",
                                                     "Relative Sum of Wos of tracks from secondary vtx vs Sum of Wos "
                                                     "of tracks for matched vertices; log10(SumWos); SecSumWos/SumWos",
                                                     50,
                                                     2.,
                                                     7.,
                                                     0.,
                                                     1.,
                                                     "s");
    meFakeTrackRelSumWosvsSumWos_ = ibook.bookProfile("FakeTrackRelSumWosvsSumWos",
                                                      "Relative Sum of Wos of fake tracks vs Sum of Wos of tracks for "
                                                      "matched vertices; log10(SumWos); FakeSumWos/SumWos",
                                                      50,
                                                      2.5,
                                                      7.5,
                                                      0.,
                                                      1.,
                                                      "s");
    mePUTrackRelSumPtvsSumPt_ = ibook.bookProfile(
        "PUTrackRelSumPtvsSumPt",
        "Relative PU Sum of Pt vs Sum of Pt of tracks for matched vertices; log10(SumPt); PUSumPt/SumPt",
        50,
        0.,
        3.,
        0.,
        1.,
        "s");
    mePUTrackRelSumPt2vsSumPt2_ = ibook.bookProfile(
        "PUTrackRelSumPt2vsSumPt2",
        "Relative PU Sum of Pt2 vs Sum of Pt2 of tracks for matched vertices; log10(SumPt2); PUSumPt2/SumPt2",
        50,
        0.,
        4.,
        0.,
        1.,
        "s");

    mePUTrackRecLVMult_ = ibook.book1D("PUTrackRecLVMult", "Number of PU tracks for matched LV; #PUTrks", 50, 0., 100.);
    mePUTrackRecLVWnt_ =
        ibook.book1D("PUTrackRecLVWnt", "Wnt of PU tracks for matched LV; PUTrkW*min(Pt, 1.)", 50, 0., 1.);
    mePUTrackRecLVSumWnt_ = ibook.book1D(
        "PUTrackRecLVSumWnt", "Sum of wnt of PU tracks for matched LV; log10(PUSumW*min(Pt, 1.))", 50, -2., 3.);
    mePUTrackRecLVSumWos_ =
        ibook.book1D("PUTrackRecLVSumWos", "Sum of wos of PU tracks for matched LV; log10(PUSumWos)", 50, -1., 7.);
    meSecTrackRecLVSumWos_ = ibook.book1D(
        "SecTrackRecLVSumWos", "Sum of wos of tracks from secondary vtx for matched LV; log10(SecSumWos)", 50, -1., 7.);
    mePUTrackRecLVSumPt_ =
        ibook.book1D("PUTrackRecLVSumPt", "Sum of Pt of PU tracks for matched LV; log10(PUSumPt)", 50, -2., 3.);
    mePUTrackRecLVSumPt2_ =
        ibook.book1D("PUTrackRecLVSumPt2", "Sum of Pt2 of PU tracks for matched LV; log10(PUSumPt2)", 50, -3., 3.);

    mePUTrackRecLVRelMultvsMult_ =
        ibook.bookProfile("PUTrackRecLVRelMultvsMult",
                          "Relative PU multiplicity vs Number of tracks for matched LV; #Trks; #PUTrks/#Trks",
                          50,
                          0.,
                          200.,
                          0.,
                          1.,
                          "s");
    meFakeTrackRecLVRelMultvsMult_ = ibook.bookProfile(
        "FakeTrackRecLVRelMultvsMult",
        "Relative multiplicity of fake tracks  vs Number of tracks for matched LV; #Trks; #FakeTrks/#Trks",
        50,
        0.,
        200.,
        0.,
        1.,
        "s");
    mePUTrackRecLVRelSumWntvsSumWnt_ =
        ibook.bookProfile("PUTrackRecLVRelSumWntvsSumWnt",
                          "Relative PU Sum of Wnt vs Sum of Wnt of tracks for matched LV; log10(SumW*min(Pt, 1.)); "
                          "PUSumW*min(Pt, 1.)/SumW*min(Pt, 1.)",
                          50,
                          1.,
                          2.3,
                          0.,
                          1.,
                          "s");
    mePUTrackRecLVRelSumWosvsSumWos_ = ibook.bookProfile(
        "PUTrackRecLVRelSumWosvsSumWos",
        "Relative PU Sum of Wos vs Sum of Wos of tracks for matched vertices; log10(SumWos); PUSumWos/SumWos",
        50,
        5.5,
        7.,
        0.,
        1.,
        "s");
    meSecTrackRecLVRelSumWosvsSumWos_ =
        ibook.bookProfile("SecTrackRecLVRelSumWosvsSumWos",
                          "Relative Sum of Wos of tracks from secondary vtx vs Sum of Wos of tracks for matched "
                          "vertices; log10(SumWos); SecSumWos/SumWos",
                          50,
                          5.5,
                          7.,
                          0.,
                          1.,
                          "s");
    meFakeTrackRecLVRelSumWosvsSumWos_ = ibook.bookProfile("FakeTrackRecLVRelSumWosvsSumWos",
                                                           "Relative Sum of Wos of fake tracks vs Sum of Wos of tracks "
                                                           "for matched vertices; log10(SumWos); FakeSumWos/SumWos",
                                                           50,
                                                           5.5,
                                                           7.,
                                                           0.,
                                                           1.,
                                                           "s");
    mePUTrackRecLVRelSumPtvsSumPt_ =
        ibook.bookProfile("PUTrackRecLVRelSumPtvsSumPt",
                          "Relative PU Sum of Pt vs Sum of Pt of tracks for matched LV; log10(SumPt); PUSumPt/SumPt",
                          50,
                          1.4,
                          3.,
                          0.,
                          1.,
                          "s");
    mePUTrackRecLVRelSumPt2vsSumPt2_ =
        ibook.bookProfile("PUTrackRecLVRelSumPt2vsSumPt2",
                          "Relative PU Sum of Pt2 vs Sum of tracks for matched LV; log10(SumPt2); PUSumPt2/SumPt2",
                          50,
                          2.,
                          4.,
                          0.,
                          1.,
                          "s");
  }

  meJetsPURelMult_ =
      ibook.book1D("JetsPURelMult",
                   "Relative contribution of PU to jet multiplicity for matched vertices; (#Jets-#JetsNoPU)/#Jets",
                   50,
                   0.,
                   1.);
  meJetsPURelHt_ =
      ibook.book1D("JetsPURelHt",
                   "Relative contribution of PU to scalar sum of Et of jets for matched vertices; (Ht-HtNoPU)/HT",
                   50,
                   0.,
                   1.);
  meJetsPURelSumPt2_ =
      ibook.book1D("JetsPURelSumPt2",
                   "Relative contribution of PU to sum of Pt2 of jets for matched vertices; (SumPt2-SumPt2NoPU)/SumPt2",
                   50,
                   0.,
                   1.);
  meJetsFakeRelSumPt2_ = ibook.book1D(
      "JetsFakeRelSumPt2",
      "Relative contribution of fake tracks to sum of Pt2 of jets for matched vertices; (SumPt2-SumPt2NoFake)/SumPt2",
      50,
      0.,
      1.);
  meJetsPURelMetPt_ =
      ibook.book1D("JetsPURelMetPt",
                   "Relative contribution of PU to Missing Transverse Energy for matched vertices; (Met-MetNoPU)/Met",
                   50,
                   -1.,
                   1.);
  meJetsPURelSumPz_ =
      ibook.book1D("JetsPURelSumPz",
                   "Relative contribution of PU to sum of Pz of jets for matched vertices; (SumPz-SumPzNoPU)/SumPz",
                   50,
                   -1.,
                   1.);

  meJetsRecLVPURelMult_ =
      ibook.book1D("JetsRecLVPURelMult",
                   "Relative contribution of PU to jet multiplicity for matched LV; (#Jets-#JetsNoPU)/#Jets",
                   50,
                   0.,
                   1.);
  meJetsRecLVPURelHt_ =
      ibook.book1D("JetsRecLVPURelHt",
                   "Relative contribution of PU to scalar sum of Et of jets for matched LV; (Ht-HtNoPU)/HT",
                   50,
                   0.,
                   1.);
  meJetsRecLVPURelSumPt2_ =
      ibook.book1D("JetsRecLVPURelSumPt2",
                   "Relative contribution of PU to sum of Pt2 of jets for matched LV; (SumPt2-SumPt2NoPU)/SumPt2",
                   50,
                   0.,
                   1.);
  meJetsRecLVFakeRelSumPt2_ = ibook.book1D(
      "JetsRecLVFakeRelSumPt2",
      "Relative contribution of fake tracks to sum of Pt2 of jets for matched LV; (SumPt2-SumPt2NoFake)/SumPt2",
      50,
      0.,
      1.);
  meJetsRecLVPURelMetPt_ =
      ibook.book1D("JetsRecLVPURelMetPt",
                   "Relative contribution of PU to Missing Transverse Energy for matched LV; (Met-MetNoPU)/Met",
                   50,
                   -1.,
                   1.);
  meJetsRecLVPURelSumPz_ =
      ibook.book1D("JetsRecLVPURelSumPz",
                   "Relative contribution of PU to sum of Pz of jets for matched LV; (SumPz-SumPzNoPU)/SumPz",
                   50,
                   -1.,
                   1.);

  if (optionalPlots_) {
    meJetsPUMult_ = ibook.book1D(
        "JetsPUMult", "Contribution of PU to jet multiplicity for matched vertices; #Jets-#JetsNoPU", 50, 0., 100.);
    meJetsPUHt_ = ibook.book1D("JetsPUHt",
                               "Contribution of PU to scalar sum of Et of jets for matched vertices; log10(Ht-HtNoPU)",
                               50,
                               -2.,
                               3.);
    meJetsPUSumPt2_ =
        ibook.book1D("JetsPUSumPt2",
                     "Contribution of PU to sum of Pt2 of jets for matched vertices; log10(sumPt2-SumPt2NoPU)",
                     50,
                     -3.,
                     3.);
    meJetsPUMetPt_ =
        ibook.book1D("JetsPUMetPt",
                     "Contribution of PU to Missing Transverse Energy for matched vertices; log10(Met-MetNoPU)",
                     50,
                     -3.,
                     2.);
    meJetsPUSumPz_ =
        ibook.book1D("JetsPUSumPz",
                     "Contribution of PU to sum of Pz of jets for matched vertices; log10(abs(SumPz-SumPzNoPU))",
                     50,
                     -4.,
                     3.);

    meJetsPURelMultvsMult_ = ibook.bookProfile("JetsPURelMultvsMult",
                                               "Relative contribution of PU to jet multiplicity vs number of jets for "
                                               "matched vertices; #Jets; (#Jets-#JetsNoPU)/#Jets",
                                               50,
                                               0.,
                                               120.,
                                               0.,
                                               1.,
                                               "s");
    meJetsPURelHtvsHt_ = ibook.bookProfile("JetsPURelHtvsHt",
                                           "Relative contribution of PU to scalar sum of Et of jets vs scalar sum of "
                                           "Et for matched vertices; log10(Ht); (Ht-HtNoPU)/HT",
                                           50,
                                           0.,
                                           3.,
                                           0.,
                                           1.,
                                           "s");
    meJetsPURelSumPt2vsSumPt2_ = ibook.bookProfile("JetsPURelSumPt2vsSumPt2",
                                                   "Relative contribution of PU to sum of Pt2 of jets vs sum of Pt2 "
                                                   "for matched vertices; log10(SumPt2); (SumPt2-SumPt2NoPU)/SumPt2",
                                                   50,
                                                   -1.,
                                                   4.,
                                                   0.,
                                                   1.,
                                                   "s");
    meJetsFakeRelSumPt2vsSumPt2_ =
        ibook.bookProfile("JetsFakeRelSumPt2vsSumPt2",
                          "Relative contribution of fake tracks to sum of Pt2 of jets vs sum of Pt2 for matched "
                          "vertices; log10(SumPt2); (SumPt2-SumPt2NoFake)/SumPt2",
                          50,
                          -1.,
                          4.,
                          0.,
                          1.,
                          "s");
    meJetsPURelMetPtvsMetPt_ = ibook.bookProfile("JetsPURelMetPtvsMetPt",
                                                 "Relative contribution of PU to Missing Transverse Energy vs MET for "
                                                 "matched vertices; log10(Met); (Met-MetNoPU)/Met",
                                                 50,
                                                 -1.,
                                                 2.,
                                                 -1.,
                                                 1.,
                                                 "s");
    meJetsPURelSumPzvsSumPz_ = ibook.bookProfile("JetsPURelSumPzvsSumPz",
                                                 "Relative contribution of PU to sum of Pz of jets vs Sum of Pz for "
                                                 "matched vertices; log10(abs SumPz); (SumPz-SumPzNoPU)/SumPz",
                                                 50,
                                                 -2.,
                                                 3.,
                                                 -1.,
                                                 1.,
                                                 "s");

    meJetsRecLVPUMult_ = ibook.book1D(
        "JetsRecLVPUMult", "Contribution of PU to jet multiplicity for matched LV; #Jets-#JetsNoPU", 50, 0., 100.);
    meJetsRecLVPUHt_ = ibook.book1D(
        "JetsRecLVPUHt", "Contribution of PU to scalar sum of Et of jets for matched LV; log10(Ht-HtNoPU)", 50, -2., 3.);
    meJetsRecLVPUSumPt2_ =
        ibook.book1D("JetsRecLVPUSumPt2",
                     "Contribution of PU to sum of Pt2 of jets for matched LV; log10(sumPt2-SumPt2NoPU)",
                     50,
                     -3.,
                     3.);
    meJetsRecLVPUMetPt_ =
        ibook.book1D("JetsRecLVPUMetPt",
                     "Contribution of PU to Missing Transverse Energy for matched LV; log10(Met-MetNoPU)",
                     50,
                     -3.,
                     2.);
    meJetsRecLVPUSumPz_ =
        ibook.book1D("JetsRecLVPUSumPz",
                     "Contribution of PU to sum of Pz of jets for matched LV; log10(abs(SumPz-SumPzNoPU))",
                     50,
                     -4.,
                     3.);

    meJetsRecLVPURelMultvsMult_ = ibook.bookProfile("JetsRecLVPURelMultvsMult",
                                                    "Relative contribution of PU to jet multiplicity vs number of jets "
                                                    "for matched vertices; #Jets; (#Jets-#JetsNoPU)/#Jets",
                                                    50,
                                                    0.,
                                                    120.,
                                                    0.,
                                                    1.,
                                                    "s");
    meJetsRecLVPURelHtvsHt_ = ibook.bookProfile("JetsRecLVPURelHtvsHt",
                                                "Relative contribution of PU to scalar sum of Et of jets vs scalar sum "
                                                "of Et for matched vertices; log10(Ht); (Ht-HtNoPU)/HT",
                                                50,
                                                1.5,
                                                3.,
                                                0.,
                                                1.,
                                                "s");
    meJetsRecLVPURelSumPt2vsSumPt2_ =
        ibook.bookProfile("JetsRecLVPURelSumPt2vsSumPt2",
                          "Relative contribution of PU to sum of Pt2 of jets vs sum of Pt2 for matched vertices; "
                          "log10(SumPt2); (SumPt2-SumPt2NoPU)/SumPt2",
                          50,
                          2.,
                          5.,
                          0.,
                          1.,
                          "s");
    meJetsRecLVFakeRelSumPt2vsSumPt2_ =
        ibook.bookProfile("JetsRecLVFakeRelSumPt2vsSumPt2",
                          "Relative contribution of fake tracks to sum of Pt2 of jets vs sum of Pt2 for matched "
                          "vertices; log10(SumPt2); (SumPt2-SumPt2NoFake)/SumPt2",
                          50,
                          2.,
                          5.,
                          0.,
                          1.,
                          "s");
    meJetsRecLVPURelMetPtvsMetPt_ = ibook.bookProfile("JetsRecLVPURelMetPtvsMetPt",
                                                      "Relative contribution of PU to Missing Transverse Energy vs MET "
                                                      "for matched vertices; log10(Met); (Met-MetNoPU)/Met",
                                                      50,
                                                      0.,
                                                      2.5,
                                                      -1.,
                                                      1.,
                                                      "s");
    meJetsRecLVPURelSumPzvsSumPz_ = ibook.bookProfile("JetsRecLVPURelSumPzvsSumPz",
                                                      "Relative contribution of PU to sum of Pz of jets vs Sum of Pz "
                                                      "for matched vertices; log10(abs SumPz); (SumPz-SumPzNoPU)/SumPz",
                                                      50,
                                                      0.5,
                                                      3.5,
                                                      -1.,
                                                      1.,
                                                      "s");
  }

  // some tests
  meTrackResLowPTot_ = ibook.book1D(
      "TrackResLowP", "t_{rec} - t_{sim} for tracks with p < 2 GeV; t_{rec} - t_{sim} [ns] ", 70, -0.15, 0.15);
  meTrackResLowP_[0] =
      ibook.book1D("TrackResLowP-LowMVA",
                   "t_{rec} - t_{sim} for tracks with MVA < 0.5 and p < 2 GeV; t_{rec} - t_{sim} [ns] ",
                   100,
                   -1.,
                   1.);
  meTrackResLowP_[1] =
      ibook.book1D("TrackResLowP-MediumMVA",
                   "t_{rec} - t_{sim} for tracks with 0.5 < MVA < 0.8 and p < 2 GeV; t_{rec} - t_{sim} [ns] ",
                   100,
                   -1.,
                   1.);
  meTrackResLowP_[2] =
      ibook.book1D("TrackResLowP-HighMVA",
                   "t_{rec} - t_{sim} for tracks with MVA > 0.8 and p < 2 GeV; t_{rec} - t_{sim} [ns] ",
                   100,
                   -1.,
                   1.);
  meTrackResHighPTot_ = ibook.book1D(
      "TrackResHighP", "t_{rec} - t_{sim} for tracks with p > 2 GeV; t_{rec} - t_{sim} [ns] ", 70, -0.15, 0.15);
  meTrackResHighP_[0] =
      ibook.book1D("TrackResHighP-LowMVA",
                   "t_{rec} - t_{sim} for tracks with MVA < 0.5 and p > 2 GeV; t_{rec} - t_{sim} [ns] ",
                   100,
                   -1.,
                   1.);
  meTrackResHighP_[1] =
      ibook.book1D("TrackResHighP-MediumMVA",
                   "t_{rec} - t_{sim} for tracks with 0.5 < MVA < 0.8 and p > 2 GeV; t_{rec} - t_{sim} [ns] ",
                   100,
                   -1.,
                   1.);
  meTrackResHighP_[2] =
      ibook.book1D("TrackResHighP-HighMVA",
                   "t_{rec} - t_{sim} for tracks with MVA > 0.8 and p > 2 GeV; t_{rec} - t_{sim} [ns] ",
                   100,
                   -1.,
                   1.);
  meTrackPullLowPTot_ =
      ibook.book1D("TrackPullLowP", "Pull for tracks with p < 2 GeV; (t_{rec}-t_{sim})/#sigma_{t}", 100, -10., 10.);
  meTrackPullLowP_[0] = ibook.book1D("TrackPullLowP-LowMVA",
                                     "Pull for tracks with MVA < 0.5 and p < 2 GeV; (t_{rec}-t_{sim})/#sigma_{t}",
                                     100,
                                     -10.,
                                     10.);
  meTrackPullLowP_[1] = ibook.book1D("TrackPullLowP-MediumMVA",
                                     "Pull for tracks with 0.5 < MVA < 0.8 and p < 2 GeV; (t_{rec}-t_{sim})/#sigma_{t}",
                                     100,
                                     -10.,
                                     10.);
  meTrackPullLowP_[2] = ibook.book1D("TrackPullLowP-HighMVA",
                                     "Pull for tracks with MVA > 0.8 and p < 2 GeV; (t_{rec}-t_{sim})/#sigma_{t}",
                                     100,
                                     -10.,
                                     10.);
  meTrackPullHighPTot_ =
      ibook.book1D("TrackPullHighP", "Pull for tracks with p > 2 GeV; (t_{rec}-t_{sim})/#sigma_{t}", 100, -10., 10.);
  meTrackPullHighP_[0] = ibook.book1D("TrackPullHighP-LowMVA",
                                      "Pull for tracks with MVA < 0.5 and p > 2 GeV; (t_{rec}-t_{sim})/#sigma_{t}",
                                      100,
                                      -10.,
                                      10.);
  meTrackPullHighP_[1] =
      ibook.book1D("TrackPullHighP-MediumMVA",
                   "Pull for tracks with 0.5 < MVA < 0.8 and p > 2 GeV; (t_{rec}-t_{sim})/#sigma_{t}",
                   100,
                   -10.,
                   10.);
  meTrackPullHighP_[2] = ibook.book1D("TrackPullHighP-HighMVA",
                                      "Pull for tracks with MVA > 0.8 and p > 2 GeV; (t_{rec}-t_{sim})/#sigma_{t}",
                                      100,
                                      -10.,
                                      10.);
  if (optionalPlots_) {
    meTrackResMassProtons_[0] =
        ibook.book1D("TrackResMass-Protons-LowMVA",
                     "t_{rec} - t_{est} for proton tracks with MVA < 0.5; t_{rec} - t_{est} [ns] ",
                     100,
                     -1.,
                     1.);
    meTrackResMassProtons_[1] =
        ibook.book1D("TrackResMass-Protons-MediumMVA",
                     "t_{rec} - t_{est} for proton tracks with 0.5 < MVA < 0.8; t_{rec} - t_{est} [ns] ",
                     100,
                     -1.,
                     1.);
    meTrackResMassProtons_[2] =
        ibook.book1D("TrackResMass-Protons-HighMVA",
                     "t_{rec} - t_{est} for proton tracks with MVA > 0.8; t_{rec} - t_{est} [ns] ",
                     100,
                     -1.,
                     1.);
    meTrackResMassTrueProtons_[0] =
        ibook.book1D("TrackResMassTrue-Protons-LowMVA",
                     "t_{est} - t_{sim} for proton tracks with MVA < 0.5; t_{est} - t_{sim} [ns] ",
                     100,
                     -1.,
                     1.);
    meTrackResMassTrueProtons_[1] =
        ibook.book1D("TrackResMassTrue-Protons-MediumMVA",
                     "t_{est} - t_{sim} for proton tracks with 0.5 < MVA < 0.8; t_{est} - t_{sim} [ns] ",
                     100,
                     -1.,
                     1.);
    meTrackResMassTrueProtons_[2] =
        ibook.book1D("TrackResMassTrue-Protons-HighMVA",
                     "t_{est} - t_{sim} for proton tracks with MVA > 0.8; t_{est} - t_{sim} [ns] ",
                     100,
                     -1.,
                     1.);

    meTrackResMassPions_[0] = ibook.book1D("TrackResMass-Pions-LowMVA",
                                           "t_{rec} - t_{est} for pion tracks with MVA < 0.5; t_{rec} - t_{est} [ns] ",
                                           100,
                                           -1.,
                                           1.);
    meTrackResMassPions_[1] =
        ibook.book1D("TrackResMass-Pions-MediumMVA",
                     "t_{rec} - t_{est} for pion tracks with 0.5 < MVA < 0.8; t_{rec} - t_{est} [ns] ",
                     100,
                     -1.,
                     1.);
    meTrackResMassPions_[2] = ibook.book1D("TrackResMass-Pions-HighMVA",
                                           "t_{rec} - t_{est} for pion tracks with MVA > 0.8; t_{rec} - t_{est} [ns] ",
                                           100,
                                           -1.,
                                           1.);
    meTrackResMassTruePions_[0] =
        ibook.book1D("TrackResMassTrue-Pions-LowMVA",
                     "t_{est} - t_{sim} for pion tracks with MVA < 0.5; t_{est} - t_{sim} [ns] ",
                     100,
                     -1.,
                     1.);
    meTrackResMassTruePions_[1] =
        ibook.book1D("TrackResMassTrue-Pions-MediumMVA",
                     "t_{est} - t_{sim} for pion tracks with 0.5 < MVA < 0.8; t_{est} - t_{sim} [ns] ",
                     100,
                     -1.,
                     1.);
    meTrackResMassTruePions_[2] =
        ibook.book1D("TrackResMassTrue-Pions-HighMVA",
                     "t_{est} - t_{sim} for pion tracks with MVA > 0.8; t_{est} - t_{sim} [ns] ",
                     100,
                     -1.,
                     1.);
  }

  meBarrelPIDp_ = ibook.book1D("BarrelPIDp", "PID track MTD momentum spectrum, |eta| < 1.5;p [GeV]", 25, 0., 10.);
  meEndcapPIDp_ = ibook.book1D("EndcapPIDp", "PID track with MTD momentum spectrum, |eta| > 1.6;p [GeV]", 25, 0., 10.);

  meBarrelNoPIDtype_ = ibook.book1D("BarrelNoPIDtype", "Barrel PID failure category", 4, 0., 4.);
  meEndcapNoPIDtype_ = ibook.book1D("EndcapNoPIDtype", "Endcap PID failure category", 4, 0., 4.);

  meBarrelTruePiNoPID_ =
      ibook.book1D("BarrelTruePiNoPID", "True pi NoPID momentum spectrum, |eta| < 1.5;p [GeV]", 25, 0., 10.);
  meBarrelTrueKNoPID_ =
      ibook.book1D("BarrelTrueKNoPID", "True k NoPID momentum spectrum, |eta| < 1.5;p [GeV]", 25, 0., 10.);
  meBarrelTruePNoPID_ =
      ibook.book1D("BarrelTruePNoPID", "True p NoPID momentum spectrum, |eta| < 1.5;p [GeV]", 25, 0., 10.);
  meEndcapTruePiNoPID_ =
      ibook.book1D("EndcapTruePiNoPID", "True pi NoPID momentum spectrum, |eta| > 1.6;p [GeV]", 25, 0., 10.);
  meEndcapTrueKNoPID_ =
      ibook.book1D("EndcapTrueKNoPID", "True k NoPID momentum spectrum, |eta| > 1.6;p [GeV]", 25, 0., 10.);
  meEndcapTruePNoPID_ =
      ibook.book1D("EndcapTruePNoPID", "True p NoPID momentum spectrum, |eta| > 1.6;p [GeV]", 25, 0., 10.);

  meBarrelTruePiAsPi_ =
      ibook.book1D("BarrelTruePiAsPi", "True pi as pi momentum spectrum, |eta| < 1.5;p [GeV]", 25, 0., 10.);
  meBarrelTruePiAsK_ =
      ibook.book1D("BarrelTruePiAsK", "True pi as k momentum spectrum, |eta| < 1.5;p [GeV]", 25, 0., 10.);
  meBarrelTruePiAsP_ =
      ibook.book1D("BarrelTruePiAsP", "True pi as p momentum spectrum, |eta| < 1.5;p [GeV]", 25, 0., 10.);
  meEndcapTruePiAsPi_ =
      ibook.book1D("EndcapTruePiAsPi", "True pi as pi momentum spectrum, |eta| > 1.6;p [GeV]", 25, 0., 10.);
  meEndcapTruePiAsK_ =
      ibook.book1D("EndcapTruePiAsK", "True pi as k momentum spectrum, |eta| > 1.6;p [GeV]", 25, 0., 10.);
  meEndcapTruePiAsP_ =
      ibook.book1D("EndcapTruePiAsP", "True pi as p momentum spectrum, |eta| > 1.6;p [GeV]", 25, 0., 10.);

  meBarrelTrueKAsPi_ =
      ibook.book1D("BarrelTrueKAsPi", "True k as pi momentum spectrum, |eta| < 1.5;p [GeV]", 25, 0., 10.);
  meBarrelTrueKAsK_ = ibook.book1D("BarrelTrueKAsK", "True k as k momentum spectrum, |eta| < 1.5;p [GeV]", 25, 0., 10.);
  meBarrelTrueKAsP_ = ibook.book1D("BarrelTrueKAsP", "True k as p momentum spectrum, |eta| < 1.5;p [GeV]", 25, 0., 10.);
  meEndcapTrueKAsPi_ =
      ibook.book1D("EndcapTrueKAsPi", "True k as pi momentum spectrum, |eta| > 1.6;p [GeV]", 25, 0., 10.);
  meEndcapTrueKAsK_ = ibook.book1D("EndcapTrueKAsK", "True k as k momentum spectrum, |eta| > 1.6;p [GeV]", 25, 0., 10.);
  meEndcapTrueKAsP_ = ibook.book1D("EndcapTrueKAsP", "True k as p momentum spectrum, |eta| > 1.6;p [GeV]", 25, 0., 10.);

  meBarrelTruePAsPi_ =
      ibook.book1D("BarrelTruePAsPi", "True p as pi momentum spectrum, |eta| < 1.5;p [GeV]", 25, 0., 10.);
  meBarrelTruePAsK_ = ibook.book1D("BarrelTruePAsK", "True p as k momentum spectrum, |eta| < 1.5;p [GeV]", 25, 0., 10.);
  meBarrelTruePAsP_ = ibook.book1D("BarrelTruePAsP", "True p as p momentum spectrum, |eta| < 1.5;p [GeV]", 25, 0., 10.);
  meEndcapTruePAsPi_ =
      ibook.book1D("EndcapTruePAsPi", "True p as pi momentum spectrum, |eta| > 1.6;p [GeV]", 25, 0., 10.);
  meEndcapTruePAsK_ = ibook.book1D("EndcapTruePAsK", "True p as k momentum spectrum, |eta| > 1.6;p [GeV]", 25, 0., 10.);
  meEndcapTruePAsP_ = ibook.book1D("EndcapTruePAsP", "True p as p momentum spectrum, |eta| > 1.6;p [GeV]", 25, 0., 10.);
}

fillDescriptions()

void Primary4DVertexValidation::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 2653 of file Primary4DVertexValidation.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                           {
  edm::ParameterSetDescription desc;

  desc.add<std::string>("folder", "MTD/Vertices");
  desc.add<edm::InputTag>("TPtoRecoTrackAssoc", edm::InputTag("trackingParticleRecoTrackAsssociation"));
  desc.add<edm::InputTag>("mtdTracks", edm::InputTag("trackExtenderWithMTD"));
  desc.add<edm::InputTag>("SimTag", edm::InputTag("mix", "MergedTrackTruth"));
  desc.add<edm::InputTag>("offlineBS", edm::InputTag("offlineBeamSpot"));
  desc.add<edm::InputTag>("offline4DPV", edm::InputTag("offlinePrimaryVertices4D"));
  desc.add<edm::InputTag>("trackAssocSrc", edm::InputTag("trackExtenderWithMTD:generalTrackassoc"))
      ->setComment("Association between General and MTD Extended tracks");
  desc.add<edm::InputTag>("pathLengthSrc", edm::InputTag("trackExtenderWithMTD:generalTrackPathLength"));
  desc.add<edm::InputTag>("momentumSrc", edm::InputTag("trackExtenderWithMTD:generalTrackp"));
  desc.add<edm::InputTag>("tmtd", edm::InputTag("trackExtenderWithMTD:generalTracktmtd"));
  desc.add<edm::InputTag>("timeSrc", edm::InputTag("trackExtenderWithMTD:generalTracktmtd"));
  desc.add<edm::InputTag>("sigmaSrc", edm::InputTag("trackExtenderWithMTD:generalTracksigmatmtd"));
  desc.add<edm::InputTag>("t0PID", edm::InputTag("tofPID:t0"));
  desc.add<edm::InputTag>("t0SafePID", edm::InputTag("tofPID:t0safe"));
  desc.add<edm::InputTag>("sigmat0SafePID", edm::InputTag("tofPID:sigmat0safe"));
  desc.add<edm::InputTag>("trackMVAQual", edm::InputTag("mtdTrackQualityMVA:mtdQualMVA"));
  desc.add<edm::InputTag>("tofPi", edm::InputTag("trackExtenderWithMTD:generalTrackTofPi"));
  desc.add<edm::InputTag>("tofK", edm::InputTag("trackExtenderWithMTD:generalTrackTofK"));
  desc.add<edm::InputTag>("tofP", edm::InputTag("trackExtenderWithMTD:generalTrackTofP"));
  desc.add<edm::InputTag>("probPi", edm::InputTag("tofPID:probPi"));
  desc.add<edm::InputTag>("probK", edm::InputTag("tofPID:probK"));
  desc.add<edm::InputTag>("probP", edm::InputTag("tofPID:probP"));
  desc.add<bool>("useOnlyChargedTracks", true);
  desc.addUntracked<bool>("debug", false);
  desc.addUntracked<bool>("optionalPlots", false);
  desc.add<bool>("use3dNoTime", false);
  desc.add<double>("trackweightTh", 0.5);
  desc.add<double>("mvaTh", 0.01);
  desc.add<double>("minProbHeavy", 0.75);

  //lineDensity parameters have been obtained by fitting the distribution of the z position of the vertices,
  //using a 200k single mu ptGun sample (gaussian fit)
  std::vector<double> lDP;
  lDP.push_back(1.87);
  lDP.push_back(0.);
  lDP.push_back(42.5);
  desc.add<std::vector<double>>("lineDensityPar", lDP);
  descriptions.add("vertices4DValid", desc);
}

getMatchedTP()

std::pair< const edm::Ref< std::vector< TrackingParticle > > *, int > Primary4DVertexValidation::getMatchedTP ( const reco::TrackBaseRef &  recoTrack, const TrackingVertexRef &  vsim  )

private

Definition at line 1367 of file Primary4DVertexValidation.cc.

References edm::AssociationMap< Tag >::end(), edm::AssociationMap< Tag >::find(), newFWLiteAna::found, r2s_, and cmsswSequenceInfo::tp.

Referenced by analyze(), and matchReco2Sim().

                                                                      {
  auto found = r2s_->find(recoTrack);

  // reco track not matched to any TP (fake tracks)
  if (found == r2s_->end())
    return std::make_pair(nullptr, -1);

  // matched TP equal to any TP of a given sim vertex
  for (const auto& tp : found->val) {
    if (std::find_if(vsim->daughterTracks_begin(), vsim->daughterTracks_end(), [&](const TrackingParticleRef& vtp) {
          return tp.first == vtp;
        }) != vsim->daughterTracks_end())
      return std::make_pair(&tp.first, 0);
    // matched TP not associated to any daughter track of a given sim vertex but having the same eventID (track from secondary vtx)
    else if (tp.first->eventId().bunchCrossing() == vsim->eventId().bunchCrossing() &&
             tp.first->eventId().event() == vsim->eventId().event()) {
      return std::make_pair(&tp.first, 1);
    }
    // matched TP not associated to any sim vertex of a given simulated event (PU track)
    else {
      return std::make_pair(&tp.first, 2);
    }
  }

  // reco track not matched to any TP from vertex
  return std::make_pair(nullptr, -1);
}

getRecoPVs()

std::vector< Primary4DVertexValidation::recoPrimaryVertex > Primary4DVertexValidation::getRecoPVs ( const edm::Handle< edm::View< reco::Vertex >> &  tVC )

private

Definition at line 1653 of file Primary4DVertexValidation.cc.

References funct::abs(), HLT_2024v10_cff::distance, edm::AssociationMap< Tag >::end(), edm::AssociationMap< Tag >::find(), muonTagProbeFilters_cff::matched, SiStripPI::min, Primary4DVertexValidation::recoPrimaryVertex::nRecoTrk, Primary4DVertexValidation::recoPrimaryVertex::num_matched_sim_tracks, Primary4DVertexValidation::recoPrimaryVertex::pt, Primary4DVertexValidation::recoPrimaryVertex::ptsq, alignCSCRings::r, r2s_, mathSSE::sqrt(), pfDeepBoostedJetPreprocessParams_cfi::sv, and findQualityFiles::v.

Referenced by analyze().

                                               {
  std::vector<Primary4DVertexValidation::recoPrimaryVertex> recopv;
  int r = -1;
  for (auto v = tVC->begin(); v != tVC->end(); ++v) {
    r++;
    // Skip junk vertices
    if (std::abs(v->z()) > 1000)
      continue;
    if (v->isFake() || !v->isValid())
      continue;

    recoPrimaryVertex sv(v->position().x(), v->position().y(), v->position().z());
    sv.recVtx = &(*v);
    sv.recVtxRef = reco::VertexBaseRef(tVC, std::distance(tVC->begin(), v));

    sv.OriginalIndex = r;
    sv.ndof = v->ndof();
    // this is a new vertex, add it to the list of reco-vertices
    recopv.push_back(sv);
    Primary4DVertexValidation::recoPrimaryVertex* vp = &recopv.back();

    // Loop over daughter track(s)
    for (auto iTrack = v->tracks_begin(); iTrack != v->tracks_end(); ++iTrack) {
      auto momentum = (*(*iTrack)).innerMomentum();
      if (momentum.mag2() == 0)
        momentum = (*(*iTrack)).momentum();
      vp->pt += std::sqrt(momentum.perp2());
      vp->ptsq += (momentum.perp2());
      vp->nRecoTrk++;

      auto matched = r2s_->find(*iTrack);
      if (matched != r2s_->end()) {
        vp->num_matched_sim_tracks++;
      }

    }  // End of for loop on daughters reconstructed tracks
  }    // End of for loop on tracking vertices

  // In case of no reco vertices, break here
  if (recopv.empty())
    return recopv;

  // Now compute the closest distance in z between all reconstructed vertex
  // first initialize
  auto prev_z = recopv.back().z;
  for (recoPrimaryVertex& vreco : recopv) {
    vreco.closest_vertex_distance_z = std::abs(vreco.z - prev_z);
    prev_z = vreco.z;
  }
  for (std::vector<recoPrimaryVertex>::iterator vreco = recopv.begin(); vreco != recopv.end(); vreco++) {
    std::vector<recoPrimaryVertex>::iterator vreco2 = vreco;
    vreco2++;
    for (; vreco2 != recopv.end(); vreco2++) {
      double distance = std::abs(vreco->z - vreco2->z);
      // need both to be complete
      vreco->closest_vertex_distance_z = std::min(vreco->closest_vertex_distance_z, distance);
      vreco2->closest_vertex_distance_z = std::min(vreco2->closest_vertex_distance_z, distance);
    }
  }
  return recopv;
}

getSimPVs()

std::vector< Primary4DVertexValidation::simPrimaryVertex > Primary4DVertexValidation::getSimPVs ( const edm::Handle< TrackingVertexCollection > &  tVC )

private

Definition at line 1550 of file Primary4DVertexValidation.cc.

References funct::abs(), cms::cuda::assert(), Primary4DVertexValidation::simPrimaryVertex::average_match_quality, debug_, HLT_2024v10_cff::distance, MillePedeFileConverter_cfg::e, edm::AssociationMap< Tag >::end(), edm::AssociationMap< Tag >::find(), SiStripPI::min, Primary4DVertexValidation::simPrimaryVertex::nGenTrk, Primary4DVertexValidation::simPrimaryVertex::num_matched_reco_tracks, Primary4DVertexValidation::simPrimaryVertex::pt, Primary4DVertexValidation::simPrimaryVertex::ptot, Primary4DVertexValidation::simPrimaryVertex::ptsq, alignCSCRings::s, s2r_, pfDeepBoostedJetPreprocessParams_cfi::sv, use_only_charged_tracks_, findQualityFiles::v, and V0Monitor_cff::v0.

Referenced by analyze().

                                                    {
  std::vector<Primary4DVertexValidation::simPrimaryVertex> simpv;
  int current_event = -1;
  int s = -1;
  for (TrackingVertexCollection::const_iterator v = tVC->begin(); v != tVC->end(); ++v) {
    // We keep only the first vertex from all the events at BX=0.
    if (v->eventId().bunchCrossing() != 0)
      continue;
    if (v->eventId().event() != current_event) {
      current_event = v->eventId().event();
    } else {
      continue;
    }
    s++;
    if (std::abs(v->position().z()) > 1000)
      continue;  // skip junk vertices

    // could be a new vertex, check  all primaries found so far to avoid multiple entries
    simPrimaryVertex sv(v->position().x(), v->position().y(), v->position().z(), v->position().t());
    sv.eventId = v->eventId();
    sv.sim_vertex = TrackingVertexRef(tVC, std::distance(tVC->begin(), v));
    sv.OriginalIndex = s;

    for (TrackingParticleRefVector::iterator iTrack = v->daughterTracks_begin(); iTrack != v->daughterTracks_end();
         ++iTrack) {
      assert((**iTrack).eventId().bunchCrossing() == 0);
    }
    simPrimaryVertex* vp = nullptr;  // will become non-NULL if a vertex is found and then point to it
    for (std::vector<simPrimaryVertex>::iterator v0 = simpv.begin(); v0 != simpv.end(); v0++) {
      if ((sv.eventId == v0->eventId) && (std::abs(sv.x - v0->x) < 1e-5) && (std::abs(sv.y - v0->y) < 1e-5) &&
          (std::abs(sv.z - v0->z) < 1e-5)) {
        vp = &(*v0);
        break;
      }
    }
    if (!vp) {
      // this is a new vertex, add it to the list of sim-vertices
      simpv.push_back(sv);
      vp = &simpv.back();
    }

    // Loop over daughter track(s) as Tracking Particles
    for (TrackingVertex::tp_iterator iTP = v->daughterTracks_begin(); iTP != v->daughterTracks_end(); ++iTP) {
      auto momentum = (*(*iTP)).momentum();
      const reco::Track* matched_best_reco_track = nullptr;
      double match_quality = -1;
      if (use_only_charged_tracks_ && (**iTP).charge() == 0)
        continue;
      if (s2r_->find(*iTP) != s2r_->end()) {
        matched_best_reco_track = (*s2r_)[*iTP][0].first.get();
        match_quality = (*s2r_)[*iTP][0].second;
      }

      vp->ptot.setPx(vp->ptot.x() + momentum.x());
      vp->ptot.setPy(vp->ptot.y() + momentum.y());
      vp->ptot.setPz(vp->ptot.z() + momentum.z());
      vp->ptot.setE(vp->ptot.e() + (**iTP).energy());
      vp->pt += (**iTP).pt();
      vp->ptsq += ((**iTP).pt() * (**iTP).pt());
      vp->nGenTrk++;

      if (matched_best_reco_track) {
        vp->num_matched_reco_tracks++;
        vp->average_match_quality += match_quality;
      }
    }  // End of for loop on daughters sim-particles
    if (vp->num_matched_reco_tracks)
      vp->average_match_quality /= static_cast<float>(vp->num_matched_reco_tracks);
    if (debug_) {
      edm::LogPrint("Primary4DVertexValidation")
          << "average number of associated tracks: " << vp->num_matched_reco_tracks / static_cast<float>(vp->nGenTrk)
          << " with average quality: " << vp->average_match_quality;
    }
  }  // End of for loop on tracking vertices

  // In case of no simulated vertices, break here
  if (simpv.empty())
    return simpv;

  // Now compute the closest distance in z between all simulated vertex
  // first initialize
  auto prev_z = simpv.back().z;
  for (simPrimaryVertex& vsim : simpv) {
    vsim.closest_vertex_distance_z = std::abs(vsim.z - prev_z);
    prev_z = vsim.z;
  }
  // then calculate
  for (std::vector<simPrimaryVertex>::iterator vsim = simpv.begin(); vsim != simpv.end(); vsim++) {
    std::vector<simPrimaryVertex>::iterator vsim2 = vsim;
    vsim2++;
    for (; vsim2 != simpv.end(); vsim2++) {
      double distance = std::abs(vsim->z - vsim2->z);
      // need both to be complete
      vsim->closest_vertex_distance_z = std::min(vsim->closest_vertex_distance_z, distance);
      vsim2->closest_vertex_distance_z = std::min(vsim2->closest_vertex_distance_z, distance);
    }
  }
  return simpv;
}

getWosWnt()

void Primary4DVertexValidation::getWosWnt ( const reco::Vertex &  recoVtx, const reco::TrackBaseRef &  recoTrk, const edm::ValueMap< float > &  sigmat0, const edm::Handle< reco::BeamSpot > &  BS, double &  wos, double &  wnt  )

private

Definition at line 1527 of file Primary4DVertexValidation.cc.

References c_, funct::cos(), reco::TrackBase::dzError(), SiStripPI::min, reco::TrackBase::phi(), funct::pow(), reco::TrackBase::pt(), beamSpotPI::sigmaZ, funct::sin(), funct::tan(), reco::TrackBase::theta(), and reco::Vertex::trackWeight().

Referenced by analyze(), and matchReco2Sim().

                                                       {
  double dz2_beam = pow((*BS).BeamWidthX() * cos(recoTrk->phi()) / tan(recoTrk->theta()), 2) +
                    pow((*BS).BeamWidthY() * sin(recoTrk->phi()) / tan(recoTrk->theta()), 2);
  double dz2 =
      pow(recoTrk->dzError(), 2) + dz2_beam + pow(0.0020, 2);  // added 20 um, some tracks have crazy small resolutions
  wos = recoVtx.trackWeight(recoTrk) / dz2;
  wnt = recoVtx.trackWeight(recoTrk) * std::min(recoTrk->pt(), 1.0);

  if (sigmat0[recoTrk] > 0) {
    double sigmaZ = (*BS).sigmaZ();
    double sigmaT = sigmaZ / c_;  // c in cm/ns
    wos = wos / erf(sigmat0[recoTrk] / sigmaT);
  }
}

isParticle()

void Primary4DVertexValidation::isParticle ( const reco::TrackBaseRef &  recoTrack, const edm::ValueMap< float > &  probPi, const edm::ValueMap< float > &  probK, const edm::ValueMap< float > &  probP, unsigned int &  no_PIDtype, bool &  no_PID, bool &  is_Pi, bool &  is_K, bool &  is_P  )

private

Definition at line 1500 of file Primary4DVertexValidation.cc.

References edm::isNotFinite(), and minProbHeavy_.

Referenced by analyze().

                                                       {
  no_PIDtype = 0;
  no_PID = false;
  is_Pi = false;
  is_K = false;
  is_P = false;
  if (probPi[recoTrack] == -1) {
    no_PIDtype = 1;
  } else if (edm::isNotFinite(probPi[recoTrack])) {
    no_PIDtype = 2;
  } else if (probPi[recoTrack] == 1 && probK[recoTrack] == 0 && probP[recoTrack] == 0) {
    no_PIDtype = 3;
  }
  no_PID = no_PIDtype > 0;
  is_Pi = !no_PID && 1. - probPi[recoTrack] < minProbHeavy_;
  is_K = !no_PID && !is_Pi && probK[recoTrack] > probP[recoTrack];
  is_P = !no_PID && !is_Pi && !is_K;
}

matchReco2Sim()

void Primary4DVertexValidation::matchReco2Sim ( std::vector< recoPrimaryVertex > &  recopv, std::vector< simPrimaryVertex > &  simpv, const edm::ValueMap< float > &  sigmat0, const edm::ValueMap< float > &  MVA, const edm::Handle< reco::BeamSpot > &  BS  )

private

Definition at line 1717 of file Primary4DVertexValidation.cc.

References funct::abs(), debug_, dqmdumpme::first, getMatchedTP(), getWosWnt(), mps_fire::i, caHitNtupletGeneratorKernels::iev, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::iv, dqmdumpme::k, maxRank_, maxTry_, mvaTh_, NOT_MATCHED, edm::second(), pfDeepBoostedJetPreprocessParams_cfi::sv, trackweightTh_, bphysicsOniaDQM_cfi::vertex, and zWosMatchMax_.

Referenced by analyze().

                                                                                 {
  for (auto vv : simpv) {
    vv.wnt.clear();
    vv.wos.clear();
  }
  for (auto rv : recopv) {
    rv.wnt.clear();
    rv.wos.clear();
  }

  for (unsigned int iv = 0; iv < recopv.size(); iv++) {
    const reco::Vertex* vertex = recopv.at(iv).recVtx;

    for (unsigned int iev = 0; iev < simpv.size(); iev++) {
      double wnt = 0;
      double wos = 0;
      double evwnt = 0;
      double evwos = 0;
      double evnt = 0;

      for (auto iTrack = vertex->tracks_begin(); iTrack != vertex->tracks_end(); ++iTrack) {
        if (vertex->trackWeight(*iTrack) < trackweightTh_)
          continue;
        if (MVA[(*iTrack)] < mvaTh_)
          continue;

        auto tp_info = getMatchedTP(*iTrack, simpv.at(iev).sim_vertex).first;
        int matchCategory = getMatchedTP(*iTrack, simpv.at(iev).sim_vertex).second;
        // matched TP equal to any TP of a given sim vertex
        if (tp_info != nullptr && matchCategory == 0) {
          getWosWnt(*vertex, *iTrack, sigmat0, BS, wos, wnt);
          simpv.at(iev).addTrack(iv, wos, wnt);
          recopv.at(iv).addTrack(iev, wos, wnt);
          evwos += wos;
          evwnt += wnt;
          evnt++;
        }
      }  // RecoTracks loop

      // require 2 tracks for a wos-match
      if ((evwos > 0) && (evwos > recopv.at(iv).maxwos) && (evnt > 1)) {
        recopv.at(iv).wosmatch = iev;
        recopv.at(iv).maxwos = evwos;
        recopv.at(iv).maxwosnt = evnt;

        simpv.at(iev).wos_dominated_recv.push_back(iv);
        simpv.at(iev).nwosmatch++;
      }

      // weighted track counting match, require at least one track
      if ((evnt > 0) && (evwnt > recopv.at(iv).maxwnt)) {
        recopv.at(iv).wntmatch = iev;
        recopv.at(iv).maxwnt = evwnt;
      }
    }  // TrackingVertex loop

  }  // RecoPrimaryVertex

  // after filling infos, goes for the sim-reco match
  for (auto& vrec : recopv) {
    vrec.sim = NOT_MATCHED;
    vrec.matchQuality = 0;
  }
  unsigned int iev = 0;
  for (auto& vv : simpv) {
    if (debug_) {
      edm::LogPrint("Primary4DVertexValidation") << "iev: " << iev;
      edm::LogPrint("Primary4DVertexValidation") << "wos_dominated_recv.size: " << vv.wos_dominated_recv.size();
    }
    for (unsigned int i = 0; i < vv.wos_dominated_recv.size(); i++) {
      auto recov = vv.wos_dominated_recv.at(i);
      if (debug_) {
        edm::LogPrint("Primary4DVertexValidation")
            << "index of reco vertex: " << recov << " that has a wos: " << vv.wos.at(recov) << " at position " << i;
      }
    }
    vv.rec = NOT_MATCHED;
    vv.matchQuality = 0;
    iev++;
  }
  // this tries a one-to-one match, taking simPV with highest wos if there are > 1 simPV candidates
  for (unsigned int rank = 1; rank < maxRank_; rank++) {
    for (unsigned int iev = 0; iev < simpv.size(); iev++) {  //loop on SimPV
      if (simpv.at(iev).rec != NOT_MATCHED)
        continue;
      if (simpv.at(iev).nwosmatch == 0)
        continue;
      if (simpv.at(iev).nwosmatch > rank)
        continue;
      unsigned int iv = NOT_MATCHED;
      for (unsigned int k = 0; k < simpv.at(iev).wos_dominated_recv.size(); k++) {
        unsigned int rec = simpv.at(iev).wos_dominated_recv.at(k);
        auto vrec = recopv.at(rec);
        if (vrec.sim != NOT_MATCHED)
          continue;  // already matched
        if (std::abs(simpv.at(iev).z - vrec.z) > zWosMatchMax_)
          continue;  // insanely far away
        if ((iv == NOT_MATCHED) || simpv.at(iev).wos.at(rec) > simpv.at(iev).wos.at(iv)) {
          iv = rec;
        }
      }
      if (iv !=
          NOT_MATCHED) {  // if the rec vertex has already been associated is possible that iv remains NOT_MATCHED at this point
        recopv.at(iv).sim = iev;
        simpv.at(iev).rec = iv;
        recopv.at(iv).matchQuality = rank;
        simpv.at(iev).matchQuality = rank;
      }
    }
  }
  //give vertices a chance that have a lot of overlap, but are still recognizably
  //caused by a specific simvertex (without being classified as dominating)
  //like a small peak sitting on the flank of a larger nearby peak
  unsigned int ntry = 0;
  while (ntry++ < maxTry_) {
    unsigned nmatch = 0;
    for (unsigned int iev = 0; iev < simpv.size(); iev++) {
      if ((simpv.at(iev).rec != NOT_MATCHED) || (simpv.at(iev).wos.empty()))
        continue;
      // find a rec vertex for the NOT_MATCHED sim vertex
      unsigned int rec = NOT_MATCHED;
      for (auto rv : simpv.at(iev).wos) {
        if ((rec == NOT_MATCHED) || (rv.second > simpv.at(iev).wos.at(rec))) {
          rec = rv.first;
        }
      }

      if (rec == NOT_MATCHED) {  // try with wnt match
        for (auto rv : simpv.at(iev).wnt) {
          if ((rec == NOT_MATCHED) || (rv.second > simpv.at(iev).wnt.at(rec))) {
            rec = rv.first;
          }
        }
      }

      if (rec == NOT_MATCHED)
        continue;
      if (recopv.at(rec).sim != NOT_MATCHED)
        continue;  // already gone

      // check if the recvertex can be  matched
      unsigned int rec2sim = NOT_MATCHED;
      for (auto sv : recopv.at(rec).wos) {
        if (simpv.at(sv.first).rec != NOT_MATCHED)
          continue;  // already used
        if ((rec2sim == NOT_MATCHED) || (sv.second > recopv.at(rec).wos.at(rec2sim))) {
          rec2sim = sv.first;
        }
      }
      if (iev == rec2sim) {
        // do the match and assign lowest quality (i.e. max rank)
        recopv.at(rec).sim = iev;
        recopv.at(rec).matchQuality = maxRank_;
        simpv.at(iev).rec = rec;
        simpv.at(iev).matchQuality = maxRank_;
        nmatch++;
      }
    }  // sim loop
    if (nmatch == 0) {
      break;
    }
  }  // ntry
}

matchRecoTrack2SimSignal()

bool Primary4DVertexValidation::matchRecoTrack2SimSignal ( const reco::TrackBaseRef &  recoTrack )

private

Definition at line 1350 of file Primary4DVertexValidation.cc.

References edm::AssociationMap< Tag >::end(), edm::AssociationMap< Tag >::find(), newFWLiteAna::found, r2s_, and cmsswSequenceInfo::tp.

                                                                                          {
  auto found = r2s_->find(recoTrack);

  // reco track not matched to any TP
  if (found == r2s_->end())
    return false;

  for (const auto& tp : found->val) {
    if (tp.first->eventId().bunchCrossing() == 0 && tp.first->eventId().event() == 0)
      return true;
  }

  // reco track not matched to any TP from signal vertex
  return false;
}

mvaRecSel()

const bool Primary4DVertexValidation::mvaRecSel ( const reco::TrackBase &  trk, const reco::Vertex &  vtx, const double &  t0, const double &  st0  )

private

Definition at line 2758 of file Primary4DVertexValidation.cc.

References funct::abs(), deltaZcut_, reco::TrackBase::eta(), etacutREC_, match(), reco::TrackBase::pt(), pTcut_, FrontierCondition_GT_autoExpress_cfi::t0, L1BJetProducer_cff::vtx, and reco::TrackBase::vz().

Referenced by analyze().

                                                                   {
  bool match = false;
  match = trk.pt() > pTcut_ && std::abs(trk.eta()) < etacutREC_ && std::abs(trk.vz() - vtx.z()) <= deltaZcut_;
  if (st0 > 0.) {
    match = match && std::abs(t0 - vtx.t()) < 3. * st0;
  }
  return match;
}

mvaTPSel()

const bool Primary4DVertexValidation::mvaTPSel ( const TrackingParticle &  tp )

private

Definition at line 2749 of file Primary4DVertexValidation.cc.

References funct::abs(), etacutGEN_, match(), pTcut_, and cmsswSequenceInfo::tp.

Referenced by analyze().

                                                                         {
  bool match = false;
  if (tp.status() != 1) {
    return match;
  }
  match = tp.charge() != 0 && tp.pt() > pTcut_ && std::abs(tp.eta()) < etacutGEN_;
  return match;
}

observablesFromJets()

void Primary4DVertexValidation::observablesFromJets ( const std::vector< reco::Track > &  reco_Tracks, const std::vector< double > &  mass_Tracks, const std::vector< int > &  category_Tracks, const std::string &  skip_Tracks, unsigned int &  n_Jets, double &  sum_EtJets, double &  sum_Pt2Jets, double &  met_Pt, double &  sum_PzJets  )

private

Definition at line 1427 of file Primary4DVertexValidation.cc.

References mps_fire::i, edm::isNotFinite(), PDWG_EXODelayedJetMET_cff::jets, EgHLTOffHistBins_cfi::mass, BTaggingMonitor_cfi::met, l1tEGammaCrystalsEmulatorProducer_cfi::scale, and mathSSE::sqrt().

Referenced by analyze().

                                                                        {
  double sum_PtJets = 0;
  n_Jets = 0;
  sum_EtJets = 0;
  sum_Pt2Jets = 0;
  met_Pt = 0;
  sum_PzJets = 0;
  auto met = LorentzVector(0, 0, 0, 0);
  std::vector<fastjet::PseudoJet> fjInputs_;
  fjInputs_.clear();
  size_t countScale0 = 0;
  for (size_t i = 0; i < reco_Tracks.size(); i++) {
    const auto recotr = reco_Tracks[i];
    const auto mass = mass_Tracks[i];
    float scale = 1.;
    if (recotr.charge() == 0) {
      continue;
    }
    // skip PU tracks in jet definition if skip_PU is required
    if (skip_Tracks == "skip_PU" && category_Tracks[i] == 2) {
      continue;
    }
    // skip fake tracks in jet definition if skip_Fake is required
    if (skip_Tracks == "skip_Fake" && category_Tracks[i] == -1) {
      continue;
    }
    if (recotr.pt() != 0) {
      scale = (recotr.pt() - recotr.ptError()) / recotr.pt();
    }
    if (edm::isNotFinite(scale)) {
      edm::LogWarning("Primary4DVertexValidation") << "Scaling is NAN ignoring this recotrack" << std::endl;
      scale = 0;
    }
    if (scale < 0) {
      scale = 0;
      countScale0++;
    }
    if (scale != 0) {
      fjInputs_.push_back(fastjet::PseudoJet(recotr.px() * scale,
                                             recotr.py() * scale,
                                             recotr.pz() * scale,
                                             std::sqrt(recotr.p() * recotr.p() + mass * mass) * scale));
    }
  }
  fastjet::ClusterSequence sequence(fjInputs_, fastjet::JetDefinition(fastjet::antikt_algorithm, 0.4));
  auto jets = fastjet::sorted_by_pt(sequence.inclusive_jets(0));
  for (const auto& pj : jets) {
    auto p4 = LorentzVector(pj.px(), pj.py(), pj.pz(), pj.e());
    sum_EtJets += std::sqrt(p4.e() * p4.e() - p4.P() * p4.P() + p4.pt() * p4.pt());
    sum_PtJets += p4.pt();
    sum_Pt2Jets += (p4.pt() * p4.pt() * 0.8 * 0.8);
    met += p4;
    sum_PzJets += p4.pz();
    n_Jets++;
  }
  met_Pt = met.pt();
  double metAbove = met_Pt - 2 * std::sqrt(sum_PtJets);
  if (metAbove > 0) {
    sum_Pt2Jets += (metAbove * metAbove);
  }
  if (countScale0 == reco_Tracks.size()) {
    sum_Pt2Jets = countScale0 * 0.01;  //leave some epsilon value to sort vertices with unknown pt
  }
}

printMatchedRecoTrackInfo()

void Primary4DVertexValidation::printMatchedRecoTrackInfo ( const reco::Vertex &  vtx, const reco::TrackBaseRef &  trk, const TrackingParticleRef &  tp, const unsigned int &  categ  )

private

Definition at line 2697 of file Primary4DVertexValidation.cc.

References reco::TrackBase::eta(), reco::TrackBase::pt(), simUnit_, AlCaHLTBitMon_QueryRunRegistry::string, cmsswSequenceInfo::tp, and L1BJetProducer_cff::vtx.

Referenced by analyze().

                                                                                     {
  std::string strTrk;
  switch (categ) {
    case 0:
      strTrk = "Reco_Track:";
      break;
    case 1:
      strTrk = "SecRecoTrk:";
      break;
    case 2:
      strTrk = "PU_RecoTrk:";
      break;
  }
  edm::LogPrint("Primary4DVertexValidation")
      << strTrk << " w =" << std::setw(6) << std::setprecision(2) << vtx.trackWeight(trk) << " pt =" << std::setw(6)
      << std::setprecision(2) << trk->pt() << " eta =" << std::setw(6) << std::setprecision(2) << trk->eta()
      << "  MatchedTP: Pt =" << std::setw(6) << std::setprecision(2) << tp->pt() << " eta =" << std::setw(6)
      << std::setprecision(2) << tp->eta() << "  Parent vtx: z =" << std::setw(8) << std::setprecision(4)
      << tp->parentVertex()->position().z() << " t =" << std::setw(8) << std::setprecision(4)
      << tp->parentVertex()->position().t() * simUnit_ << " BX =" << tp->parentVertex()->eventId().bunchCrossing()
      << " ev =" << tp->parentVertex()->eventId().event() << std::endl;
}

printSimVtxRecoVtxInfo()

void Primary4DVertexValidation::printSimVtxRecoVtxInfo ( const struct Primary4DVertexValidation::simPrimaryVertex &  simpVtx, const struct Primary4DVertexValidation::recoPrimaryVertex &  recopVtx  )

private

Definition at line 2723 of file Primary4DVertexValidation.cc.

References Primary4DVertexValidation::simPrimaryVertex::nGenTrk, Primary4DVertexValidation::recoPrimaryVertex::nRecoTrk, Primary4DVertexValidation::simPrimaryVertex::num_matched_reco_tracks, Primary4DVertexValidation::recoPrimaryVertex::num_matched_sim_tracks, Primary4DVertexValidation::simPrimaryVertex::pt, Primary4DVertexValidation::recoPrimaryVertex::pt, Primary4DVertexValidation::simPrimaryVertex::ptsq, Primary4DVertexValidation::recoPrimaryVertex::ptsq, Primary4DVertexValidation::simPrimaryVertex::sim_vertex, simUnit_, Primary4DVertexValidation::recoPrimaryVertex::sumwnt, Primary4DVertexValidation::recoPrimaryVertex::sumwos, Primary4DVertexValidation::simPrimaryVertex::t, use_only_charged_tracks_, Primary4DVertexValidation::simPrimaryVertex::x, Primary4DVertexValidation::recoPrimaryVertex::x, Primary4DVertexValidation::simPrimaryVertex::y, Primary4DVertexValidation::recoPrimaryVertex::y, Primary4DVertexValidation::simPrimaryVertex::z, and Primary4DVertexValidation::recoPrimaryVertex::z.

Referenced by analyze().

                                                                       {
  edm::LogPrint("Primary4DVertexValidation")
      << "Sim vtx (x,y,z,t) = (" << std::setprecision(4) << simpVtx.x << "," << std::setprecision(4) << simpVtx.y << ","
      << std::setprecision(4) << simpVtx.z << "," << std::setprecision(4) << simpVtx.t * simUnit_ << ")";
  edm::LogPrint("Primary4DVertexValidation")
      << "Sim vtx: pt = " << std::setprecision(4) << simpVtx.pt << " ptsq = " << std::setprecision(6) << simpVtx.ptsq
      << " nGenTrk = " << simpVtx.nGenTrk << " nmatch recotrks = " << simpVtx.num_matched_reco_tracks;
  edm::LogPrint("Primary4DVertexValidation")
      << "Reco vtx (x,y,z) = (" << std::setprecision(4) << recopVtx.x << "," << std::setprecision(4) << recopVtx.y
      << "," << std::setprecision(4) << recopVtx.z << ")";
  edm::LogPrint("Primary4DVertexValidation")
      << "Reco vtx: pt = " << std::setprecision(4) << recopVtx.pt << " ptsq = " << std::setprecision(6) << recopVtx.ptsq
      << " nrecotrks = " << recopVtx.nRecoTrk << " nmatch simtrks = " << recopVtx.num_matched_sim_tracks;
  edm::LogPrint("Primary4DVertexValidation") << "wnt " << recopVtx.sumwnt << " wos = " << recopVtx.sumwos;
  for (auto iTP = simpVtx.sim_vertex->daughterTracks_begin(); iTP != simpVtx.sim_vertex->daughterTracks_end(); ++iTP) {
    if (use_only_charged_tracks_ && (**iTP).charge() == 0) {
      continue;
    }
    edm::LogPrint("Primary4DVertexValidation")
        << "Daughter track of sim vertex: pt =" << std::setw(6) << std::setprecision(2) << (*iTP)->pt()
        << "  eta =" << std::setw(6) << std::setprecision(2) << (*iTP)->eta();
  }
}

select()

bool Primary4DVertexValidation::select ( const reco::Vertex &  v, int  level = 0  )

private

Definition at line 1408 of file Primary4DVertexValidation.cc.

References personalPlayback::level, selNdof_, and findQualityFiles::v.

Referenced by analyze().

                                                                     {
  /* level
   0  !isFake  && ndof>4  (default)
   1  !isFake  && ndof>4 && prob > 0.01
   2  !isFake  && ndof>4 && prob > 0.01 && ptmax2 > 0.4
   */
  if (v.isFake())
    return false;
  if ((level == 0) && (v.ndof() > selNdof_))
    return true;
  /*if ((level == 1) && (v.ndof() > selNdof_) && (vertex_pxy(v) > 0.01))
    return true;
  if ((level == 2) && (v.ndof() > selNdof_) && (vertex_pxy(v) > 0.01) && (vertex_ptmax2(v) > 0.4))
    return true;
  if ((level == 3) && (v.ndof() > selNdof_) && (vertex_ptmax2(v) < 0.4))
    return true;*/
  return false;
}

timeFromTrueMass()

double Primary4DVertexValidation::timeFromTrueMass ( double  mass, double  pathlength, double  momentum, double  time  )

private

Definition at line 1396 of file Primary4DVertexValidation.cc.

References c_, EgHLTOffHistBins_cfi::mass, mathSSE::sqrt(), hcalRecHitTable_cff::time, and findQualityFiles::v.

Referenced by analyze().

                                                                                                               {
  if (time > 0 && pathlength > 0 && mass > 0) {
    double gammasq = 1. + momentum * momentum / (mass * mass);
    double v = c_ * std::sqrt(1. - 1. / gammasq);  // cm / ns
    double t_est = time - (pathlength / v);

    return t_est;
  } else {
    return -1;
  }
}

Member Data Documentation

c_

constexpr double Primary4DVertexValidation::c_ = 2.99792458e1

staticprivate

Definition at line 260 of file Primary4DVertexValidation.cc.

Referenced by getWosWnt(), and timeFromTrueMass().

c_cm_ns

constexpr float Primary4DVertexValidation::c_cm_ns = geant_units::operators::convertMmToCm(CLHEP::c_light)

staticprivate

Definition at line 282 of file Primary4DVertexValidation.cc.

debug_

bool Primary4DVertexValidation::debug_

private

Definition at line 285 of file Primary4DVertexValidation.cc.

Referenced by analyze(), getSimPVs(), and matchReco2Sim().

deltaZcut_

constexpr double Primary4DVertexValidation::deltaZcut_ = 0.1

staticprivate

Definition at line 270 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and mvaRecSel().

etacutGEN_

constexpr double Primary4DVertexValidation::etacutGEN_ = 4.

staticprivate

Definition at line 267 of file Primary4DVertexValidation.cc.

Referenced by mvaTPSel().

etacutREC_

constexpr double Primary4DVertexValidation::etacutREC_ = 3.

staticprivate

Definition at line 268 of file Primary4DVertexValidation.cc.

Referenced by mvaRecSel().

folder_

const std::string Primary4DVertexValidation::folder_

private

Definition at line 257 of file Primary4DVertexValidation.cc.

Referenced by bookHistograms().

lineDensityPar_

const std::vector<double> Primary4DVertexValidation::lineDensityPar_

private

Definition at line 292 of file Primary4DVertexValidation.cc.

Referenced by analyze().

maxRank_

constexpr double Primary4DVertexValidation::maxRank_ = 8.

staticprivate

Definition at line 264 of file Primary4DVertexValidation.cc.

Referenced by matchReco2Sim().

maxTry_

constexpr double Primary4DVertexValidation::maxTry_ = 10.

staticprivate

Definition at line 265 of file Primary4DVertexValidation.cc.

Referenced by matchReco2Sim().

meBarrelNoPIDtype_

MonitorElement* Primary4DVertexValidation::meBarrelNoPIDtype_

private

Definition at line 477 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelPIDp_

MonitorElement* Primary4DVertexValidation::meBarrelPIDp_

private

Definition at line 474 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTrueKAsK_

MonitorElement* Primary4DVertexValidation::meBarrelTrueKAsK_

private

Definition at line 495 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTrueKAsP_

MonitorElement* Primary4DVertexValidation::meBarrelTrueKAsP_

private

Definition at line 496 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTrueKAsPi_

MonitorElement* Primary4DVertexValidation::meBarrelTrueKAsPi_

private

Definition at line 494 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTrueKNoPID_

MonitorElement* Primary4DVertexValidation::meBarrelTrueKNoPID_

private

Definition at line 481 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTruePAsK_

MonitorElement* Primary4DVertexValidation::meBarrelTruePAsK_

private

Definition at line 502 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTruePAsP_

MonitorElement* Primary4DVertexValidation::meBarrelTruePAsP_

private

Definition at line 503 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTruePAsPi_

MonitorElement* Primary4DVertexValidation::meBarrelTruePAsPi_

private

Definition at line 501 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTruePiAsK_

MonitorElement* Primary4DVertexValidation::meBarrelTruePiAsK_

private

Definition at line 488 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTruePiAsP_

MonitorElement* Primary4DVertexValidation::meBarrelTruePiAsP_

private

Definition at line 489 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTruePiAsPi_

MonitorElement* Primary4DVertexValidation::meBarrelTruePiAsPi_

private

Definition at line 487 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTruePiNoPID_

MonitorElement* Primary4DVertexValidation::meBarrelTruePiNoPID_

private

Definition at line 480 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTruePNoPID_

MonitorElement* Primary4DVertexValidation::meBarrelTruePNoPID_

private

Definition at line 482 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meDeltaTfakefake_

MonitorElement* Primary4DVertexValidation::meDeltaTfakefake_

private

Definition at line 357 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meDeltaTfakereal_

MonitorElement* Primary4DVertexValidation::meDeltaTfakereal_

private

Definition at line 358 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meDeltaTrealreal_

MonitorElement* Primary4DVertexValidation::meDeltaTrealreal_

private

Definition at line 356 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meDeltaZfakefake_

MonitorElement* Primary4DVertexValidation::meDeltaZfakefake_

private

Definition at line 354 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meDeltaZfakereal_

MonitorElement* Primary4DVertexValidation::meDeltaZfakereal_

private

Definition at line 355 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meDeltaZrealreal_

MonitorElement* Primary4DVertexValidation::meDeltaZrealreal_

private

Definition at line 353 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapNoPIDtype_

MonitorElement* Primary4DVertexValidation::meEndcapNoPIDtype_

private

Definition at line 478 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapPIDp_

MonitorElement* Primary4DVertexValidation::meEndcapPIDp_

private

Definition at line 475 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTrueKAsK_

MonitorElement* Primary4DVertexValidation::meEndcapTrueKAsK_

private

Definition at line 498 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTrueKAsP_

MonitorElement* Primary4DVertexValidation::meEndcapTrueKAsP_

private

Definition at line 499 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTrueKAsPi_

MonitorElement* Primary4DVertexValidation::meEndcapTrueKAsPi_

private

Definition at line 497 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTrueKNoPID_

MonitorElement* Primary4DVertexValidation::meEndcapTrueKNoPID_

private

Definition at line 484 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTruePAsK_

MonitorElement* Primary4DVertexValidation::meEndcapTruePAsK_

private

Definition at line 505 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTruePAsP_

MonitorElement* Primary4DVertexValidation::meEndcapTruePAsP_

private

Definition at line 506 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTruePAsPi_

MonitorElement* Primary4DVertexValidation::meEndcapTruePAsPi_

private

Definition at line 504 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTruePiAsK_

MonitorElement* Primary4DVertexValidation::meEndcapTruePiAsK_

private

Definition at line 491 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTruePiAsP_

MonitorElement* Primary4DVertexValidation::meEndcapTruePiAsP_

private

Definition at line 492 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTruePiAsPi_

MonitorElement* Primary4DVertexValidation::meEndcapTruePiAsPi_

private

Definition at line 490 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTruePiNoPID_

MonitorElement* Primary4DVertexValidation::meEndcapTruePiNoPID_

private

Definition at line 483 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTruePNoPID_

MonitorElement* Primary4DVertexValidation::meEndcapTruePNoPID_

private

Definition at line 485 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meFakeTrackRecLVRelMult_

MonitorElement* Primary4DVertexValidation::meFakeTrackRecLVRelMult_

private

Definition at line 397 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meFakeTrackRecLVRelMultvsMult_

MonitorElement* Primary4DVertexValidation::meFakeTrackRecLVRelMultvsMult_

private

Definition at line 412 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meFakeTrackRecLVRelSumWos_

MonitorElement* Primary4DVertexValidation::meFakeTrackRecLVRelSumWos_

private

Definition at line 404 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meFakeTrackRecLVRelSumWosvsSumWos_

MonitorElement* Primary4DVertexValidation::meFakeTrackRecLVRelSumWosvsSumWos_

private

Definition at line 416 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meFakeTrackRelMult_

MonitorElement* Primary4DVertexValidation::meFakeTrackRelMult_

private

Definition at line 372 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meFakeTrackRelMultvsMult_

MonitorElement* Primary4DVertexValidation::meFakeTrackRelMultvsMult_

private

Definition at line 387 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meFakeTrackRelSumWos_

MonitorElement* Primary4DVertexValidation::meFakeTrackRelSumWos_

private

Definition at line 379 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meFakeTrackRelSumWosvsSumWos_

MonitorElement* Primary4DVertexValidation::meFakeTrackRelSumWosvsSumWos_

private

Definition at line 390 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsFakeRelSumPt2_

MonitorElement* Primary4DVertexValidation::meJetsFakeRelSumPt2_

private

Definition at line 428 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsFakeRelSumPt2vsSumPt2_

MonitorElement* Primary4DVertexValidation::meJetsFakeRelSumPt2vsSumPt2_

private

Definition at line 435 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsPUHt_

MonitorElement* Primary4DVertexValidation::meJetsPUHt_

private

Definition at line 421 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsPUMetPt_

MonitorElement* Primary4DVertexValidation::meJetsPUMetPt_

private

Definition at line 423 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsPUMult_

MonitorElement* Primary4DVertexValidation::meJetsPUMult_

private

Definition at line 420 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsPURelHt_

MonitorElement* Primary4DVertexValidation::meJetsPURelHt_

private

Definition at line 426 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsPURelHtvsHt_

MonitorElement* Primary4DVertexValidation::meJetsPURelHtvsHt_

private

Definition at line 433 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsPURelMetPt_

MonitorElement* Primary4DVertexValidation::meJetsPURelMetPt_

private

Definition at line 429 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsPURelMetPtvsMetPt_

MonitorElement* Primary4DVertexValidation::meJetsPURelMetPtvsMetPt_

private

Definition at line 436 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsPURelMult_

MonitorElement* Primary4DVertexValidation::meJetsPURelMult_

private

Definition at line 425 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsPURelMultvsMult_

MonitorElement* Primary4DVertexValidation::meJetsPURelMultvsMult_

private

Definition at line 432 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsPURelSumPt2_

MonitorElement* Primary4DVertexValidation::meJetsPURelSumPt2_

private

Definition at line 427 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsPURelSumPt2vsSumPt2_

MonitorElement* Primary4DVertexValidation::meJetsPURelSumPt2vsSumPt2_

private

Definition at line 434 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsPURelSumPz_

MonitorElement* Primary4DVertexValidation::meJetsPURelSumPz_

private

Definition at line 430 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsPURelSumPzvsSumPz_

MonitorElement* Primary4DVertexValidation::meJetsPURelSumPzvsSumPz_

private

Definition at line 437 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsPUSumPt2_

MonitorElement* Primary4DVertexValidation::meJetsPUSumPt2_

private

Definition at line 422 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsPUSumPz_

MonitorElement* Primary4DVertexValidation::meJetsPUSumPz_

private

Definition at line 424 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsRecLVFakeRelSumPt2_

MonitorElement* Primary4DVertexValidation::meJetsRecLVFakeRelSumPt2_

private

Definition at line 447 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsRecLVFakeRelSumPt2vsSumPt2_

MonitorElement* Primary4DVertexValidation::meJetsRecLVFakeRelSumPt2vsSumPt2_

private

Definition at line 454 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsRecLVPUHt_

MonitorElement* Primary4DVertexValidation::meJetsRecLVPUHt_

private

Definition at line 440 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsRecLVPUMetPt_

MonitorElement* Primary4DVertexValidation::meJetsRecLVPUMetPt_

private

Definition at line 442 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsRecLVPUMult_

MonitorElement* Primary4DVertexValidation::meJetsRecLVPUMult_

private

Definition at line 439 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsRecLVPURelHt_

MonitorElement* Primary4DVertexValidation::meJetsRecLVPURelHt_

private

Definition at line 445 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsRecLVPURelHtvsHt_

MonitorElement* Primary4DVertexValidation::meJetsRecLVPURelHtvsHt_

private

Definition at line 452 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsRecLVPURelMetPt_

MonitorElement* Primary4DVertexValidation::meJetsRecLVPURelMetPt_

private

Definition at line 448 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsRecLVPURelMetPtvsMetPt_

MonitorElement* Primary4DVertexValidation::meJetsRecLVPURelMetPtvsMetPt_

private

Definition at line 455 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsRecLVPURelMult_

MonitorElement* Primary4DVertexValidation::meJetsRecLVPURelMult_

private

Definition at line 444 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsRecLVPURelMultvsMult_

MonitorElement* Primary4DVertexValidation::meJetsRecLVPURelMultvsMult_

private

Definition at line 451 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsRecLVPURelSumPt2_

MonitorElement* Primary4DVertexValidation::meJetsRecLVPURelSumPt2_

private

Definition at line 446 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsRecLVPURelSumPt2vsSumPt2_

MonitorElement* Primary4DVertexValidation::meJetsRecLVPURelSumPt2vsSumPt2_

private

Definition at line 453 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsRecLVPURelSumPz_

MonitorElement* Primary4DVertexValidation::meJetsRecLVPURelSumPz_

private

Definition at line 449 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsRecLVPURelSumPzvsSumPz_

MonitorElement* Primary4DVertexValidation::meJetsRecLVPURelSumPzvsSumPz_

private

Definition at line 456 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsRecLVPUSumPt2_

MonitorElement* Primary4DVertexValidation::meJetsRecLVPUSumPt2_

private

Definition at line 441 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meJetsRecLVPUSumPz_

MonitorElement* Primary4DVertexValidation::meJetsRecLVPUSumPz_

private

Definition at line 443 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMatchQual_

MonitorElement* Primary4DVertexValidation::meMatchQual_

private

Definition at line 352 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackEffEtaTot_

MonitorElement* Primary4DVertexValidation::meMVATrackEffEtaTot_

private

Definition at line 330 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackEffPtTot_

MonitorElement* Primary4DVertexValidation::meMVATrackEffPtTot_

private

Definition at line 327 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackMatchedEffEtaMtd_

MonitorElement* Primary4DVertexValidation::meMVATrackMatchedEffEtaMtd_

private

Definition at line 332 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackMatchedEffEtaTot_

MonitorElement* Primary4DVertexValidation::meMVATrackMatchedEffEtaTot_

private

Definition at line 331 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackMatchedEffPtMtd_

MonitorElement* Primary4DVertexValidation::meMVATrackMatchedEffPtMtd_

private

Definition at line 329 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackMatchedEffPtTot_

MonitorElement* Primary4DVertexValidation::meMVATrackMatchedEffPtTot_

private

Definition at line 328 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackPullTot_

MonitorElement* Primary4DVertexValidation::meMVATrackPullTot_

private

Definition at line 334 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackResTot_

MonitorElement* Primary4DVertexValidation::meMVATrackResTot_

private

Definition at line 333 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackZposResTot_

MonitorElement* Primary4DVertexValidation::meMVATrackZposResTot_

private

Definition at line 341 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackMult_

MonitorElement* Primary4DVertexValidation::mePUTrackMult_

private

Definition at line 370 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRecLVMult_

MonitorElement* Primary4DVertexValidation::mePUTrackRecLVMult_

private

Definition at line 395 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRecLVRelMult_

MonitorElement* Primary4DVertexValidation::mePUTrackRecLVRelMult_

private

Definition at line 396 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRecLVRelMultvsMult_

MonitorElement* Primary4DVertexValidation::mePUTrackRecLVRelMultvsMult_

private

Definition at line 411 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRecLVRelSumPt2_

MonitorElement* Primary4DVertexValidation::mePUTrackRecLVRelSumPt2_

private

Definition at line 409 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRecLVRelSumPt2vsSumPt2_

MonitorElement* Primary4DVertexValidation::mePUTrackRecLVRelSumPt2vsSumPt2_

private

Definition at line 418 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRecLVRelSumPt_

MonitorElement* Primary4DVertexValidation::mePUTrackRecLVRelSumPt_

private

Definition at line 407 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRecLVRelSumPtvsSumPt_

MonitorElement* Primary4DVertexValidation::mePUTrackRecLVRelSumPtvsSumPt_

private

Definition at line 417 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRecLVRelSumWnt_

MonitorElement* Primary4DVertexValidation::mePUTrackRecLVRelSumWnt_

private

Definition at line 399 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRecLVRelSumWntvsSumWnt_

MonitorElement* Primary4DVertexValidation::mePUTrackRecLVRelSumWntvsSumWnt_

private

Definition at line 413 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRecLVRelSumWos_

MonitorElement* Primary4DVertexValidation::mePUTrackRecLVRelSumWos_

private

Definition at line 401 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRecLVRelSumWosvsSumWos_

MonitorElement* Primary4DVertexValidation::mePUTrackRecLVRelSumWosvsSumWos_

private

Definition at line 414 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRecLVSumPt2_

MonitorElement* Primary4DVertexValidation::mePUTrackRecLVSumPt2_

private

Definition at line 408 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRecLVSumPt_

MonitorElement* Primary4DVertexValidation::mePUTrackRecLVSumPt_

private

Definition at line 406 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRecLVSumWnt_

MonitorElement* Primary4DVertexValidation::mePUTrackRecLVSumWnt_

private

Definition at line 398 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRecLVSumWos_

MonitorElement* Primary4DVertexValidation::mePUTrackRecLVSumWos_

private

Definition at line 400 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRecLVWnt_

MonitorElement* Primary4DVertexValidation::mePUTrackRecLVWnt_

private

Definition at line 405 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRelMult_

MonitorElement* Primary4DVertexValidation::mePUTrackRelMult_

private

Definition at line 371 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRelMultvsMult_

MonitorElement* Primary4DVertexValidation::mePUTrackRelMultvsMult_

private

Definition at line 386 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRelSumPt2_

MonitorElement* Primary4DVertexValidation::mePUTrackRelSumPt2_

private

Definition at line 384 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRelSumPt2vsSumPt2_

MonitorElement* Primary4DVertexValidation::mePUTrackRelSumPt2vsSumPt2_

private

Definition at line 393 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRelSumPt_

MonitorElement* Primary4DVertexValidation::mePUTrackRelSumPt_

private

Definition at line 382 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRelSumPtvsSumPt_

MonitorElement* Primary4DVertexValidation::mePUTrackRelSumPtvsSumPt_

private

Definition at line 392 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRelSumWnt_

MonitorElement* Primary4DVertexValidation::mePUTrackRelSumWnt_

private

Definition at line 374 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRelSumWntvsSumWnt_

MonitorElement* Primary4DVertexValidation::mePUTrackRelSumWntvsSumWnt_

private

Definition at line 388 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRelSumWos_

MonitorElement* Primary4DVertexValidation::mePUTrackRelSumWos_

private

Definition at line 376 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackRelSumWosvsSumWos_

MonitorElement* Primary4DVertexValidation::mePUTrackRelSumWosvsSumWos_

private

Definition at line 389 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackSumPt2_

MonitorElement* Primary4DVertexValidation::mePUTrackSumPt2_

private

Definition at line 383 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackSumPt_

MonitorElement* Primary4DVertexValidation::mePUTrackSumPt_

private

Definition at line 381 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackSumWnt_

MonitorElement* Primary4DVertexValidation::mePUTrackSumWnt_

private

Definition at line 373 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackSumWos_

MonitorElement* Primary4DVertexValidation::mePUTrackSumWos_

private

Definition at line 375 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUTrackWnt_

MonitorElement* Primary4DVertexValidation::mePUTrackWnt_

private

Definition at line 380 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUvsOtherFakeV_

MonitorElement* Primary4DVertexValidation::mePUvsOtherFakeV_

private

Definition at line 350 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUvsRealV_

MonitorElement* Primary4DVertexValidation::mePUvsRealV_

private

Definition at line 349 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUvsSplitV_

MonitorElement* Primary4DVertexValidation::mePUvsSplitV_

private

Definition at line 351 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meRecoPosInRecoOrigCollection_

MonitorElement* Primary4DVertexValidation::meRecoPosInRecoOrigCollection_

private

Definition at line 360 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meRecoPosInSimCollection_

MonitorElement* Primary4DVertexValidation::meRecoPosInSimCollection_

private

Definition at line 359 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meRecoPVPosSignal_

MonitorElement* Primary4DVertexValidation::meRecoPVPosSignal_

private

Definition at line 362 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meRecoPVPosSignalNotHighestPt_

MonitorElement* Primary4DVertexValidation::meRecoPVPosSignalNotHighestPt_

private

Definition at line 363 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meRecoVtxVsLineDensity_

MonitorElement* Primary4DVertexValidation::meRecoVtxVsLineDensity_

private

Definition at line 364 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meRecPVT_

MonitorElement* Primary4DVertexValidation::meRecPVT_

private

Definition at line 367 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meRecPVZ_

MonitorElement* Primary4DVertexValidation::meRecPVZ_

private

Definition at line 366 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meRecVerNumber_

MonitorElement* Primary4DVertexValidation::meRecVerNumber_

private

Definition at line 365 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meSecTrackRecLVRelSumWos_

MonitorElement* Primary4DVertexValidation::meSecTrackRecLVRelSumWos_

private

Definition at line 403 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meSecTrackRecLVRelSumWosvsSumWos_

MonitorElement* Primary4DVertexValidation::meSecTrackRecLVRelSumWosvsSumWos_

private

Definition at line 415 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meSecTrackRecLVSumWos_

MonitorElement* Primary4DVertexValidation::meSecTrackRecLVSumWos_

private

Definition at line 402 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meSecTrackRelSumWos_

MonitorElement* Primary4DVertexValidation::meSecTrackRelSumWos_

private

Definition at line 378 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meSecTrackRelSumWosvsSumWos_

MonitorElement* Primary4DVertexValidation::meSecTrackRelSumWosvsSumWos_

private

Definition at line 391 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meSecTrackSumWos_

MonitorElement* Primary4DVertexValidation::meSecTrackSumWos_

private

Definition at line 377 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meSimPosInSimOrigCollection_

MonitorElement* Primary4DVertexValidation::meSimPosInSimOrigCollection_

private

Definition at line 361 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meSimPVZ_

MonitorElement* Primary4DVertexValidation::meSimPVZ_

private

Definition at line 368 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTimePull_

MonitorElement* Primary4DVertexValidation::meTimePull_

private

Definition at line 346 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTimeRes_

MonitorElement* Primary4DVertexValidation::meTimeRes_

private

Definition at line 345 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTimeSignalPull_

MonitorElement* Primary4DVertexValidation::meTimeSignalPull_

private

Definition at line 348 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTimeSignalRes_

MonitorElement* Primary4DVertexValidation::meTimeSignalRes_

private

Definition at line 347 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrack3DposRes_

MonitorElement* Primary4DVertexValidation::meTrack3DposRes_[3]

private

Definition at line 344 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackPull_

MonitorElement* Primary4DVertexValidation::meTrackPull_[3]

private

Definition at line 338 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackPullHighP_

MonitorElement* Primary4DVertexValidation::meTrackPullHighP_[3]

private

Definition at line 467 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackPullHighPTot_

MonitorElement* Primary4DVertexValidation::meTrackPullHighPTot_

private

Definition at line 462 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackPullLowP_

MonitorElement* Primary4DVertexValidation::meTrackPullLowP_[3]

private

Definition at line 466 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackPullLowPTot_

MonitorElement* Primary4DVertexValidation::meTrackPullLowPTot_

private

Definition at line 461 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackPullTot_

MonitorElement* Primary4DVertexValidation::meTrackPullTot_

private

Definition at line 336 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackRes_

MonitorElement* Primary4DVertexValidation::meTrackRes_[3]

private

Definition at line 337 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResHighP_

MonitorElement* Primary4DVertexValidation::meTrackResHighP_[3]

private

Definition at line 465 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResHighPTot_

MonitorElement* Primary4DVertexValidation::meTrackResHighPTot_

private

Definition at line 460 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResLowP_

MonitorElement* Primary4DVertexValidation::meTrackResLowP_[3]

private

Definition at line 464 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResLowPTot_

MonitorElement* Primary4DVertexValidation::meTrackResLowPTot_

private

Definition at line 459 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResMass_

MonitorElement* Primary4DVertexValidation::meTrackResMass_[3]

private

Definition at line 339 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResMassPions_

MonitorElement* Primary4DVertexValidation::meTrackResMassPions_[3]

private

Definition at line 471 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResMassProtons_

MonitorElement* Primary4DVertexValidation::meTrackResMassProtons_[3]

private

Definition at line 469 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResMassTrue_

MonitorElement* Primary4DVertexValidation::meTrackResMassTrue_[3]

private

Definition at line 340 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResMassTruePions_

MonitorElement* Primary4DVertexValidation::meTrackResMassTruePions_[3]

private

Definition at line 472 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResMassTrueProtons_

MonitorElement* Primary4DVertexValidation::meTrackResMassTrueProtons_[3]

private

Definition at line 470 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResTot_

MonitorElement* Primary4DVertexValidation::meTrackResTot_

private

Definition at line 335 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackZposRes_

MonitorElement* Primary4DVertexValidation::meTrackZposRes_[3]

private

Definition at line 343 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackZposResTot_

MonitorElement* Primary4DVertexValidation::meTrackZposResTot_

private

Definition at line 342 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

minProbHeavy_

const double Primary4DVertexValidation::minProbHeavy_

private

Definition at line 289 of file Primary4DVertexValidation.cc.

Referenced by isParticle().

minThrMetPt_

constexpr double Primary4DVertexValidation::minThrMetPt_ = 1.e-3

staticprivate

Definition at line 279 of file Primary4DVertexValidation.cc.

Referenced by analyze().

minThrSumPt2_

constexpr double Primary4DVertexValidation::minThrSumPt2_ = 1.e-3

staticprivate

Definition at line 278 of file Primary4DVertexValidation.cc.

Referenced by analyze().

minThrSumPt_

constexpr double Primary4DVertexValidation::minThrSumPt_ = 0.01

staticprivate

Definition at line 277 of file Primary4DVertexValidation.cc.

Referenced by analyze().

minThrSumPz_

constexpr double Primary4DVertexValidation::minThrSumPz_ = 1.e-4

staticprivate

Definition at line 280 of file Primary4DVertexValidation.cc.

Referenced by analyze().

minThrSumWnt_

constexpr double Primary4DVertexValidation::minThrSumWnt_ = 0.01

staticprivate

Definition at line 275 of file Primary4DVertexValidation.cc.

Referenced by analyze().

minThrSumWos_

constexpr double Primary4DVertexValidation::minThrSumWos_ = 0.1

staticprivate

Definition at line 276 of file Primary4DVertexValidation.cc.

Referenced by analyze().

momentumToken_

edm::EDGetTokenT<edm::ValueMap<float> > Primary4DVertexValidation::momentumToken_

private

Definition at line 309 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

mvaH_

constexpr double Primary4DVertexValidation::mvaH_ = 0.8

staticprivate

Definition at line 262 of file Primary4DVertexValidation.cc.

Referenced by analyze().

mvaL_

constexpr double Primary4DVertexValidation::mvaL_ = 0.5

staticprivate

Definition at line 261 of file Primary4DVertexValidation.cc.

Referenced by analyze().

mvaTh_

const double Primary4DVertexValidation::mvaTh_

private

Definition at line 291 of file Primary4DVertexValidation.cc.

Referenced by matchReco2Sim().

NOT_MATCHED

constexpr unsigned int Primary4DVertexValidation::NOT_MATCHED = 66666

staticprivate

Definition at line 258 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and matchReco2Sim().

optionalPlots_

bool Primary4DVertexValidation::optionalPlots_

private

Definition at line 286 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

pathLengthToken_

edm::EDGetTokenT<edm::ValueMap<float> > Primary4DVertexValidation::pathLengthToken_

private

Definition at line 308 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

pdtToken_

edm::ESGetToken<HepPDT::ParticleDataTable, edm::DefaultRecord> Primary4DVertexValidation::pdtToken_

private

Definition at line 324 of file Primary4DVertexValidation.cc.

Referenced by analyze().

probKToken_

edm::EDGetTokenT<edm::ValueMap<float> > Primary4DVertexValidation::probKToken_

private

Definition at line 322 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

probPiToken_

edm::EDGetTokenT<edm::ValueMap<float> > Primary4DVertexValidation::probPiToken_

private

Definition at line 321 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

probPToken_

edm::EDGetTokenT<edm::ValueMap<float> > Primary4DVertexValidation::probPToken_

private

Definition at line 323 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

pTcut_

constexpr double Primary4DVertexValidation::pTcut_ = 0.7

staticprivate

Definition at line 269 of file Primary4DVertexValidation.cc.

Referenced by mvaRecSel(), and mvaTPSel().

r2s_

const reco::RecoToSimCollection* Primary4DVertexValidation::r2s_

private

Definition at line 293 of file Primary4DVertexValidation.cc.

Referenced by analyze(), getMatchedTP(), getRecoPVs(), and matchRecoTrack2SimSignal().

Rec4DVerToken_

edm::EDGetTokenT<edm::View<reco::Vertex> > Primary4DVertexValidation::Rec4DVerToken_

private

Definition at line 305 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

RecBeamSpotToken_

edm::EDGetTokenT<reco::BeamSpot> Primary4DVertexValidation::RecBeamSpotToken_

private

Definition at line 304 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

recoToSimAssociationToken_

edm::EDGetTokenT<reco::RecoToSimCollection> Primary4DVertexValidation::recoToSimAssociationToken_

private

Definition at line 303 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

RecTrackToken_

edm::EDGetTokenT<reco::TrackCollection> Primary4DVertexValidation::RecTrackToken_

private

Definition at line 296 of file Primary4DVertexValidation.cc.

Referenced by Primary4DVertexValidation().

s2r_

const reco::SimToRecoCollection* Primary4DVertexValidation::s2r_

private

Definition at line 294 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and getSimPVs().

selNdof_

constexpr double Primary4DVertexValidation::selNdof_ = 4.

staticprivate

Definition at line 263 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and select().

sigmat0SafePidToken_

edm::EDGetTokenT<edm::ValueMap<float> > Primary4DVertexValidation::sigmat0SafePidToken_

private

Definition at line 314 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

simToRecoAssociationToken_

edm::EDGetTokenT<reco::SimToRecoCollection> Primary4DVertexValidation::simToRecoAssociationToken_

private

Definition at line 302 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

simUnit_

constexpr double Primary4DVertexValidation::simUnit_ = 1e9

staticprivate

Definition at line 259 of file Primary4DVertexValidation.cc.

Referenced by analyze(), printMatchedRecoTrackInfo(), and printSimVtxRecoVtxInfo().

t0PidToken_

edm::EDGetTokenT<edm::ValueMap<float> > Primary4DVertexValidation::t0PidToken_

private

Definition at line 312 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

t0SafePidToken_

edm::EDGetTokenT<edm::ValueMap<float> > Primary4DVertexValidation::t0SafePidToken_

private

Definition at line 313 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

timeToken_

edm::EDGetTokenT<edm::ValueMap<float> > Primary4DVertexValidation::timeToken_

private

Definition at line 310 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

tmtdToken_

edm::EDGetTokenT<edm::ValueMap<float> > Primary4DVertexValidation::tmtdToken_

private

Definition at line 317 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

tofKToken_

edm::EDGetTokenT<edm::ValueMap<float> > Primary4DVertexValidation::tofKToken_

private

Definition at line 319 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

tofPiToken_

edm::EDGetTokenT<edm::ValueMap<float> > Primary4DVertexValidation::tofPiToken_

private

Definition at line 318 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

tofPToken_

edm::EDGetTokenT<edm::ValueMap<float> > Primary4DVertexValidation::tofPToken_

private

Definition at line 320 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

tol_

constexpr double Primary4DVertexValidation::tol_ = 1.e-4

staticprivate

Definition at line 274 of file Primary4DVertexValidation.cc.

Referenced by analyze().

trackAssocToken_

edm::EDGetTokenT<edm::ValueMap<int> > Primary4DVertexValidation::trackAssocToken_

private

Definition at line 307 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

trackingParticleCollectionToken_

edm::EDGetTokenT<TrackingParticleCollection> Primary4DVertexValidation::trackingParticleCollectionToken_

private

Definition at line 300 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

trackingVertexCollectionToken_

edm::EDGetTokenT<TrackingVertexCollection> Primary4DVertexValidation::trackingVertexCollectionToken_

private

Definition at line 301 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

trackMaxBtlEta_

constexpr double Primary4DVertexValidation::trackMaxBtlEta_ = 1.5

staticprivate

Definition at line 271 of file Primary4DVertexValidation.cc.

Referenced by analyze().

trackMaxEtlEta_

constexpr double Primary4DVertexValidation::trackMaxEtlEta_ = 3.

staticprivate

Definition at line 273 of file Primary4DVertexValidation.cc.

Referenced by analyze().

trackMinEtlEta_

constexpr double Primary4DVertexValidation::trackMinEtlEta_ = 1.6

staticprivate

Definition at line 272 of file Primary4DVertexValidation.cc.

Referenced by analyze().

trackMVAQualToken_

edm::EDGetTokenT<edm::ValueMap<float> > Primary4DVertexValidation::trackMVAQualToken_

private

Definition at line 315 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

trackweightTh_

const double Primary4DVertexValidation::trackweightTh_

private

Definition at line 290 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and matchReco2Sim().

use3dNoTime_

bool Primary4DVertexValidation::use3dNoTime_

private

Definition at line 287 of file Primary4DVertexValidation.cc.

Referenced by analyze().

use_only_charged_tracks_

bool Primary4DVertexValidation::use_only_charged_tracks_

private

Definition at line 284 of file Primary4DVertexValidation.cc.

Referenced by getSimPVs(), and printSimVtxRecoVtxInfo().

vecPileupSummaryInfoToken_

edm::EDGetTokenT<std::vector<PileupSummaryInfo> > Primary4DVertexValidation::vecPileupSummaryInfoToken_

private

Definition at line 298 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

zWosMatchMax_

constexpr double Primary4DVertexValidation::zWosMatchMax_ = 1.

staticprivate

Definition at line 266 of file Primary4DVertexValidation.cc.

Referenced by matchReco2Sim().