MtdTracksValidation Class Reference

Inheritance diagram for MtdTracksValidation:

Public Member Functions
	MtdTracksValidation (const edm::ParameterSet &)
	~MtdTracksValidation () 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 Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
const edm::Ref< std::vector< TrackingParticle > > *	getMatchedTP (const reco::TrackBaseRef &, const double &)
bool	isETL (const double eta) const
const bool	mvaGenRecMatch (const HepMC::GenParticle &, const double &, const reco::TrackBase &, const bool &)
const bool	mvaGenSel (const HepMC::GenParticle &, const float &)
const bool	mvaRecSel (const reco::TrackBase &, const reco::Vertex &, const double &, const double &)
const bool	mvaTPSel (const TrackingParticle &)
const bool	tpWithMTD (const TrackingParticle &, const std::unordered_set< unsigned long int > &)
const unsigned long int	uniqueId (const uint32_t x, const EncodedEventId &y)

Private Attributes
edm::EDGetTokenT< CrossingFrame< PSimHit > >	btlSimHitsToken_
edm::EDGetTokenT< CrossingFrame< PSimHit > >	etlSimHitsToken_
const std::string	folder_
edm::EDGetTokenT< reco::TrackCollection >	GenRecTrackToken_
edm::EDGetTokenT< edm::HepMCProduct >	HepMCProductToken_
MonitorElement *	meBTLTrackEffEtaMtd_
MonitorElement *	meBTLTrackEffEtaTot_
MonitorElement *	meBTLTrackEffPhiMtd_
MonitorElement *	meBTLTrackEffPhiTot_
MonitorElement *	meBTLTrackEffPtMtd_
MonitorElement *	meBTLTrackEffPtTot_
MonitorElement *	meBTLTrackPtRes_
MonitorElement *	meBTLTrackRPTime_
MonitorElement *	meETLTrackEffEta2Mtd_ [2]
MonitorElement *	meETLTrackEffEtaMtd_ [2]
MonitorElement *	meETLTrackEffEtaTot_ [2]
MonitorElement *	meETLTrackEffPhi2Mtd_ [2]
MonitorElement *	meETLTrackEffPhiMtd_ [2]
MonitorElement *	meETLTrackEffPhiTot_ [2]
MonitorElement *	meETLTrackEffPt2Mtd_ [2]
MonitorElement *	meETLTrackEffPtMtd_ [2]
MonitorElement *	meETLTrackEffPtTot_ [2]
MonitorElement *	meETLTrackPtRes_
MonitorElement *	meETLTrackRPTime_
MonitorElement *	meMVATrackEffEtaTot_
MonitorElement *	meMVATrackEffPtTot_
MonitorElement *	meMVATrackMatchedEffEtaMtd_
MonitorElement *	meMVATrackMatchedEffEtaTot_
MonitorElement *	meMVATrackMatchedEffPtMtd_
MonitorElement *	meMVATrackMatchedEffPtTot_
MonitorElement *	meMVATrackPullTot_
MonitorElement *	meMVATrackResTot_
MonitorElement *	meMVATrackZposResTot_
MonitorElement *	meNTrackingParticles_
MonitorElement *	meTrackMatchedTPEffEtaEtl2Mtd_
MonitorElement *	meTrackMatchedTPEffEtaMtd_
MonitorElement *	meTrackMatchedTPEffEtaTot_
MonitorElement *	meTrackMatchedTPEffPtEtl2Mtd_
MonitorElement *	meTrackMatchedTPEffPtMtd_
MonitorElement *	meTrackMatchedTPEffPtTot_
MonitorElement *	meTrackMatchedTPmtdEffEtaMtd_
MonitorElement *	meTrackMatchedTPmtdEffEtaTot_
MonitorElement *	meTrackMatchedTPmtdEffPtMtd_
MonitorElement *	meTrackMatchedTPmtdEffPtTot_
MonitorElement *	meTrackMVAQual_
MonitorElement *	meTrackNumHits_
MonitorElement *	meTrackPathLenghtvsEta_
MonitorElement *	meTrackSigmat0Pid_
MonitorElement *	meTrackSigmat0SafePid_
MonitorElement *	meTrackSigmat0Src_
MonitorElement *	meTrackt0Pid_
MonitorElement *	meTrackt0SafePid_
MonitorElement *	meTrackt0Src_
MonitorElement *	meTracktmtd_
MonitorElement *	meUnassCrysEnergy_
MonitorElement *	meUnassDeposit_
MonitorElement *	meUnassLgadsEnergy_
MonitorElement *	meUnassociatedDetId_
edm::ESGetToken< MTDGeometry, MTDDigiGeometryRecord >	mtdgeoToken_
edm::ESGetToken< MTDTopology, MTDTopologyRcd >	mtdtopoToken_
bool	optionalPlots_
edm::ESGetToken< HepPDT::ParticleDataTable, edm::DefaultRecord >	particleTableToken_
edm::EDGetTokenT< edm::ValueMap< float > >	pathLengthToken_
const reco::RecoToSimCollection *	r2s_
edm::EDGetTokenT< reco::RecoToSimCollection >	recoToSimAssociationToken_
edm::EDGetTokenT< reco::TrackCollection >	RecTrackToken_
edm::EDGetTokenT< std::vector< reco::Vertex > >	RecVertexToken_
const reco::SimToRecoCollection *	s2r_
edm::EDGetTokenT< edm::ValueMap< float > >	Sigmat0PidToken_
edm::EDGetTokenT< edm::ValueMap< float > >	Sigmat0SafePidToken_
edm::EDGetTokenT< edm::ValueMap< float > >	Sigmat0SrcToken_
edm::EDGetTokenT< edm::ValueMap< float > >	SigmatmtdToken_
edm::EDGetTokenT< reco::SimToRecoCollection >	simToRecoAssociationToken_
edm::EDGetTokenT< edm::ValueMap< float > >	t0PidToken_
edm::EDGetTokenT< edm::ValueMap< float > >	t0SafePidToken_
edm::EDGetTokenT< edm::ValueMap< float > >	t0SrcToken_
edm::EDGetTokenT< edm::ValueMap< float > >	tmtdToken_
edm::EDGetTokenT< edm::ValueMap< int > >	trackAssocToken_
edm::EDGetTokenT< TrackingParticleCollection >	trackingParticleCollectionToken_
const float	trackMaxBtlEta_
const float	trackMaxEtlEta_
const float	trackMinEtlEta_
const float	trackMinPt_
edm::EDGetTokenT< edm::ValueMap< float > >	trackMVAQualToken_

Static Private Attributes
static constexpr double	deltaDRcut_ = 0.03
static constexpr double	deltaPTcut_ = 0.05
static constexpr double	deltaZcut_ = 0.1
static constexpr double	depositBTLthreshold_ = 1
static constexpr double	depositETLthreshold_ = 0.001
static constexpr double	etacutGEN_ = 4.
static constexpr double	etacutREC_ = 3.
static constexpr double	etaMatchCut_ = 0.05
static constexpr double	pTcut_ = 0.7
static constexpr double	rBTL_ = 110.0
static constexpr double	zETL_ = 290.0

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 49 of file MtdTracksValidation.cc.

Constructor & Destructor Documentation

MtdTracksValidation()

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

explicit

Definition at line 193 of file MtdTracksValidation.cc.

References btlSimHitsToken_, etlSimHitsToken_, GenRecTrackToken_, edm::ParameterSet::getParameter(), HepMCProductToken_, mtdgeoToken_, mtdtopoToken_, particleTableToken_, pathLengthToken_, recoToSimAssociationToken_, RecTrackToken_, RecVertexToken_, Sigmat0PidToken_, Sigmat0SafePidToken_, Sigmat0SrcToken_, SigmatmtdToken_, simToRecoAssociationToken_, t0PidToken_, t0SafePidToken_, t0SrcToken_, tmtdToken_, trackAssocToken_, trackingParticleCollectionToken_, and trackMVAQualToken_.

    : folder_(iConfig.getParameter<std::string>("folder")),
      trackMinPt_(iConfig.getParameter<double>("trackMinimumPt")),
      trackMaxBtlEta_(iConfig.getParameter<double>("trackMaximumBtlEta")),
      trackMinEtlEta_(iConfig.getParameter<double>("trackMinimumEtlEta")),
      trackMaxEtlEta_(iConfig.getParameter<double>("trackMaximumEtlEta")),
      optionalPlots_(iConfig.getUntrackedParameter<bool>("optionalPlots")) {
  GenRecTrackToken_ = consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("inputTagG"));
  RecTrackToken_ = consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("inputTagT"));
  RecVertexToken_ = consumes<std::vector<reco::Vertex>>(iConfig.getParameter<edm::InputTag>("inputTagV"));
  HepMCProductToken_ = consumes<edm::HepMCProduct>(iConfig.getParameter<edm::InputTag>("inputTagH"));
  trackingParticleCollectionToken_ =
      consumes<TrackingParticleCollection>(iConfig.getParameter<edm::InputTag>("SimTag"));
  simToRecoAssociationToken_ =
      consumes<reco::SimToRecoCollection>(iConfig.getParameter<edm::InputTag>("TPtoRecoTrackAssoc"));
  recoToSimAssociationToken_ =
      consumes<reco::RecoToSimCollection>(iConfig.getParameter<edm::InputTag>("TPtoRecoTrackAssoc"));
  btlSimHitsToken_ = consumes<CrossingFrame<PSimHit>>(iConfig.getParameter<edm::InputTag>("btlSimHits"));
  etlSimHitsToken_ = consumes<CrossingFrame<PSimHit>>(iConfig.getParameter<edm::InputTag>("etlSimHits"));
  trackAssocToken_ = consumes<edm::ValueMap<int>>(iConfig.getParameter<edm::InputTag>("trackAssocSrc"));
  pathLengthToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("pathLengthSrc"));
  tmtdToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("tmtd"));
  SigmatmtdToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("sigmatmtd"));
  t0SrcToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("t0Src"));
  Sigmat0SrcToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("sigmat0Src"));
  t0PidToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("t0PID"));
  Sigmat0PidToken_ = consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("sigmat0PID"));
  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"));
  mtdgeoToken_ = esConsumes<MTDGeometry, MTDDigiGeometryRecord>();
  mtdtopoToken_ = esConsumes<MTDTopology, MTDTopologyRcd>();
  particleTableToken_ = esConsumes<HepPDT::ParticleDataTable, edm::DefaultRecord>();
}

~MtdTracksValidation()

MtdTracksValidation::~MtdTracksValidation ( )

override

Definition at line 228 of file MtdTracksValidation.cc.

{}

Member Function Documentation

analyze()

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

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 231 of file MtdTracksValidation.cc.

References funct::abs(), btlSimHitsToken_, ALCARECOTkAlJpsiMuMu_cff::charge, BTLDetId::column(), angle_units::operators::convertMmToCm(), geant_units::operators::convertUnitsTo(), MTDTopologyMode::crysLayoutFromTopoMode(), TauDecayModes::dec, deltaZcut_, depositBTLthreshold_, depositETLthreshold_, l1ctLayer1_cff::dZ, HCALHighEnergyHPDFilter_cfi::energy, etaMatchCut_, MTDTopologyMode::etlLayoutFromTopoMode(), etlSimHitsToken_, Exception, dqm::impl::MonitorElement::Fill(), GenRecTrackToken_, BTLDetId::geographicalId(), ETLDetId::geographicalId(), totem::nt2::vfat::geoId(), relativeConstraints::geom, edm::EventSetup::getHandle(), getMatchedTP(), MTDTopology::getMTDTopologyMode(), edm::EventSetup::getTransientHandle(), HepMCProductToken_, iEvent, LogDebug, edm::makeValid(), CaloTowersParam_cfi::mc, meBTLTrackEffEtaMtd_, meBTLTrackEffEtaTot_, meBTLTrackEffPhiMtd_, meBTLTrackEffPhiTot_, meBTLTrackEffPtMtd_, meBTLTrackEffPtTot_, meBTLTrackPtRes_, meBTLTrackRPTime_, meETLTrackEffEta2Mtd_, meETLTrackEffEtaMtd_, meETLTrackEffEtaTot_, meETLTrackEffPhi2Mtd_, meETLTrackEffPhiMtd_, meETLTrackEffPhiTot_, meETLTrackEffPt2Mtd_, meETLTrackEffPtMtd_, meETLTrackEffPtTot_, meETLTrackPtRes_, meETLTrackRPTime_, meMVATrackEffEtaTot_, meMVATrackEffPtTot_, meMVATrackMatchedEffEtaMtd_, meMVATrackMatchedEffEtaTot_, meMVATrackMatchedEffPtMtd_, meMVATrackMatchedEffPtTot_, meMVATrackPullTot_, meMVATrackResTot_, meMVATrackZposResTot_, meNTrackingParticles_, meTrackMatchedTPEffEtaEtl2Mtd_, meTrackMatchedTPEffEtaMtd_, meTrackMatchedTPEffEtaTot_, meTrackMatchedTPEffPtEtl2Mtd_, meTrackMatchedTPEffPtMtd_, meTrackMatchedTPEffPtTot_, meTrackMatchedTPmtdEffEtaMtd_, meTrackMatchedTPmtdEffEtaTot_, meTrackMatchedTPmtdEffPtMtd_, meTrackMatchedTPmtdEffPtTot_, meTrackMVAQual_, meTrackNumHits_, meTrackPathLenghtvsEta_, meTrackSigmat0Pid_, meTrackSigmat0SafePid_, meTrackSigmat0Src_, meTrackt0Pid_, meTrackt0SafePid_, meTrackt0Src_, meTracktmtd_, meUnassCrysEnergy_, meUnassDeposit_, meUnassLgadsEnergy_, meUnassociatedDetId_, mtdgeoToken_, mtdtopoToken_, mvaGenRecMatch(), mvaGenSel(), mvaRecSel(), mvaTPSel(), ETLDetId::nDisc(), RectangularMTDTopology::nrows(), optionalPlots_, LHEGenericFilter_cfi::ParticleID, particleTableToken_, pathLengthToken_, RectangularMTDTopology::pixelToModuleLocalPoint(), r2s_, DetId::rawId(), recoToSimAssociationToken_, RecTrackToken_, RecVertexToken_, BTLDetId::row(), s2r_, Sigmat0PidToken_, Sigmat0SafePidToken_, Sigmat0SrcToken_, SigmatmtdToken_, rpcPointValidation_cfi::simHit, simToRecoAssociationToken_, ProxyMTDTopology::specificTopology(), t0PidToken_, t0SafePidToken_, t0SrcToken_, tmtdToken_, GeomDet::toGlobal(), GeomDet::topology(), cmsswSequenceInfo::tp, ETLDetId::tp, tpWithMTD(), HLT_2023v12_cff::track, trackAssocToken_, trackingParticleCollectionToken_, trackMaxBtlEta_, trackMaxEtlEta_, trackMinEtlEta_, trackMinPt_, trackMVAQualToken_, uniqueId(), x, y, z, and MTDDetId::zside().

                                                                                     {
  using namespace edm;
  using namespace geant_units::operators;
  using namespace std;

  auto geometryHandle = iSetup.getTransientHandle(mtdgeoToken_);
  const MTDGeometry* geom = geometryHandle.product();
  auto topologyHandle = iSetup.getTransientHandle(mtdtopoToken_);
  const MTDTopology* topology = topologyHandle.product();

  bool topo1Dis = false;
  bool topo2Dis = false;
  if (MTDTopologyMode::etlLayoutFromTopoMode(topology->getMTDTopologyMode()) == ETLDetId::EtlLayout::tp) {
    topo1Dis = true;
  } else {
    topo2Dis = true;
  }

  auto GenRecTrackHandle = makeValid(iEvent.getHandle(GenRecTrackToken_));
  auto RecVertexHandle = makeValid(iEvent.getHandle(RecVertexToken_));

  std::unordered_map<uint32_t, MTDHit> m_btlHits;
  std::unordered_map<uint32_t, MTDHit> m_etlHits;
  std::unordered_map<uint32_t, std::set<unsigned long int>> m_btlTrkPerCell;
  std::unordered_map<uint32_t, std::set<unsigned long int>> m_etlTrkPerCell;

  const auto& tMtd = iEvent.get(tmtdToken_);
  const auto& SigmatMtd = iEvent.get(SigmatmtdToken_);
  const auto& t0Src = iEvent.get(t0SrcToken_);
  const auto& Sigmat0Src = iEvent.get(Sigmat0SrcToken_);
  const auto& t0Pid = iEvent.get(t0PidToken_);
  const auto& Sigmat0Pid = iEvent.get(Sigmat0PidToken_);
  const auto& t0Safe = iEvent.get(t0SafePidToken_);
  const auto& Sigmat0Safe = iEvent.get(Sigmat0SafePidToken_);
  const auto& mtdQualMVA = iEvent.get(trackMVAQualToken_);
  const auto& trackAssoc = iEvent.get(trackAssocToken_);
  const auto& pathLength = iEvent.get(pathLengthToken_);

  const auto& primRecoVtx = *(RecVertexHandle.product()->begin());

  // generator level information (HepMC format)
  auto GenEventHandle = makeValid(iEvent.getHandle(HepMCProductToken_));
  const HepMC::GenEvent* mc = GenEventHandle->GetEvent();
  double zsim = convertMmToCm((*(mc->vertices_begin()))->position().z());
  double tsim = (*(mc->vertices_begin()))->position().t() * CLHEP::mm / CLHEP::c_light;

  auto pdt = iSetup.getHandle(particleTableToken_);
  const HepPDT::ParticleDataTable* pdTable = pdt.product();

  auto simToRecoH = makeValid(iEvent.getHandle(simToRecoAssociationToken_));
  s2r_ = simToRecoH.product();

  auto recoToSimH = makeValid(iEvent.getHandle(recoToSimAssociationToken_));
  r2s_ = recoToSimH.product();

  // find all signal event trackId corresponding to an MTD simHit
  std::unordered_set<unsigned long int> mtdTrackId;

  std::unordered_set<unsigned long int> tpTrackId;
  auto tpHandle = makeValid(iEvent.getHandle(trackingParticleCollectionToken_));
  TrackingParticleCollection tpColl = *(tpHandle.product());
  for (const auto& tp : tpColl) {
    if (tp.eventId().bunchCrossing() == 0 && tp.eventId().event() == 0) {
      if (!mvaTPSel(tp))
        continue;
      if (optionalPlots_) {
        if (std::abs(tp.eta()) < trackMaxBtlEta_) {
          meNTrackingParticles_->Fill(0.5);
        } else if ((std::abs(tp.eta()) < trackMaxEtlEta_) && (std::abs(tp.eta()) > trackMinEtlEta_)) {
          meNTrackingParticles_->Fill(1.5);
        }
      }
      for (const auto& simTrk : tp.g4Tracks()) {
        auto const thisTId = uniqueId(simTrk.trackId(), simTrk.eventId());
        tpTrackId.insert(thisTId);
        LogDebug("MtdTracksValidation") << "TP simTrack id : " << thisTId;
      }
    }
  }

  //Fill maps with simhits accumulated per DetId

  auto btlSimHitsHandle = makeValid(iEvent.getHandle(btlSimHitsToken_));
  MixCollection<PSimHit> btlSimHits(btlSimHitsHandle.product());
  for (auto const& simHit : btlSimHits) {
    if (simHit.tof() < 0 || simHit.tof() > 25.)
      continue;
    DetId id = simHit.detUnitId();
    auto const thisHId = uniqueId(simHit.trackId(), simHit.eventId());
    m_btlTrkPerCell[id.rawId()].insert(thisHId);
    auto simHitIt = m_btlHits.emplace(id.rawId(), MTDHit()).first;
    // --- Accumulate the energy (in MeV) of SIM hits in the same detector cell
    (simHitIt->second).energy += convertUnitsTo(0.001_MeV, simHit.energyLoss());
  }

  uint32_t hcount(0);
  for (auto const& cell : m_btlTrkPerCell) {
    bool foundAssocTP = false;
    auto detId_key = cell.first;
    for (auto const& simtrack : cell.second) {
      if (tpTrackId.find(simtrack) != tpTrackId.end()) {
        foundAssocTP = true;
        mtdTrackId.insert(simtrack);
      }
    }
    if (foundAssocTP == false) {
      meUnassCrysEnergy_->Fill(log10(m_btlHits[detId_key].energy));
      if (m_btlHits[detId_key].energy > depositBTLthreshold_) {
        hcount++;
      }
    }
  }
  meUnassociatedDetId_->Fill(0.5, hcount);

  auto etlSimHitsHandle = makeValid(iEvent.getHandle(etlSimHitsToken_));
  MixCollection<PSimHit> etlSimHits(etlSimHitsHandle.product());
  for (auto const& simHit : etlSimHits) {
    if (simHit.tof() < 0 || simHit.tof() > 25.) {
      continue;
    }
    DetId id = simHit.detUnitId();
    auto const thisHId = uniqueId(simHit.trackId(), simHit.eventId());
    m_etlTrkPerCell[id.rawId()].insert(thisHId);
    auto simHitIt = m_etlHits.emplace(id.rawId(), MTDHit()).first;
    // --- Accumulate the energy (in MeV) of SIM hits in the same detector cell
    (simHitIt->second).energy += convertUnitsTo(0.001_MeV, simHit.energyLoss());
  }

  hcount = 0;
  for (auto const& cell : m_etlTrkPerCell) {
    bool foundAssocTP = false;
    auto detId_key = cell.first;
    for (auto const& simtrack : cell.second) {
      if (tpTrackId.find(simtrack) != tpTrackId.end()) {
        foundAssocTP = true;
        mtdTrackId.insert(simtrack);
      }
    }
    if (foundAssocTP == false) {
      meUnassLgadsEnergy_->Fill(log10(m_etlHits[detId_key].energy));
      if (m_etlHits[detId_key].energy > depositETLthreshold_) {
        hcount++;
      }
    }
  }
  meUnassociatedDetId_->Fill(1.5, hcount);

  // Search for TP without associated hits and unassociated DetIds above threshold
  if (optionalPlots_) {
    for (const auto& tp : tpColl) {
      if (tp.eventId().bunchCrossing() == 0 && tp.eventId().event() == 0) {
        if (!mvaTPSel(tp)) {
          continue;
        }

        // test BTL crystals for association

        if (std::abs(tp.eta()) < trackMaxBtlEta_) {
          bool tpIsAssoc = false;
          bool goodCell = false;
          for (auto const& cell : m_btlTrkPerCell) {
            auto detId_key = cell.first;
            if (m_btlHits[detId_key].energy < depositBTLthreshold_) {
              continue;
            }

            BTLDetId detId(detId_key);
            DetId geoId = detId.geographicalId(MTDTopologyMode::crysLayoutFromTopoMode(topology->getMTDTopologyMode()));
            const MTDGeomDet* thedet = geom->idToDet(geoId);
            if (thedet == nullptr)
              throw cms::Exception("MtdTracksValidation") << "GeographicalID: " << std::hex << geoId.rawId() << " ("
                                                          << detId.rawId() << ") is invalid!" << std::dec << std::endl;
            const ProxyMTDTopology& topoproxy = static_cast<const ProxyMTDTopology&>(thedet->topology());
            const RectangularMTDTopology& topo =
                static_cast<const RectangularMTDTopology&>(topoproxy.specificTopology());

            Local3DPoint local_point(convertMmToCm(m_btlHits[detId_key].x),
                                     convertMmToCm(m_btlHits[detId_key].y),
                                     convertMmToCm(m_btlHits[detId_key].z));

            local_point =
                topo.pixelToModuleLocalPoint(local_point, detId.row(topo.nrows()), detId.column(topo.nrows()));
            const auto& global_point = thedet->toGlobal(local_point);

            if (std::abs(tp.eta() - global_point.eta()) > etaMatchCut_) {
              continue;
            }
            goodCell = true;
            for (auto const& simtrack : cell.second) {
              for (auto const& TPsimtrack : tp.g4Tracks()) {
                auto const testId = uniqueId(TPsimtrack.trackId(), TPsimtrack.eventId());
                if (simtrack == testId) {
                  tpIsAssoc = true;
                  break;
                }
              }
            }
          }  //cell Loop
          if (!tpIsAssoc && goodCell) {
            meUnassDeposit_->Fill(0.5);
          }

        } else {
          // test ETL LGADs for association
          bool tpIsAssoc = false;
          bool goodCell = false;
          for (auto const& cell : m_etlTrkPerCell) {
            auto detId_key = cell.first;
            if (m_etlHits[detId_key].energy < depositETLthreshold_) {
              continue;
            }

            ETLDetId detId(detId_key);
            DetId geoId = detId.geographicalId();
            const MTDGeomDet* thedet = geom->idToDet(geoId);
            if (thedet == nullptr)
              throw cms::Exception("MtdTracksValidation") << "GeographicalID: " << std::hex << geoId.rawId() << " ("
                                                          << detId.rawId() << ") is invalid!" << std::dec << std::endl;

            Local3DPoint local_point(convertMmToCm(m_etlHits[detId_key].x),
                                     convertMmToCm(m_etlHits[detId_key].y),
                                     convertMmToCm(m_etlHits[detId_key].z));
            const auto& global_point = thedet->toGlobal(local_point);

            if (std::abs(tp.eta() - global_point.eta()) > etaMatchCut_) {
              continue;
            }
            goodCell = true;
            for (auto const& simtrack : cell.second) {
              for (auto const& TPsimtrack : tp.g4Tracks()) {
                auto const testId = uniqueId(TPsimtrack.trackId(), TPsimtrack.eventId());
                if (simtrack == testId) {
                  tpIsAssoc = true;
                  break;
                }
              }
            }
          }  //cell Loop
          if (!tpIsAssoc && goodCell) {
            meUnassDeposit_->Fill(1.5);
          }
        }  // tp BTL/ETL acceptance selection
      }
    }  //tp Loop
  }    //optionalPlots

  unsigned int index = 0;

  // flag to select events with reco vertex close to true simulated primary vertex, or PV fake (particle guns)
  const bool isGoodVtx = std::abs(primRecoVtx.z() - zsim) < deltaZcut_ || primRecoVtx.isFake();

  // --- Loop over all RECO tracks ---
  for (const auto& trackGen : *GenRecTrackHandle) {
    const reco::TrackRef trackref(iEvent.getHandle(GenRecTrackToken_), index);
    index++;

    if (trackAssoc[trackref] == -1) {
      LogInfo("mtdTracks") << "Extended track not associated";
      continue;
    }

    const reco::TrackRef mtdTrackref = reco::TrackRef(iEvent.getHandle(RecTrackToken_), trackAssoc[trackref]);
    const reco::Track& track = *mtdTrackref;

    bool isBTL = false;
    bool twoETLdiscs = false;

    if (track.pt() >= trackMinPt_ && std::abs(track.eta()) <= trackMaxEtlEta_) {
      meTracktmtd_->Fill(tMtd[trackref]);
      if (std::round(SigmatMtd[trackref] - Sigmat0Pid[trackref]) != 0) {
        LogWarning("mtdTracks")
            << "TimeError associated to refitted track is different from TimeError stored in tofPID "
               "sigmat0 ValueMap: this should not happen";
      }

      meTrackt0Src_->Fill(t0Src[trackref]);
      meTrackSigmat0Src_->Fill(Sigmat0Src[trackref]);

      meTrackt0Pid_->Fill(t0Pid[trackref]);
      meTrackSigmat0Pid_->Fill(Sigmat0Pid[trackref]);
      meTrackt0SafePid_->Fill(t0Safe[trackref]);
      meTrackSigmat0SafePid_->Fill(Sigmat0Safe[trackref]);
      meTrackMVAQual_->Fill(mtdQualMVA[trackref]);

      meTrackPathLenghtvsEta_->Fill(std::abs(track.eta()), pathLength[trackref]);

      if (std::abs(track.eta()) < trackMaxBtlEta_) {
        // --- all BTL tracks (with and without hit in MTD) ---
        meBTLTrackEffEtaTot_->Fill(track.eta());
        meBTLTrackEffPhiTot_->Fill(track.phi());
        meBTLTrackEffPtTot_->Fill(track.pt());

        bool MTDBtl = false;
        int numMTDBtlvalidhits = 0;
        for (const auto hit : track.recHits()) {
          if (hit->isValid() == false)
            continue;
          MTDDetId Hit = hit->geographicalId();
          if ((Hit.det() == 6) && (Hit.subdetId() == 1) && (Hit.mtdSubDetector() == 1)) {
            MTDBtl = true;
            numMTDBtlvalidhits++;
          }
        }
        meTrackNumHits_->Fill(numMTDBtlvalidhits);

        // --- keeping only tracks with last hit in MTD ---
        if (MTDBtl == true) {
          isBTL = true;
          meBTLTrackEffEtaMtd_->Fill(track.eta());
          meBTLTrackEffPhiMtd_->Fill(track.phi());
          meBTLTrackEffPtMtd_->Fill(track.pt());
          meBTLTrackRPTime_->Fill(track.t0());
          meBTLTrackPtRes_->Fill((trackGen.pt() - track.pt()) / trackGen.pt());
        }
      }  //loop over (geometrical) BTL tracks

      else {
        // --- all ETL tracks (with and without hit in MTD) ---
        if ((track.eta() < -trackMinEtlEta_) && (track.eta() > -trackMaxEtlEta_)) {
          meETLTrackEffEtaTot_[0]->Fill(track.eta());
          meETLTrackEffPhiTot_[0]->Fill(track.phi());
          meETLTrackEffPtTot_[0]->Fill(track.pt());
        }

        if ((track.eta() > trackMinEtlEta_) && (track.eta() < trackMaxEtlEta_)) {
          meETLTrackEffEtaTot_[1]->Fill(track.eta());
          meETLTrackEffPhiTot_[1]->Fill(track.phi());
          meETLTrackEffPtTot_[1]->Fill(track.pt());
        }

        bool MTDEtlZnegD1 = false;
        bool MTDEtlZnegD2 = false;
        bool MTDEtlZposD1 = false;
        bool MTDEtlZposD2 = false;
        int numMTDEtlvalidhits = 0;
        for (const auto hit : track.recHits()) {
          if (hit->isValid() == false)
            continue;
          MTDDetId Hit = hit->geographicalId();
          if ((Hit.det() == 6) && (Hit.subdetId() == 1) && (Hit.mtdSubDetector() == 2)) {
            ETLDetId ETLHit = hit->geographicalId();

            if (topo2Dis) {
              if ((ETLHit.zside() == -1) && (ETLHit.nDisc() == 1)) {
                MTDEtlZnegD1 = true;
                meETLTrackRPTime_->Fill(track.t0());
                meETLTrackPtRes_->Fill((trackGen.pt() - track.pt()) / trackGen.pt());
                numMTDEtlvalidhits++;
              }
              if ((ETLHit.zside() == -1) && (ETLHit.nDisc() == 2)) {
                MTDEtlZnegD2 = true;
                meETLTrackRPTime_->Fill(track.t0());
                meETLTrackPtRes_->Fill((trackGen.pt() - track.pt()) / trackGen.pt());
                numMTDEtlvalidhits++;
              }
              if ((ETLHit.zside() == 1) && (ETLHit.nDisc() == 1)) {
                MTDEtlZposD1 = true;
                meETLTrackRPTime_->Fill(track.t0());
                meETLTrackPtRes_->Fill((trackGen.pt() - track.pt()) / trackGen.pt());
                numMTDEtlvalidhits++;
              }
              if ((ETLHit.zside() == 1) && (ETLHit.nDisc() == 2)) {
                MTDEtlZposD2 = true;
                meETLTrackRPTime_->Fill(track.t0());
                meETLTrackPtRes_->Fill((trackGen.pt() - track.pt()) / trackGen.pt());
                numMTDEtlvalidhits++;
              }
            }

            if (topo1Dis) {
              if (ETLHit.zside() == -1) {
                MTDEtlZnegD1 = true;
                meETLTrackRPTime_->Fill(track.t0());
                numMTDEtlvalidhits++;
              }
              if (ETLHit.zside() == 1) {
                MTDEtlZposD1 = true;
                meETLTrackRPTime_->Fill(track.t0());
                numMTDEtlvalidhits++;
              }
            }
          }
        }
        meTrackNumHits_->Fill(-numMTDEtlvalidhits);

        // --- keeping only tracks with last hit in MTD ---
        if ((track.eta() < -trackMinEtlEta_) && (track.eta() > -trackMaxEtlEta_)) {
          twoETLdiscs = (MTDEtlZnegD1 == true) && (MTDEtlZnegD2 == true);
          if ((MTDEtlZnegD1 == true) || (MTDEtlZnegD2 == true)) {
            meETLTrackEffEtaMtd_[0]->Fill(track.eta());
            meETLTrackEffPhiMtd_[0]->Fill(track.phi());
            meETLTrackEffPtMtd_[0]->Fill(track.pt());
            if (twoETLdiscs) {
              meETLTrackEffEta2Mtd_[0]->Fill(track.eta());
              meETLTrackEffPhi2Mtd_[0]->Fill(track.phi());
              meETLTrackEffPt2Mtd_[0]->Fill(track.pt());
            }
          }
        }
        if ((track.eta() > trackMinEtlEta_) && (track.eta() < trackMaxEtlEta_)) {
          twoETLdiscs = (MTDEtlZposD1 == true) && (MTDEtlZposD2 == true);
          if ((MTDEtlZposD1 == true) || (MTDEtlZposD2 == true)) {
            meETLTrackEffEtaMtd_[1]->Fill(track.eta());
            meETLTrackEffPhiMtd_[1]->Fill(track.phi());
            meETLTrackEffPtMtd_[1]->Fill(track.pt());
            if (twoETLdiscs) {
              meETLTrackEffEta2Mtd_[1]->Fill(track.eta());
              meETLTrackEffPhi2Mtd_[1]->Fill(track.phi());
              meETLTrackEffPt2Mtd_[1]->Fill(track.pt());
            }
          }
        }
      }
    }

    if (isGoodVtx) {
      bool noCrack = std::abs(trackGen.eta()) < trackMaxBtlEta_ || std::abs(trackGen.eta()) > trackMinEtlEta_;
      const bool vtxFake = primRecoVtx.isFake();

      if (mvaRecSel(trackGen, primRecoVtx, t0Safe[trackref], Sigmat0Safe[trackref])) {
        // reco-gen matching used for MVA quality flag

        if (noCrack) {
          meMVATrackEffPtTot_->Fill(trackGen.pt());
        }
        meMVATrackEffEtaTot_->Fill(std::abs(trackGen.eta()));

        double dZ = trackGen.vz() - zsim;
        double dT(-9999.);
        double pullT(-9999.);
        if (Sigmat0Safe[trackref] != -1.) {
          dT = t0Safe[trackref] - tsim;
          pullT = dT / Sigmat0Safe[trackref];
        }
        for (const auto& genP : mc->particle_range()) {
          // select status 1 genParticles and match them to the reconstructed track

          float charge = pdTable->particle(HepPDT::ParticleID(genP->pdg_id())) != nullptr
                             ? pdTable->particle(HepPDT::ParticleID(genP->pdg_id()))->charge()
                             : 0.f;
          if (mvaGenSel(*genP, charge)) {
            if (mvaGenRecMatch(*genP, zsim, trackGen, vtxFake)) {
              meMVATrackZposResTot_->Fill(dZ);
              if (noCrack) {
                meMVATrackMatchedEffPtTot_->Fill(trackGen.pt());
              }
              meMVATrackMatchedEffEtaTot_->Fill(std::abs(trackGen.eta()));
              if (pullT > -9999.) {
                meMVATrackResTot_->Fill(dT);
                meMVATrackPullTot_->Fill(pullT);
                if (noCrack) {
                  meMVATrackMatchedEffPtMtd_->Fill(trackGen.pt());
                }
                meMVATrackMatchedEffEtaMtd_->Fill(std::abs(trackGen.eta()));
              }
              break;
            }
          }
        }

        // TrackingParticle based matching

        const reco::TrackBaseRef trkrefb(trackref);
        auto tp_info = getMatchedTP(trkrefb, zsim);

        if (tp_info != nullptr && mvaTPSel(**tp_info)) {
          const bool withMTD = tpWithMTD(**tp_info, mtdTrackId);
          if (noCrack) {
            meTrackMatchedTPEffPtTot_->Fill(trackGen.pt());
            if (withMTD) {
              meTrackMatchedTPmtdEffPtTot_->Fill(trackGen.pt());
            }
          }
          meTrackMatchedTPEffEtaTot_->Fill(std::abs(trackGen.eta()));
          if (withMTD) {
            meTrackMatchedTPmtdEffEtaTot_->Fill(std::abs(trackGen.eta()));
          }
          if (pullT > -9999.) {
            if (noCrack) {
              meTrackMatchedTPEffPtMtd_->Fill(trackGen.pt());
              if (isBTL || twoETLdiscs) {
                meTrackMatchedTPEffPtEtl2Mtd_->Fill(trackGen.pt());
              }
              if (withMTD) {
                meTrackMatchedTPmtdEffPtMtd_->Fill(trackGen.pt());
              }
            }
            meTrackMatchedTPEffEtaMtd_->Fill(std::abs(trackGen.eta()));
            if (isBTL || twoETLdiscs) {
              meTrackMatchedTPEffEtaEtl2Mtd_->Fill(std::abs(trackGen.eta()));
            }
            if (withMTD) {
              meTrackMatchedTPmtdEffEtaMtd_->Fill(std::abs(trackGen.eta()));
            }
          }
        }
      }
    }  // MC truth matich analysis for good PV
  }    //RECO tracks loop
}

bookHistograms()

void MtdTracksValidation::bookHistograms ( DQMStore::IBooker &  ibook, edm::Run const &  run, edm::EventSetup const &  iSetup  )

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 733 of file MtdTracksValidation.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::bookProfile(), folder_, meBTLTrackEffEtaMtd_, meBTLTrackEffEtaTot_, meBTLTrackEffPhiMtd_, meBTLTrackEffPhiTot_, meBTLTrackEffPtMtd_, meBTLTrackEffPtTot_, meBTLTrackPtRes_, meBTLTrackRPTime_, meETLTrackEffEta2Mtd_, meETLTrackEffEtaMtd_, meETLTrackEffEtaTot_, meETLTrackEffPhi2Mtd_, meETLTrackEffPhiMtd_, meETLTrackEffPhiTot_, meETLTrackEffPt2Mtd_, meETLTrackEffPtMtd_, meETLTrackEffPtTot_, meETLTrackPtRes_, meETLTrackRPTime_, meMVATrackEffEtaTot_, meMVATrackEffPtTot_, meMVATrackMatchedEffEtaMtd_, meMVATrackMatchedEffEtaTot_, meMVATrackMatchedEffPtMtd_, meMVATrackMatchedEffPtTot_, meMVATrackPullTot_, meMVATrackResTot_, meMVATrackZposResTot_, meNTrackingParticles_, meTrackMatchedTPEffEtaEtl2Mtd_, meTrackMatchedTPEffEtaMtd_, meTrackMatchedTPEffEtaTot_, meTrackMatchedTPEffPtEtl2Mtd_, meTrackMatchedTPEffPtMtd_, meTrackMatchedTPEffPtTot_, meTrackMatchedTPmtdEffEtaMtd_, meTrackMatchedTPmtdEffEtaTot_, meTrackMatchedTPmtdEffPtMtd_, meTrackMatchedTPmtdEffPtTot_, meTrackMVAQual_, meTrackNumHits_, meTrackPathLenghtvsEta_, meTrackSigmat0Pid_, meTrackSigmat0SafePid_, meTrackSigmat0Src_, meTrackt0Pid_, meTrackt0SafePid_, meTrackt0Src_, meTracktmtd_, meUnassCrysEnergy_, meUnassDeposit_, meUnassLgadsEnergy_, meUnassociatedDetId_, and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                                                                               {
  ibook.setCurrentFolder(folder_);

  // histogram booking
  meBTLTrackRPTime_ = ibook.book1D("TrackBTLRPTime", "Track t0 with respect to R.P.;t0 [ns]", 100, -1, 3);
  meBTLTrackEffEtaTot_ = ibook.book1D("TrackBTLEffEtaTot", "Track efficiency vs eta (Tot);#eta_{RECO}", 100, -1.6, 1.6);
  meBTLTrackEffPhiTot_ =
      ibook.book1D("TrackBTLEffPhiTot", "Track efficiency vs phi (Tot);#phi_{RECO} [rad]", 100, -3.2, 3.2);
  meBTLTrackEffPtTot_ = ibook.book1D("TrackBTLEffPtTot", "Track efficiency vs pt (Tot);pt_{RECO} [GeV]", 50, 0, 10);
  meBTLTrackEffEtaMtd_ = ibook.book1D("TrackBTLEffEtaMtd", "Track efficiency vs eta (Mtd);#eta_{RECO}", 100, -1.6, 1.6);
  meBTLTrackEffPhiMtd_ =
      ibook.book1D("TrackBTLEffPhiMtd", "Track efficiency vs phi (Mtd);#phi_{RECO} [rad]", 100, -3.2, 3.2);
  meBTLTrackEffPtMtd_ = ibook.book1D("TrackBTLEffPtMtd", "Track efficiency vs pt (Mtd);pt_{RECO} [GeV]", 50, 0, 10);
  meBTLTrackPtRes_ =
      ibook.book1D("TrackBTLPtRes", "Track pT resolution  ;pT_{Gentrack}-pT_{MTDtrack}/pT_{Gentrack} ", 100, -0.1, 0.1);
  meETLTrackRPTime_ = ibook.book1D("TrackETLRPTime", "Track t0 with respect to R.P.;t0 [ns]", 100, -1, 3);
  meETLTrackEffEtaTot_[0] =
      ibook.book1D("TrackETLEffEtaTotZneg", "Track efficiency vs eta (Tot) (-Z);#eta_{RECO}", 100, -3.2, -1.4);
  meETLTrackEffEtaTot_[1] =
      ibook.book1D("TrackETLEffEtaTotZpos", "Track efficiency vs eta (Tot) (+Z);#eta_{RECO}", 100, 1.4, 3.2);
  meETLTrackEffPhiTot_[0] =
      ibook.book1D("TrackETLEffPhiTotZneg", "Track efficiency vs phi (Tot) (-Z);#phi_{RECO} [rad]", 100, -3.2, 3.2);
  meETLTrackEffPhiTot_[1] =
      ibook.book1D("TrackETLEffPhiTotZpos", "Track efficiency vs phi (Tot) (+Z);#phi_{RECO} [rad]", 100, -3.2, 3.2);
  meETLTrackEffPtTot_[0] =
      ibook.book1D("TrackETLEffPtTotZneg", "Track efficiency vs pt (Tot) (-Z);pt_{RECO} [GeV]", 50, 0, 10);
  meETLTrackEffPtTot_[1] =
      ibook.book1D("TrackETLEffPtTotZpos", "Track efficiency vs pt (Tot) (+Z);pt_{RECO} [GeV]", 50, 0, 10);
  meETLTrackEffEtaMtd_[0] =
      ibook.book1D("TrackETLEffEtaMtdZneg", "Track efficiency vs eta (Mtd) (-Z);#eta_{RECO}", 100, -3.2, -1.4);
  meETLTrackEffEtaMtd_[1] =
      ibook.book1D("TrackETLEffEtaMtdZpos", "Track efficiency vs eta (Mtd) (+Z);#eta_{RECO}", 100, 1.4, 3.2);
  meETLTrackEffPhiMtd_[0] =
      ibook.book1D("TrackETLEffPhiMtdZneg", "Track efficiency vs phi (Mtd) (-Z);#phi_{RECO} [rad]", 100, -3.2, 3.2);
  meETLTrackEffPhiMtd_[1] =
      ibook.book1D("TrackETLEffPhiMtdZpos", "Track efficiency vs phi (Mtd) (+Z);#phi_{RECO} [rad]", 100, -3.2, 3.2);
  meETLTrackEffPtMtd_[0] =
      ibook.book1D("TrackETLEffPtMtdZneg", "Track efficiency vs pt (Mtd) (-Z);pt_{RECO} [GeV]", 50, 0, 10);
  meETLTrackEffPtMtd_[1] =
      ibook.book1D("TrackETLEffPtMtdZpos", "Track efficiency vs pt (Mtd) (+Z);pt_{RECO} [GeV]", 50, 0, 10);
  meETLTrackEffEta2Mtd_[0] =
      ibook.book1D("TrackETLEffEta2MtdZneg", "Track efficiency vs eta (Mtd 2 hit) (-Z);#eta_{RECO}", 100, -3.2, -1.4);
  meETLTrackEffEta2Mtd_[1] =
      ibook.book1D("TrackETLEffEta2MtdZpos", "Track efficiency vs eta (Mtd 2 hit) (+Z);#eta_{RECO}", 100, 1.4, 3.2);
  meETLTrackEffPhi2Mtd_[0] = ibook.book1D(
      "TrackETLEffPhi2MtdZneg", "Track efficiency vs phi (Mtd 2 hit) (-Z);#phi_{RECO} [rad]", 100, -3.2, 3.2);
  meETLTrackEffPhi2Mtd_[1] = ibook.book1D(
      "TrackETLEffPhi2MtdZpos", "Track efficiency vs phi (Mtd 2 hit) (+Z);#phi_{RECO} [rad]", 100, -3.2, 3.2);
  meETLTrackEffPt2Mtd_[0] =
      ibook.book1D("TrackETLEffPt2MtdZneg", "Track efficiency vs pt (Mtd 2 hit) (-Z);pt_{RECO} [GeV]", 50, 0, 10);
  meETLTrackEffPt2Mtd_[1] =
      ibook.book1D("TrackETLEffPt2MtdZpos", "Track efficiency vs pt (Mtd 2 hit) (+Z);pt_{RECO} [GeV]", 50, 0, 10);
  meETLTrackPtRes_ =
      ibook.book1D("TrackETLPtRes", "Track pT resolution;pT_{Gentrack}-pT_{MTDtrack}/pT_{Gentrack} ", 100, -0.1, 0.1);

  meTracktmtd_ = ibook.book1D("Tracktmtd", "Track time from TrackExtenderWithMTD;tmtd [ns]", 150, 1, 16);
  meTrackt0Src_ = ibook.book1D("Trackt0Src", "Track time from TrackExtenderWithMTD;t0Src [ns]", 100, -1.5, 1.5);
  meTrackSigmat0Src_ =
      ibook.book1D("TrackSigmat0Src", "Time Error from TrackExtenderWithMTD; #sigma_{t0Src} [ns]", 100, 0, 0.1);

  meTrackt0Pid_ = ibook.book1D("Trackt0Pid", "Track t0 as stored in TofPid;t0 [ns]", 100, -1, 1);
  meTrackSigmat0Pid_ = ibook.book1D("TrackSigmat0Pid", "Sigmat0 as stored in TofPid; #sigma_{t0} [ns]", 100, 0, 0.1);
  meTrackt0SafePid_ = ibook.book1D("Trackt0SafePID", "Track t0 Safe as stored in TofPid;t0 [ns]", 100, -1, 1);
  meTrackSigmat0SafePid_ =
      ibook.book1D("TrackSigmat0SafePID", "Sigmat0 Safe as stored in TofPid; #sigma_{t0} [ns]", 100, 0, 0.1);
  meTrackNumHits_ = ibook.book1D("TrackNumHits", "Number of valid MTD hits per track ; Number of hits", 10, -5, 5);
  meTrackMVAQual_ = ibook.book1D("TrackMVAQual", "Track MVA Quality as stored in Value Map ; MVAQual", 100, 0, 1);
  meTrackPathLenghtvsEta_ = ibook.bookProfile(
      "TrackPathLenghtvsEta", "MTD Track pathlength vs MTD track Eta;|#eta|;Pathlength", 100, 0, 3.2, 100.0, 400.0, "S");

  meMVATrackEffPtTot_ = ibook.book1D("MVAEffPtTot", "Pt of tracks associated to LV; track pt [GeV] ", 110, 0., 11.);
  meMVATrackMatchedEffPtTot_ =
      ibook.book1D("MVAMatchedEffPtTot", "Pt of tracks associated to LV matched to GEN; track pt [GeV] ", 110, 0., 11.);
  meMVATrackMatchedEffPtMtd_ = ibook.book1D(
      "MVAMatchedEffPtMtd", "Pt of tracks associated to LV matched to GEN with time; track pt [GeV] ", 110, 0., 11.);

  meTrackMatchedTPEffPtTot_ =
      ibook.book1D("MatchedTPEffPtTot", "Pt of tracks associated to LV matched to TP; track pt [GeV] ", 110, 0., 11.);
  meTrackMatchedTPEffPtMtd_ = ibook.book1D(
      "MatchedTPEffPtMtd", "Pt of tracks associated to LV matched to TP with time; track pt [GeV] ", 110, 0., 11.);
  meTrackMatchedTPEffPtEtl2Mtd_ =
      ibook.book1D("MatchedTPEffPtEtl2Mtd",
                   "Pt of tracks associated to LV matched to TP with time, 2 ETL hits; track pt [GeV] ",
                   110,
                   0.,
                   11.);

  meTrackMatchedTPmtdEffPtTot_ = ibook.book1D(
      "MatchedTPmtdEffPtTot", "Pt of tracks associated to LV matched to TP-mtd hit; track pt [GeV] ", 110, 0., 11.);
  meTrackMatchedTPmtdEffPtMtd_ =
      ibook.book1D("MatchedTPmtdEffPtMtd",
                   "Pt of tracks associated to LV matched to TP-mtd hit with time; track pt [GeV] ",
                   110,
                   0.,
                   11.);

  meMVATrackEffEtaTot_ = ibook.book1D("MVAEffEtaTot", "Eta of tracks associated to LV; track eta ", 66, 0., 3.3);
  meMVATrackMatchedEffEtaTot_ =
      ibook.book1D("MVAMatchedEffEtaTot", "Eta of tracks associated to LV matched to GEN; track eta ", 66, 0., 3.3);
  meMVATrackMatchedEffEtaMtd_ = ibook.book1D(
      "MVAMatchedEffEtaMtd", "Eta of tracks associated to LV matched to GEN with time; track eta ", 66, 0., 3.3);

  meTrackMatchedTPEffEtaTot_ =
      ibook.book1D("MatchedTPEffEtaTot", "Eta of tracks associated to LV matched to TP; track eta ", 66, 0., 3.3);
  meTrackMatchedTPEffEtaMtd_ = ibook.book1D(
      "MatchedTPEffEtaMtd", "Eta of tracks associated to LV matched to TP with time; track eta ", 66, 0., 3.3);
  meTrackMatchedTPEffEtaEtl2Mtd_ =
      ibook.book1D("MatchedTPEffEtaEtl2Mtd",
                   "Eta of tracks associated to LV matched to TP with time, 2 ETL hits; track eta ",
                   66,
                   0.,
                   3.3);

  meTrackMatchedTPmtdEffEtaTot_ = ibook.book1D(
      "MatchedTPmtdEffEtaTot", "Eta of tracks associated to LV matched to TP-mtd hit; track eta ", 66, 0., 3.3);
  meTrackMatchedTPmtdEffEtaMtd_ =
      ibook.book1D("MatchedTPmtdEffEtaMtd",
                   "Eta of tracks associated to LV matched to TP-mtd hit with time; track eta ",
                   66,
                   0.,
                   3.3);

  meMVATrackResTot_ = ibook.book1D(
      "MVATrackRes", "t_{rec} - t_{sim} for LV associated tracks; t_{rec} - t_{sim} [ns] ", 120, -0.15, 0.15);
  meMVATrackPullTot_ =
      ibook.book1D("MVATrackPull", "Pull for associated tracks; (t_{rec}-t_{sim})/#sigma_{t}", 50, -5., 5.);
  meMVATrackZposResTot_ = ibook.book1D(
      "MVATrackZposResTot", "Z_{PCA} - Z_{sim} for associated tracks;Z_{PCA} - Z_{sim} [cm] ", 100, -0.1, 0.1);

  meUnassociatedDetId_ = ibook.bookProfile(
      "UnassociatedDetId", "Number of MTD cell not associated to any TP per event", 2, 0., 2., 0., 100000., "S");
  meNTrackingParticles_ = ibook.book1D("NTrackingParticles", "Total #Tracking particles", 2, 0, 2);
  meUnassDeposit_ =
      ibook.book1D("UnassDeposit",
                   "#Tracking particles with deposit over threshold in MTD cell, but with no cell associated to TP;",
                   2,
                   0,
                   2);
  meUnassCrysEnergy_ =
      ibook.book1D("UnassCrysEnergy",
                   "Energy deposit in BTL crystal with no associated SimTracks;log_{10}(Energy [MeV]) ",
                   100,
                   -3.5,
                   1.5);
  meUnassLgadsEnergy_ = ibook.book1D("UnassLgadsEnergy",
                                     "Energy deposit in ETL LGADs with no associated SimTracks;log_{10}(Energy [MeV]) ",
                                     100,
                                     -3.5,
                                     1.5);
}

fillDescriptions()

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

static

Definition at line 886 of file MtdTracksValidation.cc.

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

                                                                                     {
  edm::ParameterSetDescription desc;

  desc.add<std::string>("folder", "MTD/Tracks");
  desc.add<edm::InputTag>("inputTagG", edm::InputTag("generalTracks"));
  desc.add<edm::InputTag>("inputTagT", edm::InputTag("trackExtenderWithMTD"));
  desc.add<edm::InputTag>("inputTagV", edm::InputTag("offlinePrimaryVertices4D"));
  desc.add<edm::InputTag>("inputTagH", edm::InputTag("generatorSmeared"));
  desc.add<edm::InputTag>("SimTag", edm::InputTag("mix", "MergedTrackTruth"));
  desc.add<edm::InputTag>("TPtoRecoTrackAssoc", edm::InputTag("trackingParticleRecoTrackAsssociation"));
  desc.add<edm::InputTag>("btlSimHits", edm::InputTag("mix", "g4SimHitsFastTimerHitsBarrel"));
  desc.add<edm::InputTag>("etlSimHits", edm::InputTag("mix", "g4SimHitsFastTimerHitsEndcap"));
  desc.add<edm::InputTag>("tmtd", edm::InputTag("trackExtenderWithMTD:generalTracktmtd"));
  desc.add<edm::InputTag>("sigmatmtd", edm::InputTag("trackExtenderWithMTD:generalTracksigmatmtd"));
  desc.add<edm::InputTag>("t0Src", edm::InputTag("trackExtenderWithMTD:generalTrackt0"));
  desc.add<edm::InputTag>("sigmat0Src", edm::InputTag("trackExtenderWithMTD:generalTracksigmat0"));
  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>("t0SafePID", edm::InputTag("tofPID:t0safe"));
  desc.add<edm::InputTag>("sigmat0SafePID", edm::InputTag("tofPID:sigmat0safe"));
  desc.add<edm::InputTag>("sigmat0PID", edm::InputTag("tofPID:sigmat0"));
  desc.add<edm::InputTag>("t0PID", edm::InputTag("tofPID:t0"));
  desc.add<edm::InputTag>("trackMVAQual", edm::InputTag("mtdTrackQualityMVA:mtdQualMVA"));
  desc.add<double>("trackMinimumPt", 0.7);  // [GeV]
  desc.add<double>("trackMaximumBtlEta", 1.5);
  desc.add<double>("trackMinimumEtlEta", 1.6);
  desc.add<double>("trackMaximumEtlEta", 3.);
  desc.addUntracked<bool>("optionalPlots", true);

  descriptions.add("mtdTracksValid", desc);
}

