PndSdsCalcPixel.h

Go to the documentation of this file.

#ifndef PNDSDSCALCPIXEL_H
#define PNDSDSCALCPIXEL_H

#include "PndSdsPixel.h"
#include "FairGeoVector.h"
#include <iostream>
#include <vector>
//#include <XYVector.h>
//#include "CLHEP/Vector/TwoVector.h"

class PndSdsCalcPixel {
 public:
  PndSdsCalcPixel();

  PndSdsCalcPixel(Double_t w, Double_t l, Double_t threshold = 0, Double_t noise = 0);

  std::vector<PndSdsPixel> GetPixels(Double_t inx, Double_t iny, Double_t outx, Double_t outy, Double_t energy);
  // Alternative in terms of I/O
  Int_t GetPixelsAlternative(Double_t inx, Double_t iny, Double_t outx, Double_t outy, Double_t energy, std::vector<Int_t> &cols, std::vector<Int_t> &rows,
                             std::vector<Double_t> &charges);
  void SetVerboseLevel(Int_t level) { fVerboseLevel = level; };
  std::ostream &operator<<(std::ostream &out);

 private:
  std::vector<PndSdsPixel> fPixels;
  Double_t fPixelWidth;
  Double_t fPixelLength;
  Double_t fThreshold;
  Double_t fNoise;
  Double_t fEnergy;
  Int_t fVerboseLevel;

  typedef enum { kQuadUNDEF, kUR, kUL, kDL, kDR } EQuadrantType;
  typedef enum { kPixelUNDEF, kU, kD, kL, kR } ENextPixelType;

  EQuadrantType fQuad;
  ENextPixelType fNextPixel;

  FairGeoVector fIn;  // entry point of particle in sensor
  FairGeoVector fOut; // exit point of particle out of sensor
  FairGeoVector fDir; // direction of particle inside sensor
  FairGeoVector fPos; // actual position of stepping inside sensor

  FairGeoVector fCon;

  PndSdsPixel fActivePixel;
  Double_t fCperL; // Charge per Tracklength;
  bool fStop;
  void CalcConMatrix();
  void ApplyConMatrix();
  void CalcQuadrant();
  void CalcCperL(Double_t Energy);
  void CalcStartPixel(); // pixel 0,0 is in lower left corner
  void CalcPixel();
  void ConvertPixels();

  Double_t SmearCharge(Double_t charge);
};

#endif