PndCAParameters.h

Go to the documentation of this file.

/**************************************************************************
 * This file is property of and copyright by the ALICE HLT Project        *
 * All rights reserved.                                                   *
 *                                                                        *
 * Primary Authors:                                                       *
 *     Copyright 2009       Matthias Kretz <kretz@kde.org>                *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

#ifndef PNDCAPARAMETERS_H
#define PNDCAPARAMETERS_H

#include "PndCAVector.h"

struct PndCAParameters { 
  enum {
    MaxNStations = 30,
    NMVDStations = 4,
    MinimumHitsForRecoTrack = 3,
    MaxCellLength = 6,
    //    LastCellLength = 4 // length of rightmost cell
    LastCellLength = 3
  };

  // Transform CALocal coordinates in global x y z coordinates and back
  // CALocal coor are: X0 is normal to the station module, X2 is paralel to magnetic field. X0 X1 X2 are right handed coordinates
  template <typename T>
  static void CALocalToGlobal(T x0, T x1, T angle, T &x, T &y);
  template <typename T>
  static void GlobalToCALocal(T x, T y, T angle, T &x0, T &x1);
  template <typename T>
  static void CALocalToGlobal(T x0, T x1, T x2, T angle, T &x, T &y, T &z);
  template <typename T>
  static void GlobalToCALocal(T x, T y, T z, T angle, T &x0, T &x1, T &x2);

}; // namespace PndCAParameters
// namespace Parameters{  using namespace PndCAParameters; }

/*
template< typename T >
void PndCAParameters::CALocalToGlobal( T x0, T x1, T angle, T& x, T& y ) {
  const T cA = CAMath::Cos( angle );
  const T sA = CAMath::Sin( angle );
  x =  x0*sA + x1*cA;
  y =    x0*cA - x1*sA;
}

template< typename T >
void PndCAParameters::CALocalToGlobal( T x0, T x1, T x2, T angle, T& x, T& y, T &z ) {
  CALocalToGlobal<T>(x0,x1,angle,x,y);
  z =  - x2;
}
*/
template <typename T>
void PndCAParameters::GlobalToCALocal(T x, T y, T angle, T &x0, T &x1)
{
  const T cA = CAMath::Cos(angle);
  const T sA = CAMath::Sin(angle);
  // SG!!!
  // x0 = x*cA + y*sA;
  // x1 =  x*cA - y*cA;
  x0 = x * sA + y * cA;
  x1 = x * cA - y * sA;
}

template <typename T>
void PndCAParameters::GlobalToCALocal(T x, T y, T z, T angle, T &x0, T &x1, T &x2)
{
  GlobalToCALocal<T>(x, y, angle, x0, x1);
  x2 = -z;
}

namespace PndCAPParameters { 

enum {

  MinimumHitsForMCTrack = 6, // PndCAParameters::MinimumHitsForRecoTrack,

  MinimumMCPointsForMCTrack = MinimumHitsForMCTrack, // PndCAParameters::MinimumHitsForRecoTrack,

  MinimumHitsForRecoTrack = PndCAParameters::MinimumHitsForRecoTrack
}; // enum

static const float ExtraThreshold = 0.05;

static const float RefThreshold = 1.;

static const float MinTrackPurity = .75;
} //  namespace PndCAPParameters
namespace PParameters {
using namespace PndCAPParameters;
}

#endif