AliPhysics  a34469b (a34469b)
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
AliAnalysisTaskEmcalTriggerBase.cxx
Go to the documentation of this file.
1 /**************************************************************************
2  * Copyright(c) 1998-2016, ALICE Experiment at CERN, All rights reserved. *
3  * *
4  * Author: The ALICE Off-line Project. *
5  * Contributors are mentioned in the code where appropriate. *
6  * *
7  * Permission to use, copy, modify and distribute this software and its *
8  * documentation strictly for non-commercial purposes is hereby granted *
9  * without fee, provided that the above copyright notice appears in all *
10  * copies and that both the copyright notice and this permission notice *
11  * appear in the supporting documentation. The authors make no claims *
12  * about the suitability of this software for any purpose. It is *
13  * provided "as is" without express or implied warranty. *
14  **************************************************************************/
15 #include <array>
16 #include <functional>
17 #include <iostream>
18 
19 #include <TClonesArray.h>
20 #include <TGrid.h>
21 #include <THistManager.h>
22 #include <TParameter.h>
23 
25 #include "AliAnalysisUtils.h"
26 #include "AliEMCALGeometry.h"
27 #include "AliEMCALTriggerPatchInfo.h"
28 #include "AliEMCALTriggerMapping.h"
29 #include "AliESDEvent.h"
30 #include "AliInputEventHandler.h"
31 #include "AliLog.h"
32 #include "AliOADBContainer.h"
33 #include "AliVVertex.h"
34 
38 
39 namespace EMCalTriggerPtAnalysis {
40 
41 AliAnalysisTaskEmcalTriggerBase::AliAnalysisTaskEmcalTriggerBase():
43  fTriggerSelection(nullptr),
44  fUseTriggerBits(kTRUE),
45  fRequireBunchCrossing(kTRUE),
46  fHistos(nullptr),
47  fTriggerStringFromPatches(kFALSE),
48  fSelectedTriggers(),
49  fRequireAnalysisUtils(kTRUE),
50  fVertexCut(-10., 10.),
51  fNameDownscaleOADB(""),
52  fDownscaleOADB(nullptr),
53  fDownscaleFactors(nullptr),
54  fNameMaskedFastorOADB(),
55  fMaskedFastorOADB(nullptr),
56  fMaskedFastors(),
57  fOnlineTriggerThresholds(),
58  fNameAcceptanceOADB(),
59  fSelectNoiseEvents(false),
60  fRejectNoiseEvents(false),
61  fEnableDCALTriggers(true)
62 {
63  SetNeedEmcalGeom(true);
65  SetCaloTriggerPatchInfoName("EmcalTriggers");
66 }
67 
69  AliAnalysisTaskEmcal(name, true),
70  fTriggerSelection(nullptr),
71  fUseTriggerBits(kTRUE),
72  fRequireBunchCrossing(kTRUE),
73  fHistos(nullptr),
74  fTriggerStringFromPatches(kFALSE),
75  fSelectedTriggers(),
76  fRequireAnalysisUtils(kTRUE),
77  fVertexCut(-10., 10.),
78  fNameDownscaleOADB(""),
79  fDownscaleOADB(nullptr),
80  fDownscaleFactors(nullptr),
81  fNameMaskedFastorOADB(),
82  fMaskedFastorOADB(nullptr),
83  fMaskedFastors(),
84  fOnlineTriggerThresholds(),
85  fNameAcceptanceOADB(),
86  fSelectNoiseEvents(false),
87  fRejectNoiseEvents(false),
88  fEnableDCALTriggers(true)
89 {
90  SetNeedEmcalGeom(true);
92  SetCaloTriggerPatchInfoName("EmcalTriggers");
93 }
94 
97  if(fHistos) delete fHistos;
98  if(fDownscaleOADB) delete fDownscaleOADB;
99 }
100 
103  if(fRequireAnalysisUtils && !fAliAnalysisUtils) fAliAnalysisUtils = new AliAnalysisUtils;
104 
105  fHistos = new THistManager(Form("Histos_%s", GetName()));
106 
109 
110  for(auto h : *(fHistos->GetListOfHistograms())) fOutput->Add(h);
111 
112  PostData(1, fOutput);
113 }
114 
116  // Apply trigger selection
118  if(!fSelectedTriggers.size()) return false;
119 
121 
122  const AliVVertex *vtx = fInputEvent->GetPrimaryVertex();
123  //if(!fInputEvent->IsPileupFromSPD(3, 0.8, 3., 2., 5.)) return; // reject pileup event
124  if(vtx->GetNContributors() < 1) return false;
125 
127  if(fInputEvent->IsA() == AliESDEvent::Class() && fAliAnalysisUtils->IsFirstEventInChunk(fInputEvent)) return false;
128  if(fAliAnalysisUtils->IsPileUpEvent(fInputEvent)) return false; // Apply new vertex cut
129  if(!fAliAnalysisUtils->IsVertexSelected2013pA(fInputEvent))return false; // Apply new vertex cut
130  }
131 
132  if(!fVertexCut.IsInRange(fVertex[2])) return false;
133  if(!IsUserEventSelected()) return false;
134 
135  // Do MC outlier cut
136  if(fIsPythia){
137  if(!CheckMCOutliers()){
138  AliDebugStream(1) << GetName() << ": Reject MC outliers" << std::endl;
139  return false;
140  }
141  }
142 
144  return true;
145 }
146 
147 
149  fSelectedTriggers.clear();
150  Bool_t isMC = MCEvent() != nullptr;
151 
152  TString triggerstring = "";
155  } else {
156  triggerstring = fInputEvent->GetFiredTriggerClasses();
157  }
158 
159  if(fRequireBunchCrossing && ! triggerstring.Contains("-B-")) return;
160 
161  UInt_t selectionstatus = fInputHandler->IsEventSelected();
162  Bool_t isMinBias = selectionstatus & AliVEvent::kINT7,
164  for(int itrg = 0; itrg < AliEmcalTriggerOfflineSelection::kTrgn; itrg++) emcalTriggers[itrg] = true;
165  if(!isMC){
166  // In case of data use information from the physics selection and event record
167  // Further cleanup of trigger events can be performed depending on the presence
168  // of recalc patches (after masking hot fastors in the trigger maker) above
169  // threshold
170  if(fUseTriggerBits){
171  const std::array<ULong_t, AliEmcalTriggerOfflineSelection::kTrgn> kSelectTriggerBits = {
172  AliVEvent::kEMC7, AliVEvent::kEMCEGA, AliVEvent::kEMCEGA, AliVEvent::kEMCEJE, AliVEvent::kEMCEJE,
173  AliVEvent::kEMC7, AliVEvent::kEMCEGA, AliVEvent::kEMCEGA, AliVEvent::kEMCEJE, AliVEvent::kEMCEJE,
174  };
175  for(int iclass = 0; iclass < AliEmcalTriggerOfflineSelection::kTrgn; iclass++)
176  emcalTriggers[iclass] &= bool(selectionstatus & kSelectTriggerBits[iclass]);
177  }
178 
179  // Apply cut on the trigger string - this basically discriminates high- and low-threshold
180  // triggers
181  const std::array<TString, AliEmcalTriggerOfflineSelection::kTrgn> kSelectTriggerStrings = {
182  "CEMC7", "EG1", "EG2", "EJ1", "EJ2", "CDMC7", "DG1", "DG2", "DJ1", "DJ2"
183  };
184  if(triggerstring.Contains("EMC")) AliDebugStream(1) << GetName() << ": Trigger string " << triggerstring << std::endl;
185  for(int iclass = 0; iclass < AliEmcalTriggerOfflineSelection::kTrgn; iclass++){
186  emcalTriggers[iclass] &= triggerstring.Contains(kSelectTriggerStrings[iclass]);
187  if(emcalTriggers[iclass])
188  AliDebugStream(1) << GetName() << ": Event selected as trigger " << kSelectTriggerStrings[iclass] << std::endl;
189  }
190 
191  // Online selection / rejection
193  if(fRejectNoiseEvents){
194  for(int itrg = 0; itrg < AliEmcalTriggerOfflineSelection::kTrgn; itrg++){
195  if(emcalTriggers[itrg])
197  }
198  } else {
199  for(int itrg = 0; itrg < AliEmcalTriggerOfflineSelection::kTrgn; itrg++){
200  if(emcalTriggers[itrg])
202  }
203  }
204  }
205  }
206  // Apply offline trigger selection: In this case cuts are performed on the
207  // patch energy from EMCAL cells after calibration. This method is most relevant
208  // for simulations. It can have a special use case in data in case a stronger
209  // offline selection is applied in addition to the online selection.
211  for(int itrg = 0; itrg < AliEmcalTriggerOfflineSelection::kTrgn; itrg++)
213  }
214  if(isMinBias) fSelectedTriggers.push_back("MB");
215  if(emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgEL0]){
216  fSelectedTriggers.push_back("EMC7");
217  if(!isMinBias) fSelectedTriggers.push_back("EMC7excl");
218  }
219  if(emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgEJ2]){
220  fSelectedTriggers.push_back("EJ2");
221  if(!(isMinBias || emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgEL0])) fSelectedTriggers.push_back("EJ2excl");
222  }
223  if(emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgEJ1]){
224  fSelectedTriggers.push_back("EJ1");
225  if(!(isMinBias || emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgEL0] || emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgEJ2])) fSelectedTriggers.push_back("EJ1excl");
226  }
227  if(emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgEG2]){
228  fSelectedTriggers.push_back("EG2");
229  if(!(isMinBias || emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgEL0])) fSelectedTriggers.push_back("EG2excl");
230  }
231  if(emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgEG1]){
232  fSelectedTriggers.push_back("EG1");
233  if(!(isMinBias || emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgEL0] || emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgEG2])) fSelectedTriggers.push_back("EG1excl");
234  }
235 
237  // Handle DCAL triggers only in case DCAL triggers are enabled,
238  // otherwise ignore results of the online/offline trigger selection
239  if(emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgDL0]){
240  fSelectedTriggers.push_back("DMC7");
241  if(!isMinBias) fSelectedTriggers.push_back("DMC7excl");
242  }
243  if(emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgDJ2]){
244  fSelectedTriggers.push_back("DJ2");
245  if(!(isMinBias || emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgDL0])) fSelectedTriggers.push_back("DJ2excl");
246  }
247  if(emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgDJ1]){
248  fSelectedTriggers.push_back("DJ1");
249  if(!(isMinBias || emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgDL0] || emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgDJ2])) fSelectedTriggers.push_back("DJ1excl");
250  }
251  if(emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgDG2]){
252  fSelectedTriggers.push_back("DG2");
253  if(!(isMinBias || emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgDL0])) fSelectedTriggers.push_back("DG2excl");
254  }
255  if(emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgDG1]){
256  fSelectedTriggers.push_back("DG1");
257  if(!(isMinBias || emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgDL0] || emcalTriggers[AliEmcalTriggerOfflineSelection::kTrgDG2])) fSelectedTriggers.push_back("DG1excl");
258  }
259  }
260 }
261 
264 
265  if(!fLocalInitialized){
266  return;
267  }
268 
269  // Handle OADB container with downscaling factors
270  if(fNameDownscaleOADB.Length()){
271  if(fNameDownscaleOADB.Contains("alien://") && ! gGrid) TGrid::Connect("alien://");
272  fDownscaleOADB = new AliOADBContainer("AliEmcalDownscaleFactors");
273  fDownscaleOADB->InitFromFile(fNameDownscaleOADB.Data(), "AliEmcalDownscaleFactors");
274  }
275  // Load OADB container with masked fastors (in case fastor masking is switched on)
277  if(fNameMaskedFastorOADB.Contains("alien://") && ! gGrid) TGrid::Connect("alien://");
278  fMaskedFastorOADB = new AliOADBContainer("AliEmcalMaskedFastors");
279  fMaskedFastorOADB->InitFromFile(fNameMaskedFastorOADB.Data(), "AliEmcalMaskedFastors");
280  }
281 
282  // Setting online threshold for trigger
284  if(fInputEvent->GetRunNumber() >= 15344 && fInputEvent->GetRunNumber() <= 197388){
289  SetOnlineTriggerThreshold("DG1", 0);
290  SetOnlineTriggerThreshold("DG2", 0);
291  SetOnlineTriggerThreshold("DJ1", 0);
292  SetOnlineTriggerThreshold("DJ2", 0);
293  SetOnlineTriggerThreshold("EMC7", 0);
294  SetOnlineTriggerThreshold("DMC7", 0);
295  }
296  }
297 
298  // Load acceptance OADB
299  if(fNameAcceptanceOADB.Length() && fTriggerSelection){
300  AliDebugStream(1) << GetName() << ": Loading acceptance map from OADB file " << fNameAcceptanceOADB << std::endl;
301  AliOADBContainer acceptanceCont("AliEmcalTriggerAcceptance");
302  acceptanceCont.InitFromFile(fNameAcceptanceOADB.Data(), "AliEmcalTriggerAcceptance");
303  TObjArray *acceptanceMaps = dynamic_cast<TObjArray *>(acceptanceCont.GetObject(fInputEvent->GetRunNumber()));
304  TH2 *map(nullptr);
305  if((map = dynamic_cast<TH2 *>(acceptanceMaps->FindObject("EG1")))){
306  AliDebugStream(1) << GetName() << ": Found acceptance map for trigger EG1" << std::endl;
307  map->SetDirectory(nullptr);
309  }
310  if((map = dynamic_cast<TH2 *>(acceptanceMaps->FindObject("EG2")))){
311  AliDebugStream(1) << GetName() << ": Found acceptance map for trigger EG2" << std::endl;
312  map->SetDirectory(nullptr);
314  }
315  if((map = dynamic_cast<TH2 *>(acceptanceMaps->FindObject("DG1")))){
316  AliDebugStream(1) << GetName() << ": Found acceptance map for trigger DG1" << std::endl;
317  map->SetDirectory(nullptr);
319  }
320  if((map = dynamic_cast<TH2 *>(acceptanceMaps->FindObject("DG2")))){
321  AliDebugStream(1) << GetName() << ": Found acceptance map for trigger DG2" << std::endl;
322  map->SetDirectory(nullptr);
324  }
325  if((map = dynamic_cast<TH2 *>(acceptanceMaps->FindObject("EJ1")))){
326  AliDebugStream(1) << GetName() << ": Found acceptance map for trigger EJ1" << std::endl;
327  map->SetDirectory(nullptr);
329  }
330  if((map = dynamic_cast<TH2 *>(acceptanceMaps->FindObject("EJ2")))){
331  AliDebugStream(1) << GetName() << ": Found acceptance map for trigger EJ2" << std::endl;
332  map->SetDirectory(nullptr);
334  }
335  if((map = dynamic_cast<TH2 *>(acceptanceMaps->FindObject("DJ1")))){
336  AliDebugStream(1) << GetName() << ": Found acceptance map for trigger DJ1" << std::endl;
337  map->SetDirectory(nullptr);
339  }
340  if((map = dynamic_cast<TH2 *>(acceptanceMaps->FindObject("DJ2")))){
341  AliDebugStream(1) << GetName() << ": Found acceptance map for trigger DJ2" << std::endl;
342  map->SetDirectory(nullptr);
344  }
345  }
346 }
347 
349  if(fDownscaleOADB){
350  fDownscaleFactors = static_cast<TObjArray *>(fDownscaleOADB->GetObject(runnumber));
351  }
352  if(fMaskedFastorOADB){
353  fMaskedFastors.clear();
354  TObjArray *ids = static_cast<TObjArray *>(fMaskedFastorOADB->GetObject(runnumber));
355  for(auto m : *ids){
356  TParameter<int> *id = static_cast<TParameter<int> *>(m);
357  fMaskedFastors.push_back(id->GetVal());
358  }
359  }
360 }
361 
363  // Exclusive means classes without lower trigger classes (which are downscaled) -
364  // in order to make samples statistically independent: MBExcl means MinBias && !EMCAL trigger
365  std::vector<TString> triggers = {
366  "MB", "EMC7", "EJ1", "EJ2", "EG1", "EG2", "EMC7excl", "EG2excl", "EJ2excl", "EJ1excl", "EG1excl"
367  };
368  std::vector<TString> dcaltriggers = {
369  "DMC7", "DJ1", "DJ2", "DG1", "DG2", "DMC7excl", "DG2excl", "DJ2excl", "DJ1excl", "DG1excl"
370  };
372  for(const auto &t : dcaltriggers) triggers.push_back(t);
373  }
374  return triggers;
375 }
376 
378  AliDebugStream(1) << "Using V1 online trigger selector" << std::endl;
380  int ngood(0);
381  for(auto p : *fTriggerPatchInfo){
382  AliEMCALTriggerPatchInfo *patch = static_cast<AliEMCALTriggerPatchInfo *>(p);
383  if((AliEmcalTriggerOfflineSelection::IsSingleShower(trigger) && patch->IsGammaLowRecalc()) ||
384  (!AliEmcalTriggerOfflineSelection::IsSingleShower(trigger) && patch->IsJetLowRecalc())){
385  if(patch->GetADCAmp() >= GetOnlineTriggerThresholdByIndex(trigger)) ngood++;
386  }
387  }
388  return ngood > 0;
389 }
390 
392  if(fDownscaleFactors){
393  TParameter<double> *result(nullptr);
394  // Downscaling only done on MB, L0 and the low threshold triggers
395  if(triggerclass.Contains("MB")) result = static_cast<TParameter<double> *>(fDownscaleFactors->FindObject("INT7"));
396  else if(triggerclass.Contains("EMC7")) result = static_cast<TParameter<double> *>(fDownscaleFactors->FindObject("EMC7"));
397  else if(triggerclass.Contains("EJ2")) result = static_cast<TParameter<double> *>(fDownscaleFactors->FindObject("EJ2"));
398  else if(triggerclass.Contains("EG2")) result = static_cast<TParameter<double> *>(fDownscaleFactors->FindObject("EG2"));
399  if(result) return 1./result->GetVal();
400  }
401  return 1.;
402 }
403 
405  TParameter<int> *threshpar(nullptr);
406  if((threshpar = static_cast<TParameter<int> *>(fOnlineTriggerThresholds.FindObject(triggerclass.Data())))){
407  threshpar->SetVal(threshold);
408  } else {
409  fOnlineTriggerThresholds.Add(new TParameter<int>(triggerclass.Data(), threshold));
410  }
411 }
412 
414  Int_t threshold(0);
415  TParameter<int> *val(nullptr);
416  if((val = static_cast<TParameter<int> *>(fOnlineTriggerThresholds.FindObject(name))))
417  threshold = val->GetVal();
418  return threshold;
419 }
420 
422  const TString triggernames[AliEmcalTriggerOfflineSelection::kTrgn] = {"EMC7", "EG1", "EG2", "EJ1", "EJ2", "DMC7",
423  "DG1", "DG2", "DJ1", "DJ2"};
424  return GetOnlineTriggerThresholdByName(triggernames[trigger]);
425 }
426 
428  const TString triggernames[AliEmcalTriggerOfflineSelection::kTrgn] = {"EMC7", "EG1", "EG2", "EJ1", "EJ2", "DMC7",
429  "DG1", "DG2", "DJ1", "DJ2"};
430  bool isInitialized = true;
431  for(int itrg = 0; itrg < AliEmcalTriggerOfflineSelection::kTrgn; itrg++){
432  if(!fOnlineTriggerThresholds.FindObject(triggernames[itrg].Data())) {
433  isInitialized = false;
434  break;
435  }
436  }
437  return isInitialized;
438 }
439 
441  return adc > GetOnlineTriggerThresholdByName(triggerclass);
442 }
443 
445  TString triggerstring = "";
446  Int_t nEJ1 = 0, nEJ2 = 0, nEG1 = 0, nEG2 = 0, nDJ1 = 0, nDJ2 = 0, nDG1 = 0, nDG2 = 0;
447  for(auto patchIter : *triggerpatches){
448  AliEMCALTriggerPatchInfo *patch = dynamic_cast<AliEMCALTriggerPatchInfo *>(patchIter);
449  if(!patch->IsRecalc()) continue;
450  if(patch->IsEMCal()){
451  if(patch->IsGammaLowRecalc() && SelectFiredPatch("EG1", patch->GetADCAmp())) nEG1++;
452  if(patch->IsGammaLowRecalc() && SelectFiredPatch("EG2", patch->GetADCAmp())) nEG2++;
453  if(patch->IsJetLowRecalc() && SelectFiredPatch("EJ1", patch->GetADCAmp())) nEJ1++;
454  if(patch->IsJetLowRecalc() && SelectFiredPatch("EJ2", patch->GetADCAmp())) nEJ2++;
455  } else {
456  if(patch->IsGammaLowRecalc() && SelectFiredPatch("DG1", patch->GetADCAmp())) nDG1++;
457  if(patch->IsGammaLowRecalc() && SelectFiredPatch("DG2", patch->GetADCAmp())) nDG2++;
458  if(patch->IsJetLowRecalc() && SelectFiredPatch("DJ1", patch->GetADCAmp())) nDJ1++;
459  if(patch->IsJetLowRecalc() && SelectFiredPatch("DJ2", patch->GetADCAmp())) nDJ2++;
460  }
461  }
462  if(nEJ1) triggerstring += "EJ1";
463  if(nEJ2){
464  if(triggerstring.Length()) triggerstring += ",";
465  triggerstring += "EJ2";
466  }
467  if(nEG1){
468  if(triggerstring.Length()) triggerstring += ",";
469  triggerstring += "EG1";
470  }
471  if(nEG2){
472  if(triggerstring.Length()) triggerstring += ",";
473  triggerstring += "EG2";
474  }
475  if(nDJ1){
476  if(triggerstring.Length()) triggerstring += ",";
477  triggerstring += "DJ2";
478  }
479  if(nDJ2){
480  if(triggerstring.Length()) triggerstring += ",";
481  triggerstring += "DJ2";
482  }
483  if(nDG1){
484  if(triggerstring.Length()) triggerstring += ",";
485  triggerstring += "DG1";
486  }
487  if(nDG2){
488  if(triggerstring.Length()) triggerstring += ",";
489  triggerstring += "DG2";
490  }
491  return triggerstring;
492 }
493 
494 } /* namespace EMCalTriggerPtAnalysis */
std::vector< TString > fSelectedTriggers
! Triggers selected for given event
Bool_t fIsPythia
trigger, if it is a PYTHIA production
double Double_t
Definition: External.C:58
Int_t GetOnlineTriggerThresholdByIndex(AliEmcalTriggerOfflineSelection::EmcalTriggerClass trg) const
Bool_t fRejectNoiseEvents
Reject events triggered by noisy fastors.
Base task in the EMCAL framework.
void SetOnlineTriggerThreshold(const TString &triggerclass, Int_t threshold)
Bool_t fLocalInitialized
whether or not the task has been already initialized
TString GetFiredTriggerClassesFromPatches(const TClonesArray *triggerpatches) const
AliEmcalTriggerOfflineSelection * fTriggerSelection
Offline trigger selection.
TObjArray fOnlineTriggerThresholds
Trigger thresholds applied at online level.
Bool_t fEnableDCALTriggers
Enable / Disable event selection for DCAL trigger classes.
void SetCaloTriggerPatchInfoName(const char *n)
TString fNameMaskedFastorOADB
Name of the masked fastor OADB container.
Bool_t fRequireAnalysisUtils
Switch whether to require event selection in AliAnalysisUtils.
bool SelectOnlineTrigger(AliEmcalTriggerOfflineSelection::EmcalTriggerClass trigger) const
Bool_t fRequireBunchCrossing
Require bunch-bunch events (tag -B- in trigger string)
void SetAcceptanceMap(EmcalTriggerClass trgcls, const TH2 *accmap)
int Int_t
Definition: External.C:63
unsigned int UInt_t
Definition: External.C:33
THashList * GetListOfHistograms() const
Definition: THistManager.h:504
EmcalTriggerClass
Definition of the various supported trigger types.
AliAnalysisUtils * fAliAnalysisUtils
!vertex selection (optional)
TString fNameDownscaleOADB
Name of the downscale OADB container.
TString fNameAcceptanceOADB
Name of the OADB container with the trigger acceptance.
AliOADBContainer * fMaskedFastorOADB
! Container with masked fastors
Bool_t fSelectNoiseEvents
Explicitly select events triggered only by noisy fastors.
Bool_t IsOfflineSelected(EmcalTriggerClass trgcls, const TClonesArray *const triggerpatches) const
Bool_t fUseTriggerBits
Switch whether using trigger bits (relies on physics selection)
Bool_t fTriggerStringFromPatches
Do rebuild the trigger string from trigger patches.
Bool_t isMC
AliEmcalList * fOutput
!output list
Definition: External.C:220
AliOADBContainer * fDownscaleOADB
! Container with downscale factors for different triggers
ClassImp(AliAnalysisTaskCRC) AliAnalysisTaskCRC
Double_t fVertex[3]
!event vertex
void SetMakeGeneralHistograms(Bool_t g)
TClonesArray * fTriggerPatchInfo
!trigger patch info array
void SetNeedEmcalGeom(Bool_t n)
Container class for histograms for the high- charged particle analysis.
Definition: THistManager.h:43
AliCutValueRange< double > fVertexCut
Cut on the z-position of the primary vertex.
bool Bool_t
Definition: External.C:53
Bool_t SelectFiredPatch(const TString &triggerclass, Int_t adc) const