#include <SiStripMonitorTrack.h>

Inheritance diagram for SiStripMonitorTrack:

Classes
struct	Det2MEs

struct	LayerMEs

struct	ModMEs

struct	RingMEs

struct	SubDetMEs

Public Types
typedef TrackingRecHit::ConstRecHitPointer	ConstRecHitPointer

enum	RecHitType { Single = 0, Matched = 1, Projected = 2, Null = 3 }

Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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 &run, const edm::EventSetup &es) override

	SiStripMonitorTrack (const edm::ParameterSet &)

	~SiStripMonitorTrack () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

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

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

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

	DQMEDAnalyzer ()

virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

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

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

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

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Types
enum	ClusterFlags { OffTrack, OnTrack }

Private Member Functions
void	AllClusters (const edm::Event &ev, const edm::EventSetup &es)

void	book (DQMStore::IBooker &, const TrackerTopology tTopo, const TkDetMap tkDetMap)

void	bookLayerMEs (DQMStore::IBooker &, const uint32_t, std::string &)

MonitorElement *	bookME1D (DQMStore::IBooker &, const char , const char )

MonitorElement *	bookME2D (DQMStore::IBooker &, const char , const char )

MonitorElement *	bookME3D (DQMStore::IBooker &, const char , const char )

MonitorElement *	bookMEProfile (DQMStore::IBooker &, const char , const char )

MonitorElement *	bookMETrend (DQMStore::IBooker &, const char *)

void	bookModMEs (DQMStore::IBooker &, const uint32_t)

void	bookRing (DQMStore::IBooker &, const uint32_t, std::string &)

void	bookRingMEs (DQMStore::IBooker &, const uint32_t, std::string &)

void	bookSubDetMEs (DQMStore::IBooker &, std::string &name)

bool	clusterInfos (SiStripClusterInfo cluster, const uint32_t detid, enum ClusterFlags flags, bool track_ok, LocalVector LV, const Det2MEs &MEs, const TrackerTopology tTopo, const SiStripGain stripGain, const SiStripQuality stripQuality, const edm::DetSetVector< SiStripDigi > &digilist, float clustZ, float clustPhi)

bool	fillControlViewHistos (const edm::Event &ev, const edm::EventSetup &es)

void	fillME (MonitorElement *ME, float value1)

void	fillME (MonitorElement *ME, float value1, float value2)

void	fillME (MonitorElement *ME, float value1, float value2, float value3)

void	fillME (MonitorElement *ME, float value1, float value2, float value3, float value4)

Det2MEs	findMEs (const TrackerTopology *tTopo, const uint32_t detid)

MonitorElement *	handleBookMEs (DQMStore::IBooker &, std::string &, std::string &, std::string &, std::string &)

void	hitStudy (const edm::Event &ev, const edm::EventSetup &es, const edm::DetSetVector< SiStripDigi > &digilist, const ProjectedSiStripRecHit2D projhit, const SiStripMatchedRecHit2D matchedhit, const SiStripRecHit2D hit2D, const SiStripRecHit1D hit1D, LocalVector localMomentum, const bool track_ok)

template<class T >
void	RecHitInfo (const T *tkrecHit, LocalVector LV, const edm::DetSetVector< SiStripDigi > &digilist, const edm::Event &ev, const edm::EventSetup &es, bool track_ok)

void	return2DME (MonitorElement input1, MonitorElement input2, int binx, int biny, double value)

bool	trackFilter (const reco::Track &track)

void	trackStudy (const edm::Event &ev, const edm::EventSetup &es)

void	trackStudyFromTrack (edm::Handle< reco::TrackCollection > trackCollectionHandle, const edm::DetSetVector< SiStripDigi > &digilist, const edm::Event &ev, const edm::EventSetup &es)

void	trackStudyFromTrajectory (edm::Handle< reco::TrackCollection > trackCollectionHandle, const edm::DetSetVector< SiStripDigi > &digilist, const edm::Event &ev, const edm::EventSetup &es)

void	trajectoryStudy (const reco::Track &track, const edm::DetSetVector< SiStripDigi > &digilist, const edm::Event &ev, const edm::EventSetup &es, bool track_ok)

Private Attributes
bool	applyClusterQuality_

bool	clchCMoriginTkHmap_On_

edm::InputTag	Cluster_src_

MonitorElement *	ClusterCount_OnTrack_FECCratevsFECSlot = 0

MonitorElement *	ClusterCount_OnTrack_FECSlotVsFECRing_TECM = 0

MonitorElement *	ClusterCount_OnTrack_FECSlotVsFECRing_TECP = 0

MonitorElement *	ClusterCount_OnTrack_FECSlotVsFECRing_TIBTID = 0

MonitorElement *	ClusterCount_OnTrack_FECSlotVsFECRing_TOB = 0

MonitorElement *	ClusterStoNCorr_OnTrack_FECCratevsFECSlot = 0

MonitorElement *	ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TECM = 0

MonitorElement *	ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TECP = 0

MonitorElement *	ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TIBTID = 0

MonitorElement *	ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TOB = 0

MonitorElement *	ClusterStoNCorr_OnTrack_TECM = 0

MonitorElement *	ClusterStoNCorr_OnTrack_TECP = 0

MonitorElement *	ClusterStoNCorr_OnTrack_TIBTID = 0

MonitorElement *	ClusterStoNCorr_OnTrack_TOB = 0

edm::EDGetTokenT< edmNew::DetSetVector< SiStripCluster > >	clusterToken_

edm::ParameterSet	conf_

SiStripDCSStatus *	dcsStatus_

edm::EDGetTokenT< edm::DetSetVector< SiStripDigi > >	digiToken_

edm::EventNumber_t	eventNb

int	firstEvent

SiStripFolderOrganizer	folderOrganizer_

GenericTriggerEventFlag *	genTriggerEventFlag_

std::string	histname

bool	HistoFlag_On_

float	iLumisection

float	iOrbitSec

std::map< std::string, LayerMEs >	LayerMEsMap

LocalVector	LV

bool	Mod_On_

std::map< std::string, ModMEs >	ModMEsMap

int	numTracks

bool	OffHisto_On_

edm::ParameterSet	Parameters

bool	ring_flag

std::map< std::string, RingMEs >	RingMEsMap

edm::ESHandle< SiStripDetCabling >	SiStripDetCabling_

double	sToNLowerLimit_

double	sToNUpperLimit_

const SiStripGain *	stripGain_

const SiStripQuality *	stripQuality_

std::map< std::string, SubDetMEs >	SubDetMEsMap

edm::ESHandle< TrackerGeometry >	tkgeom_

std::unique_ptr< TkHistoMap >	tkhisto_ClChPerCMfromOrigin

std::unique_ptr< TkHistoMap >	tkhisto_ClChPerCMfromTrack

std::unique_ptr< TkHistoMap >	tkhisto_ClusterWidthOffTrack

std::unique_ptr< TkHistoMap >	tkhisto_ClusterWidthOnTrack

std::unique_ptr< TkHistoMap >	tkhisto_NoiseOffTrack

std::unique_ptr< TkHistoMap >	tkhisto_NoiseOnTrack

std::unique_ptr< TkHistoMap >	tkhisto_NumberInactiveHits

std::unique_ptr< TkHistoMap >	tkhisto_NumberValidHits

std::unique_ptr< TkHistoMap >	tkhisto_NumMissingHits

std::unique_ptr< TkHistoMap >	tkhisto_NumOffTrack

std::unique_ptr< TkHistoMap >	tkhisto_NumOnTrack

std::unique_ptr< TkHistoMap >	tkhisto_StoNCorrOnTrack

bool	TkHistoMap_On_

std::string	topFolderName_

bool	trackAssociatorCollection_in_EventTree

std::string	TrackLabel_

std::string	TrackProducer_

bool	tracksCollection_in_EventTree

edm::EDGetTokenT< reco::TrackCollection >	trackToken_

bool	Trend_On_

const TrackerTopology *	tTopo_

std::unordered_set< const SiStripCluster * >	vPSiStripCluster

edm::ESWatcher< SiStripGainRcd >	watcherStripGain_

edm::ESWatcher< TrackerTopologyRcd >	watchertTopo_

double	widthLowerLimit_

double	widthUpperLimit_

Additional Inherited Members
Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 65 of file SiStripMonitorTrack.h.

Member Typedef Documentation

typedef TrackingRecHit::ConstRecHitPointer SiStripMonitorTrack::ConstRecHitPointer

Definition at line 67 of file SiStripMonitorTrack.h.

Member Enumeration Documentation

enum SiStripMonitorTrack::ClusterFlags

private

Enumerator
OffTrack
OnTrack

Definition at line 76 of file SiStripMonitorTrack.h.

76 { OffTrack, OnTrack };

SiStripMonitorTrack::OffTrack

Definition: SiStripMonitorTrack.h:76

SiStripMonitorTrack::OnTrack

Definition: SiStripMonitorTrack.h:76

enum SiStripMonitorTrack::RecHitType

Enumerator
Single
Matched
Projected
Null

Definition at line 68 of file SiStripMonitorTrack.h.

68 { Single = 0, Matched = 1, Projected = 2, Null = 3 };

SiStripMonitorTrack::Matched

Definition: SiStripMonitorTrack.h:68

SiStripMonitorTrack::Null

Definition: SiStripMonitorTrack.h:68

SiStripMonitorTrack::Single

Definition: SiStripMonitorTrack.h:68

SiStripMonitorTrack::Projected

Definition: SiStripMonitorTrack.h:68

Constructor & Destructor Documentation

SiStripMonitorTrack::SiStripMonitorTrack ( const edm::ParameterSet & conf )

explicit

Definition at line 29 of file SiStripMonitorTrack.cc.

References applyClusterQuality_, clchCMoriginTkHmap_On_, Cluster_src_, clusterToken_, conf_, edm::EDConsumerBase::consumesCollector(), dcsStatus_, digiToken_, edm::ParameterSet::getParameter(), HLT_2018_cff::InputTag, Mod_On_, sToNLowerLimit_, sToNUpperLimit_, AlCaHLTBitMon_QueryRunRegistry::string, TkHistoMap_On_, topFolderName_, TrackLabel_, TrackProducer_, trackToken_, Trend_On_, widthLowerLimit_, and widthUpperLimit_.

     : conf_(conf),
       tracksCollection_in_EventTree(true),
       firstEvent(-1),
       genTriggerEventFlag_(new GenericTriggerEventFlag(
           conf.getParameter<edm::ParameterSet>("genericTriggerEventPSet"), consumesCollector(), *this)) {
   Cluster_src_ = conf.getParameter<edm::InputTag>("Cluster_src");
   Mod_On_ = conf.getParameter<bool>("Mod_On");
   Trend_On_ = conf.getParameter<bool>("Trend_On");
   TkHistoMap_On_ = conf.getParameter<bool>("TkHistoMap_On");
   clchCMoriginTkHmap_On_ = conf.getParameter<bool>("clchCMoriginTkHmap_On");
 
   TrackProducer_ = conf_.getParameter<std::string>("TrackProducer");
   TrackLabel_ = conf_.getParameter<std::string>("TrackLabel");
 
   topFolderName_ = conf_.getParameter<std::string>("TopFolderName");
 
   digiToken_ = consumes<edm::DetSetVector<SiStripDigi>>(conf.getParameter<edm::InputTag>("ADCDigi_src"));
   clusterToken_ = consumes<edmNew::DetSetVector<SiStripCluster>>(Cluster_src_);
   trackToken_ = consumes<reco::TrackCollection>(edm::InputTag(TrackProducer_, TrackLabel_));
 
   // cluster quality conditions
   edm::ParameterSet cluster_condition = conf_.getParameter<edm::ParameterSet>("ClusterConditions");
   applyClusterQuality_ = cluster_condition.getParameter<bool>("On");
   sToNLowerLimit_ = cluster_condition.getParameter<double>("minStoN");
   sToNUpperLimit_ = cluster_condition.getParameter<double>("maxStoN");
   widthLowerLimit_ = cluster_condition.getParameter<double>("minWidth");
   widthUpperLimit_ = cluster_condition.getParameter<double>("maxWidth");
 
   // Create DCS Status
   bool checkDCS = conf_.getParameter<bool>("UseDCSFiltering");
   if (checkDCS)
     dcsStatus_ = new SiStripDCSStatus(consumesCollector());
   else
     dcsStatus_ = nullptr;
 }

SiStripMonitorTrack::~SiStripMonitorTrack ( )

override

Definition at line 67 of file SiStripMonitorTrack.cc.

References dcsStatus_, and genTriggerEventFlag_.

                                           {
   if (dcsStatus_)
     delete dcsStatus_;
   if (genTriggerEventFlag_)
     delete genTriggerEventFlag_;
 }

Member Function Documentation

void SiStripMonitorTrack::AllClusters	(	const edm::Event &	ev,
		const edm::EventSetup &	es
	)

private

Definition at line 1281 of file SiStripMonitorTrack.cc.

References clusterInfos(), clusterToken_, digiToken_, findMEs(), edm::EventSetup::get(), edm::Event::getByToken(), edm::HandleBase::isValid(), LogDebug, EcalMEFormatter_cfi::MEs, OffTrack, edm::Handle< T >::product(), edm::ESHandle< T >::product(), and vPSiStripCluster.

