Add model version of DPhi DEta CF

diff --git a/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoModelCorrFctnDEtaDPhi.cxx b/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoModelCorrFctnDEtaDPhi.cxx
new file mode 100644 (file)
index 0000000..d5aa9a1
--- /dev/null
+++ b/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoModelCorrFctnDEtaDPhi.cxx
@@ -0,0 +1,496 @@
+////////////////////////////////////////////////////////////////////////////////
+//                                                                            //
+// AliFemtoModelCorrFctnDEtaDPhi - A correlation function that analyzes            //
+// two particle correlations with respect to the azimuthal angle (phi)        //
+// and pseudorapidity (eta) difference                                        //
+//                                                                            //
+// Authors: Adam Kisiel Adam.Kisiel@cern.ch                                   //
+//                                                                            //
+////////////////////////////////////////////////////////////////////////////////
+
+#include "AliFemtoModelCorrFctnDEtaDPhi.h"
+//#include "AliFemtoHisto.hh"
+#include <cstdio>
+#include <TMath.h>
+
+#ifdef __ROOT__ 
+ClassImp(AliFemtoModelCorrFctnDEtaDPhi)
+#endif
+
+/*And some Model libraries..*/
+//1
+//#include "AliFemtoModelGausLCMSFreezeOutGenerator.h"
+//#include "AliFemtoModelHiddenInfo.h"
+//2
+#include "AliFemtoModelGausLCMSFreezeOutGenerator.h"
+#include "AliFemtoModelHiddenInfo.h"
+#include "AliFemtoPair.h"
+#include "AliFemtoModelManager.h"
+//
+
+//____________________________
+AliFemtoModelCorrFctnDEtaDPhi::AliFemtoModelCorrFctnDEtaDPhi(char* title, const int& aPhiBins=20, const int& aEtaBins=20):
+  AliFemtoModelCorrFctn(),
+  fDPhiDEtaNumeratorTrue(0),
+  fDPhiDEtaNumeratorFake(0),
+  fDPhiDEtaDenominator(0),
+  fDPhiDEtaColNumerator(0),
+  fDPhiDEtaColDenominator(0),
+  fDPhiNumeratorTrue(0),
+  fDPhiNumeratorFake(0),
+  fDPhiDenominator(0),
+  fDCosNumeratorTrue(0),
+  fDCosNumeratorFake(0),
+  fDCosDenominator(0),
+  fDPhiPtNumerator(0),
+  fDPhiPtDenominator(0),
+  fDCosPtNumerator(0),
+  fDCosPtDenominator(0)
+{
+  // set up numerator
+  char tTitNumDT[100] = "NumDPhiDEtaTrue";
+  strcat(tTitNumDT,title);
+  fDPhiDEtaNumeratorTrue = new TH2D(tTitNumDT,title,aPhiBins,-0.5*TMath::Pi(),1.5*TMath::Pi(),aEtaBins,-2.0,2.0);
+
+  char tTitNumDF[100] = "NumDPhiDEtaFake";
+  strcat(tTitNumDF,title);
+  fDPhiDEtaNumeratorFake = new TH2D(tTitNumDF,title,aPhiBins,-0.5*TMath::Pi(),1.5*TMath::Pi(),aEtaBins,-2.0,2.0);
+
+
+  // set up denominator
+  char tTitDenD[100] = "DenDPhiDEta";
+  strcat(tTitDenD,title);
+  fDPhiDEtaDenominator = new TH2D(tTitDenD,title,aPhiBins,-0.5*TMath::Pi(),1.5*TMath::Pi(),aEtaBins,-2.0,2.0);
+
+  // set up numerator
+  char tTitNumR[100] = "NumDPhiDEtaCol";
+  strcat(tTitNumR,title);
+  fDPhiDEtaColNumerator = new TH2D(tTitNumR,title,aPhiBins,-0.5*TMath::Pi(),1.5*TMath::Pi(),aEtaBins,-2.0,2.0);
+  // set up denominator
+  char tTitDenR[100] = "DenDPhiDEtaCol";
+  strcat(tTitDenR,title);
+  fDPhiDEtaColDenominator = new TH2D(tTitDenR,title,aPhiBins,-0.5*TMath::Pi(),1.5*TMath::Pi(),aEtaBins,-2.0,2.0);
+
+  // set up numerator
+  char tTitNumDPhiT[100] = "NumDPhiTrue";
+  strcat(tTitNumDPhiT,title);
+  fDPhiNumeratorTrue = new TH1D(tTitNumDPhiT,title,aPhiBins*2,-0.5*TMath::Pi(),1.5*TMath::Pi());
+  // set up numerator
+  char tTitNumDPhiF[100] = "NumDPhiFake";
+  strcat(tTitNumDPhiF,title);
+  fDPhiNumeratorFake = new TH1D(tTitNumDPhiF,title,aPhiBins*2,-0.5*TMath::Pi(),1.5*TMath::Pi());
+
+  // set up denominator
+  char tTitDenDPhi[100] = "DenDPhi";
+  strcat(tTitDenDPhi,title);
+  fDPhiDenominator = new TH1D(tTitDenDPhi,title,aPhiBins*2,-0.5*TMath::Pi(),1.5*TMath::Pi());
+
+  // set up numerator
+  char tTitNumDCosT[100] = "NumDCosTrue";
+  strcat(tTitNumDCosT,title);
+  fDCosNumeratorTrue = new TH1D(tTitNumDCosT,title,aPhiBins*2,-1.0,1.0);
+  // set up numerator
+  char tTitNumDCosF[100] = "NumDCosFake";
+  strcat(tTitNumDCosF,title);
+  fDCosNumeratorFake = new TH1D(tTitNumDCosF,title,aPhiBins*2,-1.0,1.0);
+
+  // set up denominator
+  char tTitDenDCos[100] = "DenDCos";
+  strcat(tTitDenDCos,title);
+  fDCosDenominator = new TH1D(tTitDenDCos,title,aPhiBins*2,-1.0,1.0);
+
+  // set up numerator
+  char tTitNumDPhiPt[100] = "NumDPhiPt";
+  strcat(tTitNumDPhiPt,title);
+  fDPhiPtNumerator = new TH2D(tTitNumDPhiPt,title,aPhiBins*2,-0.5*TMath::Pi(),1.5*TMath::Pi(), 30, 0.0, 3.0);
+  // set up denominator
+  char tTitDenDPhiPt[100] = "DenDPhiPt";
+  strcat(tTitDenDPhiPt,title);
+  fDPhiPtDenominator = new TH2D(tTitDenDPhiPt,title,aPhiBins*2,-0.5*TMath::Pi(),1.5*TMath::Pi(), 30, 0.0, 3.0);
+
+  // set up numerator
+  char tTitNumDCosPt[100] = "NumDCosPt";
+  strcat(tTitNumDCosPt,title);
+  fDCosPtNumerator = new TH2D(tTitNumDCosPt,title,aPhiBins*2,-1.0,1.0, 30, 0.0, 3.0);
+  // set up denominator
+  char tTitDenDCosPt[100] = "DenDCosPt";
+  strcat(tTitDenDCosPt,title);
+  fDCosPtDenominator = new TH2D(tTitDenDCosPt,title,aPhiBins*2,-1.0,1.0, 30, 0.0, 3.0);
+
+  // to enable error bar calculation...
+  fDPhiDEtaNumeratorTrue->Sumw2();
+  fDPhiDEtaNumeratorFake->Sumw2();
+  fDPhiDEtaDenominator->Sumw2();
+  fDPhiDEtaColNumerator->Sumw2();
+  fDPhiDEtaColDenominator->Sumw2();
+  fDPhiNumeratorTrue->Sumw2();
+  fDPhiNumeratorFake->Sumw2();
+  fDPhiDenominator->Sumw2();
+  fDCosNumeratorTrue->Sumw2();
+  fDCosNumeratorFake->Sumw2();
+  fDCosDenominator->Sumw2();
+  fDPhiPtNumerator->Sumw2();
+  fDPhiPtDenominator->Sumw2();
+  fDCosPtNumerator->Sumw2();
+  fDCosPtDenominator->Sumw2();
+
+}
+
+//____________________________
+AliFemtoModelCorrFctnDEtaDPhi::AliFemtoModelCorrFctnDEtaDPhi(const AliFemtoModelCorrFctnDEtaDPhi& aCorrFctn) :
+  AliFemtoModelCorrFctn(),
+  fDPhiDEtaNumeratorTrue(0),
+  fDPhiDEtaNumeratorFake(0),
+  fDPhiDEtaDenominator(0),
+  fDPhiDEtaColNumerator(0),
+  fDPhiDEtaColDenominator(0),
+  fDPhiNumeratorTrue(0),
+  fDPhiNumeratorFake(0),
+  fDPhiDenominator(0),
+  fDCosNumeratorTrue(0),
+  fDCosNumeratorFake(0),
+  fDCosDenominator(0),
+  fDPhiPtNumerator(0),
+  fDPhiPtDenominator(0),
+  fDCosPtNumerator(0),
+  fDCosPtDenominator(0)
+{
+  // copy constructor
+  if (aCorrFctn.fDPhiDEtaNumeratorTrue)
+    fDPhiDEtaNumeratorTrue = new TH2D(*aCorrFctn.fDPhiDEtaNumeratorTrue);
+  else
+    fDPhiDEtaNumeratorTrue = 0;
+  if (aCorrFctn.fDPhiDEtaNumeratorFake)
+    fDPhiDEtaNumeratorFake = new TH2D(*aCorrFctn.fDPhiDEtaNumeratorFake);
+  else
+    fDPhiDEtaNumeratorFake = 0;
+
+  if (aCorrFctn.fDPhiDEtaDenominator)
+    fDPhiDEtaDenominator = new TH2D(*aCorrFctn.fDPhiDEtaDenominator);
+  else
+    fDPhiDEtaDenominator = 0;
+
+  if (aCorrFctn.fDPhiDEtaColNumerator)
+    fDPhiDEtaColNumerator = new TH2D(*aCorrFctn.fDPhiDEtaColNumerator);
+  else
+    fDPhiDEtaColNumerator = 0;
+  if (aCorrFctn.fDPhiDEtaColDenominator)
+    fDPhiDEtaColDenominator = new TH2D(*aCorrFctn.fDPhiDEtaColDenominator);
+  else
+    fDPhiDEtaColDenominator = 0;
+
+  if (aCorrFctn.fDPhiNumeratorTrue)
+    fDPhiNumeratorTrue = new TH1D(*aCorrFctn.fDPhiNumeratorTrue);
+  else
+    fDPhiNumeratorTrue = 0;
+
+  if (aCorrFctn.fDPhiNumeratorFake)
+    fDPhiNumeratorFake = new TH1D(*aCorrFctn.fDPhiNumeratorFake);
+  else
+    fDPhiNumeratorFake = 0;
+
+  if (aCorrFctn.fDPhiDenominator)
+    fDPhiDenominator = new TH1D(*aCorrFctn.fDPhiDenominator);
+  else
+    fDPhiDenominator = 0;
+
+  if (aCorrFctn.fDCosNumeratorTrue)
+    fDCosNumeratorTrue = new TH1D(*aCorrFctn.fDCosNumeratorTrue);
+  else
+    fDCosNumeratorTrue = 0;
+  if (aCorrFctn.fDCosNumeratorFake)
+    fDCosNumeratorFake = new TH1D(*aCorrFctn.fDCosNumeratorFake);
+  else
+    fDCosNumeratorFake = 0;
+
+  if (aCorrFctn.fDCosDenominator)
+    fDCosDenominator = new TH1D(*aCorrFctn.fDCosDenominator);
+  else
+    fDCosDenominator = 0;
+
+  if (aCorrFctn.fDPhiPtNumerator)
+    fDPhiPtNumerator = new TH2D(*aCorrFctn.fDPhiPtNumerator);
+  else
+    fDPhiPtNumerator = 0;
+  if (aCorrFctn.fDPhiPtDenominator)
+    fDPhiPtDenominator = new TH2D(*aCorrFctn.fDPhiPtDenominator);
+  else
+    fDPhiPtDenominator = 0;
+
+  if (aCorrFctn.fDCosPtNumerator)
+    fDCosPtNumerator = new TH2D(*aCorrFctn.fDCosPtNumerator);
+  else
+    fDCosPtNumerator = 0;
+  if (aCorrFctn.fDCosPtDenominator)
+    fDCosPtDenominator = new TH2D(*aCorrFctn.fDCosPtDenominator);
+  else
+    fDCosPtDenominator = 0;
+}
+//____________________________
+AliFemtoModelCorrFctnDEtaDPhi::~AliFemtoModelCorrFctnDEtaDPhi(){
+  // destructor
+  if(fDPhiDEtaNumeratorTrue)      delete fDPhiDEtaNumeratorTrue;
+  if(fDPhiDEtaNumeratorFake)      delete fDPhiDEtaNumeratorFake;
+  if(fDPhiDEtaDenominator)    delete fDPhiDEtaDenominator;
+  if(fDPhiDEtaColNumerator)   delete fDPhiDEtaColNumerator;
+  if(fDPhiDEtaColDenominator) delete fDPhiDEtaColDenominator;
+  if(fDPhiNumeratorTrue)       delete fDPhiNumeratorTrue;
+  if(fDPhiNumeratorFake)       delete fDPhiNumeratorFake;
+  if(fDPhiDenominator)     delete fDPhiDenominator;
+  if(fDCosNumeratorTrue)       delete fDCosNumeratorTrue;
+  if(fDCosNumeratorFake)       delete fDCosNumeratorFake;
+  if(fDCosDenominator)     delete fDCosDenominator;
+  if(fDPhiPtNumerator)     delete fDPhiPtNumerator;
+  if(fDPhiPtDenominator)   delete fDPhiPtDenominator;
+  if(fDCosPtNumerator)     delete fDCosPtNumerator;
+  if(fDCosPtDenominator)   delete fDCosPtDenominator;
+}
+//_________________________
+AliFemtoModelCorrFctnDEtaDPhi& AliFemtoModelCorrFctnDEtaDPhi::operator=(const AliFemtoModelCorrFctnDEtaDPhi& aCorrFctn)
+{
+  // assignment operator
+  if (this == &aCorrFctn)
+    return *this;
+
+  if (aCorrFctn.fDPhiDEtaNumeratorTrue)
+    fDPhiDEtaNumeratorTrue = new TH2D(*aCorrFctn.fDPhiDEtaNumeratorTrue);
+  else
+    fDPhiDEtaNumeratorTrue = 0;
+
+  if (aCorrFctn.fDPhiDEtaNumeratorFake)
+    fDPhiDEtaNumeratorFake = new TH2D(*aCorrFctn.fDPhiDEtaNumeratorFake);
+  else
+    fDPhiDEtaNumeratorFake = 0;
+
+  if (aCorrFctn.fDPhiDEtaDenominator)
+    fDPhiDEtaDenominator = new TH2D(*aCorrFctn.fDPhiDEtaDenominator);
+  else
+    fDPhiDEtaDenominator = 0;
+
+  if (aCorrFctn.fDPhiDEtaColNumerator)
+    fDPhiDEtaColNumerator = new TH2D(*aCorrFctn.fDPhiDEtaColNumerator);
+  else
+    fDPhiDEtaColNumerator = 0;
+  if (aCorrFctn.fDPhiDEtaColDenominator)
+    fDPhiDEtaColDenominator = new TH2D(*aCorrFctn.fDPhiDEtaColDenominator);
+  else
+    fDPhiDEtaColDenominator = 0;
+
+  if (aCorrFctn.fDPhiNumeratorTrue)
+    fDPhiNumeratorTrue = new TH1D(*aCorrFctn.fDPhiNumeratorTrue);
+  else
+    fDPhiNumeratorTrue = 0;
+  if (aCorrFctn.fDPhiNumeratorFake)
+    fDPhiNumeratorFake = new TH1D(*aCorrFctn.fDPhiNumeratorFake);
+  else
+    fDPhiNumeratorFake = 0;
+
+  if (aCorrFctn.fDPhiDenominator)
+    fDPhiDenominator = new TH1D(*aCorrFctn.fDPhiDenominator);
+  else
+    fDPhiDenominator = 0;
+
+  if (aCorrFctn.fDCosNumeratorTrue)
+    fDCosNumeratorTrue = new TH1D(*aCorrFctn.fDCosNumeratorTrue);
+  else
+    fDCosNumeratorTrue = 0;
+  if (aCorrFctn.fDCosNumeratorFake)
+    fDCosNumeratorFake = new TH1D(*aCorrFctn.fDCosNumeratorFake);
+  else
+    fDCosNumeratorFake = 0;
+
+  if (aCorrFctn.fDCosDenominator)
+    fDCosDenominator = new TH1D(*aCorrFctn.fDCosDenominator);
+  else
+    fDCosDenominator = 0;
+
+  if (aCorrFctn.fDPhiPtNumerator)
+    fDPhiPtNumerator = new TH2D(*aCorrFctn.fDPhiPtNumerator);
+  else
+    fDPhiPtNumerator = 0;
+  if (aCorrFctn.fDPhiPtDenominator)
+    fDPhiPtDenominator = new TH2D(*aCorrFctn.fDPhiPtDenominator);
+  else
+    fDPhiPtDenominator = 0;
+
+  if (aCorrFctn.fDCosPtNumerator)
+    fDCosPtNumerator = new TH2D(*aCorrFctn.fDCosPtNumerator);
+  else
+    fDCosPtNumerator = 0;
+  if (aCorrFctn.fDCosPtDenominator)
+    fDCosPtDenominator = new TH2D(*aCorrFctn.fDCosPtDenominator);
+  else
+    fDCosPtDenominator = 0;
+
+  return *this;
+}
+//_________________________
+void AliFemtoModelCorrFctnDEtaDPhi::Finish(){
+  // here is where we should normalize, fit, etc...
+  // we should NOT Draw() the histos (as I had done it below),
+  // since we want to insulate ourselves from root at this level
+  // of the code.  Do it instead at root command line with browser.
+  //  mShareNumerator->Draw();
+  //mShareDenominator->Draw();
+  //mRatio->Draw();
+
+}
+
+//____________________________
+AliFemtoString AliFemtoModelCorrFctnDEtaDPhi::Report(){
+  // create report
+  string stemp = "TPC Ncls Correlation Function Report:\n";
+  char ctemp[100];
+  sprintf(ctemp,"Number of entries in numerator true:\t%E\n",fDPhiDEtaNumeratorTrue->GetEntries());
+  sprintf(ctemp,"Number of entries in numerator fake:\t%E\n",fDPhiDEtaNumeratorFake->GetEntries());
+  stemp += ctemp;
+  sprintf(ctemp,"Number of entries in denominator:\t%E\n",fDPhiDEtaDenominator->GetEntries());
+  stemp += ctemp;
+  //  stemp += mCoulombWeight->Report();
+  AliFemtoString returnThis = stemp;
+  return returnThis;
+}
+//____________________________
+void AliFemtoModelCorrFctnDEtaDPhi::AddRealPair( AliFemtoPair* pair){
+  // add real (effect) pair
+  double phi1 = pair->Track1()->Track()->P().phi();
+  double phi2 = pair->Track2()->Track()->P().phi();
+  double eta1 = pair->Track1()->Track()->P().pseudoRapidity();
+  double eta2 = pair->Track2()->Track()->P().pseudoRapidity();
+
+  double dphi = phi1 - phi2;
+  while (dphi<-TMath::Pi()/2) dphi+=TMath::Pi()*2;
+  while (dphi>3*TMath::Pi()/2) dphi-=TMath::Pi()*2;
+
+  double deta = eta1 - eta2;
+
+  double px1 = pair->Track1()->Track()->P().x();
+  double py1 = pair->Track1()->Track()->P().y();
+  double pz1 = pair->Track1()->Track()->P().z();
+
+  double px2 = pair->Track2()->Track()->P().x();
+  double py2 = pair->Track2()->Track()->P().y();
+  double pz2 = pair->Track2()->Track()->P().z();
+
+  double pt1 = TMath::Hypot(px1, py1);
+  double pt2 = TMath::Hypot(px2, py2);
+  double ptmin = pt1>pt2 ? pt2 : pt1;
+
+  double cosphi = (px1*px2 + py1*py2 + pz1*pz2)/
+    sqrt((px1*px1 + py1*py1 + pz1*pz1)*(px2*px2 + py2*py2 + pz2*pz2));
+
+/*weights*/
+  Double_t weight = fManager->GetWeight(pair);
+  fDPhiDEtaNumeratorTrue->Fill(dphi, deta,weight);
+
+  if (cosphi > 0) {
+    fDPhiDEtaColNumerator->Fill(dphi, deta);
+  }
+  else {
+    fDPhiDEtaColNumerator->Fill(dphi, -eta1-eta2);
+  }
+
+  fDPhiNumeratorTrue->Fill(dphi,weight);
+  fDCosNumeratorTrue->Fill(cosphi,weight);
+
+  fDPhiPtNumerator->Fill(dphi, ptmin);
+  fDCosPtNumerator->Fill(cosphi, ptmin);
+
+}
+//____________________________
+void AliFemtoModelCorrFctnDEtaDPhi::AddMixedPair( AliFemtoPair* pair){
+  // add mixed (background) pair
+  double phi1 = pair->Track1()->Track()->P().phi();
+  double phi2 = pair->Track2()->Track()->P().phi();
+  double eta1 = pair->Track1()->Track()->P().pseudoRapidity();
+  double eta2 = pair->Track2()->Track()->P().pseudoRapidity();
+
+  double dphi = phi1 - phi2;
+  while (dphi<-TMath::Pi()/2) dphi+=TMath::Pi()*2;
+  while (dphi>3*TMath::Pi()/2) dphi-=TMath::Pi()*2;
+
+  double deta = eta1 - eta2;
+
+  double px1 = pair->Track1()->Track()->P().x();
+  double py1 = pair->Track1()->Track()->P().y();
+  double pz1 = pair->Track1()->Track()->P().z();
+
+  double px2 = pair->Track2()->Track()->P().x();
+  double py2 = pair->Track2()->Track()->P().y();
+  double pz2 = pair->Track2()->Track()->P().z();
+
+  double pt1 = TMath::Hypot(px1, py1);
+  double pt2 = TMath::Hypot(px2, py2);
+  double ptmin = pt1>pt2 ? pt2 : pt1;
+
+  double cosphi = (px1*px2 + py1*py2 + pz1*pz2)/
+    sqrt((px1*px1 + py1*py1 + pz1*pz1)*(px2*px2 + py2*py2 + pz2*pz2));
+
+
+  Double_t weight = fManager->GetWeight(pair);
+  fDPhiDEtaNumeratorFake->Fill(dphi, deta,weight);
+
+  fDPhiDEtaDenominator->Fill(dphi, deta,1.0);
+
+  if (cosphi > 0) {
+    fDPhiDEtaColDenominator->Fill(dphi, deta);
+  }
+  else {
+    fDPhiDEtaColDenominator->Fill(dphi, -eta1-eta2);
+  }
+
+  fDPhiNumeratorFake->Fill(dphi,weight);
+  fDCosNumeratorFake->Fill(cosphi,weight);
+
+  fDPhiDenominator->Fill(dphi,1.0);
+  fDCosDenominator->Fill(cosphi,1.0);
+
+  fDPhiPtDenominator->Fill(dphi, ptmin);
+  fDCosPtDenominator->Fill(cosphi, ptmin);
+}
+
+
+void AliFemtoModelCorrFctnDEtaDPhi::WriteHistos()
+{
+  // Write out result histograms
+  fDPhiDEtaNumeratorTrue->Write();
+  fDPhiDEtaNumeratorFake->Write();
+  fDPhiDEtaDenominator->Write();
+  fDPhiDEtaColNumerator->Write();
+  fDPhiDEtaColDenominator->Write();
+  fDPhiNumeratorTrue->Write();
+  fDPhiNumeratorFake->Write();
+  fDPhiDenominator->Write();
+  fDCosNumeratorTrue->Write();
+  fDCosNumeratorFake->Write();
+  fDCosDenominator->Write();
+  fDPhiPtNumerator->Write();
+  fDPhiPtDenominator->Write();
+  fDCosPtNumerator->Write();
+  fDCosPtDenominator->Write();
+}
+
+TList* AliFemtoModelCorrFctnDEtaDPhi::GetOutputList()
+{
+  // Prepare the list of objects to be written to the output
+  TList *tOutputList = new TList();
+
+  tOutputList->Add(fDPhiDEtaNumeratorTrue);
+  tOutputList->Add(fDPhiDEtaNumeratorFake);
+  tOutputList->Add(fDPhiDEtaDenominator);
+  tOutputList->Add(fDPhiDEtaColNumerator);
+  tOutputList->Add(fDPhiDEtaColDenominator);
+  tOutputList->Add(fDPhiNumeratorTrue);
+  tOutputList->Add(fDPhiNumeratorFake);
+  tOutputList->Add(fDPhiDenominator);
+  tOutputList->Add(fDCosNumeratorTrue);
+  tOutputList->Add(fDCosNumeratorFake);
+  tOutputList->Add(fDCosDenominator);
+  tOutputList->Add(fDPhiPtNumerator);
+  tOutputList->Add(fDPhiPtDenominator);
+  tOutputList->Add(fDCosPtNumerator);
+  tOutputList->Add(fDCosPtDenominator);
+
+  return tOutputList;
+
+}

