#include <SiPixelMuonHLT.h>

Inheritance diagram for SiPixelMuonHLT:

Classes
struct	LayerMEs

Public Types
typedef edmNew::DetSet < SiPixelCluster > ::const_iterator	ClusterIterator

Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

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

virtual void	beginJob ()

virtual void	endJob ()

virtual void	Histo_init ()

	SiPixelMuonHLT (const edm::ParameterSet &conf)

	~SiPixelMuonHLT ()

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
edm::InputTag	clusterCollectionTag_

edm::EDGetTokenT < edmNew::DetSetVector < SiPixelCluster > >	clustersToken_

edm::ParameterSet	conf_

int	counterEvt_

int	eventNo

SiStripFolderOrganizer	folder_organizer

edm::InputTag	l3MuonCollectionTag_

edm::EDGetTokenT < reco::RecoChargedCandidateCollection >	l3MuonCollectionToken_

std::map< int, MonitorElement * >	MEContainerAllBarrelEta

std::map< int, MonitorElement * >	MEContainerAllBarrelEtaPhi

std::map< int, MonitorElement * >	MEContainerAllBarrelN

std::map< int, MonitorElement * >	MEContainerAllBarrelPhi

std::map< int, MonitorElement * >	MEContainerAllBarrelZ

std::map< int, MonitorElement * >	MEContainerAllBarrelZPhi

std::map< int, MonitorElement * >	MEContainerAllEndcapN

std::map< int, MonitorElement * >	MEContainerAllEndcapPhi

std::map< int, MonitorElement * >	MEContainerAllEndcapXY

std::map< int, MonitorElement * >	MEContainerOnTrackBarrelEta

std::map< int, MonitorElement * >	MEContainerOnTrackBarrelEtaPhi

std::map< int, MonitorElement * >	MEContainerOnTrackBarrelN

std::map< int, MonitorElement * >	MEContainerOnTrackBarrelPhi

std::map< int, MonitorElement * >	MEContainerOnTrackBarrelZ

std::map< int, MonitorElement * >	MEContainerOnTrackBarrelZPhi

std::map< int, MonitorElement * >	MEContainerOnTrackEndcapN

std::map< int, MonitorElement * >	MEContainerOnTrackEndcapPhi

std::map< int, MonitorElement * >	MEContainerOnTrackEndcapXY

std::string	monitorName_

std::string	outputFile_

edm::ParameterSet	parameters_

int	prescaleEvt_
	counter More...

edm::InputTag	rechitsCollectionTag_

edm::EDGetTokenT < edmNew::DetSetVector < SiPixelRecHit > >	rechitsToken_

bool	saveOUTput_

edm::InputTag	src_

DQMStore *	theDMBE

bool	verbose_
	every n events More...

Additional Inherited Members
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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)

Detailed Description

Definition at line 84 of file SiPixelMuonHLT.h.

Member Typedef Documentation

typedef edmNew::DetSet<SiPixelCluster>::const_iterator SiPixelMuonHLT::ClusterIterator

Definition at line 94 of file SiPixelMuonHLT.h.

Constructor & Destructor Documentation

SiPixelMuonHLT::SiPixelMuonHLT ( const edm::ParameterSet & conf )

explicit

Definition at line 24 of file SiPixelMuonHLT.cc.

