#include <MonitorTrackResiduals.h>

Inheritance diagram for MonitorTrackResidualsBase< pixel_or_strip >:

Classes
struct	HistoPair

struct	HistoXY

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

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

void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override

	MonitorTrackResidualsBase (const edm::ParameterSet &)

	~MonitorTrackResidualsBase () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Private Types
typedef std::map< std::pair< std::string, int32_t >, HistoXY >	HistoSet

Private Member Functions
void	createMEs (DQMStore::IBooker &, const edm::EventSetup &)

std::pair< std::string, int32_t >	findSubdetAndLayer (uint32_t ModuleID, const TrackerTopology *tTopo)

Private Attributes
TrackerValidationVariables	avalidator_

edm::ParameterSet	conf_

GenericTriggerEventFlag *	genTriggerEventFlag_

unsigned long long	m_cacheID_

HistoSet	m_ModuleResiduals

HistoSet	m_SubdetLayerResiduals

bool	ModOn

edm::EDGetTokenT< reco::VertexCollection >	offlinePrimaryVerticesToken_

edm::ParameterSet	Parameters

std::unique_ptr< TkHistoMap >	tkhisto_ResidualsMean

Detailed Description

template<TrackerType pixel_or_strip>
class MonitorTrackResidualsBase< pixel_or_strip >

Definition at line 39 of file MonitorTrackResiduals.h.

Member Typedef Documentation

template<TrackerType pixel_or_strip>

typedef std::map<std::pair<std::string, int32_t>, HistoXY> MonitorTrackResidualsBase< pixel_or_strip >::HistoSet

private

Definition at line 65 of file MonitorTrackResiduals.h.

Constructor & Destructor Documentation

template<TrackerType pixel_or_strip>

MonitorTrackResidualsBase< pixel_or_strip >::MonitorTrackResidualsBase ( const edm::ParameterSet & iConfig )

explicit

Definition at line 20 of file MonitorTrackResiduals.cc.

References MonitorTrackResidualsBase< pixel_or_strip >::conf_, edm::ParameterSet::getParameter(), MonitorTrackResidualsBase< pixel_or_strip >::ModOn, MonitorTrackResidualsBase< pixel_or_strip >::offlinePrimaryVerticesToken_, and AlCaHLTBitMon_QueryRunRegistry::string.

     : conf_(iConfig),
       m_cacheID_(0),
       genTriggerEventFlag_(new GenericTriggerEventFlag(
           iConfig.getParameter<edm::ParameterSet>("genericTriggerEventPSet"), consumesCollector(), *this)),
       avalidator_(iConfig, consumesCollector()) {
   ModOn = conf_.getParameter<bool>("Mod_On");
   offlinePrimaryVerticesToken_ = consumes<reco::VertexCollection>(std::string("offlinePrimaryVertices"));
 }

template<TrackerType pixel_or_strip>

MonitorTrackResidualsBase< pixel_or_strip >::~MonitorTrackResidualsBase ( )

override

Definition at line 31 of file MonitorTrackResiduals.cc.

References MonitorTrackResidualsBase< pixel_or_strip >::genTriggerEventFlag_.

                                                                       {
   if (genTriggerEventFlag_)
     delete genTriggerEventFlag_;
 }

Member Function Documentation

template<TrackerType pixel_or_strip>

void MonitorTrackResidualsBase< pixel_or_strip >::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 237 of file MonitorTrackResiduals.cc.

                                                                                                            {
   auto vtracks = std::vector<TrackerValidationVariables::AVTrackStruct>();
   // Filter out events if Trigger Filtering is requested
   if (genTriggerEventFlag_->on() && !genTriggerEventFlag_->accept(iEvent, iSetup))
     return;
 
   edm::Handle<reco::VertexCollection> vertices;
   iEvent.getByToken(offlinePrimaryVerticesToken_, vertices);
   if (!vertices.isValid() || vertices->empty())
     return;
   const auto primaryVertex = vertices->at(0);
 
   // Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology *const tTopo = tTopoHandle.product();
 
   avalidator_.fillTrackQuantities(iEvent,
                                   iSetup,
                                   // tell the validator to only look at good tracks
                                   [&](const reco::Track &track) -> bool {
                                     return track.pt() > 0.75 &&
                                            abs(track.dxy(primaryVertex.position())) < 5 * track.dxyError();
                                   },
                                   vtracks);
 
   for (auto &track : vtracks) {
     for (auto &it : track.hits) {
       uint RawId = it.rawDetId;
 
       auto id = DetId(RawId);
       auto isPixel = id.subdetId() == 1 || id.subdetId() == 2;
       if (isPixel != (pixel_or_strip == TRACKERTYPE_PIXEL))
         continue;
 
       if (ModOn) {
         auto &mod_histos = m_ModuleResiduals[std::make_pair("", RawId)];
         mod_histos.x.base->Fill(it.resXprime);
         mod_histos.x.normed->Fill(it.resXprime / it.resXprimeErr);
         mod_histos.y.base->Fill(it.resYprime);
         mod_histos.y.normed->Fill(it.resYprime / it.resYprimeErr);
       }
 
       auto subdetandlayer = findSubdetAndLayer(RawId, tTopo);
       auto histos = m_SubdetLayerResiduals[subdetandlayer];
       // fill if its error is not zero
       if (it.resXprimeErr != 0 && histos.x.base) {
         histos.x.base->Fill(it.resXprime);
         histos.x.normed->Fill(it.resXprime / it.resXprimeErr);
         if (!isPixel)
           tkhisto_ResidualsMean->fill(RawId, it.resXprime);
       }
       if (it.resYprimeErr != 0 && histos.y.base) {
         histos.y.base->Fill(it.resYprime);
         histos.y.normed->Fill(it.resYprime / it.resYprimeErr);
       }
     }
   }
 }

