#include <MonitorTrackResiduals.h>

Inheritance diagram for MonitorTrackResidualsBase< pixel_or_strip >:

Classes
struct	HistoPair

struct	HistoXY

Public Member Functions
virtual 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 ()

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr< dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

template<TrackerType pixel_or_strip>
class MonitorTrackResidualsBase< pixel_or_strip >

Definition at line 38 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 61 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 ( )

Definition at line 29 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
	)

overridevirtual

Definition at line 220 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->size() == 0) 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 &
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 34 of file MonitorTrackResiduals.cc.

References MonitorTrackResidualsBase< pixel_or_strip >::createMEs(), edm::EventSetup::get(), MonitorTrackResidualsBase< pixel_or_strip >::m_cacheID_, 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_);
   tkhisto_ResidualsMean = std::unique_ptr<TkHistoMap>(new TkHistoMap(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 96 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, gedGsfElectrons_cfi::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 48 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 55 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);
 }