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
virtual void	analyze (edm::Event const &, edm::EventSetup const &)
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)=0
virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &)
void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &, edm::EventSetup const &) final
void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	EDProducer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDProducerBase () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

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
void	endJob () override
Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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
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::one::EDProducerBase
typedef EDProducerBase	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 22 of file TrackerHitAnalyzer.h.

Constructor & Destructor Documentation

TrackerHitAnalyzer::TrackerHitAnalyzer

(

const edm::ParameterSet &

ps

)

Constructor.

Definition at line 40 of file TrackerHitAnalyzer.cc.

    : 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 ( )

override

Destructor.

Definition at line 336 of file TrackerHitAnalyzer.cc.

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

Member Function Documentation

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

overrideprotected

Analyze.

Definition at line 353 of file TrackerHitAnalyzer.cc.

References DEFINE_FWK_MODULE, Geom::Phi< T1, Range >::degrees(), 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(), edm::EventID::event(), dqm::impl::MonitorElement::Fill(), edm::EventSetup::get(), edm::Event::getByToken(), 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::EventBase::id(), TrackerGeometry::idToDetUnit(), edm::HandleBase::isValid(), dqm-mbProfile::log, PV3DBase< T, PVType, FrameType >::mag(), Max(), Min(), AlCaHLTBitMon_ParallelJobs::p, PV3DBase< T, PVType, FrameType >::phi(), pixelOutput, edm::EventID::run(), mathSSE::sqrt(), funct::tan(), 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

  edm::ESHandle<TrackerGeometry> tracker;
  c.get<TrackerDigiGeometryRecord>().get(tracker);

  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;
}

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

overrideprotected

Definition at line 73 of file TrackerHitAnalyzer.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), dqm::dqmstoreimpl::DQMStore::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, ZMuMuCategoriesSequences_cff::nbin, Utilities::operator, pixelOutput, g4SimHits_cfi::Region, dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), dqm::dqmstoreimpl::DQMStore::showDirStructure(), and verbose_.

                                                                                                            {
  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) {
    if (verbose_)
      fDBE->showDirStructure();
  }

  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]);
      }
    }
  }
}

void TrackerHitAnalyzer::endJob ( void  )

overrideprotectedvirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 340 of file TrackerHitAnalyzer.cc.

References fDBE, fOutputFile, runStandalone, and dqm::dqmstoreimpl::DQMStore::save().

                                {
  // According to the previous code some profile plots were created here
  // However, these profile plots are not in the final root file
  // For now we comment out these plots (since they are not created in any case)
  // Then, if needed we will consider moving the booking of the profileplots to
  // the bookHistograms function and here we will do the profile

  // Save root file only in standalone mode
  if (runStandalone && !fOutputFile.empty() && fDBE) {
    fDBE->save(fOutputFile);
  }
}

Member Data Documentation

edm::ParameterSet TrackerHitAnalyzer::conf_

private

Definition at line 53 of file TrackerHitAnalyzer.h.

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

private

Definition at line 44 of file TrackerHitAnalyzer.h.

Referenced by analyze().

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

private

Definition at line 44 of file TrackerHitAnalyzer.h.

Referenced by analyze().

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

private

Definition at line 45 of file TrackerHitAnalyzer.h.

Referenced by analyze().

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

private

Definition at line 45 of file TrackerHitAnalyzer.h.

Referenced by analyze().

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

private

Definition at line 49 of file TrackerHitAnalyzer.h.

Referenced by analyze().

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

private

Definition at line 49 of file TrackerHitAnalyzer.h.

Referenced by analyze().

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

private

Definition at line 46 of file TrackerHitAnalyzer.h.

Referenced by analyze().

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

private

Definition at line 46 of file TrackerHitAnalyzer.h.

Referenced by analyze().

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

private

Definition at line 48 of file TrackerHitAnalyzer.h.

Referenced by analyze().

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

private

Definition at line 48 of file TrackerHitAnalyzer.h.

Referenced by analyze().

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

private

Definition at line 47 of file TrackerHitAnalyzer.h.

Referenced by analyze().

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

private

Definition at line 47 of file TrackerHitAnalyzer.h.

Referenced by analyze().

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

private

Definition at line 50 of file TrackerHitAnalyzer.h.

Referenced by analyze().

DQMStore* TrackerHitAnalyzer::fDBE

private

Definition at line 52 of file TrackerHitAnalyzer.h.

Referenced by bookHistograms(), and endJob().

std::string TrackerHitAnalyzer::fOutputFile

private

Definition at line 106 of file TrackerHitAnalyzer.h.

Referenced by endJob().

MonitorElement* TrackerHitAnalyzer::h1e[12]

private

Definition at line 63 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h1ex[12]

private

Definition at line 70 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h1ey[12]

private

Definition at line 77 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h1ez[12]

private

Definition at line 84 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h1lx[12]

private

Definition at line 91 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h1ly[12]

private

Definition at line 98 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h2e[12]

private

Definition at line 64 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h2ex[12]

private

Definition at line 71 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h2ey[12]

private

Definition at line 78 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h2ez[12]

private

Definition at line 85 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h2lx[12]

private

Definition at line 92 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h2ly[12]

private

Definition at line 99 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h3e[12]

private

Definition at line 65 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h3ex[12]

private

Definition at line 72 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h3ey[12]

private

Definition at line 79 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h3ez[12]

private

Definition at line 86 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h3lx[12]

private

Definition at line 93 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h3ly[12]

private

Definition at line 100 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h4e[12]

private

Definition at line 66 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h4ex[12]

private

Definition at line 73 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h4ey[12]

private

Definition at line 80 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h4ez[12]

private

Definition at line 87 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h4lx[12]

private

Definition at line 94 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h4ly[12]

private

Definition at line 101 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h5e[12]

private

Definition at line 67 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h5ex[12]

private

Definition at line 74 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h5ey[12]

private

Definition at line 81 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h5ez[12]

private

Definition at line 88 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h5lx[12]

private

Definition at line 95 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h5ly[12]

private

Definition at line 102 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h6e[12]

private

Definition at line 68 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h6ex[12]

private

Definition at line 75 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h6ey[12]

private

Definition at line 82 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h6ez[12]

private

Definition at line 89 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h6lx[12]

private

Definition at line 96 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::h6ly[12]

private

Definition at line 103 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::htofeta

private

Definition at line 55 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::htofeta_profile

private

Definition at line 59 of file TrackerHitAnalyzer.h.

MonitorElement* TrackerHitAnalyzer::htofphi

private

Definition at line 56 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::htofphi_profile

private

Definition at line 60 of file TrackerHitAnalyzer.h.

MonitorElement* TrackerHitAnalyzer::htofr

private

Definition at line 57 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::htofr_profile

private

Definition at line 61 of file TrackerHitAnalyzer.h.

MonitorElement* TrackerHitAnalyzer::htofz

private

Definition at line 58 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* TrackerHitAnalyzer::htofz_profile

private

Definition at line 62 of file TrackerHitAnalyzer.h.

bool TrackerHitAnalyzer::pixelOutput

private

Definition at line 107 of file TrackerHitAnalyzer.h.

Referenced by analyze(), and bookHistograms().

bool TrackerHitAnalyzer::runStandalone

private

Definition at line 105 of file TrackerHitAnalyzer.h.

Referenced by endJob().

bool TrackerHitAnalyzer::verbose_

private

Definition at line 42 of file TrackerHitAnalyzer.h.

Referenced by bookHistograms().