template<TrackerType pixel_or_strip>

void MonitorTrackResidualsBase< pixel_or_strip >::bookHistograms	(	DQMStore::IBooker &	,
		edm::Run const &	,
		edm::EventSetup const &
	)

override

Definition at line 37 of file MonitorTrackResiduals.cc.

References MonitorTrackResidualsBase< pixel_or_strip >::createMEs(), edm::EventSetup::get(), MonitorTrackResidualsBase< pixel_or_strip >::m_cacheID_, edm::ESHandle< T >::product(), SiStripFolderOrganizer::setSiStripFolderName(), AlCaHLTBitMon_QueryRunRegistry::string, and MonitorTrackResidualsBase< pixel_or_strip >::tkhisto_ResidualsMean.

                                                                                             {
   unsigned long long cacheID = iSetup.get<TrackerDigiGeometryRecord>().cacheIdentifier();
   if (m_cacheID_ != cacheID) {
     m_cacheID_ = cacheID;
     this->createMEs(ibooker, iSetup);
   }
   std::string topFolderName_ = "SiStrip";
   SiStripFolderOrganizer folder_organizer;
   folder_organizer.setSiStripFolderName(topFolderName_);
   edm::ESHandle<TkDetMap> tkDetMapHandle;
   iSetup.get<TrackerTopologyRcd>().get(tkDetMapHandle);
   tkhisto_ResidualsMean = std::make_unique<TkHistoMap>(
       tkDetMapHandle.product(), ibooker, topFolderName_, "TkHMap_ResidualsMean", 0.0, true);
 }

template<TrackerType pixel_or_strip>

void MonitorTrackResidualsBase< pixel_or_strip >::createMEs	(	DQMStore::IBooker &	ibooker,
		const edm::EventSetup &	iSetup
	)

private

Definition at line 103 of file MonitorTrackResiduals.cc.

References funct::abs(), DQMStore::IBooker::book1D(), MonitorTrackResidualsBase< pixel_or_strip >::conf_, TrackerGeometry::detIds(), MonitorTrackResidualsBase< pixel_or_strip >::findSubdetAndLayer(), edm::EventSetup::get(), edm::ParameterSet::getParameter(), plotFactory::histos, GeomDetEnumerators::isBarrel(), fastTrackerRecHitType::isPixel(), MonitorTrackResidualsBase< pixel_or_strip >::m_ModuleResiduals, MonitorTrackResidualsBase< pixel_or_strip >::m_SubdetLayerResiduals, MonitorTrackResidualsBase< pixel_or_strip >::ModOn, edm::ESHandle< T >::product(), MonitorElement::setAxisTitle(), DQMStore::IBooker::setCurrentFolder(), SiStripFolderOrganizer::setDetectorFolder(), SiStripFolderOrganizer::setLayerFolder(), AlCaHLTBitMon_QueryRunRegistry::string, TRACKERTYPE_PIXEL, and geometryCSVtoXML::xy.

