CMS 3D CMS Logo

SegmentTrackAnalyzer Class Reference

DQM monitoring source for segments associated to the muon track. More...

#include <DQMOffline/Muon/src/SegmentTrackAnalyzer.h>

Inheritance diagram for SegmentTrackAnalyzer:

MuonAnalyzerBase

List of all members.

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &, const reco::Track &recoTrack)
 Get the analysis.
void beginJob (edm::EventSetup const &iSetup, DQMStore *dbe)
 Inizialize parameters for histo binning.
 SegmentTrackAnalyzer (const edm::ParameterSet &, MuonServiceProxy *theService)
 Constructor.
virtual ~SegmentTrackAnalyzer ()
 Destructor.

Private Attributes

MonitorElement * cscTrackHitPercentualVsEta
MonitorElement * cscTrackHitPercentualVsPhi
MonitorElement * cscTrackHitPercentualVsPt
MonitorElement * dtTrackHitPercentualVsEta
MonitorElement * dtTrackHitPercentualVsPhi
MonitorElement * dtTrackHitPercentualVsPt
MonitorElement * hitsNotUsed
MonitorElement * hitsNotUsedPercentual
MonitorElement * hitStaProvenance
MonitorElement * hitTkrProvenance
std::string metname
edm::ParameterSet parameters
SegmentsTrackAssociator * theSegmentsAssociator
MonitorElement * trackHitPercentualVsEta
MonitorElement * trackHitPercentualVsPhi
MonitorElement * trackHitPercentualVsPt
MonitorElement * TrackSegm

Detailed Description

DQM monitoring source for segments associated to the muon track.

Date
2008/05/23 13:48:11
Revision
1.6
Author:
G. Mila - INFN Torino

Definition at line 29 of file SegmentTrackAnalyzer.h.

Constructor & Destructor Documentation

SegmentTrackAnalyzer::SegmentTrackAnalyzer ( const edm::ParameterSet &  pSet,
MuonServiceProxy *  theService 
)

Constructor.

Definition at line 32 of file SegmentTrackAnalyzer.cc.

References edm::ParameterSet::getParameter(), parameters, and theSegmentsAssociator.

00032                                                                                                    :MuonAnalyzerBase(theService) {
00033 
00034   parameters = pSet;
00035 
00036   const ParameterSet SegmentsTrackAssociatorParameters = parameters.getParameter<ParameterSet>("SegmentsTrackAssociatorParameters");
00037   theSegmentsAssociator = new SegmentsTrackAssociator(SegmentsTrackAssociatorParameters);
00038 
00039 }

SegmentTrackAnalyzer::~SegmentTrackAnalyzer (  )  [virtual]

Destructor.

Definition at line 42 of file SegmentTrackAnalyzer.cc.

00042 { }

Member Function Documentation

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

Get the analysis.

Definition at line 112 of file SegmentTrackAnalyzer.cc.

References SegmentsTrackAssociator::associate(), MuonSubdetId::CSC, cscTrackHitPercentualVsEta, cscTrackHitPercentualVsPhi, cscTrackHitPercentualVsPt, MuonSubdetId::DT, dtTrackHitPercentualVsEta, dtTrackHitPercentualVsPhi, dtTrackHitPercentualVsPt, reco::TrackBase::eta(), MonitorElement::Fill(), hitsNotUsed, hitsNotUsedPercentual, hitStaProvenance, hitTkrProvenance, LogTrace, metname, DetId::Muon, reco::TrackBase::phi(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, reco::TrackBase::pt(), reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), MuonSubdetId::RPC, size, SiStripDetId::TEC, theSegmentsAssociator, SiStripDetId::TIB, SiStripDetId::TID, SiStripDetId::TOB, DetId::Tracker, trackHitPercentualVsEta, trackHitPercentualVsPhi, trackHitPercentualVsPt, and TrackSegm.

Referenced by MuonAnalyzer::analyze().