Referenced by analyze().

                                                                                  {
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   es.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
 
   /*edm::ESHandle<TrackerGeometry> tracker;
   es.get<TrackerDigiGeometryRecord>().get(tracker);
   const TrackerGeometry *tkgeom = &(*tracker); */
 
   edm::ESHandle<SiStripGain> gainHandle;
   es.get<SiStripGainRcd>().get(gainHandle);
   const SiStripGain* stripGain = gainHandle.product();
 
   edm::ESHandle<SiStripQuality> qualityHandle;
   es.get<SiStripQualityRcd>().get("", qualityHandle);
   const SiStripQuality* stripQuality = qualityHandle.product();
 
   edm::Handle<edm::DetSetVector<SiStripDigi>> digihandle;
   ev.getByToken(digiToken_, digihandle);
 
   edm::DetSetVector<SiStripDigi> dummy;
   auto& digilist = digihandle.isValid() ? *(digihandle.product()) : dummy;
 
   edm::Handle<edmNew::DetSetVector<SiStripCluster>> siStripClusterHandle;
   ev.getByToken(clusterToken_, siStripClusterHandle);
 
   /*edm::ESHandle<StripClusterParameterEstimator> parameterestimator;
   es.get<TkStripCPERecord>().get("StripCPEfromTrackAngle", parameterestimator); 
   const StripClusterParameterEstimator &stripcpe(*parameterestimator);*/
 
   if (siStripClusterHandle.isValid()) {
     //Loop on Dets
     for (edmNew::DetSetVector<SiStripCluster>::const_iterator DSViter = siStripClusterHandle->begin(),
                                                               DSVEnd = siStripClusterHandle->end();
          DSViter != DSVEnd;
          ++DSViter) {
       uint32_t detid = DSViter->id();
       const Det2MEs MEs = findMEs(tTopo, detid);
 
       LogDebug("SiStripMonitorTrack") << "on detid " << detid << " N Cluster= " << DSViter->size();
 
       //Loop on Clusters
       for (edmNew::DetSet<SiStripCluster>::const_iterator ClusIter = DSViter->begin(), ClusEnd = DSViter->end();
            ClusIter != ClusEnd;
            ++ClusIter) {
         if (vPSiStripCluster.find(&*ClusIter) == vPSiStripCluster.end()) {
           SiStripClusterInfo SiStripClusterInfo_(*ClusIter, es, detid);
 
           /*const StripGeomDetUnit * stripdet = (const StripGeomDetUnit*) tkgeom->idToDetUnit(detid);
           StripClusterParameterEstimator::LocalValues parameters=stripcpe.localParameters(*ClusIter, *stripdet);
           LocalPoint lp = parameters.first;
       const Surface& surface = tracker->idToDet(detid)->surface();
       GlobalPoint gp = surface.toGlobal(lp);*/
 
           clusterInfos(&SiStripClusterInfo_,
                        detid,
                        OffTrack,
                        /*track_ok*/ false,
                        LV,
                        MEs,
                        tTopo,
                        stripGain,
                        stripQuality,
                        digilist,
                        0,
                        0);
         }
       }
     }
   } else {
     edm::LogError("SiStripMonitorTrack") << "ClusterCollection is not valid!!" << std::endl;
     return;
   }
 }

void SiStripMonitorTrack::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	es
	)

override

Definition at line 108 of file SiStripMonitorTrack.cc.

References GenericTriggerEventFlag::accept(), AllClusters(), edm::ESWatcher< T >::check(), dcsStatus_, edm::EventID::event(), eventNb, fillControlViewHistos(), fillME(), genTriggerEventFlag_, edm::EventSetup::get(), SiStripDCSStatus::getStatus(), edm::EventBase::id(), iLumisection, LogDebug, SiStripMonitorTrack::SubDetMEs::nClustersOffTrack, SiStripMonitorTrack::SubDetMEs::nClustersOnTrack, SiStripMonitorTrack::SubDetMEs::nClustersTrendOffTrack, SiStripMonitorTrack::SubDetMEs::nClustersTrendOnTrack, GenericTriggerEventFlag::on(), edm::EventBase::orbitNumber(), edm::ESHandle< T >::product(), edm::EventID::run(), stripGain_, stripQuality_, SubDetMEsMap, topFolderName_, SiStripMonitorTrack::SubDetMEs::totNClustersOffTrack, SiStripMonitorTrack::SubDetMEs::totNClustersOnTrack, trackStudy(), Trend_On_, tTopo_, vPSiStripCluster, watcherStripGain_, and watchertTopo_.

                                                                             {
   // Filter out events if DCS checking is requested
   if (dcsStatus_ && !dcsStatus_->getStatus(e, es))
     return;
 
   // Filter out events if Trigger Filtering is requested
   if (genTriggerEventFlag_->on() && !genTriggerEventFlag_->accept(e, es))
     return;
 
   //initialization of global quantities
   LogDebug("SiStripMonitorTrack") << "[SiStripMonitorTrack::analyse]  "
                                   << "Run " << e.id().run() << " Event " << e.id().event() << std::endl;
   eventNb = e.id().event();
   vPSiStripCluster.clear();
 
   iLumisection = e.orbitNumber() / 262144.0;
 
   // initialise # of clusters
   for (std::map<std::string, SubDetMEs>::iterator iSubDet = SubDetMEsMap.begin(); iSubDet != SubDetMEsMap.end();
        iSubDet++) {
     iSubDet->second.totNClustersOnTrack = 0;
     iSubDet->second.totNClustersOffTrack = 0;
   }
 
   if (watchertTopo_.check(es)) {
     edm::ESHandle<TrackerTopology> tTopoHandle;
     es.get<TrackerTopologyRcd>().get(tTopoHandle);
     tTopo_ = tTopoHandle.product();
   }
 
   if (watcherStripGain_.check(es)) {
     edm::ESHandle<SiStripGain> gainHandle;
     es.get<SiStripGainRcd>().get(gainHandle);
     stripGain_ = gainHandle.product();
   }
 
   edm::ESHandle<SiStripQuality> qualityHandle;
   es.get<SiStripQualityRcd>().get("", qualityHandle);
   stripQuality_ = qualityHandle.product();
 
   //Perform track study
   trackStudy(e, es);
 
   //Perform Cluster Study (irrespectively to tracks)
 
   AllClusters(e, es);  //analyzes the off Track Clusters
 
   //Summary Counts of clusters
   std::map<std::string, MonitorElement*>::iterator iME;
   std::map<std::string, LayerMEs>::iterator iLayerME;
 
   if (!(topFolderName_.find("IsolatedBunches") != std::string::npos)) {
     fillControlViewHistos(e, es);
   }
 
   if (Trend_On_) {
     // for (std::map<std::string, SubDetMEs>::iterator iSubDet = SubDetMEsMap.begin(), iterEnd=SubDetMEsMaps.end();
     //      iSubDet != iterEnd; ++iSubDet) {
     for (auto const& iSubDet : SubDetMEsMap) {
       SubDetMEs subdet_mes = iSubDet.second;
       if (subdet_mes.totNClustersOnTrack > 0) {
         fillME(subdet_mes.nClustersOnTrack, subdet_mes.totNClustersOnTrack);
       }
       fillME(subdet_mes.nClustersOffTrack, subdet_mes.totNClustersOffTrack);
       fillME(subdet_mes.nClustersTrendOnTrack, iLumisection, subdet_mes.totNClustersOnTrack);
       fillME(subdet_mes.nClustersTrendOffTrack, iLumisection, subdet_mes.totNClustersOffTrack);
     }
   } else {
     for (auto const& iSubDet : SubDetMEsMap) {
       SubDetMEs subdet_mes = iSubDet.second;
       if (subdet_mes.totNClustersOnTrack > 0) {
         fillME(subdet_mes.nClustersOnTrack, subdet_mes.totNClustersOnTrack);
       }
       fillME(subdet_mes.nClustersOffTrack, subdet_mes.totNClustersOffTrack);
     }
   }
 }

void SiStripMonitorTrack::book	(	DQMStore::IBooker &	ibooker,
		const TrackerTopology *	tTopo,
		const TkDetMap *	tkDetMap
	)

private

Definition at line 187 of file SiStripMonitorTrack.cc.

