//
// Object to store energy distribution corrections as used in the FMD
// analysis.
+//
+//
#include "AliFMDAnaCalibEnergyDistribution.h"
#include "TAxis.h"
#include <AliLog.h>
//____________________________________________________________________
void AliFMDAnaCalibEnergyDistribution::Init()
{
+ //Init object
if(fNetaBins == 0)
AliFatal("Set Eta bins before doing Init or anything else");
AliFMDAnaCalibEnergyDistribution::GetEnergyDistribution(Int_t det,
Char_t ring,
Float_t eta) {
-
+
+ //Get Energy dist
//TAxis testaxis(fNetaBins,fEtaMin,fEtaMax);
// Int_t binnumber = testaxis.FindBin(eta);
AliFMDAnaParameters* pars = AliFMDAnaParameters::Instance();
AliFMDAnaCalibEnergyDistribution::GetEmptyEnergyDistribution(Int_t det,
Char_t ring)
{
+ //Get e dist of empty
Int_t ringNumber = (ring == 'I' ? 0 : 1);
TObjArray* detArray = (TObjArray*)fEmptyArray.At(det);
TH1F*
AliFMDAnaCalibEnergyDistribution::GetRingEnergyDistribution(Int_t det,
Char_t ring) {
-
+ // Get E dist of ring
Int_t ringNumber = (ring == 'I' ? 0 : 1);
TObjArray* detArray = (TObjArray*)fRingArray.At(det);
Float_t eta,
TH1F* edist)
{
+ //Set E dist (user)
AliFMDAnaParameters* pars = AliFMDAnaParameters::Instance();
Int_t binnumber = pars->GetEtaBin(eta);
//std::cout<<binnumber<<std::endl;
Int_t etabin,
TH1F* edist)
{
+ //Set E dist
if(!fIsInit) Init();
Int_t ringNumber = (ring == 'I' ? 0 : 1);
Char_t ring,
TH1F* edist)
{
+ //Set the empty dist
if(!fIsInit)
Init();
Char_t ring,
TH1F* edist)
{
-
+ // Set E dist of ring
if(!fIsInit) Init();
Int_t ringNumber = (ring == 'I' ? 0 : 1);
//____________________________________________________________________
void AliFMDAnaCalibEnergyDistribution::Browse(TBrowser* b)
{
+ //Browse object
for(Int_t i = 0; i<fNetaBins; i++) {
TObjArray* etaArray = static_cast<TObjArray*>(fArray.At(i));
etaArray->SetName(Form("etabin_%03d", i+1));
#ifndef ALIFMDANACALIBENERGYDISTRIBUTION_H
#define ALIFMDANACALIBENERGYDISTRIBUTION_H
-
#include <TObject.h>
#include <TObjArray.h>
#include <TH1F.h>
class TBrowser;
-
/**
* @ingroup FMD_ana
* @brief Find most-probable value of MIP peak for one or more
AliFMDAnaCalibEnergyDistribution();
void SetNetaBins(Int_t nbins) {fNetaBins = nbins;}
- Int_t GetNetaBins() { return fNetaBins;}
+ Int_t GetNetaBins() const { return fNetaBins;}
void SetEtaLimits(Float_t eta_min, Float_t eta_max) {fEtaMin = eta_min; fEtaMax = eta_max;}
void SetEnergyDistributionUser(Int_t det, Char_t ring, Float_t eta, TH1F* edist);
void SetEnergyDistribution(Int_t det, Char_t ring, Int_t etabin, TH1F* edist);
void Browse(TBrowser* b);
protected:
void Init();
- TObjArray fArray;
- TObjArray fEmptyArray;
- TObjArray fRingArray;
- Bool_t fIsInit;
- Int_t fNetaBins;
- Float_t fEtaMax;
- Float_t fEtaMin;
+ TObjArray fArray; // Overall array
+ TObjArray fEmptyArray; // Array of empty events energy dep
+ TObjArray fRingArray; // Array of rings
+ Bool_t fIsInit; // Are we init?
+ Int_t fNetaBins; // Number of eta bins (intervals)
+ Float_t fEtaMax; // Max Eta
+ Float_t fEtaMin; // Min Eta
ClassDef(AliFMDAnaCalibEnergyDistribution,3);
};