00112                                                                                                                  {
00113 
00114   LogTrace(metname)<<"[SegmentTrackAnalyzer] Filling the histos";
00115   
00116   MuonTransientTrackingRecHit::MuonRecHitContainer segments = theSegmentsAssociator->associate(iEvent, iSetup, recoTrack );
00117  
00118   LogTrace(metname)<<"[SegmentTrackAnalyzer] # of segments associated to the track: "<<(segments).size();
00119 
00120   // hit counters
00121   int hitsFromDt=0;
00122   int hitsFromCsc=0;
00123   int hitsFromRpc=0;
00124   int hitsFromTk=0;
00125   int hitsFromTrack=0;
00126   int hitsFromSegmDt=0;
00127   int hitsFromSegmCsc=0;
00128   // segment counters
00129   int segmFromDt=0;
00130   int segmFromCsc=0;
00131 
00132   for (MuonTransientTrackingRecHit::MuonRecHitContainer::const_iterator segment=segments.begin();
00133        segment!=segments.end(); segment++) {
00134    
00135     DetId id = (*segment)->geographicalId();
00136     
00137     // hits from DT segments
00138     if (id.det() == DetId::Muon && id.subdetId() == MuonSubdetId::DT ) {
00139       ++segmFromDt;
00140       const DTRecSegment4D *seg4D = dynamic_cast<const DTRecSegment4D*>((*segment)->hit());
00141       if((*seg4D).hasPhi())
00142         hitsFromSegmDt+=(*seg4D).phiSegment()->specificRecHits().size();
00143       if((*seg4D).hasZed())
00144         hitsFromSegmDt+=(*seg4D).zSegment()->specificRecHits().size();
00145       
00146     }
00147     
00148     // hits from CSC segments
00149     if (id.det() == DetId::Muon && id.subdetId() == MuonSubdetId::CSC ) {
00150       hitsFromSegmCsc+=(*segment)->recHits().size();
00151       segmFromCsc++;
00152     }
00153 
00154   }
00155 
00156 
00157   // hits from track
00158   for(trackingRecHit_iterator recHit =  recoTrack.recHitsBegin(); recHit != recoTrack.recHitsEnd(); ++recHit){
00159 
00160     hitsFromTrack++;
00161      DetId id = (*recHit)->geographicalId();
00162      // hits from DT
00163      if (id.det() == DetId::Muon && id.subdetId() == MuonSubdetId::DT ) 
00164        hitsFromDt++;   
00165      // hits from CSC
00166       if (id.det() == DetId::Muon && id.subdetId() == MuonSubdetId::CSC ) 
00167        hitsFromCsc++;
00168      // hits from RPC
00169      if (id.det() == DetId::Muon && id.subdetId() == MuonSubdetId::RPC ) 
00170        hitsFromRpc++;
00171      // hits from Tracker
00172      if (id.det() == DetId::Tracker){
00173        hitsFromTk++;
00174        if(id.subdetId() == PixelSubdetector::PixelBarrel )
00175          hitTkrProvenance->Fill(1);
00176        if(id.subdetId() == PixelSubdetector::PixelEndcap )
00177          hitTkrProvenance->Fill(2);
00178        if(id.subdetId() == SiStripDetId::TIB )
00179          hitTkrProvenance->Fill(3);
00180        if(id.subdetId() == SiStripDetId::TID )
00181          hitTkrProvenance->Fill(4);
00182        if(id.subdetId() == SiStripDetId::TOB )
00183          hitTkrProvenance->Fill(5);
00184        if(id.subdetId() == SiStripDetId::TEC )
00185          hitTkrProvenance->Fill(6);
00186      }
00187 
00188   }
00189 
00190   // fill the histos
00191   hitsNotUsed->Fill(hitsFromSegmDt+hitsFromSegmCsc+hitsFromRpc+hitsFromTk-hitsFromTrack);
00192   hitsNotUsedPercentual->Fill(double(hitsFromSegmDt+hitsFromSegmCsc+hitsFromRpc+hitsFromTk-hitsFromTrack)/hitsFromTrack);
00193 
00194   if(hitsFromDt!=0 && hitsFromCsc!=0)
00195     TrackSegm->Fill(1,segmFromDt+segmFromCsc);
00196   if(hitsFromDt!=0 && hitsFromCsc==0)
00197     TrackSegm->Fill(2,segmFromDt);
00198   if(hitsFromDt==0 && hitsFromCsc!=0)
00199     TrackSegm->Fill(3,segmFromCsc);
00200 
00201   if(hitsFromDt!=0 && hitsFromCsc==0 && hitsFromRpc==0) hitStaProvenance->Fill(1);
00202   if(hitsFromCsc!=0 && hitsFromDt==0 && hitsFromRpc==0) hitStaProvenance->Fill(2);
00203   if(hitsFromRpc!=0 && hitsFromDt==0 && hitsFromCsc==0) hitStaProvenance->Fill(3);
00204   if(hitsFromDt!=0 && hitsFromCsc!=0 && hitsFromRpc==0) hitStaProvenance->Fill(4);
00205   if(hitsFromDt!=0 && hitsFromRpc!=0 && hitsFromCsc==0) hitStaProvenance->Fill(5);
00206   if(hitsFromCsc!=0 && hitsFromRpc!=0 && hitsFromDt==0) hitStaProvenance->Fill(6);
00207   if(hitsFromDt!=0 && hitsFromCsc!=0 && hitsFromRpc!=0) hitStaProvenance->Fill(7);
00208 
00209   if(hitsFromSegmDt+hitsFromSegmCsc !=0){
00210     trackHitPercentualVsEta->Fill(recoTrack.eta(), double(hitsFromDt+hitsFromCsc)/(hitsFromSegmDt+hitsFromSegmCsc));
00211     trackHitPercentualVsPhi->Fill(recoTrack.phi(), double(hitsFromDt+hitsFromCsc)/(hitsFromSegmDt+hitsFromSegmCsc));
00212     trackHitPercentualVsPt->Fill(recoTrack.pt(), double(hitsFromDt+hitsFromCsc)/(hitsFromSegmDt+hitsFromSegmCsc));
00213   }
00214 
00215   if(hitsFromSegmDt!=0){
00216     dtTrackHitPercentualVsEta->Fill(recoTrack.eta(), double(hitsFromDt)/hitsFromSegmDt);
00217     dtTrackHitPercentualVsPhi->Fill(recoTrack.phi(), double(hitsFromDt)/hitsFromSegmDt);
00218     dtTrackHitPercentualVsPt->Fill(recoTrack.pt(), double(hitsFromDt)/hitsFromSegmDt);
00219   }
00220 
00221   if(hitsFromSegmCsc!=0){
00222     cscTrackHitPercentualVsEta->Fill(recoTrack.eta(), double(hitsFromCsc)/hitsFromSegmCsc);
00223     cscTrackHitPercentualVsPhi->Fill(recoTrack.phi(), double(hitsFromCsc)/hitsFromSegmCsc);
00224     cscTrackHitPercentualVsPt->Fill(recoTrack.pt(), double(hitsFromCsc)/hitsFromSegmCsc);
00225   }
00226 
00227 }

