X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=FASTSIM%2FAliQuenchingWeights.h;h=0f5d088afe7f7a1aff1f552ae200eeba9c6f97f5;hb=c039a4b5c88e0b45514449ad69c715aa5a4cb0d2;hp=33de5028d51f31c96b5dcf465cb8c24379e8a9bd;hpb=9d851d20e4c4fa1d0e0ace5e76f3f98b71a98e2f;p=u%2Fmrichter%2FAliRoot.git diff --git a/FASTSIM/AliQuenchingWeights.h b/FASTSIM/AliQuenchingWeights.h index 33de5028d51..0f5d088afe7 100644 --- a/FASTSIM/AliQuenchingWeights.h +++ b/FASTSIM/AliQuenchingWeights.h @@ -17,28 +17,40 @@ // A. Dainese andrea.dainese@pd.infn.it //---------------------------------------------------------------------------- +#include #include class TH1F; class AliQuenchingWeights : public TObject { public: - enum kECMethod {kDefault=0,kReweight=1}; + enum kECMethod {kDefault=0,kReweight=1,kReweightCont=2}; AliQuenchingWeights(); AliQuenchingWeights(const AliQuenchingWeights& a); + AliQuenchingWeights& operator=(const AliQuenchingWeights& a) + {a.Copy(*this); return(*this);} virtual ~AliQuenchingWeights(); - + void Init(); void Reset(); Int_t SampleEnergyLoss(); - Int_t SampleEnergyLoss(Int_t ipart, Double_t R); + Int_t SampleEnergyLoss(Int_t ipart, Double_t r); Double_t GetELossRandom(Int_t ipart, Double_t length, Double_t e=1.e10) const; Double_t CalcQuenchedEnergy(Int_t ipart, Double_t length, Double_t e) const; Double_t GetELossRandom(Int_t ipart, TH1F *hell, Double_t e=1.e10) const; Double_t CalcQuenchedEnergy(Int_t ipart, TH1F *hell, Double_t e) const; - Double_t GetELossRandomK(Int_t ipart, Double_t I0, Double_t I1, Double_t e=1.e10); Double_t CalcQuenchedEnergyK(Int_t ipart, Double_t I0, Double_t I1, Double_t e); + Double_t GetELossRandomKFast(Int_t ipart, Double_t I0, Double_t I1, Double_t e=1.e10); + Double_t GetELossRandomKFastR(Int_t ipart, Double_t r, Double_t wc, Double_t e=1.e10); + Double_t CalcQuenchedEnergyKFast(Int_t ipart, Double_t I0, Double_t I1, Double_t e); + + Double_t GetDiscreteWeight(Int_t ipart, Double_t I0, Double_t I1); + Double_t GetDiscreteWeightR(Int_t ipart, Double_t r); + void GetZeroLossProb(Double_t &p,Double_t &prw,Double_t &prwcont, + Int_t ipart,Double_t I0,Double_t I1,Double_t e=1.e10); + void GetZeroLossProbR(Double_t &p,Double_t &prw, Double_t &prwcont, + Int_t ipart,Double_t r,Double_t wc,Double_t e=1.e10); //multiple soft scattering approximation Int_t InitMult(const Char_t *contall="$(ALICE_ROOT)/FASTSIM/data/cont_mult.all", @@ -88,7 +100,10 @@ class AliQuenchingWeights : public TObject { {return CalcQk(fK,I0,I1);} Double_t CalcQk(Double_t k, Double_t I0, Double_t I1) const - {return I0*I0/2/I1/fgkConvFmToInvGeV/fgkConvFmToInvGeV*k;} + {return I0*I0/2./I1/fgkConvFmToInvGeV/fgkConvFmToInvGeV*k;} + + Double_t CalcLk(Double_t i0, Double_t i1) const + {return 2.*i1/i0;} Int_t CalcLengthMax(Double_t q) const {Double_t l3max=fgkRMax/.5/q/fgkConvFmToInvGeV/fgkConvFmToInvGeV; @@ -105,44 +120,44 @@ class AliQuenchingWeights : public TObject { Float_t GetMu() const {return fMu;} Float_t GetQTransport() const {return fQTransport;} Float_t GetK() const {return fK;} - Bool_t GetECMethod() const {return fECMethod;} + kECMethod GetECMethod() const {return fECMethod;} Bool_t GetTablesLoaded() const {return fTablesLoaded;} Bool_t GetMultSoft() const {return fMultSoft;} Int_t GetLengthMax() const {return fLengthMax;} TH1F* ComputeQWHisto (Int_t ipart,Double_t medval,Double_t length) const; TH1F* ComputeQWHistoX(Int_t ipart,Double_t medval,Double_t length) const; - TH1F* ComputeQWHistoX(Int_t ipart,Double_t R) const; + TH1F* ComputeQWHistoX(Int_t ipart,Double_t r) const; TH1F* ComputeELossHisto(Int_t ipart,Double_t medval,Double_t l,Double_t e=1.e10) const; TH1F* ComputeELossHisto(Int_t ipart,Double_t medval,TH1F *hEll,Double_t e=1.e10) const; - TH1F* ComputeELossHisto(Int_t ipart,Double_t R) const; + TH1F* ComputeELossHisto(Int_t ipart,Double_t r) const; Double_t GetMeanELoss(Int_t ipart,Double_t medval,Double_t l) const; Double_t GetMeanELoss(Int_t ipart,Double_t medval,TH1F *hEll) const; - Double_t GetMeanELoss(Int_t ipart,Double_t R) const; + Double_t GetMeanELoss(Int_t ipart,Double_t r) const; - void PlotDiscreteWeights(Double_t len=4) const; - void PlotContWeights(Int_t itype,Double_t len) const; - void PlotContWeightsVsL(Int_t itype,Double_t medval) const; - void PlotAvgELoss(Double_t len ,Double_t e=1.e10) const; - void PlotAvgELoss(TH1F *hEll,Double_t e=1.e10) const; - void PlotAvgELossVsL(Double_t e=1.e10) const; - void PlotAvgELossVsPt(Double_t medval,Double_t len) const; - void PlotAvgELossVsPt(Double_t medval,TH1F *hEll) const; + void PlotDiscreteWeights(Double_t len=4,Double_t qm=5) const; + void PlotContWeights(Int_t itype,Double_t len) const; + void PlotContWeightsVsL(Int_t itype,Double_t medval) const; + void PlotAvgELoss(Double_t len,Double_t qm=5,Double_t e=1.e10) const; + void PlotAvgELoss(TH1F *hEll,Double_t e=1.e10) const; + void PlotAvgELossVsL(Double_t e=1.e10) const; + void PlotAvgELossVsPt(Double_t medval,Double_t len) const; + void PlotAvgELossVsPt(Double_t medval,TH1F *hEll) const; protected: Int_t GetIndex(Double_t len) const; static const Double_t fgkConvFmToInvGeV; //conversion factor - static const Int_t fgBins; //number of bins for hists - static const Double_t fgMaxBin; //max. value of wc + static const Int_t fgkBins; //number of bins for hists + static const Double_t fgkMaxBin; //max. value of wc static const Double_t fgkRMax; //max. tabled value of R static Int_t fgCounter;//static instance counter Int_t fInstanceNumber; //instance number of class Bool_t fMultSoft; //approximation type - Bool_t fECMethod; //energy constraint method + kECMethod fECMethod; //energy constraint method Double_t fQTransport; //transport coefficient [GeV^2/fm]] Double_t fMu; //Debye screening mass Double_t fK; //proportional constant [fm]