EcalDrawUtils.h

Go to the documentation of this file.

#include "TH2F.h"
#include "TLine.h"
 
void DrawEB(TH2F* ebmap, float min, float max) {
  ebmap->SetXTitle("i#phi");
  ebmap->SetYTitle("i#eta");
  ebmap->GetXaxis()->SetNdivisions(-418, kFALSE);
  ebmap->GetYaxis()->SetNdivisions(-1702, kFALSE);
  ebmap->GetXaxis()->SetLabelSize(0.03);
  ebmap->GetYaxis()->SetLabelSize(0.03);
  ebmap->GetXaxis()->SetTickLength(0.01);
  ebmap->GetYaxis()->SetTickLength(0.01);
  ebmap->SetMaximum(max);
  ebmap->SetMinimum(min);
  ebmap->Draw("colz");
  TLine* l = new TLine;
  l->SetLineWidth(1);
  for (int i = 0; i < 17; i++) {
    Double_t x = 20. + (i * 20);
    l = new TLine(x, -85., x, 85.);
    l->Draw();
  }
  l = new TLine(0., 85., 360., 85.);
  l->Draw();
  l = new TLine(0., 0., 360., 0.);
  l->Draw();
}  //   DrawEB method
 
void DrawEE(TH2F* endc, float min, float max) {
  int ixSectorsEE[202] = {
      62, 62, 61, 61, 60, 60, 59, 59, 58, 58, 56, 56, 46, 46, 44, 44, 43, 43, 42, 42,  41,  41,  40,  40, 41, 41,
      42, 42, 43, 43, 44, 44, 46, 46, 56, 56, 58, 58, 59, 59, 60, 60, 61, 61, 62, 62,  0,   101, 101, 98, 98, 96,
      96, 93, 93, 88, 88, 86, 86, 81, 81, 76, 76, 66, 66, 61, 61, 41, 41, 36, 36, 26,  26,  21,  21,  16, 16, 14,
      14, 9,  9,  6,  6,  4,  4,  1,  1,  4,  4,  6,  6,  9,  9,  14, 14, 16, 16, 21,  21,  26,  26,  36, 36, 41,
      41, 61, 61, 66, 66, 76, 76, 81, 81, 86, 86, 88, 88, 93, 93, 96, 96, 98, 98, 101, 101, 0,   62,  66, 66, 71,
      71, 81, 81, 91, 91, 93, 0,  62, 66, 66, 91, 91, 98, 0,  58, 61, 61, 66, 66, 71,  71,  76,  76,  81, 81, 0,
      51, 51, 0,  44, 41, 41, 36, 36, 31, 31, 26, 26, 21, 21, 0,  40, 36, 36, 11, 11,  4,   0,   40,  36, 36, 31,
      31, 21, 21, 11, 11, 9,  0,  46, 46, 41, 41, 36, 36, 0,  56, 56, 61, 61, 66, 66};
 
  int iySectorsEE[202] = {51, 56, 56, 58, 58, 59, 59, 60, 60, 61, 61, 62, 62, 61, 61, 60, 60, 59, 59, 58,  58,  56,  56,
                          46, 46, 44, 44, 43, 43, 42, 42, 41, 41, 40, 40, 41, 41, 42, 42, 43, 43, 44, 44,  46,  46,  51,
                          0,  51, 61, 61, 66, 66, 76, 76, 81, 81, 86, 86, 88, 88, 93, 93, 96, 96, 98, 98,  101, 101, 98,
                          98, 96, 96, 93, 93, 88, 88, 86, 86, 81, 81, 76, 76, 66, 66, 61, 61, 41, 41, 36,  36,  26,  26,
                          21, 21, 16, 16, 14, 14, 9,  9,  6,  6,  4,  4,  1,  1,  4,  4,  6,  6,  9,  9,   14,  14,  16,
                          16, 21, 21, 26, 26, 36, 36, 41, 41, 51, 0,  46, 46, 41, 41, 36, 36, 31, 31, 26,  26,  0,   51,
                          51, 56, 56, 61, 61, 0,  61, 61, 66, 66, 71, 71, 76, 76, 86, 86, 88, 0,  62, 101, 0,   61,  61,
                          66, 66, 71, 71, 76, 76, 86, 86, 88, 0,  51, 51, 56, 56, 61, 61, 0,  46, 46, 41,  41,  36,  36,
                          31, 31, 26, 26, 0,  40, 31, 31, 16, 16, 6,  0,  40, 31, 31, 16, 16, 6};
 
  TLine* l = new TLine;
  l->SetLineWidth(1);
 
  endc->SetXTitle("ix");
  endc->SetYTitle("iy");
  endc->SetMaximum(max);
  endc->SetMinimum(min);
  endc->Draw("colz1");
  for (int i = 0; i < 201; i = i + 1) {
    if ((ixSectorsEE[i] != 0 || iySectorsEE[i] != 0) && (ixSectorsEE[i + 1] != 0 || iySectorsEE[i + 1] != 0)) {
      l->DrawLine(ixSectorsEE[i], iySectorsEE[i], ixSectorsEE[i + 1], iySectorsEE[i + 1]);
    }
  }
}  //   DrawEE method
 