Referenced by bookHistograms().

                                                                                                                {
   SiStripFolderOrganizer folder_organizer;
   folder_organizer.setSiStripFolderName(topFolderName_);
   //******** TkHistoMaps
   if (TkHistoMap_On_) {
     tkhisto_StoNCorrOnTrack =
         std::make_unique<TkHistoMap>(tkDetMap, topFolderName_, "TkHMap_StoNCorrOnTrack", 0.0, true);
     tkhisto_NumOnTrack =
         std::make_unique<TkHistoMap>(tkDetMap, topFolderName_, "TkHMap_NumberOfOnTrackCluster", 0.0, true);
     tkhisto_NumOffTrack =
         std::make_unique<TkHistoMap>(tkDetMap, topFolderName_, "TkHMap_NumberOfOfffTrackCluster", 0.0, true);
     tkhisto_ClChPerCMfromTrack =
         std::make_unique<TkHistoMap>(tkDetMap, topFolderName_, "TkHMap_ChargePerCMfromTrack", 0.0, true);
     tkhisto_NumMissingHits =
         std::make_unique<TkHistoMap>(tkDetMap, topFolderName_, "TkHMap_NumberMissingHits", 0.0, true, true);
     tkhisto_NumberInactiveHits =
         std::make_unique<TkHistoMap>(tkDetMap, topFolderName_, "TkHMap_NumberInactiveHits", 0.0, true, true);
     tkhisto_NumberValidHits =
         std::make_unique<TkHistoMap>(tkDetMap, topFolderName_, "TkHMap_NumberValidHits", 0.0, true, true);
     tkhisto_NoiseOnTrack = std::make_unique<TkHistoMap>(tkDetMap, topFolderName_, "TkHMap_NoiseOnTrack", 0.0, true);
     tkhisto_NoiseOffTrack = std::make_unique<TkHistoMap>(tkDetMap, topFolderName_, "TkHMap_NoiseOffTrack", 0.0, true);
     tkhisto_ClusterWidthOnTrack =
         std::make_unique<TkHistoMap>(tkDetMap, topFolderName_, "TkHMap_ClusterWidthOnTrack", 0.0, true);
     tkhisto_ClusterWidthOffTrack =
         std::make_unique<TkHistoMap>(tkDetMap, topFolderName_, "TkHMap_ClusterWidthOffTrack", 0.0, true);
   }
   if (clchCMoriginTkHmap_On_)
     tkhisto_ClChPerCMfromOrigin =
         std::make_unique<TkHistoMap>(tkDetMap, topFolderName_, "TkHMap_ChargePerCMfromOrigin", 0.0, true);
   //******** TkHistoMaps
 
   std::vector<uint32_t> vdetId_;
   SiStripDetCabling_->addActiveDetectorsRawIds(vdetId_);
   const char* tec = "TEC";
   const char* tid = "TID";
   //Histos for each detector, layer and module
   SiStripHistoId hidmanager;
 
   if (Mod_On_) {
     for (std::vector<uint32_t>::const_iterator detid_iter = vdetId_.begin(), detid_end = vdetId_.end();
          detid_iter != detid_end;
          ++detid_iter) {  //loop on all the active detid
       uint32_t detid = *detid_iter;
 
       if (detid < 1) {
         edm::LogError("SiStripMonitorTrack") << "[" << __PRETTY_FUNCTION__ << "] invalid detid " << detid << std::endl;
         continue;
       }
 
       //std::string name;
 
       // book Layer and RING plots
       std::pair<std::string, int32_t> det_layer_pair = folder_organizer.GetSubDetAndLayer(detid, tTopo, false);
       /*
       std::string thickness;
       std::pair<std::string,int32_t> det_layer_pair_test = folder_organizer.GetSubDetAndLayerThickness(detid,tTopo,thickness);
       std::cout << "[SiStripMonitorTrack::book] det_layer_pair " << det_layer_pair.first << " " << det_layer_pair.second << " " << thickness << std::endl;
       */
 
       std::string layer_id = hidmanager.getSubdetid(detid, tTopo, false);
 
       std::map<std::string, LayerMEs>::iterator iLayerME = LayerMEsMap.find(layer_id);
       if (iLayerME == LayerMEsMap.end()) {
         folder_organizer.setLayerFolder(detid, tTopo, det_layer_pair.second, false);
         bookLayerMEs(ibooker, detid, layer_id);
       }
 
       const char* subdet = det_layer_pair.first.c_str();
       if (std::strstr(subdet, tec) != nullptr || std::strstr(subdet, tid) != nullptr) {
         std::string ring_id = hidmanager.getSubdetid(detid, tTopo, true);
         std::map<std::string, RingMEs>::iterator iRingME = RingMEsMap.find(ring_id);
         if (iRingME == RingMEsMap.end()) {
           std::pair<std::string, int32_t> det_ring_pair = folder_organizer.GetSubDetAndLayer(detid, tTopo, true);
           folder_organizer.setLayerFolder(detid, tTopo, det_ring_pair.second, true);
           bookRingMEs(ibooker, detid, ring_id);
         }
       }
 
       // book sub-detector plots
       std::pair<std::string, std::string> sdet_pair = folder_organizer.getSubDetFolderAndTag(detid, tTopo);
       if (SubDetMEsMap.find(sdet_pair.second) == SubDetMEsMap.end()) {
         ibooker.setCurrentFolder(sdet_pair.first);
         bookSubDetMEs(ibooker, sdet_pair.second);
       }
       // book module plots
       folder_organizer.setDetectorFolder(detid, tTopo);
       bookModMEs(ibooker, *detid_iter);
     }  //end loop on detectors detid
   } else {
     for (std::vector<uint32_t>::const_iterator detid_iter = vdetId_.begin(), detid_end = vdetId_.end();
          detid_iter != detid_end;
          ++detid_iter) {  //loop on all the active detid
       uint32_t detid = *detid_iter;
 
       if (detid < 1) {
         edm::LogError("SiStripMonitorTrack") << "[" << __PRETTY_FUNCTION__ << "] invalid detid " << detid << std::endl;
         continue;
       }
 
       //std::string name;
 
       // book Layer and RING plots
       std::pair<std::string, int32_t> det_layer_pair = folder_organizer.GetSubDetAndLayer(detid, tTopo, false);
       /*
       std::string thickness;
       std::pair<std::string,int32_t> det_layer_pair_test = folder_organizer.GetSubDetAndLayerThickness(detid,tTopo,thickness);
       std::cout << "[SiStripMonitorTrack::book] det_layer_pair " << det_layer_pair.first << " " << det_layer_pair.second << " " << thickness << std::endl;
       */
 
       std::string layer_id = hidmanager.getSubdetid(detid, tTopo, false);
 
       std::map<std::string, LayerMEs>::iterator iLayerME = LayerMEsMap.find(layer_id);
       if (iLayerME == LayerMEsMap.end()) {
         folder_organizer.setLayerFolder(detid, tTopo, det_layer_pair.second, false);
         bookLayerMEs(ibooker, detid, layer_id);
       }
 
       const char* subdet = det_layer_pair.first.c_str();
       if (std::strstr(subdet, tec) != nullptr || std::strstr(subdet, tid) != nullptr) {
         std::string ring_id = hidmanager.getSubdetid(detid, tTopo, true);
         std::map<std::string, RingMEs>::iterator iRingME = RingMEsMap.find(ring_id);
         if (iRingME == RingMEsMap.end()) {
           std::pair<std::string, int32_t> det_ring_pair = folder_organizer.GetSubDetAndLayer(detid, tTopo, true);
           folder_organizer.setLayerFolder(detid, tTopo, det_ring_pair.second, true);
           bookRingMEs(ibooker, detid, ring_id);
         }
       }
 
       // book sub-detector plots
       std::pair<std::string, std::string> sdet_pair = folder_organizer.getSubDetFolderAndTag(detid, tTopo);
       if (SubDetMEsMap.find(sdet_pair.second) == SubDetMEsMap.end()) {
         ibooker.setCurrentFolder(sdet_pair.first);
         bookSubDetMEs(ibooker, sdet_pair.second);
       }
     }  //end loop on detectors detid
   }
 
   //book control view plots
   if (!(topFolderName_.find("IsolatedBunches") != std::string::npos)) {
     ibooker.setCurrentFolder(topFolderName_ + "/ControlView/");
 
     ClusterStoNCorr_OnTrack_TIBTID =
         ibooker.book1D("ClusterStoNCorr_OnTrack_TIBTID", "TIB/TID [FECCrate=1] (OnTrack)", 100, 0., 100.);
     ClusterStoNCorr_OnTrack_TIBTID->setAxisTitle("S/N", 1);
 
     ClusterStoNCorr_OnTrack_TOB =
         ibooker.book1D("ClusterStoNCorr_OnTrack_TOB", "TOB [FECCrate=4] (OnTrack)", 100, 0., 100.);
     ClusterStoNCorr_OnTrack_TOB->setAxisTitle("S/N", 1);
 
     ClusterStoNCorr_OnTrack_TECM =
         ibooker.book1D("ClusterStoNCorr_OnTrack_TECM", "TECM [FECCrate=3] (OnTrack)", 100, 0., 100.);
     ClusterStoNCorr_OnTrack_TECM->setAxisTitle("S/N", 1);
 
     ClusterStoNCorr_OnTrack_TECP =
         ibooker.book1D("ClusterStoNCorr_OnTrack_TECP", "TECP [FECCrate=2] (OnTrack)", 100, 0., 100.);
     ClusterStoNCorr_OnTrack_TECP->setAxisTitle("S/N", 1);
 
     ClusterStoNCorr_OnTrack_FECCratevsFECSlot =
         ibooker.book2D("ClusterStoNCorr_OnTrack_FECCratevsFECSlot", " S/N (On track)", 22, 0.5, 22.5, 4, 0.5, 4.5);
     ClusterStoNCorr_OnTrack_FECCratevsFECSlot->setAxisTitle("FEC Slot", 1);
     ClusterStoNCorr_OnTrack_FECCratevsFECSlot->setAxisTitle("FEC Crate (TTC partition)", 2);
     ClusterStoNCorr_OnTrack_FECCratevsFECSlot->setBinLabel(1, "TIB/TID", 2);
     ClusterStoNCorr_OnTrack_FECCratevsFECSlot->setBinLabel(2, "TEC+", 2);
     ClusterStoNCorr_OnTrack_FECCratevsFECSlot->setBinLabel(3, "TEC-", 2);
     ClusterStoNCorr_OnTrack_FECCratevsFECSlot->setBinLabel(4, "TOB", 2);
 
     ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TIBTID = ibooker.book2D("ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TIBTID",
                                                                      "TIB/TID [FECCrate=1] (OnTrack)",
                                                                      10,
                                                                      -0.5,
                                                                      9.5,
                                                                      22,
                                                                      0.5,
                                                                      22.5);
     ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TIBTID->setAxisTitle("FEC Ring", 1);
     ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TIBTID->setAxisTitle("FEC Slot", 2);
 
     ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TOB = ibooker.book2D(
         "ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TOB", "TOB [FECCrate=4] (OnTrack)", 10, -0.5, 9.5, 22, 0.5, 22.5);
     ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TOB->setAxisTitle("FEC Ring", 1);
     ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TOB->setAxisTitle("FEC Slot", 2);
 
     ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TECM = ibooker.book2D(
         "ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TECM", "TEC- [FECCrate=3] (OnTrack)", 10, -0.5, 9.5, 22, 0.5, 22.5);
     ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TECM->setAxisTitle("FEC Ring", 1);
     ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TECM->setAxisTitle("FEC Slot", 2);
 
     ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TECP = ibooker.book2D(
         "ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TECP", "TEC- [FECCrate=2] (OnTrack)", 10, -0.5, 9.5, 22, 0.5, 22.5);
     ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TECP->setAxisTitle("FEC Ring", 1);
     ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TECP->setAxisTitle("FEC Slot", 2);
 
     //----------------------------------------
     // for conting the number of clusters, for the mean S/N calculation
     //book control view plots
 
     ClusterCount_OnTrack_FECCratevsFECSlot =
         ibooker.book2D("ClusterCount_OnTrack_FECCratevsFECSlot", " S/N (On track)", 22, 0.5, 22.5, 4, 0.5, 4.5);
     ClusterCount_OnTrack_FECCratevsFECSlot->setAxisTitle("FEC Slot", 1);
     ClusterCount_OnTrack_FECCratevsFECSlot->setAxisTitle("FEC Crate (TTC partition)", 2);
     ClusterCount_OnTrack_FECCratevsFECSlot->setBinLabel(1, "TIB/TID", 2);
     ClusterCount_OnTrack_FECCratevsFECSlot->setBinLabel(2, "TEC+", 2);
     ClusterCount_OnTrack_FECCratevsFECSlot->setBinLabel(3, "TEC-", 2);
     ClusterCount_OnTrack_FECCratevsFECSlot->setBinLabel(4, "TOB", 2);
 
     ClusterCount_OnTrack_FECSlotVsFECRing_TIBTID = ibooker.book2D(
         "ClusterCount_OnTrack_FECSlotVsFECRing_TIBTID", "TIB/TID [FECCrate=1] (OnTrack)", 10, -0.5, 9.5, 22, 0.5, 22.5);
     ClusterCount_OnTrack_FECSlotVsFECRing_TIBTID->setAxisTitle("FEC Ring", 1);
     ClusterCount_OnTrack_FECSlotVsFECRing_TIBTID->setAxisTitle("FEC Slot", 2);
 
     ClusterCount_OnTrack_FECSlotVsFECRing_TOB = ibooker.book2D(
         "ClusterCount_OnTrack_FECSlotVsFECRing_TOB", "TOB [FECCrate=4] (OnTrack)", 10, -0.5, 9.5, 22, 0.5, 22.5);
     ClusterCount_OnTrack_FECSlotVsFECRing_TOB->setAxisTitle("FEC Ring", 1);
     ClusterCount_OnTrack_FECSlotVsFECRing_TOB->setAxisTitle("FEC Slot", 2);
 
     ClusterCount_OnTrack_FECSlotVsFECRing_TECM = ibooker.book2D(
         "ClusterCount_OnTrack_FECSlotVsFECRing_TECM", "TEC- [FECCrate=3] (OnTrack)", 10, -0.5, 9.5, 22, 0.5, 22.5);
     ClusterCount_OnTrack_FECSlotVsFECRing_TECM->setAxisTitle("FEC Ring", 1);
     ClusterCount_OnTrack_FECSlotVsFECRing_TECM->setAxisTitle("FEC Slot", 2);
 
     ClusterCount_OnTrack_FECSlotVsFECRing_TECP = ibooker.book2D(
         "ClusterCount_OnTrack_FECSlotVsFECRing_TECP", "TEC- [FECCrate=2] (OnTrack)", 10, -0.5, 9.5, 22, 0.5, 22.5);
     ClusterCount_OnTrack_FECSlotVsFECRing_TECP->setAxisTitle("FEC Ring", 1);
     ClusterCount_OnTrack_FECSlotVsFECRing_TECP->setAxisTitle("FEC Slot", 2);
   }
 }

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 96 of file SiStripMonitorTrack.cc.

References book(), edm::EventSetup::get(), and edm::ESHandle< T >::product().

                                                                                                              {
   // Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   es.get<TrackerTopologyRcd>().get(tTopoHandle);
 
   edm::ESHandle<TkDetMap> tkDetMapHandle;
   es.get<TrackerTopologyRcd>().get(tkDetMapHandle);
 
   book(ibooker, tTopoHandle.product(), tkDetMapHandle.product());
 }

void SiStripMonitorTrack::bookLayerMEs	(	DQMStore::IBooker &	ibooker,
		const uint32_t	mod_id,
		std::string &	layer_id
	)

private

Definition at line 486 of file SiStripMonitorTrack.cc.

Referenced by book().

                                                                                                            {
   std::string name = "layer";
   std::string view = "layerView";
   std::string hname;
   std::string hpar;
   SiStripHistoId hidmanager;
 
   LayerMEs theLayerMEs;
 
   // Signal/Noise (w/ cluster harge corrected)
   hname = hidmanager.createHistoLayer("Summary_ClusterStoNCorr", name, layer_id, "OnTrack");
   hpar = "TH1ClusterStoNCorr";
   theLayerMEs.ClusterStoNCorrOnTrack = handleBookMEs(ibooker, view, layer_id, hpar, hname);
 
   // Cluster Gain
   hname = hidmanager.createHistoLayer("Summary_ClusterGain", name, layer_id, "");
   hpar = "TH1ClusterGain";
   theLayerMEs.ClusterGain = handleBookMEs(ibooker, view, layer_id, hpar, hname);
 
   // Cluster Charge Corrected
   hname = hidmanager.createHistoLayer("Summary_ClusterChargeCorr", name, layer_id, "OnTrack");
   hpar = "TH1ClusterChargeCorr";
   theLayerMEs.ClusterChargeCorrOnTrack = handleBookMEs(ibooker, view, layer_id, hpar, hname);
 
   // Cluster Charge (On and Off Track)
   hname = hidmanager.createHistoLayer("Summary_ClusterCharge", name, layer_id, "OnTrack");
   hpar = "TH1ClusterCharge";
   theLayerMEs.ClusterChargeOnTrack = handleBookMEs(ibooker, view, layer_id, hpar, hname);
 
   hname = hidmanager.createHistoLayer("Summary_ClusterCharge", name, layer_id, "OffTrack");
   hpar = "TH1ClusterCharge";
   theLayerMEs.ClusterChargeOffTrack = handleBookMEs(ibooker, view, layer_id, hpar, hname);
 
   // Cluster Charge Raw (On and Off Track)
   hname = hidmanager.createHistoLayer("Summary_ClusterChargeRaw", name, layer_id, "OnTrack");
   hpar = "TH1ClusterChargeRaw";
   theLayerMEs.ClusterChargeRawOnTrack = handleBookMEs(ibooker, view, layer_id, hpar, hname);
 
   hname = hidmanager.createHistoLayer("Summary_ClusterChargeRaw", name, layer_id, "OffTrack");
   hpar = "TH1ClusterChargeRaw";
   theLayerMEs.ClusterChargeRawOffTrack = handleBookMEs(ibooker, view, layer_id, hpar, hname);
 
   // Cluster Noise (On and Off Track)
   hname = hidmanager.createHistoLayer("Summary_ClusterNoise", name, layer_id, "OnTrack");
   hpar = "TH1ClusterNoise";
   theLayerMEs.ClusterNoiseOnTrack = handleBookMEs(ibooker, view, layer_id, hpar, hname);
 
   hname = hidmanager.createHistoLayer("Summary_ClusterNoise", name, layer_id, "OffTrack");
   hpar = "TH1ClusterNoise";
   theLayerMEs.ClusterNoiseOffTrack = handleBookMEs(ibooker, view, layer_id, hpar, hname);
 
   // Cluster Width (On and Off Track)
   hname = hidmanager.createHistoLayer("Summary_ClusterWidth", name, layer_id, "OnTrack");
   hpar = "TH1ClusterWidth";
   theLayerMEs.ClusterWidthOnTrack = handleBookMEs(ibooker, view, layer_id, hpar, hname);
 
   hname = hidmanager.createHistoLayer("Summary_ClusterWidth", name, layer_id, "OffTrack");
   hpar = "TH1ClusterWidth";
   theLayerMEs.ClusterWidthOffTrack = handleBookMEs(ibooker, view, layer_id, hpar, hname);
 
   //Cluster Position
   short total_nr_strips = SiStripDetCabling_->nApvPairs(mod_id) * 2 * 128;
   if (layer_id.find("TEC") != std::string::npos)
     total_nr_strips = 3 * 2 * 128;
 
   hname = hidmanager.createHistoLayer("Summary_ClusterPosition", name, layer_id, "OnTrack");
   hpar = "TH1ClusterPos";
   if (layer_id.find("TIB") != std::string::npos || layer_id.find("TOB") != std::string::npos ||
       (conf_.getParameter<edm::ParameterSet>(hpar.c_str())).getParameter<bool>(view.c_str()))
     theLayerMEs.ClusterPosOnTrack = ibooker.book1D(hname, hname, total_nr_strips, 0.5, total_nr_strips + 0.5);
 
   //----------------------
   //add 2D z-phi map per layer
   hname = hidmanager.createHistoLayer("Summary_ClusterPosition2D", name, layer_id, "OnTrack");
   hpar = "TH2ClusterPosTOB";
   if (layer_id.find("TOB") != std::string::npos)
     theLayerMEs.ClusterPosOnTrack2D = ibooker.book2D(hname, hname, 12, -110, 110, 300, -3.2, 3.2);
 
   hname = hidmanager.createHistoLayer("Summary_ClusterPosition2D", name, layer_id, "OnTrack");
   hpar = "TH2ClusterPosTIB";
   if (layer_id.find("TIB") != std::string::npos)
     theLayerMEs.ClusterPosOnTrack2D = ibooker.book2D(hname, hname, 12, -65, 65, 450, -3.2, 3.2);
 
   hname = hidmanager.createHistoLayer("Summary_ClusterPosition2D", name, layer_id, "OnTrack");
   hpar = "TH2ClusterPosTEC";
   if (layer_id.find("TEC") != std::string::npos) {
     static constexpr int nbinR = 8;
     static constexpr float rval[9] = {0, 21.2, 30.8, 40.4, 50.0, 60.0, 75.0, 90.0, 110.0};
     static constexpr int nmodulesPhi = 40 * 6;  //max number of APV  for a ring
     float phival[nmodulesPhi];
     for (int i = 0; i < nmodulesPhi; i++)
       phival[i] = -3.2 + 2 * i * 3.2 / nmodulesPhi;
 
     TH2F* temp = new TH2F("tmp", "tmp", nbinR, rval, nmodulesPhi - 1, phival);
     theLayerMEs.ClusterPosOnTrack2D = ibooker.book2D(hname, temp);
   }
 
   hname = hidmanager.createHistoLayer("Summary_ClusterPosition2D", name, layer_id, "OnTrack");
   hpar = "TH2ClusterPosTID";
   if (layer_id.find("TID") != std::string::npos) {
     static constexpr int nbinR = 4;
     static constexpr float rval[5] = {0, 21.2, 30.8, 40.4, 50.0};
     static constexpr int nmodulesPhi = 80 * 4;  //max number of APV  for a ring
     float phival[nmodulesPhi];
     for (int i = 0; i < nmodulesPhi; i++)
       phival[i] = -3.2 + i * 2 * 3.2 / nmodulesPhi;
 
     TH2F* temp = new TH2F("tmp", "tmp", nbinR, rval, nmodulesPhi - 1, phival);
     theLayerMEs.ClusterPosOnTrack2D = ibooker.book2D(hname, temp);
   }
 
   hname = hidmanager.createHistoLayer("Summary_ClusterPosition", name, layer_id, "OffTrack");
   hpar = "TH1ClusterPos";
   if (layer_id.find("TIB") != std::string::npos || layer_id.find("TOB") != std::string::npos ||
       (conf_.getParameter<edm::ParameterSet>(hpar.c_str())).getParameter<bool>(view.c_str()))
     theLayerMEs.ClusterPosOffTrack = ibooker.book1D(hname, hname, total_nr_strips, 0.5, total_nr_strips + 0.5);
 
   // dQ/dx
   hname = hidmanager.createHistoLayer("Summary_ClusterChargePerCMfromTrack", name, layer_id, "");
   hpar = "TH1ClusterChargePerCM";
   theLayerMEs.ClusterChargePerCMfromTrack = handleBookMEs(ibooker, view, layer_id, hpar, hname);
 
   hname = hidmanager.createHistoLayer("Summary_ClusterChargePerCMfromOrigin", name, layer_id, "OnTrack");
   hpar = "TH1ClusterChargePerCM";
   theLayerMEs.ClusterChargePerCMfromOriginOnTrack = handleBookMEs(ibooker, view, layer_id, hpar, hname);
 
   hname = hidmanager.createHistoLayer("Summary_ClusterChargePerCMfromOrigin", name, layer_id, "OffTrack");
   hpar = "TH1ClusterChargePerCM";
   theLayerMEs.ClusterChargePerCMfromOriginOffTrack = handleBookMEs(ibooker, view, layer_id, hpar, hname);
 
   //bookeeping
   LayerMEsMap[layer_id] = theLayerMEs;
 }

