#include <RPCValidHistograms.h>

Public Types
typedef MonitorElement *	MEP

Public Member Functions
void	bookHistograms (DQMStore::IBooker &booker, const std::string &subDir)

	RPCValidHistograms ()

Public Attributes
MEP	avgClusterSize

MEP	avgClusterSizeBarrel

MEP	avgClusterSizeEndcap

MEP	clusterSize

MEP	clusterSizeBarrel

MEP	clusterSizeEndcap

MEP	matchOccupancyBarrel_station

MEP	matchOccupancyBarrel_wheel

MEP	matchOccupancyBarrel_wheel_station

MEP	matchOccupancyEndcap_disk

MEP	matchOccupancyEndcap_disk_ring

MEP	nMatchHitBarrel

MEP	nMatchHitEndcap

MEP	nRecHitBarrel

MEP	nRecHitEndcap

MEP	nRefHitBarrel

MEP	nRefHitEndcap

MEP	pull_disk_pull

MEP	pull_ring_pull

MEP	pull_station_pull

MEP	pull_wheel_pull

MEP	pullBarrel

MEP	pullEndcap

MEP	recHitOccupancyBarrel_station

MEP	recHitOccupancyBarrel_wheel

MEP	recHitOccupancyBarrel_wheel_station

MEP	recHitOccupancyEndcap_disk

MEP	recHitOccupancyEndcap_disk_ring

MEP	refHitOccupancyBarrel_station

MEP	refHitOccupancyBarrel_wheel

MEP	refHitOccupancyBarrel_wheel_station

MEP	refHitOccupancyEndcap_disk

MEP	refHitOccupancyEndcap_disk_ring

MEP	res_disk_res

MEP	res_ring_res

MEP	res_station_res

MEP	res_wheel_res

MEP	resBarrel

MEP	resEndcap

MEP	timeBarrel

MEP	timeCRPC

MEP	timeEndcap

MEP	timeIRPC

Private Attributes
bool	booked_

Detailed Description

Definition at line 12 of file RPCValidHistograms.h.

Member Typedef Documentation

typedef MonitorElement* RPCValidHistograms::MEP

Definition at line 14 of file RPCValidHistograms.h.

Constructor & Destructor Documentation

RPCValidHistograms::RPCValidHistograms ( )

inline

Definition at line 16 of file RPCValidHistograms.h.

References booked_, bookHistograms(), AlCaHLTBitMon_QueryRunRegistry::string, and muonValidation_cff::subDir.

   {
     booked_ = false;
   };

Member Function Documentation

void RPCValidHistograms::bookHistograms	(	DQMStore::IBooker &	booker,
		const std::string &	subDir
	)

Definition at line 5 of file RPCValidHistograms.cc.

References avgClusterSize, avgClusterSizeBarrel, avgClusterSizeEndcap, DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), booked_, clusterSize, clusterSizeBarrel, clusterSizeEndcap, MonitorElement::getTH1(), MonitorElement::getTH1F(), MonitorElement::getTH2F(), mps_fire::i, matchOccupancyBarrel_station, matchOccupancyBarrel_wheel, matchOccupancyBarrel_wheel_station, matchOccupancyEndcap_disk, matchOccupancyEndcap_disk_ring, nMatchHitBarrel, nMatchHitEndcap, nRecHitBarrel, nRecHitEndcap, nRefHitBarrel, nRefHitEndcap, pull_disk_pull, pull_ring_pull, pull_station_pull, pull_wheel_pull, pullBarrel, pullEndcap, alignCSCRings::pwd, DQMStore::IBooker::pwd(), recHitOccupancyBarrel_station, recHitOccupancyBarrel_wheel, recHitOccupancyBarrel_wheel_station, recHitOccupancyEndcap_disk, recHitOccupancyEndcap_disk_ring, refHitOccupancyBarrel_station, refHitOccupancyBarrel_wheel, refHitOccupancyBarrel_wheel_station, refHitOccupancyEndcap_disk, refHitOccupancyEndcap_disk_ring, res_disk_res, res_ring_res, res_station_res, res_wheel_res, resBarrel, resEndcap, DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, timeBarrel, timeCRPC, timeEndcap, and timeIRPC.

