Histo2DActions/TofDedxVsBetaHistoAction/TofDedxVsBetaHistoAction.cpp

/*
 * TofDedxVsBetaHistoAction.cpp
 *
 *  Created on: 25-giu-2009
 *      Author: Nicola Mori
 */

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

#include "TofDedxVsBetaHistoAction.h"

TofDedxVsBetaHistoAction::TofDedxVsBetaHistoAction(const char *actionName, unsigned int layers, bool standAlone,
    TString outFileBase, TString mode, bool outRoot, bool outText, TString title) :
  Histo2DAction<Int_t> (actionName, title, outFileBase, mode, outRoot, outText), _layers(layers), _nLayers(0),
      _badEvents(0), _standAlone(standAlone) {

  bool buildTitle = false;
  if (title == "") {
    buildTitle = true;
    title = "TOF dE/dx (";
  }
  if ((_layers & S11) == S11) {
    if (buildTitle)
      title += "S11";
    _nLayers++;
  }
  if ((_layers & S12) == S12) {
    if (buildTitle) {
      if (_nLayers > 0)
        title += ",";
      title += "S12";
    }
    _nLayers++;
  }
  if ((_layers & S21) == S21) {
    if (buildTitle) {
      if (_nLayers > 0)
        title += ",";
      title += "S21";
    }
    _nLayers++;
  }
  if ((_layers & S22) == S22) {
    if (buildTitle) {
      if (_nLayers > 0)
        title += ",";
      title += "S22";
    }
    _nLayers++;
  }
  if ((_layers & S31) == S31) {
    if (buildTitle) {
      if (_nLayers > 0)
        title += ",";
      title += "S31";
    }
    _nLayers++;
  }
  if ((_layers & S32) == S32) {
    if (buildTitle) {
      if (_nLayers > 0)
        title += ",";
      title += "S32";
    }
    _nLayers++;
  }
  if (buildTitle)
    title += ") Vs Beta";

  SetTitle(title);

}

void TofDedxVsBetaHistoAction::OnGood(PamLevel2 *event) {

  static float dEdx;
  static float dEdxLayer;
  static unsigned int badLayers;

  static int trkSeqNo;
  if (_standAlone)
    trkSeqNo = 0;
  else
    for (int i = 0; i < event->GetToFLevel2()->ntrk(); i++) {
      if (event->GetToFLevel2()->GetToFTrkVar(i)->trkseqno == event->GetTrack(0)->GetTrkTrack()->GetSeqNo())
        trkSeqNo = i;
    }

  dEdx = 0.;
  badLayers = 0;
  if ((_layers & S11) == S11) {
    dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 0, 100);
    if (dEdxLayer > 0. && dEdxLayer < 4090)
      dEdx += dEdxLayer;
    else
      badLayers++;
  }
  if ((_layers & S12) == S12) {
    dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 1, 100);
    if (dEdxLayer > 0. && dEdxLayer < 4090)
      dEdx += dEdxLayer;
    else
      badLayers++;
  }
  if ((_layers & S21) == S21) {
    dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 2, 100);
    if (dEdxLayer > 0. && dEdxLayer < 4090)
      dEdx += dEdxLayer;
    else
      badLayers++;
  }
  if ((_layers & S22) == S22) {
    dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 3, 100);
    if (dEdxLayer > 0. && dEdxLayer < 4090)
      dEdx += dEdxLayer;
    else
      badLayers++;
  }
  if ((_layers & S31) == S31) {
    dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 4, 100);
    if (dEdxLayer > 0. && dEdxLayer < 4090)
      dEdx += dEdxLayer;
    else
      badLayers++;
  }
  if ((_layers & S32) == S32) {
    dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 5, 100);
    if (dEdxLayer > 0. && dEdxLayer < 4090)
      dEdx += dEdxLayer;
    else
      badLayers++;
  }
  if (badLayers == 0) {
    dEdx /= _nLayers;
    Fill(event->GetToFLevel2()->GetToFTrkVar(trkSeqNo)->beta[12], dEdx);
  }
  else
    _badEvents++;
}