SiStripMonitorTrack::MonitorElement * SiStripMonitorTrack::bookME1D	(	DQMStore::IBooker &	ibooker,
		const char *	ParameterSetLabel,
		const char *	HistoName
	)

private

Definition at line 828 of file SiStripMonitorTrack.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), conf_, and edm::ParameterSet::getParameter().

Referenced by bookModMEs(), bookSubDetMEs(), and handleBookMEs().

                                                                                           {
   Parameters = conf_.getParameter<edm::ParameterSet>(ParameterSetLabel);
   return ibooker.book1D(HistoName,
                         HistoName,
                         Parameters.getParameter<int32_t>("Nbinx"),
                         Parameters.getParameter<double>("xmin"),
                         Parameters.getParameter<double>("xmax"));
 }

SiStripMonitorTrack::MonitorElement * SiStripMonitorTrack::bookME2D	(	DQMStore::IBooker &	ibooker,
		const char *	ParameterSetLabel,
		const char *	HistoName
	)

private

Definition at line 840 of file SiStripMonitorTrack.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), conf_, and edm::ParameterSet::getParameter().

                                                                                           {
   Parameters = conf_.getParameter<edm::ParameterSet>(ParameterSetLabel);
   return ibooker.book2D(HistoName,
                         HistoName,
                         Parameters.getParameter<int32_t>("Nbinx"),
                         Parameters.getParameter<double>("xmin"),
                         Parameters.getParameter<double>("xmax"),
                         Parameters.getParameter<int32_t>("Nbiny"),
                         Parameters.getParameter<double>("ymin"),
                         Parameters.getParameter<double>("ymax"));
 }

SiStripMonitorTrack::MonitorElement * SiStripMonitorTrack::bookME3D	(	DQMStore::IBooker &	ibooker,
		const char *	ParameterSetLabel,
		const char *	HistoName
	)

private

Definition at line 855 of file SiStripMonitorTrack.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book3D(), conf_, and edm::ParameterSet::getParameter().

                                                                                           {
   Parameters = conf_.getParameter<edm::ParameterSet>(ParameterSetLabel);
   return ibooker.book3D(HistoName,
                         HistoName,
                         Parameters.getParameter<int32_t>("Nbinx"),
                         Parameters.getParameter<double>("xmin"),
                         Parameters.getParameter<double>("xmax"),
                         Parameters.getParameter<int32_t>("Nbiny"),
                         Parameters.getParameter<double>("ymin"),
                         Parameters.getParameter<double>("ymax"),
                         Parameters.getParameter<int32_t>("Nbinz"),
                         Parameters.getParameter<double>("zmin"),
                         Parameters.getParameter<double>("zmax"));
 }

SiStripMonitorTrack::MonitorElement * SiStripMonitorTrack::bookMEProfile	(	DQMStore::IBooker &	ibooker,
		const char *	ParameterSetLabel,
		const char *	HistoName
	)

private

Definition at line 873 of file SiStripMonitorTrack.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::bookProfile(), conf_, and edm::ParameterSet::getParameter().

Referenced by bookModMEs().

                                                                                                {
   Parameters = conf_.getParameter<edm::ParameterSet>(ParameterSetLabel);
   return ibooker.bookProfile(HistoName,
                              HistoName,
                              Parameters.getParameter<int32_t>("Nbinx"),
                              Parameters.getParameter<double>("xmin"),
                              Parameters.getParameter<double>("xmax"),
                              Parameters.getParameter<int32_t>("Nbiny"),
                              Parameters.getParameter<double>("ymin"),
                              Parameters.getParameter<double>("ymax"),
                              "");
 }

SiStripMonitorTrack::MonitorElement * SiStripMonitorTrack::bookMETrend	(	DQMStore::IBooker &	ibooker,
		const char *	HistoName
	)

private

Definition at line 889 of file SiStripMonitorTrack.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::bookProfile(), conf_, edm::ParameterSet::getParameter(), dqm::impl::MonitorElement::getTH1(), hlt_dqm_clientPB-live_cfg::me, and MonitorElementData::TPROFILE.

Referenced by bookSubDetMEs().

                                                                                              {
   edm::ParameterSet ParametersTrend = conf_.getParameter<edm::ParameterSet>("Trending");
   MonitorElement* me = ibooker.bookProfile(HistoName,
                                            HistoName,
                                            ParametersTrend.getParameter<int32_t>("Nbins"),
                                            ParametersTrend.getParameter<double>("xmin"),
                                            ParametersTrend.getParameter<double>("xmax"),
                                            0,
                                            0,
                                            "");
   if (me->kind() == MonitorElement::Kind::TPROFILE)
     me->getTH1()->SetCanExtend(TH1::kAllAxes);
 
   if (!me)
     return me;
   me->setAxisTitle("Lumisection", 1);
   return me;
 }

void SiStripMonitorTrack::bookModMEs	(	DQMStore::IBooker &	ibooker,
		const uint32_t	id
	)

private

Definition at line 415 of file SiStripMonitorTrack.cc.

Referenced by book().

 {
   std::string name = "det";
   SiStripHistoId hidmanager;
   std::string hid = hidmanager.createHistoId("", name, id);
   std::map<std::string, ModMEs>::iterator iModME = ModMEsMap.find(hid);
   if (iModME == ModMEsMap.end()) {
     ModMEs theModMEs;
 
     // Cluster Width
     theModMEs.ClusterWidth =
         bookME1D(ibooker, "TH1ClusterWidth", hidmanager.createHistoId("ClusterWidth_OnTrack", name, id).c_str());
     // Cluster Gain
     theModMEs.ClusterGain =
         bookME1D(ibooker, "TH1ClusterGain", hidmanager.createHistoId("ClusterGain", name, id).c_str());
     // Cluster Charge
     theModMEs.ClusterCharge =
         bookME1D(ibooker, "TH1ClusterCharge", hidmanager.createHistoId("ClusterCharge_OnTrack", name, id).c_str());
     // Cluster Charge Raw (no gain )
     theModMEs.ClusterChargeRaw = bookME1D(
         ibooker, "TH1ClusterChargeRaw", hidmanager.createHistoId("ClusterChargeRaw_OnTrack", name, id).c_str());
     // Cluster Charge Corrected
     theModMEs.ClusterChargeCorr = bookME1D(
         ibooker, "TH1ClusterChargeCorr", hidmanager.createHistoId("ClusterChargeCorr_OnTrack", name, id).c_str());
     // Cluster StoN Corrected
     theModMEs.ClusterStoNCorr = bookME1D(
         ibooker, "TH1ClusterStoNCorrMod", hidmanager.createHistoId("ClusterStoNCorr_OnTrack", name, id).c_str());
     // Cluster Position
     short total_nr_strips = SiStripDetCabling_->nApvPairs(id) * 2 * 128;
     theModMEs.ClusterPos = ibooker.book1D(hidmanager.createHistoId("ClusterPosition_OnTrack", name, id).c_str(),
                                           hidmanager.createHistoId("ClusterPosition_OnTrack", name, id).c_str(),
                                           total_nr_strips,
                                           0.5,
                                           total_nr_strips + 0.5);
     // Cluster PGV
     theModMEs.ClusterPGV =
         bookMEProfile(ibooker, "TProfileClusterPGV", hidmanager.createHistoId("PGV_OnTrack", name, id).c_str());
     // Cluster Charge per cm
     theModMEs.ClusterChargePerCMfromTrack = bookME1D(
         ibooker, "TH1ClusterChargePerCM", hidmanager.createHistoId("ClusterChargePerCMfromTrack", name, id).c_str());
 
     theModMEs.ClusterChargePerCMfromOrigin = bookME1D(
         ibooker, "TH1ClusterChargePerCM", hidmanager.createHistoId("ClusterChargePerCMfromOrigin", name, id).c_str());
 
     ModMEsMap[hid] = theModMEs;
   }
 }

void SiStripMonitorTrack::bookRing	(	DQMStore::IBooker &	,
		const uint32_t	,
		std::string &
	)

private

void SiStripMonitorTrack::bookRingMEs	(	DQMStore::IBooker &	ibooker,
		const uint32_t	mod_id,
		std::string &	ring_id
	)

private

Definition at line 620 of file SiStripMonitorTrack.cc.

Referenced by book().

                                                                                                          {
   std::string name = "ring";
   std::string view = "ringView";
   std::string hname;
   std::string hpar;
   SiStripHistoId hidmanager;
 
   RingMEs theRingMEs;
 
   hname = hidmanager.createHistoLayer("Summary_ClusterStoNCorr", name, ring_id, "OnTrack");
   hpar = "TH1ClusterStoNCorr";
   theRingMEs.ClusterStoNCorrOnTrack = handleBookMEs(ibooker, view, ring_id, hpar, hname);
 
   // Cluster Gain
   hname = hidmanager.createHistoLayer("Summary_ClusterGain", name, ring_id, "");
   hpar = "TH1ClusterGain";
   theRingMEs.ClusterGain = handleBookMEs(ibooker, view, ring_id, hpar, hname);
 
   // Cluster Charge Corrected
   hname = hidmanager.createHistoLayer("Summary_ClusterChargeCorr", name, ring_id, "OnTrack");
   hpar = "TH1ClusterChargeCorr";
   theRingMEs.ClusterChargeCorrOnTrack = handleBookMEs(ibooker, view, ring_id, hpar, hname);
 
   // Cluster Charge (On and Off Track)
   hname = hidmanager.createHistoLayer("Summary_ClusterCharge", name, ring_id, "OnTrack");
   hpar = "TH1ClusterCharge";
   theRingMEs.ClusterChargeOnTrack = handleBookMEs(ibooker, view, ring_id, hpar, hname);
 
   hname = hidmanager.createHistoLayer("Summary_ClusterCharge", name, ring_id, "OffTrack");
   hpar = "TH1ClusterCharge";
   theRingMEs.ClusterChargeOffTrack = handleBookMEs(ibooker, view, ring_id, hpar, hname);
 
   // Cluster Charge Raw (no-gain), On and off track
   hname = hidmanager.createHistoLayer("Summary_ClusterChargeRaw", name, ring_id, "OnTrack");
   hpar = "TH1ClusterChargeRaw";
   theRingMEs.ClusterChargeRawOnTrack = handleBookMEs(ibooker, view, ring_id, hpar, hname);
 
   hname = hidmanager.createHistoLayer("Summary_ClusterChargeRaw", name, ring_id, "OffTrack");
   hpar = "TH1ClusterChargeRaw";
   theRingMEs.ClusterChargeRawOffTrack = handleBookMEs(ibooker, view, ring_id, hpar, hname);
 
   // Cluster Noise (On and Off Track)
   hname = hidmanager.createHistoLayer("Summary_ClusterNoise", name, ring_id, "OnTrack");
   hpar = "TH1ClusterNoise";
   theRingMEs.ClusterNoiseOnTrack = handleBookMEs(ibooker, view, ring_id, hpar, hname);
 
   hname = hidmanager.createHistoLayer("Summary_ClusterNoise", name, ring_id, "OffTrack");
   hpar = "TH1ClusterNoise";
   theRingMEs.ClusterNoiseOffTrack = handleBookMEs(ibooker, view, ring_id, hpar, hname);
 
   // Cluster Width (On and Off Track)
   hname = hidmanager.createHistoLayer("Summary_ClusterWidth", name, ring_id, "OnTrack");
   hpar = "TH1ClusterWidth";
   theRingMEs.ClusterWidthOnTrack = handleBookMEs(ibooker, view, ring_id, hpar, hname);
 
   hname = hidmanager.createHistoLayer("Summary_ClusterWidth", name, ring_id, "OffTrack");
   hpar = "TH1ClusterWidth";
   theRingMEs.ClusterWidthOffTrack = handleBookMEs(ibooker, view, ring_id, hpar, hname);
 
   //Cluster Position
   short total_nr_strips = SiStripDetCabling_->nApvPairs(mod_id) * 2 * 128;
   if (ring_id.find("TEC") != std::string::npos)
     total_nr_strips = 3 * 2 * 128;
 
   hname = hidmanager.createHistoLayer("Summary_ClusterPosition", name, ring_id, "OnTrack");
   hpar = "TH1ClusterPos";
   if ((conf_.getParameter<edm::ParameterSet>(hpar.c_str())).getParameter<bool>(view.c_str()))
     theRingMEs.ClusterPosOnTrack = ibooker.book1D(hname, hname, total_nr_strips, 0.5, total_nr_strips + 0.5);
 
   hname = hidmanager.createHistoLayer("Summary_ClusterPosition", name, ring_id, "OffTrack");
   hpar = "TH1ClusterPos";
   if ((conf_.getParameter<edm::ParameterSet>(hpar.c_str())).getParameter<bool>(view.c_str()))
     theRingMEs.ClusterPosOffTrack = ibooker.book1D(hname, hname, total_nr_strips, 0.5, total_nr_strips + 0.5);
 
   // dQ/dx
   hname = hidmanager.createHistoLayer("Summary_ClusterChargePerCMfromTrack", name, ring_id, "");
   hpar = "TH1ClusterChargePerCM";
   theRingMEs.ClusterChargePerCMfromTrack = handleBookMEs(ibooker, view, ring_id, hpar, hname);
 
   hname = hidmanager.createHistoLayer("Summary_ClusterChargePerCMfromOrigin", name, ring_id, "OnTrack");
   hpar = "TH1ClusterChargePerCM";
   theRingMEs.ClusterChargePerCMfromOriginOnTrack = handleBookMEs(ibooker, view, ring_id, hpar, hname);
 
   hname = hidmanager.createHistoLayer("Summary_ClusterChargePerCMfromOrigin", name, ring_id, "OffTrack");
   hpar = "TH1ClusterChargePerCM";
   theRingMEs.ClusterChargePerCMfromOriginOffTrack = handleBookMEs(ibooker, view, ring_id, hpar, hname);
 
   //bookeeping
   RingMEsMap[ring_id] = theRingMEs;
 }

