CollectionActions/CaloChi2HistoAction/CaloChi2HistoAction.cpp

/*
 * CaloChi2HistoAction.cpp
 *
 *  Created on: 04/ago/2009
 *      Author: Nicola Mori
 */

/*! @file CaloChi2HistoAction.cpp The CaloChi2HistoAction class implementation file. */

#ifndef NO_CALOAXIS

#include "CaloChi2HistoAction.h"

CaloChi2HistoAction::CaloChi2HistoAction(const char *actionName, TString outFileName, CaloAxis *xAxis, CaloAxis *yAxis,
    float min, float max, unsigned int nBins, bool logBinning, TString mode) :
    CollectionAction(actionName), _outFileName(outFileName), _xAxis(xAxis), _yAxis(yAxis), _mode(mode), _xHisto(), _yHisto(), _xyHisto() {

  Double_t *bins = new Double_t[nBins + 1];

  if (!logBinning || (logBinning && min <= 0.)) {

#ifdef DEBUGPAMCUT
    if (logbinning && rigMin <= 0.)
    cout << "Warning: logarithmic binning was chosen for X axis but min <= 0. Using linear binning."
#endif

    float step = (max - min) / nBins;
    for (unsigned int i = 0; i < nBins + 1; i++) {
      bins[i] = min + i * step;
    }

  }
  else {

    double maxExp = log10(max / min);
    for (unsigned int i = 0; i < nBins + 1; i++) {
      bins[i] = min * pow(10., (double) i / ((double) nBins) * maxExp);
    }

  }

  _xHisto.SetName(TString(GetName()) + TString("_x"));
  _xHisto.SetTitle("CALO Chi2 (X view)");
  _xHisto.SetBins(nBins, bins);
  _xHisto.GetXaxis()->SetTitle("Chi2 (X)");
  _xHisto.GetYaxis()->SetTitle("Events");

  _yHisto.SetName(TString(GetName()) + TString("_y"));
  _yHisto.SetBins(nBins, bins);
  _yHisto.SetTitle("CALO Chi2 (Y view)");
  _xHisto.GetXaxis()->SetTitle("Chi2 (Y)");
  _xHisto.GetYaxis()->SetTitle("Events");

  _xyHisto.SetName(TString(GetName()) + TString("_xy"));
  _xyHisto.SetBins(nBins, bins, nBins, bins);
  _xyHisto.SetTitle("CALO Chi2");
  _xyHisto.GetXaxis()->SetTitle("Chi2 (X)");
  _xyHisto.GetYaxis()->SetTitle("Chi2 (Y)");

}

void CaloChi2HistoAction::OnGood(PamLevel2 *event) {

  _xHisto.Fill(_xAxis->GetChi2());
  _yHisto.Fill(_yAxis->GetChi2());
  _xyHisto.Fill(_xAxis->GetChi2(), _yAxis->GetChi2());

}

void CaloChi2HistoAction::Finalize() {

  TFile outFile(_outFileName, _mode);
  outFile.cd();
  _xHisto.Write();
  _yHisto.Write();
  _xyHisto.Write();
  outFile.Close();
}

#endif /* NO_CALOAXIS */
1	/*
2	* CaloChi2HistoAction.cpp
3	*
4	* Created on: 04/ago/2009
5	* Author: Nicola Mori
6	*/
7
8	/! @file CaloChi2HistoAction.cpp The CaloChi2HistoAction class implementation file. /
9
10	#ifndef NO_CALOAXIS
11
12	#include "CaloChi2HistoAction.h"
13
14	CaloChi2HistoAction::CaloChi2HistoAction(const char actionName, TString outFileName, CaloAxis xAxis, CaloAxis *yAxis,
15	float min, float max, unsigned int nBins, bool logBinning, TString mode) :
16	CollectionAction(actionName), _outFileName(outFileName), _xAxis(xAxis), _yAxis(yAxis), _mode(mode), _xHisto(), _yHisto(), _xyHisto() {
17
18	Double_t *bins = new Double_t[nBins + 1];
19
20	if (!logBinning \|\| (logBinning && min <= 0.)) {
21
22	#ifdef DEBUGPAMCUT
23	if (logbinning && rigMin <= 0.)
24	cout << "Warning: logarithmic binning was chosen for X axis but min <= 0. Using linear binning."
25	#endif
26
27	float step = (max - min) / nBins;
28	for (unsigned int i = 0; i < nBins + 1; i++) {
29	bins[i] = min + i * step;
30	}
31
32	}
33	else {
34
35	double maxExp = log10(max / min);
36	for (unsigned int i = 0; i < nBins + 1; i++) {
37	bins[i] = min * pow(10., (double) i / ((double) nBins) * maxExp);
38	}
39
40	}
41
42	_xHisto.SetName(TString(GetName()) + TString("_x"));
43	_xHisto.SetTitle("CALO Chi2 (X view)");
44	_xHisto.SetBins(nBins, bins);
45	_xHisto.GetXaxis()->SetTitle("Chi2 (X)");
46	_xHisto.GetYaxis()->SetTitle("Events");
47
48	_yHisto.SetName(TString(GetName()) + TString("_y"));
49	_yHisto.SetBins(nBins, bins);
50	_yHisto.SetTitle("CALO Chi2 (Y view)");
51	_xHisto.GetXaxis()->SetTitle("Chi2 (Y)");
52	_xHisto.GetYaxis()->SetTitle("Events");
53
54	_xyHisto.SetName(TString(GetName()) + TString("_xy"));
55	_xyHisto.SetBins(nBins, bins, nBins, bins);
56	_xyHisto.SetTitle("CALO Chi2");
57	_xyHisto.GetXaxis()->SetTitle("Chi2 (X)");
58	_xyHisto.GetYaxis()->SetTitle("Chi2 (Y)");
59
60	}
61
62	void CaloChi2HistoAction::OnGood(PamLevel2 *event) {
63
64	_xHisto.Fill(_xAxis->GetChi2());
65	_yHisto.Fill(_yAxis->GetChi2());
66	_xyHisto.Fill(_xAxis->GetChi2(), _yAxis->GetChi2());
67
68	}
69
70	void CaloChi2HistoAction::Finalize() {
71
72	TFile outFile(_outFileName, _mode);
73	outFile.cd();
74	_xHisto.Write();
75	_yHisto.Write();
76	_xyHisto.Write();
77	outFile.Close();
78	}
79
80	#endif /* NO_CALOAXIS */