// Origin: Petr Naumenko, SPbSU-CERN, Petr.Naoumenko@cern.ch
//-------------------------------------------------------------------------
-#include "Riostream.h"
#include "AliLRCAnalysis.h"
+#include "Riostream.h"
#include "TFile.h"
#include "AliLRCFit.h"
#include "TProfile.h"
#include "TH2D.h"
#include "TPaveText.h"
#include "TF1.h"
-//#include "math.h"
+#include "TMath.h"
#include "TStyle.h"
class gStyle;
* AliLRCAnalysis class
******************************************************/
-AliLRCAnalysis::AliLRCAnalysis(): TObject(), fPrAbs(new TH1D()), fPrRel(new TH1D()), fPrf(new TH1D()), fPrb(new TH1D()), fileHist(new TFile()), fdptb(.0), fEntries(0), fSx((char*)" "), fSy((char*)" "), fxFitMin(.0), fxFitMax(.0), farel(.0), fbrel(.0), farelError(0.), fbrelError(0.), fXi2rel(.0), faabs(.0), fbabs(.0), faabsError(0.), fbabsError(0.), fXi2abs(.0){
+AliLRCAnalysis::AliLRCAnalysis(): fPrAbs(new TH1D()), fPrRel(new TH1D()), fPrf(new TH1D()), fPrb(new TH1D()), fileHist(new TFile()), fdptb(.0), fEntries(0), fSx((char*)" "), fSy((char*)" "), fxFitMin(.0), fxFitMax(.0), farel(.0), fbrel(.0), farelError(0.), fbrelError(0.), fXi2rel(.0), faabs(.0), fbabs(.0), faabsError(0.), fbabsError(0.), fXi2abs(.0){
//Empty constructor
}
-AliLRCAnalysis::AliLRCAnalysis(const AliLRCAnalysis& a): TObject(), fPrAbs(a.fPrAbs), fPrRel(a.fPrRel), fPrf(a.fPrf), fPrb(a.fPrb), fileHist(a.fileHist), fdptb(a.fdptb), fEntries(a.fEntries), fSx(a.fSx), fSy(a.fSy), fxFitMin(a.fxFitMin), fxFitMax(a.fxFitMax), farel(a.farel), fbrel(a.fbrel), farelError(a.farelError), fbrelError(a.fbrelError), fXi2rel(a.fXi2rel), faabs(a.faabs), fbabs(a.fbabs), faabsError(a.faabsError), fbabsError(a.fbabsError), fXi2abs(a.fXi2abs){
+AliLRCAnalysis::AliLRCAnalysis(const AliLRCAnalysis& a):fPrAbs(a.fPrAbs), fPrRel(a.fPrRel), fPrf(a.fPrf), fPrb(a.fPrb), fileHist(a.fileHist), fdptb(a.fdptb), fEntries(a.fEntries), fSx(a.fSx), fSy(a.fSy), fxFitMin(a.fxFitMin), fxFitMax(a.fxFitMax), farel(a.farel), fbrel(a.fbrel), farelError(a.farelError), fbrelError(a.fbrelError), fXi2rel(a.fXi2rel), faabs(a.faabs), fbabs(a.fbabs), faabsError(a.faabsError), fbabsError(a.fbabsError), fXi2abs(a.fXi2abs){
//Constructor
}
fxFitMax = a.fxFitMax;
farel = a.farel;
fbrel = a.fbrel;
+ farelError = a.farelError;
+ fbrelError = a.fbrelError;
fXi2rel = a.fXi2rel;
faabs = a.faabs;
fbabs = a.fbabs;
+ faabsError = a.faabsError;
+ fbabsError = a.fbabsError;
fXi2abs = a.fXi2abs;
}
return *this;
fSx = atitleF;
fSy = atitleB;
double mnf = fPrf->GetMean();
- fxFitMin = mnf-2*sqrt(mnf);
- fxFitMax = mnf+2*sqrt(mnf);
+ fxFitMin = mnf-2*TMath::Sqrt(mnf);
+ fxFitMax = mnf+2*TMath::Sqrt(mnf);
//delete profX;
}
void AliLRCAnalysis::SetBinsRange(int binMin, int binMax){
//Set the bin range
TH1D* h=fPrf;
- double n=h->GetNbinsX();
+ Int_t n=h->GetNbinsX();
fxFitMin = h->GetXaxis()->GetXmin();
fxFitMax = h->GetXaxis()->GetXmax();
double df = (fxFitMax-fxFitMin)/n;
}
void AliLRCAnalysis::DrawAbs() {
- //Draw abs var hist with ALL info
- int * mas = new int[N_PL_FLAGS];
- for ( int i = 0; i < N_PL_FLAGS; i++ )
- mas[i] = 1;
- DrawAbsPure( mas, 1 );
- delete [] mas;
- mas = NULL;
+//Draw abs var hist with ALL info
+ int * mas = new int [N_PL_FLAGS];
+ for ( int i = 0; i < N_PL_FLAGS; i++ )
+ mas[i] = 1;
+ DrawAbsPure( mas, 1 );
+ delete []mas;
}
void AliLRCAnalysis::DrawAbs( int * mas ) {
}
void AliLRCAnalysis::DrawRel() {
- //Draw rel var hist with ALL info
- int * mas = new int[N_PL_FLAGS];
- for ( int i = 0; i < N_PL_FLAGS; i++ )
- mas[i] = 1;
- DrawRelPure( mas, 1 );
- delete [] mas;
- mas = NULL;
+//Draw rel var hist with ALL info
+ int * mas = new int [N_PL_FLAGS];
+ for ( int i = 0; i < N_PL_FLAGS; i++ )
+ mas[i] = 1;
+ DrawRelPure( mas, 1 );
+ delete []mas;
}
void AliLRCAnalysis::DrawRel( int * mas ) {
fPrAbs->SetBinContent(i, profX->GetBinContent(i));
if(fPrf->GetBinContent(i)!=0)
{
- fPrAbs->SetBinError(i,sqrt(profX->GetBinContent(i)/fPrf->GetBinContent(i)));
+ fPrAbs->SetBinError(i,TMath::Sqrt(profX->GetBinContent(i)/fPrf->GetBinContent(i)));
}
fPrRel->SetBinContent(i, fPrAbs->GetBinContent(i)/fPrb->GetMean());
fPrRel->SetBinError(i,fPrAbs->GetBinError(i)/fPrb->GetMean());
if(fPrf->GetBinContent(i)!=0)
{
pt = profX->GetBinContent(i);
- fPrAbs->SetBinError(i,ptd*sqrt(Integral(nb,i))/fPrf->GetBinContent(i));
+ fPrAbs->SetBinError(i,ptd*TMath::Sqrt(Integral(nb,i))/fPrf->GetBinContent(i));
}
fPrRel->SetBinContent(i, fPrAbs->GetBinContent(i)/fPrb->GetMean());
fPrRel->SetBinError(i,fPrAbs->GetBinError(i)/fPrb->GetMean());
if(fPrf->GetBinContent(i)!=0)
{
pt = profX->GetBinContent(i);
- fPrAbs->SetBinError(i,ptd*sqrt(nb->GetBinContent(i)/fPrf->GetBinContent(i)));
+ fPrAbs->SetBinError(i,ptd*TMath::Sqrt(nb->GetBinContent(i)/fPrf->GetBinContent(i)));
}
fPrRel->SetBinContent(i, fPrAbs->GetBinContent(i)/fPrb->GetMean());
fPrRel->SetBinError(i,fPrAbs->GetBinError(i)/fPrb->GetMean());
if ( noError )
{
- cout << "After rounding: " << error << endl;
+ cout << "After rounding: " << error << endl;
return error;
}
else
+#ifndef ALILRCANALYSIS_H
+#define ALILRCANALYSIS_H
+
//-------------------------------------------------------------------------
-// Description:
-// This class is included into LRC library for Long-Range Correlation analysis
-// it is base class for NN, PtN, PtPt
-// implements base methods for thees classes
-// Origin: Petr Naumenko, SPbSU-CERN, Petr.Naoumenko@cern.ch
+// Description:
+// This class is included into LRC library for Long-Range Correlation analysis
+// it is base class for NN, PtN, PtPt
+// implements base methods for thees classes
+// Origin: Petr Naumenko, SPbSU-CERN, Petr.Naoumenko@cern.ch
//-------------------------------------------------------------------------
-#ifndef ALILRCANALYSIS_H
-#define ALILRCANALYSIS_H
/* See cxx source for full Copyright notice */
// Origin: Petr Naumenko, SPbSU-CERN, Petr.Naoumenko@cern.ch
//-------------------------------------------------------------------------
-#include <TObject.h>
-
#define N_PL_FLAGS 10
+#include "TObject.h"
+
class TFile;
class AliLRCFit;
class TProfile;
class math;
class TStyle;
-class AliLRCAnalysis : public TObject{
+class AliLRCAnalysis{
public:
void DrawAbs();
void DrawAbs( int * mas );
double HI2(TH1D * const h, double a, double b, double xmin, double xmax) const;
double HI2(TH1D * const h, double a, double b) const;
double Integral(TH2D* source, int nbin) const;
-
//Creating profile from histogramm
void CreateHist(char *name, char *nameAbs, char *nameRel, char *atitleF, char *atitleB,char *rtitleF, char *rtitleB,TH2D* sourceHist);
void SetGraphics() const;
void SetErrors(TH2D* source, const char *name, double ptd, TH2D* nb);
void SetErrors(TH2D* source, const char *name, double ptd, TProfile* nb);
- private:
+
+ private:
char* fSx; // Title of x axis
char* fSy; // Title of y axis
double fxFitMin; // FitMin minimum of fit baundary
double fbabsError; // b absolut error
double fXi2abs; // chi square absolut
- ClassDef(AliLRCAnalysis,0); // macro for rootcint
+ ClassDef(AliLRCAnalysis,0) // macro for rootcint
};
#endif
/* $Id$ */
#include "AliLRCFit.h"
+#include "TH1D.h"
+#include "TMath.h"
+class TH1D;
ClassImp(AliLRCFit)
//Empty constructor
}
-AliLRCFit::AliLRCFit(TH1D* h):fN (0), fTrueN(0), fNmin (0), fNmax(0), fS1(.0), fSz(.0), fSfz(.0), fSf(.0), fSf2(.0), fhi2(.0), fw(.0), fz(.0), f(.0), fnum(.0), fdf(.0), fdelta(.0), fa(.0), fb(.0), fda(.0), fdb(.0), fda1(.0), fdb1(.0), fxmin(.0), fxmax(.0){
+AliLRCFit::AliLRCFit(TH1D * const h):fN (0), fTrueN(0), fNmin (0), fNmax(0), fS1(.0), fSz(.0), fSfz(.0), fSf(.0), fSf2(.0), fhi2(.0), fw(.0), fz(.0), f(.0), fnum(.0), fdf(.0), fdelta(.0), fa(.0), fb(.0), fda(.0), fdb(.0), fda1(.0), fdb1(.0), fxmin(.0), fxmax(.0){
//Constructor make fit of 1d histogramm
fxmin = h->GetXaxis()->GetXmin();
fxmax = h->GetXaxis()->GetXmax();
for(int i=1; i<=fN; i++){
f=i*fdf-fdf/2+fnum;
fw = h->GetBinError(i);
- if(fw!=0){
+ if(fw){
fz = h->GetBinContent(i);
fS1 = fS1 + 1/(fw*fw);
fSz = fSz + fz/(fw*fw);
fdelta = fS1*fSf2 - fSf*fSf;
fb = (fS1*fSfz - fSf*fSz)/fdelta;
fa = ((fSf2-fSf)*fSz - (fSf-fS1)*fSfz)/fdelta;
- fda = sqrt((fS1+fSf2-2*fSf)/fdelta);
- fdb = sqrt(fS1/fdelta);
+ fda = TMath::Sqrt((fS1+fSf2-2*fSf)/fdelta);
+ fdb = TMath::Sqrt(fS1/fdelta);
fdb1 = 0;
fda1 = 0;
f = h->GetXaxis()->GetXmin();
for(int i=1; i<=fN; i++){
f=i*fdf-fdf/2+fnum;
fw = h->GetBinError(i);
- if(fw!=0){
+ if(fw){
fz = h->GetBinContent(i);
fdb1 = fdb1 + ((fS1*f - fSf)*(fS1*f - fSf)/(fw*fw)) * ((fz-fa-fb*(f-1))*(fz-fa-fb*(f-1))/(fw*fw));
fda1 = fda1 + (((fSf2-fSf)-(fSf-fS1)*f)*((fSf2-fSf)-(fSf-fS1)*f)/(fw*fw)) * ((fz-fa+fb-fb*f)*(fz-fa+fb-fb*f)/(fw*fw));
fTrueN++;
}
}
- fdb1 = sqrt(fdb1/(fdelta*fdelta));
- fda1 = sqrt(fda1/(fdelta*fdelta));
- fhi2 = fhi2 / (fTrueN-2);
+ fdb1 = TMath::Sqrt(fdb1/(fdelta*fdelta));
+ fda1 = TMath::Sqrt(fda1/(fdelta*fdelta));
+ fhi2 = fTrueN > 2 ? fhi2 / (fTrueN-2) : -1;
}
-AliLRCFit::AliLRCFit(TH1D *h, double xmin, double xmax):fN (0), fTrueN(0), fNmin (0), fNmax(0), fS1(.0), fSz(.0), fSfz(.0), fSf(.0), fSf2(.0), fhi2(.0), fw(.0), fz(.0), f(.0), fnum(.0), fdf(.0), fdelta(.0), fa(.0), fb(.0), fda(.0), fdb(.0), fda1(.0), fdb1(.0), fxmin(.0), fxmax(.0){
+AliLRCFit::AliLRCFit(TH1D * const h, double xmin, double xmax):fN (0), fTrueN(0), fNmin (0), fNmax(0), fS1(.0), fSz(.0), fSfz(.0), fSf(.0), fSf2(.0), fhi2(.0), fw(.0), fz(.0), f(.0), fnum(.0), fdf(.0), fdelta(.0), fa(.0), fb(.0), fda(.0), fdb(.0), fda1(.0), fdb1(.0), fxmin(.0), fxmax(.0){
//Constructor make fit of 1d histogramm between xmin and xmax
fxmin = xmin;
fxmax = xmax;
for(int i=fNmin; i<=fNmax; i++){
f=i*fdf-fdf/2+fnum;
fw = h->GetBinError(i);
- if(fw!=0){
+ if(fw){
fz = h->GetBinContent(i);
fS1 = fS1 + 1/(fw*fw);
fSz = fSz + fz/(fw*fw);
fdelta = fS1*fSf2 - fSf*fSf;
fb = (fS1*fSfz - fSf*fSz)/fdelta;
fa = ((fSf2-fSf)*fSz - (fSf-fS1)*fSfz)/fdelta;
- fda = sqrt((fS1+fSf2-2*fSf)/fdelta);
- fdb = sqrt(fS1/fdelta);
+ fda = TMath::Sqrt((fS1+fSf2-2*fSf)/fdelta);
+ fdb = TMath::Sqrt(fS1/fdelta);
fdb1 = 0;
fda1 = 0;
for(int i=fNmin; i<=fNmax; i++){
f=i*fdf-fdf/2+fnum;
fw = h->GetBinError(i);
- if(fw!=0){
+ if(fw){
fz = h->GetBinContent(i);
fdb1 = fdb1 + ((fS1*f - fSf)*(fS1*f - fSf)/(fw*fw)) * ((fz-fa-fb*(f-1))*(fz-fa-fb*(f-1))/(fw*fw));
fda1 = fda1 + (((fSf2-fSf)-(fSf-fS1)*f)*((fSf2-fSf)-(fSf-fS1)*f)/(fw*fw)) * ((fz-fa+fb-fb*f)*(fz-fa+fb-fb*f)/(fw*fw));
fTrueN++;
}
}
- fdb1 = sqrt(fdb1/(fdelta*fdelta));
- fda1 = sqrt(fda1/(fdelta*fdelta));
- fhi2 = fhi2 / (fTrueN-2);
+ fdb1 = TMath::Sqrt(fdb1/(fdelta*fdelta));
+ fda1 = TMath::Sqrt(fda1/(fdelta*fdelta));
+ fhi2 = fTrueN > 2 ? fhi2 / (fTrueN-2) : -1;
}
AliLRCFit::~AliLRCFit() {
-//-------------------------------------------------------------------------\r
-// Description: \r
-// This class include into LRC library for Long-Range Correlation analysis\r
-// it makes fit of the 1d histogramm\r
-// calculates ax+b coefficients with error and hi square\r
-// Origin: Petr Naumenko, SPbSU-CERN, Petr.Naoumenko@cern.ch\r
-//-------------------------------------------------------------------------\r
-\r
-#ifndef ALILRCFIT_H\r
-#define ALILRCFIT_H\r
-\r
-/* See cxx source for full Copyright notice */\r
-\r
-\r
-/* $Id$ */\r
-\r
-#include "TFile.h"\r
-#include "TH1D.h"\r
-#include "TH2D.h"\r
-#include "math.h"\r
-#include "TProfile.h"\r
-#include "TPaveText.h"\r
-#include "TF1.h"\r
-#include "TStyle.h"\r
-\r
-\r
-class TH1D;\r
-class TH2D;\r
-class TFile;\r
-\r
-\r
-class AliLRCFit{\r
- public:\r
- AliLRCFit();\r
- AliLRCFit(TH1D* h);\r
- AliLRCFit(TH1D* h, double xmin, double xmax);\r
- virtual ~AliLRCFit();\r
- double Geta() const;\r
- double Getb() const;\r
- double Getda() const;\r
- double Getdb() const;\r
- double Getda1() const;\r
- double Getdb1() const;\r
- double Gethi2() const;\r
- double Getf() const;\r
- double Getxmin() const;\r
- double Getxmax() const;\r
- int GetN() const;\r
-\r
- private:\r
- int fN; //Number of bins\r
- int fTrueN; //Number of bins between xmin and xmax\r
- int fNmin; //xmin bin\r
- int fNmax; //xmax bin \r
- double fS1; // work wariable\r
- double fSz; // work wariable\r
- double fSfz; // work wariable\r
- double fSf; // work wariable\r
- double fSf2; // work wariable\r
- double fhi2; // hi square\r
- double fw; // work wariable\r
- double fz; // work wariable\r
- double f; // work wariable\r
- double fnum; // work wariable\r
- double fdf; // work wariable\r
- double fdelta; // work wariable\r
- double fa; // coefficient a of ax+b\r
- double fb; // coefficient b of ax+b\r
- double fda; //error of coefficient a of ax+b\r
- double fdb; //error of coefficient b of ax+b\r
- double fda1; // work wariable\r
- double fdb1; // work wariable\r
- double fxmin; // xmin\r
- double fxmax; // xmax\r
- ClassDef(AliLRCFit,0) // macro for rootcint\r
-\r
-};\r
-\r
-#endif\r
-\r
+#ifndef ALILRCFIT_H
+#define ALILRCFIT_H
+
+//-------------------------------------------------------------------------
+// Description:
+// This class include into LRC library for Long-Range Correlation analysis
+// it makes fit of the 1d histogramm
+// calculates ax+b coefficients with error and hi square
+// Origin: Petr Naumenko, SPbSU-CERN, Petr.Naoumenko@cern.ch
+//-------------------------------------------------------------------------
+
+
+/* See cxx source for full Copyright notice */
+
+
+/* $Id$ */
+
+#include "TObject.h"
+class TH1D;
+
+class AliLRCFit {
+ public:
+ AliLRCFit();
+ AliLRCFit(TH1D * const h);
+ AliLRCFit(TH1D * const h, double xmin, double xmax);
+ virtual ~AliLRCFit();
+
+ double Geta() const;
+ double Getb() const;
+ double Getda() const;
+ double Getdb() const;
+ double Getda1() const;
+ double Getdb1() const;
+ double Gethi2() const;
+ double Getf() const;
+ double Getxmin() const;
+ double Getxmax() const;
+ int GetN() const;
+
+ private:
+ int fN; //Number of bins
+ int fTrueN; //Number of bins between xmin and xmax
+ int fNmin; //xmin bin
+ int fNmax; //xmax bin
+ double fS1; // work wariable
+ double fSz; // work wariable
+ double fSfz; // work wariable
+ double fSf; // work wariable
+ double fSf2; // work wariable
+ double fhi2; // hi square
+ double fw; // work wariable
+ double fz; // work wariable
+ double f; // work wariable
+ double fnum; // work wariable
+ double fdf; // work wariable
+ double fdelta; // work wariable
+ double fa; // coefficient a of ax+b
+ double fb; // coefficient b of ax+b
+ double fda; //error of coefficient a of ax+b
+ double fdb; //error of coefficient b of ax+b
+ double fda1; // work wariable
+ double fdb1; // work wariable
+ double fxmin; // xmin
+ double fxmax; // xmax
+
+ ClassDef(AliLRCFit,0) // macro for rootcint
+};
+
+#endif
+
#include "AliLRCNN.h"
#include "TFile.h"
+class TFile;
+class TH2D;
+class TList;
+
ClassImp(AliLRCNN)
SetErrors(sourceHist, profname);
}
-AliLRCNN::AliLRCNN(TList *LHist, char *histname, char *profname):AliLRCAnalysis() {
+AliLRCNN::AliLRCNN(TList * const LHist, char *histname, char *profname):AliLRCAnalysis() {
//Make NN from 2d histogramm from root file
SetGraphics();
TH2D* sourceHist = (TH2D*) LHist->FindObject(histname);
SetErrors(sourceHist, profname);
}
-void AliLRCNN::MakeHistogramm(TList *LHist, char *histname, char *profname) {
+void AliLRCNN::MakeHistogramm(TList * const LHist, char *histname, char *profname) {
//Make NN from 2d histogramm from root file
SetGraphics();
TH2D* sourceHist = (TH2D*) LHist->FindObject(histname);
-//-------------------------------------------------------------------------\r
-// Description: \r
-// This class include into LRC library for Long-Range Correlation analysis\r
-// it is the NN class\r
-// calculates NN correlations for abs and rel var\r
-// Origin: Petr Naumenko, SPbSU-CERN, Petr.Naoumenko@cern.ch\r
-//-------------------------------------------------------------------------\r
-\r
-#ifndef ALILRCNN_H\r
-#define ALILRCNN_H\r
-/* See cxx source for full Copyright notice */\r
-\r
-\r
-/* $Id$ */\r
-\r
-//-------------------------------------------------------------------------\r
-// LRC library for Long-Range Correlation analysis\r
-//\r
-// Origin: Petr Naumenko, SPbSU-CERN, Petr.Naoumenko@cern.ch\r
-//-------------------------------------------------------------------------\r
-#include "TProfile.h"\r
-#include "TH1D.h"\r
-#include "TH2D.h"\r
-#include "TPaveText.h"\r
-#include "TF1.h"\r
-#include "math.h"\r
-#include "TStyle.h"\r
+//-------------------------------------------------------------------------
+// Description:
+// This class include into LRC library for Long-Range Correlation analysis
+// it is the NN class
+// calculates NN correlations for abs and rel var
+// Origin: Petr Naumenko, SPbSU-CERN, Petr.Naoumenko@cern.ch
+//-------------------------------------------------------------------------
+
+#ifndef ALILRCNN_H
+#define ALILRCNN_H
+/* See cxx source for full Copyright notice */
+
+
+/* $Id$ */
+
+//-------------------------------------------------------------------------
+// LRC library for Long-Range Correlation analysis
+//
+// Origin: Petr Naumenko, SPbSU-CERN, Petr.Naoumenko@cern.ch
+//-------------------------------------------------------------------------
#include "AliLRCAnalysis.h"
-#include "TList.h"\r
-\r
-\r
-class TH1D;\r
-class TH2D;\r
-class TFile;
-class TList;\r
-\r
-class AliLRCNN : public AliLRCAnalysis{\r
- public:\r
- AliLRCNN();\r
- virtual ~AliLRCNN();\r
- AliLRCNN(char *name, TH2D* sourceHist);\r
+
+class TList;
+class TH2D;
+class AliLRCAnalysis;
+
+class AliLRCNN : public AliLRCAnalysis{
+ public:
+ AliLRCNN();
+ virtual ~AliLRCNN();
+ AliLRCNN(char *name, TH2D* sourceHist);
AliLRCNN(char *fileHistname, char *histname, char *profname);
- AliLRCNN(TList *LHist, char *histname, char *profname);\r
- void MakeHistogramm(char *name, TH2D* sourceHist);\r
+ AliLRCNN(TList *LHist, char *histname, char *profname);
+ void MakeHistogramm(char *name, TH2D* sourceHist);
void MakeHistogramm(char *fileHistname, char *histname, char *profname);
- void MakeHistogramm(TList *LHist, char *histname, char *profname);\r
- ClassDef(AliLRCNN,0) // macro for rootcint\r
-};\r
-\r
-\r
-\r
-#endif\r
-\r
+ void MakeHistogramm(TList *LHist, char *histname, char *profname);
+ ClassDef(AliLRCNN,0) // macro for rootcint
+};
+
+
+
+#endif
+
#include "AliLRCPtN.h"
#include "TFile.h"
+#include "TList.h"
+#include "TProfile.h"
+#include "TMath.h"
+
+class TFile;
+class TProfile;
+class TList;
+class TH2D;
ClassImp(AliLRCPtN)
AliLRCPtN::AliLRCPtN(char *name, TH2D* sourceHist, double ptd, TH2D* nb):AliLRCAnalysis() {
//Make PtN form 2d histogram
SetGraphics();
- CreateHist(name, (char*)"PtN_abs", (char*)"PtN_rel", (char*)"Pt_{F}", (char*)"<n_{B}>_{Pt_{F}}", (char*)"#frac{Pt_{F}}{<Pt_{F}>}", (char*)"#frac{<n_{B}>_{Pt_{F}}}{<n_{B}>}", sourceHist);
+ CreateHist(name, (char*)"PtN_abs", (char*)"PtN_rel", (char*)"n_{F}", (char*)"<Pt_{B}>_{n_{F}}", (char*)"#frac{n_{F}}{<n_{F}>}", (char*)"#frac{<Pt_{B}>_{n_{F}}}{<Pt_{B}>}", sourceHist);
SetErrors(sourceHist, name, ptd, nb);
}
fileHist = new TFile(fileHistname);
TH2D* sourceHist = (TH2D*) fileHist->Get(histname);
- CreateHist(profname, (char*)"PtN_abs", (char*)"PtN_rel", (char*)"Pt_{F}", (char*)"<n_{B}>_{Pt_{F}}", (char*)"#frac{Pt_{F}}{<Pt_{F}>}", (char*)"#frac{<n_{B}>_{Pt_{F}}}{<n_{B}>}", sourceHist);
+ CreateHist(profname, (char*)"PtN_abs", (char*)"PtN_rel", (char*)"n_{F}", (char*)"<Pt_{B}>_{n_{F}}", (char*)"#frac{n_{F}}{<n_{F}>}", (char*)"#frac{<Pt_{B}>_{n_{F}}}{<Pt_{B}>}", sourceHist);
TProfile* nbP = (TProfile*) fileHist->Get(errhistname);
SetErrors(sourceHist, profname, ptd, nbP);
}
-
-AliLRCPtN::AliLRCPtN(TList *LHist, char *histname, char *profname, char *ptdname, char *errhistname):AliLRCAnalysis() {
+//const
+AliLRCPtN::AliLRCPtN(TList * const LHist, char *histname, char *profname, char *ptdname, char *errhistname):AliLRCAnalysis() {
//Make PtN form 2d histogram from root file
SetGraphics();
TH2D* sourceHist = (TH2D*) LHist->FindObject(histname);
- CreateHist(profname, (char*)"PtN_abs", (char*)"PtN_rel", (char*)"Pt_{F}", (char*)"<n_{B}>_{Pt_{F}}", (char*)"#frac{Pt_{F}}{<Pt_{F}>}", (char*)"#frac{<n_{B}>_{Pt_{F}}}{<n_{B}>}", sourceHist);
+ CreateHist(profname, (char*)"PtN_abs", (char*)"PtN_rel", (char*)"n_{F}", (char*)"<Pt_{B}>_{n_{F}}", (char*)"#frac{n_{F}}{<n_{F}>}", (char*)"#frac{<Pt_{B}>_{n_{F}}}{<Pt_{B}>}", sourceHist);
TProfile* nbP = (TProfile*) LHist->FindObject(errhistname);
TProfile *dPtB = (TProfile*) LHist->FindObject(ptdname);
- double dptb=dPtB->GetBinError(1)*sqrt(dPtB->GetBinEntries(1));
+ double dptb=dPtB->GetBinError(1)*TMath::Sqrt(dPtB->GetBinEntries(1));
SetErrors(sourceHist, profname, dptb, nbP);
}
AliLRCPtN::AliLRCPtN(char *name, TH2D* sourceHist, double ptd, TProfile* nb):AliLRCAnalysis() {
//Make PtN form 2d histogram
SetGraphics();
- CreateHist(name, (char*)"PtN_abs", (char*)"PtN_rel", (char*)"Pt_{F}", (char*)"<n_{B}>_{Pt_{F}}", (char*)"#frac{Pt_{F}}{<Pt_{F}>}", (char*)"#frac{<n_{B}>_{Pt_{F}}}{<n_{B}>}", sourceHist);
+ CreateHist(name, (char*)"PtN_abs", (char*)"PtN_rel", (char*)"n_{F}", (char*)"<Pt_{B}>_{n_{F}}", (char*)"#frac{n_{F}}{<n_{F}>}", (char*)"#frac{<Pt_{B}>_{n_{F}}}{<Pt_{B}>}", sourceHist);
SetErrors(sourceHist, name, ptd, nb);
}
void AliLRCPtN::MakeHistogramm(char *name, TH2D* sourceHist, double ptd, TH2D* nb) {
//Make PtN form 2d histogram
SetGraphics();
- CreateHist(name, (char*)"PtN_abs", (char*)"PtN_rel", (char*)"Pt_{F}", (char*)"<n_{B}>_{Pt_{F}}", (char*)"#frac{Pt_{F}}{<Pt_{F}>}", (char*)"#frac{<n_{B}>_{Pt_{F}}}{<n_{B}>}", sourceHist);
+ CreateHist(name, (char*)"PtN_abs", (char*)"PtN_rel", (char*)"n_{F}", (char*)"<Pt_{B}>_{n_{F}}", (char*)"#frac{n_{F}}{<n_{F}>}", (char*)"#frac{<Pt_{B}>_{n_{F}}}{<Pt_{B}>}", sourceHist);
SetErrors(sourceHist, name, ptd, nb);
}
void AliLRCPtN::MakeHistogramm(char *name, TH2D* sourceHist, double ptd, TProfile* nb) {
//Make PtN form 2d histogram
SetGraphics();
- CreateHist(name, (char*)"PtN_abs", (char*)"PtN_rel", (char*)"Pt_{F}", (char*)"<n_{B}>_{Pt_{F}}", (char*)"#frac{Pt_{F}}{<Pt_{F}>}", (char*)"#frac{<n_{B}>_{Pt_{F}}}{<n_{B}>}", sourceHist);
+ CreateHist(name, (char*)"PtN_abs", (char*)"PtN_rel", (char*)"n_{F}", (char*)"<Pt_{B}>_{n_{F}}", (char*)"#frac{n_{F}}{<n_{F}>}", (char*)"#frac{<Pt_{B}>_{n_{F}}}{<Pt_{B}>}", sourceHist);
SetErrors(sourceHist, name, ptd, nb);
}
SetGraphics();
fileHist = new TFile(fileHistname);
TH2D* sourceHist = (TH2D*) fileHist->Get(histname);
- CreateHist(profname, (char*)"PtN_abs", (char*)"PtN_rel", (char*)"Pt_{F}", (char*)"<n_{B}>_{Pt_{F}}", (char*)"#frac{Pt_{F}}{<Pt_{F}>}", (char*)"#frac{<n_{B}>_{Pt_{F}}}{<n_{B}>}", sourceHist);
+ CreateHist(profname, (char*)"PtN_abs", (char*)"PtN_rel", (char*)"n_{F}", (char*)"<Pt_{B}>_{n_{F}}", (char*)"#frac{n_{F}}{<n_{F}>}", (char*)"#frac{<Pt_{B}>_{n_{F}}}{<Pt_{B}>}", sourceHist);
TProfile* nbP = (TProfile*) fileHist->Get(errhistname);
SetErrors(sourceHist, profname, ptd, nbP);
}
-void AliLRCPtN::MakeHistogramm(TList *LHist, char *histname, char *profname, char *ptdname, char *errhistname) {
+void AliLRCPtN::MakeHistogramm(TList * const LHist, char *histname, char *profname, char *ptdname, char *errhistname) {
//Make PtN form 2d histogram from root file
SetGraphics();
TH2D* sourceHist = (TH2D*) LHist->FindObject(histname);
- CreateHist(profname, (char*)"PtN_abs", (char*)"PtN_rel", (char*)"Pt_{F}", (char*)"<n_{B}>_{Pt_{F}}", (char*)"#frac{Pt_{F}}{<Pt_{F}>}", (char*)"#frac{<n_{B}>_{Pt_{F}}}{<n_{B}>}", sourceHist);
+ CreateHist(profname, (char*)"PtN_abs", (char*)"PtN_rel", (char*)"n_{F}", (char*)"<Pt_{B}>_{n_{F}}", (char*)"#frac{n_{F}}{<n_{F}>}", (char*)"#frac{<Pt_{B}>_{n_{F}}}{<Pt_{B}>}", sourceHist);
TProfile* nbP = (TProfile*) LHist->FindObject(errhistname);
TProfile *dPtB = (TProfile*) LHist->FindObject(ptdname);
- double dptb=dPtB->GetBinError(1)*sqrt(dPtB->GetBinEntries(1));
+ double dptb=dPtB->GetBinError(1)*TMath::Sqrt(dPtB->GetBinEntries(1));
SetErrors(sourceHist, profname, dptb, nbP);
}
-//-------------------------------------------------------------------------\r
-// Description: \r
-// This class include into LRC library for Long-Range Correlation analysis\r
-// it is the PtN class\r
-// calculates PtN correlations for abs and rel var\r
-// Origin: Petr Naumenko, SPbSU-CERN, Petr.Naoumenko@cern.ch\r
-//-------------------------------------------------------------------------\r
-\r
-#ifndef ALILRCPTN_H\r
-#define ALILRCPTN_H\r
-/* See cxx source for full Copyright notice */\r
-\r
-\r
-/* $Id$ */\r
-\r
-//-------------------------------------------------------------------------\r
-// LRC library for Long-Range Correlation analysis\r
-//\r
-// Origin: Petr Naumenko, SPbSU-CERN, Petr.Naoumenko@cern.ch\r
-//-------------------------------------------------------------------------\r
-#include "AliLRCAnalysis.h"\r
-#include "TH1D.h"\r
-#include "TH2D.h"\r
-#include "TFile.h"\r
-#include "TProfile.h"
-#include "TList.h"\r
-#include "math.h"\r
-#include <typeinfo>\r
-\r
-class TProfile;\r
+//-------------------------------------------------------------------------
+// Description:
+// This class include into LRC library for Long-Range Correlation analysis
+// it is the PtN class
+// calculates PtN correlations for abs and rel var
+// Origin: Petr Naumenko, SPbSU-CERN, Petr.Naoumenko@cern.ch
+//-------------------------------------------------------------------------
+
+#ifndef ALILRCPTN_H
+#define ALILRCPTN_H
+/* See cxx source for full Copyright notice */
+
+
+/* $Id$ */
+
+//-------------------------------------------------------------------------
+// LRC library for Long-Range Correlation analysis
+//
+// Origin: Petr Naumenko, SPbSU-CERN, Petr.Naoumenko@cern.ch
+//-------------------------------------------------------------------------
+
+#include "AliLRCAnalysis.h"
+
+class TProfile;
+class TList;
class TH2D;
-class TList;\r
-\r
-\r
-class AliLRCPtN : public AliLRCAnalysis {\r
- public:\r
- AliLRCPtN();\r
- ~AliLRCPtN();\r
+
+
+class AliLRCPtN : public AliLRCAnalysis {
+ public:
+ AliLRCPtN();
+ ~AliLRCPtN();
AliLRCPtN(char *name, TH2D* sourceHist, double ptd, TProfile* nb);
- AliLRCPtN(char *name, TH2D* sourceHist, double ptd, TH2D* nb);\r
+ AliLRCPtN(char *name, TH2D* sourceHist, double ptd, TH2D* nb);
AliLRCPtN(char *fileHistname, char *histname, char *profname, double ptd, char *errhistname);
- AliLRCPtN(TList *LHist, char *histname, char *profname, char *ptdname, char *errhistname); \r
- void MakeHistogramm(char *name, TH2D* sourceHist, double ptd, TH2D* nb);\r
- void MakeHistogramm(char *name, TH2D* sourceHist, double ptd, TProfile* nb);\r
+ AliLRCPtN(TList *LHist, char *histname, char *profname, char *ptdname, char *errhistname);
+ void MakeHistogramm(char *name, TH2D* sourceHist, double ptd, TH2D* nb);
+ void MakeHistogramm(char *name, TH2D* sourceHist, double ptd, TProfile* nb);
void MakeHistogramm(char *fileHistname, char *histname, char *profname, double ptd, char *errhistname);
- void MakeHistogramm(TList *LHist, char *histname, char *profname, char *ptdname, char *errhistname); \r
- ClassDef(AliLRCPtN,0) // macro for rootcint\r
-};\r
-\r
-\r
-\r
-#endif\r
-\r
+ void MakeHistogramm(TList *LHist, char *histname, char *profname, char *ptdname, char *errhistname);
+ ClassDef(AliLRCPtN,0) // macro for rootcint
+};
+
+
+
+#endif
+
//-------------------------------------------------------------------------
#include "AliLRCPtPt.h"
+#include "TFile.h"
+#include "TList.h"
+#include "TProfile.h"
+#include "TMath.h"
+
+class TFile;
+class TProfile;
+class TList;
+class TH2D;
ClassImp(AliLRCPtPt)
SetErrors(sourceHist, profname, ptd, nbP);
}
-AliLRCPtPt::AliLRCPtPt(TList *LHist, char *histname, char *profname, char *ptdname, char *errhistname):AliLRCAnalysis() {
+AliLRCPtPt::AliLRCPtPt(TList * const LHist, char *histname, char *profname, char *ptdname, char *errhistname):AliLRCAnalysis() {
//Make PtPt form 2d histogram from root file
SetGraphics();
TH2D* sourceHist = (TH2D*) LHist->FindObject(histname);
CreateHist(profname, (char*)"PtPt_abs", (char*)"PtPt_rel", (char*)"Pt_{F}", (char*)"<Pt_{B}>_{Pt_{F}}", (char*)"#frac{Pt_{F}}{<Pt_{F}>}", (char*)"#frac{<Pt_{B}>_{Pt_{F}}}{<Pt_{B}>}", sourceHist);
TProfile* nbP = (TProfile*) LHist->FindObject(errhistname);
TProfile *dPtB = (TProfile*) LHist->FindObject(ptdname);
- double dptb=dPtB->GetBinError(1)*sqrt(dPtB->GetBinEntries(1));
+ double dptb=dPtB->GetBinError(1)*TMath::Sqrt(dPtB->GetBinEntries(1));
SetErrors(sourceHist, profname, dptb, nbP);
}
SetErrors(sourceHist, profname, ptd, nbP);
}
-void AliLRCPtPt::MakeHistogramm(TList *LHist, char *histname, char *profname, char *ptdname, char *errhistname) {
+void AliLRCPtPt::MakeHistogramm(TList * const LHist, char *histname, char *profname, char *ptdname, char *errhistname) {
//Make PtPt form 2d histogram from root file
SetGraphics();
TH2D* sourceHist = (TH2D*) LHist->FindObject(histname);
CreateHist(profname, (char*)"PtPt_abs", (char*)"PtPt_rel", (char*)"Pt_{F}", (char*)"<Pt_{B}>_{Pt_{F}}", (char*)"#frac{Pt_{F}}{<Pt_{F}>}", (char*)"#frac{<Pt_{B}>_{Pt_{F}}}{<Pt_{B}>}", sourceHist);
TProfile* nbP = (TProfile*) LHist->FindObject(errhistname);
TProfile *dPtB = (TProfile*) LHist->FindObject(ptdname);
- double dptb=dPtB->GetBinError(1)*sqrt(dPtB->GetBinEntries(1));
+ double dptb=dPtB->GetBinError(1)*TMath::Sqrt(dPtB->GetBinEntries(1));
SetErrors(sourceHist, profname, dptb, nbP);
}
-//-------------------------------------------------------------------------\r
-// Description: \r
-// This class include into LRC library for Long-Range Correlation analysis\r
-// it is the PtPt class\r
-// calculates PtPt correlations for abs and rel var\r
-// Origin: Petr Naumenko, SPbSU-CERN, Petr.Naoumenko@cern.ch\r
-//-------------------------------------------------------------------------\r
-#ifndef ALILRCPTPT_H\r
-#define ALILRCPTPT_H\r
-/* See cxx source for full Copyright notice */\r
-\r
-\r
-/* $Id$ */\r
-\r
-//-------------------------------------------------------------------------\r
-// LRC library for Long-Range Correlation analysis\r
-//\r
-// Origin: Petr Naumenko, SPbSU-CERN, Petr.Naoumenko@cern.ch\r
-//-------------------------------------------------------------------------\r
-\r
-#include "AliLRCAnalysis.h"\r
-#include "TH1D.h"\r
-#include "TH2D.h"\r
-#include "TFile.h"\r
-#include "TProfile.h"
-#include "TList.h"\r
-#include "math.h"\r
-#include <typeinfo>\r
-\r
+//-------------------------------------------------------------------------
+// Description:
+// This class include into LRC library for Long-Range Correlation analysis
+// it is the PtPt class
+// calculates PtPt correlations for abs and rel var
+// Origin: Petr Naumenko, SPbSU-CERN, Petr.Naoumenko@cern.ch
+//-------------------------------------------------------------------------
+#ifndef ALILRCPTPT_H
+#define ALILRCPTPT_H
+/* See cxx source for full Copyright notice */
+
+
+/* $Id$ */
+
+//-------------------------------------------------------------------------
+// LRC library for Long-Range Correlation analysis
+//
+// Origin: Petr Naumenko, SPbSU-CERN, Petr.Naoumenko@cern.ch
+//-------------------------------------------------------------------------
+
+#include "AliLRCAnalysis.h"
+
+class TProfile;
+class TList;
class TH2D;
-class TList;\r
-\r
-class AliLRCPtPt:public AliLRCAnalysis\r
-{\r
- public:\r
- AliLRCPtPt();\r
- ~AliLRCPtPt();\r
- AliLRCPtPt(char *name, TH2D* sourceHist, double ptd, TH2D* nb);\r
- AliLRCPtPt(char *name, TH2D* sourceHist, double ptd, TProfile* nb);\r
+
+class AliLRCPtPt:public AliLRCAnalysis
+{
+ public:
+ AliLRCPtPt();
+ ~AliLRCPtPt();
+ AliLRCPtPt(char *name, TH2D* sourceHist, double ptd, TH2D* nb);
+ AliLRCPtPt(char *name, TH2D* sourceHist, double ptd, TProfile* nb);
AliLRCPtPt(char *fileHistname, char *histname, char *profname, double ptd, char *errhistname);
- AliLRCPtPt(TList *LHist, char *histname, char *profname, char *ptdname, char *errhistname);\r
- void MakeHistogramm(char *name, TH2D* sourceHist, double ptd, TH2D* nb);\r
- void MakeHistogramm(char *name, TH2D* sourceHist, double ptd, TProfile* nb);\r
+ AliLRCPtPt(TList *LHist, char *histname, char *profname, char *ptdname, char *errhistname);
+ void MakeHistogramm(char *name, TH2D* sourceHist, double ptd, TH2D* nb);
+ void MakeHistogramm(char *name, TH2D* sourceHist, double ptd, TProfile* nb);
void MakeHistogramm(char *fileHistname, char *histname, char *profname, double ptd, char *errhistname) ;
- void MakeHistogramm(TList *LHist, char *histname, char *profname, char *ptdname, char *errhistname) ; \r
- ClassDef(AliLRCPtPt,0) // macro for rootcint\r
-\r
-};\r
-\r
-#endif\r
-\r
+ void MakeHistogramm(TList *LHist, char *histname, char *profname, char *ptdname, char *errhistname) ;
+ ClassDef(AliLRCPtPt,0) // macro for rootcint
+
+};
+
+#endif
+