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
const edm::Ref< std::vector< TrackingParticle > > *	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	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 &)
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 *	meMatchQual_
MonitorElement *	meMVATrackEffEtaTot_
MonitorElement *	meMVATrackEffPtTot_
MonitorElement *	meMVATrackMatchedEffEtaMtd_
MonitorElement *	meMVATrackMatchedEffEtaTot_
MonitorElement *	meMVATrackMatchedEffPtMtd_
MonitorElement *	meMVATrackMatchedEffPtTot_
MonitorElement *	meMVATrackPullTot_
MonitorElement *	meMVATrackResTot_
MonitorElement *	meMVATrackZposResTot_
MonitorElement *	mePUvsOtherFakeV_
MonitorElement *	mePUvsRealV_
MonitorElement *	mePUvsSplitV_
MonitorElement *	meRecoPosInRecoOrigCollection_
MonitorElement *	meRecoPosInSimCollection_
MonitorElement *	meRecoPVPosSignal_
MonitorElement *	meRecoPVPosSignalNotHighestPt_
MonitorElement *	meRecoVtxVsLineDensity_
MonitorElement *	meRecPVT_
MonitorElement *	meRecPVZ_
MonitorElement *	meRecVerNumber_
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::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	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	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 54 of file Primary4DVertexValidation.cc.

Member Typedef Documentation

LorentzVector

typedef math::XYZTLorentzVector Primary4DVertexValidation::LorentzVector

private

Definition at line 55 of file Primary4DVertexValidation.cc.

Constructor & Destructor Documentation

Primary4DVertexValidation()

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

explicit

Definition at line 369 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")),
      minProbHeavy_(iConfig.getParameter<double>("minProbHeavy")),
      trackweightTh_(iConfig.getParameter<double>("trackweightTh")),
      mvaTh_(iConfig.getParameter<double>("mvaTh")),
      lineDensityPar_(iConfig.getParameter<std::vector<double>>("lineDensityPar")) {
  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 406 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 1159 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(), getMatchedTP(), getRecoPVs(), getSimPVs(), caHitNtupletGeneratorKernels::iev, iEvent, CommonMethods::isnan(), edm::HandleBase::isValid(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::iv, lineDensityPar_, LogTrace, EgHLTOffHistBins_cfi::mass, matchReco2Sim(), 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_, meMatchQual_, meMVATrackEffEtaTot_, meMVATrackEffPtTot_, meMVATrackMatchedEffEtaMtd_, meMVATrackMatchedEffEtaTot_, meMVATrackMatchedEffPtMtd_, meMVATrackMatchedEffPtTot_, meMVATrackPullTot_, meMVATrackResTot_, meMVATrackZposResTot_, mePUvsOtherFakeV_, mePUvsRealV_, mePUvsSplitV_, meRecoPosInRecoOrigCollection_, meRecoPosInSimCollection_, meRecoPVPosSignal_, meRecoPVPosSignalNotHighestPt_, meRecoVtxVsLineDensity_, meRecPVT_, meRecPVZ_, meRecVerNumber_, 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_, minProbHeavy_, momentumToken_, mvaH_, mvaL_, mvaRecSel(), mvaTPSel(), NOT_MATCHED, optionalPlots_, pathLengthToken_, probKToken_, probPiToken_, probPToken_, edm::Handle< T >::product(), Primary4DVertexValidation::simPrimaryVertex::ptsq, r2s_, Rec4DVerToken_, RecBeamSpotToken_, recoToSimAssociationToken_, s2r_, select(), selNdof_, sigmat0SafePidToken_, simToRecoAssociationToken_, simUnit_, submitPVValidationJobs::split(), mathSSE::sqrt(), t0PidToken_, t0SafePidToken_, hcalRecHitTable_cff::time, timeFromTrueMass(), timeToken_, tmtdToken_, tofKToken_, tofPiToken_, tofPToken_, tol_, trackAssocToken_, trackingParticleCollectionToken_, trackingVertexCollectionToken_, trackMaxBtlEta_, trackMaxEtlEta_, trackMinEtlEta_, trackMVAQualToken_, trackweightTh_, 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_);

  //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 &&
                                   simpv.at(iev).eventId.bunchCrossing() == 0 && simpv.at(iev).eventId.event() == 0 &&
                                   recopv.at(iv).OriginalIndex == 0;
        if (selectedVtxMatching && !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();
        }
        double vzsim = simpv.at(iev).z;
        double vtsim = simpv.at(iev).t * simUnit_;

        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 (selectedVtxMatching && selectRecoTrk) {
            if (noCrack) {
              meMVATrackEffPtTot_->Fill((*iTrack)->pt());
            }
            meMVATrackEffEtaTot_->Fill(std::abs((*iTrack)->eta()));
          }

          auto tp_info = getMatchedTP(*iTrack, vsim);
          if (tp_info != nullptr) {
            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 (selectedVtxMatching && 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 (selectedVtxMatching && 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;
              if (probPi[*iTrack] == -1) {
                noPIDtype = 1;
              } else if (std::isnan(probPi[*iTrack])) {
                noPIDtype = 2;
              } else if (probPi[*iTrack] == 1 && probK[*iTrack] == 0 && probP[*iTrack] == 0) {
                noPIDtype = 3;
              }
              bool noPID = noPIDtype > 0;
              bool isPi = !noPID && 1. - probPi[*iTrack] < minProbHeavy_;
              bool isK = !noPID && !isPi && probK[*iTrack] > probP[*iTrack];
              bool isP = !noPID && !isPi && !isK;

              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
        }
      }
    }  // ndof
  }

  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 (std::isinf(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 411 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_, meMatchQual_, meMVATrackEffEtaTot_, meMVATrackEffPtTot_, meMVATrackMatchedEffEtaMtd_, meMVATrackMatchedEffEtaTot_, meMVATrackMatchedEffPtMtd_, meMVATrackMatchedEffPtTot_, meMVATrackPullTot_, meMVATrackResTot_, meMVATrackZposResTot_, mePUvsOtherFakeV_, mePUvsRealV_, mePUvsSplitV_, meRecoPosInRecoOrigCollection_, meRecoPosInSimCollection_, meRecoPVPosSignal_, meRecoPVPosSignalNotHighestPt_, meRecoVtxVsLineDensity_, meRecPVT_, meRecPVZ_, meRecVerNumber_, 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", "Pt 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", "Pt of tracks associated to LV matched to TP; track eta ", 66, 0., 3.3);
  meMVATrackMatchedEffEtaMtd_ = ibook.book1D(
      "MVAMatchedEffEtaMtd", "Pt 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", 200, 0, 200);
  meRecoPVPosSignalNotHighestPt_ =
      ibook.book1D("RecoPVPosSignalNotHighestPt",
                   "Position in reco collection of PV associated to sim signal not highest Pt",
                   200,
                   0,
                   200);
  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.);

  //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.);

  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 NoPIDpi 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 1717 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<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()

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

private

Definition at line 761 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
  if (found == r2s_->end())
    return nullptr;

  //matched TP equal to any TP of 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 &tp.first;
  }

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

