#include <TkAlCaRecoMonitor.h>

Inheritance diagram for TkAlCaRecoMonitor:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

	TkAlCaRecoMonitor (const edm::ParameterSet &)

	~TkAlCaRecoMonitor () override=default

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

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	fillHitmaps (const reco::Track &track, const TrackerGeometry &geometry)

void	fillRawIdMap (const TrackerGeometry &geometry)

Private Attributes
MonitorElement *	AlCaRecoTrackEfficiency_

std::map< int, int >	binByRawId_

edm::ParameterSet	conf_

const double	daughterMass_

const bool	fillInvariantMass_

const bool	fillRawIdMap_

MonitorElement *	Hits_perDetId_

MonitorElement *	Hits_XvsY_

MonitorElement *	Hits_ZvsR_

MonitorElement *	invariantMass_

const edm::EDGetTokenT< reco::CaloJetCollection >	jetCollection_

MonitorElement *	jetPt_

const double	maxJetPt_

const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	mfToken_

MonitorElement *	minJetDeltaR_

MonitorElement *	minTrackDeltaR_

const edm::EDGetTokenT< reco::TrackCollection >	referenceTrackProducer_

const bool	runsOnReco_

MonitorElement *	sumCharge_

const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tkGeomToken_

MonitorElement *	TrackCurvature_

const edm::EDGetTokenT< reco::TrackCollection >	trackProducer_

MonitorElement *	TrackPtNegative_

MonitorElement *	TrackPtPositive_

MonitorElement *	TrackQuality_

const bool	useSignedR_

Static Private Attributes
static constexpr const double	kMuonMass_ = 0.1056583755

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

DQM/TrackerMonitorTrack/src/TkAlCaRecoMonitor.cc Monitoring special quantities related to Tracker Alignment AlCaReco Production.

Definition at line 34 of file TkAlCaRecoMonitor.h.

Constructor & Destructor Documentation

◆ TkAlCaRecoMonitor()

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

explicit

Definition at line 19 of file TkAlCaRecoMonitor.cc.

References conf_.

     : tkGeomToken_(esConsumes()),
       mfToken_(esConsumes()),
       trackProducer_(consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("TrackProducer"))),
       referenceTrackProducer_(
           consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("ReferenceTrackProducer"))),
       jetCollection_(mayConsume<reco::CaloJetCollection>(iConfig.getParameter<edm::InputTag>("CaloJetCollection"))),
       daughterMass_(iConfig.getParameter<double>("daughterMass")),
       maxJetPt_(iConfig.getParameter<double>("maxJetPt")),
       fillInvariantMass_(iConfig.getParameter<bool>("fillInvariantMass")),
       fillRawIdMap_(iConfig.getParameter<bool>("fillRawIdMap")),
       runsOnReco_(iConfig.getParameter<bool>("runsOnReco")),
       useSignedR_(iConfig.getParameter<bool>("useSignedR")) {
   // copy configuration object to use it in bookHistograms
   conf_ = iConfig;
 }

◆ ~TkAlCaRecoMonitor()

TkAlCaRecoMonitor::~TkAlCaRecoMonitor ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 185 of file TkAlCaRecoMonitor.cc.

