#include <SiStripMonitorTrack.h>

Inheritance diagram for SiStripMonitorTrack:

Classes
struct	Det2MEs

struct	LayerMEs

struct	ModMEs

struct	RingMEs

struct	SubDetMEs

Public Types
typedef TrackingRecHit::ConstRecHitPointer	ConstRecHitPointer

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

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

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T...>

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T...>

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

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

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

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

	SiStripMonitorTrack (const edm::ParameterSet &)

	~SiStripMonitorTrack () override

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

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

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

void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer ()

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

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

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Private Types
enum	ClusterFlags { OffTrack, OnTrack }

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

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

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

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

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

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

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

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

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

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

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

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

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

bool	fillControlViewHistos (const edm::Event &ev)

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::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, 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)

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

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

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

Private Attributes
bool	applyClusterQuality_

bool	clchCMoriginTkHmap_On_

edm::InputTag	Cluster_src_

MonitorElement *	ClusterCount_OnTrack_FECCratevsFECSlot = nullptr

MonitorElement *	ClusterCount_OnTrack_FECSlotVsFECRing_TECM = nullptr

MonitorElement *	ClusterCount_OnTrack_FECSlotVsFECRing_TECP = nullptr

MonitorElement *	ClusterCount_OnTrack_FECSlotVsFECRing_TIBTID = nullptr

MonitorElement *	ClusterCount_OnTrack_FECSlotVsFECRing_TOB = nullptr

MonitorElement *	ClusterStoNCorr_OnTrack_FECCratevsFECSlot = nullptr

MonitorElement *	ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TECM = nullptr

MonitorElement *	ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TECP = nullptr

MonitorElement *	ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TIBTID = nullptr

MonitorElement *	ClusterStoNCorr_OnTrack_FECSlotVsFECRing_TOB = nullptr

MonitorElement *	ClusterStoNCorr_OnTrack_TECM = nullptr

MonitorElement *	ClusterStoNCorr_OnTrack_TECP = nullptr

MonitorElement *	ClusterStoNCorr_OnTrack_TIBTID = nullptr

MonitorElement *	ClusterStoNCorr_OnTrack_TOB = nullptr

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

SiStripClusterInfo	siStripClusterInfo_

const SiStripDetCabling *	siStripDetCabling_ = nullptr

edm::ESGetToken < SiStripDetCabling, SiStripDetCablingRcd >	siStripDetCablingToken_

double	sToNLowerLimit_

double	sToNUpperLimit_

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

edm::ESGetToken< TkDetMap, TrackerTopologyRcd >	tkDetMapToken_

const TrackerGeometry *	tkgeom_ = nullptr

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

edm::ESGetToken < TrackerGeometry, TrackerDigiGeometryRecord >	trackerGeometryToken_

const TrackerTopology *	trackerTopology_ = nullptr

edm::ESGetToken < TrackerTopology, TrackerTopologyRcd >	trackerTopologyEventToken_

edm::ESGetToken < TrackerTopology, TrackerTopologyRcd >	trackerTopologyRunToken_

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 DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr < DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

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

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 65 of file SiStripMonitorTrack.h.

Member Typedef Documentation

typedef TrackingRecHit::ConstRecHitPointer SiStripMonitorTrack::ConstRecHitPointer

Definition at line 67 of file SiStripMonitorTrack.h.

Member Enumeration Documentation

enum SiStripMonitorTrack::ClusterFlags

private

Enumerator
OffTrack
OnTrack

Definition at line 76 of file SiStripMonitorTrack.h.

76 { OffTrack, OnTrack };

SiStripMonitorTrack::OffTrack

Definition: SiStripMonitorTrack.h:76

SiStripMonitorTrack::OnTrack

Definition: SiStripMonitorTrack.h:76

enum SiStripMonitorTrack::RecHitType

Enumerator
Single
Matched
Projected
Null

Definition at line 68 of file SiStripMonitorTrack.h.

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

SiStripMonitorTrack::Matched

Definition: SiStripMonitorTrack.h:68

SiStripMonitorTrack::Null

Definition: SiStripMonitorTrack.h:68

SiStripMonitorTrack::Single

Definition: SiStripMonitorTrack.h:68

SiStripMonitorTrack::Projected

Definition: SiStripMonitorTrack.h:68

Constructor & Destructor Documentation

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

explicit

Definition at line 24 of file SiStripMonitorTrack.cc.

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

     : conf_(conf),
       siStripClusterInfo_(consumesCollector()),
       trackerGeometryToken_(esConsumes<TrackerGeometry, TrackerDigiGeometryRecord, edm::Transition::BeginRun>()),
       siStripDetCablingToken_(esConsumes<SiStripDetCabling, SiStripDetCablingRcd, edm::Transition::BeginRun>()),
       trackerTopologyRunToken_(esConsumes<TrackerTopology, TrackerTopologyRcd, edm::Transition::BeginRun>()),
       tkDetMapToken_(esConsumes<TkDetMap, TrackerTopologyRcd, edm::Transition::BeginRun>()),
       trackerTopologyEventToken_(esConsumes<TrackerTopology, TrackerTopologyRcd>()),
       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 68 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 )