getRecoPVs()

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

private

Definition at line 916 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 815 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::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->ptsq += ((**iTP).pt() * (**iTP).pt());

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

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 980 of file Primary4DVertexValidation.cc.

References funct::abs(), c_, funct::cos(), debug_, getMatchedTP(), mps_fire::i, caHitNtupletGeneratorKernels::iev, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::iv, dqmdumpme::k, maxRank_, maxTry_, SiStripPI::min, mvaTh_, NOT_MATCHED, funct::pow(), DiDispStaMuonMonitor_cfi::pt, beamSpotPI::sigmaZ, funct::sin(), pfDeepBoostedJetPreprocessParams_cfi::sv, funct::tan(), 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) {
        double pt = (*iTrack)->pt();

        if (vertex->trackWeight(*iTrack) < trackweightTh_)
          continue;
        if (MVA[(*iTrack)] < mvaTh_)
          continue;

        auto tp_info = getMatchedTP(*iTrack, simpv.at(iev).sim_vertex);
        if (tp_info != nullptr) {
          double dz2_beam = pow((*BS).BeamWidthX() * cos((*iTrack)->phi()) / tan((*iTrack)->theta()), 2) +
                            pow((*BS).BeamWidthY() * sin((*iTrack)->phi()) / tan((*iTrack)->theta()), 2);
          double dz2 = pow((*iTrack)->dzError(), 2) + dz2_beam +
                       pow(0.0020, 2);  // added 20 um, some tracks have crazy small resolutions
          wos = vertex->trackWeight(*iTrack) / dz2;
          wnt = vertex->trackWeight(*iTrack) * std::min(pt, 1.0);

          if (sigmat0[(*iTrack)] > 0) {
            double sigmaZ = (*BS).sigmaZ();
            double sigmaT = sigmaZ / c_;  // c in cm/ns
            wos = wos / erf(sigmat0[(*iTrack)] / sigmaT);
          }
          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 744 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 1769 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 1760 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;
}

select()

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

private

Definition at line 793 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 781 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 215 of file Primary4DVertexValidation.cc.

Referenced by matchReco2Sim(), and timeFromTrueMass().

c_cm_ns

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

staticprivate

Definition at line 231 of file Primary4DVertexValidation.cc.

debug_

bool Primary4DVertexValidation::debug_

private

Definition at line 234 of file Primary4DVertexValidation.cc.

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

deltaZcut_

constexpr double Primary4DVertexValidation::deltaZcut_ = 0.1

staticprivate

Definition at line 225 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and mvaRecSel().

etacutGEN_

constexpr double Primary4DVertexValidation::etacutGEN_ = 4.

staticprivate

Definition at line 222 of file Primary4DVertexValidation.cc.

Referenced by mvaTPSel().

etacutREC_

constexpr double Primary4DVertexValidation::etacutREC_ = 3.

staticprivate

Definition at line 223 of file Primary4DVertexValidation.cc.

Referenced by mvaRecSel().

folder_

const std::string Primary4DVertexValidation::folder_

private

Definition at line 212 of file Primary4DVertexValidation.cc.

Referenced by bookHistograms().

lineDensityPar_

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

private

Definition at line 240 of file Primary4DVertexValidation.cc.

Referenced by analyze().

maxRank_

constexpr double Primary4DVertexValidation::maxRank_ = 8.

staticprivate

Definition at line 219 of file Primary4DVertexValidation.cc.

Referenced by matchReco2Sim().

maxTry_

constexpr double Primary4DVertexValidation::maxTry_ = 10.

staticprivate

Definition at line 220 of file Primary4DVertexValidation.cc.

Referenced by matchReco2Sim().

meBarrelNoPIDtype_

MonitorElement* Primary4DVertexValidation::meBarrelNoPIDtype_

private

Definition at line 336 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelPIDp_

MonitorElement* Primary4DVertexValidation::meBarrelPIDp_

private

Definition at line 333 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTrueKAsK_

MonitorElement* Primary4DVertexValidation::meBarrelTrueKAsK_

private

Definition at line 354 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTrueKAsP_

MonitorElement* Primary4DVertexValidation::meBarrelTrueKAsP_

private

Definition at line 355 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTrueKAsPi_

MonitorElement* Primary4DVertexValidation::meBarrelTrueKAsPi_

private

Definition at line 353 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTrueKNoPID_

MonitorElement* Primary4DVertexValidation::meBarrelTrueKNoPID_

private

Definition at line 340 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTruePAsK_

MonitorElement* Primary4DVertexValidation::meBarrelTruePAsK_

private

Definition at line 361 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTruePAsP_

MonitorElement* Primary4DVertexValidation::meBarrelTruePAsP_

private

Definition at line 362 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTruePAsPi_

MonitorElement* Primary4DVertexValidation::meBarrelTruePAsPi_

private

Definition at line 360 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTruePiAsK_

MonitorElement* Primary4DVertexValidation::meBarrelTruePiAsK_

private

Definition at line 347 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTruePiAsP_

MonitorElement* Primary4DVertexValidation::meBarrelTruePiAsP_

private

Definition at line 348 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTruePiAsPi_

MonitorElement* Primary4DVertexValidation::meBarrelTruePiAsPi_

private

Definition at line 346 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTruePiNoPID_

MonitorElement* Primary4DVertexValidation::meBarrelTruePiNoPID_

private

Definition at line 339 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meBarrelTruePNoPID_

MonitorElement* Primary4DVertexValidation::meBarrelTruePNoPID_

private

Definition at line 341 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meDeltaTfakefake_

MonitorElement* Primary4DVertexValidation::meDeltaTfakefake_

private

Definition at line 304 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meDeltaTfakereal_

MonitorElement* Primary4DVertexValidation::meDeltaTfakereal_

private

Definition at line 305 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meDeltaTrealreal_

MonitorElement* Primary4DVertexValidation::meDeltaTrealreal_

private

Definition at line 303 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meDeltaZfakefake_

MonitorElement* Primary4DVertexValidation::meDeltaZfakefake_

private

Definition at line 301 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meDeltaZfakereal_

MonitorElement* Primary4DVertexValidation::meDeltaZfakereal_

private

Definition at line 302 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meDeltaZrealreal_

MonitorElement* Primary4DVertexValidation::meDeltaZrealreal_

private

Definition at line 300 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapNoPIDtype_

MonitorElement* Primary4DVertexValidation::meEndcapNoPIDtype_

private

Definition at line 337 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapPIDp_

MonitorElement* Primary4DVertexValidation::meEndcapPIDp_

private

Definition at line 334 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTrueKAsK_

MonitorElement* Primary4DVertexValidation::meEndcapTrueKAsK_

private

Definition at line 357 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTrueKAsP_

MonitorElement* Primary4DVertexValidation::meEndcapTrueKAsP_

private

Definition at line 358 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTrueKAsPi_

MonitorElement* Primary4DVertexValidation::meEndcapTrueKAsPi_

private

Definition at line 356 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTrueKNoPID_

MonitorElement* Primary4DVertexValidation::meEndcapTrueKNoPID_

private

Definition at line 343 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTruePAsK_

MonitorElement* Primary4DVertexValidation::meEndcapTruePAsK_

private

Definition at line 364 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTruePAsP_

MonitorElement* Primary4DVertexValidation::meEndcapTruePAsP_

private

Definition at line 365 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTruePAsPi_

MonitorElement* Primary4DVertexValidation::meEndcapTruePAsPi_

private

Definition at line 363 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTruePiAsK_

MonitorElement* Primary4DVertexValidation::meEndcapTruePiAsK_

private

Definition at line 350 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTruePiAsP_

MonitorElement* Primary4DVertexValidation::meEndcapTruePiAsP_

private

Definition at line 351 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTruePiAsPi_

MonitorElement* Primary4DVertexValidation::meEndcapTruePiAsPi_

private

Definition at line 349 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTruePiNoPID_

MonitorElement* Primary4DVertexValidation::meEndcapTruePiNoPID_

private

Definition at line 342 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meEndcapTruePNoPID_

MonitorElement* Primary4DVertexValidation::meEndcapTruePNoPID_

private

Definition at line 344 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMatchQual_

MonitorElement* Primary4DVertexValidation::meMatchQual_

private

Definition at line 299 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackEffEtaTot_

MonitorElement* Primary4DVertexValidation::meMVATrackEffEtaTot_

private

Definition at line 277 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackEffPtTot_

MonitorElement* Primary4DVertexValidation::meMVATrackEffPtTot_

private

Definition at line 274 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackMatchedEffEtaMtd_

MonitorElement* Primary4DVertexValidation::meMVATrackMatchedEffEtaMtd_

private

Definition at line 279 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackMatchedEffEtaTot_

MonitorElement* Primary4DVertexValidation::meMVATrackMatchedEffEtaTot_

private

Definition at line 278 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackMatchedEffPtMtd_

MonitorElement* Primary4DVertexValidation::meMVATrackMatchedEffPtMtd_

private

Definition at line 276 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackMatchedEffPtTot_

MonitorElement* Primary4DVertexValidation::meMVATrackMatchedEffPtTot_

private

Definition at line 275 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackPullTot_

MonitorElement* Primary4DVertexValidation::meMVATrackPullTot_

private

Definition at line 281 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackResTot_

MonitorElement* Primary4DVertexValidation::meMVATrackResTot_

private

Definition at line 280 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackZposResTot_

MonitorElement* Primary4DVertexValidation::meMVATrackZposResTot_

private

Definition at line 288 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUvsOtherFakeV_

MonitorElement* Primary4DVertexValidation::mePUvsOtherFakeV_

private

Definition at line 297 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUvsRealV_

MonitorElement* Primary4DVertexValidation::mePUvsRealV_

private

Definition at line 296 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

mePUvsSplitV_

MonitorElement* Primary4DVertexValidation::mePUvsSplitV_

private

Definition at line 298 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meRecoPosInRecoOrigCollection_

MonitorElement* Primary4DVertexValidation::meRecoPosInRecoOrigCollection_

private

Definition at line 307 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meRecoPosInSimCollection_

MonitorElement* Primary4DVertexValidation::meRecoPosInSimCollection_

private

Definition at line 306 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meRecoPVPosSignal_

MonitorElement* Primary4DVertexValidation::meRecoPVPosSignal_

private

Definition at line 309 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meRecoPVPosSignalNotHighestPt_

MonitorElement* Primary4DVertexValidation::meRecoPVPosSignalNotHighestPt_

private

Definition at line 310 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meRecoVtxVsLineDensity_

MonitorElement* Primary4DVertexValidation::meRecoVtxVsLineDensity_

private

Definition at line 311 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meRecPVT_

MonitorElement* Primary4DVertexValidation::meRecPVT_

private

Definition at line 314 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meRecPVZ_

MonitorElement* Primary4DVertexValidation::meRecPVZ_

private

Definition at line 313 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meRecVerNumber_

MonitorElement* Primary4DVertexValidation::meRecVerNumber_

private

Definition at line 312 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meSimPosInSimOrigCollection_

MonitorElement* Primary4DVertexValidation::meSimPosInSimOrigCollection_

private

Definition at line 308 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meSimPVZ_

MonitorElement* Primary4DVertexValidation::meSimPVZ_

private

Definition at line 315 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTimePull_

MonitorElement* Primary4DVertexValidation::meTimePull_

private

Definition at line 293 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTimeRes_

MonitorElement* Primary4DVertexValidation::meTimeRes_

private

Definition at line 292 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTimeSignalPull_

MonitorElement* Primary4DVertexValidation::meTimeSignalPull_

private

Definition at line 295 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTimeSignalRes_

MonitorElement* Primary4DVertexValidation::meTimeSignalRes_

private

Definition at line 294 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrack3DposRes_

MonitorElement* Primary4DVertexValidation::meTrack3DposRes_[3]

private

Definition at line 291 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackPull_

MonitorElement* Primary4DVertexValidation::meTrackPull_[3]

private

Definition at line 285 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackPullHighP_

MonitorElement* Primary4DVertexValidation::meTrackPullHighP_[3]

private

Definition at line 326 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackPullHighPTot_

MonitorElement* Primary4DVertexValidation::meTrackPullHighPTot_

private

Definition at line 321 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackPullLowP_

MonitorElement* Primary4DVertexValidation::meTrackPullLowP_[3]

private

Definition at line 325 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackPullLowPTot_

MonitorElement* Primary4DVertexValidation::meTrackPullLowPTot_

private

Definition at line 320 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackPullTot_

MonitorElement* Primary4DVertexValidation::meTrackPullTot_

private

Definition at line 283 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackRes_

MonitorElement* Primary4DVertexValidation::meTrackRes_[3]

private

Definition at line 284 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResHighP_

MonitorElement* Primary4DVertexValidation::meTrackResHighP_[3]

private

Definition at line 324 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResHighPTot_

MonitorElement* Primary4DVertexValidation::meTrackResHighPTot_

private

Definition at line 319 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResLowP_

MonitorElement* Primary4DVertexValidation::meTrackResLowP_[3]

private

Definition at line 323 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResLowPTot_

MonitorElement* Primary4DVertexValidation::meTrackResLowPTot_

private

Definition at line 318 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResMass_

MonitorElement* Primary4DVertexValidation::meTrackResMass_[3]

private

Definition at line 286 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResMassPions_

MonitorElement* Primary4DVertexValidation::meTrackResMassPions_[3]

private

Definition at line 330 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResMassProtons_

MonitorElement* Primary4DVertexValidation::meTrackResMassProtons_[3]

private

Definition at line 328 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResMassTrue_

MonitorElement* Primary4DVertexValidation::meTrackResMassTrue_[3]

private

Definition at line 287 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResMassTruePions_

MonitorElement* Primary4DVertexValidation::meTrackResMassTruePions_[3]

private

Definition at line 331 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResMassTrueProtons_

MonitorElement* Primary4DVertexValidation::meTrackResMassTrueProtons_[3]

private

Definition at line 329 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackResTot_

MonitorElement* Primary4DVertexValidation::meTrackResTot_

private

Definition at line 282 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackZposRes_

MonitorElement* Primary4DVertexValidation::meTrackZposRes_[3]

private

Definition at line 290 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackZposResTot_

MonitorElement* Primary4DVertexValidation::meTrackZposResTot_

private

Definition at line 289 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

minProbHeavy_

const double Primary4DVertexValidation::minProbHeavy_

private

Definition at line 237 of file Primary4DVertexValidation.cc.

Referenced by analyze().

momentumToken_

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

private

Definition at line 257 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

mvaH_

constexpr double Primary4DVertexValidation::mvaH_ = 0.8

staticprivate

Definition at line 217 of file Primary4DVertexValidation.cc.

Referenced by analyze().

mvaL_

constexpr double Primary4DVertexValidation::mvaL_ = 0.5

staticprivate

Definition at line 216 of file Primary4DVertexValidation.cc.

Referenced by analyze().

mvaTh_

const double Primary4DVertexValidation::mvaTh_

private

Definition at line 239 of file Primary4DVertexValidation.cc.

Referenced by matchReco2Sim().

NOT_MATCHED

constexpr unsigned int Primary4DVertexValidation::NOT_MATCHED = 66666

staticprivate

Definition at line 213 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and matchReco2Sim().

optionalPlots_

bool Primary4DVertexValidation::optionalPlots_

private

Definition at line 235 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and bookHistograms().

pathLengthToken_

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

private

Definition at line 256 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

probKToken_

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

private

Definition at line 270 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

probPiToken_

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

private

Definition at line 269 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

probPToken_

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

private

Definition at line 271 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

pTcut_

constexpr double Primary4DVertexValidation::pTcut_ = 0.7

staticprivate

Definition at line 224 of file Primary4DVertexValidation.cc.

Referenced by mvaRecSel(), and mvaTPSel().

r2s_

const reco::RecoToSimCollection* Primary4DVertexValidation::r2s_

private

Definition at line 241 of file Primary4DVertexValidation.cc.

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

Rec4DVerToken_

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

private

Definition at line 253 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

RecBeamSpotToken_

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

private

Definition at line 252 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

recoToSimAssociationToken_

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

private

Definition at line 251 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

RecTrackToken_

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

private

Definition at line 244 of file Primary4DVertexValidation.cc.

Referenced by Primary4DVertexValidation().

s2r_

const reco::SimToRecoCollection* Primary4DVertexValidation::s2r_

private

Definition at line 242 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and getSimPVs().

selNdof_

constexpr double Primary4DVertexValidation::selNdof_ = 4.

staticprivate

Definition at line 218 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and select().

sigmat0SafePidToken_

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

private

Definition at line 262 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

simToRecoAssociationToken_

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

private

Definition at line 250 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

simUnit_

constexpr double Primary4DVertexValidation::simUnit_ = 1e9

staticprivate

Definition at line 214 of file Primary4DVertexValidation.cc.

Referenced by analyze().

t0PidToken_

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

private

Definition at line 260 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

t0SafePidToken_

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

private

Definition at line 261 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

timeToken_

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

private

Definition at line 258 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

tmtdToken_

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

private

Definition at line 265 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

tofKToken_

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

private

Definition at line 267 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

tofPiToken_

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

private

Definition at line 266 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

tofPToken_

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

private

Definition at line 268 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

tol_

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

staticprivate

Definition at line 229 of file Primary4DVertexValidation.cc.

Referenced by analyze().

trackAssocToken_

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

private

Definition at line 255 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

trackingParticleCollectionToken_

edm::EDGetTokenT<TrackingParticleCollection> Primary4DVertexValidation::trackingParticleCollectionToken_

private

Definition at line 248 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

trackingVertexCollectionToken_

edm::EDGetTokenT<TrackingVertexCollection> Primary4DVertexValidation::trackingVertexCollectionToken_

private

Definition at line 249 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

trackMaxBtlEta_

constexpr double Primary4DVertexValidation::trackMaxBtlEta_ = 1.5

staticprivate

Definition at line 226 of file Primary4DVertexValidation.cc.

Referenced by analyze().

trackMaxEtlEta_

constexpr double Primary4DVertexValidation::trackMaxEtlEta_ = 3.

staticprivate

Definition at line 228 of file Primary4DVertexValidation.cc.

Referenced by analyze().

trackMinEtlEta_

constexpr double Primary4DVertexValidation::trackMinEtlEta_ = 1.6

staticprivate

Definition at line 227 of file Primary4DVertexValidation.cc.

Referenced by analyze().

trackMVAQualToken_

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

private

Definition at line 263 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

trackweightTh_

const double Primary4DVertexValidation::trackweightTh_

private

Definition at line 238 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and matchReco2Sim().

use_only_charged_tracks_

bool Primary4DVertexValidation::use_only_charged_tracks_

private

Definition at line 233 of file Primary4DVertexValidation.cc.

Referenced by getSimPVs().

vecPileupSummaryInfoToken_

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

private

Definition at line 246 of file Primary4DVertexValidation.cc.

Referenced by analyze(), and Primary4DVertexValidation().

zWosMatchMax_

constexpr double Primary4DVertexValidation::zWosMatchMax_ = 1.

staticprivate

Definition at line 221 of file Primary4DVertexValidation.cc.

Referenced by matchReco2Sim().