void DrawEE_Tower(TH2F* endc, TLine* l, double minScale, double maxScale) {
  endc->SetStats(false);
  endc->SetMinimum(minScale);
  endc->SetMaximum(maxScale);
  endc->Draw("colz");
 
  int ixSectorsEE[136] = {8,  14, 14, 17, 17, 18, 18, 19, 19, 20, 20, 21, 21, 20, 20, 19, 19, 18, 18, 17, 17, 14, 14,
                          8,  8,  5,  5,  4,  4,  3,  3,  2,  2,  1,  1,  2,  2,  3,  3,  4,  4,  5,  5,  8,  8,  8,
                          9,  9,  10, 10, 12, 12, 13, 13, 12, 12, 10, 10, 9,  9,  10, 10, 0,  11, 11, 0,  10, 9,  9,
                          8,  8,  7,  7,  6,  6,  5,  5,  0,  12, 13, 13, 14, 14, 15, 15, 16, 16, 17, 17, 0,  9,  8,
                          8,  3,  3,  1,  0,  13, 14, 14, 19, 19, 21, 0,  9,  8,  8,  7,  7,  5,  5,  3,  3,  2,  0,
                          13, 14, 14, 15, 15, 17, 17, 19, 19, 20, 0,  14, 14, 13, 13, 12, 12, 0};
 
  int iySectorsEE[136] = {1,  1,  2,  2,  3,  3,  4,  4,  5,  5,  8,  8,  14, 14, 17, 17, 18, 18, 19, 19, 20, 20, 21,
                          21, 20, 20, 19, 19, 18, 18, 17, 17, 14, 14, 8,  8,  5,  5,  4,  4,  3,  3,  2,  2,  1,  4,
                          4,  7,  7,  9,  9,  10, 10, 12, 12, 13, 13, 12, 12, 10, 10, 9,  0,  13, 21, 0,  13, 13, 14,
                          14, 15, 15, 16, 16, 18, 18, 19, 0,  13, 13, 14, 14, 15, 15, 16, 16, 18, 18, 19, 0,  11, 11,
                          12, 12, 13, 13, 0,  11, 11, 12, 12, 13, 13, 0,  10, 10, 9,  9,  8,  8,  7,  7,  6,  6,  0,
                          10, 10, 9,  9,  8,  8,  7,  7,  6,  6,  0,  2,  4,  4,  7,  7,  9,  0};
 
  for (int i = 0; i < 136; i = i + 1)
    if ((ixSectorsEE[i] != 0 || iySectorsEE[i] != 0) && (ixSectorsEE[i + 1] != 0 || iySectorsEE[i + 1] != 0))
      l->DrawLine(ixSectorsEE[i], iySectorsEE[i], ixSectorsEE[i + 1], iySectorsEE[i + 1]);
 
}  //draw EE in case of a tower
 
void drawTable(int nbRows, int nbColumns) {
  TLine* l = new TLine;
  l->SetLineWidth(1);
  for (int i = 1; i < nbRows; i++) {
    double y = (double)i;
    l = new TLine(0., y, nbColumns, y);
    l->Draw();
  }
 
  for (int i = 1; i < nbColumns; i++) {
    double x = (double)i;
    double y = (double)nbRows;
    l = new TLine(x, 0., x, y);
    l->Draw();
  }
}