CollectionActions/CaloDedxVsBetaHistoAction/CaloDedxVsBetaHistoAction.cpp

/*
 * CaloDedxVsBetaHistoAction.cpp
 *
 *  Created on: 03/ago/2009
 *      Author: Nicola Mori
 */

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

#include "CaloDedxVsBetaHistoAction.h"

CaloDedxVsBetaHistoAction::CaloDedxVsBetaHistoAction(const char *actionName, TString outFileName, unsigned int plane,
    unsigned int view, TString mode, float minDedx, float maxDedx, unsigned int nBinsDedx, float minBeta,
    float maxBeta, unsigned int nBinsBeta) :
  CollectionAction(actionName), _outFileName(outFileName), _plane(plane), _view(view), _histo(), _mode(mode) {

  _histo.SetName(GetName());
  TString title = "CALO dE/dx (plane ";
  title += _plane;
  title += ", ";
  if (_view == 0)
    title += "X";
  if (_view == 1)
    title += "Y";
  title += " view) Vs Beta";
  _histo.SetTitle(title);
  _histo.GetXaxis()->SetTitle("Beta");
  _histo.GetYaxis()->SetTitle("dE/dx (MIP)");

  _histo.SetBins(nBinsBeta, minBeta, maxBeta, nBinsDedx, minDedx, maxDedx);
}

void CaloDedxVsBetaHistoAction::OnGood(PamLevel2 *event) {

  _histo.Fill(event->GetTrack(0)->GetToFTrack()->beta[12], event->GetCaloLevel1()->qtotpl(_view, _plane));
  //cout << event->GetTrack(0)->GetToFTrack()->beta[12] << "  "<< event->GetCaloLevel1()->qtotpl(_view, _plane) << endl;
}

void CaloDedxVsBetaHistoAction::Finalize() {

  TFile outFile(_outFileName, _mode);
  outFile.cd();
  _histo.Write();
  outFile.Close();
}
1	pam-fi	1.1	/*
2			* CaloDedxVsBetaHistoAction.cpp
3			*
4			* Created on: 03/ago/2009
5			* Author: Nicola Mori
6			*/
7
8			/! @file CaloDedxVsBetaHistoAction.cpp The CaloDedxVsBetaHistoAction class implementation file. /
9
10			#include "CaloDedxVsBetaHistoAction.h"
11
12			CaloDedxVsBetaHistoAction::CaloDedxVsBetaHistoAction(const char *actionName, TString outFileName, unsigned int plane,
13			unsigned int view, TString mode, float minDedx, float maxDedx, unsigned int nBinsDedx, float minBeta,
14			float maxBeta, unsigned int nBinsBeta) :
15			CollectionAction(actionName), _outFileName(outFileName), _plane(plane), _view(view), _histo(), _mode(mode) {
16
17			_histo.SetName(GetName());
18			TString title = "CALO dE/dx (plane ";
19			title += _plane;
20			title += ", ";
21			if (_view == 0)
22			title += "X";
23			if (_view == 1)
24			title += "Y";
25			title += " view) Vs Beta";
26			_histo.SetTitle(title);
27			_histo.GetXaxis()->SetTitle("Beta");
28			_histo.GetYaxis()->SetTitle("dE/dx (MIP)");
29
30			_histo.SetBins(nBinsBeta, minBeta, maxBeta, nBinsDedx, minDedx, maxDedx);
31			}
32
33			void CaloDedxVsBetaHistoAction::OnGood(PamLevel2 *event) {
34
35			_histo.Fill(event->GetTrack(0)->GetToFTrack()->beta[12], event->GetCaloLevel1()->qtotpl(_view, _plane));
36			//cout << event->GetTrack(0)->GetToFTrack()->beta[12] << " "<< event->GetCaloLevel1()->qtotpl(_view, _plane) << endl;
37			}
38
39			void CaloDedxVsBetaHistoAction::Finalize() {
40
41			TFile outFile(_outFileName, _mode);
42			outFile.cd();
43			_histo.Write();
44			outFile.Close();
45			}