void SegmentTrackAnalyzer::beginJob ( edm::EventSetup const &  iSetup,
DQMStore *  dbe 
) [virtual]

Inizialize parameters for histo binning.

Implements MuonAnalyzerBase.

Definition at line 45 of file SegmentTrackAnalyzer.cc.

References DQMStore::book1D(), DQMStore::book2D(), cscTrackHitPercentualVsEta, cscTrackHitPercentualVsPhi, cscTrackHitPercentualVsPt, dtTrackHitPercentualVsEta, dtTrackHitPercentualVsPhi, dtTrackHitPercentualVsPt, gen_hydjet_hard_muon_trigger_cfg::etaMax, gen_hydjet_hard_muon_trigger_cfg::etaMin, edm::ParameterSet::getParameter(), hitsNotUsed, hitsNotUsedPercentual, hitStaProvenance, hitTkrProvenance, label, LogTrace, metname, parameters, RegionalCKFTracksForL3Isolation_cfi::ptMin, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), trackHitPercentualVsEta, trackHitPercentualVsPhi, trackHitPercentualVsPt, and TrackSegm.

Referenced by MuonAnalyzer::beginJob().

00045                                                                               {
00046 
00047 
00048   metname = "segmTrackAnalyzer";
00049   string trackCollection = parameters.getParameter<edm::InputTag>("MuTrackCollection").label() + parameters.getParameter<edm::InputTag>("MuTrackCollection").instance();
00050   LogTrace(metname)<<"[SegmentTrackAnalyzer] Parameters initialization";
00051   dbe->setCurrentFolder("Muons/SegmentTrackAnalyzer");
00052   
00053   // histograms initalization
00054   hitsNotUsed = dbe->book1D("HitsNotUsedForGlobalTracking_"+trackCollection, "recHits not used for GLB    ["+trackCollection+"]", 50, -0.5, 49.5);
00055   hitsNotUsedPercentual = dbe->book1D("HitsNotUsedForGlobalTrackingDvHitUsed_"+trackCollection, "(recHits_{notUsedForGLB}) / (recHits_{GLB})    ["+trackCollection+"]", 100, 0, 1.);
00056 
00057   TrackSegm = dbe->book2D("trackSegments_"+trackCollection, "Number of segments associated to the track    ["+trackCollection+"]", 3, 0.5, 3.5, 8, 0, 8);
00058   TrackSegm->setBinLabel(1,"DT+CSC",1);
00059   TrackSegm->setBinLabel(2,"DT",1);
00060   TrackSegm->setBinLabel(3,"CSC",1);
00061   
00062   hitStaProvenance = dbe->book1D("trackHitStaProvenance_"+trackCollection, "Number of recHits_{STAinTrack}    ["+trackCollection+"]", 7, 0.5, 7.5);
00063   hitStaProvenance->setBinLabel(1,"DT");
00064   hitStaProvenance->setBinLabel(2,"CSC");
00065   hitStaProvenance->setBinLabel(3,"RPC");
00066   hitStaProvenance->setBinLabel(4,"DT+CSC");
00067   hitStaProvenance->setBinLabel(5,"DT+RPC");
00068   hitStaProvenance->setBinLabel(6,"CSC+RPC");
00069   hitStaProvenance->setBinLabel(7,"DT+CSC+RPC");
00070 
00071 
00072   if(trackCollection!="standAloneMuons"){
00073     hitTkrProvenance = dbe->book1D("trackHitTkrProvenance_"+trackCollection, "Number of recHits_{TKinTrack}    ["+trackCollection+"]", 6, 0.5, 6.5);
00074     hitTkrProvenance->setBinLabel(1,"PixBarrel");
00075     hitTkrProvenance->setBinLabel(2,"PixEndCap");
00076     hitTkrProvenance->setBinLabel(3,"TIB");
00077     hitTkrProvenance->setBinLabel(4,"TID");
00078     hitTkrProvenance->setBinLabel(5,"TOB");
00079     hitTkrProvenance->setBinLabel(6,"TEC");
00080   }
00081 
00082   int etaBin = parameters.getParameter<int>("etaBin");
00083   double etaMin = parameters.getParameter<double>("etaMin");
00084   double etaMax = parameters.getParameter<double>("etaMax");
00085   trackHitPercentualVsEta = dbe->book2D("trackHitDivtrackSegmHitVsEta_"+trackCollection, "(recHits_{Track} / recHits_{associatedSegm}) vs #eta    [" +trackCollection+"]", etaBin, etaMin, etaMax, 20, 0, 1);
00086   dtTrackHitPercentualVsEta = dbe->book2D("dtTrackHitDivtrackSegmHitVsEta_"+trackCollection, "(recHits_{DTinTrack} / recHits_{associatedSegm}) vs #eta    [" +trackCollection+"]", etaBin, etaMin, etaMax, 20, 0, 1);
00087   cscTrackHitPercentualVsEta = dbe->book2D("cscTrackHitDivtrackSegmHitVsEta_"+trackCollection, "(recHits_{CSCinTrack} / recHits_{associatedSegm}) vs #eta    [" +trackCollection+"]", etaBin, etaMin, etaMax, 20, 0, 1);
00088 
00089   int phiBin = parameters.getParameter<int>("phiBin");
00090   double phiMin = parameters.getParameter<double>("phiMin");
00091   double phiMax = parameters.getParameter<double>("phiMax");
00092   trackHitPercentualVsPhi = dbe->book2D("trackHitDivtrackSegmHitVsPhi_"+trackCollection, "(recHits_{Track} / recHits_{associatedSegm}) vs #phi    [" +trackCollection+"]", phiBin, phiMin, phiMax, 20, 0, 1);
00093   trackHitPercentualVsPhi->setAxisTitle("rad",2);
00094   dtTrackHitPercentualVsPhi = dbe->book2D("dtTrackHitDivtrackSegmHitVsPhi_"+trackCollection, "(recHits_{DTinTrack} / recHits_{associatedSegm}) vs #phi    [" +trackCollection+"]", phiBin, phiMin, phiMax, 20, 0, 1);
00095   dtTrackHitPercentualVsPhi->setAxisTitle("rad",2);
00096   cscTrackHitPercentualVsPhi = dbe->book2D("cscTrackHitDivtrackSegmHitVsPhi_"+trackCollection, "(recHits_{CSCinTrack} / recHits_{associatedSegm}) vs #phi    [" +trackCollection+"]", phiBin, phiMin, phiMax, 20, 0, 1);
00097   cscTrackHitPercentualVsPhi->setAxisTitle("rad",2);
00098 
00099   int ptBin = parameters.getParameter<int>("ptBin");
00100   double ptMin = parameters.getParameter<double>("ptMin");
00101   double ptMax = parameters.getParameter<double>("ptMax");
00102   trackHitPercentualVsPt = dbe->book2D("trackHitDivtrackSegmHitVsPt_"+trackCollection, "(recHits_{Track} / recHits_{associatedSegm}) vs 1/p_{t}    [" +trackCollection+"]", ptBin, ptMin, ptMax, 20, 0, 1);
00103   trackHitPercentualVsPt->setAxisTitle("GeV",2);
00104   dtTrackHitPercentualVsPt = dbe->book2D("dtTrackHitDivtrackSegmHitVsPt_"+trackCollection, "(recHits_{DTinTrack} / recHits_{associatedSegm}) vs 1/p_{t}    [" +trackCollection+"]", ptBin, ptMin, ptMax, 20, 0, 1);
00105   dtTrackHitPercentualVsPt->setAxisTitle("GeV",2);
00106   cscTrackHitPercentualVsPt = dbe->book2D("cscTrackHitDivtrackSegmHitVsPt_"+trackCollection, "(recHits_{CSCinTrack} / recHits_{associatedSegm}) vs 1/p_{t}    [" +trackCollection+"]", ptBin, ptMin, ptMax, 20, 0, 1);
00107   cscTrackHitPercentualVsPt->setAxisTitle("GeV",2);
00108 
00109 }

