PndCAParam.h

Go to the documentation of this file.

//-*- Mode: C++ -*-
// @(#) $Id: PndCAParam.h,v 1.4 2011/10/01 00:23:44 perev Exp $
// ************************************************************************
// This file is property of and copyright by the ALICE HLT Project        *
// ALICE Experiment at CERN, All rights reserved.                         *
// See cxx source for full Copyright notice                               *
//                                                                        *
//*************************************************************************

#ifndef PNDCAPARAM_H
#define PNDCAPARAM_H

#include "PndCADef.h"
#include "PndCAVector.h"
#include "PndCAMath.h"
//#include "PndCATrackParam.h";
class PndCATrackParam;
class PndCATrackParamVector;
#include "PndCAParameters.h"
#include "PndCAStation.h"
#include "PndCAFieldValue.h"
#include <cstdio>
#include <sstream>

class PndCAParam {
  friend std::istringstream &operator>>(std::istringstream &, PndCAParam &);
  //   friend std::ostream &operator<<( std::ostream &, const PndCAParam & );
 public:
  PndCAParam();
  ~PndCAParam()
  {
    if (fStations)
      delete[] fStations;
  }
  int NStations() const { return fNStations; }

  const PndCAStation &Station(short i) const { return fStations[i]; }

  float Bz() const { return fBz; }
  float cBz() const { return fBz * 0.000299792458; }

  void SetBz(float v) { fBz = v; }

  void SetBoundaries(float z, float r)
  {
    fMaxZ = z;
    fMaxR = r;
  }

  float GetBz() const { return fBz; }
  float GetBz(float x, float y, float z) const;
  float_v GetBz(float_v x, float_v y, float_v z) const;
  float GetBz(const PndCATrackParam &t) const;
  float_v GetBz(const PndCATrackParamVector &t) const;

  float GetXOverX0(short iSt) const { return fStations[iSt].xOverX0; };
  float GetXTimesRho(short iSt) const { return fStations[iSt].xTimesRho; };
  float GetR(short iSt) const { return fStations[iSt].r; };
  float_v GetXOverX0(int_v iSt, float_m mask) const
  {
    float_v r;
    r.gather(fStations, &PndCAStation::xOverX0, static_cast<uint_v>(iSt), mask);
    return r;
  };
  float_v GetXTimesRho(int_v iSt, float_m mask) const
  {
    float_v r;
    r.gather(fStations, &PndCAStation::xTimesRho, static_cast<uint_v>(iSt), mask);
    return r;
  };
  float_v GetR(int_v iSt, const float_m &mask) const
  {
    float_v r;
    r.gather(fStations, &PndCAStation::r, static_cast<uint_v>(iSt), mask);
    return r;
  }

  float MinZ() const { return fMinZ; }
  float MaxZ() const { return fMaxZ; }
  float MinR() const { return fMinR; }
  float MaxR() const { return fMaxR; }

  const PndCAFieldValue &VtxFieldValue() const { return fVtxFieldValue; }

 protected:
  int fNStations;
  PndCAStation *fStations;
  float fBz;

  PndCAFieldValue fVtxFieldValue; // field at the vertex position.

  float fMinZ, fMaxZ, fMinR, fMaxR;
};

std::istringstream &operator>>(std::istringstream &, PndCAParam &);
// std::ostream &operator<<( std::ostream &, const PndCAParam & );

#endif