AliAnalysisTaskEmcalJetCDF.h

Go to the documentation of this file.

#ifndef ALIANALYSISTASKEMCALJETCDF_H
#define ALIANALYSISTASKEMCALJETCDF_H

/* Copyright(c) 1998-2016, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

#include "AliAnalysisTaskEmcalJet.h"
#include <THistManager.h>

class AliAnalysisTaskEmcalJetCDF : public AliAnalysisTaskEmcalJet
  {
  public:

    AliAnalysisTaskEmcalJetCDF();
    AliAnalysisTaskEmcalJetCDF ( const char *name );
    virtual ~AliAnalysisTaskEmcalJetCDF();

    void                        UserCreateOutputObjects();
    void                        Terminate ( Option_t *option );

    THistManager                fHistManager   ;

  protected:
    void                        ExecOnce();
    Bool_t                      Run() ;

    Bool_t                      FillHistograms()   ;

    TObject* GetHistogram ( const char* histName );

  private:
    AliAnalysisTaskEmcalJetCDF ( const AliAnalysisTaskEmcalJetCDF& );           // not implemented
    AliAnalysisTaskEmcalJetCDF &operator= ( const AliAnalysisTaskEmcalJetCDF& ); // not implemented

    ClassDef ( AliAnalysisTaskEmcalJetCDF, 7 );

  };

namespace NS_AliAnalysisTaskEmcalJetCDF {
  typedef std::pair<Double_t, Int_t> ptidx_pair;

  struct sort_descend
    {
    bool operator () ( const ptidx_pair &p1, const ptidx_pair &p2 )  { return p1.first > p2.first ; }
    };

  std::vector<Int_t>          SortTracksPt ( AliVEvent *event );

  std::vector<Int_t>          SortTracksPt ( AliParticleContainer *track_container );

  Double_t                    Z_ptot ( const AliEmcalJet* jet, const AliVParticle* trk ); // Get Z of constituent trk ; p total

  Double_t                    Z_pt   ( const AliEmcalJet* jet, const AliVParticle* trk ); // Get Z of constituent trk ; pt

  inline Double_t              Xi ( Double_t z ) { return TMath::Log ( 1/z ); } // Get Xi of value z

  Double_t                    DeltaR ( const AliVParticle* part1, const AliVParticle* part2 );

  Int_t                      IdxInArray ( Int_t index, TArrayI &array );

  TObject* AddTaskEmcalJetCDF (
                                const char* ntracks    = "usedefault",
                                const char* nclusters  = "usedefault",
                                const char* ncells     = "usedefault",
                                const char* tag        = "CDF"
                              );

  TObject* jetContSetParams (
                             AliJetContainer* jetCont,
                             Float_t jetptmin = 1.,
                             Float_t jetptmax = 500.,
                             Float_t jetareacutperc = 0.,
                             Int_t leadhadtype = 2,
                             Int_t nLeadJets = 1,
                             Float_t mintrackpt = 0.15,
                             Float_t maxtrackpt = 1000.
                            );

}

#endif

// kate: indent-mode none; indent-width 2; replace-tabs on;