/[PAMELA software]/PamCut/CaloCuts/CaloGeomCut/CaloGeomCut.cpp
ViewVC logotype

Diff of /PamCut/CaloCuts/CaloGeomCut/CaloGeomCut.cpp

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

revision 1.5 by pam-fi, Thu Jul 28 13:26:31 2011 UTC revision 1.6 by pam-fi, Fri Sep 2 16:25:38 2011 UTC
# Line 14  const Int_t CaloGeomCut::_nPoint = TrkPa Line 14  const Int_t CaloGeomCut::_nPoint = TrkPa
14  const Float_t CaloGeomCut::_zCaloTop = -26.181;  const Float_t CaloGeomCut::_zCaloTop = -26.181;
15    
16  CaloGeomCut::CaloGeomCut(const char *cutName, CaloAxis *xCaloAxis, CaloAxis *yCaloAxis, UInt_t iMethod,  CaloGeomCut::CaloGeomCut(const char *cutName, CaloAxis *xCaloAxis, CaloAxis *yCaloAxis, UInt_t iMethod,
17      Float_t xTolCaloTrack, Float_t yTolCaloTrack, Float_t mass, Float_t resMax, Float_t qualCut, Float_t chi2Cut) :      Float_t xTolCaloTrack, Float_t yTolCaloTrack, Float_t mass, Float_t zeta, Float_t resMax, Float_t qualCut, Float_t chi2Cut) :
18      PamCut(cutName), _xCaloAxis(xCaloAxis), _yCaloAxis(yCaloAxis), _iMethod(iMethod), _xTolCaloTrack(xTolCaloTrack), _yTolCaloTrack(      PamCut(cutName), _xCaloAxis(xCaloAxis), _yCaloAxis(yCaloAxis), _iMethod(iMethod), _xTolCaloTrack(xTolCaloTrack), _yTolCaloTrack(
19          yTolCaloTrack), _mass(mass), _resMax(resMax), _qualCut(qualCut), _chi2Cut(chi2Cut), _traj(NULL), _zIn(_nPoint) {          yTolCaloTrack), _mass(mass),  _zeta(zeta), _resMax(resMax), _qualCut(qualCut), _chi2Cut(chi2Cut), _traj(NULL), _zIn(_nPoint) {
20    
21    if (_iMethod == 1) {    if (_iMethod == 1) {
22      _traj = (Trajectory*) new Trajectory(_nPoint);      _traj = (Trajectory*) new Trajectory(_nPoint);
# Line 56  int CaloGeomCut::Check(PamLevel2 *event) Line 56  int CaloGeomCut::Check(PamLevel2 *event)
56      alCaloTop[3] = TMath::Pi() + atan2(tany, tanx); // [phi]      alCaloTop[3] = TMath::Pi() + atan2(tany, tanx); // [phi]
57    
58      beta = event->GetToFLevel2()->CalcBeta(0, _resMax, _qualCut, _chi2Cut); // TOF stand-alone beta      beta = event->GetToFLevel2()->CalcBeta(0, _resMax, _qualCut, _chi2Cut); // TOF stand-alone beta
59      alCaloTop[4] = 1. / (_mass * beta / sqrt(1. - pow(beta, 2))); // [etaP] [rho = (m/Z)*(beta/(sqrt(1-beta^2)))]      alCaloTop[4] = 1. / ((_mass/_zeta) * beta / sqrt(1. - pow(beta, 2))); // [etaP] [rho = (m/Z)*(beta/(sqrt(1-beta^2)))]
60    
61      _traj->DoTrack(alCaloTop, _zCaloTop);      _traj->DoTrack(alCaloTop, _zCaloTop);
62    
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.23