TrackerHitAnalyzer Class Reference

#include <TrackerHitAnalyzer.h>

Inheritance diagram for TrackerHitAnalyzer:

Public Member Functions
	TrackerHitAnalyzer (const edm::ParameterSet &ps)
	Constructor. More...
	~TrackerHitAnalyzer () 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
	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

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override
	Analyze. More...
void	bookHistograms (DQMStore::IBooker &ibooker, const edm::Run &run, const edm::EventSetup &es) override
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Attributes
edm::ParameterSet	conf_
edm::EDGetTokenT< edm::PSimHitContainer >	edmPSimHitContainer_pxlBrlHigh_Token_
edm::EDGetTokenT< edm::PSimHitContainer >	edmPSimHitContainer_pxlBrlLow_Token_
edm::EDGetTokenT< edm::PSimHitContainer >	edmPSimHitContainer_pxlFwdHigh_Token_
edm::EDGetTokenT< edm::PSimHitContainer >	edmPSimHitContainer_pxlFwdLow_Token_
edm::EDGetTokenT< edm::PSimHitContainer >	edmPSimHitContainer_siTECHigh_Token_
edm::EDGetTokenT< edm::PSimHitContainer >	edmPSimHitContainer_siTECLow_Token_
edm::EDGetTokenT< edm::PSimHitContainer >	edmPSimHitContainer_siTIBHigh_Token_
edm::EDGetTokenT< edm::PSimHitContainer >	edmPSimHitContainer_siTIBLow_Token_
edm::EDGetTokenT< edm::PSimHitContainer >	edmPSimHitContainer_siTIDHigh_Token_
edm::EDGetTokenT< edm::PSimHitContainer >	edmPSimHitContainer_siTIDLow_Token_
edm::EDGetTokenT< edm::PSimHitContainer >	edmPSimHitContainer_siTOBHigh_Token_
edm::EDGetTokenT< edm::PSimHitContainer >	edmPSimHitContainer_siTOBLow_Token_
edm::EDGetTokenT< edm::SimTrackContainer >	edmSimTrackContainerToken_
DQMStore *	fDBE
std::string	fOutputFile
MonitorElement *	h1e [12]
MonitorElement *	h1ex [12]
MonitorElement *	h1ey [12]
MonitorElement *	h1ez [12]
MonitorElement *	h1lx [12]
MonitorElement *	h1ly [12]
MonitorElement *	h2e [12]
MonitorElement *	h2ex [12]
MonitorElement *	h2ey [12]
MonitorElement *	h2ez [12]
MonitorElement *	h2lx [12]
MonitorElement *	h2ly [12]
MonitorElement *	h3e [12]
MonitorElement *	h3ex [12]
MonitorElement *	h3ey [12]
MonitorElement *	h3ez [12]
MonitorElement *	h3lx [12]
MonitorElement *	h3ly [12]
MonitorElement *	h4e [12]
MonitorElement *	h4ex [12]
MonitorElement *	h4ey [12]
MonitorElement *	h4ez [12]
MonitorElement *	h4lx [12]
MonitorElement *	h4ly [12]
MonitorElement *	h5e [12]
MonitorElement *	h5ex [12]
MonitorElement *	h5ey [12]
MonitorElement *	h5ez [12]
MonitorElement *	h5lx [12]
MonitorElement *	h5ly [12]
MonitorElement *	h6e [12]
MonitorElement *	h6ex [12]
MonitorElement *	h6ey [12]
MonitorElement *	h6ez [12]
MonitorElement *	h6lx [12]
MonitorElement *	h6ly [12]
MonitorElement *	htofeta
MonitorElement *	htofeta_profile
MonitorElement *	htofphi
MonitorElement *	htofphi_profile
MonitorElement *	htofr
MonitorElement *	htofr_profile
MonitorElement *	htofz
MonitorElement *	htofz_profile
bool	pixelOutput
bool	runStandalone
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tGeomEsToken_
bool	verbose_

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 >
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
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 Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 24 of file TrackerHitAnalyzer.h.

Constructor & Destructor Documentation

TrackerHitAnalyzer()

TrackerHitAnalyzer::TrackerHitAnalyzer

(

const edm::ParameterSet &

ps

)

Constructor.

Definition at line 38 of file TrackerHitAnalyzer.cc.

    : tGeomEsToken_(esConsumes()),
      verbose_(ps.getUntrackedParameter<bool>("Verbosity", false)),
      edmPSimHitContainer_pxlBrlLow_Token_(
          consumes<edm::PSimHitContainer>(ps.getParameter<edm::InputTag>("PxlBrlLowSrc"))),
      edmPSimHitContainer_pxlBrlHigh_Token_(
          consumes<edm::PSimHitContainer>(ps.getParameter<edm::InputTag>("PxlBrlHighSrc"))),
      edmPSimHitContainer_pxlFwdLow_Token_(
          consumes<edm::PSimHitContainer>(ps.getParameter<edm::InputTag>("PxlFwdLowSrc"))),
      edmPSimHitContainer_pxlFwdHigh_Token_(
          consumes<edm::PSimHitContainer>(ps.getParameter<edm::InputTag>("PxlFwdHighSrc"))),
      edmPSimHitContainer_siTIBLow_Token_(
          consumes<edm::PSimHitContainer>(ps.getParameter<edm::InputTag>("SiTIBLowSrc"))),
      edmPSimHitContainer_siTIBHigh_Token_(
          consumes<edm::PSimHitContainer>(ps.getParameter<edm::InputTag>("SiTIBHighSrc"))),
      edmPSimHitContainer_siTOBLow_Token_(
          consumes<edm::PSimHitContainer>(ps.getParameter<edm::InputTag>("SiTOBLowSrc"))),
      edmPSimHitContainer_siTOBHigh_Token_(
          consumes<edm::PSimHitContainer>(ps.getParameter<edm::InputTag>("SiTOBHighSrc"))),
      edmPSimHitContainer_siTIDLow_Token_(
          consumes<edm::PSimHitContainer>(ps.getParameter<edm::InputTag>("SiTIDLowSrc"))),
      edmPSimHitContainer_siTIDHigh_Token_(
          consumes<edm::PSimHitContainer>(ps.getParameter<edm::InputTag>("SiTIDHighSrc"))),
      edmPSimHitContainer_siTECLow_Token_(
          consumes<edm::PSimHitContainer>(ps.getParameter<edm::InputTag>("SiTECLowSrc"))),
      edmPSimHitContainer_siTECHigh_Token_(
          consumes<edm::PSimHitContainer>(ps.getParameter<edm::InputTag>("SiTECHighSrc"))),
      edmSimTrackContainerToken_(consumes<edm::SimTrackContainer>(ps.getParameter<edm::InputTag>("G4TrkSrc"))),
      fDBE(nullptr),
      conf_(ps),
      runStandalone(ps.getParameter<bool>("runStandalone")),
      fOutputFile(ps.getUntrackedParameter<std::string>("outputFile", "TrackerHitHisto.root")),
      pixelOutput(ps.getParameter<bool>("pixelOutput")) {}

