#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 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 &ev, edm::EventSetup const &es) final

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

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

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

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

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

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

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

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

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

	~DQMEDAnalyzer () override=default

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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

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

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)

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)

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_

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_

std::unordered_set< const SiStripCluster * >	vPSiStripCluster

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::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

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

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

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

void	consumesMany (const TypeToGet &id)

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

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

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

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

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

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 61 of file SiStripMonitorTrack.h.

Member Typedef Documentation

typedef TrackingRecHit::ConstRecHitPointer SiStripMonitorTrack::ConstRecHitPointer

Definition at line 63 of file SiStripMonitorTrack.h.

Member Enumeration Documentation

enum SiStripMonitorTrack::ClusterFlags

private

Enumerator
OffTrack
OnTrack

Definition at line 72 of file SiStripMonitorTrack.h.

                     {
     OffTrack,
     OnTrack
   };

enum SiStripMonitorTrack::RecHitType

Enumerator
Single
Matched
Projected
Null

Definition at line 64 of file SiStripMonitorTrack.h.

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

SiStripMonitorTrack::Matched

Definition: SiStripMonitorTrack.h:64

SiStripMonitorTrack::Null

Definition: SiStripMonitorTrack.h:64

SiStripMonitorTrack::Single

Definition: SiStripMonitorTrack.h:64

SiStripMonitorTrack::Projected

Definition: SiStripMonitorTrack.h:64

Constructor & Destructor Documentation

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

explicit

Definition at line 26 of file SiStripMonitorTrack.cc.

References applyClusterQuality_, clchCMoriginTkHmap_On_, Cluster_src_, clusterToken_, conf_, edm::EDConsumerBase::consumesCollector(), dcsStatus_, digiToken_, edm::ParameterSet::getParameter(), 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 63 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 1179 of file SiStripMonitorTrack.cc.

References clusterInfos(), clusterToken_, digiToken_, findMEs(), edm::EventSetup::get(), edm::Event::getByToken(), edm::HandleBase::isValid(), LogDebug, 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<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);
   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);
           clusterInfos(&SiStripClusterInfo_,detid,OffTrack, /*track_ok*/ false,LV,MEs, tTopo,stripGain,stripQuality,digilist);
         }
       }
     }
   } 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 95 of file SiStripMonitorTrack.cc.

References GenericTriggerEventFlag::accept(), AllClusters(), dcsStatus_, edm::EventID::event(), eventNb, fillControlViewHistos(), fillME(), genTriggerEventFlag_, 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::EventID::run(), SubDetMEsMap, topFolderName_, SiStripMonitorTrack::SubDetMEs::totNClustersOffTrack, SiStripMonitorTrack::SubDetMEs::totNClustersOnTrack, trackStudy(), Trend_On_, and vPSiStripCluster.

 {
   // 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;
   }
 
   //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 154 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 82 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 435 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);
 
   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;
 
 }

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

inlineprivate

Definition at line 732 of file SiStripMonitorTrack.cc.

References 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")
             );
 }

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

inlineprivate

Definition at line 743 of file SiStripMonitorTrack.cc.

References 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")
              );
 }

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

inlineprivate

Definition at line 757 of file SiStripMonitorTrack.cc.

References 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")
              );
 }

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

inlineprivate

Definition at line 774 of file SiStripMonitorTrack.cc.

References 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"),
                             "" );
 }

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

private

Definition at line 788 of file SiStripMonitorTrack.cc.

References DQMStore::IBooker::bookProfile(), conf_, MonitorElement::DQM_KIND_TPROFILE, edm::ParameterSet::getParameter(), and MonitorElement::getTH1().

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::DQM_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 370 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());
     ibooker.tag(theModMEs.ClusterWidth,id);
     // Cluster Gain
     theModMEs.ClusterGain=bookME1D(ibooker , "TH1ClusterGain", hidmanager.createHistoId("ClusterGain",name,id).c_str());
     ibooker.tag(theModMEs.ClusterGain,id);
     // Cluster Charge
     theModMEs.ClusterCharge=bookME1D(ibooker , "TH1ClusterCharge", hidmanager.createHistoId("ClusterCharge_OnTrack",name,id).c_str());
     ibooker.tag(theModMEs.ClusterCharge,id);
     // Cluster Charge Raw (no gain )
     theModMEs.ClusterChargeRaw=bookME1D(ibooker , "TH1ClusterChargeRaw", hidmanager.createHistoId("ClusterChargeRaw_OnTrack",name,id).c_str());
     ibooker.tag(theModMEs.ClusterChargeRaw,id);
     // Cluster Charge Corrected
     theModMEs.ClusterChargeCorr=bookME1D(ibooker , "TH1ClusterChargeCorr", hidmanager.createHistoId("ClusterChargeCorr_OnTrack",name,id).c_str());
     ibooker.tag(theModMEs.ClusterChargeCorr,id);
     // Cluster StoN Corrected
     theModMEs.ClusterStoNCorr=bookME1D(ibooker , "TH1ClusterStoNCorrMod", hidmanager.createHistoId("ClusterStoNCorr_OnTrack",name,id).c_str());
     ibooker.tag(theModMEs.ClusterStoNCorr,id);
     // 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);
     ibooker.tag(theModMEs.ClusterPos,id);
     // Cluster PGV
     theModMEs.ClusterPGV=bookMEProfile(ibooker , "TProfileClusterPGV", hidmanager.createHistoId("PGV_OnTrack",name,id).c_str());
     ibooker.tag(theModMEs.ClusterPGV,id);
     // Cluster Charge per cm
     theModMEs.ClusterChargePerCMfromTrack = bookME1D(ibooker , "TH1ClusterChargePerCM", hidmanager.createHistoId("ClusterChargePerCMfromTrack",name,id).c_str());
     ibooker.tag(theModMEs.ClusterChargePerCMfromTrack,id);
 
     theModMEs.ClusterChargePerCMfromOrigin = bookME1D(ibooker , "TH1ClusterChargePerCM", hidmanager.createHistoId("ClusterChargePerCMfromOrigin",name,id).c_str());
     ibooker.tag(theModMEs.ClusterChargePerCMfromOrigin,id);
 
     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 526 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 620 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->getTH1()->StatOverflows(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->getTH1()->StatOverflows(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
	)

private

Definition at line 1377 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);
       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) == SiStripDetId::ModuleGeometry::W5 || tTopo->moduleGeometry(detid) == SiStripDetId::ModuleGeometry::W6 || tTopo->moduleGeometry(detid) == SiStripDetId::ModuleGeometry::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 69 of file SiStripMonitorTrack.cc.