References clusterCollectionTag_, clustersToken_, edm::ParameterSet::getUntrackedParameter(), Histo_init(), l3MuonCollectionTag_, l3MuonCollectionToken_, monitorName_, NULL, cppFunctionSkipper::operator, outputFile_, parameters_, rechitsCollectionTag_, rechitsToken_, saveOUTput_, AlCaHLTBitMon_QueryRunRegistry::string, theDMBE, and verbose_.

                                                              :
   conf_(iConfig)
 {
 
   parameters_ = iConfig;
 
   verbose_ = parameters_.getUntrackedParameter < bool > ("verbose", false);
   monitorName_ = parameters_.getUntrackedParameter < std::string > ("monitorName", "HLT/HLTMonMuon");
   saveOUTput_  = parameters_.getUntrackedParameter < bool > ("saveOUTput", true);
 
   //tags
   clusterCollectionTag_ = parameters_.getUntrackedParameter < edm::InputTag > ("clusterCollectionTag", edm::InputTag ("hltSiPixelClusters"));
   rechitsCollectionTag_ = parameters_.getUntrackedParameter < edm::InputTag > ("rechitsCollectionTag", edm::InputTag ("hltSiPixelRecHits"));
   l3MuonCollectionTag_  = parameters_.getUntrackedParameter < edm::InputTag > ("l3MuonCollectionTag", edm::InputTag ("hltL3MuonCandidates"));
 
    theDMBE = edm::Service<DQMStore>().operator->();
    edm::LogInfo ("PixelHLTDQM") << "SiPixelMuonHLT::SiPixelMuonHLT: Got DQM BackEnd interface"<<std::endl;
    outputFile_ = parameters_.getUntrackedParameter < std::string > ("outputFile", "");
    if (outputFile_.size () != 0)
      edm::LogWarning ("HLTMuonDQMSource") << "Muon HLT Monitoring histograms will be saved to " << outputFile_ << std::endl;
    else
      outputFile_ = "PixelHLTDQM.root";
    if (theDMBE != NULL) theDMBE->setCurrentFolder (monitorName_); 
    SiPixelMuonHLT::Histo_init();
 
    //set Token(-s)
    clustersToken_ = consumes<edmNew::DetSetVector<SiPixelCluster> >(std::string("hltSiPixelClusters"));
    rechitsToken_ = consumes<edmNew::DetSetVector<SiPixelRecHit> >(std::string("hltSiPixelRecHits"));
    l3MuonCollectionToken_ = consumes<reco::RecoChargedCandidateCollection>(parameters_.getUntrackedParameter < edm::InputTag > ("l3MuonCollectionTag", edm::InputTag ("hltL3MuonCandidates")));
 }

SiPixelMuonHLT::~SiPixelMuonHLT ( )

Definition at line 57 of file SiPixelMuonHLT.cc.

 {
    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)
   edm::LogInfo ("PixelHLTDQM") << "SiPixelMuonHLT::~SiPixelMuonHLT: Destructor"<<std::endl;
 
 }

Member Function Documentation

void SiPixelMuonHLT::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

virtual

Implements edm::EDAnalyzer.

Definition at line 83 of file SiPixelMuonHLT.cc.