~TrackerHitAnalyzer()

TrackerHitAnalyzer::~TrackerHitAnalyzer ( )

override

Destructor.

Definition at line 330 of file TrackerHitAnalyzer.cc.

                                        {
  // don't try to delete any pointers - they're handled by DQM machinery
}

Member Function Documentation

analyze()

void TrackerHitAnalyzer::analyze ( const edm::Event &  e, const edm::EventSetup &  c  )

overrideprotectedvirtual

Analyze.

Reimplemented from DQMEDAnalyzer.

Definition at line 334 of file TrackerHitAnalyzer.cc.

References DummyCfis::c, Geom::Phi< T1, Range >::degrees(), MillePedeFileConverter_cfg::e, edmPSimHitContainer_pxlBrlHigh_Token_, edmPSimHitContainer_pxlBrlLow_Token_, edmPSimHitContainer_pxlFwdHigh_Token_, edmPSimHitContainer_pxlFwdLow_Token_, edmPSimHitContainer_siTECHigh_Token_, edmPSimHitContainer_siTECLow_Token_, edmPSimHitContainer_siTIBHigh_Token_, edmPSimHitContainer_siTIBLow_Token_, edmPSimHitContainer_siTIDHigh_Token_, edmPSimHitContainer_siTIDLow_Token_, edmPSimHitContainer_siTOBHigh_Token_, edmPSimHitContainer_siTOBLow_Token_, edmSimTrackContainerToken_, PVValHelper::eta, PV3DBase< T, PVType, FrameType >::eta(), dqm::impl::MonitorElement::Fill(), h1e, h1ex, h1ey, h1ez, h1lx, h1ly, h2e, h2ex, h2ey, h2ez, h2lx, h2ly, h3e, h3ex, h3ey, h3ez, h3lx, h3ly, h4e, h4ex, h4ey, h4ez, h4lx, h4ly, h5e, h5ex, h5ey, h5ez, h5lx, h5ly, h6e, h6ex, h6ey, h6ez, h6lx, h6ly, htofeta, htofphi, htofr, htofz, edm::HandleBase::isValid(), dqm-mbProfile::log, PV3DBase< T, PVType, FrameType >::mag(), METSkim_cff::Max, METSkim_cff::Min, AlCaHLTBitMon_ParallelJobs::p, PV3DBase< T, PVType, FrameType >::phi(), pixelOutput, mathSSE::sqrt(), funct::tan(), tGeomEsToken_, theta(), PbPb_ZMuSkimMuonDPG_cff::tracker, and PV3DBase< T, PVType, FrameType >::z().

                                                                          {
  edm::LogInfo("EventInfo") << " Run = " << e.id().run() << " Event = " << e.id().event();

  // iterator to access containers
  edm::PSimHitContainer::const_iterator itHit;
  // get Pixel Barrel information
  // extract low container
  edm::Handle<edm::PSimHitContainer> PxlBrlLowContainer;
  e.getByToken(edmPSimHitContainer_pxlBrlLow_Token_, PxlBrlLowContainer);
  if (!PxlBrlLowContainer.isValid()) {
    edm::LogError("TrackerHitAnalyzer::analyze") << "Unable to find TrackerHitsPixelBarrelLowTof in event!";
    return;
  }
  // extract high container
  edm::Handle<edm::PSimHitContainer> PxlBrlHighContainer;
  e.getByToken(edmPSimHitContainer_pxlBrlHigh_Token_, PxlBrlHighContainer);
  if (!PxlBrlHighContainer.isValid()) {
    edm::LogError("TrackerHitAnalyzer::analyze") << "Unable to find TrackerHitsPixelBarrelHighTof in event!";
    return;
  }
  // get Pixel Forward information
  // extract low container
  edm::Handle<edm::PSimHitContainer> PxlFwdLowContainer;
  e.getByToken(edmPSimHitContainer_pxlFwdLow_Token_, PxlFwdLowContainer);
  if (!PxlFwdLowContainer.isValid()) {
    edm::LogError("TrackerHitAnalyzer::analyze") << "Unable to find TrackerHitsPixelEndcapLowTof in event!";
    return;
  }
  // extract high container
  edm::Handle<edm::PSimHitContainer> PxlFwdHighContainer;
  e.getByToken(edmPSimHitContainer_pxlFwdHigh_Token_, PxlFwdHighContainer);
  if (!PxlFwdHighContainer.isValid()) {
    edm::LogError("TrackerHitAnalyzer::analyze") << "Unable to find TrackerHitsPixelEndcapHighTof in event!";
    return;
  }
  // get Silicon TIB information
  // extract TIB low container
  edm::Handle<edm::PSimHitContainer> SiTIBLowContainer;
  e.getByToken(edmPSimHitContainer_siTIBLow_Token_, SiTIBLowContainer);
  if (!SiTIBLowContainer.isValid()) {
    edm::LogError("TrackerHitProducer::analyze") << "Unable to find TrackerHitsTIBLowTof in event!";
    return;
  }
  // extract TIB high container
  edm::Handle<edm::PSimHitContainer> SiTIBHighContainer;
  e.getByToken(edmPSimHitContainer_siTIBHigh_Token_, SiTIBHighContainer);
  if (!SiTIBHighContainer.isValid()) {
    edm::LogError("TrackerHitProducer::analyze") << "Unable to find TrackerHitsTIBHighTof in event!";
    return;
  }
  // get Silicon TOB information
  // extract TOB low container
  edm::Handle<edm::PSimHitContainer> SiTOBLowContainer;
  e.getByToken(edmPSimHitContainer_siTOBLow_Token_, SiTOBLowContainer);
  if (!SiTOBLowContainer.isValid()) {
    edm::LogError("TrackerHitProducer::analyze") << "Unable to find TrackerHitsTOBLowTof in event!";
    return;
  }
  // extract TOB high container
  edm::Handle<edm::PSimHitContainer> SiTOBHighContainer;
  e.getByToken(edmPSimHitContainer_siTOBHigh_Token_, SiTOBHighContainer);
  if (!SiTOBHighContainer.isValid()) {
    edm::LogError("TrackerHitProducer::analyze") << "Unable to find TrackerHitsTOBHighTof in event!";
    return;
  }

  // get Silicon TID information
  // extract TID low container
  edm::Handle<edm::PSimHitContainer> SiTIDLowContainer;
  e.getByToken(edmPSimHitContainer_siTIDLow_Token_, SiTIDLowContainer);
  if (!SiTIDLowContainer.isValid()) {
    edm::LogError("TrackerHitProducer::analyze") << "Unable to find TrackerHitsTIDLowTof in event!";
    return;
  }
  // extract TID high container
  edm::Handle<edm::PSimHitContainer> SiTIDHighContainer;
  e.getByToken(edmPSimHitContainer_siTIDHigh_Token_, SiTIDHighContainer);
  if (!SiTIDHighContainer.isValid()) {
    edm::LogError("TrackerHitProducer::analyze") << "Unable to find TrackerHitsTIDHighTof in event!";
    return;
  }
  // get Silicon TEC information
  // extract TEC low container
  edm::Handle<edm::PSimHitContainer> SiTECLowContainer;
  e.getByToken(edmPSimHitContainer_siTECLow_Token_, SiTECLowContainer);
  if (!SiTECLowContainer.isValid()) {
    edm::LogError("TrackerHitProducer::analyze") << "Unable to find TrackerHitsTECLowTof in event!";
    return;
  }
  // extract TEC high container
  edm::Handle<edm::PSimHitContainer> SiTECHighContainer;
  e.getByToken(edmPSimHitContainer_siTECHigh_Token_, SiTECHighContainer);
  if (!SiTECHighContainer.isValid()) {
    edm::LogError("TrackerHitProducer::analyze") << "Unable to find TrackerHitsTECHighTof in event!";
    return;
  }

  // get G4Track information

  edm::Handle<edm::SimTrackContainer> G4TrkContainer;
  e.getByToken(edmSimTrackContainerToken_, G4TrkContainer);
  if (!G4TrkContainer.isValid()) {
    edm::LogError("TrackerHitAnalyzer::analyze") << "Unable to find SimTrack in event!";
    return;
  }

  // Get geometry information
  const auto &tracker = &c.getData(tGeomEsToken_);

  int ir = -100;
  edm::SimTrackContainer::const_iterator itTrk;
  for (itTrk = G4TrkContainer->begin(); itTrk != G4TrkContainer->end(); ++itTrk) {
    //    std::cout << "itTrk = "<< itTrk << std::endl;
    double eta = 0, p = 0;
    const CLHEP::HepLorentzVector &G4Trk = CLHEP::HepLorentzVector(
        itTrk->momentum().x(), itTrk->momentum().y(), itTrk->momentum().z(), itTrk->momentum().e());
    p = sqrt(G4Trk[0] * G4Trk[0] + G4Trk[1] * G4Trk[1] + G4Trk[2] * G4Trk[2]);
    if (p == 0)
      edm::LogError("TrackerHitAnalyzer::analyze") << "TrackerTest::INFO: Primary has p = 0 ";
    else {
      double costheta = G4Trk[2] / p;
      double theta = acos(TMath::Min(TMath::Max(costheta, -1.), 1.));
      eta = -log(tan(theta / 2));

      if (eta > 0.0 && eta <= 0.5)
        ir = 0;
      if (eta > 0.5 && eta <= 1.0)
        ir = 1;
      if (eta > 1.0 && eta <= 1.5)
        ir = 2;
      if (eta > 1.5 && eta <= 2.0)
        ir = 3;
      if (eta > 2.0 && eta <= 2.5)
        ir = 4;
      if (eta > 2.5)
        ir = 5;

      if (eta > -0.5 && eta <= 0.0)
        ir = 6;
      if (eta > -1.0 && eta <= -0.5)
        ir = 7;
      if (eta > -1.5 && eta <= -1.0)
        ir = 8;
      if (eta > -2.0 && eta <= -1.5)
        ir = 9;
      if (eta > -2.5 && eta <= -2.0)
        ir = 10;
      if (eta <= -2.5)
        ir = 11;
      //          edm::LogInfo("EventInfo") << " eta = " << eta << " ir = " <<
      //          ir;
      //      std::cout << " " << std::endl;
      //          std::cout << "eta " << eta << " ir = " << ir << std::endl;
      //      std::cout << " " << std::endl;
    }
  }
  // get Pixel information
  // If Phase 1, do not run - will run in new Phase 1 module

  if (pixelOutput) {
    for (itHit = PxlBrlLowContainer->begin(); itHit != PxlBrlLowContainer->end(); ++itHit) {
      DetId detid = DetId(itHit->detUnitId());
      const GeomDetUnit *det = (const GeomDetUnit *)tracker->idToDetUnit(detid);
      GlobalPoint gpos = det->toGlobal(itHit->localPosition());
      htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
      htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
      htofr->Fill(gpos.mag(), itHit->timeOfFlight());
      htofz->Fill(gpos.z(), itHit->timeOfFlight());

      h5e[ir]->Fill(itHit->energyLoss());
      h5ex[ir]->Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
      h5ey[ir]->Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
      h5ez[ir]->Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
      h5lx[ir]->Fill(itHit->localPosition().x());
      h5ly[ir]->Fill(itHit->localPosition().y());
    }
    for (itHit = PxlBrlHighContainer->begin(); itHit != PxlBrlHighContainer->end(); ++itHit) {
      DetId detid = DetId(itHit->detUnitId());
      const GeomDetUnit *det = (const GeomDetUnit *)tracker->idToDetUnit(detid);
      GlobalPoint gpos = det->toGlobal(itHit->localPosition());
      htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
      htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
      htofr->Fill(gpos.mag(), itHit->timeOfFlight());
      htofz->Fill(gpos.z(), itHit->timeOfFlight());

      h5e[ir]->Fill(itHit->energyLoss());
      h5ex[ir]->Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
      h5ey[ir]->Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
      h5ez[ir]->Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
      h5lx[ir]->Fill(itHit->localPosition().x());
      h5ly[ir]->Fill(itHit->localPosition().y());
    }
    for (itHit = PxlFwdLowContainer->begin(); itHit != PxlFwdLowContainer->end(); ++itHit) {
      DetId detid = DetId(itHit->detUnitId());
      const GeomDetUnit *det = (const GeomDetUnit *)tracker->idToDetUnit(detid);
      GlobalPoint gpos = det->toGlobal(itHit->localPosition());
      htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
      htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
      htofr->Fill(gpos.mag(), itHit->timeOfFlight());
      htofz->Fill(gpos.z(), itHit->timeOfFlight());

      h6e[ir]->Fill(itHit->energyLoss());
      h6ex[ir]->Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
      h6ey[ir]->Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
      h6ez[ir]->Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
      h6lx[ir]->Fill(itHit->localPosition().x());
      h6ly[ir]->Fill(itHit->localPosition().y());
    }
    for (itHit = PxlFwdHighContainer->begin(); itHit != PxlFwdHighContainer->end(); ++itHit) {
      DetId detid = DetId(itHit->detUnitId());
      const GeomDetUnit *det = (const GeomDetUnit *)tracker->idToDetUnit(detid);
      GlobalPoint gpos = det->toGlobal(itHit->localPosition());
      htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
      htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
      htofr->Fill(gpos.mag(), itHit->timeOfFlight());
      htofz->Fill(gpos.z(), itHit->timeOfFlight());

      h6e[ir]->Fill(itHit->energyLoss());
      h6ex[ir]->Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
      h6ey[ir]->Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
      h6ez[ir]->Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
      h6lx[ir]->Fill(itHit->localPosition().x());
      h6ly[ir]->Fill(itHit->localPosition().y());
    }
  }
  // get TIB information
  for (itHit = SiTIBLowContainer->begin(); itHit != SiTIBLowContainer->end(); ++itHit) {
    DetId detid = DetId(itHit->detUnitId());
    const GeomDetUnit *det = (const GeomDetUnit *)tracker->idToDetUnit(detid);
    GlobalPoint gpos = det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());

    h1e[ir]->Fill(itHit->energyLoss());
    h1ex[ir]->Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
    h1ey[ir]->Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
    h1ez[ir]->Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
    h1lx[ir]->Fill(itHit->localPosition().x());
    h1ly[ir]->Fill(itHit->localPosition().y());
  }
  for (itHit = SiTIBHighContainer->begin(); itHit != SiTIBHighContainer->end(); ++itHit) {
    DetId detid = DetId(itHit->detUnitId());
    const GeomDetUnit *det = (const GeomDetUnit *)tracker->idToDetUnit(detid);
    GlobalPoint gpos = det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());

    h1e[ir]->Fill(itHit->energyLoss());
    h1ex[ir]->Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
    h1ey[ir]->Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
    h1ez[ir]->Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
    h1lx[ir]->Fill(itHit->localPosition().x());
    h1ly[ir]->Fill(itHit->localPosition().y());
  }
  // get TOB information
  for (itHit = SiTOBLowContainer->begin(); itHit != SiTOBLowContainer->end(); ++itHit) {
    DetId detid = DetId(itHit->detUnitId());
    const GeomDetUnit *det = (const GeomDetUnit *)tracker->idToDetUnit(detid);
    GlobalPoint gpos = det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());

    h2e[ir]->Fill(itHit->energyLoss());
    h2ex[ir]->Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
    h2ey[ir]->Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
    h2ez[ir]->Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
    h2lx[ir]->Fill(itHit->localPosition().x());
    h2ly[ir]->Fill(itHit->localPosition().y());
  }
  for (itHit = SiTOBHighContainer->begin(); itHit != SiTOBHighContainer->end(); ++itHit) {
    DetId detid = DetId(itHit->detUnitId());
    const GeomDetUnit *det = (const GeomDetUnit *)tracker->idToDetUnit(detid);
    GlobalPoint gpos = det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());

    h2e[ir]->Fill(itHit->energyLoss());
    h2ex[ir]->Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
    h2ey[ir]->Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
    h2ez[ir]->Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
    h2lx[ir]->Fill(itHit->localPosition().x());
    h2ly[ir]->Fill(itHit->localPosition().y());
  }
  // get TID information
  for (itHit = SiTIDLowContainer->begin(); itHit != SiTIDLowContainer->end(); ++itHit) {
    DetId detid = DetId(itHit->detUnitId());
    const GeomDetUnit *det = (const GeomDetUnit *)tracker->idToDetUnit(detid);
    GlobalPoint gpos = det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());

    h3e[ir]->Fill(itHit->energyLoss());
    h3ex[ir]->Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
    h3ey[ir]->Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
    h3ez[ir]->Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
    h3lx[ir]->Fill(itHit->localPosition().x());
    h3ly[ir]->Fill(itHit->localPosition().y());
  }
  for (itHit = SiTIDHighContainer->begin(); itHit != SiTIDHighContainer->end(); ++itHit) {
    DetId detid = DetId(itHit->detUnitId());
    const GeomDetUnit *det = (const GeomDetUnit *)tracker->idToDetUnit(detid);
    GlobalPoint gpos = det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());

    h3e[ir]->Fill(itHit->energyLoss());
    h3ex[ir]->Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
    h3ey[ir]->Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
    h3ez[ir]->Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
    h3lx[ir]->Fill(itHit->localPosition().x());
    h3ly[ir]->Fill(itHit->localPosition().y());
  }
  // get TEC information
  for (itHit = SiTECLowContainer->begin(); itHit != SiTECLowContainer->end(); ++itHit) {
    DetId detid = DetId(itHit->detUnitId());
    const GeomDetUnit *det = (const GeomDetUnit *)tracker->idToDetUnit(detid);
    GlobalPoint gpos = det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());

    h4e[ir]->Fill(itHit->energyLoss());
    h4ex[ir]->Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
    h4ey[ir]->Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
    h4ez[ir]->Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
    h4lx[ir]->Fill(itHit->localPosition().x());
    h4ly[ir]->Fill(itHit->localPosition().y());
  }
  for (itHit = SiTECHighContainer->begin(); itHit != SiTECHighContainer->end(); ++itHit) {
    DetId detid = DetId(itHit->detUnitId());
    const GeomDetUnit *det = (const GeomDetUnit *)tracker->idToDetUnit(detid);
    GlobalPoint gpos = det->toGlobal(itHit->localPosition());
    htofeta->Fill(gpos.eta(), itHit->timeOfFlight());
    htofphi->Fill(gpos.phi().degrees(), itHit->timeOfFlight());
    htofr->Fill(gpos.mag(), itHit->timeOfFlight());
    htofz->Fill(gpos.z(), itHit->timeOfFlight());

    h4e[ir]->Fill(itHit->energyLoss());
    h4ex[ir]->Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
    h4ey[ir]->Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
    h4ez[ir]->Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
    h4lx[ir]->Fill(itHit->localPosition().x());
    h4ly[ir]->Fill(itHit->localPosition().y());
  }

  return;
}

