PndSdsCalcStrip.h

Go to the documentation of this file.

//****************************************************************************
//*                   This file is part of PandaRoot.                        *
//*                                                                          *
//*            PandaRoot is distributed under the terms of the               *
//*              GNU General Public License (GPL) version 3,                 *
//*                 copied verbatim in the file "LICENSE".                   *
//*                                                                          *
//*  Copyright (C) 2006 - 2024 FAIR GmbH and copyright holders of PandaRoot  *
//*     The copyright holders are listed in the file "COPYRIGHTHOLDERS".     *
//*               The authors are listed in the file "AUTHORS".              *
//****************************************************************************

//
// C++ Interface: PndSdsCalcStrip
//
// Description:
//
//
// Authors: HG Zaunick <hg.zaunick@physik.tu-dresden.de>, (C) 2007
//          R Kliemt (Tu Dresden,Uni Bonn)
//
// Copyright: See COPYING file that comes with this distribution
//
//
#ifndef PNDSDSCALCSTRIP_H
#define PNDSDSCALCSTRIP_H

#include "PndSdsStrip.h"
//#include "FairGeoVector.h"
#include "PndSdsStripDigiPar.h"

#include <TVector2.h>

#include <iostream>
#include <vector>

#include "PndDetectorList.h"


class PndSdsCalcStrip {
 public:
  PndSdsCalcStrip();

  PndSdsCalcStrip(Double_t pitch, Double_t orient, Int_t nrStrips, Int_t nrFeChannels, const TVector2 &firstStripAnchor, Double_t threshold, Double_t noise, Double_t csigma);

  PndSdsCalcStrip(const PndSdsStripDigiPar *digipar, SensorSide side = SensorSide::kTOP);

  std::vector<PndSdsStrip> GetStrips(Double_t inx, Double_t iny, Double_t inz, Double_t outx, Double_t outy, Double_t outz, Double_t eLoss);

  // Calculate without diffusion
  std::vector<PndSdsStrip> GetStripsNoDif(Double_t nuIn, Double_t nuOut, Double_t Q);

  // Calculate with diffusion
  std::vector<PndSdsStrip> GetStripsDif(Double_t nuIn, Double_t nuOut, Double_t Q);
  Double_t CalcFk(Double_t strip, Double_t x, Double_t sig);

  // for debugging; call one of the GetStrips* but returns the number of strips and int arrays to their channel/charge
  Int_t GetStripsAlternative(Double_t nuIn, Double_t nuOut, Double_t Q, Int_t mode, std::vector<Int_t> &indice, std::vector<Double_t> &charges);

  Int_t CalcFEfromStrip(Int_t stripNr) const;

  Int_t CalcChannelfromStrip(Int_t stripNr) const;

  void CalcFeChToStrip(Int_t fe, Int_t channel, Int_t &strip, enum SensorSide &side) const;
  void CalcStripPointOnLine(const Double_t strip, TVector2 &point) const;

  const TVector2 GetStripDirection() const { return fStripDir; }

  void SetVerboseLevel(Int_t level) { fVerboseLevel = level; }
  // std::ostream& operator<<(std::ostream& out);  // not implemented
  Double_t GetThreshold() const { return fThreshold; }
  Double_t GetNoise() const { return fNoise; }

  void Print() const;

 private:
  Double_t fPitch;     
  Double_t fOrient;    
  Int_t fNrStrips;     
  Int_t fNrFeChannels; 
  TVector2 fAnchor;    
  Double_t fThreshold; 
  Double_t fNoise;     
  Double_t fCSigma;    
  TVector2 fStripDir;  
  TVector2 fOrthoDir;  
  Int_t fVerboseLevel;

  Double_t CalcStripFromPoint(Double_t x, Double_t y);

  Double_t ChargeFromEloss(Double_t eloss) const { return eloss / (3.61e-9); }
  void InjectStripCharge(std::vector<PndSdsStrip> &array, Int_t istrip, Double_t charge);
  Double_t SmearCharge(Double_t charge);
};

#endif /* PNDSDSCALCSTRIP_H */