FTSCANPletsV.h

Go to the documentation of this file.

//-*- Mode: C++ -*-
// *****************************************************************************
//                                                                             *
// @Autors: I.Kulakov; M.Pugach; M.Zyzak; I.Kisel                                        *
// @e-mail: I.Kulakov@gsi.de; M.Pugach@gsi.de; M.Zyzak@gsi.de; I.Kisel@compeng.uni-frankfurt.de *
//                                                                             *
// *****************************************************************************

#ifndef FTSCANPletsV_H
#define FTSCANPletsV_H

#include <vector>
using std::vector;
#include "FTSCAStationArray.h"
#include "PndFTSCATrackParamVector.h"

class FTSCANPletV {

  //* Parameters of N-plet in vectors

 public:
  FTSCANPletV(int n = 0) : fParam(), fIsValid(false) { fIHit.resize(n); }

  FTSCANPletV(const FTSCANPletV &s1, const PndFTSCATrackParamVector &param, const float_m isValid) : fParam(param), fIsValid(isValid) { fNHits = s1.N() + 1; }

  FTSCANPletV(const FTSCANPletV &s1, int ista, unsigned int ihit, const PndFTSCATrackParamVector &param, const float_m isValid) : fParam(param), fIsValid(isValid)
  {
    fNHits = s1.N() + 1;
    fIHit.resize(s1.N() + 1);
    for (int i = 0; i < N() - 1; i++)
      fIHit[i] = s1.IHit(i);
    fIHit[N() - 1] = TESV(ista, ihit);
  }

  FTSCANPletV(const uint_v &ihit, const int_v &ista, const PndFTSCATrackParamVector &param, const float_m &isValid) : fParam(param), fIsValid(isValid)
  {
    fIHit.resize(1, TESV(ista, ihit));
  }
  FTSCANPletV(const FTSCANPletV &s1, const FTSCANPletV &s2, const int iV, const PndFTSCATrackParamVector &param, const float_m isValid) : fParam(param), fIsValid(isValid)
  {
    fIHit.resize(s1.N() + 1);
    assert(s1.N() == s2.N());

    for (int i = 0; i < N() - 1; i++)
      fIHit[i] = s1.IHit(i);
    fIHit[N() - 1] = s2.IHit(N() - 2)[iV];
  }

  int N() const { return fNHits; /*fIHit.size();*/ }

  const TESV &IHit(int IH) const { return fIHit[IH]; }
  const vector<TESV> &IHit() const { return fIHit; }
  const PndFTSCATrackParamVector &Param() const { return fParam; }
  PndFTSCATrackParamVector &Param() { return fParam; }
  PndFTSCATrackParamVector &ParamRef() { return fParam; } // dbg

  float_m IsValid() const { return fIsValid; }

  // check wether a is neighbour from the right to this
  float_m IsRightNeighbour(const FTSCANPletV &a, unsigned int iV) const
  {
    float_m active = fIsValid;
    assert(a.N() >= N() - 1);
    for (int i = 0; i < N() - 1; i++) {
      active &= IHit(i + 1) == a.IHit(i)[iV];
    }
    return active;
  }

  void CopyOne(int iV, FTSCANPletV a, int iVa)
  {
    for (int i = 0; i < N(); i++) {
      fIHit[i].s[iV] = a.fIHit[i].s[iVa];
      fIHit[i].e[iV] = a.fIHit[i].e[iVa];
    }
    fParam.SetTrackParamOne(iV, a.fParam, iVa);
    fIsValid = (float_m(uint_v(Vc::IndexesFromZero) == uint_v(iV)) && float_m(a.fIsValid[iVa])) || (!float_m(uint_v(Vc::IndexesFromZero) == uint_v(iV)) && fIsValid);
  }

  int_v fLastHit;
  int fNHits;
  // private:
  vector<TESV> fIHit; // index of hit on station
  PndFTSCATrackParamVector fParam;

  float_m fIsValid;
};

class FTSCANPletsV : public FTSCAStationArray<FTSCANPletV> {
 public:
  FTSCANPletsV(){};
  FTSCANPletsV(int nSta, const FTSCAHits *hits) : FTSCAStationArray<FTSCANPletV>(nSta, hits){};
};

#endif