#include <SegmentTrackAnalyzer.h>

Inheritance diagram for SegmentTrackAnalyzer:

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

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

	SegmentTrackAnalyzer (const edm::ParameterSet &)
	Constructor. More...

virtual	~SegmentTrackAnalyzer ()
	Destructor. More...

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

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

	DQMEDAnalyzer (void)

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

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

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
MonitorElement *	cscTrackHitPercentualVsEta

MonitorElement *	cscTrackHitPercentualVsPhi

MonitorElement *	cscTrackHitPercentualVsPt

MonitorElement *	dtTrackHitPercentualVsEta

MonitorElement *	dtTrackHitPercentualVsPhi

MonitorElement *	dtTrackHitPercentualVsPt

int	etaBin

double	etaMax

double	etaMin

MonitorElement *	hitsNotUsed

MonitorElement *	hitsNotUsedPercentual

MonitorElement *	hitStaProvenance

MonitorElement *	hitTkrProvenance

std::string	metname

edm::ParameterSet	parameters

int	phiBin

double	phiMax

double	phiMin

int	ptBin

double	ptMax

double	ptMin

edm::EDGetTokenT < reco::TrackCollection >	theMuTrackCollectionLabel_

SegmentsTrackAssociator *	theSegmentsAssociator

MuonServiceProxy *	theService

std::string	trackCollection

MonitorElement *	trackHitPercentualVsEta

MonitorElement *	trackHitPercentualVsPhi

MonitorElement *	trackHitPercentualVsPt

MonitorElement *	TrackSegm

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

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

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

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

void	consumesMany (const TypeToGet &id)

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

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

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

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

Detailed Description

DQM monitoring source for segments associated to the muon track

Author: G. Mila - INFN Torino

Definition at line 32 of file SegmentTrackAnalyzer.h.

Constructor & Destructor Documentation

SegmentTrackAnalyzer::SegmentTrackAnalyzer ( const edm::ParameterSet & pSet )

Constructor.

Definition at line 27 of file SegmentTrackAnalyzer.cc.

References instance, diffTwoXMLs::label, MuonServiceProxy_cff::MuonServiceProxy, Parameters::parameters, and PtMinSelector_cfg::ptMin.

                                                                       {
   parameters = pSet;
   
   // MuonService
   theService = new MuonServiceProxy(parameters.getParameter<ParameterSet>("ServiceParameters"));
   
   // Read Data:
   theMuTrackCollectionLabel_ = consumes<reco::TrackCollection>(parameters.getParameter<edm::InputTag>("MuTrackCollection"));
 
   // SegmentTrackAssociatior Parameters
   const ParameterSet SegmentsTrackAssociatorParameters = parameters.getParameter<ParameterSet>("SegmentsTrackAssociatorParameters");
 
   edm::ConsumesCollector iC   = consumesCollector();
   theSegmentsAssociator = new SegmentsTrackAssociator(SegmentsTrackAssociatorParameters,iC);
 
   trackCollection = parameters.getParameter<edm::InputTag>("MuTrackCollection").label() + parameters.getParameter<edm::InputTag>("MuTrackCollection").instance();
 
   etaBin = parameters.getParameter<int>("etaBin");
   etaMin = parameters.getParameter<double>("etaMin");
   etaMax = parameters.getParameter<double>("etaMax");
   phiBin = parameters.getParameter<int>("phiBin");
   phiMin = parameters.getParameter<double>("phiMin");
   phiMax = parameters.getParameter<double>("phiMax");
   ptBin  = parameters.getParameter<int>("ptBin");
   ptMin  = parameters.getParameter<double>("ptMin");
   ptMax  = parameters.getParameter<double>("ptMax");
 }

virtual SegmentTrackAnalyzer::~SegmentTrackAnalyzer ( )

inlinevirtual

Destructor.

Definition at line 39 of file SegmentTrackAnalyzer.h.

39 {};

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 110 of file SegmentTrackAnalyzer.cc.