References SiPixelCluster::charge(), clustersToken_, PV3DBase< T, PVType, FrameType >::eta(), edm::EventAuxiliary::event(), edm::Event::eventAuxiliary(), eventNo, edm::HandleBase::failedToGet(), HcalObjRepresent::Fill(), GeomDet::geographicalId(), edm::EventSetup::get(), edm::Ref< C, T, F >::get(), edm::Event::getByToken(), i, edm::HandleBase::id(), TrackerGeometry::idToDet(), edm::HandleBase::isValid(), l3MuonCollectionToken_, MEContainerAllBarrelEta, MEContainerAllBarrelEtaPhi, MEContainerAllBarrelN, MEContainerAllBarrelPhi, MEContainerAllBarrelZ, MEContainerAllBarrelZPhi, MEContainerAllEndcapN, MEContainerAllEndcapPhi, MEContainerAllEndcapXY, MEContainerOnTrackBarrelEta, MEContainerOnTrackBarrelEtaPhi, MEContainerOnTrackBarrelN, MEContainerOnTrackBarrelPhi, MEContainerOnTrackBarrelZ, MEContainerOnTrackBarrelZPhi, MEContainerOnTrackEndcapN, MEContainerOnTrackEndcapPhi, MEContainerOnTrackEndcapXY, PV3DBase< T, PVType, FrameType >::phi(), edm::ESHandle< class >::product(), TrackerTopology::pxbLayer(), TrackerTopology::pxfDisk(), TrackerTopology::pxfSide(), HI_PhotonSkim_cff::rechits, rechitsToken_, edm::EventAuxiliary::run(), PixelGeomDetUnit::specificTopology(), DetId::subdetId(), GeomDet::surface(), Surface::toGlobal(), DetId::Tracker, PV3DBase< T, PVType, FrameType >::x(), SiPixelCluster::x(), PV3DBase< T, PVType, FrameType >::y(), SiPixelCluster::y(), and PV3DBase< T, PVType, FrameType >::z().

 {
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   iSetup.get<IdealGeometryRecord>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
 
 
   eventNo++;
 
   edm::ESHandle < TrackerGeometry > TG;
   iSetup.get < TrackerDigiGeometryRecord > ().get (TG);
   const TrackerGeometry *theTrackerGeometry = TG.product ();
   const TrackerGeometry & theTracker (*theTrackerGeometry);
   edm::Handle < reco::RecoChargedCandidateCollection > l3mucands;
   reco::RecoChargedCandidateCollection::const_iterator cand;
   edm::Handle <edmNew::DetSetVector<SiPixelCluster> > clusters;
   edm::Handle <edmNew::DetSetVector<SiPixelRecHit> > rechits;
 
   bool GotClusters = true;
   bool GotRecHits  = true;
   bool GotL3Muons  = true;
   
   iEvent.getByToken(clustersToken_, clusters);
   if(!clusters.isValid()){
     edm::LogInfo("PixelHLTDQM") << "No pix clusters, cannot run for event " << iEvent.eventAuxiliary ().event() <<" run: "<<iEvent.eventAuxiliary ().run()  << std::endl;
     GotClusters = false;
   }
   iEvent.getByToken(rechitsToken_, rechits);
   if(!rechits.isValid()){
     edm::LogInfo("PixelHLTDQM") << "No pix rechits, cannot run for event " << iEvent.eventAuxiliary ().event() <<" run: "<<iEvent.eventAuxiliary ().run()  << std::endl;
     GotRecHits = false;
   }
   iEvent.getByToken(l3MuonCollectionToken_, l3mucands);
   if(!l3mucands.isValid()){
     edm::LogInfo("PixelHLTDQM") << "No L3 Muons, cannot run for event " << iEvent.eventAuxiliary ().event() <<" run: "<<iEvent.eventAuxiliary ().run()  << std::endl;
     GotL3Muons = false;
   }
   
   if (GotClusters){
     if(!clusters.failedToGet ())
       {
     int NBarrel[4] = {0,0,0,0};
     int NEndcap[5] = {0,0,0,0,0};
     for (size_t i = 0; i < clusters->size(); ++i){ 
       const SiPixelCluster* clust = clusters->data(i);
       clust->charge();
       uint detID = clusters->id(i);
       const PixelGeomDetUnit *PixGeom = dynamic_cast < const PixelGeomDetUnit * >(theTracker.idToDet (detID));
       const PixelTopology *topol = dynamic_cast < const PixelTopology * >(&(PixGeom->specificTopology ()));
       // get the cluster position in local coordinates (cm)
       LocalPoint clustlp = topol->localPosition (MeasurementPoint(clust->x(),clust->y()));
       GlobalPoint clustgp = PixGeom->surface ().toGlobal (clustlp);
       if(PixGeom->geographicalId().subdetId() == 1){ //1 Defines a barrel hit
         int clustLay = tTopo->pxbLayer(detID);
         //Eta-Phi
         MEContainerAllBarrelEtaPhi[0]->Fill(clustgp.eta(),clustgp.phi());
         MEContainerAllBarrelZPhi[0]->Fill(clustgp.z(),clustgp.phi());
         MEContainerAllBarrelEtaPhi[clustLay]->Fill(clustgp.eta(),clustgp.phi());
         MEContainerAllBarrelZPhi[clustLay]->Fill(clustgp.z(),clustgp.phi());
         //Eta
         MEContainerAllBarrelEta[0]->Fill(clustgp.eta());
         MEContainerAllBarrelZ[0]->Fill(clustgp.z());
         MEContainerAllBarrelEta[clustLay]->Fill(clustgp.eta());
         MEContainerAllBarrelZ[clustLay]->Fill(clustgp.z());
         //Phi
         MEContainerAllBarrelPhi[0]->Fill(clustgp.phi());
         MEContainerAllBarrelPhi[clustLay]->Fill(clustgp.phi());
         ++NBarrel[0]; //N clusters all layers
         ++NBarrel[clustLay]; //N clusters all layers
       }
       if(PixGeom->geographicalId().subdetId() == 2){ //2 Defines a Endcap hit
         int clustDisk  = tTopo->pxfDisk(detID);
         if( tTopo->pxfSide(detID) == 2)
           clustDisk = clustDisk +2;//neg z disks have ID 3 and 4
         MEContainerAllEndcapXY[0]->Fill(clustgp.x(),clustgp.y());
         MEContainerAllEndcapXY[clustDisk]->Fill(clustgp.x(),clustgp.y());
         MEContainerAllEndcapPhi[0]->Fill(clustgp.phi());
         MEContainerAllEndcapPhi[clustDisk]->Fill(clustgp.phi());
         ++NEndcap[0];
         ++NEndcap[clustDisk];
       }
       
     }
     MEContainerAllBarrelN[0]->Fill(NBarrel[0]);
     for (int lay = 1; lay < 4; ++lay)
       MEContainerAllBarrelN[lay]->Fill(NBarrel[lay]);
     MEContainerAllEndcapN[0]->Fill(NEndcap[0]);
     for (int disk = 1; disk < 5; ++disk)
       MEContainerAllEndcapN[disk]->Fill(NEndcap[disk]); 
       }//if clust (!failedToGet)
   }
   bool doRecHits = false;
   
   if (GotRecHits && doRecHits){
     if(!rechits.failedToGet ())
       {
     for (size_t i = 0; i < rechits->size(); ++i){ 
       const SiPixelRecHit* myhit = rechits->data(i);
       uint detID = rechits->id(i);
       const PixelGeomDetUnit *PixGeom = dynamic_cast < const PixelGeomDetUnit * >(theTracker.idToDet (detID));
       //edm::LogInfo("PixelHLTDQM") << "" << PixGeom->geographicalId().subdetId() << std::endl;
       //const PixelTopology *topol = dynamic_cast < const PixelTopology * >(&(PixGeom->specificTopology ()));
       // get the hit position in local coordinates (cm)
       //LocalPoint hitlp = topol->localPosition (MeasurementPoint(myhit->x(),myhit->y()));
       if(PixGeom->geographicalId().subdetId() == 1 && myhit->hasPositionAndError()){
         GlobalPoint hitgp = PixGeom->surface ().toGlobal (myhit->localPosition());
         edm::LogInfo("PixelHLTDQM") << " (From SiPixelRecHit) Hit Eta: " << hitgp.eta()   << " Hit Phi: " << hitgp.phi()  << std::endl;
       }
     }      
       }
   }
   if(GotL3Muons){
     if(!l3mucands.failedToGet ())
       {
     int NBarrel[4] = {0,0,0,0};
     int NEndcap[5] = {0,0,0,0,0};
     for (cand = l3mucands->begin (); cand != l3mucands->end (); ++cand){
       reco::TrackRef l3tk = cand->get < reco::TrackRef > ();
       for (size_t hit = 0; hit < l3tk->recHitsSize (); hit++){
         if (l3tk->recHit (hit)->isValid () == true && l3tk->recHit (hit)->geographicalId ().det () == DetId::Tracker){
           int detID = l3tk->recHit(hit)->geographicalId().rawId();
           //if hit is in pixel detector say true
           bool IdMatch = typeid(*(l3tk->recHit(hit))) == typeid(SiPixelRecHit);
           if (IdMatch){
         const SiPixelRecHit *pixhit = dynamic_cast < const SiPixelRecHit * >(l3tk->recHit(hit).get());
         if((*pixhit).isValid() == true){
           edm::Ref<edmNew::DetSetVector<SiPixelCluster>, SiPixelCluster> const& pixclust = (*pixhit).cluster();
           if (!(*pixhit).cluster().isAvailable()) 
             {continue;}
           const PixelGeomDetUnit *PixGeom = dynamic_cast < const PixelGeomDetUnit * >(theTracker.idToDet (detID));
           const PixelTopology *topol = dynamic_cast < const PixelTopology * >(&(PixGeom->specificTopology ()));
           LocalPoint clustlp = topol->localPosition (MeasurementPoint(pixclust->x(),pixclust->y()));
           GlobalPoint clustgp = PixGeom->surface ().toGlobal (clustlp);
           if(l3tk->recHit(hit)->geographicalId().subdetId() == 1){ //1 Defines a barrel hit
             //get the cluster position in local coordinates (cm)      
             int clustLay = tTopo->pxbLayer(detID);
             MEContainerOnTrackBarrelEtaPhi[0]->Fill(clustgp.eta(),clustgp.phi());
             MEContainerOnTrackBarrelZPhi[0]->Fill(clustgp.z(),clustgp.phi());
             MEContainerOnTrackBarrelEtaPhi[clustLay]->Fill(clustgp.eta(),clustgp.phi());
             MEContainerOnTrackBarrelZPhi[clustLay]->Fill(clustgp.z(),clustgp.phi());
             MEContainerOnTrackBarrelEta[0]->Fill(clustgp.eta());
             MEContainerOnTrackBarrelZ[0]->Fill(clustgp.z());
             MEContainerOnTrackBarrelEta[clustLay]->Fill(clustgp.eta());
             MEContainerOnTrackBarrelZ[clustLay]->Fill(clustgp.z());
             MEContainerOnTrackBarrelPhi[0]->Fill(clustgp.phi());
             MEContainerOnTrackBarrelPhi[clustLay]->Fill(clustgp.phi());
             ++NBarrel[0];
             ++NBarrel[clustLay];
           }//subdet ==1
           if(l3tk->recHit(hit)->geographicalId().subdetId() == 2){ //2 Defines a Endcap hit
             int clustDisk  = tTopo->pxfDisk(detID);
             if( tTopo->pxfDisk(detID) == 2)
               clustDisk = clustDisk +2;
             MEContainerOnTrackEndcapXY[0]->Fill(clustgp.x(),clustgp.y());
             MEContainerOnTrackEndcapXY[clustDisk]->Fill(clustgp.x(),clustgp.y());
             MEContainerOnTrackEndcapPhi[0]->Fill(clustgp.phi());
             MEContainerOnTrackEndcapPhi[clustDisk]->Fill(clustgp.phi());
             ++NEndcap[0];
             ++NEndcap[clustDisk];
           }//subdet ==2
         }//pixhit valid
           }//typeid match
         }//l3tk->recHit (hit)->isValid () == true
       }//loop over RecHits
     }//loop over l3mucands
     MEContainerOnTrackBarrelN[0]->Fill(NBarrel[0]);
     for (int lay = 1; lay < 4; ++lay)
       MEContainerOnTrackBarrelN[lay]->Fill(NBarrel[lay]);
     MEContainerOnTrackEndcapN[0]->Fill(NEndcap[0]);
     for (int disk = 1; disk < 5; ++disk)
       MEContainerOnTrackEndcapN[disk]->Fill(NEndcap[disk]);
     
       }//if l3mucands  
   }
 }

void SiPixelMuonHLT::beginJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 65 of file SiPixelMuonHLT.cc.

References eventNo.

                              {
 
   edm::LogInfo ("PixelHLTDQM") << " SiPixelMuonHLT::beginJob - Initialisation ... " << std::endl;
   eventNo = 0;
 
 }

void SiPixelMuonHLT::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 73 of file SiPixelMuonHLT.cc.

References outputFile_, DQMStore::save(), saveOUTput_, and theDMBE.

                                {
   if(saveOUTput_){
     edm::LogInfo ("PixelHLTDQM") << " SiPixelMuonHLT::endJob - Saving Root File " << std::endl;
     theDMBE->save( outputFile_.c_str() );
   }
 }

void SiPixelMuonHLT::Histo_init ( )

virtual

Definition at line 262 of file SiPixelMuonHLT.cc.

References DQMStore::book1D(), DQMStore::book2D(), Gflash::EtaMax, i, MEContainerAllBarrelEta, MEContainerAllBarrelEtaPhi, MEContainerAllBarrelN, MEContainerAllBarrelPhi, MEContainerAllBarrelZ, MEContainerAllBarrelZPhi, MEContainerAllEndcapN, MEContainerAllEndcapPhi, MEContainerAllEndcapXY, MEContainerOnTrackBarrelEta, MEContainerOnTrackBarrelEtaPhi, MEContainerOnTrackBarrelN, MEContainerOnTrackBarrelPhi, MEContainerOnTrackBarrelZ, MEContainerOnTrackBarrelZPhi, MEContainerOnTrackEndcapN, MEContainerOnTrackEndcapPhi, MEContainerOnTrackEndcapXY, monitorName_, DQMStore::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, theDMBE, and indexGen::title.