private

Definition at line 1250 of file SiStripMonitorTrack.cc.

References clusterInfos(), clusterToken_, digiToken_, findMEs(), edm::Event::getByToken(), edm::HandleBase::isValid(), LogDebug, OffTrack, edm::Handle< T >::product(), SiStripClusterInfo::setCluster(), siStripClusterInfo_, SiStripClusterInfo::siStripGain(), SiStripClusterInfo::siStripQuality(), trackerTopology_, and vPSiStripCluster.

Referenced by analyze().

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

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 98 of file SiStripMonitorTrack.cc.

References GenericTriggerEventFlag::accept(), AllClusters(), dcsStatus_, edm::EventID::event(), eventNb, fillControlViewHistos(), fillME(), genTriggerEventFlag_, edm::EventSetup::getData(), SiStripDCSStatus::getStatus(), edm::EventBase::id(), iLumisection, SiStripClusterInfo::initEvent(), LogDebug, SiStripMonitorTrack::SubDetMEs::nClustersOffTrack, SiStripMonitorTrack::SubDetMEs::nClustersOnTrack, SiStripMonitorTrack::SubDetMEs::nClustersTrendOffTrack, SiStripMonitorTrack::SubDetMEs::nClustersTrendOnTrack, GenericTriggerEventFlag::on(), edm::EventBase::orbitNumber(), edm::EventID::run(), siStripClusterInfo_, SubDetMEsMap, topFolderName_, SiStripMonitorTrack::SubDetMEs::totNClustersOffTrack, SiStripMonitorTrack::SubDetMEs::totNClustersOnTrack, trackerTopology_, trackerTopologyEventToken_, trackStudy(), Trend_On_, and vPSiStripCluster.

                                                                                 {
   siStripClusterInfo_.initEvent(iSetup);
 
   // Filter out events if DCS checking is requested
   if (dcsStatus_ && !dcsStatus_->getStatus(e, iSetup))
     return;
 
   // Filter out events if Trigger Filtering is requested
   if (genTriggerEventFlag_->on() && !genTriggerEventFlag_->accept(e, iSetup))
     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;
   }
 
   trackerTopology_ = &iSetup.getData(trackerTopologyEventToken_);
 
   //Perform track study
   trackStudy(e);
 
   //Perform Cluster Study (irrespectively to tracks)
 
   AllClusters(e);  //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);
   }
 
   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 165 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 &	ibooker,
		edm::Run const &	run,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 88 of file SiStripMonitorTrack.cc.

References book(), edm::EventSetup::getData(), tkDetMapToken_, and trackerTopologyRunToken_.

                                                                       {
   const TrackerTopology* trackerTopology = &iSetup.getData(trackerTopologyRunToken_);
   const TkDetMap* tkDetMap = &iSetup.getData(tkDetMapToken_);
 
   book(ibooker, trackerTopology, tkDetMap);
 }

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

private

Definition at line 464 of file SiStripMonitorTrack.cc.

Referenced by book().

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

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

private

Definition at line 806 of file SiStripMonitorTrack.cc.

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

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

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

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

private

Definition at line 818 of file SiStripMonitorTrack.cc.

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

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

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

private

Definition at line 833 of file SiStripMonitorTrack.cc.

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

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

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

private

Definition at line 851 of file SiStripMonitorTrack.cc.

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

Referenced by bookModMEs().

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

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

private

Definition at line 867 of file SiStripMonitorTrack.cc.

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