References gather_cfg::cout, MuonSubdetId::CSC, MuonSubdetId::DT, edm::Event::getByToken(), DetId::Muon, DTRecSegment4D::phiSegment(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, MuonSubdetId::RPC, findQualityFiles::size, DTRecSegment2D::specificRecHits(), SiStripDetId::TEC, SiStripDetId::TIB, SiStripDetId::TID, SiStripDetId::TOB, and DetId::Tracker.

                                                                                      {
   theService->update(iSetup);
   
   Handle<reco::TrackCollection> glbTracks;
   iEvent.getByToken(theMuTrackCollectionLabel_,glbTracks);
   
   for (reco::TrackCollection::const_iterator recoTrack = glbTracks->begin(); recoTrack!=glbTracks->end(); ++recoTrack){
     MuonTransientTrackingRecHit::MuonRecHitContainer segments = theSegmentsAssociator->associate(iEvent, iSetup, *recoTrack );
   
 #ifdef DEBUG
     cout << "[SegmentTrackAnalyzer] # of segments associated to the track: "<<(segments).size() <<endl;
 #endif
 
     // hit counters
     int hitsFromDt=0;
     int hitsFromCsc=0;
     int hitsFromRpc=0;
     int hitsFromTk=0;
     int hitsFromTrack=0;
     int hitsFromSegmDt=0;
     int hitsFromSegmCsc=0;
     // segment counters
     int segmFromDt=0;
     int segmFromCsc=0;
 
     for (MuonTransientTrackingRecHit::MuonRecHitContainer::const_iterator segment=segments.begin();
      segment!=segments.end(); segment++) {
    
       DetId id = (*segment)->geographicalId();
     
       // hits from DT segments
       if (id.det() == DetId::Muon && id.subdetId() == MuonSubdetId::DT ) {
     ++segmFromDt;
     const DTRecSegment4D *seg4D = dynamic_cast<const DTRecSegment4D*>((*segment)->hit());
     if((*seg4D).hasPhi())
       hitsFromSegmDt+=(*seg4D).phiSegment()->specificRecHits().size();
     if((*seg4D).hasZed())
       hitsFromSegmDt+=(*seg4D).zSegment()->specificRecHits().size();
       
       }
     
       // hits from CSC segments
       if (id.det() == DetId::Muon && id.subdetId() == MuonSubdetId::CSC ) {
     hitsFromSegmCsc+=(*segment)->recHits().size();
     segmFromCsc++;
       }
 
     }
 
 
     // hits from track
     for(trackingRecHit_iterator recHit =  recoTrack->recHitsBegin(); recHit != recoTrack->recHitsEnd(); ++recHit){
 
       hitsFromTrack++;
       DetId id = (*recHit)->geographicalId();
       // hits from DT
       if (id.det() == DetId::Muon && id.subdetId() == MuonSubdetId::DT ) 
     hitsFromDt++;   
       // hits from CSC
       if (id.det() == DetId::Muon && id.subdetId() == MuonSubdetId::CSC ) 
     hitsFromCsc++;
       // hits from RPC
       if (id.det() == DetId::Muon && id.subdetId() == MuonSubdetId::RPC ) 
     hitsFromRpc++;
       // hits from Tracker
       if (id.det() == DetId::Tracker){
     hitsFromTk++;
     if(id.subdetId() == PixelSubdetector::PixelBarrel )
       hitTkrProvenance->Fill(1);
     if(id.subdetId() == PixelSubdetector::PixelEndcap )
       hitTkrProvenance->Fill(2);
     if(id.subdetId() == SiStripDetId::TIB )
       hitTkrProvenance->Fill(3);
     if(id.subdetId() == SiStripDetId::TID )
       hitTkrProvenance->Fill(4);
     if(id.subdetId() == SiStripDetId::TOB )
       hitTkrProvenance->Fill(5);
     if(id.subdetId() == SiStripDetId::TEC )
       hitTkrProvenance->Fill(6);
       }
 
     }
 
     // fill the histos
     hitsNotUsed->Fill(hitsFromSegmDt+hitsFromSegmCsc+hitsFromRpc+hitsFromTk-hitsFromTrack);
     hitsNotUsedPercentual->Fill(double(hitsFromSegmDt+hitsFromSegmCsc+hitsFromRpc+hitsFromTk-hitsFromTrack)/hitsFromTrack);
 
     if(hitsFromDt!=0 && hitsFromCsc!=0)
       TrackSegm->Fill(1,segmFromDt+segmFromCsc);
     if(hitsFromDt!=0 && hitsFromCsc==0)
       TrackSegm->Fill(2,segmFromDt);
     if(hitsFromDt==0 && hitsFromCsc!=0)
       TrackSegm->Fill(3,segmFromCsc);
 
     if(hitsFromDt!=0 && hitsFromCsc==0 && hitsFromRpc==0) hitStaProvenance->Fill(1);
     if(hitsFromCsc!=0 && hitsFromDt==0 && hitsFromRpc==0) hitStaProvenance->Fill(2);
     if(hitsFromRpc!=0 && hitsFromDt==0 && hitsFromCsc==0) hitStaProvenance->Fill(3);
     if(hitsFromDt!=0 && hitsFromCsc!=0 && hitsFromRpc==0) hitStaProvenance->Fill(4);
     if(hitsFromDt!=0 && hitsFromRpc!=0 && hitsFromCsc==0) hitStaProvenance->Fill(5);
     if(hitsFromCsc!=0 && hitsFromRpc!=0 && hitsFromDt==0) hitStaProvenance->Fill(6);
     if(hitsFromDt!=0 && hitsFromCsc!=0 && hitsFromRpc!=0) hitStaProvenance->Fill(7);
 
     if(hitsFromSegmDt+hitsFromSegmCsc !=0){
       trackHitPercentualVsEta->Fill(recoTrack->eta(), double(hitsFromDt+hitsFromCsc)/(hitsFromSegmDt+hitsFromSegmCsc));
       trackHitPercentualVsPhi->Fill(recoTrack->phi(), double(hitsFromDt+hitsFromCsc)/(hitsFromSegmDt+hitsFromSegmCsc));
       trackHitPercentualVsPt->Fill(recoTrack->pt(), double(hitsFromDt+hitsFromCsc)/(hitsFromSegmDt+hitsFromSegmCsc));
     }
 
     if(hitsFromSegmDt!=0){
       dtTrackHitPercentualVsEta->Fill(recoTrack->eta(), double(hitsFromDt)/hitsFromSegmDt);
       dtTrackHitPercentualVsPhi->Fill(recoTrack->phi(), double(hitsFromDt)/hitsFromSegmDt);
       dtTrackHitPercentualVsPt->Fill(recoTrack->pt(), double(hitsFromDt)/hitsFromSegmDt);
     }
 
     if(hitsFromSegmCsc!=0){
       cscTrackHitPercentualVsEta->Fill(recoTrack->eta(), double(hitsFromCsc)/hitsFromSegmCsc);
       cscTrackHitPercentualVsPhi->Fill(recoTrack->phi(), double(hitsFromCsc)/hitsFromSegmCsc);
       cscTrackHitPercentualVsPt->Fill(recoTrack->pt(), double(hitsFromCsc)/hitsFromSegmCsc);
     }
   }
 
 } 

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

overridevirtual

Implements thread_unsafe::DQMEDAnalyzer.

Definition at line 54 of file SegmentTrackAnalyzer.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::cd(), PtMinSelector_cfg::ptMin, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), and DQMStore::IBooker::setCurrentFolder().

                                              {
   ibooker.cd();
   ibooker.setCurrentFolder("Muons/SegmentTrackAnalyzer");
   
   
   // histograms initalization
   hitsNotUsed = ibooker.book1D("HitsNotUsedForGlobalTracking_"+trackCollection, "recHits not used for GLB    ["+trackCollection+"]", 50, -0.5, 49.5);
   hitsNotUsedPercentual = ibooker.book1D("HitsNotUsedForGlobalTrackingDvHitUsed_"+trackCollection, "(recHits_{notUsedForGLB}) / (recHits_{GLB})    ["+trackCollection+"]", 100, 0, 1.);
 
   TrackSegm = ibooker.book2D("trackSegments_"+trackCollection, "Number of segments associated to the track    ["+trackCollection+"]", 3, 0.5, 3.5, 8, 0, 8);
   TrackSegm->setBinLabel(1,"DT+CSC",1);
   TrackSegm->setBinLabel(2,"DT",1);
   TrackSegm->setBinLabel(3,"CSC",1);
   
   hitStaProvenance = ibooker.book1D("trackHitStaProvenance_"+trackCollection, "Number of recHits_{STAinTrack}    ["+trackCollection+"]", 7, 0.5, 7.5);
   hitStaProvenance->setBinLabel(1,"DT");
   hitStaProvenance->setBinLabel(2,"CSC");
   hitStaProvenance->setBinLabel(3,"RPC");
   hitStaProvenance->setBinLabel(4,"DT+CSC");
   hitStaProvenance->setBinLabel(5,"DT+RPC");
   hitStaProvenance->setBinLabel(6,"CSC+RPC");
   hitStaProvenance->setBinLabel(7,"DT+CSC+RPC");
 
 
   if(trackCollection!="standAloneMuons"){
     hitTkrProvenance = ibooker.book1D("trackHitTkrProvenance_"+trackCollection, "Number of recHits_{TKinTrack}    ["+trackCollection+"]", 6, 0.5, 6.5);
     hitTkrProvenance->setBinLabel(1,"PixBarrel");
     hitTkrProvenance->setBinLabel(2,"PixEndCap");
     hitTkrProvenance->setBinLabel(3,"TIB");
     hitTkrProvenance->setBinLabel(4,"TID");
     hitTkrProvenance->setBinLabel(5,"TOB");
     hitTkrProvenance->setBinLabel(6,"TEC");
   }
 
   trackHitPercentualVsEta = ibooker.book2D("trackHitDivtrackSegmHitVsEta_"+trackCollection, "(recHits_{Track} / recHits_{associatedSegm}) vs #eta    [" +trackCollection+"]", etaBin, etaMin, etaMax, 20, 0, 1);
   dtTrackHitPercentualVsEta = ibooker.book2D("dtTrackHitDivtrackSegmHitVsEta_"+trackCollection, "(recHits_{DTinTrack} / recHits_{associatedSegm}) vs #eta    [" +trackCollection+"]", etaBin, etaMin, etaMax, 20, 0, 1);
   cscTrackHitPercentualVsEta = ibooker.book2D("cscTrackHitDivtrackSegmHitVsEta_"+trackCollection, "(recHits_{CSCinTrack} / recHits_{associatedSegm}) vs #eta    [" +trackCollection+"]", etaBin, etaMin, etaMax, 20, 0, 1);
 
   trackHitPercentualVsPhi = ibooker.book2D("trackHitDivtrackSegmHitVsPhi_"+trackCollection, "(recHits_{Track} / recHits_{associatedSegm}) vs #phi    [" +trackCollection+"]", phiBin, phiMin, phiMax, 20, 0, 1);
   trackHitPercentualVsPhi->setAxisTitle("rad",2);
   dtTrackHitPercentualVsPhi = ibooker.book2D("dtTrackHitDivtrackSegmHitVsPhi_"+trackCollection, "(recHits_{DTinTrack} / recHits_{associatedSegm}) vs #phi    [" +trackCollection+"]", phiBin, phiMin, phiMax, 20, 0, 1);
   dtTrackHitPercentualVsPhi->setAxisTitle("rad",2);
   cscTrackHitPercentualVsPhi = ibooker.book2D("cscTrackHitDivtrackSegmHitVsPhi_"+trackCollection, "(recHits_{CSCinTrack} / recHits_{associatedSegm}) vs #phi    [" +trackCollection+"]", phiBin, phiMin, phiMax, 20, 0, 1);
   cscTrackHitPercentualVsPhi->setAxisTitle("rad",2);
 
   trackHitPercentualVsPt = ibooker.book2D("trackHitDivtrackSegmHitVsPt_"+trackCollection, "(recHits_{Track} / recHits_{associatedSegm}) vs 1/p_{t}    [" +trackCollection+"]", ptBin, ptMin, ptMax, 20, 0, 1);
   trackHitPercentualVsPt->setAxisTitle("GeV",2);
   dtTrackHitPercentualVsPt = ibooker.book2D("dtTrackHitDivtrackSegmHitVsPt_"+trackCollection, "(recHits_{DTinTrack} / recHits_{associatedSegm}) vs 1/p_{t}    [" +trackCollection+"]", ptBin, ptMin, ptMax, 20, 0, 1);
   dtTrackHitPercentualVsPt->setAxisTitle("GeV",2);
   cscTrackHitPercentualVsPt = ibooker.book2D("cscTrackHitDivtrackSegmHitVsPt_"+trackCollection, "(recHits_{CSCinTrack} / recHits_{associatedSegm}) vs 1/p_{t}    [" +trackCollection+"]", ptBin, ptMin, ptMax, 20, 0, 1);
   cscTrackHitPercentualVsPt->setAxisTitle("GeV",2);
 }

Member Data Documentation

MonitorElement* SegmentTrackAnalyzer::cscTrackHitPercentualVsEta

private

Definition at line 78 of file SegmentTrackAnalyzer.h.

MonitorElement* SegmentTrackAnalyzer::cscTrackHitPercentualVsPhi

private

Definition at line 79 of file SegmentTrackAnalyzer.h.

MonitorElement* SegmentTrackAnalyzer::cscTrackHitPercentualVsPt

private

Definition at line 80 of file SegmentTrackAnalyzer.h.

MonitorElement* SegmentTrackAnalyzer::dtTrackHitPercentualVsEta

private

Definition at line 75 of file SegmentTrackAnalyzer.h.

MonitorElement* SegmentTrackAnalyzer::dtTrackHitPercentualVsPhi

private

Definition at line 76 of file SegmentTrackAnalyzer.h.

MonitorElement* SegmentTrackAnalyzer::dtTrackHitPercentualVsPt

private

Definition at line 77 of file SegmentTrackAnalyzer.h.

int SegmentTrackAnalyzer::etaBin

private

Definition at line 56 of file SegmentTrackAnalyzer.h.

double SegmentTrackAnalyzer::etaMax

private

Definition at line 58 of file SegmentTrackAnalyzer.h.

double SegmentTrackAnalyzer::etaMin