Referenced by SiPixelMuonHLT().

 {
    monitorName_ = monitorName_+"/SiPixel";
    int NBinsEta = 100;
    int NBinsPhi = 80;
    float EtaMax = 3.0;
    float ZMax   = 28.0;
    int NBinsZ   = 112;
    float PhiMax = 3.142;
    int   NBinsN = 800;
    float NMax   = 800.;
    int   NBinsX = 100;
    int   NBinsY = 100;
    float XMax   = 20.;
    float YMax   = 20.;
    std::string histoname;
    std::string title;
 
    theDMBE->setCurrentFolder (monitorName_ + "/Barrel");   
    std::string layerLabel[4] = {"All_Layers", "Layer1", "Layer2", "Layer3"};
    for (unsigned int i = 0; i < 4; i++)
      {
        //Eta-Phi
        histoname = "EtaPhiAllBarrelMap_" + layerLabel[i];
        title     = "#eta-#phi Pixel Barrel Cluster Occupancy Map for " + layerLabel[i];
        MEContainerAllBarrelEtaPhi[i] = theDMBE->book2D (histoname, title, NBinsEta, -EtaMax, EtaMax, NBinsPhi, -PhiMax, PhiMax);
        //Z-Phi
        histoname = "ZPhiAllBarrelMap_" + layerLabel[i];
        title     = "Z-#phi Pixel Barrel Cluster Occupancy Map for " + layerLabel[i];
        MEContainerAllBarrelZPhi[i] = theDMBE->book2D (histoname, title, NBinsZ, -ZMax, ZMax, NBinsPhi, -PhiMax, PhiMax);
        //Eta
        histoname = "EtaAllBarrelMap_" + layerLabel[i];
        title     = "#eta Pixel Barrel Cluster Occupancy Map for " + layerLabel[i];
        MEContainerAllBarrelEta[i] = theDMBE->book1D (histoname, title, NBinsEta, -EtaMax, EtaMax);
        //Z
        histoname = "ZAllBarrelMap_" + layerLabel[i];
        title     = "Z Pixel Barrel Cluster Occupancy Map for " + layerLabel[i];
        MEContainerAllBarrelZ[i] = theDMBE->book1D (histoname, title, NBinsZ, -ZMax, ZMax);
        //Phi
        histoname = "PhiAllBarrelMap_" + layerLabel[i];
        title     = "#phi Pixel Barrel Cluster Occupancy Map for " + layerLabel[i];
        MEContainerAllBarrelPhi[i] = theDMBE->book1D (histoname, title, NBinsPhi, -PhiMax, PhiMax);
        //N clusters
        histoname = "NAllBarrelMap_" + layerLabel[i];
        title     = "#phi Pixel Barrel Cluster Occupancy Map for " + layerLabel[i];
        MEContainerAllBarrelN[i] = theDMBE->book1D (histoname, title, NBinsN, 0, NMax);
        //Eta-Phi
        histoname = "EtaPhiOnTrackBarrelMap_" + layerLabel[i];
        title = "#eta-#phi On Track Pixel Barrel Cluster Occupancy Map for " + layerLabel[i];
        MEContainerOnTrackBarrelEtaPhi[i] = theDMBE->book2D (histoname, title, NBinsEta, -EtaMax, EtaMax, NBinsPhi, -PhiMax, PhiMax);
        //Z-Phi
        histoname = "ZPhiOnTrackBarrelMap_" + layerLabel[i];
        title = "Z-#phi On Track Pixel Barrel Cluster Occupancy Map for " + layerLabel[i];
        MEContainerOnTrackBarrelZPhi[i] = theDMBE->book2D (histoname, title, NBinsZ, -ZMax, ZMax, NBinsPhi, -PhiMax, PhiMax);
        //Eta
        histoname = "EtaOnTrackBarrelMap_" + layerLabel[i];
        title     = "#eta On Track Pixel Barrel Cluster Occupancy Map for " + layerLabel[i];
        MEContainerOnTrackBarrelEta[i] = theDMBE->book1D (histoname, title, NBinsEta, -EtaMax, EtaMax);
        //Z
        histoname = "ZOnTrackBarrelMap_" + layerLabel[i];
        title     = "Z On Track Pixel Barrel Cluster Occupancy Map for " + layerLabel[i];
        MEContainerOnTrackBarrelZ[i] = theDMBE->book1D (histoname, title, NBinsZ, -ZMax, ZMax);
        //Phi
        histoname = "PhiOnTrackBarrelMap_" + layerLabel[i];
        title     = "#phi On Track Pixel Barrel Cluster Occupancy Map for " + layerLabel[i];
        MEContainerOnTrackBarrelPhi[i] = theDMBE->book1D (histoname, title, NBinsPhi, -PhiMax, PhiMax);
        //N clusters
        histoname = "NOnTrackBarrelMap_" + layerLabel[i];
        title     = "N_{Clusters} On Track Pixel Barrel Cluster Occupancy Map for " + layerLabel[i];
        MEContainerOnTrackBarrelN[i] = theDMBE->book1D (histoname, title, NBinsN, 0, NMax);
      }
 
    theDMBE->setCurrentFolder (monitorName_ + "/EndCap");
    std::string diskLabel[5] = {"All_Disks", "InnerPosZ", "OuterPosZ", "InnerNegZ", "OuterNegZ"};
    for (int i = 0;i < 5; ++i) 
      {
        //XY
        histoname = "XYAllEndcapMap_" + diskLabel[i];
        title     = "X-Y Pixel Endcap Cluster Occupancy Map for " + diskLabel[i];
        MEContainerAllEndcapXY[i] = theDMBE->book2D (histoname, title, NBinsX, -XMax, XMax, NBinsY, -YMax, YMax);
        //Phi
        histoname = "PhiAllEndcapMap_" + diskLabel[i];
        title     = "#phi Pixel Endcap Cluster Occupancy Map for " + diskLabel[i];
        MEContainerAllEndcapPhi[i] = theDMBE->book1D (histoname, title, NBinsPhi, -PhiMax, PhiMax);
        //N clusters
        histoname = "NAllEndcapMap_" + diskLabel[i];
        title     = "#phi Pixel Endcap Cluster Occupancy Map for " + diskLabel[i];
        MEContainerAllEndcapN[i] = theDMBE->book1D (histoname, title, NBinsN, 0, NMax);
        //XY
        histoname = "XYOnTrackEndcapMap_" + diskLabel[i];
        title     = "X-Y Pixel Endcap On Track Cluster Occupancy Map for " + diskLabel[i];
        MEContainerOnTrackEndcapXY[i] = theDMBE->book2D (histoname, title, NBinsX, -XMax, XMax, NBinsY, -YMax, YMax);
        //Phi
        histoname = "PhiOnTrackEndcapMap_" + diskLabel[i];
        title     = "#phi Pixel Endcap On Track Cluster Occupancy Map for " + diskLabel[i];
        MEContainerOnTrackEndcapPhi[i] = theDMBE->book1D (histoname, title, NBinsPhi, -PhiMax, PhiMax);
        //N clusters
        histoname = "NOnTrackEndcapMap_" + diskLabel[i];
        title     = "#phi Pixel Endcap On Track Cluster Occupancy Map for " + diskLabel[i];
        MEContainerOnTrackEndcapN[i] = theDMBE->book1D (histoname, title, NBinsN, 0, NMax);
    }
    return;
 }