#include <DQM/TrackerMonitorTrack/interface/MonitorTrackResiduals.cc>

Inheritance diagram for MonitorTrackResiduals:

Classes
struct	HistoPair

struct	HistoXY

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

virtual void	beginJob (void)

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

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

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

virtual void	endJob (void)

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

	MonitorTrackResiduals (const edm::ParameterSet &)

	~MonitorTrackResiduals ()

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

virtual void	beginStream (edm::StreamID id) final

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

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

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

Private Attributes
TrackerValidationVariables	avalidator_

edm::ParameterSet	conf_

DQMStore *	dqmStore_

GenericTriggerEventFlag *	genTriggerEventFlag_

unsigned long long	m_cacheID_

HistoSet	m_ModuleResiduals

HistoSet	m_SubdetLayerResiduals

bool	ModOn

edm::ParameterSet	Parameters

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

Monitoring source for track residuals on each detector module

Definition at line 31 of file MonitorTrackResiduals.h.

Member Typedef Documentation

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

private

Definition at line 61 of file MonitorTrackResiduals.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 16 of file MonitorTrackResiduals.cc.

References conf_, edm::ParameterSet::getParameter(), and ModOn.

    : dqmStore_( edm::Service<DQMStore>().operator->() )
    , conf_(iConfig), m_cacheID_(0)
    , genTriggerEventFlag_(new GenericTriggerEventFlag(iConfig, consumesCollector(), *this))
    , avalidator_(iConfig, consumesCollector()) {
   ModOn = conf_.getParameter<bool>("Mod_On");
 }

MonitorTrackResiduals::~MonitorTrackResiduals ( )

Definition at line 24 of file MonitorTrackResiduals.cc.

References genTriggerEventFlag_.

                                               {
   if (genTriggerEventFlag_) delete genTriggerEventFlag_;
 }

Member Function Documentation

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

virtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 222 of file MonitorTrackResiduals.cc.

References GenericTriggerEventFlag::accept(), avalidator_, TrackerValidationVariables::fillHitQuantities(), findSubdetAndLayer(), genTriggerEventFlag_, edm::EventSetup::get(), mergeVDriftHistosByStation::histos, m_ModuleResiduals, m_SubdetLayerResiduals, ModOn, GenericTriggerEventFlag::on(), and edm::ESHandle< class >::product().

                                                                                        {
 
   // Filter out events if Trigger Filtering is requested
   if (genTriggerEventFlag_->on()&& ! genTriggerEventFlag_->accept( iEvent, iSetup) ) return;
 
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
 
   std::vector<TrackerValidationVariables::AVHitStruct> v_hitstruct;
   avalidator_.fillHitQuantities(iEvent,v_hitstruct);
   for (auto it : v_hitstruct) {
     uint RawId = it.rawDetId;
 
     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(it.resYprimeErr != 0 && histos.y.base) {
       histos.y.base->Fill(it.resYprime);
       histos.y.normed->Fill(it.resYprime/it.resYprimeErr);
     }
   }
 }

void MonitorTrackResiduals::beginJob ( void )

virtual

Definition at line 29 of file MonitorTrackResiduals.cc.

29 {

30 }

void MonitorTrackResiduals::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	run,
		edm::EventSetup const &	iSetup
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 32 of file MonitorTrackResiduals.cc.

References createMEs(), edm::EventSetup::get(), and m_cacheID_.

 {
   unsigned long long cacheID = iSetup.get<SiStripDetCablingRcd>().cacheIdentifier();
   if (m_cacheID_ != cacheID) {
     m_cacheID_ = cacheID;
     this->createMEs( ibooker , iSetup);
   }
 }

void MonitorTrackResiduals::createMEs	(	DQMStore::IBooker &	ibooker,
		const edm::EventSetup &	iSetup
	)

Definition at line 87 of file MonitorTrackResiduals.cc.

References funct::abs(), DQMStore::IBooker::book1D(), conf_, findSubdetAndLayer(), edm::EventSetup::get(), edm::ParameterSet::getParameter(), mergeVDriftHistosByStation::histos, GeomDetEnumerators::isBarrel(), fastTrackerRecHitType::isPixel(), m_ModuleResiduals, m_SubdetLayerResiduals, ModOn, edm::ESHandle< class >::product(), MonitorElement::setAxisTitle(), SiStripFolderOrganizer::setDetectorFolder(), SiStripFolderOrganizer::setLayerFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and create_public_lumi_plots::xy.