void SiStripMonitorTrack::bookSubDetMEs	(	DQMStore::IBooker &	ibooker,
		std::string &	name
	)

private

Definition at line 714 of file SiStripMonitorTrack.cc.

Referenced by book().

                                                                                  {
   std::string subdet_tag;
   subdet_tag = "__" + name;
   std::string completeName;
   std::string axisName;
 
   SubDetMEs theSubDetMEs;
 
   // TotalNumber of Cluster OnTrack
   completeName = "Summary_TotalNumberOfClusters_OnTrack" + subdet_tag;
   axisName = "Number of on-track clusters in " + name;
   theSubDetMEs.nClustersOnTrack = bookME1D(ibooker, "TH1nClustersOn", completeName.c_str());
   theSubDetMEs.nClustersOnTrack->setAxisTitle(axisName);
   theSubDetMEs.nClustersOnTrack->setStatOverflows(kTRUE);
 
   // TotalNumber of Cluster OffTrack
   completeName = "Summary_TotalNumberOfClusters_OffTrack" + subdet_tag;
   axisName = "Number of off-track clusters in " + name;
   theSubDetMEs.nClustersOffTrack = bookME1D(ibooker, "TH1nClustersOff", completeName.c_str());
   theSubDetMEs.nClustersOffTrack->setAxisTitle(axisName);
 
   double xmaximum = 0;
   if (name.find("TIB") != std::string::npos) {
     xmaximum = 40000.0;
     theSubDetMEs.nClustersOffTrack->setAxisRange(0.0, xmaximum, 1);
   }
   if (name.find("TOB") != std::string::npos) {
     xmaximum = 40000.0;
     theSubDetMEs.nClustersOffTrack->setAxisRange(0.0, xmaximum, 1);
   }
   if (name.find("TID") != std::string::npos) {
     xmaximum = 10000.0;
     theSubDetMEs.nClustersOffTrack->setAxisRange(0.0, xmaximum, 1);
   }
   if (name.find("TEC") != std::string::npos) {
     xmaximum = 40000.0;
     theSubDetMEs.nClustersOffTrack->setAxisRange(0.0, xmaximum, 1);
   }
 
   theSubDetMEs.nClustersOffTrack->setStatOverflows(kTRUE);
 
   // Cluster Gain
   completeName = "Summary_ClusterGain" + subdet_tag;
   theSubDetMEs.ClusterGain = bookME1D(ibooker, "TH1ClusterGain", completeName.c_str());
 
   // Cluster StoN On Track
   completeName = "Summary_ClusterStoNCorr_OnTrack" + subdet_tag;
   theSubDetMEs.ClusterStoNCorrOnTrack = bookME1D(ibooker, "TH1ClusterStoNCorr", completeName.c_str());
 
   completeName = "Summary_ClusterStoNCorrThin_OnTrack" + subdet_tag;
   if (subdet_tag.find("TEC") != std::string::npos)
     theSubDetMEs.ClusterStoNCorrThinOnTrack = bookME1D(ibooker, "TH1ClusterStoNCorr", completeName.c_str());
 
   completeName = "Summary_ClusterStoNCorrThick_OnTrack" + subdet_tag;
   if (subdet_tag.find("TEC") != std::string::npos)
     theSubDetMEs.ClusterStoNCorrThickOnTrack = bookME1D(ibooker, "TH1ClusterStoNCorr", completeName.c_str());
 
   // Cluster Charge Corrected
   completeName = "Summary_ClusterChargeCorr_OnTrack" + subdet_tag;
   theSubDetMEs.ClusterChargeCorrOnTrack = bookME1D(ibooker, "TH1ClusterChargeCorr", completeName.c_str());
 
   completeName = "Summary_ClusterChargeCorrThin_OnTrack" + subdet_tag;
   if (subdet_tag.find("TEC") != std::string::npos)
     theSubDetMEs.ClusterChargeCorrThinOnTrack = bookME1D(ibooker, "TH1ClusterChargeCorr", completeName.c_str());
 
   completeName = "Summary_ClusterChargeCorrThick_OnTrack" + subdet_tag;
   if (subdet_tag.find("TEC") != std::string::npos)
     theSubDetMEs.ClusterChargeCorrThickOnTrack = bookME1D(ibooker, "TH1ClusterChargeCorr", completeName.c_str());
 
   // Cluster Charge On Track
   completeName = "Summary_ClusterCharge_OnTrack" + subdet_tag;
   theSubDetMEs.ClusterChargeOnTrack = bookME1D(ibooker, "TH1ClusterCharge", completeName.c_str());
 
   // Cluster Charge On Track, Raw (no-gain)
   completeName = "Summary_ClusterChargeRaw_OnTrack" + subdet_tag;
   theSubDetMEs.ClusterChargeRawOnTrack = bookME1D(ibooker, "TH1ClusterChargeRaw", completeName.c_str());
 
   // Cluster Charge Off Track
   completeName = "Summary_ClusterCharge_OffTrack" + subdet_tag;
   theSubDetMEs.ClusterChargeOffTrack = bookME1D(ibooker, "TH1ClusterCharge", completeName.c_str());
 
   // Cluster Charge Off Track, Raw (no-gain)
   completeName = "Summary_ClusterChargeRaw_OffTrack" + subdet_tag;
   theSubDetMEs.ClusterChargeRawOffTrack = bookME1D(ibooker, "TH1ClusterChargeRaw", completeName.c_str());
 
   // Cluster Charge StoN Off Track
   completeName = "Summary_ClusterStoN_OffTrack" + subdet_tag;
   theSubDetMEs.ClusterStoNOffTrack = bookME1D(ibooker, "TH1ClusterStoN", completeName.c_str());
 
   // cluster charge per cm on track
   completeName = "Summary_ClusterChargePerCMfromTrack" + subdet_tag;
   theSubDetMEs.ClusterChargePerCMfromTrack = bookME1D(ibooker, "TH1ClusterChargePerCM", completeName.c_str());
 
   // cluster charge per cm on track
   completeName = "Summary_ClusterChargePerCMfromOrigin_OnTrack" + subdet_tag;
   theSubDetMEs.ClusterChargePerCMfromOriginOnTrack = bookME1D(ibooker, "TH1ClusterChargePerCM", completeName.c_str());
 
   // cluster charge per cm off track
   completeName = "Summary_ClusterChargePerCMfromOrigin_OffTrack" + subdet_tag;
   theSubDetMEs.ClusterChargePerCMfromOriginOffTrack = bookME1D(ibooker, "TH1ClusterChargePerCM", completeName.c_str());
 
   if (Trend_On_) {
     // TotalNumber of Cluster
     completeName = "Trend_TotalNumberOfClusters_OnTrack" + subdet_tag;
     theSubDetMEs.nClustersTrendOnTrack = bookMETrend(ibooker, completeName.c_str());
     completeName = "Trend_TotalNumberOfClusters_OffTrack" + subdet_tag;
     theSubDetMEs.nClustersTrendOffTrack = bookMETrend(ibooker, completeName.c_str());
   }
 
   //bookeeping
   SubDetMEsMap[name] = theSubDetMEs;
 }

bool SiStripMonitorTrack::clusterInfos	(	SiStripClusterInfo *	cluster,
		const uint32_t	detid,
		enum ClusterFlags	flags,
		bool	track_ok,
		LocalVector	LV,
		const Det2MEs &	MEs,
		const TrackerTopology *	tTopo,
		const SiStripGain *	stripGain,
		const SiStripQuality *	stripQuality,
		const edm::DetSetVector< SiStripDigi > &	digilist,
		float	clustZ,
		float	clustPhi
	)

private

Definition at line 1520 of file SiStripMonitorTrack.cc.