getMatchedTP()

const edm::Ref< std::vector< TrackingParticle > > * MtdTracksValidation::getMatchedTP ( const reco::TrackBaseRef &  recoTrack, const double &  zsim  )

private

Definition at line 968 of file MtdTracksValidation.cc.

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

Referenced by analyze().

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

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

  //matched TP equal to any TP associated to signal sim vertex
  for (const auto& tp : found->val) {
    if (tp.first->eventId().bunchCrossing() == 0 && tp.first->eventId().event() == 0 &&
        std::abs(tp.first->parentVertex()->position().z() - zsim) < deltaZcut_) {
      return &tp.first;
    }
  }

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

isETL()

bool MtdTracksValidation::isETL ( const double  eta ) const

inlineprivate

Definition at line 78 of file MtdTracksValidation.cc.

References funct::abs(), PVValHelper::eta, trackMaxEtlEta_, and trackMinEtlEta_.

{ return (std::abs(eta) > trackMinEtlEta_) && (std::abs(eta) < trackMaxEtlEta_); }

mvaGenRecMatch()

const bool MtdTracksValidation::mvaGenRecMatch ( const HepMC::GenParticle &  genP, const double &  zsim, const reco::TrackBase &  trk, const bool &  vtxFake  )

private

Definition at line 956 of file MtdTracksValidation.cc.

