#include <SiPixelMuonHLT.h>

Inheritance diagram for SiPixelMuonHLT:

Classes
struct	LayerMEs
Public Types
typedef edmNew::DetSet < SiPixelCluster > ::const_iterator	ClusterIterator
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 ()
Private Attributes
edm::InputTag	clusterCollectionTag_
edm::ParameterSet	conf_
int	counterEvt_
int	eventNo
SiStripFolderOrganizer	folder_organizer
edm::InputTag	l3MuonCollectionTag_
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
edm::InputTag	rechitsCollectionTag_
bool	saveOUTput_
edm::InputTag	src_
DQMStore *	theDMBE
bool	verbose_
	every n events

Detailed Description

Definition at line 86 of file SiPixelMuonHLT.h.

Member Typedef Documentation

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

Definition at line 96 of file SiPixelMuonHLT.h.

Constructor & Destructor Documentation

SiPixelMuonHLT::SiPixelMuonHLT ( const edm::ParameterSet & conf ) [explicit]

Definition at line 24 of file SiPixelMuonHLT.cc.

References clusterCollectionTag_, edm::ParameterSet::getUntrackedParameter(), Histo_init(), l3MuonCollectionTag_, monitorName_, NULL, cmsCodeRules::cppFunctionSkipper::operator, outputFile_, parameters_, rechitsCollectionTag_, saveOUTput_, 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();
   
}

SiPixelMuonHLT::~SiPixelMuonHLT ( )

Definition at line 53 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 79 of file SiPixelMuonHLT.cc.

References SiPixelCluster::charge(), PXFDetId::disk(), PV3DBase< T, PVType, FrameType >::eta(), edm::Event::eventAuxiliary(), eventNo, edm::HandleBase::failedToGet(), GeomDet::geographicalId(), edm::EventSetup::get(), edm::Ref< C, T, F >::get(), edm::Event::getByLabel(), BaseSiTrackerRecHit2DLocalPos::hasPositionAndError(), i, edm::HandleBase::id(), TrackerGeometry::idToDet(), edm::HandleBase::isValid(), l3MuonCollectionTag_, PXBDetId::layer(), BaseSiTrackerRecHit2DLocalPos::localPosition(), 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< T >::product(), PixelGeomDetUnit::specificTopology(), DetId::subdetId(), GeomDet::surface(), Surface::toGlobal(), align::Tracker, SiPixelCluster::x(), PV3DBase< T, PVType, FrameType >::x(), SiPixelCluster::y(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

{
  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.getByLabel("hltSiPixelClusters", 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.getByLabel("hltSiPixelRecHits", 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.getByLabel (l3MuonCollectionTag_, 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 = PXBDetId(detID).layer();
            //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  = PXFDetId(detID).disk();
            if( PXFDetId(detID).side() == 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 = PXBDetId(detID).layer();
                    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  = PXFDetId(detID).disk();
                    if( PXFDetId(detID).disk() == 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 61 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 69 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 252 of file SiPixelMuonHLT.cc.

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

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