CMS 3D CMS Logo
List of all members | Public Member Functions | Private Attributes
PhysicsObjectsMonitor Class Reference

#include <PhysicsObjectsMonitor.h>

Inheritance diagram for PhysicsObjectsMonitor:
one::DQMEDAnalyzer< T > one::dqmimplementation::DQMBaseClass< T... >

Public Member Functions

void analyze (const edm::Event &event, const edm::EventSetup &eventSetup) override
 
void bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
 
 PhysicsObjectsMonitor (const edm::ParameterSet &pset)
 Constructor. More...
 
 ~PhysicsObjectsMonitor () override
 Destructor. More...
 
- Public Member Functions inherited from one::DQMEDAnalyzer< T >
 DQMEDAnalyzer ()=default
 
 DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
 
 DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
 
 ~DQMEDAnalyzer () override=default
 

Private Attributes

MonitorElement * charge
 
MonitorElement * CSCvsRPC
 
MonitorElement * DTvsCSC
 
MonitorElement * DTvsRPC
 
MonitorElement * h1_Pres
 
MonitorElement * hPres
 
MonitorElement * NCSChits
 
MonitorElement * NDThits
 
MonitorElement * Nmuon
 
MonitorElement * Nrechits
 
MonitorElement * NRPChits
 
MonitorElement * pt
 
MonitorElement * ptot
 
MonitorElement * px
 
MonitorElement * py
 
MonitorElement * pz
 
std::string theDataType
 
std::string theSeedCollectionLabel
 
std::string theSTAMuonLabel
 
edm::EDGetTokenT< reco::TrackCollection > theSTAMuonToken_
 

Detailed Description

For now: Analyzer of StandAlone muon tracks Later: Add other detectors and more Reco

Author
M. Mulders - CERN marti.nosp@m.jn.m.nosp@m.ulder.nosp@m.s@ce.nosp@m.rn.ch Based on STAMuonAnalyzer by R. Bellan - INFN Torino ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch

Analyzer of the StandAlone muon tracks

Author
M. Mulders - CERN marti.nosp@m.jn.m.nosp@m.ulder.nosp@m.s@ce.nosp@m.rn.ch Based on STAMuonAnalyzer by R. Bellan - INFN Torino ricca.nosp@m.rdo..nosp@m.bella.nosp@m.n@ce.nosp@m.rn.ch

Definition at line 31 of file PhysicsObjectsMonitor.h.

Constructor & Destructor Documentation

PhysicsObjectsMonitor::PhysicsObjectsMonitor ( const edm::ParameterSet &  pset)

Constructor.

Definition at line 32 of file PhysicsObjectsMonitor.cc.

References edm::ParameterSet::getUntrackedParameter().

32  {
33  theSTAMuonLabel = pset.getUntrackedParameter<string>("StandAloneTrackCollectionLabel");
34  theSeedCollectionLabel = pset.getUntrackedParameter<string>("MuonSeedCollectionLabel");
35  theDataType = pset.getUntrackedParameter<string>("DataType");
36 
37  if (theDataType != "RealData" && theDataType != "SimData")
38  edm::LogInfo("PhysicsObjectsMonitor") << "Error in Data Type!!" << endl;
39 
40  if (theDataType == "SimData") {
41  edm::LogInfo("PhysicsObjectsMonitor") << "Sorry! Running this package on simulation is no longer supported! ";
42  }
43 
44  // set Token(-s)
45  theSTAMuonToken_ =
46  consumes<reco::TrackCollection>(pset.getUntrackedParameter<string>("StandAloneTrackCollectionLabel"));
47 }
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
PhysicsObjectsMonitor::theSTAMuonToken_
edm::EDGetTokenT< reco::TrackCollection > theSTAMuonToken_
Definition: PhysicsObjectsMonitor.h:68
PhysicsObjectsMonitor::~PhysicsObjectsMonitor ( )
override

Destructor.

Definition at line 50 of file PhysicsObjectsMonitor.cc.

50 {}

Member Function Documentation

void PhysicsObjectsMonitor::analyze ( const edm::Event &  event,
const edm::EventSetup &  eventSetup 
)
override

Definition at line 76 of file PhysicsObjectsMonitor.cc.

References ALCARECOTkAlJpsiMuMu_cff::charge, TrajectoryStateClosestToPoint::charge(), MuonSubdetId::CSC, debug, DEFINE_FWK_MODULE, MuonSubdetId::DT, MuonPatternRecoDumper::dumpFTS(), edm::EventSetup::get(), reco::TransientTrack::impactPointTSCP(), PV3DBase< T, PVType, FrameType >::mag(), TrajectoryStateClosestToPoint::momentum(), PV3DBase< T, PVType, FrameType >::perp(), EnergyCorrector::pt, reco::TransientTrack::recHitsBegin(), reco::TransientTrack::recHitsEnd(), recPt, MuonSubdetId::RPC, TrajectoryStateClosestToPoint::theState(), HiIsolationCommonParameters_cff::track, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

