#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 &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &)=0

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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

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)

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::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 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 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 24 of file TrackerHitAnalyzer.h.

Constructor & Destructor Documentation

TrackerHitAnalyzer::TrackerHitAnalyzer ( const edm::ParameterSet & ps )

Constructor.

Definition at line 41 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 326 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 345 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
 
   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 63 of file TrackerHitAnalyzer.cc.

References DQMStore::IBooker::book1D(), 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, htofeta, htofphi, htofr, htofz, mps_fire::i, Utilities::operator, pixelOutput, DQMStore::IBooker::setCurrentFolder(), DQMStore::setVerbose(), 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->setVerbose(1);
      } else {
        fDBE->setVerbose(0);
      }
   }
                                                                         
   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 333 of file TrackerHitAnalyzer.cc.

References fDBE, fOutputFile, runStandalone, and 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);}
 }