Referenced by AllClusters(), and RecHitInfo().

                                                    {
   if (cluster == nullptr)
     return false;
   // if one imposes a cut on the clusters, apply it
   if ((applyClusterQuality_) &&
       (cluster->signalOverNoise() < sToNLowerLimit_ || cluster->signalOverNoise() > sToNUpperLimit_ ||
        cluster->width() < widthLowerLimit_ || cluster->width() > widthUpperLimit_))
     return false;
   // start of the analysis
 
   float cosRZ = -2;
   LogDebug("SiStripMonitorTrack") << "\n\tLV " << LV.x() << " " << LV.y() << " " << LV.z() << " " << LV.mag()
                                   << std::endl;
   if (LV.mag()) {
     cosRZ = fabs(LV.z()) / LV.mag();
     LogDebug("SiStripMonitorTrack") << "\n\t cosRZ " << cosRZ << std::endl;
   }
 
   // Filling SubDet/Layer Plots (on Track + off Track)
   float StoN = cluster->signalOverNoise();
   float noise = cluster->noiseRescaledByGain();
   uint16_t charge = cluster->charge();
   uint16_t width = cluster->width();
   float position = cluster->baryStrip();
 
   // Getting raw charge with strip gain.
   double chargeraw = 0;
   double clustergain = 0;
   auto digi_it = digilist.find(detid);  //(digilist.isValid() ? digilist.find(detid) : digilist.end());
   // SiStripClusterInfo.stripCharges() <==> SiStripCluster.amplitudes()
   for (size_t chidx = 0; chidx < cluster->stripCharges().size(); ++chidx) {
     if (cluster->stripCharges().at(chidx) <= 0) {
       continue;
     }  // nonzero amplitude
     if (stripQuality->IsStripBad(stripQuality->getRange(detid), cluster->firstStrip() + chidx)) {
       continue;
     }
     clustergain += stripGain->getStripGain(cluster->firstStrip() + chidx, stripGain->getRange(detid));
     // Getting raw adc charge from digi collections
     if (digi_it == digilist.end()) {
       continue;
     }  // skipping if not found
     for (const auto& digiobj : *digi_it) {
       if (digiobj.strip() == cluster->firstStrip() + chidx) {
         chargeraw += digiobj.adc();
       }
     }
   }
 
   clustergain /= double(cluster->stripCharges().size());  // calculating average gain inside cluster
 
   // new dE/dx (chargePerCM)
   // https://indico.cern.ch/event/342236/session/5/contribution/10/material/slides/0.pdf
   float dQdx_fromTrack = siStripClusterTools::chargePerCM(detid, *cluster, LV);
   // from straigth line origin-sensor centre
   const StripGeomDetUnit* DetUnit = static_cast<const StripGeomDetUnit*>(tkgeom_->idToDetUnit(DetId(detid)));
   LocalPoint locVtx = DetUnit->toLocal(GlobalPoint(0.0, 0.0, 0.0));
   LocalVector locDir(locVtx.x(), locVtx.y(), locVtx.z());
   float dQdx_fromOrigin = siStripClusterTools::chargePerCM(detid, *cluster, locDir);
 
   if (TkHistoMap_On_ && (flag == OnTrack)) {
     uint32_t adet = cluster->detId();
     if (track_ok)
       tkhisto_ClChPerCMfromTrack->fill(adet, dQdx_fromTrack);
   }
   if (clchCMoriginTkHmap_On_ && (flag == OffTrack)) {
     uint32_t adet = cluster->detId();
     if (track_ok)
       tkhisto_ClChPerCMfromOrigin->fill(adet, dQdx_fromOrigin);
   }
 
   if (flag == OnTrack) {
     if (MEs.iSubdet != nullptr)
       MEs.iSubdet->totNClustersOnTrack++;
     // layerMEs
     if (MEs.iLayer != nullptr) {
       if (noise > 0.0)
         fillME(MEs.iLayer->ClusterStoNCorrOnTrack, StoN * cosRZ);
       if (noise == 0.0)
         LogDebug("SiStripMonitorTrack") << "Module " << detid << " in Event " << eventNb << " noise "
                                         << cluster->noiseRescaledByGain() << std::endl;
       fillME(MEs.iLayer->ClusterGain, clustergain);
       fillME(MEs.iLayer->ClusterChargeCorrOnTrack, charge * cosRZ);
       fillME(MEs.iLayer->ClusterChargeOnTrack, charge);
       fillME(MEs.iLayer->ClusterChargeRawOnTrack, chargeraw);
       fillME(MEs.iLayer->ClusterNoiseOnTrack, noise);
       fillME(MEs.iLayer->ClusterWidthOnTrack, width);
       fillME(MEs.iLayer->ClusterPosOnTrack, position);
       //auto clustgp = cluster->globalPosition();
 
       fillME(MEs.iLayer->ClusterPosOnTrack2D, valX, valY);
 
       if (track_ok)
         fillME(MEs.iLayer->ClusterChargePerCMfromTrack, dQdx_fromTrack);
       if (track_ok)
         fillME(MEs.iLayer->ClusterChargePerCMfromOriginOnTrack, dQdx_fromOrigin);
     }
     // ringMEs
     if (MEs.iRing != nullptr) {
       if (noise > 0.0)
         fillME(MEs.iRing->ClusterStoNCorrOnTrack, StoN * cosRZ);
       if (noise == 0.0)
         LogDebug("SiStripMonitorTrack") << "Module " << detid << " in Event " << eventNb << " noise "
                                         << cluster->noiseRescaledByGain() << std::endl;
       fillME(MEs.iRing->ClusterGain, clustergain);
       fillME(MEs.iRing->ClusterChargeCorrOnTrack, charge * cosRZ);
       fillME(MEs.iRing->ClusterChargeOnTrack, charge);
       fillME(MEs.iRing->ClusterChargeRawOnTrack, chargeraw);
       fillME(MEs.iRing->ClusterNoiseOnTrack, noise);
       fillME(MEs.iRing->ClusterWidthOnTrack, width);
       fillME(MEs.iRing->ClusterPosOnTrack, position);
       if (track_ok)
         fillME(MEs.iRing->ClusterChargePerCMfromTrack, dQdx_fromTrack);
       if (track_ok)
         fillME(MEs.iRing->ClusterChargePerCMfromOriginOnTrack, dQdx_fromOrigin);
     }
     // subdetMEs
     if (MEs.iSubdet != nullptr) {
       fillME(MEs.iSubdet->ClusterGain, clustergain);
       fillME(MEs.iSubdet->ClusterChargeOnTrack, charge);
       fillME(MEs.iSubdet->ClusterChargeRawOnTrack, chargeraw);
       if (noise > 0.0)
         fillME(MEs.iSubdet->ClusterStoNCorrOnTrack, StoN * cosRZ);
       fillME(MEs.iSubdet->ClusterChargeCorrOnTrack, charge * cosRZ);
       if (track_ok)
         fillME(MEs.iSubdet->ClusterChargePerCMfromTrack, dQdx_fromTrack);
       if (track_ok)
         fillME(MEs.iSubdet->ClusterChargePerCMfromOriginOnTrack, dQdx_fromOrigin);
       if (tTopo->moduleGeometry(detid) == SiStripModuleGeometry::W5 ||
           tTopo->moduleGeometry(detid) == SiStripModuleGeometry::W6 ||
           tTopo->moduleGeometry(detid) == SiStripModuleGeometry::W7) {
         if (noise > 0.0)
           fillME(MEs.iSubdet->ClusterStoNCorrThickOnTrack, StoN * cosRZ);
         fillME(MEs.iSubdet->ClusterChargeCorrThickOnTrack, charge * cosRZ);
       } else {
         if (noise > 0.0)
           fillME(MEs.iSubdet->ClusterStoNCorrThinOnTrack, StoN * cosRZ);
         fillME(MEs.iSubdet->ClusterChargeCorrThinOnTrack, charge * cosRZ);
       }
     }
     //******** TkHistoMaps
     if (TkHistoMap_On_) {
       uint32_t adet = cluster->detId();
       tkhisto_NumOnTrack->add(adet, 1.);
       if (noise > 0.0)
         tkhisto_StoNCorrOnTrack->fill(adet, cluster->signalOverNoise() * cosRZ);
       if (noise == 0.0)
         LogDebug("SiStripMonitorTrack") << "Module " << detid << " in Event " << eventNb << " noise " << noise
                                         << std::endl;
     }
     if (TkHistoMap_On_ && (flag == OnTrack)) {
       uint32_t adet = cluster->detId();
       tkhisto_ClusterWidthOnTrack->fill(adet, cluster->width());
       if (noise > 0.0)
         tkhisto_NoiseOnTrack->fill(adet, cluster->noiseRescaledByGain() * cosRZ);
     }
 
     // Module plots filled only for onTrack Clusters
     if (Mod_On_) {
       SiStripHistoId hidmanager2;
       std::string name = hidmanager2.createHistoId("", "det", detid);
       //fillModMEs
       std::map<std::string, ModMEs>::iterator iModME = ModMEsMap.find(name);
       if (iModME != ModMEsMap.end()) {
         if (noise > 0.0)
           fillME(iModME->second.ClusterStoNCorr, StoN * cosRZ);
         if (noise == 0.0)
           LogDebug("SiStripMonitorTrack")
               << "Module " << name << " in Event " << eventNb << " noise " << noise << std::endl;
         fillME(iModME->second.ClusterGain, clustergain);
         fillME(iModME->second.ClusterCharge, charge);
         fillME(iModME->second.ClusterChargeRaw, chargeraw);
 
         fillME(iModME->second.ClusterChargeCorr, charge * cosRZ);
 
         fillME(iModME->second.ClusterWidth, width);
         fillME(iModME->second.ClusterPos, position);
 
         if (track_ok)
           fillME(iModME->second.ClusterChargePerCMfromTrack, dQdx_fromTrack);
         if (track_ok)
           fillME(iModME->second.ClusterChargePerCMfromOrigin, dQdx_fromOrigin);
 
         //fill the PGV histo
         float PGVmax = cluster->maxCharge();
         int PGVposCounter = cluster->maxIndex();
         for (int i = int(conf_.getParameter<edm::ParameterSet>("TProfileClusterPGV").getParameter<double>("xmin"));
              i < PGVposCounter;
              ++i)
           fillME(iModME->second.ClusterPGV, i, 0.);
         for (auto it = cluster->stripCharges().begin(); it < cluster->stripCharges().end(); ++it) {
           fillME(iModME->second.ClusterPGV, PGVposCounter++, (*it) / PGVmax);
         }
         for (int i = PGVposCounter;
              i < int(conf_.getParameter<edm::ParameterSet>("TProfileClusterPGV").getParameter<double>("xmax"));
              ++i)
           fillME(iModME->second.ClusterPGV, i, 0.);
         //end fill the PGV histo
       }
     }
   } else {
     if (flag == OffTrack) {
       if (MEs.iSubdet != nullptr)
         MEs.iSubdet->totNClustersOffTrack++;
       //******** TkHistoMaps
       if (TkHistoMap_On_) {
         uint32_t adet = cluster->detId();
         tkhisto_NumOffTrack->add(adet, 1.);
         if (charge > 250) {
           LogDebug("SiStripMonitorTrack") << "Module firing " << detid << " in Event " << eventNb << std::endl;
         }
       }
       if (TkHistoMap_On_) {
         uint32_t adet = cluster->detId();
         tkhisto_ClusterWidthOffTrack->fill(adet, cluster->width());
         if (noise > 0.0)
           tkhisto_NoiseOffTrack->fill(adet, cluster->noiseRescaledByGain());
       }
     }
     // layerMEs
     if (MEs.iLayer != nullptr) {
       fillME(MEs.iLayer->ClusterGain, clustergain);
       fillME(MEs.iLayer->ClusterChargeOffTrack, charge);
       fillME(MEs.iLayer->ClusterChargeRawOffTrack, chargeraw);
       fillME(MEs.iLayer->ClusterNoiseOffTrack, noise);
       fillME(MEs.iLayer->ClusterWidthOffTrack, width);
       fillME(MEs.iLayer->ClusterPosOffTrack, position);
       fillME(MEs.iLayer->ClusterChargePerCMfromOriginOffTrack, dQdx_fromOrigin);
     }
     // ringMEs
     if (MEs.iRing != nullptr) {
       fillME(MEs.iRing->ClusterGain, clustergain);
       fillME(MEs.iRing->ClusterChargeOffTrack, charge);
       fillME(MEs.iRing->ClusterChargeRawOffTrack, chargeraw);
       fillME(MEs.iRing->ClusterNoiseOffTrack, noise);
       fillME(MEs.iRing->ClusterWidthOffTrack, width);
       fillME(MEs.iRing->ClusterPosOffTrack, position);
       fillME(MEs.iRing->ClusterChargePerCMfromOriginOffTrack, dQdx_fromOrigin);
     }
     // subdetMEs
     if (MEs.iSubdet != nullptr) {
       fillME(MEs.iSubdet->ClusterGain, clustergain);
       fillME(MEs.iSubdet->ClusterChargeOffTrack, charge);
       fillME(MEs.iSubdet->ClusterChargeRawOffTrack, chargeraw);
       if (noise > 0.0)
         fillME(MEs.iSubdet->ClusterStoNOffTrack, StoN);
       fillME(MEs.iSubdet->ClusterChargePerCMfromOriginOffTrack, dQdx_fromOrigin);
     }
   }
   return true;
 }

void SiStripMonitorTrack::dqmBeginRun	(	const edm::Run &	run,
		const edm::EventSetup &	es
	)

override

Definition at line 75 of file SiStripMonitorTrack.cc.

References TrackerGeometry::detUnits(), genTriggerEventFlag_, edm::EventSetup::get(), GenericTriggerEventFlag::initRun(), LogDebug, GenericTriggerEventFlag::on(), edm::ESHandle< T >::product(), SiStripDetCabling_, stripGain_, tkgeom_, and tTopo_.

                                                                                 {
   //get geom
   es.get<TrackerDigiGeometryRecord>().get(tkgeom_);
   LogDebug("SiStripMonitorTrack") << "[SiStripMonitorTrack::beginRun] There are " << tkgeom_->detUnits().size()
                                   << " detectors instantiated in the geometry" << std::endl;
   es.get<SiStripDetCablingRcd>().get(SiStripDetCabling_);
 
   edm::ESHandle<TrackerTopology> tTopoHandle;
   es.get<TrackerTopologyRcd>().get(tTopoHandle);
   tTopo_ = tTopoHandle.product();
 
   edm::ESHandle<SiStripGain> gainHandle;
   es.get<SiStripGainRcd>().get(gainHandle);
   stripGain_ = gainHandle.product();
 
   // Initialize the GenericTriggerEventFlag
   if (genTriggerEventFlag_->on())
     genTriggerEventFlag_->initRun(run, es);
 }

bool SiStripMonitorTrack::fillControlViewHistos	(	const edm::Event &	ev,
		const edm::EventSetup &	es
	)

private

Definition at line 1388 of file SiStripMonitorTrack.cc.

Referenced by analyze().

                                                                                            {
   edm::Handle<reco::TrackCollection> tracks;
   ev.getByToken(trackToken_, tracks);  //takes the track collection
 
   //check that tracks are valid
   if (!tracks.isValid())
     return false;
 
   // loop over the tracks
   for (const auto& track : *tracks) {
     // loop over the rechits of this track
     for (trackingRecHit_iterator hit = track.recHitsBegin(), ehit = track.recHitsEnd(); hit != ehit; ++hit) {
       uint32_t thedetid = (*hit)->rawId();
       if (!(DetId(thedetid).subdetId() >= 3 && DetId(thedetid).subdetId() <= 6)) {
         continue;
       }
 
       if (!(*hit)->isValid())
         continue;
 
       const TrackingRecHit* theHit = (*hit);
       if (theHit == nullptr) {
         continue;
       }
 
       edm::ESHandle<TrackerTopology> tTopoHandle;
       es.get<TrackerTopologyRcd>().get(tTopoHandle);
       const TrackerTopology* const tTopo = tTopoHandle.product();
 
       const SiStripRecHit1D* hit1D = dynamic_cast<const SiStripRecHit1D*>(theHit);
       const SiStripRecHit2D* hit2D = dynamic_cast<const SiStripRecHit2D*>(theHit);
 
       float sovn = -1.;
       if (hit1D && !hit2D) {
         const SiStripCluster* SiStripCluster_ = &*(hit1D->cluster());
         SiStripClusterInfo SiStripClusterInfo_(*SiStripCluster_, es, thedetid);
         sovn = SiStripClusterInfo_.signalOverNoise();
       }
 
       else if (!hit1D && hit2D) {
         const SiStripCluster* SiStripCluster_ = &*(hit2D->cluster());
         SiStripClusterInfo SiStripClusterInfo_(*SiStripCluster_, es, thedetid);
         sovn = SiStripClusterInfo_.signalOverNoise();
       }
 
       std::vector<const FedChannelConnection*> getFedChanConnections;
       getFedChanConnections = SiStripDetCabling_->getConnections(thedetid);
 
       //      SiStripFolderOrganizer folder_organizer;
       //      std::string sistripsubdet = folder_organizer.getSubDetFolderAndTag(thedetid, tTopo).second;
 
       // loop over the fed chan connections
       for (const auto& getFedChanConnection : getFedChanConnections) {
         if (getFedChanConnection == nullptr) {
           continue;
         }
 
         int binfeccrate = getFedChanConnection->fecCrate();
         int binfecslot = getFedChanConnection->fecSlot();
         int binfecring = getFedChanConnection->fecRing();
         //int binccuchan  = getFedChanConnections[i0]->ccuChan(); //will be used in a new PR
         //int binccuadd   = getFedChanConnections[i0]->ccuAddr(); //will be used in a new PR
 
         return2DME(ClusterStoNCorr_OnTrack_FECCratevsFECSlot,
                    ClusterCount_OnTrack_FECCratevsFECSlot,
                    binfecslot,
                    binfeccrate,
                    sovn);
 
         // TIB/TID
         //        if ((sistripsubdet.find("TIB")) || (sistripsubdet.find("TID"))) {
         if ((DetId(thedetid).subdetId() == SiStripDetId::TIB) || ((DetId(thedetid).subdetId() == SiStripDetId::TID))) {
           ClusterStoNCorr_OnTrack_TIBTID->Fill(sovn);
           return2DME(ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TIBTID,
                      ClusterCount_OnTrack_FECSlotVsFECRing_TIBTID,
                      binfecring,
                      binfecslot,
                      sovn);
         }
 
         // TOB
         if (DetId(thedetid).subdetId() == SiStripDetId::TOB) {
           ClusterStoNCorr_OnTrack_TOB->Fill(sovn);
           return2DME(ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TOB,
                      ClusterCount_OnTrack_FECSlotVsFECRing_TOB,
                      binfecring,
                      binfecslot,
                      sovn);
         }
 
         // TECM
         if ((DetId(thedetid).subdetId() == SiStripDetId::TEC) && (tTopo->tecSide(thedetid) == 1)) {
           ClusterStoNCorr_OnTrack_TECM->Fill(sovn);
           return2DME(ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TECM,
                      ClusterCount_OnTrack_FECSlotVsFECRing_TECM,
                      binfecring,
                      binfecslot,
                      sovn);
         }
 
         // TECP
         if ((DetId(thedetid).subdetId() == SiStripDetId::TEC) && (tTopo->tecSide(thedetid) == 2)) {
           ClusterStoNCorr_OnTrack_TECP->Fill(sovn);
           return2DME(ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TECP,
                      ClusterCount_OnTrack_FECSlotVsFECRing_TECP,
                      binfecring,
                      binfecslot,
                      sovn);
         }
 
       }  // end of looping over the fed chan connections
     }    // end of looping over the rechits of the track
   }      // end of looping over the tracks
 
   return true;
 }

