/[PAMELA software]/PamCut/CollectionActions/TrkDedxVsRigHistoAction/TrkDedxVsRigHistoAction.cpp
ViewVC logotype

Diff of /PamCut/CollectionActions/TrkDedxVsRigHistoAction/TrkDedxVsRigHistoAction.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1 by pam-fi, Wed Jun 10 12:41:41 2009 UTC revision 1.2 by pam-fi, Wed Aug 5 14:01:39 2009 UTC
# Line 10  Line 10 
10  #include "TrkDedxVsRigHistoAction.h"  #include "TrkDedxVsRigHistoAction.h"
11    
12  TrkDedxVsRigHistoAction::TrkDedxVsRigHistoAction(const char *actionName, TString outFileName, TString mode,  TrkDedxVsRigHistoAction::TrkDedxVsRigHistoAction(const char *actionName, TString outFileName, TString mode,
13      float minDedx, float maxDedx, unsigned int nBinsDedx, float minR, float maxR, unsigned int nBinsR) :      bool logBinningR, float minDedx, float maxDedx, unsigned int nBinsDedx, float minR, float maxR,
14        unsigned int nBinsR) :
15    CollectionAction(actionName), _outFileName(outFileName), _histo(), _mode(mode) {    CollectionAction(actionName), _outFileName(outFileName), _histo(), _mode(mode) {
16    
17    _histo.SetName(GetName());    _histo.SetName(GetName());
# Line 18  TrkDedxVsRigHistoAction::TrkDedxVsRigHis Line 19  TrkDedxVsRigHistoAction::TrkDedxVsRigHis
19    _histo.GetXaxis()->SetTitle("R (GV)");    _histo.GetXaxis()->SetTitle("R (GV)");
20    _histo.GetYaxis()->SetTitle("dE/dx (MIP)");    _histo.GetYaxis()->SetTitle("dE/dx (MIP)");
21    
22    _histo.SetBins(nBinsR, minR, maxR, nBinsDedx, minDedx, maxDedx);    if (!logBinningR || minR <= 0.)
23        _histo.SetBins(nBinsR, minR, maxR, nBinsDedx, minDedx, maxDedx);
24      else {
25        double rigBins[nBinsR + 1];
26        double maxExp = log10(maxR/minR);
27        for (unsigned int i = 0; i < nBinsR + 1; i++) {
28          rigBins[i] = minR * pow(10., (double)i / ((double)nBinsR) * maxExp);
29        }
30        double dedxBins[nBinsDedx + 1];
31        double dedxStep = (maxDedx - minDedx) / (double) nBinsDedx;
32        for (unsigned int i = 0; i < nBinsDedx + 1; i++) {
33          dedxBins[i] = minDedx + i * dedxStep;
34        }
35        _histo.SetBins(nBinsR, rigBins, nBinsDedx, dedxBins);
36      }
37  }  }
38    
39  void TrkDedxVsRigHistoAction::OnGood(PamLevel2 *event) {  void TrkDedxVsRigHistoAction::OnGood(PamLevel2 *event) {
40    
41    _histo.Fill(1/event->GetTrack(0)->GetTrkTrack()->GetDeflection(), event->GetTrack(0)->GetTrkTrack()->GetDEDX());    _histo.Fill(1 / event->GetTrack(0)->GetTrkTrack()->GetDeflection(), event->GetTrack(0)->GetTrkTrack()->GetDEDX());
42  }  }
43    
44  void TrkDedxVsRigHistoAction::Finalize() {  void TrkDedxVsRigHistoAction::Finalize() {
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.23