LitMath.h

Go to the documentation of this file.

/* Copyright (C) 2009-2013 GSI/JINR-LIT, Darmstadt/Dubna
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Andrey Lebedev [committer] */
 
#ifndef LITMATH_H_
#define LITMATH_H_
 
#include "LitScalTrack.h"
#include "LitTypes.h"
 
namespace lit
{
  namespace parallel
  {
 
    inline fscal rcp(const fscal& a) { return 1. / a; }
 
    inline fscal sgn(const fscal& a) { return a < 0 ? -1 : (a > 0 ? 1 : 0); }
 
    inline unsigned short NDF(const LitScalTrack& track)
    {
      short ndf = 2 * track.GetNofHits() - 5;
      if (ndf > 0) {
        return ndf;
      }
      else {
        return 1;
      }
 
      // int ndf = 0;
      // for (int i = 0; i < track->GetNofHits(); i++){
      //    if (track->GetHit(i)->GetType() == kLITPIXELHIT) ndf += 2;
      //    else if (track->GetHit(i)->GetType() == kLITSTRIPHIT) ndf++;
      // }
      // ndf -= 5;
      // if (ndf > 0) return ndf; else return 1;
    }
 
  }  // namespace parallel
}  // namespace lit
#endif /* LITMATH_H_ */