PndEmcXClMoments.h

Go to the documentation of this file.

//--------------------------------------------------------------------------
// File and Version Information:
//
// Description:
//  Class PndEmcXClMoments.
//
// Environment:
//  Software developed for the BaBar Detector at the SLAC B-Factory.
//
// Author List:
//  Thorsten Brandt     Originator
//
// Copyright Information:
//
// Dima Melnichuk, adaption for pandaroot
//
//------------------------------------------------------------------------
//#pragma once
#ifndef PNDEMCXCLMOMENTS_H
#define PNDEMCXCLMOMENTS_H

#include "PndEmcAbsClusterProperty.h"
#include "TObject.h"
#include <iostream>
#include <vector>

class TVector3;
class PndEmcDigi;
class PndEmcCluster;

struct PndEmcClEnergyDeposition {
  Double_t deposited_energy;
  Double_t r;
  Double_t phi;
};

class PndEmcXClMoments : public PndEmcAbsClusterProperty {

 public:
  // Constructors

  PndEmcXClMoments(const PndEmcCluster &cluster, const TClonesArray *digiArray);
  PndEmcXClMoments(const PndEmcXClMoments &);

  // Destructor

  virtual ~PndEmcXClMoments();

  // Operators

  // Selectors (const)

  // Second moment of energy about centroid
  virtual Double_t SecondMoment() const;
  virtual Double_t SecondMomentPhi() const;
  virtual Double_t SecondMomentTheta() const;

  // Absolute value Zernike-Moment of index(m,n); according to the
  // definition of the Zernike moments, n and m must satisfy
  // (a) n>=m  (b) n-m even  (c) n,m>=0
  // Otherwise this function returns -1;
  virtual Double_t AbsZernikeMoment(int n, int m, Double_t R0 = 15) const;

  // Parameter for description of the LATeral energydeposition
  // within the cluster; if cluster contains less than 3 digis,
  // the function return 0
  virtual Double_t Lat() const;

  virtual void Print(const Option_t *opt = "") const;

  // Modifiers

 protected:
  // Helper functions

  void Init();

  // explicit implementation of polynomial part of
  // Zernike-Functions for n<=5;

  Double_t f00(Double_t r) const;
  Double_t f11(Double_t r) const;
  Double_t f20(Double_t r) const;
  Double_t f22(Double_t r) const;
  Double_t f31(Double_t r) const;
  Double_t f33(Double_t r) const;
  Double_t f40(Double_t r) const;
  Double_t f42(Double_t r) const;
  Double_t f44(Double_t r) const;
  Double_t f51(Double_t r) const;
  Double_t f53(Double_t r) const;
  Double_t f55(Double_t r) const;

  // Calculation of Zernike-Moments for n<=5 :
  Double_t Fast_AbsZernikeMoment(int n, int m, Double_t R0) const;

  // Calculation of Zernike-Moments for general values of (n,m)
  Double_t Calc_AbsZernikeMoment(int n, int m, Double_t R0) const;

  Double_t Fak(int n) const;

 private:
  // These should never get used
  PndEmcXClMoments();
  PndEmcXClMoments &operator=(const PndEmcXClMoments &);

  // Data members
  std::vector<PndEmcClEnergyDeposition> *fEnergyDistribution;
  Int_t fClusterSize;
  Double_t (PndEmcXClMoments::*fFcn[12])(Double_t) const;
  //    Double_t (PndEmcXClMoments::*fFcn2)( Double_t ); // STE: Is it needed??

  ClassDef(PndEmcXClMoments, 1)
};

#endif // PNDEMCXCLMOMENTS_HH