void SiStripMonitorTrack::fillME	(	MonitorElement *	ME,
		float	value1
	)

inlineprivate

Definition at line 147 of file SiStripMonitorTrack.h.

References dqm::impl::MonitorElement::Fill().

Referenced by analyze(), and clusterInfos().

                                                        {
     if (ME != nullptr)
       ME->Fill(value1);
   }

void SiStripMonitorTrack::fillME	(	MonitorElement *	ME,
		float	value1,
		float	value2
	)

inlineprivate

Definition at line 151 of file SiStripMonitorTrack.h.

References dqm::impl::MonitorElement::Fill().

                                                                      {
     if (ME != nullptr)
       ME->Fill(value1, value2);
   }

void SiStripMonitorTrack::fillME	(	MonitorElement *	ME,
		float	value1,
		float	value2,
		float	value3
	)

inlineprivate

Definition at line 155 of file SiStripMonitorTrack.h.

References dqm::impl::MonitorElement::Fill().

                                                                                    {
     if (ME != nullptr)
       ME->Fill(value1, value2, value3);
   }

void SiStripMonitorTrack::fillME	(	MonitorElement *	ME,
		float	value1,
		float	value2,
		float	value3,
		float	value4
	)

inlineprivate

Definition at line 159 of file SiStripMonitorTrack.h.

References dqm::impl::MonitorElement::Fill(), and findMEs().

                                                                                                  {
     if (ME != nullptr)
       ME->Fill(value1, value2, value3, value4);
   }

SiStripMonitorTrack::Det2MEs SiStripMonitorTrack::findMEs	(	const TrackerTopology *	tTopo,
		const uint32_t	detid
	)

private

Definition at line 1358 of file SiStripMonitorTrack.cc.

References folderOrganizer_, SiStripFolderOrganizer::getSubDetFolderAndTag(), SiStripHistoId::getSubdetid(), SiStripMonitorTrack::Det2MEs::iLayer, SiStripMonitorTrack::Det2MEs::iRing, SiStripMonitorTrack::Det2MEs::iSubdet, LayerMEsMap, hlt_dqm_clientPB-live_cfg::me, alignCSCRings::ring_id, RingMEsMap, AlCaHLTBitMon_QueryRunRegistry::string, and SubDetMEsMap.

Referenced by AllClusters(), fillME(), and RecHitInfo().

                                                                                                           {
   SiStripHistoId hidmanager1;
 
   std::string layer_id = hidmanager1.getSubdetid(detid, tTopo, false);
   std::string ring_id = hidmanager1.getSubdetid(detid, tTopo, true);
   std::string sdet_tag = folderOrganizer_.getSubDetFolderAndTag(detid, tTopo).second;
 
   Det2MEs me;
   me.iLayer = nullptr;
   me.iRing = nullptr;
   me.iSubdet = nullptr;
 
   std::map<std::string, LayerMEs>::iterator iLayer = LayerMEsMap.find(layer_id);
   if (iLayer != LayerMEsMap.end()) {
     me.iLayer = &(iLayer->second);
   }
 
   std::map<std::string, RingMEs>::iterator iRing = RingMEsMap.find(ring_id);
   if (iRing != RingMEsMap.end()) {
     me.iRing = &(iRing->second);
   }
 
   std::map<std::string, SubDetMEs>::iterator iSubdet = SubDetMEsMap.find(sdet_tag);
   if (iSubdet != SubDetMEsMap.end()) {
     me.iSubdet = &(iSubdet->second);
   }
 
   return me;
 }

SiStripMonitorTrack::MonitorElement * SiStripMonitorTrack::handleBookMEs	(	DQMStore::IBooker &	ibooker,
		std::string &	viewParameter,
		std::string &	id,
		std::string &	histoParameters,
		std::string &	histoName
	)

private

Definition at line 463 of file SiStripMonitorTrack.cc.

References bookME1D(), conf_, spr::find(), edm::ParameterSet::getParameter(), and hlt_dqm_clientPB-live_cfg::me.

Referenced by bookLayerMEs(), and bookRingMEs().

                                                                                               {
   MonitorElement* me = nullptr;
   bool view = false;
   view = (conf_.getParameter<edm::ParameterSet>(histoParameters.c_str())).getParameter<bool>(viewParameter.c_str());
   if (id.find("TEC") == std::string::npos && id.find("TID") == std::string::npos) {
     me = bookME1D(ibooker, histoParameters.c_str(), histoName.c_str());
   } else {
     if (view) {
       //      histoName = histoName + "__" + thickness;
       me = bookME1D(ibooker, histoParameters.c_str(), histoName.c_str());
     }
   }
   return me;
 }

void SiStripMonitorTrack::hitStudy	(	const edm::Event &	ev,
		const edm::EventSetup &	es,
		const edm::DetSetVector< SiStripDigi > &	digilist,
		const ProjectedSiStripRecHit2D *	projhit,
		const SiStripMatchedRecHit2D *	matchedhit,
		const SiStripRecHit2D *	hit2D,
		const SiStripRecHit1D *	hit1D,
		LocalVector	localMomentum,
		const bool	track_ok
	)

private

Definition at line 1003 of file SiStripMonitorTrack.cc.

