Fixing warning and adding new functionality (Andrey)
authorpchrist <pchrist@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 27 Mar 2010 09:11:38 +0000 (09:11 +0000)
committerpchrist <pchrist@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 27 Mar 2010 09:11:38 +0000 (09:11 +0000)
PWG2/EBYE/LRC/AliLRCAnalysis.cxx
PWG2/EBYE/LRC/AliLRCAnalysis.h
PWG2/EBYE/LRC/AliLRCFit.cxx
PWG2/EBYE/LRC/AliLRCFit.h
PWG2/EBYE/LRC/AliLRCNN.cxx
PWG2/EBYE/LRC/AliLRCNN.h
PWG2/EBYE/LRC/AliLRCPtN.cxx
PWG2/EBYE/LRC/AliLRCPtN.h
PWG2/EBYE/LRC/AliLRCPtPt.cxx
PWG2/EBYE/LRC/AliLRCPtPt.h

index 32b28e7..86df979 100755 (executable)
@@ -26,8 +26,8 @@
 //    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"
@@ -35,7 +35,7 @@
 #include "TH2D.h"
 #include "TPaveText.h"
 #include "TF1.h"
-//#include "math.h"
+#include "TMath.h"
 #include "TStyle.h"
 
 class gStyle;
@@ -47,11 +47,11 @@ ClassImp(AliLRCAnalysis)
  * 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
 }
 
@@ -71,9 +71,13 @@ AliLRCAnalysis& AliLRCAnalysis::operator= (const AliLRCAnalysis& a){
                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;
@@ -155,8 +159,8 @@ void AliLRCAnalysis::CreateHist(char *name, char *nameAbs, char *nameRel, char *
     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;
 
 }
@@ -164,7 +168,7 @@ void AliLRCAnalysis::CreateHist(char *name, char *nameAbs, char *nameRel, char *
 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;
@@ -280,13 +284,12 @@ void AliLRCAnalysis::Calculate(){
 }
 
 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 ) {
@@ -367,13 +370,12 @@ void AliLRCAnalysis::DrawAbsPure( const int * const mDrawArray, bool drawPaveLab
 }
 
 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 ) {
@@ -480,7 +482,7 @@ void AliLRCAnalysis::SetErrors(TH2D* source, const char *name){
                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());  
@@ -498,7 +500,7 @@ void AliLRCAnalysis::SetErrors(TH2D* source, const char *name, double ptd, TH2D*
                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());  
@@ -517,7 +519,7 @@ void AliLRCAnalysis::SetErrors(TH2D* source, const char *name, double ptd, TProf
                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());  
@@ -575,7 +577,7 @@ double AliLRCAnalysis::GetRoundValueErrorPrecision( double value, double error,
        
        if ( noError )
        {
-         cout << "After rounding: " << error << endl;
+               cout << "After rounding: " << error << endl;
                return error;
        }
        else
index 894217d..069dcb7 100755 (executable)
@@ -1,13 +1,14 @@
+#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;
@@ -34,7 +35,7 @@ class TF1;
 class math;
 class TStyle;
 
-class AliLRCAnalysis : public TObject{
+class AliLRCAnalysis{
  public:
   void DrawAbs();
   void DrawAbs( int * mas );
@@ -83,7 +84,6 @@ class AliLRCAnalysis : public TObject{
   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;
@@ -91,7 +91,8 @@ class AliLRCAnalysis : public TObject{
   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
@@ -107,7 +108,7 @@ class AliLRCAnalysis : public TObject{
   double fbabsError;   // b absolut error
   double fXi2abs;      // chi square absolut  
   
-  ClassDef(AliLRCAnalysis,0);                 // macro for rootcint
+  ClassDef(AliLRCAnalysis,0)                 // macro for rootcint
 };
 
 #endif
index 4c44a97..b190342 100755 (executable)
 /* $Id$ */
 
 #include "AliLRCFit.h"
+#include "TH1D.h"
+#include "TMath.h"
 
+class TH1D;
 
 ClassImp(AliLRCFit) 
 
@@ -31,7 +34,7 @@ AliLRCFit::AliLRCFit():fN (0), fTrueN(0), fNmin (0), fNmax(0), fS1(.0), fSz(.0),
 //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();
@@ -43,7 +46,7 @@ AliLRCFit::AliLRCFit(TH1D* h):fN (0), fTrueN(0), fNmin (0), fNmax(0), fS1(.0), f
     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);
@@ -55,15 +58,15 @@ AliLRCFit::AliLRCFit(TH1D* h):fN (0), fTrueN(0), fNmin (0), fNmax(0), fS1(.0), f
     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));
@@ -71,12 +74,12 @@ AliLRCFit::AliLRCFit(TH1D* h):fN (0), fTrueN(0), fNmin (0), fNmax(0), fS1(.0), f
             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;
@@ -88,7 +91,7 @@ AliLRCFit::AliLRCFit(TH1D *h, double xmin, double xmax):fN (0), fTrueN(0), fNmin
     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);
@@ -100,14 +103,14 @@ AliLRCFit::AliLRCFit(TH1D *h, double xmin, double xmax):fN (0), fTrueN(0), fNmin
     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));
@@ -115,9 +118,9 @@ AliLRCFit::AliLRCFit(TH1D *h, double xmin, double xmax):fN (0), fTrueN(0), fNmin
             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() {
index 3f69ad0..560d141 100755 (executable)
@@ -1,80 +1,70 @@
-//-------------------------------------------------------------------------\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
+
index e142350..45376d3 100755 (executable)
 #include "AliLRCNN.h"
 #include "TFile.h"
 
+class TFile;
+class TH2D;
+class TList;
+
 ClassImp(AliLRCNN) 
 
 
@@ -55,7 +59,7 @@ AliLRCNN::AliLRCNN(char *fileHistname, char *histname, char *profname):AliLRCAna
     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);
@@ -80,7 +84,7 @@ void AliLRCNN::MakeHistogramm(char *fileHistname, char *histname, char *profname
     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);
index f9be156..2d5ab47 100755 (executable)
@@ -1,53 +1,43 @@
-//-------------------------------------------------------------------------\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
+
index 9c27161..c322b73 100755 (executable)
 
 #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) 
 
@@ -48,7 +56,7 @@ AliLRCPtN::~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);
 }
 
@@ -58,41 +66,41 @@ AliLRCPtN::AliLRCPtN(char *fileHistname, char *histname, char *profname, double
     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);
 }
 
@@ -101,18 +109,18 @@ void AliLRCPtN::MakeHistogramm(char *fileHistname, char *histname, char *profnam
     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);
 }
index 60b4f79..cf02621 100755 (executable)
@@ -1,53 +1,47 @@
-//-------------------------------------------------------------------------\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
+
index f5e3bba..5acc9f8 100755 (executable)
 //-------------------------------------------------------------------------
 
 #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)
 
@@ -65,14 +74,14 @@ AliLRCPtPt::AliLRCPtPt(char *fileHistname, char *histname, char *profname, doubl
        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);
 }
 
@@ -101,13 +110,13 @@ void AliLRCPtPt::MakeHistogramm(char *fileHistname, char *histname, char *profna
        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);
 }
index c7b4534..0a74462 100755 (executable)
@@ -1,51 +1,45 @@
-//-------------------------------------------------------------------------\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
+