Referenced by 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
   std::vector<uint32_t> activeDets;
   auto ids = TG->detIds(); // or detUnitIds?
   for (DetId id : ids) {
     activeDets.push_back(id.rawId());
   }
 
   // book histo per each detector module
   for(auto ModuleID : activeDets)
     {
       auto id = DetId(ModuleID);
 
       // TODO: Not yet implemented for Pixel.
       // 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
       case 1:   pixel_organizer.setModuleFolder(ibooker, ModuleID, 2); break;
       case 2:   pixel_organizer.setModuleFolder(ibooker, ModuleID, 6); break;
       // All strip
       default:  strip_organizer.setLayerFolder(ModuleID,tTopo,subdetandlayer.second);
     }
     
     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
       auto isBarrel = subdetandlayer.first.find("B") != std::string::npos;
       auto isPixel = subdetandlayer.first.find("X") != std::string::npos;
       // TODO: We use a legacy name to stay compatible with other code. 
       // Check if this is necessary.
       std::string histoname = Form("HitResiduals_%s__%s__%d%s",
         subdetandlayer.first.c_str(),
         isBarrel ? "Layer" : "wheel",
         std::abs(subdetandlayer.second),
         histopair.second[0] == 'X' ? "" : "_Y");
       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
 }

void MonitorTrackResiduals::dqmBeginRun	(	const edm::Run &	run,
		const edm::EventSetup &	iSetup
	)

virtual

Reimplemented from DQMEDAnalyzer.

Definition at line 41 of file MonitorTrackResiduals.cc.

References genTriggerEventFlag_, GenericTriggerEventFlag::initRun(), and GenericTriggerEventFlag::on().

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

void MonitorTrackResiduals::endJob ( void )

virtual

Definition at line 211 of file MonitorTrackResiduals.cc.

References conf_, dqmStore_, edm::ParameterSet::getParameter(), dumpDBToFile_GT_ttrig_cfg::outputFileName, DQMStore::save(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                       {
 
   //dqmStore_->showDirStructure();
   bool outputMEsInRootFile = conf_.getParameter<bool>("OutputMEsInRootFile");
   std::string outputFileName = conf_.getParameter<std::string>("OutputFileName");
   if(outputMEsInRootFile){
     dqmStore_->save(outputFileName);
   }
 }

void MonitorTrackResiduals::endRun	(	const edm::Run &	,
		const edm::EventSetup &
	)

virtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 207 of file MonitorTrackResiduals.cc.

207 {

208 }

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

private

Definition at line 47 of file MonitorTrackResiduals.cc.

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

Referenced by analyze(), and 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->tecSide(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);
 }

Member Data Documentation

TrackerValidationVariables MonitorTrackResiduals::avalidator_

private

Definition at line 69 of file MonitorTrackResiduals.h.

Referenced by analyze().

edm::ParameterSet MonitorTrackResiduals::conf_

private

Definition at line 47 of file MonitorTrackResiduals.h.

Referenced by createMEs(), endJob(), and MonitorTrackResiduals().

DQMStore* MonitorTrackResiduals::dqmStore_

private

Definition at line 46 of file MonitorTrackResiduals.h.

Referenced by endJob().

GenericTriggerEventFlag* MonitorTrackResiduals::genTriggerEventFlag_

private

Definition at line 68 of file MonitorTrackResiduals.h.

Referenced by analyze(), dqmBeginRun(), and ~MonitorTrackResiduals().

unsigned long long MonitorTrackResiduals::m_cacheID_

private

Definition at line 66 of file MonitorTrackResiduals.h.

Referenced by bookHistograms().

HistoSet MonitorTrackResiduals::m_ModuleResiduals

private

Definition at line 64 of file MonitorTrackResiduals.h.

Referenced by analyze(), and createMEs().

HistoSet MonitorTrackResiduals::m_SubdetLayerResiduals

private

Definition at line 63 of file MonitorTrackResiduals.h.

Referenced by analyze(), and createMEs().

bool MonitorTrackResiduals::ModOn

private

Definition at line 67 of file MonitorTrackResiduals.h.

Referenced by analyze(), createMEs(), and MonitorTrackResiduals().

edm::ParameterSet MonitorTrackResiduals::Parameters

private

Definition at line 48 of file MonitorTrackResiduals.h.

Classes

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation