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 bool	mvaGenRecMatch (const HepMC::GenParticle &, const double &, const reco::TrackBase &)

const bool	mvaGenSel (const HepMC::GenParticle &, const float &)

const bool	mvaRecSel (const reco::TrackBase &, const reco::Vertex &, const double &, const double &)

Private Attributes
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 *	meETLTrackEffEtaMtd_ [2]

MonitorElement *	meETLTrackEffEtaTot_ [2]

MonitorElement *	meETLTrackEffPhiMtd_ [2]

MonitorElement *	meETLTrackEffPhiTot_ [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 *	meTrackMVAQual_

MonitorElement *	meTrackNumHits_

MonitorElement *	meTrackPathLenghtvsEta_

MonitorElement *	meTrackSigmat0Pid_

MonitorElement *	meTrackSigmat0SafePid_

MonitorElement *	meTrackSigmat0Src_

MonitorElement *	meTrackt0Pid_

MonitorElement *	meTrackt0SafePid_

MonitorElement *	meTrackt0Src_

MonitorElement *	meTracktmtd_

edm::ESGetToken< MTDTopology, MTDTopologyRcd >	mtdtopoToken_

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

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

edm::EDGetTokenT< reco::TrackCollection >	RecTrackToken_

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

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

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	etacutGEN_ = 4.

static constexpr double	etacutREC_ = 3.

static constexpr double	pTcut_ = 0.7

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

Constructor & Destructor Documentation

◆ MtdTracksValidation()

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

explicit

Definition at line 135 of file MtdTracksValidation.cc.

References GenRecTrackToken_, edm::ParameterSet::getParameter(), HepMCProductToken_, mtdtopoToken_, particleTableToken_, pathLengthToken_, RecTrackToken_, RecVertexToken_, Sigmat0PidToken_, Sigmat0SafePidToken_, Sigmat0SrcToken_, SigmatmtdToken_, t0PidToken_, t0SafePidToken_, t0SrcToken_, tmtdToken_, trackAssocToken_, 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")) {
   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"));
   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"));
   mtdtopoToken_ = esConsumes<MTDTopology, MTDTopologyRcd>();
   particleTableToken_ = esConsumes<HepPDT::ParticleDataTable, edm::DefaultRecord>();
 }

◆ ~MtdTracksValidation()

MtdTracksValidation::~MtdTracksValidation ( )

override

Definition at line 160 of file MtdTracksValidation.cc.

160 {}

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 163 of file MtdTracksValidation.cc.

