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

Inheritance diagram for MonitorTrackResiduals:

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 thread_unsafe::DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

	DQMEDAnalyzer (void)

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	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 (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
edm::ParameterSet	conf_

DQMStore *	dqmStore_

SiStripFolderOrganizer	folder_organizer

GenericTriggerEventFlag *	genTriggerEventFlag_

HistoClass	HitResidual

unsigned long long	m_cacheID_

std::map< std::pair < std::string, int32_t > , MonitorElement * >	m_SubdetLayerNormedResiduals

std::map< std::pair < std::string, int32_t > , MonitorElement * >	m_SubdetLayerResiduals

bool	ModOn

HistoClass	NormedHitResiduals

edm::ParameterSet	Parameters

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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.

Constructor & Destructor Documentation

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

explicit

Definition at line 24 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())) {
   ModOn = conf_.getParameter<bool>("Mod_On");
 }

MonitorTrackResiduals::~MonitorTrackResiduals ( )

Definition at line 31 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::EDAnalyzer.

Definition at line 150 of file MonitorTrackResiduals.cc.

References GenericTriggerEventFlag::accept(), conf_, TrackerValidationVariables::fillHitQuantities(), folder_organizer, genTriggerEventFlag_, edm::EventSetup::get(), SiStripFolderOrganizer::GetSubDetAndLayer(), HitResidual, m_SubdetLayerNormedResiduals, m_SubdetLayerResiduals, ModOn, NormedHitResiduals, GenericTriggerEventFlag::on(), edm::ESHandle< class >::product(), and SiStripDetId::subDetector().

                                                                                        {
 
   // 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<IdealGeometryRecord>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
 
   TrackerValidationVariables avalidator_(iSetup,conf_);
   std::vector<TrackerValidationVariables::AVHitStruct> v_hitstruct;
   avalidator_.fillHitQuantities(iEvent,v_hitstruct);
   for (std::vector<TrackerValidationVariables::AVHitStruct>::const_iterator it = v_hitstruct.begin(),
        itEnd = v_hitstruct.end(); it != itEnd; ++it) {
     uint RawId = it->rawDetId;
 
     // fill if hit belongs to StripDetector and its error is not zero
     if( it->resXprimeErr != 0 && SiStripDetId(RawId).subDetector()  != 0 )  {
       if (ModOn && HitResidual[RawId]) {
     HitResidual[RawId]->Fill(it->resXprime);
     NormedHitResiduals[RawId]->Fill(it->resXprime/it->resXprimeErr);
       }
       std::pair<std::string, int32_t> subdetandlayer = folder_organizer.GetSubDetAndLayer(RawId, tTopo);
       if(m_SubdetLayerResiduals[subdetandlayer]) {
     m_SubdetLayerResiduals[subdetandlayer]->Fill(it->resXprime);
     m_SubdetLayerNormedResiduals[subdetandlayer]->Fill(it->resXprime/it->resXprimeErr);
       }
     }
   }
 
 }

void MonitorTrackResiduals::beginJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 36 of file MonitorTrackResiduals.cc.

36 {

37 }

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

overridevirtual

Implements thread_unsafe::DQMEDAnalyzer.

Definition at line 39 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 54 of file MonitorTrackResiduals.cc.

References funct::abs(), DQMStore::IBooker::book1D(), conf_, folder_organizer, edm::EventSetup::get(), edm::ParameterSet::getParameter(), SiStripFolderOrganizer::GetSubDetAndLayer(), HitResidual, m_SubdetLayerNormedResiduals, m_SubdetLayerResiduals, ModOn, NormedHitResiduals, edm::ESHandle< class >::product(), MonitorElement::setAxisTitle(), SiStripFolderOrganizer::setDetectorFolder(), SiStripFolderOrganizer::setLayerFolder(), SiStripFolderOrganizer::setSiStripFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and SiStripDetId::subDetector().

