cyclops/FillCalHist.C at master · tkhemmick/cyclops · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
////////////////////////////////////
// FillCalHist.C                    //
////////////////////////////////////

//  This makes a histogram that can show the calibrated hits.

#include "FillCalHist.h"
#include "TF1.h"
#include "TH1D.h"
#include "TH2D.h"
#include "APad.h"
#include "AZig.h"
#include "TSpectrum.h"

#include "cyclops.h"

#include <iostream>
#include <cmath>

TH2D *CalHist=0;
TH1D *CalDist[NumberOfPads];

TH2D *CalZigHist=0;
TH1D *CalZigDist[NumberOfZigs];

TH1D *slice=0;

using namespace std;

void FillCalHist()
{

  int nTime = APad::Cal[0].size();
  double left = -0.5;
  double right = (double)nTime - 0.5;

  if (!CalHist)
    {
      CalHist = new TH2D("CalHist","CalHist",NumberOfPads,-0.5,(double)NumberOfPads-0.5,nTime,left,right);
      cout << "Making the histogram" << endl;

      char name[500];
      for (int i=0; i<NumberOfPads; i++)
	{
	  sprintf(name,"CalDist%03d",i);
	  CalDist[i] = new TH1D(name,name,3301,-300.5,3000.5);
	}

    }
  CalHist->Reset();

  for (int i=0; i<NumberOfPads; i++)
    {
      for (int j=0; j<APad::Cal[i].size(); j++)
	{
	  CalHist->SetBinContent(i+1,j+1,APad::Cal[i][j]);
	  CalDist[i]->Fill(APad::Cal[i][j]);
	}
    }

}


void FillCalZigHist()
{

  int nTime = AZig::Cal[0].size();
  double left = -0.5;
  double right = (double)nTime - 0.5;

  if (!CalZigHist)
    {
      CalZigHist = new TH2D("CalZigHist","CalZigHist",NumberOfZigs,-0.5,(double)NumberOfZigs-0.5,nTime,left,right);
      cout << "Making the histogram" << endl;

      char name[500];
      for (int i=0; i<NumberOfZigs; i++)
	{
	  sprintf(name,"CalZigDist%03d",i);
	  CalZigDist[i] = new TH1D(name,name,3301,-300.5,3000.5);
	}

    }
  CalZigHist->Reset();

  for (int i=0; i<NumberOfZigs; i++)
    {
      for (int j=0; j<AZig::Cal[i].size(); j++)
	{
	  CalZigHist->SetBinContent(i+1,j+1,AZig::Cal[i][j]);
	  CalZigDist[i]->Fill(AZig::Cal[i][j]);
	}
    }

}