AliAnalysisTaskEmcalClustersRef.h

Go to the documentation of this file.

#ifndef ALIANALYSISTASKEMCALCLUSTERSREF_H
#define ALIANALYSISTASKEMCALCLUSTERSREF_H
/* Copyright(c) 1998-2015, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

#include "AliAnalysisTaskEmcalTriggerBase.h"
#include "AliCutValueRange.h"
#include <TCustomBinning.h>
#include <TString.h>

class TClonesArray;

namespace EMCalTriggerPtAnalysis {

class AliAnalysisTaskEmcalClustersRef : public AliAnalysisTaskEmcalTriggerBase {
public:

  enum EnergyDefinition_t {
    kDefaultEnergy,       
    kNonLinCorrEnergy,    
    kHadCorrEnergy        
  };

  AliAnalysisTaskEmcalClustersRef();

  AliAnalysisTaskEmcalClustersRef(const char *name);

  virtual ~AliAnalysisTaskEmcalClustersRef();

  void EnableSumw2(Bool_t doEnable) { fEnableSumw2 = doEnable; }

  void SetEnergyDefinition(EnergyDefinition_t edef) { fEnergyDefinition = edef; }

  void SetClusterTimeRange(double mintime, double maxtime) { fClusterTimeRange.SetLimits(mintime, maxtime); }

  void SetCentralityRange(double min, double max) { fCentralityRange.SetLimits(min, max); fRequestCentrality = true; }

  void SetUserCentralityEstimator(TString centest) { fCentralityEstimator = centest; }

  void SetBunchCrossingIndex(Int_t bunchCrossingIndex) { fBunchCrossingIndex = bunchCrossingIndex; };

  static AliAnalysisTaskEmcalClustersRef *AddTaskEmcalClustersRef(const TString &nClusters = "usedefault", const TString &suffix = "");

  static AliAnalysisTaskEmcalClustersRef *AddTaskEmcalClustersRefDefault(const TString &nClusters = "usedefault");

protected:

  virtual void CreateUserObjects() {}

  virtual void CreateUserHistos();

  virtual bool IsUserEventSelected();

  virtual bool Run();

  virtual void UserFillHistosAfterEventSelection();

  void GetPatchBoundaries(TObject *o, Double_t *boundaries) const;

  void FillClusterHistograms(const TString &triggerclass, double energy, double transversenergy, double eta, double phi, double clustertime, int ncell, TList *triggerpatches);

  void FindPatchesForTrigger(TString triggerclass, const TClonesArray * triggerpatches, TList &foundpatches) const;

  Bool_t CorrelateToTrigger(Double_t etaclust, Double_t phiclust, TList *triggerpatches) const;


  AliCutValueRange<double>            fCentralityRange;           
  Bool_t                              fRequestCentrality;         
  Double_t                            fEventCentrality;           
  TString                             fCentralityEstimator;       
  Char_t                              fBunchCrossingIndex;        

  EnergyDefinition_t                  fEnergyDefinition;          
  Bool_t                              fEnableSumw2;               
  AliCutValueRange<double>            fClusterTimeRange;          

private:

  class EnergyBinning : public TCustomBinning {
  public:
    EnergyBinning();
    virtual ~EnergyBinning() {}
  };

  int CountTracklets(double etamin, double etamax, double phimin, double phimax);
  int CountEmcalClusters(double ecut);
  int GetEMCALCellOccupancy(double ecut);

  AliAnalysisTaskEmcalClustersRef(const AliAnalysisTaskEmcalClustersRef &);
  AliAnalysisTaskEmcalClustersRef &operator=(const AliAnalysisTaskEmcalClustersRef &);

  ClassDef(AliAnalysisTaskEmcalClustersRef, 1);
};

} /* namespace EMCalTriggerPtAnalysis */

#endif /* ALIANALYSISTASKEMCALCLUSTERSREF_H */