bookHistograms()

void TrackerHitAnalyzer::bookHistograms ( DQMStore::IBooker &  ibooker, const edm::Run &  run, const edm::EventSetup &  es  )

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 72 of file TrackerHitAnalyzer.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), fDBE, h1e, h1ex, h1ey, h1ez, h1lx, h1ly, h2e, h2ex, h2ey, h2ez, h2lx, h2ly, h3e, h3ex, h3ey, h3ez, h3lx, h3ly, h4e, h4ex, h4ey, h4ez, h4lx, h4ly, h5e, h5ex, h5ey, h5ez, h5lx, h5ly, h6e, h6ex, h6ey, h6ez, h6lx, h6ly, LaserClient_cfi::high, htofeta, htofphi, htofr, htofz, mps_fire::i, LaserClient_cfi::low, Utilities::operator, pixelOutput, g4SimHits_cfi::Region, and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                                                                            {
  fDBE = edm::Service<DQMStore>().operator->();

  Char_t hname1[50], htitle1[80];
  Char_t hname2[50], htitle2[80];
  Char_t hname3[50], htitle3[80];
  Char_t hname4[50], htitle4[80];
  Char_t hname5[50], htitle5[80];
  Char_t hname6[50], htitle6[80];

  if (fDBE != nullptr) {
    //   fDBE->setCurrentFolder("TrackerHitsV/TrackerHitTask");

    // is there any way to record CPU Info ???
    // if so, it can be done once - via beginJob()
    int nbin = 5000;

    ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/");
    htofeta = ibooker.book2D("tof_eta", "Time of flight vs eta", nbin, -3.0, 3.0, 200, -100, 100);
    htofphi = ibooker.book2D("tof_phi", "Time of flight vs phi", nbin, -180, 180, 200, -100, 100);
    htofr = ibooker.book2D("tof_r", "Time of flight vs r", nbin, 0, 300, 200, -100, 100);
    htofz = ibooker.book2D("tof_z", "Time of flight vs z", nbin, -280, 280, 200, -100, 100);

    const float E2NEL = 1.;

    const char *Region[] = {"005", "051", "115", "152", "225", "253", "-050", "-105", "-151", "-215", "-252", "-325"};
    nbin = 200;

    // Energy loss histograms
    for (int i = 0; i < 12; i++) {
      sprintf(htitle1, "Energy loss in TIB %s", Region[i]);
      sprintf(htitle2, "Energy loss in TOB %s", Region[i]);
      sprintf(htitle3, "Energy loss in TID %s", Region[i]);
      sprintf(htitle4, "Energy loss in TEC %s", Region[i]);
      if (pixelOutput)
        sprintf(htitle5, "Energy loss in BPIX %s", Region[i]);
      if (pixelOutput)
        sprintf(htitle6, "Energy loss in FPIX %s", Region[i]);

      sprintf(hname1, "Eloss_TIB_%i", i + 1);
      sprintf(hname2, "Eloss_TOB_%i", i + 1);
      sprintf(hname3, "Eloss_TID_%i", i + 1);
      sprintf(hname4, "Eloss_TEC_%i", i + 1);
      if (pixelOutput)
        sprintf(hname5, "Eloss_BPIX_%i", i + 1);
      if (pixelOutput)
        sprintf(hname6, "Eloss_FPIX_%i", i + 1);

      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TIBHit");
      h1e[i] = ibooker.book1D(hname1, htitle1, nbin, 0.0, 0.001 * E2NEL);
      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TOBHit");
      h2e[i] = ibooker.book1D(hname2, htitle2, nbin, 0.0, 0.001 * E2NEL);
      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TIDHit");
      h3e[i] = ibooker.book1D(hname3, htitle3, nbin, 0.0, 0.001 * E2NEL);
      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TECHit");
      h4e[i] = ibooker.book1D(hname4, htitle4, nbin, 0.0, 0.001 * E2NEL);
      if (pixelOutput) {
        ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/BPIXHit");
        h5e[i] = ibooker.book1D(hname5, htitle5, nbin, 0.0, 0.001 * E2NEL);
        ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/FPIXHit");
        h6e[i] = ibooker.book1D(hname6, htitle6, nbin, 0.0, 0.001 * E2NEL);
      }
    }

    // limits
    const float high[] = {0.03, 0.03, 0.02, 0.03, 0.03, 0.03};
    const float low[] = {-0.03, -0.03, -0.02, -0.03, -0.03, -0.03};

    for (int i = 0; i < 12; i++) {
      sprintf(htitle1, "Entryx-Exitx in TIB %s", Region[i]);
      sprintf(htitle2, "Entryx-Exitx in TOB %s", Region[i]);
      sprintf(htitle3, "Entryx-Exitx in TID %s", Region[i]);
      sprintf(htitle4, "Entryx-Exitx in TEC %s", Region[i]);
      if (pixelOutput)
        sprintf(htitle5, "Entryx-Exitx in BPIX %s", Region[i]);
      if (pixelOutput)
        sprintf(htitle6, "Entryx-Exitx in FPIX %s", Region[i]);

      sprintf(hname1, "Entryx-Exitx_TIB_%i", i + 1);
      sprintf(hname2, "Entryx-Exitx_TOB_%i", i + 1);
      sprintf(hname3, "Entryx-Exitx_TID_%i", i + 1);
      sprintf(hname4, "Entryx-Exitx_TEC_%i", i + 1);
      if (pixelOutput)
        sprintf(hname5, "Entryx-Exitx_BPIX_%i", i + 1);
      if (pixelOutput)
        sprintf(hname6, "Entryx-Exitx_FPIX_%i", i + 1);

      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TIBHit");
      h1ex[i] = ibooker.book1D(hname1, htitle1, nbin, low[0], high[0]);
      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TOBHit");
      h2ex[i] = ibooker.book1D(hname2, htitle2, nbin, low[1], high[1]);
      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TIDHit");
      h3ex[i] = ibooker.book1D(hname3, htitle3, nbin, low[2], high[2]);
      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TECHit");
      h4ex[i] = ibooker.book1D(hname4, htitle4, nbin, low[3], high[3]);
      if (pixelOutput) {
        ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/BPIXHit");
        h5ex[i] = ibooker.book1D(hname5, htitle5, nbin, low[4], high[4]);
        ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/FPIXHit");
        h6ex[i] = ibooker.book1D(hname6, htitle6, nbin, low[5], high[5]);
      }
    }

    const float high0[] = {0.05, 0.06, 0.03, 0.03, 0.03, 0.03};
    const float low0[] = {-0.05, -0.06, -0.03, -0.03, -0.03, -0.03};

    for (int i = 0; i < 12; i++) {
      sprintf(htitle1, "Entryy-Exity in TIB %s", Region[i]);
      sprintf(htitle2, "Entryy-Exity in TOB %s", Region[i]);
      sprintf(htitle3, "Entryy-Exity in TID %s", Region[i]);
      sprintf(htitle4, "Entryy-Exity in TEC %s", Region[i]);
      if (pixelOutput)
        sprintf(htitle5, "Entryy-Exity in BPIX %s", Region[i]);
      if (pixelOutput)
        sprintf(htitle6, "Entryy-Exity in FPIX %s", Region[i]);

      sprintf(hname1, "Entryy-Exity_TIB_%i", i + 1);
      sprintf(hname2, "Entryy-Exity_TOB_%i", i + 1);
      sprintf(hname3, "Entryy-Exity_TID_%i", i + 1);
      sprintf(hname4, "Entryy-Exity_TEC_%i", i + 1);
      if (pixelOutput)
        sprintf(hname5, "Entryy-Exity_BPIX_%i", i + 1);
      if (pixelOutput)
        sprintf(hname6, "Entryy-Exity_FPIX_%i", i + 1);

      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TIBHit");
      h1ey[i] = ibooker.book1D(hname1, htitle1, nbin, low0[0], high0[0]);
      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TOBHit");
      h2ey[i] = ibooker.book1D(hname2, htitle2, nbin, low0[1], high0[1]);
      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TIDHit");
      h3ey[i] = ibooker.book1D(hname3, htitle3, nbin, low0[2], high0[2]);
      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TECHit");
      h4ey[i] = ibooker.book1D(hname4, htitle4, nbin, low0[3], high0[3]);
      if (pixelOutput) {
        ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/BPIXHit");
        h5ey[i] = ibooker.book1D(hname5, htitle5, nbin, low0[4], high0[4]);
        ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/FPIXHit");
        h6ey[i] = ibooker.book1D(hname6, htitle6, nbin, low0[5], high0[5]);
      }
    }

    const float high1[] = {0.05, 0.06, 0.05, 0.06, 0.05, 0.05};
    const float low1[] = {0., 0., 0., 0., 0., 0.};

    for (int i = 0; i < 12; i++) {
      sprintf(htitle1, "abs(Entryz-Exitz) in TIB %s", Region[i]);
      sprintf(htitle2, "abs(Entryz-Exitz) in TOB %s", Region[i]);
      sprintf(htitle3, "abs(Entryz-Exitz) in TID %s", Region[i]);
      sprintf(htitle4, "abs(Entryz-Exitz) in TEC %s", Region[i]);
      if (pixelOutput)
        sprintf(htitle5, "abs(Entryz-Exitz) in BPIX %s", Region[i]);
      if (pixelOutput)
        sprintf(htitle6, "abs(Entryz-Exitz) in FPIX %s", Region[i]);

      sprintf(hname1, "Entryz-Exitz_TIB_%i", i + 1);
      sprintf(hname2, "Entryz-Exitz_TOB_%i", i + 1);
      sprintf(hname3, "Entryz-Exitz_TID_%i", i + 1);
      sprintf(hname4, "Entryz-Exitz_TEC_%i", i + 1);
      if (pixelOutput)
        sprintf(hname5, "Entryz-Exitz_BPIX_%i", i + 1);
      if (pixelOutput)
        sprintf(hname6, "Entryz-Exitz_FPIX_%i", i + 1);

      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TIBHit");
      h1ez[i] = ibooker.book1D(hname1, htitle1, nbin, low1[0], high1[0]);
      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TOBHit");
      h2ez[i] = ibooker.book1D(hname2, htitle2, nbin, low1[1], high1[1]);
      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TIDHit");
      h3ez[i] = ibooker.book1D(hname3, htitle3, nbin, low1[2], high1[2]);
      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TECHit");
      h4ez[i] = ibooker.book1D(hname4, htitle4, nbin, low1[3], high1[3]);
      if (pixelOutput) {
        ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/BPIXHit");
        h5ez[i] = ibooker.book1D(hname5, htitle5, nbin, low1[4], high1[4]);
        ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/FPIXHit");
        h6ez[i] = ibooker.book1D(hname6, htitle6, nbin, low1[5], high1[5]);
      }
    }

    const float high2[] = {3.2, 5.0, 5.5, 6.2, 0.85, 0.5};
    const float low2[] = {-3.2, -5.0, -5.5, -6.2, -0.85, -0.5};

    for (int i = 0; i < 12; i++) {
      sprintf(htitle1, "Localx in TIB %s", Region[i]);
      sprintf(htitle2, "Localx in TOB %s", Region[i]);
      sprintf(htitle3, "Localx in TID %s", Region[i]);
      sprintf(htitle4, "Localx in TEC %s", Region[i]);
      if (pixelOutput)
        sprintf(htitle5, "Localx in BPIX %s", Region[i]);
      if (pixelOutput)
        sprintf(htitle6, "Localx in FPIX %s", Region[i]);

      sprintf(hname1, "Localx_TIB_%i", i + 1);
      sprintf(hname2, "Localx_TOB_%i", i + 1);
      sprintf(hname3, "Localx_TID_%i", i + 1);
      sprintf(hname4, "Localx_TEC_%i", i + 1);
      if (pixelOutput)
        sprintf(hname5, "Localx_BPIX_%i", i + 1);
      if (pixelOutput)
        sprintf(hname6, "Localx_FPIX_%i", i + 1);

      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TIBHit");
      h1lx[i] = ibooker.book1D(hname1, htitle1, nbin, low2[0], high2[0]);
      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TOBHit");
      h2lx[i] = ibooker.book1D(hname2, htitle2, nbin, low2[1], high2[1]);
      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TIDHit");
      h3lx[i] = ibooker.book1D(hname3, htitle3, nbin, low2[2], high2[2]);
      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TECHit");
      h4lx[i] = ibooker.book1D(hname4, htitle4, nbin, low2[3], high2[3]);
      if (pixelOutput) {
        ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/BPIXHit");
        h5lx[i] = ibooker.book1D(hname5, htitle5, nbin, low2[4], high2[4]);
        ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/FPIXHit");
        h6lx[i] = ibooker.book1D(hname6, htitle6, nbin, low2[5], high2[5]);
      }
    }

    const float high3[] = {6.0, 10., 5.6, 10.5, 3.4, 0.52};
    const float low3[] = {-6.0, -10., -5.6, -10.5, -3.4, -0.52};

    for (int i = 0; i < 12; i++) {
      sprintf(htitle1, "Localy in TIB %s", Region[i]);
      sprintf(htitle2, "Localy in TOB %s", Region[i]);
      sprintf(htitle3, "Localy in TID %s", Region[i]);
      sprintf(htitle4, "Localy in TEC %s", Region[i]);
      if (pixelOutput)
        sprintf(htitle5, "Localy in BPIX %s", Region[i]);
      if (pixelOutput)
        sprintf(htitle6, "Localy in FPIX %s", Region[i]);

      sprintf(hname1, "Localy_TIB_%i", i + 1);
      sprintf(hname2, "Localy_TOB_%i", i + 1);
      sprintf(hname3, "Localy_TID_%i", i + 1);
      sprintf(hname4, "Localy_TEC_%i", i + 1);
      if (pixelOutput)
        sprintf(hname5, "Localy_BPIX_%i", i + 1);
      if (pixelOutput)
        sprintf(hname6, "Localy_FPIX_%i", i + 1);

      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TIBHit");
      h1ly[i] = ibooker.book1D(hname1, htitle1, nbin, low3[0], high3[0]);
      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TOBHit");
      h2ly[i] = ibooker.book1D(hname2, htitle2, nbin, low3[1], high3[1]);
      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TIDHit");
      h3ly[i] = ibooker.book1D(hname3, htitle3, nbin, low3[2], high3[2]);
      ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/TECHit");
      h4ly[i] = ibooker.book1D(hname4, htitle4, nbin, low3[3], high3[3]);
      if (pixelOutput) {
        ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/BPIXHit");
        h5ly[i] = ibooker.book1D(hname5, htitle5, nbin, low3[4], high3[4]);
        ibooker.setCurrentFolder("TrackerHitsV/TrackerHit/FPIXHit");
        h6ly[i] = ibooker.book1D(hname6, htitle6, nbin, low3[5], high3[5]);
      }
    }
  }
}

