#include <RPCValidHistograms.h>

Public Types
typedef MonitorElement *	MEP
Public Member Functions
void	bookHistograms (DQMStore *dbe, 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	umOccupancyBarrel_station
MEP	umOccupancyBarrel_wheel
MEP	umOccupancyBarrel_wheel_station
MEP	umOccupancyEndcap_disk
MEP	umOccupancyEndcap_disk_ring
Private Attributes
bool	booked_

Detailed Description

Definition at line 10 of file RPCValidHistograms.h.

Member Typedef Documentation

typedef MonitorElement* RPCValidHistograms::MEP

Definition at line 12 of file RPCValidHistograms.h.

Constructor & Destructor Documentation

RPCValidHistograms::RPCValidHistograms ( ) [inline]

Definition at line 14 of file RPCValidHistograms.h.

References booked_.

  {
    booked_ = false;
  };

Member Function Documentation

void RPCValidHistograms::bookHistograms	(	DQMStore *	dbe,
		const std::string	subDir
	)

Definition at line 5 of file RPCValidHistograms.cc.

{
  if ( !dbe ) 
  {
    edm::LogError("RPCValidHistograms") << "DBE not initialized\n";
    return;
  }
  if ( booked_ )
  {
    edm::LogError("RPCValidHistograms") << "Histogram is already booked\n";
    return;
  }

  const std::string pwd = dbe->pwd();
  dbe->setCurrentFolder(subDir);

  // Book histograms
  dbe->setCurrentFolder(subDir+"/HitProperty");
  clusterSize = dbe->book1D("ClusterSize", "Cluster size;Cluster size", 11, -0.5, 10.5);
  clusterSizeBarrel = dbe->book1D("ClusterSizeBarrel", "Cluster size in Barrel;Cluster size", 11, -0.5, 10.5);
  clusterSizeEndcap = dbe->book1D("ClusterSizeEndcap", "Cluster size in Endcap;Cluster size", 11, -0.5, 10.5);

  avgClusterSize = dbe->book1D("AverageClusterSize", "Average cluster size;Average clsuter size", 11, -0.5, 10.5);
  avgClusterSizeBarrel = dbe->book1D("AverageClusterSizeBarrel", "Average cluster size in Barrel;Average clsuter size", 11, -0.5, 10.5);
  avgClusterSizeEndcap = dbe->book1D("AverageClusterSizeEndcap", "Average cluster size in Endcap;Average clsuter size", 11, -0.5, 10.5);

  nRecHitBarrel = dbe->book1D("NRecHitBarrel", "Number of RPC recHits per event in Barrel;Number of RPC hits", 25, 0, 25);
  nRecHitEndcap = dbe->book1D("NRecHitEndcap", "Number of RPC recHits per event in Endcap;Number of RPC hits", 25, 0, 25);

  nRefHitBarrel = dbe->book1D("NRefHitBarrel", "Number of reference hits per event in Barrel;Number of RPC hits", 25, 0, 25);
  nRefHitEndcap = dbe->book1D("NRefHitEndcap", "Number of reference hits per event in Endcap;Number of RPC hits", 25, 0, 25);

  nMatchHitBarrel = dbe->book1D("nMatchBarrel", "Number of matched reference hits per event in Barrel;Number of RPC hits", 25, 0, 25);
  nMatchHitEndcap = dbe->book1D("nMatchEndcap", "Number of matched reference hits per event in Endcap;Number of RPC hits", 25, 0, 25);

  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
  dbe->setCurrentFolder(subDir+"/Occupancy");
  refHitOccupancyBarrel_wheel   = dbe->book1D("RefHitOccupancyBarrel_wheel"  , "Reference Hit occupancy", 5, -2.5, 2.5);
  refHitOccupancyEndcap_disk    = dbe->book1D("RefHitOccupancyEndcap_disk"   , "Reference Hit occupancy", 9, -4.5, 4.5);
  refHitOccupancyBarrel_station = dbe->book1D("RefHitOccupancyBarrel_station", "Reference Hit occupancy", 4,  0.5, 4.5);

  recHitOccupancyBarrel_wheel   = dbe->book1D("RecHitOccupancyBarrel_wheel"  , "RecHit occupancy", 5, -2.5, 2.5);
  recHitOccupancyEndcap_disk    = dbe->book1D("RecHitOccupancyEndcap_disk"   , "RecHit occupancy", 9, -4.5, 4.5);
  recHitOccupancyBarrel_station = dbe->book1D("RecHitOccupancyBarrel_station", "RecHit occupancy", 4,  0.5, 4.5);

  matchOccupancyBarrel_wheel   = dbe->book1D("MatchOccupancyBarrel_wheel"  , "Matched hit occupancy", 5, -2.5, 2.5);
  matchOccupancyEndcap_disk    = dbe->book1D("MatchOccupancyEndcap_disk"   , "Matched hit occupancy", 9, -4.5, 4.5);
  matchOccupancyBarrel_station = dbe->book1D("MatchOccupancyBarrel_station", "Matched hit occupancy", 4,  0.5, 4.5);

  umOccupancyBarrel_wheel   = dbe->book1D("UmOccupancyBarrel_wheel"  , "Un-matched hit occupancy", 5, -2.5, 2.5);
  umOccupancyEndcap_disk    = dbe->book1D("UmOccupancyEndcap_disk"   , "Un-matched hit occupancy", 9, -4.5, 4.5);
  umOccupancyBarrel_station = dbe->book1D("UmOccupancyBarrel_station", "Un-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);
                               
  umOccupancyBarrel_wheel  ->getTH1()->SetMinimum(0);
  umOccupancyEndcap_disk   ->getTH1()->SetMinimum(0);
  umOccupancyBarrel_station->getTH1()->SetMinimum(0);

  // Occupancy 2D
  refHitOccupancyBarrel_wheel_station = dbe->book2D("RefHitOccupancyBarrel_wheel_station", "Reference hit occupancy", 5, -2.5, 2.5, 4, 0.5, 4.5);
  refHitOccupancyEndcap_disk_ring     = dbe->book2D("RefHitOccupancyEndcap_disk_ring"    , "Reference hit occupancy", 9, -4.5, 4.5, 4, 0.5, 4.5);

  recHitOccupancyBarrel_wheel_station = dbe->book2D("RecHitOccupancyBarrel_wheel_station", "RecHit occupancy", 5, -2.5, 2.5, 4, 0.5, 4.5);
  recHitOccupancyEndcap_disk_ring     = dbe->book2D("RecHitOccupancyEndcap_disk_ring"    , "RecHit occupancy", 9, -4.5, 4.5, 4, 0.5, 4.5);

  matchOccupancyBarrel_wheel_station = dbe->book2D("MatchOccupancyBarrel_wheel_station", "Matched hit occupancy", 5, -2.5, 2.5, 4, 0.5, 4.5);
  matchOccupancyEndcap_disk_ring     = dbe->book2D("MatchOccupancyEndcap_disk_ring"    , "Matched hit occupancy", 9, -4.5, 4.5, 4, 0.5, 4.5);

  umOccupancyBarrel_wheel_station = dbe->book2D("UmOccupancyBarrel_wheel_station", "Un-matched hit occupancy", 5, -2.5, 2.5, 4, 0.5, 4.5);
  umOccupancyEndcap_disk_ring     = dbe->book2D("UmOccupancyEndcap_disk_ring"    , "Un-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);

  umOccupancyBarrel_wheel_station->getTH2F()->SetMinimum(0);
  umOccupancyEndcap_disk_ring    ->getTH2F()->SetMinimum(0);

  // Residuals
  dbe->setCurrentFolder(subDir+"/Residual");
  resBarrel = dbe->book1D("ResBarrel", "Global Residuals for Barrel;Residual [cm]"  , 100, -8, 8);
  resEndcap = dbe->book1D("ResEndcap", "Global Residuals for Endcap;Residual [cm]"   , 100, -8, 8);

  resBarrel->getTH1()->SetMinimum(0);
  resEndcap->getTH1()->SetMinimum(0);

  res_wheel_res   = dbe->book2D("Res_wheel_res"  , "Residuals vs Wheel;;Residual [cm]", 5, -2.5, 2.5, 50, -8, 8);
  res_disk_res    = dbe->book2D("Res_disk_res"   , "Residuals vs Disk;;Residual [cm]", 9, -4.5, 4.5, 50, -8, 8);
  res_station_res = dbe->book2D("Res_station_res", "Redisuals vs Station;;Residual [cm]", 4, 0.5, 4.5, 50, -8, 8);
  res_ring_res    = dbe->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 = dbe->book1D("PullBarrel", "Global Pull for Barrel;Pull", 100, -3, 3);
  pullEndcap = dbe->book1D("PullEndcap", "Global Pull for Endcap;Pull", 100, -3, 3);

  pullBarrel->getTH1()->SetMinimum(0);
  pullEndcap->getTH1()->SetMinimum(0);

  pull_wheel_pull   = dbe->book2D("Pull_wheel_pull"  , "Pull vs Wheel;;Pull"  , 5, -2.5, 2.5, 50, -3, 3);
  pull_disk_pull    = dbe->book2D("Pull_disk_pull"   , "Pull vs Disk;;Pull"   , 9, -4.5, 4.5, 50, -3, 3);
  pull_station_pull = dbe->book2D("Pull_station_pull", "Pull vs Station;;Pull", 4,  0.5, 4.5, 50, -3, 3);
  pull_ring_pull    = dbe->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");
  umOccupancyBarrel_wheel_station    ->getTH2F()->SetOption("COLZ");
  umOccupancyEndcap_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);
  umOccupancyBarrel_wheel_station    ->getTH2F()->SetContour(10);
  umOccupancyEndcap_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);
  umOccupancyBarrel_wheel_station    ->getTH2F()->SetStats(0);
  umOccupancyEndcap_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);
    umOccupancyBarrel_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);
    umOccupancyBarrel_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);
    umOccupancyEndcap_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);
    umOccupancyEndcap_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);
    umOccupancyBarrel_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);
    umOccupancyBarrel_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);
    umOccupancyEndcap_disk_ring    ->getTH2F()->GetYaxis()->SetBinLabel(i, binLabel);
  }

  dbe->setCurrentFolder(pwd);
  booked_ = true;
}