private

Definition at line 57 of file SegmentTrackAnalyzer.h.

MonitorElement* SegmentTrackAnalyzer::hitsNotUsed

private

Definition at line 67 of file SegmentTrackAnalyzer.h.

MonitorElement* SegmentTrackAnalyzer::hitsNotUsedPercentual

private

Definition at line 68 of file SegmentTrackAnalyzer.h.

MonitorElement* SegmentTrackAnalyzer::hitStaProvenance

private

Definition at line 70 of file SegmentTrackAnalyzer.h.

MonitorElement* SegmentTrackAnalyzer::hitTkrProvenance

private

Definition at line 71 of file SegmentTrackAnalyzer.h.

std::string SegmentTrackAnalyzer::metname

private

Definition at line 51 of file SegmentTrackAnalyzer.h.

edm::ParameterSet SegmentTrackAnalyzer::parameters

private

Definition at line 47 of file SegmentTrackAnalyzer.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties().

int SegmentTrackAnalyzer::phiBin

private

Definition at line 59 of file SegmentTrackAnalyzer.h.

double SegmentTrackAnalyzer::phiMax

private

Definition at line 61 of file SegmentTrackAnalyzer.h.

double SegmentTrackAnalyzer::phiMin

private

Definition at line 60 of file SegmentTrackAnalyzer.h.