Member Data Documentation

MonitorElement* SegmentTrackAnalyzer::cscTrackHitPercentualVsEta [private]

Definition at line 66 of file SegmentTrackAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SegmentTrackAnalyzer::cscTrackHitPercentualVsPhi [private]

Definition at line 67 of file SegmentTrackAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SegmentTrackAnalyzer::cscTrackHitPercentualVsPt [private]

Definition at line 68 of file SegmentTrackAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SegmentTrackAnalyzer::dtTrackHitPercentualVsEta [private]

Definition at line 63 of file SegmentTrackAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SegmentTrackAnalyzer::dtTrackHitPercentualVsPhi [private]

Definition at line 64 of file SegmentTrackAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SegmentTrackAnalyzer::dtTrackHitPercentualVsPt [private]

Definition at line 65 of file SegmentTrackAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SegmentTrackAnalyzer::hitsNotUsed [private]

Definition at line 55 of file SegmentTrackAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SegmentTrackAnalyzer::hitsNotUsedPercentual [private]

Definition at line 56 of file SegmentTrackAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SegmentTrackAnalyzer::hitStaProvenance [private]

Definition at line 58 of file SegmentTrackAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SegmentTrackAnalyzer::hitTkrProvenance [private]

Definition at line 59 of file SegmentTrackAnalyzer.h.