diff --git a/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoModelCorrFctnDEtaDPhi.h b/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoModelCorrFctnDEtaDPhi.h
new file mode 100644 (file)
index 0000000..6e2a880
--- /dev/null
+++ b/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoModelCorrFctnDEtaDPhi.h
@@ -0,0 +1,67 @@
+////////////////////////////////////////////////////////////////////////////////
+//                                                                            //
+// AliFemtoCorrFctnDEtaDPhi - A correlation function that analyzes            //
+// two particle correlations with respect to the azimuthal angle (phi)        //
+// and pseudorapidity (eta) difference                                        //
+//                                                                            //
+// Authors: Adam Kisiel Adam.Kisiel@cern.ch                                   //
+//                                                                            //
+////////////////////////////////////////////////////////////////////////////////
+
+#ifndef ALIFEMTOMODELCORRFCTNDETADPHI_H
+#define ALIFEMTOMODELCORRFCTNDETADPHI_H
+
+#include "TH1D.h"
+#include "TH2D.h"
+#include "AliFemtoCorrFctn.h"
+#include "AliFemtoPair.h"
+#include "AliFemtoModelManager.h"
+#include "AliFemtoModelCorrFctn.h"
+
+class AliFemtoModelCorrFctnDEtaDPhi : public AliFemtoModelCorrFctn {
+public:
+  AliFemtoModelCorrFctnDEtaDPhi(char* title, const int& aPhiBins, const int& aEtaBins);
+  AliFemtoModelCorrFctnDEtaDPhi(const AliFemtoModelCorrFctnDEtaDPhi& aCorrFctn);
+  virtual ~AliFemtoModelCorrFctnDEtaDPhi();
+
+  AliFemtoModelCorrFctnDEtaDPhi& operator=(const AliFemtoModelCorrFctnDEtaDPhi& aCorrFctn);
+
+  virtual AliFemtoString Report();
+  virtual void AddRealPair(AliFemtoPair* aPair);
+  virtual void AddMixedPair(AliFemtoPair* aPair);
+
+  virtual void Finish();
+
+  void WriteHistos();
+  virtual TList* GetOutputList();
+private:
+  
+  TH2D *fDPhiDEtaNumeratorTrue;      // Numerator of dEta dPhi true function
+  TH2D *fDPhiDEtaNumeratorFake;      // Numerator of dEta dPhi fake function
+  TH2D *fDPhiDEtaDenominator;        // Denominator of dEta dPhi function
+
+  TH2D *fDPhiDEtaColNumerator;       // Numerator of colinear dEta dPhi function 
+  TH2D *fDPhiDEtaColDenominator;     // Denominator of colinear dEta dPhi function
+
+  TH1D *fDPhiNumeratorTrue;          // Numerator of dPhi true correlation
+  TH1D *fDPhiNumeratorFake;          // Numerator of dPhi fake correlation
+  TH1D *fDPhiDenominator;            // Denominator of dPhi correlation
+
+  TH1D *fDCosNumeratorTrue;           // Numerator of colinearity true correlation
+  TH1D *fDCosNumeratorFake;           // Numerator of colinearity fake correlation
+  TH1D *fDCosDenominator;            // Denominator of colinearity correlation
+
+  TH2D *fDPhiPtNumerator;            // Numerator of dPhi correlation vs. Pt min
+  TH2D *fDPhiPtDenominator;          // Denominator of dPhi correlation vs. Pt min
+
+  TH2D *fDCosPtNumerator;            // Numerator of colinearity correlation vs. Pt min
+  TH2D *fDCosPtDenominator;          // Denominator of colinearity correlation vs. Pt min
+
+#ifdef __ROOT__
+  ClassDef(AliFemtoModelCorrFctnDEtaDPhi, 1)
+#endif
+};
+
+
+#endif
+

diff --git a/PWG2/PWG2femtoscopyUserLinkDef.h b/PWG2/PWG2femtoscopyUserLinkDef.h
index e6b79ac93538662613a9b90d34019c4f8e17c28e..6157178205dcc4f54ac0645e317fa6129ff7ed31 100644 (file)
--- a/PWG2/PWG2femtoscopyUserLinkDef.h
+++ b/PWG2/PWG2femtoscopyUserLinkDef.h
@@ -11,6 +11,7 @@
 #pragma link C++ class AliFemtoTPCInnerCorrFctn;
 #pragma link C++ class AliFemtoChi2CorrFctn;
 #pragma link C++ class AliFemtoModelCorrFctnSource;
+#pragma link C++ class AliFemtoModelCorrFctnDEtaDPhi;
 #pragma link C++ class AliFemtoModelBPLCMSCorrFctn;
 #pragma link C++ class AliFemtoModelGausRinvFreezeOutGenerator;
 #pragma link C++ class AliFemtoModelCorrFctnTrueQ;
@@ -34,3 +35,4 @@
 #pragma link C++ class AliFemtoCorrFctnDEtaDPhi;
 #pragma link C++ class AliFemtoCutMonitorParticlePtPDG;
 #pragma link C++ class AliFemtoCutMonitorParticleEtCorr;