References funct::abs(), MTDTopologyMode::barphiflat, ALCARECOTkAlJpsiMuMu_cff::charge, angle_units::operators::convertMmToCm(), deltaZcut_, dqm::impl::MonitorElement::Fill(), GenRecTrackToken_, edm::EventSetup::getHandle(), MTDTopology::getMTDTopologyMode(), edm::EventSetup::getTransientHandle(), HepMCProductToken_, iEvent, edm::makeValid(), CaloTowersParam_cfi::mc, meBTLTrackEffEtaMtd_, meBTLTrackEffEtaTot_, meBTLTrackEffPhiMtd_, meBTLTrackEffPhiTot_, meBTLTrackEffPtMtd_, meBTLTrackEffPtTot_, meBTLTrackPtRes_, meBTLTrackRPTime_, meETLTrackEffEtaMtd_, meETLTrackEffEtaTot_, meETLTrackEffPhiMtd_, meETLTrackEffPhiTot_, meETLTrackEffPtMtd_, meETLTrackEffPtTot_, meETLTrackPtRes_, meETLTrackRPTime_, meMVATrackEffEtaTot_, meMVATrackEffPtTot_, meMVATrackMatchedEffEtaMtd_, meMVATrackMatchedEffEtaTot_, meMVATrackMatchedEffPtMtd_, meMVATrackMatchedEffPtTot_, meMVATrackPullTot_, meMVATrackResTot_, meMVATrackZposResTot_, meTrackMVAQual_, meTrackNumHits_, meTrackPathLenghtvsEta_, meTrackSigmat0Pid_, meTrackSigmat0SafePid_, meTrackSigmat0Src_, meTrackt0Pid_, meTrackt0SafePid_, meTrackt0Src_, meTracktmtd_, mtdtopoToken_, mvaGenRecMatch(), mvaGenSel(), mvaRecSel(), ETLDetId::nDisc(), LHEGenericFilter_cfi::ParticleID, particleTableToken_, pathLengthToken_, RecTrackToken_, RecVertexToken_, Sigmat0PidToken_, Sigmat0SafePidToken_, Sigmat0SrcToken_, SigmatmtdToken_, t0PidToken_, t0SafePidToken_, t0SrcToken_, tmtdToken_, HLT_2022v12_cff::track, trackAssocToken_, trackMaxBtlEta_, trackMaxEtlEta_, trackMinEtlEta_, trackMinPt_, trackMVAQualToken_, and MTDDetId::zside().

                                                                                      {
   using namespace edm;
   using namespace geant_units::operators;
   using namespace std;
 
   auto topologyHandle = iSetup.getTransientHandle(mtdtopoToken_);
   const MTDTopology* topology = topologyHandle.product();
 
   bool topo1Dis = false;
   bool topo2Dis = false;
   if (topology->getMTDTopologyMode() <= static_cast<int>(MTDTopologyMode::Mode::barphiflat)) {
     topo1Dis = true;
   }
   if (topology->getMTDTopologyMode() > static_cast<int>(MTDTopologyMode::Mode::barphiflat)) {
     topo2Dis = true;
   }
 
   auto GenRecTrackHandle = makeValid(iEvent.getHandle(GenRecTrackToken_));
   auto RecVertexHandle = makeValid(iEvent.getHandle(RecVertexToken_));
 
   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_);
 
   unsigned int index = 0;
   // --- 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;
 
     if (track.pt() < trackMinPt_ || std::abs(track.eta()) > trackMaxEtlEta_)
       continue;
 
     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) {
         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_)) {
         if ((MTDEtlZnegD1 == true) || (MTDEtlZnegD2 == true)) {
           meETLTrackEffEtaMtd_[0]->Fill(track.eta());
           meETLTrackEffPhiMtd_[0]->Fill(track.phi());
           meETLTrackEffPtMtd_[0]->Fill(track.pt());
         }
       }
       if ((track.eta() > trackMinEtlEta_) && (track.eta() < trackMaxEtlEta_)) {
         if ((MTDEtlZposD1 == true) || (MTDEtlZposD2 == true)) {
           meETLTrackEffEtaMtd_[1]->Fill(track.eta());
           meETLTrackEffPhiMtd_[1]->Fill(track.phi());
           meETLTrackEffPtMtd_[1]->Fill(track.pt());
         }
       }
     }
   }  //RECO tracks loop
 
   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();
 
   // select events with reco vertex close to true simulated primary vertex
   if (std::abs(primRecoVtx.z() - zsim) < deltaZcut_) {
     index = 0;
     for (const auto& trackGen : *GenRecTrackHandle) {
       const reco::TrackRef trackref(iEvent.getHandle(GenRecTrackToken_), index);
       index++;
 
       // select the reconstructed track
 
       if (trackAssoc[trackref] == -1) {
         continue;
       }
 
       bool noCrack = std::abs(trackGen.eta()) < trackMaxBtlEta_ || std::abs(trackGen.eta()) > trackMinEtlEta_;
 
       // reco-gen matching used for MVA quality flag
       if (mvaRecSel(trackGen, primRecoVtx, t0Safe[trackref], Sigmat0Safe[trackref])) {
         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)) {
               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;
             }
           }
         }
       }
     }
   }
 }

◆ bookHistograms()

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

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 416 of file MtdTracksValidation.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::bookProfile(), folder_, meBTLTrackEffEtaMtd_, meBTLTrackEffEtaTot_, meBTLTrackEffPhiMtd_, meBTLTrackEffPhiTot_, meBTLTrackEffPtMtd_, meBTLTrackEffPtTot_, meBTLTrackPtRes_, meBTLTrackRPTime_, meETLTrackEffEtaMtd_, meETLTrackEffEtaTot_, meETLTrackEffPhiMtd_, meETLTrackEffPhiTot_, meETLTrackEffPtMtd_, meETLTrackEffPtTot_, meETLTrackPtRes_, meETLTrackRPTime_, meMVATrackEffEtaTot_, meMVATrackEffPtTot_, meMVATrackMatchedEffEtaMtd_, meMVATrackMatchedEffEtaTot_, meMVATrackMatchedEffPtMtd_, meMVATrackMatchedEffPtTot_, meMVATrackPullTot_, meMVATrackResTot_, meMVATrackZposResTot_, meTrackMVAQual_, meTrackNumHits_, meTrackPathLenghtvsEta_, meTrackSigmat0Pid_, meTrackSigmat0SafePid_, meTrackSigmat0Src_, meTrackt0Pid_, meTrackt0SafePid_, meTrackt0Src_, meTracktmtd_, 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);
   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.);
   meMVATrackEffEtaTot_ = ibook.book1D("MVAEffEtaTot", "Pt of tracks associated to LV; track eta ", 66, 0., 3.3);
   meMVATrackMatchedEffEtaTot_ =
       ibook.book1D("MVAMatchedEffEtaTot", "Pt of tracks associated to LV matched to GEN; track eta ", 66, 0., 3.3);
   meMVATrackMatchedEffEtaMtd_ = ibook.book1D(
       "MVAMatchedEffEtaMtd", "Pt of tracks associated to LV matched to GEN 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);
 }

◆ fillDescriptions()

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

static

Definition at line 493 of file MtdTracksValidation.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, HLT_2022v12_cff::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>("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.);
 
   descriptions.add("mtdTracksValid", desc);
 }

◆ mvaGenRecMatch()

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

private

Definition at line 542 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_;
   return match;
 }

◆ mvaGenSel()

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

private

Definition at line 521 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 530 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_;
   if (st0 > 0.) {
     match = match && std::abs(t0 - vtx.t()) < 3. * st0;
   }
   return match;
 }