Referenced by analyze(), and beginJob().

std::string SegmentTrackAnalyzer::metname [private]

Definition at line 50 of file SegmentTrackAnalyzer.h.

Referenced by analyze(), and beginJob().

edm::ParameterSet SegmentTrackAnalyzer::parameters [private]

Definition at line 48 of file SegmentTrackAnalyzer.h.

Referenced by beginJob(), and SegmentTrackAnalyzer().

SegmentsTrackAssociator* SegmentTrackAnalyzer::theSegmentsAssociator [private]

Definition at line 52 of file SegmentTrackAnalyzer.h.

Referenced by analyze(), and SegmentTrackAnalyzer().

MonitorElement* SegmentTrackAnalyzer::trackHitPercentualVsEta [private]

Definition at line 60 of file SegmentTrackAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SegmentTrackAnalyzer::trackHitPercentualVsPhi [private]

Definition at line 61 of file SegmentTrackAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SegmentTrackAnalyzer::trackHitPercentualVsPt [private]

Definition at line 62 of file SegmentTrackAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* SegmentTrackAnalyzer::TrackSegm [private]

Definition at line 57 of file SegmentTrackAnalyzer.h.

Referenced by analyze(), and beginJob().

The documentation for this class was generated from the following files:
Generated on Tue Jun 9 18:31:16 2009 for CMSSW by  doxygen 1.5.4