+#pragma link C++ class AliFemtoPhiPairCut;

diff --git a/PWG2/libPWG2femtoscopyUser.pkg b/PWG2/libPWG2femtoscopyUser.pkg
index 37aadae1223c2dc7b2b39b3eb5d11c65964e12cf..203f31014df20ffa1c50957de80360b7d74422bc 100644 (file)
--- a/PWG2/libPWG2femtoscopyUser.pkg
+++ b/PWG2/libPWG2femtoscopyUser.pkg
@@ -8,6 +8,7 @@ SRCS= FEMTOSCOPY/AliFemtoUser/AliFemtoQPairCut.cxx \
       FEMTOSCOPY/AliFemtoUser/AliFemtoTPCInnerCorrFctn.cxx \
       FEMTOSCOPY/AliFemtoUser/AliFemtoChi2CorrFctn.cxx \
       FEMTOSCOPY/AliFemtoUser/AliFemtoModelCorrFctnSource.cxx \
+      FEMTOSCOPY/AliFemtoUser/AliFemtoModelCorrFctnDEtaDPhi.cxx \
       FEMTOSCOPY/AliFemtoUser/AliFemtoModelBPLCMSCorrFctn.cxx \
       FEMTOSCOPY/AliFemtoUser/AliFemtoModelGausRinvFreezeOutGenerator.cxx \
       FEMTOSCOPY/AliFemtoUser/AliFemtoModelCorrFctnTrueQ.cxx \
@@ -24,6 +25,7 @@ SRCS= FEMTOSCOPY/AliFemtoUser/AliFemtoQPairCut.cxx \
       FEMTOSCOPY/AliFemtoUser/AliFemtoShareQualityQAPairCut.cxx \
       FEMTOSCOPY/AliFemtoUser/AliFemtoShareQualityCorrFctn.cxx \
       FEMTOSCOPY/AliFemtoUser/AliFemtoYlm.cxx \
+      FEMTOSCOPY/AliFemtoUser/AliFemtoPhiPairCut.cxx \
       FEMTOSCOPY/AliFemtoUser/AliFemtoCorrFctnDirectYlm.cxx \
       FEMTOSCOPY/AliFemtoUser/AliFemtoModelCorrFctnDirectYlm.cxx \
       FEMTOSCOPY/AliFemtoUser/AliFemtoCutMonitorParticlePtPDG.cxx \