PndSdsChargeWeightingAlgorithms.h

Go to the documentation of this file.

#ifndef _PNDSDSCHARGEWEIGHTINGALGORITHMS_H
#define _PNDSDSCHARGEWEIGHTINGALGORITHMS_H

#include <vector>
#include "TObject.h"
#include "TClonesArray.h"
#include "TH2F.h"
#include "PndSdsDigiStrip.h"
#include "PndSdsChargeConversion.h"
class PndSdsCluster;
//class RunInfo;
class PndSdsCalcStrip;


class PndSdsChargeWeightingAlgorithms : public TObject
{
public:
        PndSdsChargeWeightingAlgorithms(TClonesArray* arr);
  PndSdsChargeWeightingAlgorithms(PndSdsChargeWeightingAlgorithms& other):
    TObject(other),
          fDigiArray(other.fDigiArray),
          fCalcStrip(other.fCalcStrip),
          fChargeConverter(other.fChargeConverter),
    fVerbose(other.fVerbose)
  {
  };
        virtual ~PndSdsChargeWeightingAlgorithms();
  PndSdsChargeWeightingAlgorithms& operator=(PndSdsChargeWeightingAlgorithms& other)
  {
          fDigiArray=other.fDigiArray;
          fCalcStrip=other.fCalcStrip;
          fChargeConverter=other.fChargeConverter;
    fVerbose=other.fVerbose;
    return *this;
  };

  std::pair<Double_t,Double_t> CenterOfGravity(const PndSdsCluster* Cluster);

        std::pair<Double_t,Double_t> HeadTail(const PndSdsCluster* Cluster);

        std::pair<Double_t,Double_t> Eta(const PndSdsCluster* Cluster, const TH2F* PosVsEta);

        std::pair<Double_t,Double_t> EtaValue(const PndSdsCluster* Cluster, Double_t &stripno, Int_t &NmbOfStrips); // calculates the eta value for the given cluster.
        std::pair<Double_t,Double_t> Binary(const PndSdsCluster* Cluster);

        std::pair<Double_t,Double_t> Median(const PndSdsCluster* Cluster);

        
        std::pair<Double_t,Double_t> AutoSelect(const PndSdsCluster* Cluster);

//      void MakedNdEta(const std::vector<PndSdsCluster>& Cluster, RunInfo& info);
        void SetDigiArray(TClonesArray* darray){fDigiArray = darray;};
        void SetCalcStrip(PndSdsCalcStrip* calc){fCalcStrip = calc;};
        void SetChargeConverter(PndSdsChargeConversion* ChargeConverter){fChargeConverter = ChargeConverter;};
  void SetVerbose(Int_t level=0){fVerbose=level;};

private:

        Double_t DigiCharge(Int_t digiIndex);
        Double_t DigiChargeError(Int_t digiIndex);
        Int_t DigiStripno(Int_t digiIndex);
        Double_t Erfmod(Double_t x, Double_t p0, Double_t p1, Double_t p2, Double_t p3);
        TClonesArray* fDigiArray;
        PndSdsCalcStrip* fCalcStrip;
        PndSdsChargeConversion* fChargeConverter;
  Int_t fVerbose;

        ClassDef(PndSdsChargeWeightingAlgorithms,1);
};
#endif