Referenced by MonitorTrackResidualsBase< pixel_or_strip >::bookHistograms().

                                                                                                                {
   // Retrieve tracker topology and geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology *const tTopo = tTopoHandle.product();
 
   edm::ESHandle<TrackerGeometry> TG;
   iSetup.get<TrackerDigiGeometryRecord>().get(TG);
 
   Parameters = conf_.getParameter<edm::ParameterSet>("TH1ResModules");
   int32_t i_residuals_Nbins = Parameters.getParameter<int32_t>("Nbinx");
   double d_residual_xmin = Parameters.getParameter<double>("xmin");
   double d_residual_xmax = Parameters.getParameter<double>("xmax");
   Parameters = conf_.getParameter<edm::ParameterSet>("TH1NormResModules");
   int32_t i_normres_Nbins = Parameters.getParameter<int32_t>("Nbinx");
   double d_normres_xmin = Parameters.getParameter<double>("xmin");
   double d_normres_xmax = Parameters.getParameter<double>("xmax");
 
   // use SistripHistoId for producing histogram id (and title)
   SiStripHistoId hidmanager;
 
   SiStripFolderOrganizer strip_organizer;
   auto pixel_organizer = SiPixelFolderOrganizer(false);
 
   // Collect list of modules from Tracker Geometry
   // book histo per each detector module
   auto ids = TG->detIds();  // or detUnitIds?
   for (DetId id : ids) {
     auto ModuleID = id.rawId();
     auto isPixel = id.subdetId() == 1 || id.subdetId() == 2;
     if (isPixel != (pixel_or_strip == TRACKERTYPE_PIXEL))
       continue;
 
     // Book module histogramms?
     if (ModOn) {
       switch (id.subdetId()) {
         case 1:
           pixel_organizer.setModuleFolder(ibooker, ModuleID, 0);
           break;
         case 2:
           pixel_organizer.setModuleFolder(ibooker, ModuleID, 0);
           break;
         default:
           strip_organizer.setDetectorFolder(ModuleID, tTopo);
       }
       {
         // this sounds strip specific but also works for pixel
         std::string hid = hidmanager.createHistoId("HitResidualsX", "det", ModuleID);
         std::string normhid = hidmanager.createHistoId("NormalizedHitResidualsX", "det", ModuleID);
         auto &histos = m_ModuleResiduals[std::make_pair("", ModuleID)];
         histos.x.base = ibooker.book1D(hid, hid, i_residuals_Nbins, d_residual_xmin, d_residual_xmax);
         histos.x.base->setAxisTitle("(x_{pred} - x_{rec})' [cm]");
         histos.x.normed = ibooker.book1D(normhid, normhid, i_normres_Nbins, d_normres_xmin, d_normres_xmax);
         histos.x.normed->setAxisTitle("(x_{pred} - x_{rec})'/#sigma");
       }
       {
         std::string hid = hidmanager.createHistoId("HitResidualsY", "det", ModuleID);
         std::string normhid = hidmanager.createHistoId("NormalizedHitResidualsY", "det", ModuleID);
         auto &histos = m_ModuleResiduals[std::make_pair("", ModuleID)];
         histos.y.base = ibooker.book1D(hid, hid, i_residuals_Nbins, d_residual_xmin, d_residual_xmax);
         histos.y.base->setAxisTitle("(y_{pred} - y_{rec})' [cm]");
         histos.y.normed = ibooker.book1D(normhid, normhid, i_normres_Nbins, d_normres_xmin, d_normres_xmax);
         histos.y.normed->setAxisTitle("(y_{pred} - y_{rec})'/#sigma");
       }
     }
 
     auto subdetandlayer = findSubdetAndLayer(ModuleID, tTopo);
     if (m_SubdetLayerResiduals.find(subdetandlayer) == m_SubdetLayerResiduals.end()) {
       // add new histograms
       auto &histos = m_SubdetLayerResiduals[subdetandlayer];
       switch (id.subdetId()) {
         // Pixel Barrel, Endcap
         // We can't use the folder organizer here (SiPixelActionExecutor.cc#1638
         // does the same)
         case 1:
           ibooker.setCurrentFolder("Pixel/Barrel");
           break;
         case 2:
           ibooker.setCurrentFolder("Pixel/Endcap");
           break;
         // All strip
         default:
           strip_organizer.setLayerFolder(ModuleID, tTopo, subdetandlayer.second);
       }
 
       auto isBarrel = subdetandlayer.first.find("B") != std::string::npos;
 
       auto xy = std::vector<std::pair<HistoPair &, const char *>>{std::make_pair(std::ref(histos.x), "X"),
                                                                   std::make_pair(std::ref(histos.y), "Y")};
       for (auto &histopair : xy) {
         // book histogramms on layer level, check for barrel/pixel only for
         // correct labeling
 
         // Skip the Y plots for strips.
         if (!isPixel && histopair.second[0] == 'Y')
           continue;
 
         std::string histoname = isPixel ? (  // Pixel name
                                               Form("HitResiduals%s_%s%d",
                                                    histopair.second,
                                                    isBarrel ? "L" : (subdetandlayer.second > 0 ? "Dp" : "Dm"),
                                                    std::abs(subdetandlayer.second)))
                                         : (Form("HitResiduals_%s__%s__%d",  // Strip TODO: We use a
                                                                             // legacy name.
                                                 subdetandlayer.first.c_str(),
                                                 isBarrel ? "Layer" : "wheel",
                                                 std::abs(subdetandlayer.second)));
 
         std::string histotitle = Form("HitResiduals %s on %s%s full %s %d",
                                       histopair.second,
                                       subdetandlayer.first.c_str(),
                                       isBarrel ? "" : (subdetandlayer.second > 0 ? "+" : "-"),
                                       isBarrel ? "Layer" : (isPixel ? "Disk" : "Wheel"),
                                       std::abs(subdetandlayer.second));
 
         std::string normhistoname = Form("Normalized%s", histoname.c_str());
         std::string normhistotitle = Form("Normalized%s", histotitle.c_str());
 
         // std::cout << "##### Booking: " << ibooker.pwd() << " title " <<
         // histoname << std::endl;
 
         histopair.first.base =
             ibooker.book1D(histoname.c_str(), histotitle.c_str(), i_residuals_Nbins, d_residual_xmin, d_residual_xmax);
         histopair.first.base->setAxisTitle("(x_{pred} - x_{rec})' [cm]");
 
         histopair.first.normed = ibooker.book1D(
             normhistoname.c_str(), normhistotitle.c_str(), i_normres_Nbins, d_normres_xmin, d_normres_xmax);
         histopair.first.normed->setAxisTitle("(x_{pred} - x_{rec})'/#sigma");
       }
     }
   }  // end loop over activeDets
 }