References funct::abs(), deltaDRcut_, deltaPTcut_, reco::deltaR(), deltaZcut_, HGC3DClusterGenMatchSelector_cfi::dR, match(), reco::TrackBase::momentum(), reco::TrackBase::pt(), and reco::TrackBase::vz().

Referenced by analyze().

                                                                    {
  bool match = false;
  double dR = reco::deltaR(genP.momentum(), trk.momentum());
  double genPT = genP.momentum().perp();
  match = std::abs(genPT - trk.pt()) < trk.pt() * deltaPTcut_ && dR < deltaDRcut_ &&
          (std::abs(trk.vz() - zsim) < deltaZcut_ || vtxFake);
  return match;
}

mvaGenSel()

const bool MtdTracksValidation::mvaGenSel ( const HepMC::GenParticle &  gp, const float &  charge  )

private

Definition at line 919 of file MtdTracksValidation.cc.

References funct::abs(), ALCARECOTkAlJpsiMuMu_cff::charge, etacutGEN_, runTauDisplay::gp, match(), and pTcut_.

Referenced by analyze().

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

mvaRecSel()

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

private

Definition at line 943 of file MtdTracksValidation.cc.

References funct::abs(), deltaZcut_, reco::TrackBase::eta(), etacutREC_, match(), reco::TrackBase::pt(), pTcut_, FrontierCondition_GT_autoExpress_cfi::t0, extraflags_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_ || vtx.isFake());
  if (st0 > 0.) {
    match = match && std::abs(t0 - vtx.t()) < 3. * st0;
  }
  return match;
}