Referenced by bookSubDetMEs().

                                                                                              {
   edm::ParameterSet ParametersTrend = conf_.getParameter<edm::ParameterSet>("Trending");
   MonitorElement* me = ibooker.bookProfile(HistoName,
                                            HistoName,
                                            ParametersTrend.getParameter<int32_t>("Nbins"),
                                            ParametersTrend.getParameter<double>("xmin"),
                                            ParametersTrend.getParameter<double>("xmax"),
                                            0,
                                            0,
                                            "");
   if (me->kind() == MonitorElement::Kind::TPROFILE)
     me->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 393 of file SiStripMonitorTrack.cc.

Referenced by book().

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

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

private

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

private

Definition at line 598 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 692 of file SiStripMonitorTrack.cc.

Referenced by book().

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

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

private

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

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 76 of file SiStripMonitorTrack.cc.

References TrackerGeometry::detUnits(), genTriggerEventFlag_, edm::EventSetup::getData(), GenericTriggerEventFlag::initRun(), LogDebug, GenericTriggerEventFlag::on(), siStripDetCabling_, siStripDetCablingToken_, tkgeom_, and trackerGeometryToken_.

                                                                                     {
   tkgeom_ = &iSetup.getData(trackerGeometryToken_);
   LogDebug("SiStripMonitorTrack") << "[SiStripMonitorTrack::beginRun] There are " << tkgeom_->detUnits().size()
                                   << " detectors instantiated in the geometry" << std::endl;
   siStripDetCabling_ = &iSetup.getData(siStripDetCablingToken_);
 
   // Initialize the GenericTriggerEventFlag
   if (genTriggerEventFlag_->on())
     genTriggerEventFlag_->initRun(run, iSetup);
 }

bool SiStripMonitorTrack::fillControlViewHistos ( const edm::Event & ev )

private

Definition at line 1342 of file SiStripMonitorTrack.cc.

Referenced by analyze().

                                                                   {
   const TrackerTopology* tTopo = trackerTopology_;
 
   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;
       }
 
       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_.setCluster(*SiStripCluster_, thedetid);
         sovn = siStripClusterInfo_.signalOverNoise();
       }
 
       else if (!hit1D && hit2D) {
         const SiStripCluster* SiStripCluster_ = &*(hit2D->cluster());
         siStripClusterInfo_.setCluster(*SiStripCluster_, 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 142 of file SiStripMonitorTrack.h.

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

Referenced by analyze(), and clusterInfos().

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

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

inlineprivate

Definition at line 146 of file SiStripMonitorTrack.h.

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

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

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

inlineprivate

Definition at line 150 of file SiStripMonitorTrack.h.

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

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

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

inlineprivate

Definition at line 154 of file SiStripMonitorTrack.h.

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

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

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

private

Definition at line 1312 of file SiStripMonitorTrack.cc.

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

Referenced by AllClusters(), and RecHitInfo().

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

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

private

Definition at line 441 of file SiStripMonitorTrack.cc.

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

Referenced by bookLayerMEs(), and bookRingMEs().

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

void SiStripMonitorTrack::hitStudy	(	const edm::Event &	ev,
		const edm::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 980 of file SiStripMonitorTrack.cc.

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

Referenced by trackStudyFromTrack().

                                                         {
   LocalVector statedirection;
   if (matchedhit) {  // type=Matched;
     LogTrace("SiStripMonitorTrack") << "\nMatched recHit found" << std::endl;
 
     const GluedGeomDet* gdet = static_cast<const GluedGeomDet*>(tkgeom_->idToDet(matchedhit->geographicalId()));
 
     GlobalVector gtrkdirup = gdet->toGlobal(localMomentum);
 
     //mono side
     const GeomDetUnit* monodet = gdet->monoDet();
     statedirection = monodet->toLocal(gtrkdirup);
     SiStripRecHit2D m = matchedhit->monoHit();
     if (statedirection.mag())
       RecHitInfo<SiStripRecHit2D>(&m, statedirection, digilist, ev, 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, 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, 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, track_ok);
     }
   } else if (hit2D) {  // type=2D
     statedirection = localMomentum;
     if (statedirection.mag())
       RecHitInfo<SiStripRecHit2D>(hit2D, statedirection, digilist, ev, track_ok);
   } else if (hit1D) {  // type=1D
     statedirection = localMomentum;
     if (statedirection.mag())
       RecHitInfo<SiStripRecHit1D>(hit1D, statedirection, digilist, ev, 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,
		bool	track_ok
	)

private

Definition at line 1191 of file SiStripMonitorTrack.cc.

References clusterInfos(), findMEs(), TrackerGeometry::idToDet(), LogTrace, OnTrack, PV3DBase< T, PVType, FrameType >::phi(), funct::pow(), SiStripClusterInfo::setCluster(), siStripClusterInfo_, SiStripClusterInfo::siStripGain(), SiStripClusterInfo::siStripQuality(), DetId::subdetId(), SiStripDetId::TIB, tkgeom_, SiStripDetId::TOB, GeomDet::toGlobal(), trackerTopology_, vPSiStripCluster, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

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

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

private

Definition at line 1457 of file SiStripMonitorTrack.cc.

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

Referenced by fillControlViewHistos().

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

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

private

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

private

Definition at line 1047 of file SiStripMonitorTrack.cc.

References digiToken_, edm::Event::getByToken(), dt_dqm_sourceclient_common_cff::reco, 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);
   } 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
	)

private

Definition at line 1087 of file SiStripMonitorTrack.cc.

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

                                                                   {
   //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, 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
	)

private

Definition at line 1164 of file SiStripMonitorTrack.cc.

References mps_fire::i, LogDebug, numTracks, HLT_FULL_cff::track, HLT_FULL_cff::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);
     bool track_ok = trackFilter(*track);
     trajectoryStudy(*track, digilist, ev, track_ok);
   }
 }

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

private

Definition at line 888 of file SiStripMonitorTrack.cc.

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

Referenced by trackStudyFromTrajectory().

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