template<TrackerType pixel_or_strip>

void MonitorTrackResidualsBase< pixel_or_strip >::dqmBeginRun	(	const edm::Run &	,
		const edm::EventSetup &
	)

override

Definition at line 55 of file MonitorTrackResiduals.cc.

References MonitorTrackResidualsBase< pixel_or_strip >::genTriggerEventFlag_, GenericTriggerEventFlag::initRun(), and GenericTriggerEventFlag::on().

                                                                                                           {
   // Initialize the GenericTriggerEventFlag
   if (genTriggerEventFlag_->on())
     genTriggerEventFlag_->initRun(run, iSetup);
 }

template<TrackerType pixel_or_strip>

std::pair< std::string, int32_t > MonitorTrackResidualsBase< pixel_or_strip >::findSubdetAndLayer	(	uint32_t	ModuleID,
		const TrackerTopology *	tTopo
	)

private

Definition at line 62 of file MonitorTrackResiduals.cc.

References TrackerTopology::pxbLayer(), TrackerTopology::pxfDisk(), TrackerTopology::pxfSide(), AlCaHLTBitMon_QueryRunRegistry::string, TrackerTopology::tecSide(), TrackerTopology::tecWheel(), TrackerTopology::tibLayer(), TrackerTopology::tidSide(), TrackerTopology::tidWheel(), and TrackerTopology::tobLayer().

Referenced by MonitorTrackResidualsBase< pixel_or_strip >::analyze(), and MonitorTrackResidualsBase< pixel_or_strip >::createMEs().

                                                      {
   std::string subdet = "";
   int32_t layer = 0;
   auto id = DetId(ModuleID);
   switch (id.subdetId()) {
     // Pixel Barrel, Endcap
     case 1:
       subdet = "BPIX";
       layer = tTopo->pxbLayer(id);
       break;
     case 2:
       subdet = "FPIX";
       layer = tTopo->pxfDisk(id) * (tTopo->pxfSide(ModuleID) == 1 ? -1 : +1);
       break;
     // Strip TIB, TID, TOB, TEC
     case 3:
       subdet = "TIB";
       layer = tTopo->tibLayer(id);
       break;
     case 4:
       subdet = "TID";
       layer = tTopo->tidWheel(id) * (tTopo->tidSide(ModuleID) == 1 ? -1 : +1);
       break;
     case 5:
       subdet = "TOB";
       layer = tTopo->tobLayer(id);
       break;
     case 6:
       subdet = "TEC";
       layer = tTopo->tecWheel(id) * (tTopo->tecSide(ModuleID) == 1 ? -1 : +1);
       break;
     default:
       // TODO: Fail loudly.
       subdet = "UNKNOWN";
       layer = 0;
   }
   return std::make_pair(subdet, layer);
 }