Update of B->Jpsi->ee analysis classes

[u/mrichter/AliRoot.git] / PWG3 / vertexingHF / AliBtoJPSItoEleCDFfitHandler.h

#ifndef ALIBTOJPSITOELECDFFITHANDLER_H\r
#define ALIBTOJPSITOELECDFFITHANDLER_H\r
/* Copyright(c) 1998-2009, ALICE Experiment at CERN, All rights reserved. *\r
 * See cxx source for full Copyright notice                               */\r
\r
//___________________________________________________________________________\r
//                  Class AliBtoJPSItoEleCDFfitHandler\r
//            Class to perform unbinned log-likelihood fit\r
//      \r
//                         Origin: C. Di Giglio\r
//     Contact: Carmelo.Digiglio@ba.infn.it; Giuseppe.Bruno@ba.infn.it\r
//____________________________________________________________________________\r
\r
#include <TNamed.h>\r
#include <TBits.h>\r
#include "TMath.h"\r
#include "AliBtoJPSItoEleCDFfitFCN.h"\r
#include "AliLog.h"\r
\r
class TBits; \r
\r
class AliBtoJPSItoEleCDFfitHandler : public TNamed {\r
 public:\r
  //\r
  AliBtoJPSItoEleCDFfitHandler();\r
  AliBtoJPSItoEleCDFfitHandler& operator= (const AliBtoJPSItoEleCDFfitHandler& c);\r
  AliBtoJPSItoEleCDFfitHandler(const AliBtoJPSItoEleCDFfitHandler& c);\r
  AliBtoJPSItoEleCDFfitHandler(Double_t* decaytime, Double_t* invariantmass, Int_t ncand);\r
  ~AliBtoJPSItoEleCDFfitHandler(); \r
  //Double_t Up() const { return fUp*TMath::Sqrt(fLikely->GetIntegral()); }\r
  Double_t Up() const { return fUp; }\r
  void SetErrorDef(Double_t up) {fUp = up;}\r
  void SetPrintStatus(Bool_t okPrintStatus) { fPrintStatus = okPrintStatus; } \r
  Bool_t GetPrintStatus() { return fPrintStatus ; }\r
  void SetParamStartValues (Double_t*);\r
  Double_t* GetStartParamValues() { return fParamStartValues; }\r
  TBits GetFixedParamList() { return fIsParamFixed; }\r
  void FixParam(UInt_t param, Bool_t value) { fIsParamFixed.SetBitNumber(param,value); }\r
  void FixAllParam(Bool_t value) { for(UInt_t par=0;par<16;par++) fIsParamFixed.SetBitNumber(par,value); }\r
  Bool_t IsParamFixed(UInt_t param) { return fIsParamFixed.TestBitNumber(param); }\r
  void SetResolutionConstants();\r
  void SetCrystalBallFunction(Bool_t okCB);\r
  void SetMassWndHigh(Double_t limit);\r
  void SetMassWndLow(Double_t limit);\r
\r
  Double_t operator()(const Double_t* par) const ;\r
  void CdfFCN(Int_t & /* npar */, Double_t * /* gin */, Double_t &f, Double_t *par, Int_t /* iflag */);\r
\r
  Double_t* Decaytime() const         { return fX; }\r
  Double_t* InvariantMass() const     { return fM; }\r
  AliBtoJPSItoEleCDFfitFCN* LikelihoodPointer() const { return fLikely; }\r
  Int_t DoMinimization();\r
\r
 private:\r
  //\r
  TBits fIsParamFixed;                               //array of bits: 0 = param free; 1 = param fixed;\r
  Bool_t fPrintStatus;                               //flag to enable the prit out of the algorithm at each step\r
  Double_t fParamStartValues[16];                    //array of parameters input value\r
  Double_t fUp;                                      //error definition \r
  Double_t* fX; 	                     	     //pseudo-proper decay time X\r
  Double_t* fM;                                      //invariant mass M\r
  AliBtoJPSItoEleCDFfitFCN* fLikely;                 //Log likelihood function\r
  Int_t fNcand;                                      //number of candidates\r
  //\r
  ClassDef(AliBtoJPSItoEleCDFfitHandler,1);\r
\r
}; \r
#endif\r