CaloCuts/CaloGeomCut/CaloGeomCut.cpp

/*
 * CaloGeomCut.cpp
 *
 *  Created on: 19-mar-2009
 *      Author: Sergio Ricciarini, Nicola Mori
 */

/*! @file CaloGeomCut.cpp The CaloCrossCut class implementation file */

#ifndef NO_CALOAXIS

#include "CaloGeomCut.h"
const Int_t CaloGeomCut::_nPoint = TrkParams::nGF + 1;
const Float_t CaloGeomCut::_zCaloTop = -26.181;

CaloGeomCut::CaloGeomCut(const char *cutName, CaloAxis *xCaloAxis, CaloAxis *yCaloAxis, UInt_t iMethod,
    Float_t xTolCaloTrack, Float_t yTolCaloTrack, Float_t mass, Float_t zeta, Float_t resMax, Float_t qualCut, Float_t chi2Cut) :
    PamCut(cutName), _xCaloAxis(xCaloAxis), _yCaloAxis(yCaloAxis), _iMethod(iMethod), _xTolCaloTrack(xTolCaloTrack), _yTolCaloTrack(
        yTolCaloTrack), _mass(mass),  _zeta(zeta), _resMax(resMax), _qualCut(qualCut), _chi2Cut(chi2Cut), _traj(NULL), _zIn(_nPoint) {

  if (_iMethod == 1) {
    _traj = (Trajectory*) new Trajectory(_nPoint);

    for (Int_t igf = 0; igf < TrkParams::nGF; igf++) {
      _zIn[igf] = TrkParams::zGF[igf];
    }
    _zIn[TrkParams::nGF] = _zCaloTop;

  }
}

/*! @brief Destructor. */
CaloGeomCut::~CaloGeomCut() {

  delete _traj;
}

int CaloGeomCut::Check(PamLevel2 *event) {

  Float_t xCaloTrack = 0.;
  Float_t yCaloTrack = 0.;

  Double_t tanx, tany, beta;
  Float_t alCaloTop[5]; // MUST be Float_t

  tanx = _xCaloAxis->par[1];
  tany = _yCaloAxis->par[1];

  if (_iMethod == 0) {
  }
  else if (_iMethod == 1) { // curved track in the Pamela acceptance (back-propagated starting from the straight track in the calorimeter, evaluated at the top of calorimeter, which is of course better than at higher z coordinates)

    alCaloTop[0] = _xCaloAxis->GetYfit(_zCaloTop); // Yfit   [xCaloTop]
    alCaloTop[1] = _yCaloAxis->GetYfit(_zCaloTop); // Yfit  [yCaloTop]
    alCaloTop[2] = sin(atan(sqrt(pow(tanx, 2) + pow(tany, 2)))); // [sintheta]
    alCaloTop[3] = TMath::Pi() + atan2(tany, tanx); // [phi]

    beta = event->GetToFLevel2()->CalcBeta(0, _resMax, _qualCut, _chi2Cut); // TOF stand-alone beta
    alCaloTop[4] = 1. / ((_mass/_zeta) * beta / sqrt(1. - pow(beta, 2))); // [etaP] [rho = (m/Z)*(beta/(sqrt(1-beta^2)))]

    _traj->DoTrack(alCaloTop, _zCaloTop);

  }
  else {
    return 1;
  }

  for (Int_t igf = 0; igf < TrkParams::nGF; igf++) {

    if (_iMethod == 0) { // straight track
      // NOTE: par[0] is the x (y) coordinate on plane z=0 in PAMELA reference (magnet centre)
      xCaloTrack = _xCaloAxis->par[0] + tanx * TrkParams::zGF[igf]; // cm
      yCaloTrack = _yCaloAxis->par[0] + tany * TrkParams::zGF[igf]; // cm
    }
    else if (_iMethod == 1) {
      xCaloTrack = _traj->x[igf];
      yCaloTrack = _traj->y[igf];
    }

    if (!(TrkParams::xGF_min[igf] + _xTolCaloTrack < xCaloTrack
        && xCaloTrack < TrkParams::xGF_max[igf] - _xTolCaloTrack)) {
      return 0;
    }

    if (!(TrkParams::yGF_min[igf] + _yTolCaloTrack < yCaloTrack
        && yCaloTrack < TrkParams::yGF_max[igf] - _yTolCaloTrack)) {
      return 0;
    }

  }

  return CUTOK;

}