References AlCaRecoTrackEfficiency_, B, binByRawId_, daughterMass_, HLTMuonOfflineAnalyzer_cfi::deltaR2, MillePedeFileConverter_cfg::e, dqm::impl::MonitorElement::Fill(), fillHitmaps(), fillInvariantMass_, fillRawIdMap(), fillRawIdMap_, edm::EventSetup::getHandle(), mps_fire::i, iEvent, invariantMass_, jetCollection_, jetPt_, PDWG_EXODelayedJetMET_cff::jets, HLT_2023v12_cff::magneticField, maxJetPt_, mfToken_, minJetDeltaR_, minTrackDeltaR_, reco::TrackBase::qualitySize, referenceTrackProducer_, runsOnReco_, mathSSE::sqrt(), sumCharge_, tkGeomToken_, HLT_2023v12_cff::track, JetHT_cfg::trackCollection, TrackCurvature_, trackProducer_, TrackPtNegative_, TrackPtPositive_, and TrackQuality_.

 {
   edm::Handle<reco::TrackCollection> trackCollection;
   iEvent.getByToken(trackProducer_, trackCollection);
   if (!trackCollection.isValid()) {
     edm::LogError("Alignment") << "invalid trackcollection encountered!";
     return;
   }
 
   edm::Handle<reco::TrackCollection> referenceTrackCollection;
   iEvent.getByToken(referenceTrackProducer_, referenceTrackCollection);
   if (!trackCollection.isValid()) {
     edm::LogError("Alignment") << "invalid reference track-collection encountered!";
     return;
   }
 
   const auto &geometry = iSetup.getHandle(tkGeomToken_);
   if (!geometry.isValid()) {
     edm::LogError("Alignment") << "invalid geometry found in event setup!";
   }
 
   const auto &magneticField = iSetup.getHandle(mfToken_);
   if (!magneticField.isValid()) {
     edm::LogError("Alignment") << "invalid magnetic field configuration encountered!";
     return;
   }
 
   edm::Handle<reco::CaloJetCollection> jets;
   if (runsOnReco_) {
     iEvent.getByToken(jetCollection_, jets);
     if (!jets.isValid()) {
       edm::LogError("Alignment") << "no jet collection found in event!";
     }
   }
   // fill only once - not yet in beginJob since no access to geometry
   if (fillRawIdMap_ && binByRawId_.empty())
     this->fillRawIdMap(*geometry);
 
   AlCaRecoTrackEfficiency_->Fill(static_cast<double>((*trackCollection).size()) / (*referenceTrackCollection).size());
 
   double sumOfCharges = 0;
   for (const auto &track : *trackCollection) {
     double dR2 = 0;
     if (runsOnReco_) {
       double minJetDeltaR2 = 10 * 10;  // some number > 2pi
       for (const auto &itJet : *jets) {
         jetPt_->Fill(itJet.pt());
         dR2 = deltaR2(track, itJet);
         if (itJet.pt() > maxJetPt_ && dR2 < minJetDeltaR2)
           minJetDeltaR2 = dR2;
 
         // edm::LogInfo("Alignment") <<">  isolated: "<< isolated << " jetPt "<<
         // (*itJet).pt() <<" deltaR: "<< deltaR(*(*it),(*itJet)) ;
       }
       minJetDeltaR_->Fill(std::sqrt(minJetDeltaR2));
     }
 
     double minTrackDeltaR2 = 10 * 10;  // some number > 2pi
     for (const auto &track2 : *trackCollection) {
       dR2 = deltaR2(track, track2);
       if (dR2 < minTrackDeltaR2 && dR2 > 1e-12)
         minTrackDeltaR2 = dR2;
     }
 
     for (int i = 0; i < reco::TrackBase::qualitySize; ++i) {
       if (track.quality(reco::TrackBase::TrackQuality(i))) {
         TrackQuality_->Fill(i);
       }
     }
 
     GlobalPoint gPoint(track.vx(), track.vy(), track.vz());
     double B = magneticField->inTesla(gPoint).z();
     double curv = -track.charge() * 0.002998 * B / track.pt();
     TrackCurvature_->Fill(curv);
 
     track.charge() > 0 ? TrackPtPositive_->Fill(track.pt()) : TrackPtNegative_->Fill(track.pt());
 
     minTrackDeltaR_->Fill(std::sqrt(minTrackDeltaR2));
     fillHitmaps(track, *geometry);
     sumOfCharges += track.charge();
   }
 
   sumCharge_->Fill(sumOfCharges);
 
   if (fillInvariantMass_) {
     if ((*trackCollection).size() == 2) {
       TLorentzVector track0(
           (*trackCollection).at(0).px(),
           (*trackCollection).at(0).py(),
           (*trackCollection).at(0).pz(),
           sqrt(((*trackCollection).at(0).p() * (*trackCollection).at(0).p()) + daughterMass_ * daughterMass_));
       TLorentzVector track1(
           (*trackCollection).at(1).px(),
           (*trackCollection).at(1).py(),
           (*trackCollection).at(1).pz(),
           sqrt(((*trackCollection).at(1).p() * (*trackCollection).at(1).p()) + daughterMass_ * daughterMass_));
       TLorentzVector mother = track0 + track1;
 
       invariantMass_->Fill(mother.M());
     } else {
       edm::LogInfo("Alignment") << "wrong number of tracks trackcollection encountered: " << (*trackCollection).size();
     }
   }
 }

◆ bookHistograms()

void TkAlCaRecoMonitor::bookHistograms	(	DQMStore::IBooker &	iBooker,
		edm::Run const &	,
		edm::EventSetup const &
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 36 of file TkAlCaRecoMonitor.cc.

                                                                                                       {
   std::string histname;  // for naming the histograms according to algorithm used
 
   std::string AlgoName = conf_.getParameter<std::string>("AlgoName");
   std::string MEFolderName = conf_.getParameter<std::string>("FolderName");
   iBooker.setCurrentFolder(MEFolderName + "/TkAlignmentSpecific");
 
   //
   unsigned int MassBin = conf_.getParameter<unsigned int>("MassBin");
   double MassMin = conf_.getParameter<double>("MassMin");
   double MassMax = conf_.getParameter<double>("MassMax");
 
   if (fillInvariantMass_) {
     histname = "InvariantMass_";
     invariantMass_ = iBooker.book1D(histname + AlgoName, histname + AlgoName, MassBin, MassMin, MassMax);
     invariantMass_->setAxisTitle("invariant Mass / GeV");
   } else {
     invariantMass_ = nullptr;
   }
 
   unsigned int TrackPtPositiveBin = conf_.getParameter<unsigned int>("TrackPtBin");
   double TrackPtPositiveMin = conf_.getParameter<double>("TrackPtMin");
   double TrackPtPositiveMax = conf_.getParameter<double>("TrackPtMax");
 
   histname = "TrackPtPositive_";
   TrackPtPositive_ = iBooker.book1D(
       histname + AlgoName, histname + AlgoName, TrackPtPositiveBin, TrackPtPositiveMin, TrackPtPositiveMax);
   TrackPtPositive_->setAxisTitle("p_{T} of tracks charge > 0");
 
   unsigned int TrackPtNegativeBin = conf_.getParameter<unsigned int>("TrackPtBin");
   double TrackPtNegativeMin = conf_.getParameter<double>("TrackPtMin");
   double TrackPtNegativeMax = conf_.getParameter<double>("TrackPtMax");
 
   histname = "TrackPtNegative_";
   TrackPtNegative_ = iBooker.book1D(
       histname + AlgoName, histname + AlgoName, TrackPtNegativeBin, TrackPtNegativeMin, TrackPtNegativeMax);
   TrackPtNegative_->setAxisTitle("p_{T} of tracks charge < 0");
 
   histname = "TrackQuality_";
   TrackQuality_ = iBooker.book1D(
       histname + AlgoName, histname + AlgoName, reco::TrackBase::qualitySize, -0.5, reco::TrackBase::qualitySize - 0.5);
   TrackQuality_->setAxisTitle("quality");
   for (int i = 0; i < reco::TrackBase::qualitySize; ++i) {
     TrackQuality_->getTH1()->GetXaxis()->SetBinLabel(
         i + 1, reco::TrackBase::qualityName(reco::TrackBase::TrackQuality(i)).c_str());
   }
 
   unsigned int SumChargeBin = conf_.getParameter<unsigned int>("SumChargeBin");
   double SumChargeMin = conf_.getParameter<double>("SumChargeMin");
   double SumChargeMax = conf_.getParameter<double>("SumChargeMax");
 
   histname = "SumCharge_";
   sumCharge_ = iBooker.book1D(histname + AlgoName, histname + AlgoName, SumChargeBin, SumChargeMin, SumChargeMax);
   sumCharge_->setAxisTitle("#SigmaCharge");
 
   unsigned int TrackCurvatureBin = conf_.getParameter<unsigned int>("TrackCurvatureBin");
   double TrackCurvatureMin = conf_.getParameter<double>("TrackCurvatureMin");
   double TrackCurvatureMax = conf_.getParameter<double>("TrackCurvatureMax");
 
   histname = "TrackCurvature_";
   TrackCurvature_ =
       iBooker.book1D(histname + AlgoName, histname + AlgoName, TrackCurvatureBin, TrackCurvatureMin, TrackCurvatureMax);
   TrackCurvature_->setAxisTitle("#kappa track");
 
   if (runsOnReco_) {
     unsigned int JetPtBin = conf_.getParameter<unsigned int>("JetPtBin");
     double JetPtMin = conf_.getParameter<double>("JetPtMin");
     double JetPtMax = conf_.getParameter<double>("JetPtMax");
 
     histname = "JetPt_";
     jetPt_ = iBooker.book1D(histname + AlgoName, histname + AlgoName, JetPtBin, JetPtMin, JetPtMax);
     jetPt_->setAxisTitle("jet p_{T} / GeV");
 
     unsigned int MinJetDeltaRBin = conf_.getParameter<unsigned int>("MinJetDeltaRBin");
     double MinJetDeltaRMin = conf_.getParameter<double>("MinJetDeltaRMin");
     double MinJetDeltaRMax = conf_.getParameter<double>("MinJetDeltaRMax");
 
     histname = "MinJetDeltaR_";
     minJetDeltaR_ =
         iBooker.book1D(histname + AlgoName, histname + AlgoName, MinJetDeltaRBin, MinJetDeltaRMin, MinJetDeltaRMax);
     minJetDeltaR_->setAxisTitle("minimal Jet #DeltaR / rad");
   } else {
     jetPt_ = nullptr;
     minJetDeltaR_ = nullptr;
   }
 
   unsigned int MinTrackDeltaRBin = conf_.getParameter<unsigned int>("MinTrackDeltaRBin");
   double MinTrackDeltaRMin = conf_.getParameter<double>("MinTrackDeltaRMin");
   double MinTrackDeltaRMax = conf_.getParameter<double>("MinTrackDeltaRMax");
 
   histname = "MinTrackDeltaR_";
   minTrackDeltaR_ =
       iBooker.book1D(histname + AlgoName, histname + AlgoName, MinTrackDeltaRBin, MinTrackDeltaRMin, MinTrackDeltaRMax);
   minTrackDeltaR_->setAxisTitle("minimal Track #DeltaR / rad");
 
   unsigned int TrackEfficiencyBin = conf_.getParameter<unsigned int>("TrackEfficiencyBin");
   double TrackEfficiencyMin = conf_.getParameter<double>("TrackEfficiencyMin");
   double TrackEfficiencyMax = conf_.getParameter<double>("TrackEfficiencyMax");
 
   histname = "AlCaRecoTrackEfficiency_";
   AlCaRecoTrackEfficiency_ = iBooker.book1D(
       histname + AlgoName, histname + AlgoName, TrackEfficiencyBin, TrackEfficiencyMin, TrackEfficiencyMax);
   Labels l_tp, l_rtp;
   labelsForToken(referenceTrackProducer_, l_rtp);
   labelsForToken(trackProducer_, l_tp);
   AlCaRecoTrackEfficiency_->setAxisTitle("n(" + std::string(l_tp.module) + ") / n(" + std::string(l_rtp.module) + ")");
 
   int zBin = conf_.getParameter<unsigned int>("HitMapsZBin");  // 300
   double zMax = conf_.getParameter<double>("HitMapZMax");      // 300.0; //cm
 
   int rBin = conf_.getParameter<unsigned int>("HitMapsRBin");  // 120;
   double rMax = conf_.getParameter<double>("HitMapRMax");      // 120.0; //cm
 
   histname = "Hits_ZvsR_";
   double rMin = 0.0;
   if (useSignedR_)
     rMin = -rMax;
 
   Hits_ZvsR_ = iBooker.book2D(histname + AlgoName, histname + AlgoName, zBin, -zMax, zMax, rBin, rMin, rMax);
 
   histname = "Hits_XvsY_";
   Hits_XvsY_ = iBooker.book2D(histname + AlgoName, histname + AlgoName, rBin, -rMax, rMax, rBin, -rMax, rMax);
 
   if (fillRawIdMap_) {
     histname = "Hits_perDetId_";
 
     // leads to differences in axsis between samples??
     // int nModules = binByRawId_.size();
     // Hits_perDetId_ = iBooker.book1D(histname+AlgoName, histname+AlgoName,
     // nModules, static_cast<double>(nModules) -0.5,
     // static_cast<double>(nModules) -0.5);
     Hits_perDetId_ = iBooker.book1D(histname + AlgoName, histname + AlgoName, 16601, -0.5, 16600.5);
     Hits_perDetId_->setAxisTitle("rawId Bins");
 
     //  std::stringstream binLabel;
     //  for( std::map<int,int>::iterator it = binByRawId_.begin(); it !=
     //  binByRawId_.end(); ++it ){
     //    binLabel.str() = "";
     //    binLabel << (*it).first;
     //    Hits_perDetId_->getTH1()->GetXaxis()->SetBinLabel( (*it).second +1,
     //    binLabel.str().c_str());
     //  }
   }
 }

◆ fillDescriptions()

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

static

Definition at line 335 of file TkAlCaRecoMonitor.cc.

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

 {
   edm::ParameterSetDescription desc;
   desc.setComment("Generic track analyzer to check ALCARECO Tracker Alignment specific sample quantities");
   desc.add<edm::InputTag>("TrackProducer", edm::InputTag("generalTracks"));
   desc.add<edm::InputTag>("ReferenceTrackProducer", edm::InputTag("generalTrakcs"));
   desc.add<edm::InputTag>("CaloJetCollection", edm::InputTag("ak4CaloJets"));
   desc.add<std::string>("AlgoName", "testTkAlCaReco");
   desc.add<std::string>("FolderName", "TkAlCaRecoMonitor");
   desc.add<double>("daughterMass", kMuonMass_)->setComment("GeV");
   desc.add<double>("maxJetPt", 10.)->setComment("GeV");
   desc.add<bool>("fillInvariantMass", false);
   desc.add<bool>("runsOnReco", false);
   desc.add<bool>("useSignedR", false);
   desc.add<bool>("fillRawIdMap", false);
   desc.add<unsigned int>("MassBin", 100);
   desc.add<double>("MassMin", 0.0);
   desc.add<double>("MassMax", 100.0);
   desc.add<unsigned int>("TrackPtBin", 110);
   desc.add<double>("TrackPtMin", 0.0);
   desc.add<double>("TrackPtMax", 110.);
   desc.add<unsigned int>("TrackCurvatureBin", 2000);
   desc.add<double>("TrackCurvatureMin", -0.01)->setComment("1/GeV");
   desc.add<double>("TrackCurvatureMax", 0.01)->setComment("1/GeV");
   desc.add<unsigned int>("SumChargeBin", 11);
   desc.add<double>("SumChargeMin", -5.5);
   desc.add<double>("SumChargeMax", 5.5);
   desc.add<unsigned int>("JetPtBin", 100);
   desc.add<double>("JetPtMin", 0.0);
   desc.add<double>("JetPtMax", 50.0);
   desc.add<unsigned int>("MinJetDeltaRBin", 100);
   desc.add<double>("MinJetDeltaRMin", 0);
   desc.add<double>("MinJetDeltaRMax", 10);
   desc.add<unsigned int>("MinTrackDeltaRBin", 100);
   desc.add<double>("MinTrackDeltaRMin", 0);
   desc.add<double>("MinTrackDeltaRMax", 3.2);
   desc.add<unsigned int>("TrackEfficiencyBin", 102);
   desc.add<double>("TrackEfficiencyMin", -0.01);
   desc.add<double>("TrackEfficiencyMax", 1.01);
   desc.add<unsigned int>("HitMapsZBin", 300);
   desc.add<double>("HitMapZMax", 300.)->setComment("cm");
   desc.add<unsigned int>("HitMapsRBin", 120);
   desc.add<double>("HitMapRMax", 120.)->setComment("cm");
   descriptions.addWithDefaultLabel(desc);
 }

◆ fillHitmaps()

void TkAlCaRecoMonitor::fillHitmaps	(	const reco::Track &	track,
		const TrackerGeometry &	geometry
	)

private

Definition at line 292 of file TkAlCaRecoMonitor.cc.

References binByRawId_, dqm::impl::MonitorElement::Fill(), fillRawIdMap_, Hits_perDetId_, Hits_XvsY_, Hits_ZvsR_, alignCSCRings::r, mathSSE::sqrt(), HLT_2023v12_cff::track, and useSignedR_.

Referenced by analyze().

 {
   for (auto const &iHit : track.recHits()) {
     if (iHit->isValid()) {
       const DetId geoId(iHit->geographicalId());
       const GeomDet *gd = geometry.idToDet(geoId);
       // since 2_1_X local hit positions are transient. taking center of the hit
       // module for now. The alternative would be the coarse estimation or a
       // refit.
       // const GlobalPoint globP( gd->toGlobal( hit->localPosition() ) );
       const GlobalPoint globP(gd->toGlobal(Local3DPoint(0., 0., 0.)));
       double r = sqrt(globP.x() * globP.x() + globP.y() * globP.y());
       if (useSignedR_)
         r *= globP.y() / fabs(globP.y());
 
       Hits_ZvsR_->Fill(globP.z(), r);
       Hits_XvsY_->Fill(globP.x(), globP.y());
 
       if (fillRawIdMap_)
         Hits_perDetId_->Fill(binByRawId_[geoId.rawId()]);
     }
   }
 }

◆ fillRawIdMap()

void TkAlCaRecoMonitor::fillRawIdMap ( const TrackerGeometry & geometry )

private

Definition at line 318 of file TkAlCaRecoMonitor.cc.

References binByRawId_, mps_fire::i, and jetUpdater_cfi::sort.

Referenced by analyze().

 {
   std::vector<int> sortedRawIds;
   for (const auto &iDetId : geometry.detUnitIds()) {
     sortedRawIds.push_back(iDetId.rawId());
   }
   std::sort(sortedRawIds.begin(), sortedRawIds.end());
 
   int i = 0;
   for (const auto &iRawId : sortedRawIds) {
     binByRawId_[iRawId] = i;
     ++i;
   }
 }

Member Data Documentation

◆ AlCaRecoTrackEfficiency_

MonitorElement* TkAlCaRecoMonitor::AlCaRecoTrackEfficiency_

private

Definition at line 74 of file TkAlCaRecoMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ binByRawId_

std::map<int, int> TkAlCaRecoMonitor::binByRawId_

private

Definition at line 83 of file TkAlCaRecoMonitor.h.

Referenced by analyze(), fillHitmaps(), and fillRawIdMap().

◆ conf_

edm::ParameterSet TkAlCaRecoMonitor::conf_

private

Definition at line 65 of file TkAlCaRecoMonitor.h.

Referenced by bookHistograms(), and TkAlCaRecoMonitor().

◆ daughterMass_

const double TkAlCaRecoMonitor::daughterMass_

private

Definition at line 58 of file TkAlCaRecoMonitor.h.

Referenced by analyze().

◆ fillInvariantMass_

const bool TkAlCaRecoMonitor::fillInvariantMass_

private

Definition at line 60 of file TkAlCaRecoMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ fillRawIdMap_

const bool TkAlCaRecoMonitor::fillRawIdMap_

private

Definition at line 61 of file TkAlCaRecoMonitor.h.

Referenced by analyze(), bookHistograms(), and fillHitmaps().

◆ Hits_perDetId_

MonitorElement* TkAlCaRecoMonitor::Hits_perDetId_

private

Definition at line 75 of file TkAlCaRecoMonitor.h.

Referenced by bookHistograms(), and fillHitmaps().

◆ Hits_XvsY_

MonitorElement* TkAlCaRecoMonitor::Hits_XvsY_

private

Definition at line 81 of file TkAlCaRecoMonitor.h.

Referenced by bookHistograms(), and fillHitmaps().

◆ Hits_ZvsR_

MonitorElement* TkAlCaRecoMonitor::Hits_ZvsR_

private

Definition at line 80 of file TkAlCaRecoMonitor.h.

Referenced by bookHistograms(), and fillHitmaps().

◆ invariantMass_

MonitorElement* TkAlCaRecoMonitor::invariantMass_

private

Definition at line 68 of file TkAlCaRecoMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ jetCollection_

const edm::EDGetTokenT<reco::CaloJetCollection> TkAlCaRecoMonitor::jetCollection_

private

Definition at line 56 of file TkAlCaRecoMonitor.h.

Referenced by analyze().

◆ jetPt_

MonitorElement* TkAlCaRecoMonitor::jetPt_

private

Definition at line 71 of file TkAlCaRecoMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ kMuonMass_

constexpr const double TkAlCaRecoMonitor::kMuonMass_ = 0.1056583755

staticprivate

Definition at line 45 of file TkAlCaRecoMonitor.h.

Referenced by fillDescriptions().

◆ maxJetPt_

const double TkAlCaRecoMonitor::maxJetPt_

private

Definition at line 59 of file TkAlCaRecoMonitor.h.

Referenced by analyze().

◆ mfToken_

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> TkAlCaRecoMonitor::mfToken_

private

Definition at line 52 of file TkAlCaRecoMonitor.h.

Referenced by analyze().

◆ minJetDeltaR_

MonitorElement* TkAlCaRecoMonitor::minJetDeltaR_

private

Definition at line 72 of file TkAlCaRecoMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ minTrackDeltaR_

MonitorElement* TkAlCaRecoMonitor::minTrackDeltaR_

private

Definition at line 73 of file TkAlCaRecoMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ referenceTrackProducer_

const edm::EDGetTokenT<reco::TrackCollection> TkAlCaRecoMonitor::referenceTrackProducer_

private

Definition at line 55 of file TkAlCaRecoMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ runsOnReco_

const bool TkAlCaRecoMonitor::runsOnReco_

private

Definition at line 62 of file TkAlCaRecoMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ sumCharge_

MonitorElement* TkAlCaRecoMonitor::sumCharge_

private

Definition at line 69 of file TkAlCaRecoMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ tkGeomToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> TkAlCaRecoMonitor::tkGeomToken_

private

Definition at line 51 of file TkAlCaRecoMonitor.h.

Referenced by analyze().

◆ TrackCurvature_

MonitorElement* TkAlCaRecoMonitor::TrackCurvature_

private

Definition at line 78 of file TkAlCaRecoMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ trackProducer_

const edm::EDGetTokenT<reco::TrackCollection> TkAlCaRecoMonitor::trackProducer_

private

Definition at line 54 of file TkAlCaRecoMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ TrackPtNegative_

MonitorElement* TkAlCaRecoMonitor::TrackPtNegative_

private

Definition at line 77 of file TkAlCaRecoMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ TrackPtPositive_

MonitorElement* TkAlCaRecoMonitor::TrackPtPositive_

private

Definition at line 76 of file TkAlCaRecoMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ TrackQuality_

MonitorElement* TkAlCaRecoMonitor::TrackQuality_

private

Definition at line 70 of file TkAlCaRecoMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ useSignedR_

const bool TkAlCaRecoMonitor::useSignedR_

private

Definition at line 63 of file TkAlCaRecoMonitor.h.

Referenced by bookHistograms(), and fillHitmaps().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ TkAlCaRecoMonitor()

◆ ~TkAlCaRecoMonitor()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ fillDescriptions()

◆ fillHitmaps()

◆ fillRawIdMap()

Member Data Documentation

◆ AlCaRecoTrackEfficiency_

◆ binByRawId_

◆ conf_

◆ daughterMass_

◆ fillInvariantMass_

◆ fillRawIdMap_

◆ Hits_perDetId_

◆ Hits_XvsY_

◆ Hits_ZvsR_

◆ invariantMass_

◆ jetCollection_

◆ jetPt_

◆ kMuonMass_

◆ maxJetPt_

◆ mfToken_

◆ minJetDeltaR_

◆ minTrackDeltaR_

◆ referenceTrackProducer_

◆ runsOnReco_

◆ sumCharge_

◆ tkGeomToken_

◆ TrackCurvature_

◆ trackProducer_

◆ TrackPtNegative_

◆ TrackPtPositive_

◆ TrackQuality_

◆ useSignedR_