Member Data Documentation

conf_

edm::ParameterSet TrackerHitAnalyzer::conf_

private

Definition at line 53 of file TrackerHitAnalyzer.h.

edmPSimHitContainer_pxlBrlHigh_Token_

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_pxlBrlHigh_Token_

private

Definition at line 44 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edmPSimHitContainer_pxlBrlLow_Token_

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_pxlBrlLow_Token_

private

Definition at line 44 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edmPSimHitContainer_pxlFwdHigh_Token_

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_pxlFwdHigh_Token_

private

Definition at line 45 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edmPSimHitContainer_pxlFwdLow_Token_

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_pxlFwdLow_Token_

private

Definition at line 45 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edmPSimHitContainer_siTECHigh_Token_

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_siTECHigh_Token_

private

Definition at line 49 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edmPSimHitContainer_siTECLow_Token_

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_siTECLow_Token_

private

Definition at line 49 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edmPSimHitContainer_siTIBHigh_Token_

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_siTIBHigh_Token_

private

Definition at line 46 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edmPSimHitContainer_siTIBLow_Token_

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_siTIBLow_Token_

private

Definition at line 46 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edmPSimHitContainer_siTIDHigh_Token_

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_siTIDHigh_Token_

private

Definition at line 48 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edmPSimHitContainer_siTIDLow_Token_

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_siTIDLow_Token_

