PhysicsObjectsMonitor Class Reference

#include <PhysicsObjectsMonitor.h>

Inheritance diagram for PhysicsObjectsMonitor:

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 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
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final

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_

Additional Inherited Members
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 >
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
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

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 30 of file PhysicsObjectsMonitor.h.

Constructor & Destructor Documentation

PhysicsObjectsMonitor()

PhysicsObjectsMonitor::PhysicsObjectsMonitor

(

const edm::ParameterSet &

pset

)

Constructor.

Definition at line 31 of file PhysicsObjectsMonitor.cc.

                                                                     {
  theSTAMuonLabel = pset.getUntrackedParameter<string>("StandAloneTrackCollectionLabel");
  theSeedCollectionLabel = pset.getUntrackedParameter<string>("MuonSeedCollectionLabel");
  theDataType = pset.getUntrackedParameter<string>("DataType");
 
  if (theDataType != "RealData" && theDataType != "SimData")
    edm::LogInfo("PhysicsObjectsMonitor") << "Error in Data Type!!" << endl;
 
  if (theDataType == "SimData") {
    edm::LogInfo("PhysicsObjectsMonitor") << "Sorry! Running this package on simulation is no longer supported! ";
  }
 
  // set Token(-s)
  theSTAMuonToken_ =
      consumes<reco::TrackCollection>(pset.getUntrackedParameter<string>("StandAloneTrackCollectionLabel"));
}

References muonDTDigis_cfi::pset.

~PhysicsObjectsMonitor()

PhysicsObjectsMonitor::~PhysicsObjectsMonitor ( )

override

Destructor.

Definition at line 49 of file PhysicsObjectsMonitor.cc.

{}

Member Function Documentation

analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 75 of file PhysicsObjectsMonitor.cc.

                                                                                    {
  edm::LogInfo("PhysicsObjectsMonitor") << "Run: " << event.id().run() << " Event: " << event.id().event();
  MuonPatternRecoDumper debug;
 
  // Get the RecTrack collection from the event
  Handle<reco::TrackCollection> staTracks;
  event.getByToken(theSTAMuonToken_, staTracks);
 
  ESHandle<MagneticField> theMGField;
  eventSetup.get<IdealMagneticFieldRecord>().get(theMGField);
 
  double recPt = 0.;
  double simPt = 0.;
 
  reco::TrackCollection::const_iterator staTrack;
 
  edm::LogInfo("PhysicsObjectsMonitor") << "Reconstructed tracks: " << staTracks->size() << endl;
  Nmuon->Fill(staTracks->size());
  for (staTrack = staTracks->begin(); staTrack != staTracks->end(); ++staTrack) {
    reco::TransientTrack track(*staTrack, &*theMGField);
 
    int nrechits = 0;
    int nDThits = 0;
    int nCSChits = 0;
    int nRPChits = 0;
 
    for (trackingRecHit_iterator it = track.recHitsBegin(); it != track.recHitsEnd(); it++) {
      if ((*it)->isValid()) {
        edm::LogInfo("PhysicsObjectsMonitor") << "Analyzer:  Aha this looks like a Rechit!" << std::endl;
        if ((*it)->geographicalId().subdetId() == MuonSubdetId::DT) {
          nDThits++;
        } else if ((*it)->geographicalId().subdetId() == MuonSubdetId::CSC) {
          nCSChits++;
        } else if ((*it)->geographicalId().subdetId() == MuonSubdetId::RPC) {
          nRPChits++;
        } else {
          edm::LogInfo("PhysicsObjectsMonitor") << "This is an UNKNOWN hit !! " << std::endl;
        }
        nrechits++;
      }
    }
 
    Nrechits->Fill(nrechits);
    NDThits->Fill(nDThits);
    NCSChits->Fill(nCSChits);
    DTvsCSC->Fill(nDThits, nCSChits);
    NRPChits->Fill(nRPChits);
 
    debug.dumpFTS(track.impactPointTSCP().theState());
 
    recPt = track.impactPointTSCP().momentum().perp();
    edm::LogInfo("PhysicsObjectsMonitor")
        << " p: " << track.impactPointTSCP().momentum().mag() << " pT: " << recPt << endl;
    pt->Fill(recPt);
    ptot->Fill(track.impactPointTSCP().momentum().mag());
    charge->Fill(track.impactPointTSCP().charge());
    px->Fill(track.impactPointTSCP().momentum().x());
    py->Fill(track.impactPointTSCP().momentum().y());
    pz->Fill(track.impactPointTSCP().momentum().z());
  }
  edm::LogInfo("PhysicsObjectsMonitor") << "---" << endl;
  if (recPt && theDataType == "SimData") {
    hPres->Fill((recPt - simPt) / simPt);
    h1_Pres->Fill((1 / recPt - 1 / simPt) / (1 / simPt));
  }
}

References ALCARECOTkAlJpsiMuMu_cff::charge, MuonSubdetId::CSC, debug, MuonSubdetId::DT, edm::EventSetup::get(), get, DiDispStaMuonMonitor_cfi::pt, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, recPt, MuonSubdetId::RPC, muonTagProbeFilters_cff::staTracks, and HLT_FULL_cff::track.

bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 51 of file PhysicsObjectsMonitor.cc.

                                                                                                          {
  iBooker.setCurrentFolder("PhysicsObjects/MuonReconstruction");
 
  hPres = iBooker.book1D("pTRes", "pT Resolution", 100, -2, 2);
  h1_Pres = iBooker.book1D("invPTRes", "1/pT Resolution", 100, -2, 2);
 
  charge = iBooker.book1D("charge", "track charge", 5, -2.5, 2.5);
  ptot = iBooker.book1D("ptot", "track momentum", 50, 0, 50);
  pt = iBooker.book1D("pt", "track pT", 100, 0, 50);
  px = iBooker.book1D("px ", "track px", 100, -50, 50);
  py = iBooker.book1D("py", "track py", 100, -50, 50);
  pz = iBooker.book1D("pz", "track pz", 100, -50, 50);
  Nmuon = iBooker.book1D("Nmuon", "Number of muon tracks", 11, -.5, 10.5);
  Nrechits = iBooker.book1D("Nrechits", "Number of RecHits/Segments on track", 21, -.5, 21.5);
  NDThits = iBooker.book1D("NDThits", "Number of DT Hits/Segments on track", 31, -.5, 31.5);
  NCSChits = iBooker.book1D("NCSChits", "Number of CSC Hits/Segments on track", 31, -.5, 31.5);
  NRPChits = iBooker.book1D("NRPChits", "Number of RPC hits on track", 11, -.5, 11.5);
 
  DTvsCSC = iBooker.book2D("DTvsCSC", "Number of DT vs CSC hits on track", 29, -.5, 28.5, 29, -.5, 28.5);
  TH2F *root_ob = DTvsCSC->getTH2F();
  root_ob->SetXTitle("Number of DT hits");
  root_ob->SetYTitle("Number of CSC hits");
}

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), ALCARECOTkAlJpsiMuMu_cff::charge, dqm::legacy::MonitorElement::getTH2F(), DiDispStaMuonMonitor_cfi::pt, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, and dqm::implementation::NavigatorBase::setCurrentFolder().

Member Data Documentation

charge

MonitorElement* PhysicsObjectsMonitor::charge

private

Definition at line 49 of file PhysicsObjectsMonitor.h.

CSCvsRPC

MonitorElement* PhysicsObjectsMonitor::CSCvsRPC

private

Definition at line 61 of file PhysicsObjectsMonitor.h.

DTvsCSC

MonitorElement* PhysicsObjectsMonitor::DTvsCSC

private

Definition at line 59 of file PhysicsObjectsMonitor.h.

DTvsRPC

MonitorElement* PhysicsObjectsMonitor::DTvsRPC

private

Definition at line 60 of file PhysicsObjectsMonitor.h.

h1_Pres

MonitorElement* PhysicsObjectsMonitor::h1_Pres

private

Definition at line 46 of file PhysicsObjectsMonitor.h.

hPres

MonitorElement* PhysicsObjectsMonitor::hPres

private

Definition at line 45 of file PhysicsObjectsMonitor.h.

NCSChits

MonitorElement* PhysicsObjectsMonitor::NCSChits

private

Definition at line 58 of file PhysicsObjectsMonitor.h.

NDThits

MonitorElement* PhysicsObjectsMonitor::NDThits

private

Definition at line 57 of file PhysicsObjectsMonitor.h.

Nmuon

MonitorElement* PhysicsObjectsMonitor::Nmuon

private

Definition at line 55 of file PhysicsObjectsMonitor.h.

Nrechits

MonitorElement* PhysicsObjectsMonitor::Nrechits

private

Definition at line 56 of file PhysicsObjectsMonitor.h.

NRPChits

MonitorElement* PhysicsObjectsMonitor::NRPChits

private

Definition at line 62 of file PhysicsObjectsMonitor.h.

pt

MonitorElement* PhysicsObjectsMonitor::pt

private

Definition at line 51 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().

ptot

MonitorElement* PhysicsObjectsMonitor::ptot

private

Definition at line 50 of file PhysicsObjectsMonitor.h.

px

MonitorElement* PhysicsObjectsMonitor::px

private

Definition at line 52 of file PhysicsObjectsMonitor.h.

py

MonitorElement* PhysicsObjectsMonitor::py

private

Definition at line 53 of file PhysicsObjectsMonitor.h.

pz

MonitorElement* PhysicsObjectsMonitor::pz

private

Definition at line 54 of file PhysicsObjectsMonitor.h.

theDataType

std::string PhysicsObjectsMonitor::theDataType

private

Definition at line 64 of file PhysicsObjectsMonitor.h.

theSeedCollectionLabel

std::string PhysicsObjectsMonitor::theSeedCollectionLabel

private

Definition at line 42 of file PhysicsObjectsMonitor.h.

theSTAMuonLabel

std::string PhysicsObjectsMonitor::theSTAMuonLabel

private

Definition at line 41 of file PhysicsObjectsMonitor.h.

theSTAMuonToken_

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

private

Definition at line 67 of file PhysicsObjectsMonitor.h.