Line data    Source code

       1             : #ifndef ALIMATHBASE_H
       2             : #define ALIMATHBASE_H
       3             : /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
       4             :  * See cxx source for full Copyright notice                               */
       5             : 
       6             : 
       7             :  
       8             : #include "TObject.h"
       9             : #include "TVectorD.h"
      10             : #include "TMatrixD.h"
      11             : #include "TGraph2D.h"
      12             : #include "TGraph.h"
      13             : 
      14             : class TH1F;
      15             : class TH3;
      16             : 
      17             :  
      18           0 : class AliMathBase : public TObject
      19             : {
      20             :  public:
      21             :   AliMathBase();
      22             :   virtual ~AliMathBase();
      23             :   static void    EvaluateUni(const Int_t nvectors, const Double_t *data, Double_t &mean, Double_t &sigma, const Int_t hSub);
      24             :   static void    EvaluateUniExternal(Int_t nvectors, Double_t *data, Double_t &mean, Double_t &sigma, Int_t hh, Float_t externalfactor=1);
      25             :   static Int_t  Freq(Int_t n, const Int_t *inlist, Int_t *outlist, Bool_t down);    
      26             :   static void TruncatedMean(TH1F * his, TVectorD *param, Float_t down=0, Float_t up=1.0, Bool_t verbose=kFALSE);
      27             :   static void LTM(TH1F * his, TVectorD *param=0 , Float_t fraction=1,  Bool_t verbose=kFALSE);
      28             :   static Double_t  FitGaus(TH1F* his, TVectorD *param=0, TMatrixD *matrix=0, Float_t xmin=0, Float_t xmax=0,  Bool_t verbose=kFALSE);
      29             :   static Double_t  FitGaus(Float_t *arr, Int_t nBins, Float_t xMin, Float_t xMax, TVectorD *param=0, TMatrixD *matrix=0, Bool_t verbose=kFALSE);
      30             :   static Float_t  GetCOG(Short_t *arr, Int_t nBins, Float_t xMin, Float_t xMax, Float_t *rms=0, Float_t *sum=0);
      31             : 
      32             :   static Double_t TruncatedGaus(Double_t mean, Double_t sigma, Double_t cutat);
      33             :   static Double_t TruncatedGaus(Double_t mean, Double_t sigma, Double_t leftCut, Double_t rightCut);
      34             : 
      35             :   static TGraph2D *  MakeStat2D(TH3 * his, Int_t delta0, Int_t delta1, Int_t type);
      36             :   static TGraph *  MakeStat1D(TH3 * his, Int_t delta1, Int_t type);
      37             : 
      38             :   static Double_t ErfcFast(Double_t x);                           // Complementary error function erfc(x)
      39    20544784 :   static Double_t ErfFast(Double_t x) {return 1-ErfcFast(x);}     // Error function erf(x)
      40             : 
      41             :   //
      42             :   // TestFunctions:
      43             :   //
      44             :  static  void TestGausFit(Int_t nhistos=5000);
      45             : 
      46             :   //
      47             :   // Bethe-Bloch formula parameterizations
      48             :   //
      49             :   static Double_t BetheBlochAleph(Double_t bg,
      50             :                                   Double_t kp1=0.76176e-1,
      51             :                                   Double_t kp2=10.632,
      52             :                                   Double_t kp3=0.13279e-4,
      53             :                                   Double_t kp4=1.8631,
      54             :                                   Double_t kp5=1.9479
      55             :                                   );
      56             : 
      57             :   //
      58             :   // Distributions
      59             :   //
      60             : 
      61             :   static Double_t Gamma(Double_t k=0);
      62             : 
      63         176 :  ClassDef(AliMathBase,0) // Various mathematical tools for physics analysis - which are not included in ROOT TMath
      64             :  
      65             : };
      66             : #endif