76  {
77  edm::LogInfo("PhysicsObjectsMonitor") << "Run: " << event.id().run() << " Event: " << event.id().event();
78  MuonPatternRecoDumper debug;
79 
80  // Get the RecTrack collection from the event
81  Handle<reco::TrackCollection> staTracks;
82  event.getByToken(theSTAMuonToken_, staTracks);
83 
84  ESHandle<MagneticField> theMGField;
85  eventSetup.get<IdealMagneticFieldRecord>().get(theMGField);
86 
87  double recPt = 0.;
88  double simPt = 0.;
89 
90  reco::TrackCollection::const_iterator staTrack;
91 
92  edm::LogInfo("PhysicsObjectsMonitor") << "Reconstructed tracks: " << staTracks->size() << endl;
93  Nmuon->Fill(staTracks->size());
94  for (staTrack = staTracks->begin(); staTrack != staTracks->end(); ++staTrack) {
95  reco::TransientTrack track(*staTrack, &*theMGField);
96 
97  int nrechits = 0;
98  int nDThits = 0;
99  int nCSChits = 0;
100  int nRPChits = 0;
101 
102  for (trackingRecHit_iterator it = track.recHitsBegin(); it != track.recHitsEnd(); it++) {
103  if ((*it)->isValid()) {
104  edm::LogInfo("PhysicsObjectsMonitor") << "Analyzer: Aha this looks like a Rechit!" << std::endl;
105  if ((*it)->geographicalId().subdetId() == MuonSubdetId::DT) {
106  nDThits++;
107  } else if ((*it)->geographicalId().subdetId() == MuonSubdetId::CSC) {
108  nCSChits++;
109  } else if ((*it)->geographicalId().subdetId() == MuonSubdetId::RPC) {
110  nRPChits++;
111  } else {
112  edm::LogInfo("PhysicsObjectsMonitor") << "This is an UNKNOWN hit !! " << std::endl;
113  }
114  nrechits++;
115  }
116  }
117 
118  Nrechits->Fill(nrechits);
119  NDThits->Fill(nDThits);
120  NCSChits->Fill(nCSChits);
121  DTvsCSC->Fill(nDThits, nCSChits);
122  NRPChits->Fill(nRPChits);
123 
124  debug.dumpFTS(track.impactPointTSCP().theState());
125 
126  recPt = track.impactPointTSCP().momentum().perp();
127  edm::LogInfo("PhysicsObjectsMonitor")
128  << " p: " << track.impactPointTSCP().momentum().mag() << " pT: " << recPt << endl;
129  pt->Fill(recPt);
130  ptot->Fill(track.impactPointTSCP().momentum().mag());
131  charge->Fill(track.impactPointTSCP().charge());
132  px->Fill(track.impactPointTSCP().momentum().x());
133  py->Fill(track.impactPointTSCP().momentum().y());
134  pz->Fill(track.impactPointTSCP().momentum().z());
135  }
136  edm::LogInfo("PhysicsObjectsMonitor") << "---" << endl;
137  if (recPt && theDataType == "SimData") {
138  hPres->Fill((recPt - simPt) / simPt);
139  h1_Pres->Fill((1 / recPt - 1 / simPt) / (1 / simPt));
140  }
141 }
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:166
MuonPatternRecoDumper::dumpFTS
std::string dumpFTS(const FreeTrajectoryState &fts) const
Definition: MuonPatternRecoDumper.cc:49
debug
#define debug
Definition: HDRShower.cc:19
HiIsolationCommonParameters_cff.track
track
Definition: HiIsolationCommonParameters_cff.py:9
MuonSubdetId::RPC
static constexpr int RPC
Definition: MuonSubdetId.h:14
PhysicsObjectsMonitor::theSTAMuonToken_
edm::EDGetTokenT< reco::TrackCollection > theSTAMuonToken_
Definition: PhysicsObjectsMonitor.h:68
edm::EventSetup::get
T get() const
Definition: EventSetup.h:71
edm::OwnVector::const_iterator
Definition: OwnVector.h:41
recPt
double recPt
Definition: PFJetBenchmarkAnalyzer.cc:73
MuonSubdetId::DT
static constexpr int DT
Definition: MuonSubdetId.h:12
MuonSubdetId::CSC
static constexpr int CSC
Definition: MuonSubdetId.h:13
void PhysicsObjectsMonitor::bookHistograms ( DQMStore::IBooker &  iBooker,
edm::Run const &  ,
edm::EventSetup const &   
)
override

Definition at line 52 of file PhysicsObjectsMonitor.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), ALCARECOTkAlJpsiMuMu_cff::charge, MonitorElement::getTH2F(), EnergyCorrector::pt, and DQMStore::IBooker::setCurrentFolder().

