AliAnalysisTaskEmcalJetSpectra8TeVTriggerQA.h

Go to the documentation of this file.

#ifndef ALIANALYSISTASKEMCALJETSPECTRA8TEVTRIGGERQA_H
#define ALIANALYSISTASKEMCALJETSPECTRA8TEVTRIGGERQA_H

/**********************************************************************************
 * Copyright (C) 2018, Copyright Holders of the ALICE Collaboration                *
 * All rights reserved.                                                            *
 *                                                                                 *
 * Redistribution and use in source and binary forms, with or without              *
 * modification, are permitted provided that the following conditions are met:     *
 *   * Redistributions of source code must retain the above copyright              *
 *     notice, this list of conditions and the following disclaimer.               *
 *   * Redistributions in binary form must reproduce the above copyright           *
 *     notice, this list of conditions and the following disclaimer in the         *
 *     documentation and/or other materials provided with the distribution.        *
 *   * Neither the name of the <organization> nor the                              *
 *     names of its contributors may be used to endorse or promote products        *
 *     derived from this software without specific prior written permission.       *
 *                                                                                 *
 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND *
 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED   *
 * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE          *
 * DISCLAIMED. IN NO EVENT SHALL ALICE COLLABORATION BE LIABLE FOR ANY             *
 * DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES      *
 * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;    *
 * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND     *
 * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT      *
 * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS   *
 * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.                    *
 * *********************************************************************************/

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



#include "AliAnalysisTaskEmcalJet.h"
#include "THistManager.h"
#include <string>
#include <vector>
#include "Tracks/AliAnalysisTaskEmcalTriggerBase.h"
#include "Tracks/AliCutValueRange.h"
#include "AliEventCuts.h"
#include <TCustomBinning.h>
#include <TString.h>

class AliOADBContainer;
class AliEMCALTriggerPatchInfo;
class THistManager;
class TObjArray;
class TFormula;
class TH2F;
class TH1F;
class TF1;
class THnSparse;
class TRandom3;
class TObjArray;
class TClonesArray;
class TObject;
class TString;
class TProfile2D;
class AliAODEvent;
class AliESDEvent;
class AliMCEvent;
class AliVEvent;
class AliStack;
class AliEMCALGeometry;
//class AliEMCalJet;
class AliEMCALRecoUtils;
class AliESDCaloCluster;
class AliVTrack;
class AliMCEvent;
class AliAODEvent;
class AliESDEvent;

class AliAnalysisTaskEmcalJetSpectra8TeVTriggerQA : public AliAnalysisTaskEmcalJet {
public:

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

    void                        UserCreateOutputObjects()                         ;
    void                        Terminate(Option_t *option)                       ;
    void                        ExtractMainPatch()                                ;
    
    //setters
    void                        SetUseSumw2(Bool_t b) { fUseSumw2          = b   ; }

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

    void                        AllocateJetHistograms()                           ;
    void                        AllocateTrackHistograms()                         ;
    void                        AllocateClusterHistograms()                       ;
    void                        AllocateCellHistograms()                          ;
    void                        AllocateParticleHistograms()                      ;
    //void                        AllocateMCJetHistograms()                         ;///< MC truth Jet Histograms
    
    

    void                        DoJetLoop()                                       ;
    void                        DoTrackLoop()                                     ;
    void                        DoClusterLoop()                                   ;
    void                        DoCellLoop()                                      ;
    void                        DoParticleLoop()                                  ;
    //void                        DoMCJetLoop()                                     ;
    //Bool_t                      IsLEDEvent() const                                ;
    
    Bool_t                      fUseRecalcPatches                                 ;
    Bool_t                      SelectSingleShowerPatch(const AliEMCALTriggerPatchInfo *patch) const;
    Bool_t                      SelectJetPatch(const AliEMCALTriggerPatchInfo  *patch) const;
    THistManager                fHistManager                                      ;
    
    // Sparse Definition
    virtual THnSparse*     NewTHnSparseF(const char* name, UInt_t entries);
    virtual void           GetDimParams(Int_t iEntry,TString &label, Int_t &nbins, Double_t &xmin, Double_t &xmax);
    // Binning helper functions
    
    //  Utilities
    //Double_t                    GetZ(const Double_t trkPx, const Double_t trkPy, const Double_t trkPz, const Double_t jetPx, const Double_t jetPy, const Double_t jetPz) const;
    //Double_t                    GetSmearedTrackPt(AliVTrack *track)                ; //!<!               smeared Pt resolution with gaussian
    //Double_t                    GetFcross(const AliVCluster *cluster, AliVCaloCells *cells);
    
    //AliVEvent                   *fRecevent                                         ;//!<!                Reconstructed event
    //AliMCEvent                  *fMCevent                                          ;//!<!                Monte-Carlo event
    //AliMCParticleContainer*     fGeneratorLevel                                    ;//!<!                generator level container
   // AliJetContainer*            fMCJetContainer                                    ;//!<!                truth-level jet container

private:
    
  //AliEMCALGeometry              *fGeom                     ;//!<!          EMCal goemetry utility
   // AliEMCALRecoUtils             *fRecoUtil                 ;//!<!          Reco utility
    //TF1                           *fClusterEResolution       ;//!<!          Parameterization of cluster energy resolution from 2010 test beam results a = 4.35 b = 9.07 c = 1.63
  //  Double_t                      fVaryTrkPtRes              ;//!<!          Variation of tracking momentum resolution
    Bool_t                        fUseSumw2                  ;

    TH1F                          *fHistNumbJets             ;
    TH1F                          *fHistJetPt                ;
    TH1F                          *fHistJetJetPatchE         ;
    TH1F                          *fHistJetGammaPatchE       ;
    TH1F                          *fHistJetJetPatchPt        ;
    TH1F                          *fHistJetGammaPatchPt      ;
    TH1F                          *fHistTriggerPatchE        ;
    
    //THnSparse                     *fhnMBJetSpectra           ;//!<!          MB Jet Spectra sparse
    //THnSparse                     *fhnTrkQA                  ;//!<!          Charged Track QA Sparse
   // THnSparse                     *fhnClusQA                 ;//!<!          EMCal Cluster QA
    
    
    //TH1F                          *fHistEMCalTowerMult[9]    ;//!<!          EMCal Tower Multiplicity by SM

    AliAnalysisTaskEmcalJetSpectra8TeVTriggerQA(const AliAnalysisTaskEmcalJetSpectra8TeVTriggerQA&)           ; // not implemented
    AliAnalysisTaskEmcalJetSpectra8TeVTriggerQA &operator=(const AliAnalysisTaskEmcalJetSpectra8TeVTriggerQA&); // not implemented

    ClassDef(AliAnalysisTaskEmcalJetSpectra8TeVTriggerQA, 7);
};
#endif