PndSdsChargeWeightingAlgorithms.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".              *
//****************************************************************************

#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