private

Definition at line 48 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edmPSimHitContainer_siTOBHigh_Token_

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_siTOBHigh_Token_

private

Definition at line 47 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edmPSimHitContainer_siTOBLow_Token_

edm::EDGetTokenT<edm::PSimHitContainer> TrackerHitAnalyzer::edmPSimHitContainer_siTOBLow_Token_

private

Definition at line 47 of file TrackerHitAnalyzer.h.

Referenced by analyze().

edmSimTrackContainerToken_

edm::EDGetTokenT<edm::SimTrackContainer> TrackerHitAnalyzer::edmSimTrackContainerToken_

private

Definition at line 50 of file TrackerHitAnalyzer.h.

Referenced by analyze().

fDBE

DQMStore* TrackerHitAnalyzer::fDBE

private

Definition at line 52 of file TrackerHitAnalyzer.h.

Referenced by bookHistograms().

fOutputFile

std::string TrackerHitAnalyzer::fOutputFile

private

Definition at line 106 of file TrackerHitAnalyzer.h.

h1e

MonitorElement* TrackerHitAnalyzer::h1e[12]

private

Definition at line 63 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h1ex

MonitorElement* TrackerHitAnalyzer::h1ex[12]

private

Definition at line 70 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h1ey

MonitorElement* TrackerHitAnalyzer::h1ey[12]

