AliPhysics  59e0e03 (59e0e03)
AliAnalysisTaskSoftDropResponse.h
Go to the documentation of this file.
1 #ifndef ALIANALYSISTASKSOFTDROPRESPONSE_H
2 #define ALIANALYSISTASKSOFTDROPRESPONSE_H
3 /**************************************************************************
4 * Copyright(c) 1998-2009, ALICE Experiment at CERN, All rights reserved. *
5 * *
6 * Author: The ALICE Off-line Project. *
7 * Contributors are mentioned in the code where appropriate. *
8 * *
9 * Permission to use, copy, modify and distribute this software and its *
10 * documentation strictly for non-commercial purposes is hereby granted *
11 * without fee, provided that the above copyright notice appears in all *
12 * copies and that both the copyright notice and this permission notice *
13 * appear in the supporting documentation. The authors make no claims *
14 * about the suitability of this software for any purpose. It is *
15 * provided "as is" without express or implied warranty. *
16 **************************************************************************/
17 
18 //-----------------------------------------------------------------------
19 // Author : Kirill Lapidus, Yale University, kirill.lapidus@cern.ch
20 //-----------------------------------------------------------------------
21 
22 #include "AliJetResponseMaker.h"
23 #include "AliEmcalEmbeddingQA.h"
24 
25 #include "FJ_includes.h"
26 
28  public:
30  AliAnalysisTaskSoftDropResponse(const char *name);
32 
34 
36  const char *ntracks1 = "Tracks",
37  const char *nclusters1 = "CaloClusters",
38  const char *njets1 = "Jets",
39  const char *nrho1 = "Rho",
40  const Double_t jetradius1 = 0.4,
41  const char *ntracks2 = "MCParticles",
42  const char *nclusters2 = "",
43  const char *njets2 = "MCJets",
44  const char *nrho2 = "",
45  const Double_t jetradius2 = 0.4,
46  const Double_t jetptcut = 1,
47  const Double_t jetareacut = 0.557,
48  const Double_t jetBias = 5,
49  const Int_t biasType = 0, // 0 = charged, 1 = neutral, 2 = both
51  const Double_t maxDistance1 = 0.25,
52  const Double_t maxDistance2 = 0.25,
53  const char *cutType = "TPC",
54  const Int_t ptHardBin = -999,
55  const Double_t minCent = -999,
56  const Double_t maxCent = -999,
57  const char *taskname = "AliAnalysisTaskSoftDropResponse",
58  const Bool_t biggerMatrix = kFALSE,
60  const Double_t maxTrackPt = 100);
61 
62  void SetZ2gAxis(Int_t b) { fZ2gAxis = b; }
63 
64  protected:
65 
67  void FillMatchingHistos(AliEmcalJet* jet1, AliEmcalJet* jet2, Double_t d, Double_t CE1, Double_t CE2);
68  void AllocateTHnSparse();
69 
70  void CalculateZg(AliEmcalJet* jet, const Float_t zcut, const Float_t beta, Float_t zg);
71 
73 
74  private:
75 
76  void Decluster(const fastjet::PseudoJet& jet);
77  fastjet::ClusterSequence* Recluster(const AliEmcalJet* jet);
78 
79  void FillZgRgVectors(const AliEmcalJet* jet);
80 
82  std::vector<Float_t> fZg_values;
83  std::vector<Float_t> fRg_values;
84  std::vector<Int_t> fSDsteps_values;
85 
88 
90 };
91 #endif
Bool_t FillHistograms()
Function filling histograms.
double Double_t
Definition: External.C:58
std::vector< Float_t > fRg_values
! n_SD ordering
AliAnalysisTaskSoftDropResponse & operator=(const AliAnalysisTaskSoftDropResponse &)
int Int_t
Definition: External.C:63
std::vector< Float_t > fZg_values
! n_SD ordering
float Float_t
Definition: External.C:68
void CalculateZg(AliEmcalJet *jet, const Float_t zcut, const Float_t beta, Float_t zg)
void Decluster(const fastjet::PseudoJet &jet)
static AliAnalysisTaskSoftDropResponse * AddTaskSoftDropResponse(const char *ntracks1="Tracks", const char *nclusters1="CaloClusters", const char *njets1="Jets", const char *nrho1="Rho", const Double_t jetradius1=0.4, const char *ntracks2="MCParticles", const char *nclusters2="", const char *njets2="MCJets", const char *nrho2="", const Double_t jetradius2=0.4, const Double_t jetptcut=1, const Double_t jetareacut=0.557, const Double_t jetBias=5, const Int_t biasType=0, const AliAnalysisTaskSoftDropResponse::MatchingType matching=AliAnalysisTaskSoftDropResponse::kGeometrical, const Double_t maxDistance1=0.25, const Double_t maxDistance2=0.25, const char *cutType="TPC", const Int_t ptHardBin=-999, const Double_t minCent=-999, const Double_t maxCent=-999, const char *taskname="AliAnalysisTaskSoftDropResponse", const Bool_t biggerMatrix=kFALSE, AliAnalysisTaskSoftDropResponse *address=0, const Double_t maxTrackPt=100)
Represent a jet reconstructed using the EMCal jet framework.
Definition: AliEmcalJet.h:51
bool Bool_t
Definition: External.C:53
Int_t fZ2gAxis
add Z2g axis in matching THnSparse (default=0)
void FillMatchingHistos(AliEmcalJet *jet1, AliEmcalJet *jet2, Double_t d, Double_t CE1, Double_t CE2)
fastjet::ClusterSequence * Recluster(const AliEmcalJet *jet)