52  {
53  iBooker.setCurrentFolder("PhysicsObjects/MuonReconstruction");
54 
55  hPres = iBooker.book1D("pTRes", "pT Resolution", 100, -2, 2);
56  h1_Pres = iBooker.book1D("invPTRes", "1/pT Resolution", 100, -2, 2);
57 
58  charge = iBooker.book1D("charge", "track charge", 5, -2.5, 2.5);
59  ptot = iBooker.book1D("ptot", "track momentum", 50, 0, 50);
60  pt = iBooker.book1D("pt", "track pT", 100, 0, 50);
61  px = iBooker.book1D("px ", "track px", 100, -50, 50);
62  py = iBooker.book1D("py", "track py", 100, -50, 50);
63  pz = iBooker.book1D("pz", "track pz", 100, -50, 50);
64  Nmuon = iBooker.book1D("Nmuon", "Number of muon tracks", 11, -.5, 10.5);
65  Nrechits = iBooker.book1D("Nrechits", "Number of RecHits/Segments on track", 21, -.5, 21.5);
66  NDThits = iBooker.book1D("NDThits", "Number of DT Hits/Segments on track", 31, -.5, 31.5);
67  NCSChits = iBooker.book1D("NCSChits", "Number of CSC Hits/Segments on track", 31, -.5, 31.5);
68  NRPChits = iBooker.book1D("NRPChits", "Number of RPC hits on track", 11, -.5, 11.5);
69 
70  DTvsCSC = iBooker.book2D("DTvsCSC", "Number of DT vs CSC hits on track", 29, -.5, 28.5, 29, -.5, 28.5);
71  TH2F *root_ob = DTvsCSC->getTH2F();
72  root_ob->SetXTitle("Number of DT hits");
73  root_ob->SetYTitle("Number of CSC hits");
74 }
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:268
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:106
MonitorElement::getTH2F
TH2F * getTH2F() const
Definition: MonitorElement.cc:1507
DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:109

Member Data Documentation

MonitorElement* PhysicsObjectsMonitor::charge
private

Definition at line 50 of file PhysicsObjectsMonitor.h.

Referenced by Vispa.Plugins.EdmBrowser.ParticleDataList.ParticleData::__repr__().

MonitorElement* PhysicsObjectsMonitor::CSCvsRPC
private

Definition at line 62 of file PhysicsObjectsMonitor.h.

MonitorElement* PhysicsObjectsMonitor::DTvsCSC
private

Definition at line 60 of file PhysicsObjectsMonitor.h.

MonitorElement* PhysicsObjectsMonitor::DTvsRPC
private

Definition at line 61 of file PhysicsObjectsMonitor.h.

MonitorElement* PhysicsObjectsMonitor::h1_Pres
private

Definition at line 47 of file PhysicsObjectsMonitor.h.

MonitorElement* PhysicsObjectsMonitor::hPres
private

Definition at line 46 of file PhysicsObjectsMonitor.h.

MonitorElement* PhysicsObjectsMonitor::NCSChits
private

Definition at line 59 of file PhysicsObjectsMonitor.h.

MonitorElement* PhysicsObjectsMonitor::NDThits
private

Definition at line 58 of file PhysicsObjectsMonitor.h.

MonitorElement* PhysicsObjectsMonitor::Nmuon
private

Definition at line 56 of file PhysicsObjectsMonitor.h.

MonitorElement* PhysicsObjectsMonitor::Nrechits
private

Definition at line 57 of file PhysicsObjectsMonitor.h.

MonitorElement* PhysicsObjectsMonitor::NRPChits
private

Definition at line 63 of file PhysicsObjectsMonitor.h.

MonitorElement* PhysicsObjectsMonitor::pt
private

Definition at line 52 of file PhysicsObjectsMonitor.h.

Referenced by Particle.Particle::__str__(), ZMuMuRochCorAnalyzer.DiMuon::__str__(), DiObject.DiMuon::__str__(), Photon.Photon::calScaledIsoValueExp(), Photon.Photon::calScaledIsoValueLin(), Photon.Photon::calScaledIsoValueQuadr(), Electron.Electron::mvaIDLoose(), Electron.Electron::mvaIDRun2(), Electron.Electron::mvaIDTight(), Electron.Electron::ptErr(), ntupleDataFormat.Track::ptPull(), Lepton.Lepton::relIso(), Lepton.Lepton::relIsoFromEA(), Lepton.Lepton::relIsoR(), and Jet.Jet::setCorrP4().

MonitorElement* PhysicsObjectsMonitor::ptot
private

Definition at line 51 of file PhysicsObjectsMonitor.h.

MonitorElement* PhysicsObjectsMonitor::px
private

Definition at line 53 of file PhysicsObjectsMonitor.h.

MonitorElement* PhysicsObjectsMonitor::py
private

Definition at line 54 of file PhysicsObjectsMonitor.h.

MonitorElement* PhysicsObjectsMonitor::pz
private

Definition at line 55 of file PhysicsObjectsMonitor.h.

std::string PhysicsObjectsMonitor::theDataType
private

Definition at line 65 of file PhysicsObjectsMonitor.h.

std::string PhysicsObjectsMonitor::theSeedCollectionLabel
private

Definition at line 43 of file PhysicsObjectsMonitor.h.

std::string PhysicsObjectsMonitor::theSTAMuonLabel
private

Definition at line 42 of file PhysicsObjectsMonitor.h.

edm::EDGetTokenT<reco::TrackCollection> PhysicsObjectsMonitor::theSTAMuonToken_
private

Definition at line 68 of file PhysicsObjectsMonitor.h.

Generated for CMSSW Reference Manual by doxygen 1.8.11