Referenced by bookHistograms().

                                                                                              {
 
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   iSetup.get<IdealGeometryRecord>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
 
   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;
   folder_organizer.setSiStripFolder(); // top SiStrip folder
 
   // take from eventSetup the SiStripDetCabling object
   edm::ESHandle<SiStripDetCabling> tkmechstruct;
   iSetup.get<SiStripDetCablingRcd>().get(tkmechstruct);
 
   // get list of active detectors from SiStripDetCabling
   std::vector<uint32_t> activeDets;
   activeDets.clear(); // just in case
   tkmechstruct->addActiveDetectorsRawIds(activeDets);
 
   // use SiStripSubStructure for selecting certain regions
   SiStripSubStructure substructure;
   std::vector<uint32_t> DetIds = activeDets;
 
   // book histo per each detector module
   for (std::vector<uint32_t>::const_iterator DetItr=activeDets.begin(),
      DetItrEnd = activeDets.end(); DetItr!=DetItrEnd; ++DetItr)
     {
       uint ModuleID = (*DetItr);
 
       // is this a StripModule?
       if( SiStripDetId(ModuleID).subDetector() != 0 ) {
 
     folder_organizer.setDetectorFolder(ModuleID, tTopo); //  detid sets appropriate detector folder
     // Book module histogramms?
     if (ModOn) {
       std::string hid = hidmanager.createHistoId("HitResiduals","det",ModuleID);
       std::string normhid = hidmanager.createHistoId("NormalizedHitResiduals","det",ModuleID);
       HitResidual[ModuleID] = ibooker.book1D(hid, hid,
                          i_residuals_Nbins,d_residual_xmin,d_residual_xmax);
       HitResidual[ModuleID]->setAxisTitle("(x_{pred} - x_{rec})' [cm]");
       NormedHitResiduals[ModuleID] = ibooker.book1D(normhid, normhid,
                             i_normres_Nbins,d_normres_xmin,d_normres_xmax);
       NormedHitResiduals[ModuleID]->setAxisTitle("(x_{pred} - x_{rec})'/#sigma");
     }
     // book layer level histogramms
     std::pair<std::string,int32_t> subdetandlayer = folder_organizer.GetSubDetAndLayer(ModuleID, tTopo);
     folder_organizer.setLayerFolder(ModuleID,tTopo,subdetandlayer.second);
     if(! m_SubdetLayerResiduals[subdetandlayer ] ) {
       // book histogramms on layer level, check for barrel for correct labeling
       std::string histoname = Form(subdetandlayer.first.find("B") != std::string::npos ?
                        "HitResiduals_%s__Layer__%d" : "HitResiduals_%s__wheel__%d" ,
                        subdetandlayer.first.c_str(),std::abs(subdetandlayer.second));
       std::string normhistoname =
         Form(subdetandlayer.first.find("B") != std::string::npos ?
          "NormalizedHitResidual_%s__Layer__%d" : "NormalizedHitResidual_%s__wheel__%d" ,
          subdetandlayer.first.c_str(),std::abs(subdetandlayer.second));
       m_SubdetLayerResiduals[subdetandlayer] =
         ibooker.book1D(histoname.c_str(),histoname.c_str(),
                i_residuals_Nbins,d_residual_xmin,d_residual_xmax);
       m_SubdetLayerResiduals[subdetandlayer]->setAxisTitle("(x_{pred} - x_{rec})' [cm]");
       m_SubdetLayerNormedResiduals[subdetandlayer] =
         ibooker.book1D(normhistoname.c_str(),normhistoname.c_str(),
                i_normres_Nbins,d_normres_xmin,d_normres_xmax);
       m_SubdetLayerNormedResiduals[subdetandlayer]->setAxisTitle("(x_{pred} - x_{rec})'/#sigma");
     }
       } // end 'is strip module'
     } // end loop over activeDets
 }

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

virtual

Reimplemented from thread_unsafe::DQMEDAnalyzer.

Definition at line 48 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

Reimplemented from edm::EDAnalyzer.

Definition at line 139 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::EDAnalyzer.

Definition at line 135 of file MonitorTrackResiduals.cc.

135 {

136 }

Member Data Documentation

edm::ParameterSet MonitorTrackResiduals::conf_

private

Definition at line 47 of file MonitorTrackResiduals.h.

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

DQMStore* MonitorTrackResiduals::dqmStore_

private

Definition at line 46 of file MonitorTrackResiduals.h.

Referenced by endJob().

SiStripFolderOrganizer MonitorTrackResiduals::folder_organizer

private

Definition at line 53 of file MonitorTrackResiduals.h.

Referenced by analyze(), and createMEs().

GenericTriggerEventFlag* MonitorTrackResiduals::genTriggerEventFlag_

private

Definition at line 56 of file MonitorTrackResiduals.h.

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

HistoClass MonitorTrackResiduals::HitResidual

private

Definition at line 51 of file MonitorTrackResiduals.h.

Referenced by analyze(), and createMEs().

unsigned long long MonitorTrackResiduals::m_cacheID_

private

Definition at line 54 of file MonitorTrackResiduals.h.

Referenced by bookHistograms().

std::map< std::pair<std::string,int32_t>, MonitorElement* > MonitorTrackResiduals::m_SubdetLayerNormedResiduals

private

Definition at line 50 of file MonitorTrackResiduals.h.

Referenced by analyze(), and createMEs().

std::map< std::pair<std::string,int32_t>, MonitorElement* > MonitorTrackResiduals::m_SubdetLayerResiduals

private

Definition at line 49 of file MonitorTrackResiduals.h.

Referenced by analyze(), and createMEs().

bool MonitorTrackResiduals::ModOn

private

Definition at line 55 of file MonitorTrackResiduals.h.

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

HistoClass MonitorTrackResiduals::NormedHitResiduals

private

Definition at line 52 of file MonitorTrackResiduals.h.

Referenced by analyze(), and createMEs().

edm::ParameterSet MonitorTrackResiduals::Parameters

private

Definition at line 48 of file MonitorTrackResiduals.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation