AliPhysics  master (3d17d9d)
AliAnalysisTaskEmcalSoftDropResponse.h
Go to the documentation of this file.
1 /************************************************************************************
2  * Copyright (C) 2019, Copyright Holders of the ALICE Collaboration *
3  * All rights reserved. *
4  * *
5  * Redistribution and use in source and binary forms, with or without *
6  * modification, are permitted provided that the following conditions are met: *
7  * * Redistributions of source code must retain the above copyright *
8  * notice, this list of conditions and the following disclaimer. *
9  * * Redistributions in binary form must reproduce the above copyright *
10  * notice, this list of conditions and the following disclaimer in the *
11  * documentation and/or other materials provided with the distribution. *
12  * * Neither the name of the <organization> nor the *
13  * names of its contributors may be used to endorse or promote products *
14  * derived from this software without specific prior written permission. *
15  * *
16  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND *
17  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED *
18  * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE *
19  * DISCLAIMED. IN NO EVENT SHALL ALICE COLLABORATION BE LIABLE FOR ANY *
20  * DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES *
21  * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; *
22  * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND *
23  * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT *
24  * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS *
25  * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. *
26  ************************************************************************************/
27 #ifndef ALIANALYSISTASKEMCALSOFTDROPRESPONSE_H
28 #define ALIANALYSISTASKEMCALSOFTDROPRESPONSE_H
29 
30 #include <vector>
32 #include "THistManager.h"
33 
34 class RooUnfoldResponse;
35 class TBinning;
36 class TH2;
37 class TRandom;
38 
39 namespace PWGJE{
40 
41 namespace EMCALJetTasks {
42 
44 public:
49  };
51  kCAAlgo = 0,
52  kKTAlgo = 1,
54  };
55 
57  AliAnalysisTaskEmcalSoftDropResponse(const char *name);
59 
62  void SetBinningMode(EBinningMode_t binmode) { fBinningMode = binmode; }
65  void SetFractionResponseClosure(Double_t fracClosure) { fFractionResponseClosure = fracClosure; }
66  void SetBeta(double beta) { fBeta = beta; }
67  void SetZcut(double zcut) { fZcut = zcut; }
69  void SetSampleFraction(Double_t samplefraction) { fSampleFraction = samplefraction; }
70  void SetMinFractionShared(Double_t minsharedfraction) { fMinFractionShared = minsharedfraction; }
71  void SetUseChargedConstituents(bool doUse) { fUseChargedConstituents = doUse; }
72  void SetUseNeutralConstituents(bool doUse) { fUseNeutralConstituents = doUse; }
73  void SetNameMCParticleContainer(const char *name) { fNameMCParticles = name; }
74  void SetNamePartLevelJetContainer(const char *name) { fNamePartLevelJetContainer = name; }
75  void SetNameDetLevelJetContainer(const char *name) { fNameDetLevelJetContainer = name; }
77  void SetIsEmbeddedEvent(bool isEmbedded) {fIsEmbeddedEvent = isEmbedded; }
78 
79  static AliAnalysisTaskEmcalSoftDropResponse *AddTaskEmcalSoftDropResponse(Double_t jetradius, AliJetContainer::EJetType_t jettype, AliJetContainer::ERecoScheme_t recombinationScheme, bool ifembed, const char *namepartcont, const char *trigger);
80 
81 protected:
82  virtual void UserCreateOutputObjects();
83  virtual Bool_t CheckMCOutliers();
84  virtual bool Run();
85 
88  TBinning *GetZgBinning() const;
89  TBinning *GetRgBinning(double R) const;
90 
91  std::vector<double> MakeSoftdrop(const AliEmcalJet &jet, double jetradius, const AliParticleContainer *tracks, const AliClusterContainer *clusters) const;
92 
93 private:
94 
100  Double32_t fSampleFraction;
102 
108  TRandom *fSampleSplitter;
109  TRandom *fSampleTrimmer;
116  std::vector<RooUnfoldResponse*> fZgResponse;
117  std::vector<RooUnfoldResponse*> fZgResponseClosure;
118  std::vector<RooUnfoldResponse*> fRgResponse;
119  std::vector<RooUnfoldResponse*> fRgResponseClosure;
120  std::vector<RooUnfoldResponse*> fNsdResponse;
121  std::vector<RooUnfoldResponse*> fNsdResponseClosure;
122  std::vector<RooUnfoldResponse*> fThetagResponse;
123  std::vector<RooUnfoldResponse*> fThetagResponseClosure;
125 
128 
130 
131 };
132 
133 }
134 
135 }
136 
137 
138 #endif
Bool_t fHasResponseMatrixSparse
Fill also THnSparse representation of response matrix.
double Double_t
Definition: External.C:58
std::vector< RooUnfoldResponse * > fNsdResponse
! RooUnfold response for n_sd
Float_t fMinFractionShared
only fill histos for jets if shared fraction larger than X
AliAnalysisTaskEmcalSoftDropResponse & operator=(const AliAnalysisTaskEmcalSoftDropResponse &)
std::vector< RooUnfoldResponse * > fThetagResponseClosure
! RooUnfold response for n_sd for the closure test
TString fNameDetLevelJetContainer
Name of the detector (or hybrid if embedding) level jet container.
std::vector< RooUnfoldResponse * > fZgResponseClosure
! RooUnfold response for z_g for the closure test
Interface for binnings used by the histogram handler.
Definition: TBinning.h:23
Container for particles within the EMCAL framework.
Double_t fFractionResponseClosure
Fraction of events used for the response matrix in the closure test.
std::vector< RooUnfoldResponse * > fNsdResponseClosure
! RooUnfold response for n_sd for the closure test
virtual Bool_t CheckMCOutliers()
Filter the mc tails in pt-hard distributions.
float Float_t
Definition: External.C:68
std::vector< RooUnfoldResponse * > fRgResponse
! RooUnfold response for r_g
static AliAnalysisTaskEmcalSoftDropResponse * AddTaskEmcalSoftDropResponse(Double_t jetradius, AliJetContainer::EJetType_t jettype, AliJetContainer::ERecoScheme_t recombinationScheme, bool ifembed, const char *namepartcont, const char *trigger)
virtual bool Run()
Run function. This is the core function of the analysis and contains the user code. Therefore users have to implement this function.
std::vector< RooUnfoldResponse * > fThetagResponse
! RooUnfold response for theta_g
Definition: External.C:220
TString fNameUnSubLevelJetContainer
Name of the unsubtracted hybrid level jet container.
Base task in the EMCAL jet framework.
Represent a jet reconstructed using the EMCal jet framework.
Definition: AliEmcalJet.h:51
Container class for histograms.
Definition: THistManager.h:99
std::vector< RooUnfoldResponse * > fRgResponseClosure
! RooUnfold response for r_g for the closure test
Double32_t fSampleFraction
Fraction of statistics used for the analysis.
TString fNamePartLevelJetContainer
Name of the particle level jet container.
bool Bool_t
Definition: External.C:53
std::vector< RooUnfoldResponse * > fZgResponse
! RooUnfold response for z_g
Container structure for EMCAL clusters.
std::vector< double > MakeSoftdrop(const AliEmcalJet &jet, double jetradius, const AliParticleContainer *tracks, const AliClusterContainer *clusters) const