int SegmentTrackAnalyzer::ptBin

private

Definition at line 62 of file SegmentTrackAnalyzer.h.

double SegmentTrackAnalyzer::ptMax

private

Definition at line 64 of file SegmentTrackAnalyzer.h.

double SegmentTrackAnalyzer::ptMin

private

Definition at line 63 of file SegmentTrackAnalyzer.h.

edm::EDGetTokenT<reco::TrackCollection> SegmentTrackAnalyzer::theMuTrackCollectionLabel_

private

Definition at line 48 of file SegmentTrackAnalyzer.h.

SegmentsTrackAssociator* SegmentTrackAnalyzer::theSegmentsAssociator

private

Definition at line 54 of file SegmentTrackAnalyzer.h.

MuonServiceProxy* SegmentTrackAnalyzer::theService

private

Definition at line 46 of file SegmentTrackAnalyzer.h.

std::string SegmentTrackAnalyzer::trackCollection

private

Definition at line 52 of file SegmentTrackAnalyzer.h.

MonitorElement* SegmentTrackAnalyzer::trackHitPercentualVsEta

private

Definition at line 72 of file SegmentTrackAnalyzer.h.

MonitorElement* SegmentTrackAnalyzer::trackHitPercentualVsPhi

private

Definition at line 73 of file SegmentTrackAnalyzer.h.

MonitorElement* SegmentTrackAnalyzer::trackHitPercentualVsPt

private

Definition at line 74 of file SegmentTrackAnalyzer.h.

MonitorElement* SegmentTrackAnalyzer::TrackSegm

private

Definition at line 69 of file SegmentTrackAnalyzer.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation