#include <SiPixelPhase1MonitorTrackSoA.cc>

Inheritance diagram for SiPixelPhase1MonitorTrackSoA:

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

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

	SiPixelPhase1MonitorTrackSoA (const edm::ParameterSet &)

	~SiPixelPhase1MonitorTrackSoA () override=default

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Attributes
MonitorElement *	hchi2

MonitorElement *	hChi2VsEta

MonitorElement *	hChi2VsPhi

MonitorElement *	heta

MonitorElement *	hnHits

MonitorElement *	hnHitsVsEta

MonitorElement *	hnHitsVsPhi

MonitorElement *	hnLayers

MonitorElement *	hnLayersVsEta

MonitorElement *	hnLayersVsPhi

MonitorElement *	hnLooseAndAboveTracks

MonitorElement *	hnTracks

MonitorElement *	hphi

MonitorElement *	hpt

MonitorElement *	hquality

MonitorElement *	htip

MonitorElement *	hz

pixelTrack::Quality	minQuality_

edm::EDGetTokenT < PixelTrackHeterogeneous >	tokenSoATrack_

std::string	topFolderName_

bool	useQualityCut_

Detailed Description

bookLayer

Definition at line 28 of file SiPixelPhase1MonitorTrackSoA.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 64 of file SiPixelPhase1MonitorTrackSoA.cc.

References edm::ParameterSet::getParameter(), minQuality_, pixelTrack::qualityByName(), AlCaHLTBitMon_QueryRunRegistry::string, tokenSoATrack_, topFolderName_, and useQualityCut_.

                                                                                          {
   tokenSoATrack_ = consumes<PixelTrackHeterogeneous>(iConfig.getParameter<edm::InputTag>("pixelTrackSrc"));
   topFolderName_ = iConfig.getParameter<std::string>("topFolderName");  //"SiPixelHeterogeneous/PixelTrackSoA";
   useQualityCut_ = iConfig.getParameter<bool>("useQualityCut");
   minQuality_ = pixelTrack::qualityByName(iConfig.getParameter<std::string>("minQuality"));
 }

SiPixelPhase1MonitorTrackSoA::~SiPixelPhase1MonitorTrackSoA ( )

overridedefault

Member Function Documentation

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

override

Definition at line 74 of file SiPixelPhase1MonitorTrackSoA.cc.

References HLT_FULL_cff::chi2, PVValHelper::eta, dqm::impl::MonitorElement::Fill(), edm::Event::getHandle(), hchi2, hChi2VsEta, hChi2VsPhi, heta, hnHits, hnHitsVsEta, hnHitsVsPhi, hnLayers, hnLayersVsEta, hnLayersVsPhi, hnLooseAndAboveTracks, hnTracks, hphi, hpt, hquality, htip, hz, HLT_FULL_cff::maxTracks, minQuality_, nHits, mkfit::Config::nLayers, BeamSpotPI::nTracks, phi, DiDispStaMuonMonitor_cfi::pt, quality, HLT_FULL_cff::tip, tokenSoATrack_, useQualityCut_, and ComparisonHelper::zip().

                                                                                               {
   const auto& tsoaHandle = iEvent.getHandle(tokenSoATrack_);
   if (!tsoaHandle.isValid()) {
     edm::LogWarning("SiPixelPhase1MonitorTrackSoA") << "No Track SoA found \n returning!" << std::endl;
     return;
   }
 
   auto const& tsoa = *((tsoaHandle.product())->get());
   auto maxTracks = tsoa.stride();
   auto const* quality = tsoa.qualityData();
   int32_t nTracks = 0;
   int32_t nLooseAndAboveTracks = 0;
 
   for (int32_t it = 0; it < maxTracks; ++it) {
     auto nHits = tsoa.nHits(it);
     auto nLayers = tsoa.nLayers(it);
     if (nHits == 0)
       break;  // this is a guard
     float pt = tsoa.pt(it);
     if (!(pt > 0.))
       continue;
 
     // fill the quality for all tracks
     pixelTrack::Quality qual = tsoa.quality(it);
     hquality->Fill(int(qual));
     nTracks++;
 
     if (useQualityCut_ && quality[it] < minQuality_)
       continue;
 
     // fill parameters only for quality >= loose
     float chi2 = tsoa.chi2(it);
     float phi = tsoa.phi(it);
     float zip = tsoa.zip(it);
     float eta = tsoa.eta(it);
     float tip = tsoa.tip(it);
 
     hchi2->Fill(chi2);
     hChi2VsPhi->Fill(phi, chi2);
     hChi2VsEta->Fill(eta, chi2);
     hnHits->Fill(nHits);
     hnLayers->Fill(nLayers);
     hnHitsVsPhi->Fill(phi, nHits);
     hnHitsVsEta->Fill(eta, nHits);
     hnLayersVsPhi->Fill(phi, nLayers);
     hnLayersVsEta->Fill(eta, nLayers);
     hpt->Fill(pt);
     heta->Fill(eta);
     hphi->Fill(phi);
     hz->Fill(zip);
     htip->Fill(tip);
     nLooseAndAboveTracks++;
   }
   hnTracks->Fill(nTracks);
   hnLooseAndAboveTracks->Fill(nLooseAndAboveTracks);
 }

void SiPixelPhase1MonitorTrackSoA::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	iRun,
		edm::EventSetup const &	iSetup
	)

override

Definition at line 134 of file SiPixelPhase1MonitorTrackSoA.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::bookProfile(), dqm::implementation::NavigatorBase::cd(), hchi2, hChi2VsEta, hChi2VsPhi, heta, hnHits, hnHitsVsEta, hnHitsVsPhi, hnLayers, hnLayersVsEta, hnLayersVsPhi, hnLooseAndAboveTracks, hnTracks, hphi, hpt, hquality, htip, hz, mps_fire::i, M_PI, submitPVResolutionJobs::q, pixelTrack::qualityName, dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, topFolderName_, and parallelization::uint().

                                                                                {
   iBook.cd();
   iBook.setCurrentFolder(topFolderName_);
 
   // clang-format off
   std::string toRep = "Number of tracks";
   hnTracks = iBook.book1D("nTracks", fmt::sprintf(";%s per event;#events",toRep), 1001, -0.5, 1000.5);
   hnLooseAndAboveTracks = iBook.book1D("nLooseAndAboveTracks", fmt::sprintf(";%s (quality #geq loose) per event;#events",toRep), 1001, -0.5, 1000.5);
 
   toRep = "Number of all RecHits per track (quality #geq loose)";
   hnHits = iBook.book1D("nRecHits", fmt::sprintf(";%s;#tracks",toRep), 15, -0.5, 14.5);
   hnHitsVsPhi = iBook.bookProfile("nHitsPerTrackVsPhi", fmt::sprintf("%s vs track #eta;Track #eta;%s",toRep,toRep), 30, -3., 3., 0., 15.);
   hnHitsVsEta = iBook.bookProfile("nHitsPerTrackVsEta", fmt::sprintf("%s vs track #phi;Track #phi;%s",toRep,toRep), 30,-M_PI,-M_PI,0.,15.);
 
   toRep = "Number of all layers per track (quality #geq loose)";
   hnLayers = iBook.book1D("nLayers", fmt::sprintf(";%s;#tracks",toRep), 15, -0.5, 14.5);
   hnLayersVsPhi = iBook.bookProfile("nLayersPerTrackVsPhi", fmt::sprintf("%s vs track #eta;Track #eta;%s",toRep,toRep), 30, -3., 3., 0., 15.);
   hnLayersVsEta = iBook.bookProfile("nLayersPerTrackVsEta", fmt::sprintf("%s vs track #phi;Track #phi;%s",toRep,toRep), 30,-M_PI,-M_PI,0.,15.);
 
   toRep = "Track (quality #geq loose) #chi^{2}/ndof";
   hchi2 = iBook.book1D("nChi2ndof", fmt::sprintf(";%s;#tracks",toRep), 40, 0., 20.);
   hChi2VsPhi = iBook.bookProfile("nChi2ndofVsPhi", fmt::sprintf("%s vs track #eta;Track #eta;%s",toRep,toRep), 30, -3., 3., 0., 20.);
   hChi2VsEta = iBook.bookProfile("nChi2ndofVsEta", fmt::sprintf("%s vs track #phi;Track #phi;%s",toRep,toRep), 30, -M_PI, -M_PI, 0., 20.);
   // clang-format on
 
   hpt = iBook.book1D("pt", ";Track (quality #geq loose) p_{T} [GeV];#tracks", 200, 0., 200.);
   heta = iBook.book1D("eta", ";Track (quality #geq loose) #eta;#tracks", 30, -3., 3.);
   hphi = iBook.book1D("phi", ";Track (quality #geq loose) #phi;#tracks", 30, -M_PI, M_PI);
   hz = iBook.book1D("z", ";Track (quality #geq loose) z [cm];#tracks", 30, -30., 30.);
   htip = iBook.book1D("tip", ";Track (quality #geq loose) TIP [cm];#tracks", 100, -0.5, 0.5);
   hquality = iBook.book1D("quality", ";Track Quality;#tracks", 7, -0.5, 6.5);
   uint i = 1;
   for (const auto& q : pixelTrack::qualityName) {
     hquality->setBinLabel(i, q, 1);
     i++;
   }
 }

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

static

Definition at line 174 of file SiPixelPhase1MonitorTrackSoA.cc.

References edm::ParameterSetDescription::add(), edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, HLT_FULL_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                               {
   // monitorpixelTrackSoA
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("pixelTrackSrc", edm::InputTag("pixelTracksSoA"));
   desc.add<std::string>("topFolderName", "SiPixelHeterogeneous/PixelTrackSoA");
   desc.add<bool>("useQualityCut", true);
   desc.add<std::string>("minQuality", "loose");
   descriptions.addWithDefaultLabel(desc);
 }