cavata_prescription.h

 
#ifndef __KVIMPACTPARAMETER_H
#define __KVIMPACTPARAMETER_H
#include "KVBase.h"
#include "TH1.h"
#include "TString.h"
#include "TGraph.h"
#include "TF1.h"
#include "KVHistoManipulator.h"
#include <vector>
 
namespace KVImpactParameters {
   class cavata_prescription : public KVBase {
      TH1* fData; 
      TString fEvol; 
      TGraph* fIPScale; 
      TGraph* fXSecScale; 
      TF1* fObsTransform; 
      TF1* fObsTransformXSec; 
      KVHistoManipulator HM; 
      Double_t Bmax; 
      Double_t Smax; 
 
      void make_scale(Int_t npoints);
#ifndef WITH_CPP11
      KVImpactParameter(const KVImpactParameter&) : KVBase() {}; // copying is not possible
#endif
   public:
#ifdef WITH_CPP11
      cavata_prescription(const cavata_prescription&) = delete; // copying is not possible
#endif
      cavata_prescription()
      {
         fData = nullptr;
         fIPScale = nullptr;
         fXSecScale = nullptr;
         fObsTransform = nullptr;
         fObsTransformXSec = nullptr;
         Bmax = 1.0;
      }
      cavata_prescription(TH1* h, Option_t* evol = "D");
      virtual ~cavata_prescription();
 
      void MakeScale(Int_t npoints = 100, Double_t bmax = 1.0);
      void MakeAbsoluteScale(Int_t npoints = 100, Double_t bmax = 1.0);
      TGraph* GetScale() const
      {
         return fIPScale;
      }
      TGraph* GetXSecScale() const
      {
         return fXSecScale;
      }
      Double_t BTransform(Double_t*, Double_t*);
      Double_t XTransform(Double_t*, Double_t*);
      TF1* GetTransFunc() const
      {
         return fObsTransform;
      }
      TF1* GetXSecTransFunc() const
      {
         return fObsTransformXSec;
      }
      Double_t GetImpactParameter(Double_t obs)
      {
         return fObsTransform->Eval(obs);
      }
      Double_t GetCrossSection(Double_t obs)
      {
         return fObsTransformXSec->Eval(obs);
      }
      Double_t GetObservable(Double_t b)
      {
         return fObsTransform->GetX(b);
      }
      Double_t GetObservableXSec(Double_t sigma)
      {
         return fObsTransformXSec->GetX(sigma);
      }
      TH1* GetIPDistribution(TH1* obs, Int_t nbinx = 100, Option_t* norm = "");
      TGraph* GetIPEvolution(TH2* obscor, TString moment, TString axis = "Y");
      TH1* GetXSecDistribution(TH1* obs, Int_t nbinx = 100, Option_t* norm = "");
      TGraph* GetXSecEvolution(TH2* obscor, TString moment, TString axis = "Y");
 
      static Double_t GetXSecFromIP(Double_t bmax)
      {
         return 10.*TMath::Pi() * pow(bmax, 2);
      }
      static Double_t GetIPFromXSec(Double_t xsec)
      {
         return pow(xsec / 10. / TMath::Pi(), 0.5);
      }
      std::vector<Double_t> SliceXSec(Int_t nslices, Double_t totXsec);
      double GetMeanBForSCA(double bmin, double bmax) const;
      double GetSigmaBForSCA(double bmin, double bmax) const;
 
      ClassDef(cavata_prescription, 2) //Impact parameter analysis tools
   };
}
#endif