#endif /* NO_CALOAXIS */
1	pam-fi	1.1	/*
2			* CaloGeomCut.cpp
3			*
4			* Created on: 19-mar-2009
5			* Author: Sergio Ricciarini, Nicola Mori
6			*/
7
8			/! @file CaloGeomCut.cpp The CaloCrossCut class implementation file /
9
10	pam-fi	1.4	#ifndef NO_CALOAXIS
11
12	pam-fi	1.1	#include "CaloGeomCut.h"
13	pam-fi	1.5	const Int_t CaloGeomCut::_nPoint = TrkParams::nGF + 1;
14			const Float_t CaloGeomCut::_zCaloTop = -26.181;
15
16			CaloGeomCut::CaloGeomCut(const char cutName, CaloAxis xCaloAxis, CaloAxis *yCaloAxis, UInt_t iMethod,
17	pam-fi	1.6	Float_t xTolCaloTrack, Float_t yTolCaloTrack, Float_t mass, Float_t zeta, Float_t resMax, Float_t qualCut, Float_t chi2Cut) :
18	pam-fi	1.5	PamCut(cutName), _xCaloAxis(xCaloAxis), _yCaloAxis(yCaloAxis), _iMethod(iMethod), _xTolCaloTrack(xTolCaloTrack), _yTolCaloTrack(
19	pam-fi	1.6	yTolCaloTrack), _mass(mass), _zeta(zeta), _resMax(resMax), _qualCut(qualCut), _chi2Cut(chi2Cut), _traj(NULL), _zIn(_nPoint) {
20	pam-fi	1.5
21			if (_iMethod == 1) {
22			_traj = (Trajectory*) new Trajectory(_nPoint);
23
24			for (Int_t igf = 0; igf < TrkParams::nGF; igf++) {
25			_zIn[igf] = TrkParams::zGF[igf];
26			}
27			_zIn[TrkParams::nGF] = _zCaloTop;
28
29			}
30			}
31
32			/! @brief Destructor. /
33			CaloGeomCut::~CaloGeomCut() {
34
35			delete _traj;
36			}
37	pam-fi	1.1
38			int CaloGeomCut::Check(PamLevel2 *event) {
39
40	pam-fi	1.3	Float_t xCaloTrack = 0.;
41			Float_t yCaloTrack = 0.;
42	pam-fi	1.5
43	pam-fi	1.2	Double_t tanx, tany, beta;
44			Float_t alCaloTop[5]; // MUST be Float_t
45	pam-fi	1.5
46			tanx = _xCaloAxis->par[1];
47			tany = _yCaloAxis->par[1];
48
49			if (_iMethod == 0) {
50	pam-fi	1.2	}
51	pam-fi	1.5	else if (_iMethod == 1) { // curved track in the Pamela acceptance (back-propagated starting from the straight track in the calorimeter, evaluated at the top of calorimeter, which is of course better than at higher z coordinates)
52
53			alCaloTop[0] = _xCaloAxis->GetYfit(_zCaloTop); // Yfit [xCaloTop]
54			alCaloTop[1] = _yCaloAxis->GetYfit(_zCaloTop); // Yfit [yCaloTop]
55			alCaloTop[2] = sin(atan(sqrt(pow(tanx, 2) + pow(tany, 2)))); // [sintheta]
56			alCaloTop[3] = TMath::Pi() + atan2(tany, tanx); // [phi]
57
58			beta = event->GetToFLevel2()->CalcBeta(0, _resMax, _qualCut, _chi2Cut); // TOF stand-alone beta
59	pam-fi	1.6	alCaloTop[4] = 1. / ((_mass/_zeta) * beta / sqrt(1. - pow(beta, 2))); // [etaP] [rho = (m/Z)*(beta/(sqrt(1-beta^2)))]
60	pam-fi	1.2
61	pam-fi	1.5	_traj->DoTrack(alCaloTop, _zCaloTop);
62	pam-fi	1.2
63			}
64			else {
65			return 1;
66			}
67	pam-fi	1.1
68			for (Int_t igf = 0; igf < TrkParams::nGF; igf++) {
69	pam-fi	1.5
70			if (_iMethod == 0) { // straight track
71	pam-fi	1.2	// NOTE: par[0] is the x (y) coordinate on plane z=0 in PAMELA reference (magnet centre)
72			xCaloTrack = _xCaloAxis->par[0] + tanx * TrkParams::zGF[igf]; // cm
73			yCaloTrack = _yCaloAxis->par[0] + tany * TrkParams::zGF[igf]; // cm
74			}
75	pam-fi	1.5	else if (_iMethod == 1) {
76			xCaloTrack = _traj->x[igf];
77			yCaloTrack = _traj->y[igf];
78	pam-fi	1.2	}
79	pam-fi	1.5
80			if (!(TrkParams::xGF_min[igf] + _xTolCaloTrack < xCaloTrack
81			&& xCaloTrack < TrkParams::xGF_max[igf] - _xTolCaloTrack)) {
82	pam-fi	1.1	return 0;
83			}
84
85	pam-fi	1.5	if (!(TrkParams::yGF_min[igf] + _yTolCaloTrack < yCaloTrack
86			&& yCaloTrack < TrkParams::yGF_max[igf] - _yTolCaloTrack)) {
87	pam-fi	1.1	return 0;
88			}
89	pam-fi	1.5
90	pam-fi	1.1	}
91
92			return CUTOK;
93	pam-fi	1.5
94	pam-fi	1.1	}
95	pam-fi	1.4
96			#endif /* NO_CALOAXIS */