mvaTPSel()

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

private

Definition at line 928 of file MtdTracksValidation.cc.

References funct::abs(), etacutGEN_, match(), pTcut_, rBTL_, mathSSE::sqrt(), cmsswSequenceInfo::tp, and zETL_.

Referenced by analyze().

                                                                   {
  bool match = false;
  if (tp.status() != 1) {
    return match;
  }
  auto x_pv = tp.parentVertex()->position().x();
  auto y_pv = tp.parentVertex()->position().y();
  auto z_pv = tp.parentVertex()->position().z();

  auto r_pv = std::sqrt(x_pv * x_pv + y_pv * y_pv);

  match = tp.charge() != 0 && tp.pt() > pTcut_ && std::abs(tp.eta()) < etacutGEN_ && r_pv < rBTL_ && z_pv < zETL_;
  return match;
}

tpWithMTD()

const bool MtdTracksValidation::tpWithMTD ( const TrackingParticle &  tp, const std::unordered_set< unsigned long int > &  trkList  )

private

Definition at line 988 of file MtdTracksValidation.cc.

References cmsswSequenceInfo::tp, and uniqueId().

Referenced by analyze().

                                                                                              {
  for (const auto& simTrk : tp.g4Tracks()) {
    for (const auto& mtdTrk : trkList) {
      if (uniqueId(simTrk.trackId(), simTrk.eventId()) == mtdTrk) {
        return true;
      }
    }
  }
  return false;
}

uniqueId()

const unsigned long int MtdTracksValidation::uniqueId ( const uint32_t  x, const EncodedEventId &  y  )

inlineprivate

Definition at line 68 of file MtdTracksValidation.cc.

References a, b, x, and y.

Referenced by analyze(), and tpWithMTD().

                                                                              {
    const uint64_t a = static_cast<uint64_t>(x);
    const uint64_t b = static_cast<uint64_t>(y.rawId());

    if (x < y.rawId())
      return (b << 32) | a;
    else
      return (a << 32) | b;
  }

Member Data Documentation

btlSimHitsToken_

edm::EDGetTokenT<CrossingFrame<PSimHit> > MtdTracksValidation::btlSimHitsToken_

private

Definition at line 113 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

deltaDRcut_

constexpr double MtdTracksValidation::deltaDRcut_ = 0.03

staticprivate

Definition at line 93 of file MtdTracksValidation.cc.

Referenced by mvaGenRecMatch().

deltaPTcut_

constexpr double MtdTracksValidation::deltaPTcut_ = 0.05

staticprivate

Definition at line 92 of file MtdTracksValidation.cc.

Referenced by mvaGenRecMatch().

deltaZcut_

constexpr double MtdTracksValidation::deltaZcut_ = 0.1

staticprivate

Definition at line 91 of file MtdTracksValidation.cc.

Referenced by analyze(), getMatchedTP(), mvaGenRecMatch(), and mvaRecSel().

depositBTLthreshold_

constexpr double MtdTracksValidation::depositBTLthreshold_ = 1

staticprivate

Definition at line 94 of file MtdTracksValidation.cc.

Referenced by analyze().

depositETLthreshold_

constexpr double MtdTracksValidation::depositETLthreshold_ = 0.001

staticprivate

Definition at line 95 of file MtdTracksValidation.cc.

Referenced by analyze().

etacutGEN_

constexpr double MtdTracksValidation::etacutGEN_ = 4.

staticprivate

Definition at line 88 of file MtdTracksValidation.cc.

Referenced by mvaGenSel(), and mvaTPSel().

etacutREC_

constexpr double MtdTracksValidation::etacutREC_ = 3.

staticprivate

Definition at line 89 of file MtdTracksValidation.cc.

Referenced by mvaRecSel().

etaMatchCut_

constexpr double MtdTracksValidation::etaMatchCut_ = 0.05

staticprivate

Definition at line 98 of file MtdTracksValidation.cc.

Referenced by analyze().

etlSimHitsToken_

edm::EDGetTokenT<CrossingFrame<PSimHit> > MtdTracksValidation::etlSimHitsToken_

private

Definition at line 114 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

folder_

const std::string MtdTracksValidation::folder_

private

Definition at line 82 of file MtdTracksValidation.cc.

Referenced by bookHistograms().

GenRecTrackToken_

edm::EDGetTokenT<reco::TrackCollection> MtdTracksValidation::GenRecTrackToken_

private

Definition at line 105 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

HepMCProductToken_

edm::EDGetTokenT<edm::HepMCProduct> MtdTracksValidation::HepMCProductToken_

private

Definition at line 109 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

meBTLTrackEffEtaMtd_

MonitorElement* MtdTracksValidation::meBTLTrackEffEtaMtd_

private

Definition at line 137 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meBTLTrackEffEtaTot_

MonitorElement* MtdTracksValidation::meBTLTrackEffEtaTot_

private

Definition at line 134 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meBTLTrackEffPhiMtd_

MonitorElement* MtdTracksValidation::meBTLTrackEffPhiMtd_

private

Definition at line 138 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meBTLTrackEffPhiTot_

MonitorElement* MtdTracksValidation::meBTLTrackEffPhiTot_

private

Definition at line 135 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meBTLTrackEffPtMtd_

MonitorElement* MtdTracksValidation::meBTLTrackEffPtMtd_

private

Definition at line 139 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meBTLTrackEffPtTot_

MonitorElement* MtdTracksValidation::meBTLTrackEffPtTot_

private

Definition at line 136 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meBTLTrackPtRes_

MonitorElement* MtdTracksValidation::meBTLTrackPtRes_

private

Definition at line 140 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meBTLTrackRPTime_

MonitorElement* MtdTracksValidation::meBTLTrackRPTime_

private

Definition at line 133 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meETLTrackEffEta2Mtd_

MonitorElement* MtdTracksValidation::meETLTrackEffEta2Mtd_[2]

private

Definition at line 149 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meETLTrackEffEtaMtd_

MonitorElement* MtdTracksValidation::meETLTrackEffEtaMtd_[2]

private

Definition at line 146 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meETLTrackEffEtaTot_

MonitorElement* MtdTracksValidation::meETLTrackEffEtaTot_[2]

private

Definition at line 143 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meETLTrackEffPhi2Mtd_

MonitorElement* MtdTracksValidation::meETLTrackEffPhi2Mtd_[2]

private

Definition at line 150 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meETLTrackEffPhiMtd_

MonitorElement* MtdTracksValidation::meETLTrackEffPhiMtd_[2]

private

Definition at line 147 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meETLTrackEffPhiTot_

MonitorElement* MtdTracksValidation::meETLTrackEffPhiTot_[2]

private

Definition at line 144 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meETLTrackEffPt2Mtd_

MonitorElement* MtdTracksValidation::meETLTrackEffPt2Mtd_[2]

private

Definition at line 151 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meETLTrackEffPtMtd_

MonitorElement* MtdTracksValidation::meETLTrackEffPtMtd_[2]

private

Definition at line 148 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meETLTrackEffPtTot_

