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

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

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

overrideprotectedvirtual

Analyze.

Reimplemented from DQMEDAnalyzer.

Definition at line 334 of file TrackerHitAnalyzer.cc.

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

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, htofeta, htofphi, htofr, htofz, mps_fire::i, Utilities::operator, pixelOutput, 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]);
       }
     }
   }
 }