]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - HBTAN/AliHBTashbtCorrFctn.cxx
Updates in GRP Preprocessor (Ernesto)
[u/mrichter/AliRoot.git] / HBTAN / AliHBTashbtCorrFctn.cxx
index 5c83d9a3e174106f008089ecbdfc228022b5bce7..d85ec710b623a4b4a0feab3b31961766aea8f709 100644 (file)
@@ -1,6 +1,4 @@
 #include "AliHBTashbtCorrFctn.h"
-#include <TH1.h>
-#include <Riostream.h>
 
 ///////////////////////////////////////////////////////
 //                                                   //
 
 ClassImp(AliHBTashbtCorrFctn)
 
-AliHBTashbtCorrFctn::AliHBTashbtCorrFctn(const char* name, const char* title):
- AliHBTOnePairFctn1D(name,title),
-    fNumOut1(0x0),
-    fNumOut2(0x0),
-    fNumOut3(0x0),
-    fNumOut4(0x0),
-    fNumOut5(0x0),
-    fNumOut6(0x0),
-    fNumOut7(0x0),
-    fNumOut8(0x0),
-
-    fDenOut1(0x0),
-    fDenOut2(0x0),
-    fDenOut3(0x0),
-    fDenOut4(0x0),
-    fDenOut5(0x0),
-    fDenOut6(0x0),
-    fDenOut7(0x0),
-    fDenOut8(0x0),
-   
-    fRatOut1(0x0),
-    fRatOut2(0x0),
-    fRatOut3(0x0),
-    fRatOut4(0x0),
-    fRatOut5(0x0),
-    fRatOut6(0x0),
-    fRatOut7(0x0),
-    fRatOut8(0x0),
-   
-
-    fNumSide1(0x0),
-    fNumSide2(0x0),
-    fNumSide3(0x0),
-    fNumSide4(0x0),
-    fNumSide5(0x0),
-    fNumSide6(0x0),
-    fNumSide7(0x0),
-    fNumSide8(0x0),
-
-    fDenSide1(0x0),
-    fDenSide2(0x0),
-    fDenSide3(0x0),
-    fDenSide4(0x0),
-    fDenSide5(0x0),
-    fDenSide6(0x0),
-    fDenSide7(0x0),
-    fDenSide8(0x0),
-   
-    fRatSide1(0x0),
-    fRatSide2(0x0),
-    fRatSide3(0x0),
-    fRatSide4(0x0),
-    fRatSide5(0x0),
-    fRatSide6(0x0),
-    fRatSide7(0x0),
-    fRatSide8(0x0),
-   
-    fNumLong1(0x0),
-    fNumLong2(0x0),
-    fNumLong3(0x0),
-    fNumLong4(0x0),
-    fNumLong5(0x0),
-    fNumLong6(0x0),
-    fNumLong7(0x0),
-    fNumLong8(0x0),
-
-    fDenLong1(0x0),
-    fDenLong2(0x0),
-    fDenLong3(0x0),
-    fDenLong4(0x0),
-    fDenLong5(0x0),
-    fDenLong6(0x0),
-    fDenLong7(0x0),
-    fDenLong8(0x0),
-   
-    fRatLong1(0x0),
-    fRatLong2(0x0),
-    fRatLong3(0x0),
-    fRatLong4(0x0),
-    fRatLong5(0x0),
-    fRatLong6(0x0),
-    fRatLong7(0x0),
-    fRatLong8(0x0)
-
-
-{
-//ctor
-}
-/******************************************************************/
-AliHBTashbtCorrFctn::AliHBTashbtCorrFctn(const char* name, const char* title, Int_t nbins, Float_t maxXval, Float_t minXval):
- AliHBTOnePairFctn1D(name,title,nbins,maxXval,minXval),
-
-    fNumOut1(0x0),
-    fNumOut2(0x0),
-    fNumOut3(0x0),
-    fNumOut4(0x0),
-    fNumOut5(0x0),
-    fNumOut6(0x0),
-    fNumOut7(0x0),
-    fNumOut8(0x0),
-
-    fDenOut1(0x0),
-    fDenOut2(0x0),
-    fDenOut3(0x0),
-    fDenOut4(0x0),
-    fDenOut5(0x0),
-    fDenOut6(0x0),
-    fDenOut7(0x0),
-    fDenOut8(0x0),
-   
-    fRatOut1(0x0),
-    fRatOut2(0x0),
-    fRatOut3(0x0),
-    fRatOut4(0x0),
-    fRatOut5(0x0),
-    fRatOut6(0x0),
-    fRatOut7(0x0),
-    fRatOut8(0x0),
-   
-
-    fNumSide1(0x0),
-    fNumSide2(0x0),
-    fNumSide3(0x0),
-    fNumSide4(0x0),
-    fNumSide5(0x0),
-    fNumSide6(0x0),
-    fNumSide7(0x0),
-    fNumSide8(0x0),
-
-    fDenSide1(0x0),
-    fDenSide2(0x0),
-    fDenSide3(0x0),
-    fDenSide4(0x0),
-    fDenSide5(0x0),
-    fDenSide6(0x0),
-    fDenSide7(0x0),
-    fDenSide8(0x0),
-   
-    fRatSide1(0x0),
-    fRatSide2(0x0),
-    fRatSide3(0x0),
-    fRatSide4(0x0),
-    fRatSide5(0x0),
-    fRatSide6(0x0),
-    fRatSide7(0x0),
-    fRatSide8(0x0),
-   
-    fNumLong1(0x0),
-    fNumLong2(0x0),
-    fNumLong3(0x0),
-    fNumLong4(0x0),
-    fNumLong5(0x0),
-    fNumLong6(0x0),
-    fNumLong7(0x0),
-    fNumLong8(0x0),
-
-    fDenLong1(0x0),
-    fDenLong2(0x0),
-    fDenLong3(0x0),
-    fDenLong4(0x0),
-    fDenLong5(0x0),
-    fDenLong6(0x0),
-    fDenLong7(0x0),
-    fDenLong8(0x0),
-   
-    fRatLong1(0x0),
-    fRatLong2(0x0),
-    fRatLong3(0x0),
-    fRatLong4(0x0),
-    fRatLong5(0x0),
-    fRatLong6(0x0),
-    fRatLong7(0x0),
-    fRatLong8(0x0)
-
-
-{
-BuildHistos(nbins, maxXval, minXval);
-//ctor
-}
-
-/******************************************************************/
-AliHBTashbtCorrFctn::AliHBTashbtCorrFctn(const AliHBTashbtCorrFctn& in):
- AliHBTOnePairFctn1D(in),
-    fNumOut1(0x0),
-    fNumOut2(0x0),
-    fNumOut3(0x0),
-    fNumOut4(0x0),
-    fNumOut5(0x0),
-    fNumOut6(0x0),
-    fNumOut7(0x0),
-    fNumOut8(0x0),
-
-    fDenOut1(0x0),
-    fDenOut2(0x0),
-    fDenOut3(0x0),
-    fDenOut4(0x0),
-    fDenOut5(0x0),
-    fDenOut6(0x0),
-    fDenOut7(0x0),
-    fDenOut8(0x0),
-   
-    fRatOut1(0x0),
-    fRatOut2(0x0),
-    fRatOut3(0x0),
-    fRatOut4(0x0),
-    fRatOut5(0x0),
-    fRatOut6(0x0),
-    fRatOut7(0x0),
-    fRatOut8(0x0),
-   
-
-    fNumSide1(0x0),
-    fNumSide2(0x0),
-    fNumSide3(0x0),
-    fNumSide4(0x0),
-    fNumSide5(0x0),
-    fNumSide6(0x0),
-    fNumSide7(0x0),
-    fNumSide8(0x0),
-
-    fDenSide1(0x0),
-    fDenSide2(0x0),
-    fDenSide3(0x0),
-    fDenSide4(0x0),
-    fDenSide5(0x0),
-    fDenSide6(0x0),
-    fDenSide7(0x0),
-    fDenSide8(0x0),
-   
-    fRatSide1(0x0),
-    fRatSide2(0x0),
-    fRatSide3(0x0),
-    fRatSide4(0x0),
-    fRatSide5(0x0),
-    fRatSide6(0x0),
-    fRatSide7(0x0),
-    fRatSide8(0x0),
-   
-    fNumLong1(0x0),
-    fNumLong2(0x0),
-    fNumLong3(0x0),
-    fNumLong4(0x0),
-    fNumLong5(0x0),
-    fNumLong6(0x0),
-    fNumLong7(0x0),
-    fNumLong8(0x0),
-
-    fDenLong1(0x0),
-    fDenLong2(0x0),
-    fDenLong3(0x0),
-    fDenLong4(0x0),
-    fDenLong5(0x0),
-    fDenLong6(0x0),
-    fDenLong7(0x0),
-    fDenLong8(0x0),
-   
-    fRatLong1(0x0),
-    fRatLong2(0x0),
-    fRatLong3(0x0),
-    fRatLong4(0x0),
-    fRatLong5(0x0),
-    fRatLong6(0x0),
-    fRatLong7(0x0),
-    fRatLong8(0x0)
-  
-
- {
-//ctor
-}
-
-/******************************************************************/
-
-AliHBTashbtCorrFctn::~AliHBTashbtCorrFctn()
-{
- //dtor
-    delete fNumOut1;
-    delete fNumOut2;
-    delete fNumOut3;
-    delete fNumOut4;
-    delete fNumOut5;
-    delete fNumOut6;
-    delete fNumOut7;
-    delete fNumOut8;
-
-    delete fDenOut1;
-    delete fDenOut2;
-    delete fDenOut3;
-    delete fDenOut4;
-    delete fDenOut5;
-    delete fDenOut6;
-    delete fDenOut7;
-    delete fDenOut8;
-   
-    delete fRatOut1;
-    delete fRatOut2;
-    delete fRatOut3;
-    delete fRatOut4;
-    delete fRatOut5;
-    delete fRatOut6;
-    delete fRatOut7;
-    delete fRatOut8;
-   
-
-    delete fNumSide1;
-    delete fNumSide2;
-    delete fNumSide3;
-    delete fNumSide4;
-    delete fNumSide5;
-    delete fNumSide6;
-    delete fNumSide7;
-    delete fNumSide8;
-
-    delete fDenSide1;
-    delete fDenSide2;
-    delete fDenSide3;
-    delete fDenSide4;
-    delete fDenSide5;
-    delete fDenSide6;
-    delete fDenSide7;
-    delete fDenSide8;
-   
-    delete fRatSide1;
-    delete fRatSide2;
-    delete fRatSide3;
-    delete fRatSide4;
-    delete fRatSide5;
-    delete fRatSide6;
-    delete fRatSide7;
-    delete fRatSide8;
-   
-    delete fNumLong1;
-    delete fNumLong2;
-    delete fNumLong3;
-    delete fNumLong4;
-    delete fNumLong5;
-    delete fNumLong6;
-    delete fNumLong7;
-    delete fNumLong8;
-
-    delete fDenLong1;
-    delete fDenLong2;
-    delete fDenLong3;
-    delete fDenLong4;
-    delete fDenLong5;
-    delete fDenLong6;
-    delete fDenLong7;
-    delete fDenLong8;
-   
-    delete fRatLong1;
-    delete fRatLong2;
-    delete fRatLong3;
-    delete fRatLong4;
-    delete fRatLong5;
-    delete fRatLong6;
-    delete fRatLong7;
-    delete fRatLong8; 
-
-}
-
-/******************************************************************/
-void AliHBTashbtCorrFctn::WriteFunction()
-{
-//out    
-   Double_t out1scale = Scale(fNumOut1,fDenOut1);
-   cout <<"out1scale = "<<out1scale<<endl;
-   fRatOut1->Divide(fNumOut1,fDenOut1,out1scale);
-   
-   Double_t out2scale = Scale(fNumOut2,fDenOut2);
-   cout <<"out2scale = "<<out2scale<<endl;
-   fRatOut2->Divide(fNumOut2,fDenOut2,out2scale);
-   
-   Double_t out3scale = Scale(fNumOut3,fDenOut3);
-   cout <<"out3scale = "<<out3scale<<endl;
-   fRatOut3->Divide(fNumOut3,fDenOut3,out3scale);
-   
-   Double_t out4scale = Scale(fNumOut4,fDenOut4);
-   cout <<"out4scale = "<<out4scale<<endl;
-   fRatOut4->Divide(fNumOut4,fDenOut4,out4scale);
-   
-   Double_t out5scale = Scale(fNumOut5,fDenOut5);
-   cout <<"out5scale = "<<out5scale<<endl;
-   fRatOut5->Divide(fNumOut5,fDenOut5,out5scale);
-   
-   Double_t out6scale = Scale(fNumOut6,fDenOut6);
-   cout <<"out6scale = "<<out6scale<<endl;
-   fRatOut6->Divide(fNumOut6,fDenOut6,out6scale);
-   
-   Double_t out7scale = Scale(fNumOut7,fDenOut7);
-   cout <<"out7scale = "<<out7scale<<endl;
-   fRatOut7->Divide(fNumOut7,fDenOut7,out7scale);
-   
-   Double_t out8scale = Scale(fNumOut8,fDenOut8);
-   cout <<"out8scale = "<<out8scale<<endl;
-   fRatOut8->Divide(fNumOut8,fDenOut8,out8scale);
-
-   //side
-   Double_t side1scale = Scale(fNumSide1,fDenSide1);
-   cout <<"side1scale = "<<side1scale<<endl;
-   fRatSide1->Divide(fNumSide1,fDenSide1,side1scale);
-   
-   Double_t side2scale = Scale(fNumSide2,fDenSide2);
-   cout <<"side2scale = "<<side2scale<<endl;
-   fRatSide2->Divide(fNumSide2,fDenSide2,side2scale);
-   
-   Double_t side3scale = Scale(fNumSide3,fDenSide3);
-   cout <<"side3scale = "<<side3scale<<endl;
-   fRatSide3->Divide(fNumSide3,fDenSide3,side3scale);
-   
-   Double_t side4scale = Scale(fNumSide4,fDenSide4);
-   cout <<"side4scale = "<<side4scale<<endl;
-   fRatSide4->Divide(fNumSide4,fDenSide4,side4scale);
-   
-   Double_t side5scale = Scale(fNumSide5,fDenSide5);
-   cout <<"side5scale = "<<side5scale<<endl;
-   fRatSide5->Divide(fNumSide5,fDenSide5,side5scale);
-   
-   Double_t side6scale = Scale(fNumSide6,fDenSide6);
-   cout <<"side6scale = "<<side6scale<<endl;
-   fRatSide6->Divide(fNumSide6,fDenSide6,side6scale);
-   
-   Double_t side7scale = Scale(fNumSide7,fDenSide7);
-   cout <<"side7scale = "<<side7scale<<endl;
-   fRatSide7->Divide(fNumSide7,fDenSide7,side7scale);
-   
-   Double_t side8scale = Scale(fNumSide8,fDenSide8);
-   cout <<"side8scale = "<<side8scale<<endl;
-   fRatSide8->Divide(fNumSide8,fDenSide8,side8scale);
-   
-//long
-   Double_t long1scale = Scale(fNumLong1,fDenLong1);
-   cout <<"long1scale = "<<long1scale<<endl;
-   fRatLong1->Divide(fNumLong1,fDenLong1,long1scale);
-   
-   Double_t long2scale = Scale(fNumLong2,fDenLong2);
-   cout <<"long2scale = "<<long2scale<<endl;
-   fRatLong2->Divide(fNumLong2,fDenLong2,long2scale);
-   
-   Double_t long3scale = Scale(fNumLong3,fDenLong3);
-   cout <<"long3scale = "<<long3scale<<endl;
-   fRatLong3->Divide(fNumLong3,fDenLong3,long3scale);
-   
-   Double_t long4scale = Scale(fNumLong4,fDenLong4);
-   cout <<"long4scale = "<<long4scale<<endl;
-   fRatLong4->Divide(fNumLong4,fDenLong4,long4scale);
-   
-   Double_t long5scale = Scale(fNumLong5,fDenLong5);
-   cout <<"long5scale = "<<long5scale<<endl;
-   fRatLong5->Divide(fNumLong5,fDenLong5,long5scale);
-   
-   Double_t long6scale = Scale(fNumLong6,fDenLong6);
-   cout <<"long6scale = "<<long6scale<<endl;
-   fRatLong6->Divide(fNumLong6,fDenLong6,long6scale);
-   
-   Double_t long7scale = Scale(fNumLong7,fDenLong7);
-   cout <<"long7scale = "<<long7scale<<endl;
-   fRatLong7->Divide(fNumLong7,fDenLong7,long7scale);
-   
-   Double_t long8scale = Scale(fNumLong8,fDenLong8);
-   cout <<"long8scale = "<<long8scale<<endl;
-   fRatLong8->Divide(fNumLong8,fDenLong8,long8scale);
-
-        fNumOut1->Write();
-     fNumOut2->Write();
-     fNumOut3->Write();
-     fNumOut4->Write();
-     fNumOut5->Write();
-     fNumOut6->Write();
-     fNumOut7->Write();
-     fNumOut8->Write();
-
-     fDenOut1->Write();
-     fDenOut2->Write();
-     fDenOut3->Write();
-     fDenOut4->Write();
-     fDenOut5->Write();
-     fDenOut6->Write();
-     fDenOut7->Write();
-     fDenOut8->Write();
-   
-     fRatOut1->Write();
-     fRatOut2->Write();
-     fRatOut3->Write();
-     fRatOut4->Write();
-     fRatOut5->Write();
-     fRatOut6->Write();
-     fRatOut7->Write();
-     fRatOut8->Write();
-   
-
-     fNumSide1->Write();
-     fNumSide2->Write();
-     fNumSide3->Write();
-     fNumSide4->Write();
-     fNumSide5->Write();
-     fNumSide6->Write();
-     fNumSide7->Write();
-     fNumSide8->Write();
-
-     fDenSide1->Write();
-     fDenSide2->Write();
-     fDenSide3->Write();
-     fDenSide4->Write();
-     fDenSide5->Write();
-     fDenSide6->Write();
-     fDenSide7->Write();
-     fDenSide8->Write();
-   
-     fRatSide1->Write();
-     fRatSide2->Write();
-     fRatSide3->Write();
-     fRatSide4->Write();
-     fRatSide5->Write();
-     fRatSide6->Write();
-     fRatSide7->Write();
-     fRatSide8->Write();
-   
-     fNumLong1->Write();
-     fNumLong2->Write();
-     fNumLong3->Write();
-     fNumLong4->Write();
-     fNumLong5->Write();
-     fNumLong6->Write();
-     fNumLong7->Write();
-     fNumLong8->Write();
-
-     fDenLong1->Write();
-     fDenLong2->Write();
-     fDenLong3->Write();
-     fDenLong4->Write();
-     fDenLong5->Write();
-     fDenLong6->Write();
-     fDenLong7->Write();
-     fDenLong8->Write();
-   
-     fRatLong1->Write();
-     fRatLong2->Write();
-     fRatLong3->Write();
-     fRatLong4->Write();
-     fRatLong5->Write();
-     fRatLong6->Write();
-     fRatLong7->Write();
-     fRatLong8->Write(); 
-
-}
-
-//-------------------------------------
-void AliHBTashbtCorrFctn::ProcessSameEventParticles(AliHBTPair* pair)
-{
-    //Fills the numerator using pair from the same event
-    pair = CheckPair(pair);
-    if(pair == 0x0) return;
-   Double_t rplane=0.;   //reaction plane angle - 2 B determined
-   Double_t phi=(pair->Particle1()->Phi()+pair->Particle2()->Phi())/2.-rplane; //deltaphi bo nie mam nic innego pod reka
-   phi=phi*360/(2*TMath::Pi());
-   Double_t qout=pair->GetQOutLCMS();
-   Double_t qside=pair->GetQSideLCMS();
-   Double_t qlong=pair->GetQLongLCMS();
-   
-   if((phi>=0.) && (phi <45.))
-     {
-       fNumOut1->Fill(qout);
-       fNumSide1->Fill(qside);
-       fNumLong1->Fill(qlong);
-     }
-   else if((phi>=45.) && (phi <90.))
-     {
-       fNumOut2->Fill(qout);
-       fNumSide2->Fill(qside);
-       fNumLong2->Fill(qlong);
-     }
-   else if((phi>=90.) && (phi <135.))
-     {
-       fNumOut3->Fill(qout);
-       fNumSide3->Fill(qside);
-       fNumLong3->Fill(qlong);
-     }
-   else if((phi>=135.) && (phi <180.))
-     {
-       fNumOut4->Fill(qout);
-       fNumSide4->Fill(qside);
-       fNumLong4->Fill(qlong);
-     }
-   else if((phi>=180.) && (phi <225.))
-     {
-       fNumOut5->Fill(qout);
-       fNumSide5->Fill(qside);
-       fNumLong5->Fill(qlong);
-     }
-   else if((phi>=225.) && (phi <270.))
-     {
-       fNumOut6->Fill(qout);
-       fNumSide6->Fill(qside);
-       fNumLong6->Fill(qlong);
-     }
-   else if((phi>=270.) && (phi <315.))
-     {
-       fNumOut7->Fill(qout);
-       fNumSide7->Fill(qside);
-       fNumLong7->Fill(qlong);
-     }
-   else if((phi>=315.) && (phi <360.))
-     {
-       fNumOut8->Fill(qout);
-       fNumSide8->Fill(qside);
-       fNumLong8->Fill(qlong);
-     }
-
-   
-}
-
-/****************************************************************/
-void AliHBTashbtCorrFctn::Init()
-{
-     fNumOut1->Reset();
-     fNumOut2->Reset();
-     fNumOut3->Reset();
-     fNumOut4->Reset();
-     fNumOut5->Reset();
-     fNumOut6->Reset();
-     fNumOut7->Reset();
-     fNumOut8->Reset();
-
-     fDenOut1->Reset();
-     fDenOut2->Reset();
-     fDenOut3->Reset();
-     fDenOut4->Reset();
-     fDenOut5->Reset();
-     fDenOut6->Reset();
-     fDenOut7->Reset();
-     fDenOut8->Reset();
-   
-     fRatOut1->Reset();
-     fRatOut2->Reset();
-     fRatOut3->Reset();
-     fRatOut4->Reset();
-     fRatOut5->Reset();
-     fRatOut6->Reset();
-     fRatOut7->Reset();
-     fRatOut8->Reset();
-   
-
-     fNumSide1->Reset();
-     fNumSide2->Reset();
-     fNumSide3->Reset();
-     fNumSide4->Reset();
-     fNumSide5->Reset();
-     fNumSide6->Reset();
-     fNumSide7->Reset();
-     fNumSide8->Reset();
-
-     fDenSide1->Reset();
-     fDenSide2->Reset();
-     fDenSide3->Reset();
-     fDenSide4->Reset();
-     fDenSide5->Reset();
-     fDenSide6->Reset();
-     fDenSide7->Reset();
-     fDenSide8->Reset();
-   
-     fRatSide1->Reset();
-     fRatSide2->Reset();
-     fRatSide3->Reset();
-     fRatSide4->Reset();
-     fRatSide5->Reset();
-     fRatSide6->Reset();
-     fRatSide7->Reset();
-     fRatSide8->Reset();
-   
-     fNumLong1->Reset();
-     fNumLong2->Reset();
-     fNumLong3->Reset();
-     fNumLong4->Reset();
-     fNumLong5->Reset();
-     fNumLong6->Reset();
-     fNumLong7->Reset();
-     fNumLong8->Reset();
-
-     fDenLong1->Reset();
-     fDenLong2->Reset();
-     fDenLong3->Reset();
-     fDenLong4->Reset();
-     fDenLong5->Reset();
-     fDenLong6->Reset();
-     fDenLong7->Reset();
-     fDenLong8->Reset();
-   
-     fRatLong1->Reset();
-     fRatLong2->Reset();
-     fRatLong3->Reset();
-     fRatLong4->Reset();
-     fRatLong5->Reset();
-     fRatLong6->Reset();
-     fRatLong7->Reset();
-     fRatLong8->Reset(); 
-   
-
- }
-/****************************************************************/
-
-
-void AliHBTashbtCorrFctn::ProcessDiffEventParticles(AliHBTPair* pair)
-{
-
-      Double_t rplane=0.;   //reaction plane angle - 2 B determined
-   Double_t phi=(pair->Particle1()->Phi()+pair->Particle2()->Phi())/2.-rplane; //deltaphi bo nie mam nic innego pod reka
-   phi=phi*360/(2*TMath::Pi());
-   Double_t qout=pair->GetQOutLCMS();
-   Double_t qside=pair->GetQSideLCMS();
-   Double_t qlong=pair->GetQLongLCMS();
-   
-   if((phi>=0.) && (phi <45.))
-     {
-       fDenOut1->Fill(qout);
-       fDenSide1->Fill(qside);
-       fDenLong1->Fill(qlong);
-     }
-   else if((phi>=45.) && (phi <90.))
-     {
-       fDenOut2->Fill(qout);
-       fDenSide2->Fill(qside);
-       fDenLong2->Fill(qlong);
-     }
-   else if((phi>=90.) && (phi <135.))
-     {
-       fDenOut3->Fill(qout);
-       fDenSide3->Fill(qside);
-       fDenLong3->Fill(qlong);
-     }
-   else if((phi>=135.) && (phi <180.))
-     {
-       fDenOut4->Fill(qout);
-       fDenSide4->Fill(qside);
-       fDenLong4->Fill(qlong);
-     }
-   else if((phi>=180.) && (phi <225.))
-     {
-       fDenOut5->Fill(qout);
-       fDenSide5->Fill(qside);
-       fDenLong5->Fill(qlong);
-     }
-   else if((phi>=225.) && (phi <270.))
-     {
-       fDenOut6->Fill(qout);
-       fDenSide6->Fill(qside);
-       fDenLong6->Fill(qlong);
-     }
-   else if((phi>=270.) && (phi <315.))
-     {
-       fDenOut7->Fill(qout);
-       fDenSide7->Fill(qside);
-       fDenLong7->Fill(qlong);
-     }
-   else if((phi>=315.) && (phi <360.))
-     {
-       fDenOut8->Fill(qout);
-       fDenSide8->Fill(qside);
-       fDenLong8->Fill(qlong);
-     }
-
-}
-
-
-/******************************************************************/
-
-void AliHBTashbtCorrFctn::BuildHistos(Int_t nbins, Float_t max, Float_t min)
-{
-    
-   
-   TString nameNumOut1 = "NumOut0-45deg";
-   TString nameNumOut2 = "NumOut45-90deg";
-   TString nameNumOut3 = "NumOut90-135deg";
-   TString nameNumOut4 = "NumOut135-180deg";
-   TString nameNumOut5 = "NumOut180-225deg";
-   TString nameNumOut6 = "NumOut225-270deg";
-   TString nameNumOut7 = "NumOut270-315deg";
-   TString nameNumOut8 = "NumOut315-360deg";
-   
-   
-   TString nameDenOut1 = "DenOut0-45deg";
-   TString nameDenOut2 = "DenOut45-90deg";
-   TString nameDenOut3 = "DenOut90-135deg";
-   TString nameDenOut4 = "DenOut135-180deg";
-   TString nameDenOut5 = "DenOut180-225deg";
-   TString nameDenOut6 = "DenOut225-270deg";
-   TString nameDenOut7 = "DenOut270-315deg";
-   TString nameDenOut8 = "DenOut315-360deg";
-   
-   
-   TString nameRatOut1 = "RatOut0-45deg";
-   TString nameRatOut2 = "RatOut45-90deg";
-   TString nameRatOut3 = "RatOut90-135deg";
-   TString nameRatOut4 = "RatOut135-180deg";
-   TString nameRatOut5 = "RatOut180-225deg";
-   TString nameRatOut6 = "RatOut225-270deg";
-   TString nameRatOut7 = "RatOut270-315deg";
-   TString nameRatOut8 = "RatOut315-360deg";
-   
-   TString nameNumSide1 = "NumSide0-45deg";
-   TString nameNumSide2 = "NumSide45-90deg";
-   TString nameNumSide3 = "NumSide90-135deg";
-   TString nameNumSide4 = "NumSide135-180deg";
-   TString nameNumSide5 = "NumSide180-225deg";
-   TString nameNumSide6 = "NumSide225-270deg";
-   TString nameNumSide7 = "NumSide270-315deg";
-   TString nameNumSide8 = "NumSide315-360deg";
-   
-   
-   TString nameDenSide1 = "DenSide0-45deg";
-   TString nameDenSide2 = "DenSide45-90deg";
-   TString nameDenSide3 = "DenSide90-135deg";
-   TString nameDenSide4 = "DenSide135-180deg";
-   TString nameDenSide5 = "DenSide180-225deg";
-   TString nameDenSide6 = "DenSide225-270deg";
-   TString nameDenSide7 = "DenSide270-315deg";
-   TString nameDenSide8 = "DenSide315-360deg";
-   
-   
-   TString nameRatSide1 = "RatSide0-45deg";
-   TString nameRatSide2 = "RatSide45-90deg";
-   TString nameRatSide3 = "RatSide90-135deg";
-   TString nameRatSide4 = "RatSide135-180deg";
-   TString nameRatSide5 = "RatSide180-225deg";
-   TString nameRatSide6 = "RatSide225-270deg";
-   TString nameRatSide7 = "RatSide270-315deg";
-   TString nameRatSide8 = "RatSide315-360deg";
-   
-   TString nameNumLong1 = "NumLong0-45deg";
-   TString nameNumLong2 = "NumLong45-90deg";
-   TString nameNumLong3 = "NumLong90-135deg";
-   TString nameNumLong4 = "NumLong135-180deg";
-   TString nameNumLong5 = "NumLong180-225deg";
-   TString nameNumLong6 = "NumLong225-270deg";
-   TString nameNumLong7 = "NumLong270-315deg";
-   TString nameNumLong8 = "NumLong315-360deg";
-   
-   
-   TString nameDenLong1 = "DenLong0-45deg";
-   TString nameDenLong2 = "DenLong45-90deg";
-   TString nameDenLong3 = "DenLong90-135deg";
-   TString nameDenLong4 = "DenLong135-180deg";
-   TString nameDenLong5 = "DenLong180-225deg";
-   TString nameDenLong6 = "DenLong225-270deg";
-   TString nameDenLong7 = "DenLong270-315deg";
-   TString nameDenLong8 = "DenLong315-360deg";
-   
-   
-   TString nameRatLong1 = "RatLong0-45deg";
-   TString nameRatLong2 = "RatLong45-90deg";
-   TString nameRatLong3 = "RatLong90-135deg";
-   TString nameRatLong4 = "RatLong135-180deg";
-   TString nameRatLong5 = "RatLong180-225deg";
-   TString nameRatLong6 = "RatLong225-270deg";
-   TString nameRatLong7 = "RatLong270-315deg";
-   TString nameRatLong8 = "RatLong315-360deg";
-
-   fNumOut1 = new TH1D(nameNumOut1.Data(),nameNumOut1.Data(),nbins,min,max);
-   fNumOut2 = new TH1D(nameNumOut2.Data(),nameNumOut2.Data(),nbins,min,max);
-   fNumOut3 = new TH1D(nameNumOut3.Data(),nameNumOut3.Data(),nbins,min,max);
-   fNumOut4 = new TH1D(nameNumOut4.Data(),nameNumOut4.Data(),nbins,min,max);
-   fNumOut5 = new TH1D(nameNumOut5.Data(),nameNumOut5.Data(),nbins,min,max);
-   fNumOut6 = new TH1D(nameNumOut6.Data(),nameNumOut6.Data(),nbins,min,max);
-   fNumOut7 = new TH1D(nameNumOut7.Data(),nameNumOut7.Data(),nbins,min,max);
-   fNumOut8 = new TH1D(nameNumOut8.Data(),nameNumOut8.Data(),nbins,min,max);
-   
-   fDenOut1 = new TH1D(nameDenOut1.Data(),nameDenOut1.Data(),nbins,min,max);
-   fDenOut2 = new TH1D(nameDenOut2.Data(),nameDenOut2.Data(),nbins,min,max);
-   fDenOut3 = new TH1D(nameDenOut3.Data(),nameDenOut3.Data(),nbins,min,max);
-   fDenOut4 = new TH1D(nameDenOut4.Data(),nameDenOut4.Data(),nbins,min,max);
-   fDenOut5 = new TH1D(nameDenOut5.Data(),nameDenOut5.Data(),nbins,min,max);
-   fDenOut6 = new TH1D(nameDenOut6.Data(),nameDenOut6.Data(),nbins,min,max);
-   fDenOut7 = new TH1D(nameDenOut7.Data(),nameDenOut7.Data(),nbins,min,max);
-   fDenOut8 = new TH1D(nameDenOut8.Data(),nameDenOut8.Data(),nbins,min,max);
-   
-   fRatOut1 = new TH1D(nameRatOut1.Data(),nameRatOut1.Data(),nbins,min,max);
-   fRatOut2 = new TH1D(nameRatOut2.Data(),nameRatOut2.Data(),nbins,min,max);
-   fRatOut3 = new TH1D(nameRatOut3.Data(),nameRatOut3.Data(),nbins,min,max);
-   fRatOut4 = new TH1D(nameRatOut4.Data(),nameRatOut4.Data(),nbins,min,max);
-   fRatOut5 = new TH1D(nameRatOut5.Data(),nameRatOut5.Data(),nbins,min,max);
-   fRatOut6 = new TH1D(nameRatOut6.Data(),nameRatOut6.Data(),nbins,min,max);
-   fRatOut7 = new TH1D(nameRatOut7.Data(),nameRatOut7.Data(),nbins,min,max);
-   fRatOut8 = new TH1D(nameRatOut8.Data(),nameRatOut8.Data(),nbins,min,max);
-   
-   fNumSide1 = new TH1D(nameNumSide1.Data(),nameNumSide1.Data(),nbins,min,max);
-   fNumSide2 = new TH1D(nameNumSide2.Data(),nameNumSide2.Data(),nbins,min,max);
-   fNumSide3 = new TH1D(nameNumSide3.Data(),nameNumSide3.Data(),nbins,min,max);
-   fNumSide4 = new TH1D(nameNumSide4.Data(),nameNumSide4.Data(),nbins,min,max);
-   fNumSide5 = new TH1D(nameNumSide5.Data(),nameNumSide5.Data(),nbins,min,max);
-   fNumSide6 = new TH1D(nameNumSide6.Data(),nameNumSide6.Data(),nbins,min,max);
-   fNumSide7 = new TH1D(nameNumSide7.Data(),nameNumSide7.Data(),nbins,min,max);
-   fNumSide8 = new TH1D(nameNumSide8.Data(),nameNumSide8.Data(),nbins,min,max);
-   
-   fDenSide1 = new TH1D(nameDenSide1.Data(),nameDenSide1.Data(),nbins,min,max);
-   fDenSide2 = new TH1D(nameDenSide2.Data(),nameDenSide2.Data(),nbins,min,max);
-   fDenSide3 = new TH1D(nameDenSide3.Data(),nameDenSide3.Data(),nbins,min,max);
-   fDenSide4 = new TH1D(nameDenSide4.Data(),nameDenSide4.Data(),nbins,min,max);
-   fDenSide5 = new TH1D(nameDenSide5.Data(),nameDenSide5.Data(),nbins,min,max);
-   fDenSide6 = new TH1D(nameDenSide6.Data(),nameDenSide6.Data(),nbins,min,max);
-   fDenSide7 = new TH1D(nameDenSide7.Data(),nameDenSide7.Data(),nbins,min,max);
-   fDenSide8 = new TH1D(nameDenSide8.Data(),nameDenSide8.Data(),nbins,min,max);
-   
-   fRatSide1 = new TH1D(nameRatSide1.Data(),nameRatSide1.Data(),nbins,min,max);
-   fRatSide2 = new TH1D(nameRatSide2.Data(),nameRatSide2.Data(),nbins,min,max);
-   fRatSide3 = new TH1D(nameRatSide3.Data(),nameRatSide3.Data(),nbins,min,max);
-   fRatSide4 = new TH1D(nameRatSide4.Data(),nameRatSide4.Data(),nbins,min,max);
-   fRatSide5 = new TH1D(nameRatSide5.Data(),nameRatSide5.Data(),nbins,min,max);
-   fRatSide6 = new TH1D(nameRatSide6.Data(),nameRatSide6.Data(),nbins,min,max);
-   fRatSide7 = new TH1D(nameRatSide7.Data(),nameRatSide7.Data(),nbins,min,max);
-   fRatSide8 = new TH1D(nameRatSide8.Data(),nameRatSide8.Data(),nbins,min,max);
-   
-   fNumLong1 = new TH1D(nameNumLong1.Data(),nameNumLong1.Data(),nbins,min,max);
-   fNumLong2 = new TH1D(nameNumLong2.Data(),nameNumLong2.Data(),nbins,min,max);
-   fNumLong3 = new TH1D(nameNumLong3.Data(),nameNumLong3.Data(),nbins,min,max);
-   fNumLong4 = new TH1D(nameNumLong4.Data(),nameNumLong4.Data(),nbins,min,max);
-   fNumLong5 = new TH1D(nameNumLong5.Data(),nameNumLong5.Data(),nbins,min,max);
-   fNumLong6 = new TH1D(nameNumLong6.Data(),nameNumLong6.Data(),nbins,min,max);
-   fNumLong7 = new TH1D(nameNumLong7.Data(),nameNumLong7.Data(),nbins,min,max);
-   fNumLong8 = new TH1D(nameNumLong8.Data(),nameNumLong8.Data(),nbins,min,max);
-   
-   fDenLong1 = new TH1D(nameDenLong1.Data(),nameDenLong1.Data(),nbins,min,max);
-   fDenLong2 = new TH1D(nameDenLong2.Data(),nameDenLong2.Data(),nbins,min,max);
-   fDenLong3 = new TH1D(nameDenLong3.Data(),nameDenLong3.Data(),nbins,min,max);
-   fDenLong4 = new TH1D(nameDenLong4.Data(),nameDenLong4.Data(),nbins,min,max);
-   fDenLong5 = new TH1D(nameDenLong5.Data(),nameDenLong5.Data(),nbins,min,max);
-   fDenLong6 = new TH1D(nameDenLong6.Data(),nameDenLong6.Data(),nbins,min,max);
-   fDenLong7 = new TH1D(nameDenLong7.Data(),nameDenLong7.Data(),nbins,min,max);
-   fDenLong8 = new TH1D(nameDenLong8.Data(),nameDenLong8.Data(),nbins,min,max);
-   
-   fRatLong1 = new TH1D(nameRatLong1.Data(),nameRatLong1.Data(),nbins,min,max);
-   fRatLong2 = new TH1D(nameRatLong2.Data(),nameRatLong2.Data(),nbins,min,max);
-   fRatLong3 = new TH1D(nameRatLong3.Data(),nameRatLong3.Data(),nbins,min,max);
-   fRatLong4 = new TH1D(nameRatLong4.Data(),nameRatLong4.Data(),nbins,min,max);
-   fRatLong5 = new TH1D(nameRatLong5.Data(),nameRatLong5.Data(),nbins,min,max);
-   fRatLong6 = new TH1D(nameRatLong6.Data(),nameRatLong6.Data(),nbins,min,max);
-   fRatLong7 = new TH1D(nameRatLong7.Data(),nameRatLong7.Data(),nbins,min,max);
-   fRatLong8 = new TH1D(nameRatLong8.Data(),nameRatLong8.Data(),nbins,min,max);
-   
-   fNumOut1->Sumw2();
-   fNumOut2->Sumw2();
-     fNumOut3->Sumw2();
-     fNumOut4->Sumw2();
-     fNumOut5->Sumw2();
-     fNumOut6->Sumw2();
-     fNumOut7->Sumw2();
-     fNumOut8->Sumw2();
-
-     fDenOut1->Sumw2();
-     fDenOut2->Sumw2();
-     fDenOut3->Sumw2();
-     fDenOut4->Sumw2();
-     fDenOut5->Sumw2();
-     fDenOut6->Sumw2();
-     fDenOut7->Sumw2();
-     fDenOut8->Sumw2();
-   
-     fRatOut1->Sumw2();
-     fRatOut2->Sumw2();
-     fRatOut3->Sumw2();
-     fRatOut4->Sumw2();
-     fRatOut5->Sumw2();
-     fRatOut6->Sumw2();
-     fRatOut7->Sumw2();
-     fRatOut8->Sumw2();
-   
-
-     fNumSide1->Sumw2();
-     fNumSide2->Sumw2();
-     fNumSide3->Sumw2();
-     fNumSide4->Sumw2();
-     fNumSide5->Sumw2();
-     fNumSide6->Sumw2();
-     fNumSide7->Sumw2();
-     fNumSide8->Sumw2();
-
-     fDenSide1->Sumw2();
-     fDenSide2->Sumw2();
-     fDenSide3->Sumw2();
-     fDenSide4->Sumw2();
-     fDenSide5->Sumw2();
-     fDenSide6->Sumw2();
-     fDenSide7->Sumw2();
-     fDenSide8->Sumw2();
-   
-     fRatSide1->Sumw2();
-     fRatSide2->Sumw2();
-     fRatSide3->Sumw2();
-     fRatSide4->Sumw2();
-     fRatSide5->Sumw2();
-     fRatSide6->Sumw2();
-     fRatSide7->Sumw2();
-     fRatSide8->Sumw2();
-   
-     fNumLong1->Sumw2();
-     fNumLong2->Sumw2();
-     fNumLong3->Sumw2();
-     fNumLong4->Sumw2();
-     fNumLong5->Sumw2();
-     fNumLong6->Sumw2();
-     fNumLong7->Sumw2();
-     fNumLong8->Sumw2();
-
-     fDenLong1->Sumw2();
-     fDenLong2->Sumw2();
-     fDenLong3->Sumw2();
-     fDenLong4->Sumw2();
-     fDenLong5->Sumw2();
-     fDenLong6->Sumw2();
-     fDenLong7->Sumw2();
-     fDenLong8->Sumw2();
-   
-     fRatLong1->Sumw2();
-     fRatLong2->Sumw2();
-     fRatLong3->Sumw2();
-     fRatLong4->Sumw2();
-     fRatLong5->Sumw2();
-     fRatLong6->Sumw2();
-     fRatLong7->Sumw2();
-     fRatLong8->Sumw2();
-
- }
-
- TH1* AliHBTashbtCorrFctn::GetResult()
- {
-     return fRatOut1;
- }