/*
 * TrkDedxVsRigHistoAction.cpp
 *
 *  Created on: 9-giu-2009
 *      Author: Nicola Mori
 */

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

#include "TrkDedxVsRigHistoAction.h"

TrkDedxVsRigHistoAction::TrkDedxVsRigHistoAction(const char *actionName, TString outFileName, TString mode,
    float minDedx, float maxDedx, unsigned int nBinsDedx, float minR, float maxR, unsigned int nBinsR) :
  CollectionAction(actionName), _outFileName(outFileName), _histo(), _mode(mode) {

  _histo.SetName(GetName());
  _histo.SetTitle("TRK dE/dx Vs Rigidity");
  _histo.GetXaxis()->SetTitle("R (GV)");
  _histo.GetYaxis()->SetTitle("dE/dx (MIP)");

  _histo.SetBins(nBinsR, minR, maxR, nBinsDedx, minDedx, maxDedx);
}

void TrkDedxVsRigHistoAction::OnGood(PamLevel2 *event) {

  _histo.Fill(1/event->GetTrack(0)->GetTrkTrack()->GetDeflection(), event->GetTrack(0)->GetTrkTrack()->GetDEDX());
}

void TrkDedxVsRigHistoAction::Finalize() {

  TFile outFile(_outFileName, _mode);
  outFile.cd();
  _histo.Write();
  outFile.Close();
}