References TrackerGeometry::detUnits(), genTriggerEventFlag_, edm::EventSetup::get(), GenericTriggerEventFlag::initRun(), LogDebug, GenericTriggerEventFlag::on(), SiStripDetCabling_, and tkgeom_.

 {
   //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_ );
 
 
   // 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 1260 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 151 of file SiStripMonitorTrack.h.

References MonitorElement::Fill().

Referenced by analyze(), and clusterInfos().

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

MonitorElement::Fill

void Fill(long long x)

Definition: MonitorElement.h:166

hiDetachedQuadStep_cff.value1

value1

Definition: hiDetachedQuadStep_cff.py:70

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

inlineprivate

Definition at line 152 of file SiStripMonitorTrack.h.

References MonitorElement::Fill().

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

hiDetachedQuadStep_cff.value2

value2

Definition: hiDetachedQuadStep_cff.py:70

MonitorElement::Fill

void Fill(long long x)

Definition: MonitorElement.h:166

hiDetachedQuadStep_cff.value1

value1

Definition: hiDetachedQuadStep_cff.py:70

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

inlineprivate

Definition at line 153 of file SiStripMonitorTrack.h.

References MonitorElement::Fill().

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

hiDetachedQuadStep_cff.value2

value2

Definition: hiDetachedQuadStep_cff.py:70

MonitorElement::Fill

void Fill(long long x)

Definition: MonitorElement.h:166

hiDetachedQuadStep_cff.value1

value1

Definition: hiDetachedQuadStep_cff.py:70

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

inlineprivate

Definition at line 154 of file SiStripMonitorTrack.h.

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

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

hiDetachedQuadStep_cff.value2

value2

Definition: hiDetachedQuadStep_cff.py:70

MonitorElement::Fill

void Fill(long long x)

Definition: MonitorElement.h:166

hiDetachedQuadStep_cff.value1

value1

Definition: hiDetachedQuadStep_cff.py:70

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

private

Definition at line 1230 of file SiStripMonitorTrack.cc.

References folderOrganizer_, SiStripFolderOrganizer::getSubDetFolderAndTag(), SiStripHistoId::getSubdetid(), SiStripMonitorTrack::Det2MEs::iLayer, SiStripMonitorTrack::Det2MEs::iRing, SiStripMonitorTrack::Det2MEs::iSubdet, LayerMEsMap, 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;
 }

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

private

Definition at line 415 of file SiStripMonitorTrack.cc.

References bookME1D(), conf_, spr::find(), and edm::ParameterSet::getParameter().

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 899 of file SiStripMonitorTrack.cc.

References ev, TrackingRecHit::geographicalId(), TrackerGeometry::idToDet(), LogDebug, LogTrace, funct::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 1120 of file SiStripMonitorTrack.cc.

References clusterInfos(), findMEs(), edm::EventSetup::get(), TrackerGeometry::idToDet(), LogTrace, OnTrack, edm::ESHandle< T >::product(), tkgeom_, GeomDet::toGlobal(), 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;
 
     // FIXME: MOVE ALL THE EV AND ES ACCESS OUTSIDE THE LOOP!!!!
 
     //Retrieve tracker topology from geometry
     edm::ESHandle<TrackerTopology> tTopoHandle;
     es.get<TrackerTopologyRcd>().get(tTopoHandle);
     const TrackerTopology* const tTopo = tTopoHandle.product();
 
     // Getting SiStrip Gain settings
     edm::ESHandle<SiStripGain>  gainHandle;
     es.get<SiStripGainRcd>().get( gainHandle );
     const SiStripGain* const 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 );
 
     //Get SiStripCluster from SiStripRecHit
     if ( tkrecHit != nullptr ){
       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))
       {
         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 1359 of file SiStripMonitorTrack.cc.

References MonitorElement::getBinContent(), 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 997 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 967 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 1005 of file SiStripMonitorTrack.cc.

References DetId::det(), hitStudy(), tkhisto_NumberInactiveHits, tkhisto_NumberValidHits, tkhisto_NumMissingHits, TkHistoMap_On_, HiIsolationCommonParameters_cff::track, findElectronsInSiStrips_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 1089 of file SiStripMonitorTrack.cc.

References mps_fire::i, LogDebug, numTracks, HiIsolationCommonParameters_cff::track, findElectronsInSiStrips_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 804 of file SiStripMonitorTrack.cc.

References ev, reco::Track::extra(), TrackingRecHit::geographicalId(), h, hcalSimParameters_cfi::hb, TrackerGeometry::idToDet(), LogDebug, LogTrace, funct::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());
     }
 
   }
 
 }