private

Definition at line 77 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h1ez

MonitorElement* TrackerHitAnalyzer::h1ez[12]

private

Definition at line 84 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h1lx

MonitorElement* TrackerHitAnalyzer::h1lx[12]

private

Definition at line 91 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h1ly

MonitorElement* TrackerHitAnalyzer::h1ly[12]

private

Definition at line 98 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h2e

MonitorElement* TrackerHitAnalyzer::h2e[12]

private

Definition at line 64 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h2ex

MonitorElement* TrackerHitAnalyzer::h2ex[12]

private

Definition at line 71 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h2ey

MonitorElement* TrackerHitAnalyzer::h2ey[12]

private

Definition at line 78 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h2ez

MonitorElement* TrackerHitAnalyzer::h2ez[12]

private

Definition at line 85 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h2lx

MonitorElement* TrackerHitAnalyzer::h2lx[12]

private

Definition at line 92 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h2ly

MonitorElement* TrackerHitAnalyzer::h2ly[12]

private

Definition at line 99 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h3e

MonitorElement* TrackerHitAnalyzer::h3e[12]

private

Definition at line 65 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h3ex

MonitorElement* TrackerHitAnalyzer::h3ex[12]

private

Definition at line 72 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h3ey

MonitorElement* TrackerHitAnalyzer::h3ey[12]

