10 #include "AliQuarkoniaAcceptance.h"
15 //_______________________________________________________________________
16 AliQuarkoniaAcceptance::AliQuarkoniaAcceptance(Int_t quarkoniaResonance, Int_t decayChannel):
17 fAcceptanceFileName("$ALICE_ROOT/PWG3/QuarkoniaAcceptance.root"),
18 fQuarkoniaResonance(quarkoniaResonance),
19 fDecayChannel(decayChannel),
25 //_______________________________________________________________________
26 AliQuarkoniaAcceptance::~AliQuarkoniaAcceptance()
31 //_______________________________________________________________________
32 void AliQuarkoniaAcceptance::Init()
34 switch (fQuarkoniaResonance) {
48 SetTitle("UpsilonPP");
58 switch ( fDecayChannel) {
63 SetName("Dielectron");
67 if(!fAcceptance) delete fAcceptance;
69 TFile acceptanceFile(fAcceptanceFileName);
70 if ( acceptanceFile.IsOpen() ) {
72 sprintf(histoname,"h%s%sAccep",GetTitle(),GetName());
73 fAcceptance = (TH2F*) acceptanceFile.Get(histoname);
75 AliError(Form("Acceptance data for quarkonia %s and channel %s not found", GetTitle(), GetName() ));
78 fAcceptance->SetDirectory(0);
80 acceptanceFile.Close();
83 AliError(Form("File %s not found",fAcceptanceFileName.Data()));
86 //_______________________________________________________________________
87 TH2F* AliQuarkoniaAcceptance::GetAcceptanceHisto() const
89 if (fAcceptance) return fAcceptance;
91 AliError(Form("Acceptance data for quarkonia %s and channel %s not found",GetTitle(),GetName()));
95 //_______________________________________________________________________
96 void AliQuarkoniaAcceptance::GetAcceptance(Float_t rap, Float_t pT, Double_t &accep, Double_t &error)
102 AliError(Form("Acceptance data for quarkonia %s and channel %s not found",GetTitle(),GetName()));
105 if ( rap < (fAcceptance->GetXaxis())->GetXmin() ||
106 rap > (fAcceptance->GetXaxis())->GetXmax() ||
107 pT < (fAcceptance->GetYaxis())->GetXmin() ||
108 pT > (fAcceptance->GetYaxis())->GetXmax() ) {
109 AliInfo("Values out of range");
114 binx = fAcceptance->GetXaxis()->FindBin(rap);
115 biny = fAcceptance->GetYaxis()->FindBin(pT);
116 accep = fAcceptance->GetBinContent(binx,biny);
117 error = fAcceptance->GetBinError(binx,biny);