KVImpactParameter.h

Go to the documentation of this file.

 
#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>
 
class KVImpactParameter : 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);
 
public:
   KVImpactParameter(TH1*, Option_t* evol = "D");
   virtual ~KVImpactParameter();
 
   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);
 
   ClassDef(KVImpactParameter, 2) //Impact parameter analysis tools
};
 
#endif