private

Definition at line 79 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h3ez

MonitorElement* TrackerHitAnalyzer::h3ez[12]

private

Definition at line 86 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h3lx

MonitorElement* TrackerHitAnalyzer::h3lx[12]

private

Definition at line 93 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h3ly

MonitorElement* TrackerHitAnalyzer::h3ly[12]

private

Definition at line 100 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h4e

MonitorElement* TrackerHitAnalyzer::h4e[12]

private

Definition at line 66 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h4ex

MonitorElement* TrackerHitAnalyzer::h4ex[12]

private

Definition at line 73 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h4ey

MonitorElement* TrackerHitAnalyzer::h4ey[12]

private

Definition at line 80 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h4ez

MonitorElement* TrackerHitAnalyzer::h4ez[12]

private

Definition at line 87 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h4lx

MonitorElement* TrackerHitAnalyzer::h4lx[12]

private

Definition at line 94 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h4ly

MonitorElement* TrackerHitAnalyzer::h4ly[12]

private

Definition at line 101 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h5e

MonitorElement* TrackerHitAnalyzer::h5e[12]

private

Definition at line 67 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h5ex

MonitorElement* TrackerHitAnalyzer::h5ex[12]

private

Definition at line 74 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h5ey

MonitorElement* TrackerHitAnalyzer::h5ey[12]