MonitorElement* MtdTracksValidation::meETLTrackEffPtTot_[2]

private

Definition at line 145 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meETLTrackPtRes_

MonitorElement* MtdTracksValidation::meETLTrackPtRes_

private

Definition at line 152 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meETLTrackRPTime_

MonitorElement* MtdTracksValidation::meETLTrackRPTime_

private

Definition at line 142 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackEffEtaTot_

MonitorElement* MtdTracksValidation::meMVATrackEffEtaTot_

private

Definition at line 173 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackEffPtTot_

MonitorElement* MtdTracksValidation::meMVATrackEffPtTot_

private

Definition at line 165 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackMatchedEffEtaMtd_

MonitorElement* MtdTracksValidation::meMVATrackMatchedEffEtaMtd_

private

Definition at line 175 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackMatchedEffEtaTot_

MonitorElement* MtdTracksValidation::meMVATrackMatchedEffEtaTot_

private

Definition at line 174 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackMatchedEffPtMtd_

MonitorElement* MtdTracksValidation::meMVATrackMatchedEffPtMtd_

private

Definition at line 167 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackMatchedEffPtTot_

MonitorElement* MtdTracksValidation::meMVATrackMatchedEffPtTot_

private

Definition at line 166 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackPullTot_

MonitorElement* MtdTracksValidation::meMVATrackPullTot_

private

Definition at line 182 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackResTot_

MonitorElement* MtdTracksValidation::meMVATrackResTot_

private

Definition at line 181 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meMVATrackZposResTot_

MonitorElement* MtdTracksValidation::meMVATrackZposResTot_

private

Definition at line 183 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meNTrackingParticles_

MonitorElement* MtdTracksValidation::meNTrackingParticles_

private

Definition at line 189 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackMatchedTPEffEtaEtl2Mtd_

MonitorElement* MtdTracksValidation::meTrackMatchedTPEffEtaEtl2Mtd_

private

Definition at line 178 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackMatchedTPEffEtaMtd_

MonitorElement* MtdTracksValidation::meTrackMatchedTPEffEtaMtd_

private

Definition at line 177 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackMatchedTPEffEtaTot_

MonitorElement* MtdTracksValidation::meTrackMatchedTPEffEtaTot_

private

Definition at line 176 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackMatchedTPEffPtEtl2Mtd_

MonitorElement* MtdTracksValidation::meTrackMatchedTPEffPtEtl2Mtd_

private

Definition at line 170 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackMatchedTPEffPtMtd_

MonitorElement* MtdTracksValidation::meTrackMatchedTPEffPtMtd_

private

Definition at line 169 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackMatchedTPEffPtTot_

MonitorElement* MtdTracksValidation::meTrackMatchedTPEffPtTot_

private

Definition at line 168 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackMatchedTPmtdEffEtaMtd_

MonitorElement* MtdTracksValidation::meTrackMatchedTPmtdEffEtaMtd_

private

Definition at line 180 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackMatchedTPmtdEffEtaTot_

MonitorElement* MtdTracksValidation::meTrackMatchedTPmtdEffEtaTot_

private

Definition at line 179 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackMatchedTPmtdEffPtMtd_

MonitorElement* MtdTracksValidation::meTrackMatchedTPmtdEffPtMtd_

private

Definition at line 172 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackMatchedTPmtdEffPtTot_

MonitorElement* MtdTracksValidation::meTrackMatchedTPmtdEffPtTot_

private

Definition at line 171 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackMVAQual_

MonitorElement* MtdTracksValidation::meTrackMVAQual_

private

Definition at line 162 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackNumHits_

MonitorElement* MtdTracksValidation::meTrackNumHits_

private

Definition at line 161 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackPathLenghtvsEta_

MonitorElement* MtdTracksValidation::meTrackPathLenghtvsEta_

private

Definition at line 163 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackSigmat0Pid_

MonitorElement* MtdTracksValidation::meTrackSigmat0Pid_

private

Definition at line 158 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackSigmat0SafePid_

MonitorElement* MtdTracksValidation::meTrackSigmat0SafePid_

private

Definition at line 160 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackSigmat0Src_

MonitorElement* MtdTracksValidation::meTrackSigmat0Src_

private

Definition at line 156 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackt0Pid_

MonitorElement* MtdTracksValidation::meTrackt0Pid_

private

Definition at line 157 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackt0SafePid_

MonitorElement* MtdTracksValidation::meTrackt0SafePid_

private

Definition at line 159 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTrackt0Src_

MonitorElement* MtdTracksValidation::meTrackt0Src_

private

Definition at line 155 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meTracktmtd_

MonitorElement* MtdTracksValidation::meTracktmtd_

private

Definition at line 154 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meUnassCrysEnergy_

MonitorElement* MtdTracksValidation::meUnassCrysEnergy_

private

Definition at line 186 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meUnassDeposit_

MonitorElement* MtdTracksValidation::meUnassDeposit_

private

Definition at line 188 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meUnassLgadsEnergy_

MonitorElement* MtdTracksValidation::meUnassLgadsEnergy_

private

Definition at line 187 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

meUnassociatedDetId_

MonitorElement* MtdTracksValidation::meUnassociatedDetId_

private

Definition at line 185 of file MtdTracksValidation.cc.

Referenced by analyze(), and bookHistograms().

mtdgeoToken_

edm::ESGetToken<MTDGeometry, MTDDigiGeometryRecord> MtdTracksValidation::mtdgeoToken_

private

Definition at line 129 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

mtdtopoToken_

edm::ESGetToken<MTDTopology, MTDTopologyRcd> MtdTracksValidation::mtdtopoToken_

private

Definition at line 130 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

optionalPlots_

bool MtdTracksValidation::optionalPlots_

private

Definition at line 100 of file MtdTracksValidation.cc.

Referenced by analyze().

particleTableToken_

edm::ESGetToken<HepPDT::ParticleDataTable, edm::DefaultRecord> MtdTracksValidation::particleTableToken_

private

Definition at line 131 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

pathLengthToken_

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

private

Definition at line 117 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

pTcut_

constexpr double MtdTracksValidation::pTcut_ = 0.7

staticprivate

Definition at line 90 of file MtdTracksValidation.cc.

Referenced by mvaGenSel(), mvaRecSel(), and mvaTPSel().

r2s_

const reco::RecoToSimCollection* MtdTracksValidation::r2s_

private

Definition at line 102 of file MtdTracksValidation.cc.

Referenced by analyze(), and getMatchedTP().

rBTL_

constexpr double MtdTracksValidation::rBTL_ = 110.0

staticprivate

Definition at line 96 of file MtdTracksValidation.cc.

Referenced by mvaTPSel().

recoToSimAssociationToken_

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

private

Definition at line 112 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

RecTrackToken_

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

private

Definition at line 106 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

RecVertexToken_

edm::EDGetTokenT<std::vector<reco::Vertex> > MtdTracksValidation::RecVertexToken_

private

Definition at line 107 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

s2r_

const reco::SimToRecoCollection* MtdTracksValidation::s2r_

private

Definition at line 103 of file MtdTracksValidation.cc.

Referenced by analyze().

Sigmat0PidToken_

edm::EDGetTokenT<edm::ValueMap<float> > MtdTracksValidation::Sigmat0PidToken_

private

Definition at line 124 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

Sigmat0SafePidToken_

edm::EDGetTokenT<edm::ValueMap<float> > MtdTracksValidation::Sigmat0SafePidToken_

private

Definition at line 126 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

Sigmat0SrcToken_

edm::EDGetTokenT<edm::ValueMap<float> > MtdTracksValidation::Sigmat0SrcToken_

private

Definition at line 122 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

SigmatmtdToken_

edm::EDGetTokenT<edm::ValueMap<float> > MtdTracksValidation::SigmatmtdToken_

private

Definition at line 120 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

simToRecoAssociationToken_

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

private

Definition at line 111 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

t0PidToken_

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

private

Definition at line 123 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

t0SafePidToken_

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

private

Definition at line 125 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

t0SrcToken_

edm::EDGetTokenT<edm::ValueMap<float> > MtdTracksValidation::t0SrcToken_

private

Definition at line 121 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

tmtdToken_

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

private

Definition at line 119 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

trackAssocToken_

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

private

Definition at line 116 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

trackingParticleCollectionToken_

edm::EDGetTokenT<TrackingParticleCollection> MtdTracksValidation::trackingParticleCollectionToken_

private

Definition at line 110 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

trackMaxBtlEta_

const float MtdTracksValidation::trackMaxBtlEta_

private

Definition at line 84 of file MtdTracksValidation.cc.

Referenced by analyze().

trackMaxEtlEta_

const float MtdTracksValidation::trackMaxEtlEta_

private

Definition at line 86 of file MtdTracksValidation.cc.

Referenced by analyze(), and isETL().

trackMinEtlEta_

const float MtdTracksValidation::trackMinEtlEta_

private

Definition at line 85 of file MtdTracksValidation.cc.

Referenced by analyze(), and isETL().

trackMinPt_

const float MtdTracksValidation::trackMinPt_

private

Definition at line 83 of file MtdTracksValidation.cc.

Referenced by analyze().

trackMVAQualToken_

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

private

Definition at line 127 of file MtdTracksValidation.cc.

Referenced by analyze(), and MtdTracksValidation().

zETL_

constexpr double MtdTracksValidation::zETL_ = 290.0

staticprivate

Definition at line 97 of file MtdTracksValidation.cc.

Referenced by mvaTPSel().