References ev, TrackingRecHit::geographicalId(), TrackerGeometry::idToDet(), LogDebug, LogTrace, visualization-live-secondInstance_cfg::m, PV3DBase< T, PVType, FrameType >::mag(), GluedGeomDet::monoDet(), SiStripMatchedRecHit2D::monoHit(), ProjectedSiStripRecHit2D::originalHit(), DetId::rawId(), alignCSCRings::s, GluedGeomDet::stereoDet(), SiStripMatchedRecHit2D::stereoHit(), tkgeom_, GeomDet::toGlobal(), GeomDet::toLocal(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by trackStudyFromTrack().

                                                         {
   LocalVector statedirection;
   if (matchedhit) {  // type=Matched;
     LogTrace("SiStripMonitorTrack") << "\nMatched recHit found" << std::endl;
 
     const GluedGeomDet* gdet = static_cast<const GluedGeomDet*>(tkgeom_->idToDet(matchedhit->geographicalId()));
 
     GlobalVector gtrkdirup = gdet->toGlobal(localMomentum);
 
     //mono side
     const GeomDetUnit* monodet = gdet->monoDet();
     statedirection = monodet->toLocal(gtrkdirup);
     SiStripRecHit2D m = matchedhit->monoHit();
     if (statedirection.mag())
       RecHitInfo<SiStripRecHit2D>(&m, statedirection, digilist, ev, es, track_ok);
 
     //stereo side
     const GeomDetUnit* stereodet = gdet->stereoDet();
     statedirection = stereodet->toLocal(gtrkdirup);
     SiStripRecHit2D s = matchedhit->stereoHit();
     if (statedirection.mag())
       RecHitInfo<SiStripRecHit2D>(&s, statedirection, digilist, ev, es, track_ok);
   } else if (projhit) {  // type=Projected;
     LogTrace("SiStripMonitorTrack") << "\nProjected recHit found" << std::endl;
 
     const GluedGeomDet* gdet = static_cast<const GluedGeomDet*>(tkgeom_->idToDet(projhit->geographicalId()));
 
     GlobalVector gtrkdirup = gdet->toGlobal(localMomentum);
     const SiStripRecHit2D originalhit = projhit->originalHit();
 
     const GeomDetUnit* det;
     if (!StripSubdetector(originalhit.geographicalId().rawId()).stereo()) {
       //mono side
       LogTrace("SiStripMonitorTrack") << "\nProjected recHit found  MONO" << std::endl;
       det = gdet->monoDet();
       statedirection = det->toLocal(gtrkdirup);
       if (statedirection.mag())
         RecHitInfo<SiStripRecHit2D>(&(originalhit), statedirection, digilist, ev, es, track_ok);
     } else {
       LogTrace("SiStripMonitorTrack") << "\nProjected recHit found STEREO" << std::endl;
       //stereo side
       det = gdet->stereoDet();
       statedirection = det->toLocal(gtrkdirup);
       if (statedirection.mag())
         RecHitInfo<SiStripRecHit2D>(&(originalhit), statedirection, digilist, ev, es, track_ok);
     }
   } else if (hit2D) {  // type=2D
     statedirection = localMomentum;
     if (statedirection.mag())
       RecHitInfo<SiStripRecHit2D>(hit2D, statedirection, digilist, ev, es, track_ok);
   } else if (hit1D) {  // type=1D
     statedirection = localMomentum;
     if (statedirection.mag())
       RecHitInfo<SiStripRecHit1D>(hit1D, statedirection, digilist, ev, es, track_ok);
   } else {
     LogDebug("SiStripMonitorTrack") << " LocalMomentum: " << statedirection
                                     << "\nLocal x-z plane angle: " << atan2(statedirection.x(), statedirection.z());
   }
 }

template<class T >

void SiStripMonitorTrack::RecHitInfo	(	const T *	tkrecHit,
		LocalVector	LV,
		const edm::DetSetVector< SiStripDigi > &	digilist,
		const edm::Event &	ev,
		const edm::EventSetup &	es,
		bool	track_ok
	)

private

Definition at line 1221 of file SiStripMonitorTrack.cc.

References clusterInfos(), findMEs(), TrackerGeometry::idToDet(), LogTrace, EcalMEFormatter_cfi::MEs, OnTrack, PV3DBase< T, PVType, FrameType >::phi(), funct::pow(), stripGain_, stripQuality_, DetId::subdetId(), SiStripDetId::TIB, tkgeom_, SiStripDetId::TOB, GeomDet::toGlobal(), tTopo_, vPSiStripCluster, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                     {
   if (!tkrecHit->isValid()) {
     LogTrace("SiStripMonitorTrack") << "\t\t Invalid Hit " << std::endl;
     return;
   }
 
   const uint32_t& detid = tkrecHit->geographicalId().rawId();
 
   LogTrace("SiStripMonitorTrack")
       << "\n\t\tRecHit on det " << detid << "\n\t\tRecHit in LP " << tkrecHit->localPosition() << "\n\t\tRecHit in GP "
       << tkgeom_->idToDet(tkrecHit->geographicalId())->surface().toGlobal(tkrecHit->localPosition())
       << "\n\t\tRecHit trackLocal vector " << LV.x() << " " << LV.y() << " " << LV.z() << std::endl;
 
   //Get SiStripCluster from SiStripRecHit
   if (tkrecHit != nullptr && tkrecHit->isValid()) {
     const DetId detid = tkrecHit->geographicalId();
     int subDet = detid.subdetId();
     //std::cout << __LINE__ << std::endl;
     float clust_Pos1 = -1000;
     float clust_Pos2 = -1000;
 
     GlobalPoint theGlobalPos =
         tkgeom_->idToDet(tkrecHit->geographicalId())->surface().toGlobal(tkrecHit->localPosition());
     if (subDet == SiStripDetId::TIB || subDet == SiStripDetId::TOB) {
       clust_Pos1 = theGlobalPos.z();
       clust_Pos2 = theGlobalPos.phi();
     } else {
       clust_Pos1 = pow(theGlobalPos.x() * theGlobalPos.x() + theGlobalPos.y() * theGlobalPos.y(), 0.5);
       clust_Pos2 = theGlobalPos.phi();
     }
 
     const SiStripCluster* SiStripCluster_ = &*(tkrecHit->cluster());
     SiStripClusterInfo SiStripClusterInfo_(*SiStripCluster_, es, detid);
 
     const Det2MEs MEs = findMEs(tTopo_, detid);
     if (clusterInfos(&SiStripClusterInfo_,
                      detid,
                      OnTrack,
                      track_ok,
                      LV,
                      MEs,
                      tTopo_,
                      stripGain_,
                      stripQuality_,
                      digilist,
                      clust_Pos1,
                      clust_Pos2)) {
       vPSiStripCluster.insert(SiStripCluster_);
     }
   } else {
     edm::LogError("SiStripMonitorTrack") << "NULL hit" << std::endl;
   }
 }

void SiStripMonitorTrack::return2DME	(	MonitorElement *	input1,
		MonitorElement *	input2,
		int	binx,
		int	biny,
		double	value
	)

private

Definition at line 1505 of file SiStripMonitorTrack.cc.

References dqm::impl::MonitorElement::getBinContent(), dqm::impl::MonitorElement::setBinContent(), and relativeConstraints::value.

Referenced by fillControlViewHistos().

                                                                                                                      {
   if (input1->getBinContent(binx, biny) == 0.) {
     input1->setBinContent(binx, biny, value);
     input2->setBinContent(binx, biny, 1);
   } else {
     double nentries = input2->getBinContent(binx, biny);
     double theMeanSoN = (input1->getBinContent(binx, biny) * nentries + value) / (nentries + 1);
     //input1->setBinContent(binx,biny,((input1->getBinContent(binx,biny)+value)/2.));
     input1->setBinContent(binx, biny, theMeanSoN);
     input2->setBinContent(binx, biny, input2->getBinContent(binx, biny) + 1);
   }
 }

bool SiStripMonitorTrack::trackFilter ( const reco::Track & track )

private

Definition at line 1099 of file SiStripMonitorTrack.cc.

References reco::TrackBase::hitPattern(), reco::TrackBase::normalizedChi2(), reco::HitPattern::numberOfValidTrackerHits(), reco::TrackBase::p(), and reco::TrackBase::pt().

Referenced by trackStudyFromTrack(), and trackStudyFromTrajectory().

                                                             {
   if (track.pt() < 0.8)
     return false;
   if (track.p() < 2.0)
     return false;
   if (track.hitPattern().numberOfValidTrackerHits() <= 6)
     return false;
   if (track.normalizedChi2() > 10.0)
     return false;
   return true;
 }

void SiStripMonitorTrack::trackStudy	(	const edm::Event &	ev,
		const edm::EventSetup &	es
	)

private

Definition at line 1071 of file SiStripMonitorTrack.cc.

References digiToken_, edm::Event::getByToken(), edm::HandleBase::isValid(), edm::Handle< T >::product(), TrackLabel_, trackStudyFromTrajectory(), and trackToken_.

Referenced by analyze().

                                                                                 {
   using namespace std;
   using namespace edm;
   using namespace reco;
 
   //    edm::Handle<std::vector<Trajectory> > trajectories;
   //    ev.getByToken(trajectoryToken_, trajectories);
 
   // track input
   edm::Handle<reco::TrackCollection> trackCollectionHandle;
   ev.getByToken(trackToken_, trackCollectionHandle);  //takes the track collection
 
   // digis list
   edm::Handle<edm::DetSetVector<SiStripDigi>> digihandle;
   ev.getByToken(digiToken_, digihandle);
   edm::DetSetVector<SiStripDigi> dummy;
   auto& digilist = digihandle.isValid() ? *(digihandle.product()) : dummy;
 
   if (trackCollectionHandle.isValid()) {
     trackStudyFromTrajectory(trackCollectionHandle, digilist, ev, es);
   } else {
     edm::LogError("SiStripMonitorTrack") << "also Track Collection is not valid !! " << TrackLabel_ << std::endl;
     return;
   }
 }

void SiStripMonitorTrack::trackStudyFromTrack	(	edm::Handle< reco::TrackCollection >	trackCollectionHandle,
		const edm::DetSetVector< SiStripDigi > &	digilist,
		const edm::Event &	ev,
		const edm::EventSetup &	es
	)

private

Definition at line 1111 of file SiStripMonitorTrack.cc.

References DetId::det(), hitStudy(), tkhisto_NumberInactiveHits, tkhisto_NumberValidHits, tkhisto_NumMissingHits, TkHistoMap_On_, HLT_2018_cff::track, duplicaterechits_cfi::trackCollection, DetId::Tracker, and trackFilter().

                                                                        {
   //  edm::ESHandle<TransientTrackBuilder> builder;
   //  es.get<TransientTrackRecord>().get("TransientTrackBuilder",builder);
   //  const TransientTrackBuilder* transientTrackBuilder = builder.product();
 
   //numTracks = trackCollectionHandle->size();
   reco::TrackCollection trackCollection = *trackCollectionHandle;
   for (reco::TrackCollection::const_iterator track = trackCollection.begin(), etrack = trackCollection.end();
        track != etrack;
        ++track) {
     bool track_ok = trackFilter(*track);
     //    const reco::TransientTrack transientTrack = transientTrackBuilder->build(track);
 
     for (trackingRecHit_iterator hit = track->recHitsBegin(), ehit = track->recHitsEnd(); hit != ehit; ++hit) {
       if (TkHistoMap_On_) {
         uint32_t thedetid = (*hit)->rawId();
         if (SiStripDetId(thedetid).subDetector() >= 3 &&
             SiStripDetId(thedetid).subDetector() <= 6) {  //TIB/TID + TOB + TEC only
           if (((*hit)->getType() == 1))
             tkhisto_NumMissingHits->add(thedetid, 1.);
           if (((*hit)->getType() == 2))
             tkhisto_NumberInactiveHits->add(thedetid, 1.);
           if (((*hit)->getType() == 0))
             tkhisto_NumberValidHits->add(thedetid, 1.);
         }
       }
 
       if (!(*hit)->isValid())
         continue;
       DetId detID = (*hit)->geographicalId();
       if (detID.det() != DetId::Tracker)
         continue;
       const TrackingRecHit* theHit = (*hit);
       const ProjectedSiStripRecHit2D* projhit = dynamic_cast<const ProjectedSiStripRecHit2D*>((theHit));
       const SiStripMatchedRecHit2D* matchedhit = dynamic_cast<const SiStripMatchedRecHit2D*>((theHit));
       const SiStripRecHit2D* hit2D = dynamic_cast<const SiStripRecHit2D*>((theHit));
       const SiStripRecHit1D* hit1D = dynamic_cast<const SiStripRecHit1D*>((theHit));
 
       //      GlobalPoint globalPoint = hit->globalPosition();
       //      reco::TrajectoryStateOnSurface stateOnSurface = transientTrack->stateOnSurface(globalPoint);
 
       LocalVector localMomentum;
       hitStudy(ev, es, digilist, projhit, matchedhit, hit2D, hit1D, localMomentum, track_ok);
     }
 
     // hit pattern of the track
     // const reco::HitPattern & hitsPattern = track->hitPattern();
     // loop over the hits of the track
     //    for (int i=0; i<hitsPattern.numberOfAllHits(); i++) {
     // for (int i=0; i<hitsPattern.numberOfAllHits(reco::HitPattern::TRACK_HITS); i++) {
     //   uint32_t hit = hitsPattern.getHitPattern(reco::HitPattern::TRACK_HITS,i);
 
     // if the hit is valid and in pixel barrel, print out the layer
     //   if (hitsPattern.validHitFilter(hit) && hitsPattern.pixelBarrelHitFilter(hit))
 
     //   if (!hitsPattern.validHitFilter(hit)) continue;
     //      if (hitsPattern.pixelHitFilter(hit))       std::cout << "pixel"        << std::endl;       // pixel
     //      if (hitsPattern.pixelBarrelHitFilter(hit)) std::cout << "pixel barrel" << std::endl; // pixel barrel
     //      if (hitsPattern.pixelEndcapHitFilter(hit)) std::cout << "pixel endcap" << std::endl; // pixel endcap
     //      if (hitsPattern.stripHitFilter(hit))       std::cout << "strip" << std::endl;       // strip
     //      if (hitsPattern.stripTIBHitFilter(hit))    std::cout << "TIB" << std::endl;    // strip TIB
     //      if (hitsPattern.stripTIDHitFilter(hit))    std::cout << "TID" << std::endl;    // strip TID
     //      if (hitsPattern.stripTOBHitFilter(hit))    std::cout << "TOB" << std::endl;    // strip TOB
     //      if (hitsPattern.stripTECHitFilter(hit))    std::cout << "TEC" << std::endl;    // strip TEC
     //      if (hitsPattern.muonDTHitFilter(hit))      std::cout << "DT"  << std::endl;      // muon DT
     //      if (hitsPattern.muonCSCHitFilter(hit))     std::cout << "CSC" << std::endl;     // muon CSC
     //      if (hitsPattern.muonRPCHitFilter(hit))     std::cout << "RPC" << std::endl;     // muon RPC
     //
     //      // expert level: printout the hit in 10-bit binary format
     //      std::cout << "hit in 10-bit binary format = ";
     //      for (int j=9; j>=0; j--) {
     //        int bit = (hit >> j) & 0x1;
     //        std::cout << bit;
     //      }
     //      std::cout << std::endl;
     // }
   }
 }

void SiStripMonitorTrack::trackStudyFromTrajectory	(	edm::Handle< reco::TrackCollection >	trackCollectionHandle,
		const edm::DetSetVector< SiStripDigi > &	digilist,
		const edm::Event &	ev,
		const edm::EventSetup &	es
	)

private

Definition at line 1193 of file SiStripMonitorTrack.cc.

References mps_fire::i, LogDebug, numTracks, HLT_2018_cff::track, duplicaterechits_cfi::trackCollection, trackFilter(), and trajectoryStudy().

Referenced by trackStudy().

                                                                             {
   //Perform track study
   int i = 0;
   reco::TrackCollection trackCollection = *trackCollectionHandle;
   numTracks = trackCollection.size();
   for (reco::TrackCollection::const_iterator track = trackCollection.begin(), etrack = trackCollection.end();
        track != etrack;
        ++track) {
     LogDebug("SiStripMonitorTrack") << "Track number " << ++i << std::endl;
     //      << "\n\tmomentum: " << trackref->momentum()
     //      << "\n\tPT: " << trackref->pt()
     //      << "\n\tvertex: " << trackref->vertex()
     //      << "\n\timpact parameter: " << trackref->d0()
     //      << "\n\tcharge: " << trackref->charge()
     //      << "\n\tnormalizedChi2: " << trackref->normalizedChi2()
     //      <<"\n\tFrom EXTRA : "
     //      <<"\n\t\touter PT "<< trackref->outerPt()<<std::endl;
 
     //    trajectoryStudy(traj_iterator,trackref,es);
     bool track_ok = trackFilter(*track);
     trajectoryStudy(*track, digilist, ev, es, track_ok);
   }
 }

void SiStripMonitorTrack::trajectoryStudy	(	const reco::Track &	track,
		const edm::DetSetVector< SiStripDigi > &	digilist,
		const edm::Event &	ev,
		const edm::EventSetup &	es,
		bool	track_ok
	)

private

Definition at line 910 of file SiStripMonitorTrack.cc.

References ev, reco::Track::extra(), TrackingRecHit::geographicalId(), h, hcalSimParameters_cfi::hb, TrackerGeometry::idToDet(), LogDebug, LogTrace, visualization-live-secondInstance_cfg::m, GluedGeomDet::monoDet(), SiStripMatchedRecHit2D::monoHit(), ProjectedSiStripRecHit2D::originalHit(), DetId::rawId(), reco::Track::recHitsBegin(), reco::Track::recHitsSize(), alignCSCRings::s, GluedGeomDet::stereoDet(), SiStripMatchedRecHit2D::stereoHit(), tkgeom_, tkhisto_NumberInactiveHits, tkhisto_NumberValidHits, tkhisto_NumMissingHits, TkHistoMap_On_, GeomDet::toGlobal(), and GeomDet::toLocal().

Referenced by trackStudyFromTrajectory().

                                                          {
   auto const& trajParams = track.extra()->trajParams();
   assert(trajParams.size() == track.recHitsSize());
   auto hb = track.recHitsBegin();
   for (unsigned int h = 0; h < track.recHitsSize(); h++) {
     auto ttrh = *(hb + h);
 
     if (TkHistoMap_On_) {
       uint32_t thedetid = ttrh->rawId();
       if (SiStripDetId(thedetid).subDetector() >= 3 &&
           SiStripDetId(thedetid).subDetector() <= 6) {  //TIB/TID + TOB + TEC only
         if ((ttrh->getType() == 1))
           tkhisto_NumMissingHits->fill(thedetid, 1.);
         if ((ttrh->getType() == 2))
           tkhisto_NumberInactiveHits->fill(thedetid, 1.);
         if ((ttrh->getType() == 0))
           tkhisto_NumberValidHits->fill(thedetid, 1.);
       }
     }
 
     if (!ttrh->isValid())
       continue;
 
     //trajectory local direction and position on detector
     auto statedirection = trajParams[h].momentum();
 
     const ProjectedSiStripRecHit2D* projhit = dynamic_cast<const ProjectedSiStripRecHit2D*>(ttrh->hit());
     const SiStripMatchedRecHit2D* matchedhit = dynamic_cast<const SiStripMatchedRecHit2D*>(ttrh->hit());
     const SiStripRecHit2D* hit2D = dynamic_cast<const SiStripRecHit2D*>(ttrh->hit());
     const SiStripRecHit1D* hit1D = dynamic_cast<const SiStripRecHit1D*>(ttrh->hit());
 
     //      RecHitType type=Single;
 
     if (matchedhit) {
       LogTrace("SiStripMonitorTrack") << "\nMatched recHit found" << std::endl;
       //    type=Matched;
 
       const GluedGeomDet* gdet = static_cast<const GluedGeomDet*>(tkgeom_->idToDet(matchedhit->geographicalId()));
       GlobalVector gtrkdirup = gdet->toGlobal(statedirection);
       //mono side
       const GeomDetUnit* monodet = gdet->monoDet();
       statedirection = monodet->toLocal(gtrkdirup);
       SiStripRecHit2D m = matchedhit->monoHit();
       //      if(statedirection.mag() != 0)   RecHitInfo<SiStripRecHit2D>(&m,statedirection,trackref,es);
       if (statedirection.mag())
         RecHitInfo<SiStripRecHit2D>(&m, statedirection, digilist, ev, es, track_ok);
       //stereo side
       const GeomDetUnit* stereodet = gdet->stereoDet();
       statedirection = stereodet->toLocal(gtrkdirup);
       SiStripRecHit2D s = matchedhit->stereoHit();
       //      if(statedirection.mag() != 0)   RecHitInfo<SiStripRecHit2D>(&s,statedirection,trackref,es);
       if (statedirection.mag())
         RecHitInfo<SiStripRecHit2D>(&s, statedirection, digilist, ev, es, track_ok);
     } else if (projhit) {
       LogTrace("SiStripMonitorTrack") << "\nProjected recHit found" << std::endl;
       //    type=Projected;
       const GluedGeomDet* gdet = static_cast<const GluedGeomDet*>(tkgeom_->idToDet(projhit->geographicalId()));
 
       GlobalVector gtrkdirup = gdet->toGlobal(statedirection);
       const SiStripRecHit2D originalhit = projhit->originalHit();
       const GeomDetUnit* det;
       if (!StripSubdetector(originalhit.geographicalId().rawId()).stereo()) {
         //mono side
         LogTrace("SiStripMonitorTrack") << "\nProjected recHit found  MONO" << std::endl;
         det = gdet->monoDet();
         statedirection = det->toLocal(gtrkdirup);
         if (statedirection.mag())
           RecHitInfo<SiStripRecHit2D>(&(originalhit), statedirection, digilist, ev, es, track_ok);
       } else {
         LogTrace("SiStripMonitorTrack") << "\nProjected recHit found STEREO" << std::endl;
         //stereo side
         det = gdet->stereoDet();
         statedirection = det->toLocal(gtrkdirup);
         if (statedirection.mag())
           RecHitInfo<SiStripRecHit2D>(&(originalhit), statedirection, digilist, ev, es, track_ok);
       }
     } else if (hit2D) {
       if (statedirection.mag())
         RecHitInfo<SiStripRecHit2D>(hit2D, statedirection, digilist, ev, es, track_ok);
     } else if (hit1D) {
       if (statedirection.mag())
         RecHitInfo<SiStripRecHit1D>(hit1D, statedirection, digilist, ev, es, track_ok);
     } else {
       LogDebug("SiStripMonitorTrack") << " LocalMomentum: " << statedirection
                                       << "\nLocal x-z plane angle: " << atan2(statedirection.x(), statedirection.z());
     }
   }
 }