AliPhysics  fffcdf3 (fffcdf3)
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
AliAnalysisTaskEMCALClusterize.h
Go to the documentation of this file.
1 #ifndef ALIANALYSISTASKEMCALCLUSTERIZE_H
2 #define ALIANALYSISTASKEMCALCLUSTERIZE_H
3 
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5  * See cxx source for full Copyright notice */
6 
7 //_________________________________________________________________________
18 //_________________________________________________________________________
19 
20 // Root
21 class TTree;
22 class TClonesArray;
23 #include <TRandom3.h>
24 
25 // EMCAL
26 class AliEMCALGeometry;
27 class AliEMCALClusterizer;
28 class AliEMCALAfterBurnerUF;
29 class AliEMCALRecPoint;
30 class AliAODCaloCluster;
31 class AliCentrality;
32 class AliMultSelection;
33 class AliVCaloCells;
34 
35 #include "AliEMCALRecParam.h"
36 #include "AliEMCALRecoUtils.h"
37 
38 #include "AliAnalysisTaskSE.h"
39 
41 
42  public:
43 
45  AliAnalysisTaskEMCALClusterize(const char *name);
47 
48  virtual void UserCreateOutputObjects();
49  virtual void UserExec(Option_t *option);
50  virtual void Init();
51  virtual void LocalInit() { Init() ; }
52  void PrintParam();
53 
54  // Event methods, settings
55 
56  Bool_t AcceptCell(Int_t absID);
61  void SetEMCALNcellsCut(Int_t cut) { fEMCALNcellsCut = cut ; }
62 
65 
66  void CheckAndGetEvent();
67 
71 
72  Bool_t IsLEDEvent(const Int_t run);
75 
76  // OCDB, avoid acessing OCDB!
77 
78  void AccessOCDB();
79  void SwitchOnAccessOCDB() { fAccessOCDB = kTRUE ; }
80  void SwitchOffAccessOCDB() { fAccessOCDB = kFALSE ; }
81  void SetOCDBPath(const char *path) { fOCDBpath = path ; }
82 
83  // Geometry methods
84 
85  void InitGeometry();
86  void SetGeometryName(TString &name) { fGeomName = name ; }
87  TString GeometryName() const { return fGeomName ; }
90  void SetGeometryMatrixInSM(TGeoHMatrix* m, Int_t i) { fGeomMatrix[i] = m ; }
91 
93  TString pa = "") { fImportGeometryFromFile = im ;
94  fImportGeometryFilePath = pa ; }
95  // Outout AOD branch methods
96 
98  void SetAODCellsName(TString &name) { fOutputAODCellsName = name ; }
99  void FillAODFile(Bool_t yesno) { fFillAODFile = yesno ; }
100  void FillAODCaloCells();
101  void FillAODHeader();
103  void SwitchOffFillAODHeader() { fFillAODHeader = kFALSE ; }
106 
109 
110  // Algorithms settings
111 
112  AliEMCALRecParam * GetRecParam() { if(!fRecParam) fRecParam = new AliEMCALRecParam ;
113  return fRecParam ; }
114 
115  AliEMCALRecoUtils* GetRecoUtils() { if(!fRecoUtils) fRecoUtils = new AliEMCALRecoUtils ;
116  return fRecoUtils ; }
117  void ConfigureEMCALRecoUtils(Bool_t bMC = kFALSE, Bool_t bExotic= kTRUE, Bool_t bNonLin= kFALSE,
118  Bool_t bRecalE= kTRUE , Bool_t bBad = kTRUE, Bool_t bRecalT= kTRUE, Int_t debug = -1);
119 
120  void InitClusterization();
121  void ClusterizeCells();
122  void ClusterUnfolding();
123  void JustUnfold(Bool_t yesno) { fJustUnfold = yesno ; }
124  void UpdateCells();
125 
126  void SetConfigFileName(TString name) { fConfigName = name ; }
127  void SetMaxEvent(Int_t max) { fMaxEvent = max ; }
128 
131 
132  void SwitchOnUpdateCell() { fUpdateCell = kTRUE ; }
133  void SwitchOffUpdateCell() { fUpdateCell = kFALSE ; }
134 
135  // Cell selection after unfolding
136 
140  fSelectCellMinFrac = frac ; }
142 
143  // OADB options settings
144 
145  void AccessOADB() ;
146 
147  TString GetPass() ;
148 
149  void SwitchOnEMCALOADB() { fAccessOADB = kTRUE ; }
150  void SwitchOffEMCALOADB() { fAccessOADB = kFALSE ; }
151 
152  void SetOADBFilePath(TString path) { fOADBFilePath = path ; }
153 
155 
156  // Centrality selection
157 
158  AliCentrality* GetCentrality() const { return InputEvent()->GetCentrality() ; }
159  AliMultSelection* GetMultSelCen() const { return (AliMultSelection * ) InputEvent()->FindListObject("MultSelection") ; }
162 
165  Float_t GetEventCentrality() const ;
166  void SetCentralityBin(Int_t min, Int_t max) //Set the centrality bin to select the event. If used, then need to get percentile
167  { fCentralityBin[0]=min ; fCentralityBin[1]=max ; }
168  Float_t GetCentralityBin(Int_t i) const { if(i < 0 || i > 1) return -1 ;
169  else return fCentralityBin[i] ; }
170 
171  // MC label properly assignation methods
172 
173  void RemapMCLabelForAODs(Int_t &label);
176 
177  void SetClustersMCLabelFrom2SelectedLabels(AliEMCALRecPoint* recPoint, AliAODCaloCluster *clus) ;
178  void SetClustersMCLabelFromOriginalClusters(AliAODCaloCluster * clus) ;
179 
182 
186 
187  //-----------------------------------------
188  // T-Card correlation emulation, do on MC
189 
190  void MakeCellTCardCorrelation() ;
192 
195  void SwitchOnTCardCorrelation(Bool_t conservEnergy = kFALSE) { fTCardCorrEmulation = kTRUE ; fTCardCorrClusEnerConserv = conservEnergy ; }
196 
199 
200 
206  for(Int_t ism = 0; ism < 22; ism++) {
207  fTCardCorrInduceEner[0][ism] = ud; fTCardCorrInduceEner[1][ism] = udlr;
208  fTCardCorrInduceEner[2][ism] = lr; fTCardCorrInduceEner[3][ism] = sec ; } }
209 
215  for(Int_t ism = 0; ism < 22; ism++) {
216  fTCardCorrInduceEnerFrac[0][ism] = ud; fTCardCorrInduceEnerFrac[1][ism] = udlr;
217  fTCardCorrInduceEnerFrac[2][ism] = lr; fTCardCorrInduceEnerFrac[3][ism] = sec ; } }
218 
224  for(Int_t ism = 0; ism < 22; ism++) {
225  fTCardCorrInduceEnerFracP1[0][ism] = ud; fTCardCorrInduceEnerFracP1[1][ism] = udlr;
226  fTCardCorrInduceEnerFracP1[2][ism] = lr; fTCardCorrInduceEnerFracP1[3][ism] = sec ; } }
227 
234  if ( sm < 22 && sm >= 0 ) {
235  fTCardCorrInduceEner[0][sm] = ud; fTCardCorrInduceEner[1][sm] = udlr;
236  fTCardCorrInduceEner[2][sm] = lr; fTCardCorrInduceEner[3][sm] = sec ; } }
237 
244  if ( sm < 22 && sm >= 0 ) {
245  fTCardCorrInduceEnerFrac[0][sm] = ud; fTCardCorrInduceEnerFrac[1][sm] = udlr;
246  fTCardCorrInduceEnerFrac[2][sm] = lr; fTCardCorrInduceEnerFrac[3][sm] = sec ; } }
247 
254  if ( sm < 22 && sm >= 0 ) {
255  fTCardCorrInduceEnerFracP1[0][sm] = ud; fTCardCorrInduceEnerFracP1[1][sm] = udlr;
256  fTCardCorrInduceEnerFracP1[2][sm] = lr; fTCardCorrInduceEnerFracP1[3][sm] = sec ; } }
257 
258 
264  for(Int_t ism = 0; ism < 22; ism++) {
265  fTCardCorrInduceEnerFracWidth[0][ism] = ud; fTCardCorrInduceEnerFracWidth[1][ism] = udlr;
266  fTCardCorrInduceEnerFracWidth[2][ism] = lr; fTCardCorrInduceEnerFracWidth[3][ism] = sec ; } }
267 
274  if ( sm < 22 && sm >= 0 ) {
276  fTCardCorrInduceEnerFracWidth[2][sm] = lr; fTCardCorrInduceEnerFracWidth[3][sm] = sec ; } }
277 
282  if ( sm < 22 && sm >= 0 ) fTCardCorrInduceEnerProb[sm] = prob ; }
283 
286 
289 
290  void PrintTCardParam();
291 
292  //------------------------------------------
293 
294 private:
295 
296  virtual void FillCaloClusterInEvent();
297 
298  virtual void RecPoints2Clusters();
299 
300  virtual void ResetArrays();
301 
302  AliVEvent *fEvent;
303 
304  // Geometry
305  AliEMCALGeometry *fGeom;
307  TGeoHMatrix *fGeomMatrix[22];
310 
311  // OCDB
314 
315  // Temporal arrays
316  TClonesArray *fDigitsArr;
319  AliVCaloCells *fCaloCells;
320 
321  // Clusterizers
322  AliEMCALRecParam *fRecParam;
323  AliEMCALClusterizer *fClusterizer;
324  AliEMCALAfterBurnerUF *fUnfolder;
326 
327  // AOD
328  TClonesArray *fOutputAODBranch;
330  AliAODCaloCells *fOutputAODCells;
337 
339 
340  AliEMCALRecoUtils* fRecoUtils;
342 
343  static const Int_t fgkNEMCalCells = 17664;
344 
351 
353 
355 
364 
367 
372 
373  // Centrality
377 
378  // Event selection with some signal in EMCAL
382 
388 
392 
394 
396 
397 
398  // T-Card correlation emulation, do on MC
403 
409 
410  TRandom3 fRandom ;
412 
415 
417 
420 
423 
425  ClassDef(AliAnalysisTaskEMCALClusterize, 40) ;
427 
428 };
429 
430 #endif //ALIANALYSISTASKEMCALCLUSTERIZE_H
TString fOutputAODCellsName
New of output cells AOD branch.
void FillAODCaloCells()
Put calo cells in standard branch.
Float_t fTCardCorrMinAmp
Minimum cell energy to induce signal on adjacent cells.
AliVCaloCells * fCaloCells
! CaloCells container
Bool_t fTCardCorrEmulation
Activate T-Card cells energy correlation.
Float_t fTCardCorrInduceEnerFracP1[4][22]
Induced energy loss gauss fraction param1 on 0-same row, diff col, 1-up/down cells left/right col 2-l...
double Double_t
Definition: External.C:58
virtual void ResetArrays()
Reset arrays containing information for all possible cells.
const int debug
Definition: scanAll.C:15
void SetCellCuts(Float_t e, Float_t frac)
TString GetPass()
Get or guess pass number/string from path of filename.
Bool_t fFillAODHeader
Copy header to standard branch.
TClonesArray * fDigitsArr
! Digits array
Bool_t fTCardCorrClusEnerConserv
When making correlation, subtract from the reference cell the induced energy on the neighbour cells...
void SetInducedEnergyLossFraction(Float_t ud, Float_t udlr, Float_t lr, Float_t sec)
AliEMCALClusterizer * fClusterizer
! EMCAL clusterizer
Bool_t fJustUnfold
Just unfold, do not recluster.
Bool_t fImportGeometryFromFile
Import geometry settings in geometry.root file.
Bool_t fSelectCell
Reject cells from cluster if energy is too low and recalculate position/energy and other...
TString fConfigName
Name of analysis configuration file.
void SetInducedEnergyLossFractionWidth(Float_t ud, Float_t udlr, Float_t lr, Float_t sec)
void PrintTCardParam()
Print parameters for T-Card correlation emulation.
Int_t fCellSecondLabels[fgkNEMCalCells]
Array with Second MC label to be passed to digit.
void SwitchOffTCardCorrelation()
De-activate T-Card cells correlation,.
void SetClustersMCLabelFrom2SelectedLabels(AliEMCALRecPoint *recPoint, AliAODCaloCluster *clus)
Float_t fCellMatchdEta[fgkNEMCalCells]
Array with cluster-track dPhi.
TClonesArray * fOutputAODBranch
! AOD Branch with output clusters
TString fOADBFilePath
Default path $ALICE_PHYSICS/OADB/EMCAL, if needed change.
Float_t fTCardCorrCellsEner[fgkNEMCalCells]
Array with induced cell energy in T-Card neighbour cells.
TString fImportGeometryFilePath
path fo geometry.root file
Float_t fTCardCorrInduceEnerFrac[4][22]
Induced energy loss gauss fraction param0 on 0-same row, diff col, 1-up/down cells left/right col 2-l...
void ConfigureEMCALRecoUtils(Bool_t bMC=kFALSE, Bool_t bExotic=kTRUE, Bool_t bNonLin=kFALSE, Bool_t bRecalE=kTRUE, Bool_t bBad=kTRUE, Bool_t bRecalT=kTRUE, Int_t debug=-1)
Bool_t fPrintOnce
Print once analysis parameters.
Bool_t fSelectEMCALEvent
Process the event if there is some high energy cluster.
Float_t fTCardCorrMaxInduced
Maximum induced energy signal on adjacent cells.
Float_t GetEventCentrality() const
Get centrality/multiplicity percentile.
Bool_t fInputFromFilter
Get the input from AODs from the filter.
TString fOutputAODBranchName
New of output clusters AOD branch.
AliEMCALGeometry * fGeom
EMCAL geometry.
void ClusterUnfolding()
Take the event clusters and unfold them.
int Int_t
Definition: External.C:63
Bool_t fAccessOADB
Get calibration from OADB for EMCAL.
Bool_t fRemoveLEDEvents
Remove LED events, use only for LHC11a.
void SetInducedEnergyLossConstantPerSM(Int_t sm, Float_t ud, Float_t udlr, Float_t lr, Float_t sec)
AliEMCALRecoUtils * fRecoUtils
Access to factorized reconstruction algorithms.
float Float_t
Definition: External.C:68
Bool_t fRejectBelowThreshold
split (false-default) or reject (true) cell energy below threshold after UF
Float_t fConstantTimeShift
Apply a 600 ns time shift in case of simulation, shift in ns.
void SetClustersMCLabelFromOriginalClusters(AliAODCaloCluster *clus)
Bool_t fUpdateCell
On/Off the upate of the CaloCells container.
Float_t fSelectCellMinE
Min energy cell threshold, after unfolding.
Bool_t fOADBSet
AODB parameters already set.
AliAODCaloCells * fOutputAODCells
! AOD Branch with output cells
Float_t fTCardCorrInduceEner[4][22]
Induced energy loss gauss constant on 0-same row, diff col, 1-up/down cells left/right col 2-left/rig...
Int_t fCellLabels[fgkNEMCalCells]
Array with MC label to be passed to digit.
Bool_t fRemoveExoticEvents
Remove exotic events.
void SwitchOnTCardCorrelation(Bool_t conservEnergy=kFALSE)
Bool_t fDoTrackMatching
On/Off the matching recalculation to speed up analysis in PbPb.
void SetInducedEnergyLossFractionP1PerSM(Int_t sm, Float_t ud, Float_t udlr, Float_t lr, Float_t sec)
Float_t fCellMatchdPhi[fgkNEMCalCells]
Array with cluster-track dEta.
Float_t fTCardCorrInduceEnerProb[22]
Probability to induce energy loss per SM.
TGeoHMatrix * fGeomMatrix[22]
Geometry matrices with alignments.
Bool_t IsLEDEvent(const Int_t run)
Check if event is LED, is so remove it. Affected LHC11a runs.
AliEMCALRecParam * fRecParam
Reconstruction parameters container.
TString fGeomName
Name of geometry to use.
Bool_t fOutputAODBranchSet
Set the AOD clusters branch in the input event once.
void PrintParam()
Print clusterization task parameters.
TString fOCDBpath
Path with OCDB location.
AliEMCALAfterBurnerUF * fUnfolder
! Unfolding procedure
Bool_t fGeomMatrixSet
Set geometry matrices only once, for the first event.
Float_t fCentralityBin[2]
Minimum and maximum value of the centrality for the analysis.
Bool_t fFillAODCaloCells
Copy calocells to standard branch.
TObjArray * fCaloClusterArr
! CaloClusters array
Double_t fCellTime[fgkNEMCalCells]
Array with cluster time to be passed to digit in case of AODs.
TString fCentralityClass
Name of selected centrality class.
void SetImportGeometryFromFile(Bool_t im, TString pa="")
void SetInducedEnergyLossFractionP1(Float_t ud, Float_t udlr, Float_t lr, Float_t sec)
Bool_t fAccessOCDB
Need to access info from OCDB (not really)
void SetInducedEnergyLossProbabilityPerSM(Float_t prob, Int_t sm)
AliAnalysisTaskEMCALClusterize & operator=(const AliAnalysisTaskEMCALClusterize &)
Assignment operator not implemented.
TObjArray * fClusterArr
! Recpoints array
Bool_t fRecalibrateWithClusterTime
Use fCellTime to store time of cells in cluster.
const char Option_t
Definition: External.C:48
void SetInducedEnergyLossFractionWidthPerSM(Int_t sm, Float_t ud, Float_t udlr, Float_t lr, Float_t sec)
void SetInducedEnergyLossConstant(Float_t ud, Float_t udlr, Float_t lr, Float_t sec)
bool Bool_t
Definition: External.C:53
void FillAODHeader()
Put event header information in standard AOD branch.
static const Int_t fgkNEMCalCells
Total number of cells in the calorimeter, 10*48*24 (EMCal) + 4*48*8 (EMCal/DCal 1/3) + 6*32*24 (DCal)...
Int_t fEMCALNcellsCut
At least an EMCAL cluster with fNCellsCut cells over fEnergyCut.
void SetGeometryMatrixInSM(TGeoHMatrix *m, Int_t i)
Reclusterize EMCal clusters, put them in a new branch for other following analysis.
Bool_t fTCardCorrCellsNew[fgkNEMCalCells]
Array with induced cell energy in T-Card neighbour cells, that before had no signal.
Bool_t fLoadGeomMatrices
Matrices set from configuration, not get from geometry.root or from ESDs/AODs.
Float_t fEMCALEnergyCut
At least an EMCAL cluster with this energy in the event.
Bool_t fRandomizeTCard
Use random induced energy.
Bool_t fUseAliCentrality
Use the centrality estimator from AliCentrality or AliMultSelection.
void SetInducedEnergyLossFractionPerSM(Int_t sm, Float_t ud, Float_t udlr, Float_t lr, Float_t sec)
Float_t fTCardCorrInduceEnerFracWidth[4][22]
Induced energy loss gauss witdth on 0-same row, diff col, 1-up/down cells left/right col 2-left/righ ...
Int_t fOrgClusterCellId[fgkNEMCalCells]
Array ID of cluster to wich the cell belongs in unmodified clusters.
Float_t fSelectCellMinFrac
Min fraction of cell energy after unfolding cut.
Bool_t reject
Bool_t fRemapMCLabelForAODs
Remap AOD cells MC label. Needed in old AOD productions.