private

Definition at line 81 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h5ez

MonitorElement* TrackerHitAnalyzer::h5ez[12]

private

Definition at line 88 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h5lx

MonitorElement* TrackerHitAnalyzer::h5lx[12]

private

Definition at line 95 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h5ly

MonitorElement* TrackerHitAnalyzer::h5ly[12]

private

Definition at line 102 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h6e

MonitorElement* TrackerHitAnalyzer::h6e[12]

private

Definition at line 68 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h6ex

MonitorElement* TrackerHitAnalyzer::h6ex[12]

private

Definition at line 75 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h6ey

MonitorElement* TrackerHitAnalyzer::h6ey[12]

private

Definition at line 82 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h6ez

MonitorElement* TrackerHitAnalyzer::h6ez[12]

private

Definition at line 89 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h6lx

MonitorElement* TrackerHitAnalyzer::h6lx[12]

private

Definition at line 96 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

h6ly

MonitorElement* TrackerHitAnalyzer::h6ly[12]

private

Definition at line 103 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

htofeta

MonitorElement* TrackerHitAnalyzer::htofeta

private

Definition at line 55 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

htofeta_profile

MonitorElement* TrackerHitAnalyzer::htofeta_profile

private

Definition at line 59 of file TrackerHitAnalyzer.h.

htofphi

MonitorElement* TrackerHitAnalyzer::htofphi

private

Definition at line 56 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

htofphi_profile

MonitorElement* TrackerHitAnalyzer::htofphi_profile

private

Definition at line 60 of file TrackerHitAnalyzer.h.

htofr

MonitorElement* TrackerHitAnalyzer::htofr

private

Definition at line 57 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

htofr_profile

MonitorElement* TrackerHitAnalyzer::htofr_profile

private

Definition at line 61 of file TrackerHitAnalyzer.h.

htofz

MonitorElement* TrackerHitAnalyzer::htofz

private

Definition at line 58 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

htofz_profile

MonitorElement* TrackerHitAnalyzer::htofz_profile

private

Definition at line 62 of file TrackerHitAnalyzer.h.

pixelOutput

bool TrackerHitAnalyzer::pixelOutput

private

Definition at line 107 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

runStandalone

bool TrackerHitAnalyzer::runStandalone

private

Definition at line 105 of file TrackerHitAnalyzer.h.

tGeomEsToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> TrackerHitAnalyzer::tGeomEsToken_

private

Definition at line 41 of file TrackerHitAnalyzer.h.

Referenced by analyze().

verbose_

bool TrackerHitAnalyzer::verbose_

private

Definition at line 42 of file TrackerHitAnalyzer.h.