]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWG3/base/AliQuarkoniaAcceptance.cxx
Transition PWG3/base -> PWGDQ/base
[u/mrichter/AliRoot.git] / PWG3 / base / AliQuarkoniaAcceptance.cxx
diff --git a/PWG3/base/AliQuarkoniaAcceptance.cxx b/PWG3/base/AliQuarkoniaAcceptance.cxx
deleted file mode 100644 (file)
index f73e058..0000000
+++ /dev/null
@@ -1,120 +0,0 @@
-/* $Id$ */
-
-// Root 
-#include "TAxis.h"
-#include "TFile.h"
-#include "TH2.h"
-#include "TString.h"
-
-// AliRoot includes
-#include "AliQuarkoniaAcceptance.h"
-#include "AliLog.h"
-
-
-
-//_______________________________________________________________________
-AliQuarkoniaAcceptance::AliQuarkoniaAcceptance(Int_t quarkoniaResonance, Int_t decayChannel):
-  fAcceptanceFileName("$ALICE_ROOT/PWG3/QuarkoniaAcceptance.root"),
-  fQuarkoniaResonance(quarkoniaResonance),
-  fDecayChannel(decayChannel),
-  fAcceptance(0x0)
-{
-
-
-}
-//_______________________________________________________________________
-AliQuarkoniaAcceptance::~AliQuarkoniaAcceptance()
-{
-  delete fAcceptance; 
-}
-
-//_______________________________________________________________________
-void AliQuarkoniaAcceptance::Init()
-{
-  switch (fQuarkoniaResonance)  {  
-  case kJpsi:
-    SetTitle("Jpsi");
-    break;
-  case kPsiP:
-    SetTitle("PsiP");
-    break;
-  case kUpsilon:
-    SetTitle("Upsilon");
-    break;
-  case kUpsilonP:
-    SetTitle("UpsilonP");
-    break;
-  case kUpsilonPP:
-    SetTitle("UpsilonPP");
-    break;
-  case kOmega:
-    SetTitle("Omega");
-    break;
-  case kPhi:
-    SetTitle("Phi");
-    break;
-  }
-
-  switch ( fDecayChannel) {
-  case kDimuon:
-    SetName("Dimuon");
-    break;
-  case kDielectron:
-    SetName("Dielectron");
-    break;
-  }
-
-  if(!fAcceptance) delete fAcceptance; 
-
-  TFile acceptanceFile(fAcceptanceFileName);
-  if ( acceptanceFile.IsOpen() ) {
-    char histoname[30];
-    snprintf(histoname,30,"h%s%sAccep",GetTitle(),GetName());
-    fAcceptance = (TH2F*) acceptanceFile.Get(histoname);
-    if ( !fAcceptance ) {
-      AliError(Form("Acceptance data for quarkonia %s and channel %s not found", GetTitle(), GetName() ));
-    }
-    else {
-      fAcceptance->SetDirectory(0);
-    }
-    acceptanceFile.Close();
-  }
-  else {
-    AliError(Form("File %s not found",fAcceptanceFileName.Data()));
-  }
-}
-//_______________________________________________________________________  
-TH2F*  AliQuarkoniaAcceptance::GetAcceptanceHisto() const
-{
-  if (fAcceptance) return fAcceptance;
-  else {
-    AliError(Form("Acceptance data for quarkonia %s and channel %s not found",GetTitle(),GetName()));
-    return 0x0;
-  }
-}
-//_______________________________________________________________________  
-void  AliQuarkoniaAcceptance::GetAcceptance(Float_t rap, Float_t pT, Double_t &accep, Double_t &error)
-{
-  Int_t binx=0;
-  Int_t biny=0;
-
-  if (!fAcceptance) {
-    AliError(Form("Acceptance data for quarkonia %s and channel %s not found",GetTitle(),GetName()));
-  }
-  else {
-    if (   rap < (fAcceptance->GetXaxis())->GetXmin()  ||  
-          rap > (fAcceptance->GetXaxis())->GetXmax()  ||
-          pT  < (fAcceptance->GetYaxis())->GetXmin()  ||  
-          pT  > (fAcceptance->GetYaxis())->GetXmax()     ) {
-      AliInfo("Values out of range");
-      accep = 0.;
-      error = 0.;
-    }
-    else  { 
-      binx  = fAcceptance->GetXaxis()->FindBin(rap);  
-      biny  = fAcceptance->GetYaxis()->FindBin(pT);
-      accep = fAcceptance->GetBinContent(binx,biny);
-      error = fAcceptance->GetBinError(binx,biny);
-    }
-  }
-}