1	/*
2	* TofDedxVsBetaHistoAction.cpp
3	*
4	* Created on: 25-giu-2009
5	* Author: Nicola Mori
6	*/
7
8	/! @file TofDedxVsBetaHistoAction.cpp The TofDedxVsBetaHistoAction class implementation file. /
9
10	#include "TofDedxVsBetaHistoAction.h"
11
12	TofDedxVsBetaHistoAction::TofDedxVsBetaHistoAction(const char *actionName, unsigned int layers, bool standAlone,
13	TString outFileBase, TString mode, bool outRoot, bool outText, TString title) :
14	Histo2DAction<Int_t> (actionName, title, outFileBase, mode, outRoot, outText), _layers(layers), _nLayers(0),
15	_badEvents(0), _standAlone(standAlone) {
16
17	bool buildTitle = false;
18	if (title == "") {
19	buildTitle = true;
20	title = "TOF dE/dx (";
21	}
22	if ((_layers & S11) == S11) {
23	if (buildTitle)
24	title += "S11";
25	_nLayers++;
26	}
27	if ((_layers & S12) == S12) {
28	if (buildTitle) {
29	if (_nLayers > 0)
30	title += ",";
31	title += "S12";
32	}
33	_nLayers++;
34	}
35	if ((_layers & S21) == S21) {
36	if (buildTitle) {
37	if (_nLayers > 0)
38	title += ",";
39	title += "S21";
40	}
41	_nLayers++;
42	}
43	if ((_layers & S22) == S22) {
44	if (buildTitle) {
45	if (_nLayers > 0)
46	title += ",";
47	title += "S22";
48	}
49	_nLayers++;
50	}
51	if ((_layers & S31) == S31) {
52	if (buildTitle) {
53	if (_nLayers > 0)
54	title += ",";
55	title += "S31";
56	}
57	_nLayers++;
58	}
59	if ((_layers & S32) == S32) {
60	if (buildTitle) {
61	if (_nLayers > 0)
62	title += ",";
63	title += "S32";
64	}
65	_nLayers++;
66	}
67	if (buildTitle)
68	title += ") Vs Beta";
69
70	SetTitle(title);
71
72	}
73
74	void TofDedxVsBetaHistoAction::OnGood(PamLevel2 *event) {
75
76	static float dEdx;
77	static float dEdxLayer;
78	static unsigned int badLayers;
79
80	static int trkSeqNo;
81	if (_standAlone)
82	trkSeqNo = 0;
83	else
84	for (int i = 0; i < event->GetToFLevel2()->ntrk(); i++) {
85	if (event->GetToFLevel2()->GetToFTrkVar(i)->trkseqno == event->GetTrack(0)->GetTrkTrack()->GetSeqNo())
86	trkSeqNo = i;
87	}
88
89	dEdx = 0.;
90	badLayers = 0;
91	if ((_layers & S11) == S11) {
92	dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 0, 100);
93	if (dEdxLayer > 0. && dEdxLayer < 4090)
94	dEdx += dEdxLayer;
95	else
96	badLayers++;
97	}
98	if ((_layers & S12) == S12) {
99	dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 1, 100);
100	if (dEdxLayer > 0. && dEdxLayer < 4090)
101	dEdx += dEdxLayer;
102	else
103	badLayers++;
104	}
105	if ((_layers & S21) == S21) {
106	dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 2, 100);
107	if (dEdxLayer > 0. && dEdxLayer < 4090)
108	dEdx += dEdxLayer;
109	else
110	badLayers++;
111	}
112	if ((_layers & S22) == S22) {
113	dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 3, 100);
114	if (dEdxLayer > 0. && dEdxLayer < 4090)
115	dEdx += dEdxLayer;
116	else
117	badLayers++;
118	}
119	if ((_layers & S31) == S31) {
120	dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 4, 100);
121	if (dEdxLayer > 0. && dEdxLayer < 4090)
122	dEdx += dEdxLayer;
123	else
124	badLayers++;
125	}
126	if ((_layers & S32) == S32) {
127	dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 5, 100);
128	if (dEdxLayer > 0. && dEdxLayer < 4090)
129	dEdx += dEdxLayer;
130	else
131	badLayers++;
132	}
133	if (badLayers == 0) {
134	dEdx /= _nLayers;
135	Fill(event->GetToFLevel2()->GetToFTrkVar(trkSeqNo)->beta[12], dEdx);
136	}
137	else
138	_badEvents++;
139	}
140