PndSdsCalcStrip.h

Go to the documentation of this file.

//
// 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 = 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 */