Referenced by RPCValidHistograms().

 {
   if ( booked_ )
   {
     edm::LogError("RPCValidHistograms") << "Histogram is already booked\n";
     return;
   }
 
   const std::string pwd = booker.pwd();
   booker.setCurrentFolder(subDir);
 
   // Book histograms
   booker.setCurrentFolder(subDir+"/HitProperty");
   clusterSize = booker.book1D("ClusterSize", "Cluster size;Cluster size", 11, -0.5, 10.5);
   clusterSizeBarrel = booker.book1D("ClusterSizeBarrel", "Cluster size in Barrel;Cluster size", 11, -0.5, 10.5);
   clusterSizeEndcap = booker.book1D("ClusterSizeEndcap", "Cluster size in Endcap;Cluster size", 11, -0.5, 10.5);
 
   avgClusterSize = booker.book1D("AverageClusterSize", "Average cluster size;Average clsuter size", 11, -0.5, 10.5);
   avgClusterSizeBarrel = booker.book1D("AverageClusterSizeBarrel", "Average cluster size in Barrel;Average clsuter size", 11, -0.5, 10.5);
   avgClusterSizeEndcap = booker.book1D("AverageClusterSizeEndcap", "Average cluster size in Endcap;Average clsuter size", 11, -0.5, 10.5);
 
   nRecHitBarrel = booker.book1D("NRecHitBarrel", "Number of RPC recHits per event in Barrel;Number of RPC hits", 25, 0, 25);
   nRecHitEndcap = booker.book1D("NRecHitEndcap", "Number of RPC recHits per event in Endcap;Number of RPC hits", 25, 0, 25);
 
   nRefHitBarrel = booker.book1D("NRefHitBarrel", "Number of reference hits per event in Barrel;Number of RPC hits", 25, 0, 25);
   nRefHitEndcap = booker.book1D("NRefHitEndcap", "Number of reference hits per event in Endcap;Number of RPC hits", 25, 0, 25);
 
   nMatchHitBarrel = booker.book1D("nMatchBarrel", "Number of matched reference hits per event in Barrel;Number of RPC hits", 25, 0, 25);
   nMatchHitEndcap = booker.book1D("nMatchEndcap", "Number of matched reference hits per event in Endcap;Number of RPC hits", 25, 0, 25);
 
   timeBarrel = booker.book1D("RecHitTimeBarrel", "RecHit time in Barrel;Time (ns)", 100, -12.5, 12.5);
   timeEndcap = booker.book1D("RecHitTimeEndcap", "RecHit time in Endcap;Time (ns)", 100, -12.5, 12.5);
   timeIRPC = booker.book1D("RecHitTimeIRPC", "RecHit time of iRPC;Time (ns)", 100, -12.5, 12.5);
   timeCRPC = booker.book1D("RecHitTimeCPRC", "RecHit time in current RPC;Time (ns)", 100, -12.5, 12.5);
 
   clusterSize->getTH1()->SetMinimum(0);
   clusterSizeBarrel->getTH1()->SetMinimum(0);
   clusterSizeEndcap->getTH1()->SetMinimum(0);
                       
   avgClusterSize->getTH1()->SetMinimum(0);
   avgClusterSizeBarrel->getTH1()->SetMinimum(0);
   avgClusterSizeEndcap->getTH1()->SetMinimum(0);
                       
   nRecHitBarrel->getTH1()->SetMinimum(0);
   nRecHitEndcap->getTH1()->SetMinimum(0);
                
   nRefHitBarrel->getTH1()->SetMinimum(0);
   nRefHitEndcap->getTH1()->SetMinimum(0);
                       
   nMatchHitBarrel->getTH1()->SetMinimum(0);
   nMatchHitEndcap->getTH1()->SetMinimum(0);
 
   // Occupancy 1D
   booker.setCurrentFolder(subDir+"/Occupancy");
   refHitOccupancyBarrel_wheel   = booker.book1D("RefHitOccupancyBarrel_wheel"  , "Reference Hit occupancy", 5, -2.5, 2.5);
   refHitOccupancyEndcap_disk    = booker.book1D("RefHitOccupancyEndcap_disk"   , "Reference Hit occupancy", 9, -4.5, 4.5);
   refHitOccupancyBarrel_station = booker.book1D("RefHitOccupancyBarrel_station", "Reference Hit occupancy", 4,  0.5, 4.5);
 
   recHitOccupancyBarrel_wheel   = booker.book1D("RecHitOccupancyBarrel_wheel"  , "RecHit occupancy", 5, -2.5, 2.5);
   recHitOccupancyEndcap_disk    = booker.book1D("RecHitOccupancyEndcap_disk"   , "RecHit occupancy", 9, -4.5, 4.5);
   recHitOccupancyBarrel_station = booker.book1D("RecHitOccupancyBarrel_station", "RecHit occupancy", 4,  0.5, 4.5);
 
   matchOccupancyBarrel_wheel   = booker.book1D("MatchOccupancyBarrel_wheel"  , "Matched hit occupancy", 5, -2.5, 2.5);
   matchOccupancyEndcap_disk    = booker.book1D("MatchOccupancyEndcap_disk"   , "Matched hit occupancy", 9, -4.5, 4.5);
   matchOccupancyBarrel_station = booker.book1D("MatchOccupancyBarrel_station", "Matched hit occupancy", 4,  0.5, 4.5);
 
   refHitOccupancyBarrel_wheel  ->getTH1()->SetMinimum(0);
   refHitOccupancyEndcap_disk   ->getTH1()->SetMinimum(0);
   refHitOccupancyBarrel_station->getTH1()->SetMinimum(0);
                                
   recHitOccupancyBarrel_wheel  ->getTH1()->SetMinimum(0);
   recHitOccupancyEndcap_disk   ->getTH1()->SetMinimum(0);
   recHitOccupancyBarrel_station->getTH1()->SetMinimum(0);
                                
   matchOccupancyBarrel_wheel  ->getTH1()->SetMinimum(0);
   matchOccupancyEndcap_disk   ->getTH1()->SetMinimum(0);
   matchOccupancyBarrel_station->getTH1()->SetMinimum(0);
                                
   // Occupancy 2D
   refHitOccupancyBarrel_wheel_station = booker.book2D("RefHitOccupancyBarrel_wheel_station", "Reference hit occupancy", 5, -2.5, 2.5, 4, 0.5, 4.5);
   refHitOccupancyEndcap_disk_ring     = booker.book2D("RefHitOccupancyEndcap_disk_ring"    , "Reference hit occupancy", 9, -4.5, 4.5, 4, 0.5, 4.5);
 
   recHitOccupancyBarrel_wheel_station = booker.book2D("RecHitOccupancyBarrel_wheel_station", "RecHit occupancy", 5, -2.5, 2.5, 4, 0.5, 4.5);
   recHitOccupancyEndcap_disk_ring     = booker.book2D("RecHitOccupancyEndcap_disk_ring"    , "RecHit occupancy", 9, -4.5, 4.5, 4, 0.5, 4.5);
 
   matchOccupancyBarrel_wheel_station = booker.book2D("MatchOccupancyBarrel_wheel_station", "Matched hit occupancy", 5, -2.5, 2.5, 4, 0.5, 4.5);
   matchOccupancyEndcap_disk_ring     = booker.book2D("MatchOccupancyEndcap_disk_ring"    , "Matched hit occupancy", 9, -4.5, 4.5, 4, 0.5, 4.5);
 
   refHitOccupancyBarrel_wheel_station->getTH2F()->SetMinimum(0);
   refHitOccupancyEndcap_disk_ring    ->getTH2F()->SetMinimum(0);
 
   recHitOccupancyBarrel_wheel_station->getTH2F()->SetMinimum(0);
   recHitOccupancyEndcap_disk_ring    ->getTH2F()->SetMinimum(0);
 
   matchOccupancyBarrel_wheel_station->getTH2F()->SetMinimum(0);
   matchOccupancyEndcap_disk_ring    ->getTH2F()->SetMinimum(0);
 
   // Residuals
   booker.setCurrentFolder(subDir+"/Residual");
   resBarrel = booker.book1D("ResBarrel", "Global Residuals for Barrel;Residual [cm]"  , 100, -8, 8);
   resEndcap = booker.book1D("ResEndcap", "Global Residuals for Endcap;Residual [cm]"   , 100, -8, 8);
 
   resBarrel->getTH1()->SetMinimum(0);
   resEndcap->getTH1()->SetMinimum(0);
 
   res_wheel_res   = booker.book2D("Res_wheel_res"  , "Residuals vs Wheel;;Residual [cm]", 5, -2.5, 2.5, 50, -8, 8);
   res_disk_res    = booker.book2D("Res_disk_res"   , "Residuals vs Disk;;Residual [cm]", 9, -4.5, 4.5, 50, -8, 8);
   res_station_res = booker.book2D("Res_station_res", "Redisuals vs Station;;Residual [cm]", 4, 0.5, 4.5, 50, -8, 8);
   res_ring_res    = booker.book2D("Res_ring_res"   , "Redisuals vs Ring;;Residual [cm]", 4, 0.5, 4.5, 50, -8, 8);
 
   res_wheel_res  ->getTH2F()->SetMinimum(0);
   res_disk_res   ->getTH2F()->SetMinimum(0);
   res_station_res->getTH2F()->SetMinimum(0);
   res_ring_res   ->getTH2F()->SetMinimum(0);
 
   // Pulls
   pullBarrel = booker.book1D("PullBarrel", "Global Pull for Barrel;Pull", 100, -3, 3);
   pullEndcap = booker.book1D("PullEndcap", "Global Pull for Endcap;Pull", 100, -3, 3);
 
   pullBarrel->getTH1()->SetMinimum(0);
   pullEndcap->getTH1()->SetMinimum(0);
 
   pull_wheel_pull   = booker.book2D("Pull_wheel_pull"  , "Pull vs Wheel;;Pull"  , 5, -2.5, 2.5, 50, -3, 3);
   pull_disk_pull    = booker.book2D("Pull_disk_pull"   , "Pull vs Disk;;Pull"   , 9, -4.5, 4.5, 50, -3, 3);
   pull_station_pull = booker.book2D("Pull_station_pull", "Pull vs Station;;Pull", 4,  0.5, 4.5, 50, -3, 3);
   pull_ring_pull    = booker.book2D("Pull_ring_pull"   , "Pull vs Ring;;Pull"   , 4,  0.5, 4.5, 50, -3, 3);
 
   pull_wheel_pull  ->getTH2F()->SetMinimum(0);
   pull_disk_pull   ->getTH2F()->SetMinimum(0);
   pull_station_pull->getTH2F()->SetMinimum(0);
   pull_ring_pull   ->getTH2F()->SetMinimum(0);
 
   // Set plot options
   refHitOccupancyBarrel_wheel_station->getTH2F()->SetOption("COLZ");
   refHitOccupancyEndcap_disk_ring    ->getTH2F()->SetOption("COLZ");
   recHitOccupancyBarrel_wheel_station->getTH2F()->SetOption("COLZ");
   recHitOccupancyEndcap_disk_ring    ->getTH2F()->SetOption("COLZ");
   matchOccupancyBarrel_wheel_station ->getTH2F()->SetOption("COLZ");
   matchOccupancyEndcap_disk_ring     ->getTH2F()->SetOption("COLZ");
 
   res_wheel_res  ->getTH2F()->SetOption("COLZ");
   res_disk_res   ->getTH2F()->SetOption("COLZ");
   res_station_res->getTH2F()->SetOption("COLZ");
   res_ring_res   ->getTH2F()->SetOption("COLZ");
 
   pull_wheel_pull  ->getTH2F()->SetOption("COLZ");
   pull_disk_pull   ->getTH2F()->SetOption("COLZ");
   pull_station_pull->getTH2F()->SetOption("COLZ");
   pull_ring_pull   ->getTH2F()->SetOption("COLZ");
 
   refHitOccupancyBarrel_wheel_station->getTH2F()->SetContour(10);
   refHitOccupancyEndcap_disk_ring    ->getTH2F()->SetContour(10);
   recHitOccupancyBarrel_wheel_station->getTH2F()->SetContour(10);
   recHitOccupancyEndcap_disk_ring    ->getTH2F()->SetContour(10);
   matchOccupancyBarrel_wheel_station ->getTH2F()->SetContour(10);
   matchOccupancyEndcap_disk_ring     ->getTH2F()->SetContour(10);
 
   res_wheel_res  ->getTH2F()->SetContour(10);
   res_disk_res   ->getTH2F()->SetContour(10);
   res_station_res->getTH2F()->SetContour(10);
   res_ring_res   ->getTH2F()->SetContour(10);
 
   pull_wheel_pull  ->getTH2F()->SetContour(10);
   pull_disk_pull   ->getTH2F()->SetContour(10);
   pull_station_pull->getTH2F()->SetContour(10);
   pull_ring_pull   ->getTH2F()->SetContour(10);
  
   refHitOccupancyBarrel_wheel_station->getTH2F()->SetStats(0);
   refHitOccupancyEndcap_disk_ring    ->getTH2F()->SetStats(0);
   recHitOccupancyBarrel_wheel_station->getTH2F()->SetStats(0);
   recHitOccupancyEndcap_disk_ring    ->getTH2F()->SetStats(0);
   matchOccupancyBarrel_wheel_station ->getTH2F()->SetStats(0);
   matchOccupancyEndcap_disk_ring     ->getTH2F()->SetStats(0);
 
   res_wheel_res  ->getTH2F()->SetStats(0);
   res_disk_res   ->getTH2F()->SetStats(0);
   res_station_res->getTH2F()->SetStats(0);
   res_ring_res   ->getTH2F()->SetStats(0);
 
   pull_wheel_pull  ->getTH2F()->SetStats(0);
   pull_disk_pull   ->getTH2F()->SetStats(0);
   pull_station_pull->getTH2F()->SetStats(0);
   pull_ring_pull   ->getTH2F()->SetStats(0);
  
   // Set bin labels
   for ( int i=1; i<=5; ++i )
   {
     TString binLabel = Form("Wheel %d", i-3);
 
     refHitOccupancyBarrel_wheel->getTH1F()->GetXaxis()->SetBinLabel(i, binLabel);
     recHitOccupancyBarrel_wheel->getTH1F()->GetXaxis()->SetBinLabel(i, binLabel);
     matchOccupancyBarrel_wheel ->getTH1F()->GetXaxis()->SetBinLabel(i, binLabel);
 
     refHitOccupancyBarrel_wheel_station->getTH2F()->GetXaxis()->SetBinLabel(i, binLabel);
     recHitOccupancyBarrel_wheel_station->getTH2F()->GetXaxis()->SetBinLabel(i, binLabel);
     matchOccupancyBarrel_wheel_station ->getTH2F()->GetXaxis()->SetBinLabel(i, binLabel);
 
     res_wheel_res  ->getTH2F()->GetXaxis()->SetBinLabel(i, binLabel);
     pull_wheel_pull->getTH2F()->GetXaxis()->SetBinLabel(i, binLabel);
   }
 
   for ( int i=1; i<=9; ++i )
   {
     TString binLabel = Form("Disk %d", i-5);
 
     refHitOccupancyEndcap_disk->getTH1F()->GetXaxis()->SetBinLabel(i, binLabel);
     recHitOccupancyEndcap_disk->getTH1F()->GetXaxis()->SetBinLabel(i, binLabel);
     matchOccupancyEndcap_disk ->getTH1F()->GetXaxis()->SetBinLabel(i, binLabel);
 
     refHitOccupancyEndcap_disk_ring->getTH2F()->GetXaxis()->SetBinLabel(i, binLabel);
     recHitOccupancyEndcap_disk_ring->getTH2F()->GetXaxis()->SetBinLabel(i, binLabel);
     matchOccupancyEndcap_disk_ring ->getTH2F()->GetXaxis()->SetBinLabel(i, binLabel);
 
     res_disk_res  ->getTH2F()->GetXaxis()->SetBinLabel(i, binLabel);
     pull_disk_pull->getTH2F()->GetXaxis()->SetBinLabel(i, binLabel);
   }
 
   for ( int i=1; i<=4; ++i )
   {
     TString binLabel = Form("Station %d", i);
 
     refHitOccupancyBarrel_station->getTH1F()->GetXaxis()->SetBinLabel(i, binLabel);
     recHitOccupancyBarrel_station->getTH1F()->GetXaxis()->SetBinLabel(i, binLabel);
     matchOccupancyBarrel_station ->getTH1F()->GetXaxis()->SetBinLabel(i, binLabel);
 
     refHitOccupancyBarrel_wheel_station->getTH2F()->GetYaxis()->SetBinLabel(i, binLabel);
     recHitOccupancyBarrel_wheel_station->getTH2F()->GetYaxis()->SetBinLabel(i, binLabel);
     matchOccupancyBarrel_wheel_station ->getTH2F()->GetYaxis()->SetBinLabel(i, binLabel);
 
     res_station_res  ->getTH2F()->GetXaxis()->SetBinLabel(i, binLabel);
     pull_station_pull->getTH2F()->GetXaxis()->SetBinLabel(i, binLabel);
   }
 
   for ( int i=1; i<=4; ++i )
   {
     TString binLabel = Form("Ring %d", i);
 
     refHitOccupancyEndcap_disk_ring->getTH2F()->GetYaxis()->SetBinLabel(i, binLabel);
     recHitOccupancyEndcap_disk_ring->getTH2F()->GetYaxis()->SetBinLabel(i, binLabel);
     matchOccupancyEndcap_disk_ring ->getTH2F()->GetYaxis()->SetBinLabel(i, binLabel);
   }
 
   booked_ = true;
 
   booker.setCurrentFolder(pwd);
 }