AliPhysics  aaf9c62 (aaf9c62)
AddTaskGammaHadronCorrelation.C
Go to the documentation of this file.
1 
19 // Set includes for compilation
20 
21 #if !defined(__CINT__) || defined(__MAKECINT__)
22 
23 #include <TString.h>
24 #include <TSystem.h>
25 #include <TROOT.h>
26 
27 #include "AliLog.h"
29 #include "AliCaloTrackESDReader.h"
30 #include "AliCaloTrackAODReader.h"
31 #include "AliCalorimeterUtils.h"
32 #include "AliAnaPhoton.h"
33 #include "AliAnaPi0EbE.h"
34 #include "AliHistogramRanges.h"
37 #include "AliAnaChargedParticles.h"
38 #include "AliAnaCalorimeterQA.h"
39 #include "AliAnaGeneratorKine.h"
42 #include "AliAnalysisManager.h"
43 #include "AliInputEventHandler.h"
44 #include "AliVTrack.h"
46 #include "AliESDtrackCuts.h"
47 #include "CreateTrackCutsPWGJE.C"
48 //#include "ConfigureEMCALRecoUtils.C"
50 #endif
51 
52 // Declare methods for compilation
53 
55  TString clustersArray, Bool_t tender,
56  TString calorimeter, Bool_t nonLinOn,
57  TString trigger, Bool_t rejectEMCTrig,
58  Int_t minCen, Int_t maxCen,
59  Bool_t printSettings, Int_t debug );
61  Bool_t tender, TString calorimeter,
62  Bool_t nonLinOn, Int_t year,
63  Bool_t printSettings, Int_t debug );
65  TString calorimeter, Int_t year, Int_t tm,
66  Bool_t printSettings, Int_t debug );
68  Bool_t useSSIso, Bool_t useAsy,
69  TString col, Bool_t simulation,
70  TString calorimeter, Int_t year, Int_t tm,
71  Bool_t printSettings, Int_t debug );
73  (TString particle, Int_t leading,
74  Int_t partInCone, Int_t thresType,
75  Float_t cone, Float_t pth, Bool_t multi,
76  TString col, Bool_t simulation,
77  TString calorimeter, Int_t year, Int_t tm,
78  Bool_t printSettings, Int_t debug );
80  (TString particle, Int_t leading,
81  Bool_t bIsolated, Float_t shshMax,
82  Int_t partInCone, Int_t thresType,
83  Float_t cone, Float_t pth, Bool_t mixOn,
84  TString col, Bool_t simulation,
85  TString calorimeter, Int_t year, Int_t tm,
86  Bool_t printSettings, Int_t debug );
88  (Bool_t simulation, Bool_t printSettings, Int_t debug);
91  Bool_t printSettings, Int_t debug );
93  TString col, Bool_t simulation,
95  Bool_t printSettings, Int_t debug );
96 
98  TString calorimeter , Int_t year,
99  TString col , Bool_t simulation,
100  Bool_t printSettings, Int_t debug);
101 
102 
105 
137 (
138  TString calorimeter = "EMCAL", // "DCAL", "PHOS"
139  Bool_t simulation = kFALSE,
140  Int_t year = 2011,
141  TString col = "pp",
142  Int_t rejectEMCTrig = 0,
143  TString clustersArray = "",
144  Bool_t tender = kFALSE,
145  Bool_t nonLinOn = kFALSE,
146  Float_t shshMax = 0.27,
147  Float_t isoCone = 0.4,
148  Float_t isoPtTh = 0.5,
151  Int_t leading = 0,
152  Int_t tm = 1,
153  Int_t minCen = -1,
154  Int_t maxCen = -1,
155  Bool_t mixOn = kTRUE,
156  Bool_t qaAn = kFALSE,
157  Bool_t chargedAn = kFALSE,
158  TString outputfile = "",
159  Bool_t printSettings = kFALSE,
160  Int_t debug = 0,
161  const char *trigSuffix = "EMC7"
162 )
163 {
164  // Check the global variables, and reset the provided ones if empty.
165  //
166  TString trigger = trigSuffix;
167 
168  gROOT->LoadMacro("$ALICE_PHYSICS/PWGGA/CaloTrackCorrelations/macros/GetAlienGlobalProductionVariables.C");
169 
170  TString period = "LHC17";
171  Bool_t printGlobal = kTRUE;
172 
173  GetAlienGlobalProductionVariables(simulation,col,period,year,printGlobal);
174 
175  printf("Passed settings:\n calorimeter <%s>, simulation <%d>, year <%d>,\n col <%s>, trigger <%s>, reject EMC <%d>, clustersArray <%s>, tender <%d>, non linearity <%d>\n shshMax <%2.2f>, isoCone <%2.2f>, isoPtTh <%2.2f>, isoMethod <%d>,isoContent <%d>,\n leading <%d>, tm <%d>, minCen <%d>, maxCen <%d>, mixOn <%d>,\n qaAn <%d>, chargedAn <%d>, outputfile <%s>, printSettings <%d>, debug <%d>\n",
176  calorimeter.Data(),simulation,year,col.Data(),trigger.Data(), rejectEMCTrig, clustersArray.Data(),tender, nonLinOn, shshMax,
177  isoCone,isoPtTh,isoMethod,isoContent,leading,tm,
178  minCen,maxCen,mixOn,qaAn,chargedAn,outputfile.Data(),printSettings,debug);
179 
180  // Get the pointer to the existing analysis manager via the static access method.
181  AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
182  if (!mgr)
183  {
184  ::Error("AddTask", "No analysis manager to connect to.");
185  return NULL;
186  }
187 
188  // Check the analysis type using the event handlers connected to the analysis manager.
189 
190  if (!mgr->GetInputEventHandler())
191  {
192  ::Error("AddTask", "This task requires an input event handler");
193  return NULL;
194  }
195 
196  // Make sure the B field is enabled for track selection, some cuts need it
197  ((AliInputEventHandler*)mgr->GetInputEventHandler())->SetNeedField(kTRUE);
198 
199  // Name for containers
200 
201  //kAnaGammaHadronCorr = Form("%s_Trig%s_Cl%s_TM%d_R%1.1f_Pt%1.1f",calorimeter.Data(), trigger.Data(),kClusterArray.Data(),tm,cone,pth);
202  kAnaGammaHadronCorr = Form("GammaHadron_%s_Trig%s_Col_%s_Year%d_Cl%s_Ten%d_TM%d_M02_%1.2f_IsoParam_C%1.2fPt%1.2fM%dPa%d_Lead%d_Mix%d",
203  calorimeter.Data(),trigger.Data(),col.Data(),year,clustersArray.Data(),tender,
204  tm, shshMax,isoCone,isoPtTh,isoMethod,isoContent,leading,mixOn);
205 
206  if(col=="PbPb" && maxCen>=0) kAnaGammaHadronCorr+=Form("Cen%d_%d",minCen,maxCen);
207 
208  printf("<<<< NAME: %s >>>>>\n",kAnaGammaHadronCorr.Data());
209 
210  // #### Configure analysis ####
211 
213 
214  // General frame setting and configuration
215  maker->SetReader ( ConfigureReader (col,simulation,clustersArray,tender,calorimeter,nonLinOn,trigger,rejectEMCTrig,minCen,maxCen,printSettings,debug) );
216  maker->SetCaloUtils( ConfigureCaloUtils(col,simulation, tender,calorimeter,nonLinOn,year,printSettings,debug) );
217 
218  // Analysis tasks setting and configuration
219  Int_t n = 0;//Analysis number, order is important
220 
221  //
222  // Photon analysis
223  //
224  maker->AddAnalysis(ConfigurePhotonAnalysis(col,simulation,calorimeter,year,tm,printSettings,debug), n++); // Photon cluster selection
225 
226 // maker->AddAnalysis(ConfigurePi0EbEAnalysis("Pi0" , AliAnaPi0EbE::kIMCalo,kFALSE,kFALSE,
227 // col,simulation,calorimeter,year,tm,printSettings,debug), n++); // Pi0 event by event selection, invariant mass and photon tagging from decay
228 //
229 // maker->AddAnalysis(ConfigurePi0EbEAnalysis("Eta" , AliAnaPi0EbE::kIMCalo,kFALSE,kFALSE,
230 // col,simulation,calorimeter,year,tm,printSettings,debug), n++); // Eta event by event selection, invariant mass and photon tagging from decay
231 //
232 // maker->AddAnalysis(ConfigurePi0EbEAnalysis("Pi0SideBand", AliAnaPi0EbE::kIMCalo,kFALSE,kFALSE,
233 // col,simulation,calorimeter,year,tm,printSettings,debug), n++); // Pi0 out of peak event by event selection, and photon tagging from decay
234 //
235 // maker->AddAnalysis(ConfigurePi0EbEAnalysis("EtaSideBand", AliAnaPi0EbE::kIMCalo,kFALSE,kFALSE,
236 // col,simulation,calorimeter,year,tm,printSettings,debug), n++); // Eta out of peak event by event selection, and photon tagging from decay
237 
238 
239  maker->AddAnalysis(ConfigureIsolationAnalysis("Photon", leading, isoContent,isoMethod,isoCone,isoPtTh, kFALSE,
240  col,simulation,calorimeter,year,tm,printSettings,debug), n++); // Photon isolation
241 
242 
243  maker->AddAnalysis(ConfigureHadronCorrelationAnalysis("Photon", leading, kFALSE, shshMax, isoContent,isoMethod,isoCone,isoPtTh, mixOn,
244  col,simulation,calorimeter,year,tm,printSettings,debug), n++); // Gamma-hadron correlation
245 
246  maker->AddAnalysis(ConfigureHadronCorrelationAnalysis("Photon", leading, kTRUE, shshMax, isoContent,isoMethod,isoCone,isoPtTh, mixOn,
247  col,simulation,calorimeter,year,tm,printSettings,debug) , n++); // Isolated gamma hadron correlation
248 
249  //
250  // Merged pi0 analysis
251  //
252  maker->AddAnalysis(ConfigurePi0EbEAnalysis("Pi0", AliAnaPi0EbE::kSSCalo,kTRUE,kTRUE,
253  col,simulation,calorimeter,year,tm,printSettings,debug), n++); // Pi0 event by event selection, cluster splitting
254 
255  maker->AddAnalysis(ConfigureIsolationAnalysis("Pi0SS", leading, isoContent,isoMethod,isoCone,isoPtTh, kFALSE,
256  col,simulation,calorimeter,year,tm,printSettings,debug), n++); // Pi0 isolation, cluster splits
257 
258 
259  maker->AddAnalysis(ConfigureHadronCorrelationAnalysis("Pi0SS", leading, kFALSE, shshMax, isoContent,isoMethod,isoCone,isoPtTh, mixOn,
260  col,simulation,calorimeter,year,tm,printSettings,debug), n++); // Pi0-hadron correlation
261 
262  maker->AddAnalysis(ConfigureHadronCorrelationAnalysis("Pi0SS", leading, kTRUE, shshMax, isoContent,isoMethod,isoCone,isoPtTh, mixOn,
263  col,simulation,calorimeter,year,tm,printSettings,debug) , n++); // Isolated pi0-hadron correlation
264 
265  // Check the generated kinematics
266 // if(simulation) maker->AddAnalysis(ConfigureGenKineAnalysis(isoMethod,isoCone,isoPtTh,
267 // col,simulation,calorimeter,year,printSettings,debug), n++);
268 
269  // Charged analysis
270  if(chargedAn) maker->AddAnalysis(ConfigureChargedAnalysis(simulation,printSettings,debug), n++); // track selection checks
271 
272  // Calo QA
273  if(qaAn) maker->AddAnalysis(ConfigureQAAnalysis(col,simulation,calorimeter,year,printSettings,debug) , n++);
274 
275  maker->SetAnaDebug(debug) ;
276 
277  maker->SwitchOnHistogramsMaker() ;
278  maker->SwitchOnAODsMaker() ;
279 
280  if( simulation || !trigger.Contains("EMC") ) maker->SwitchOffDataControlHistograms();
281 
282  if(simulation)
283  {
284  // Calculate the cross section weights, apply them to all histograms
285  // and fill xsec and trial histo. Sumw2 must be activated.
286  //maker->GetReader()->GetWeightUtils()->SwitchOnMCCrossSectionCalculation();
287  //maker->SwitchOnSumw2Histograms();
288 
289  // For recent productions where the cross sections and trials are not stored in separate file
290  //maker->GetReader()->GetWeightUtils()->SwitchOnMCCrossSectionFromEventHeader() ;
291 
292  // Just fill cross section and trials histograms.
293  maker->GetReader()->GetWeightUtils()->SwitchOnMCCrossSectionHistoFill();
294 
295  // Add control histogram with pT hard to control aplication of weights
296  maker->SwitchOnPtHardHistogram();
297  }
298 
299  if(printSettings) maker->Print("");
300 
301  printf("<< End Configuration of %d analysis for calorimeter %s >>\n",n, calorimeter.Data());
302 
303  //
304  // Create task, pass the maker and add it to the manager
305  //
306  AliAnalysisTaskCaloTrackCorrelation * task = new AliAnalysisTaskCaloTrackCorrelation (Form("%s",kAnaGammaHadronCorr.Data()));
307 
308  task->SetDebugLevel(debug);
309 
310  //task->SetBranches("ESD:AliESDRun.,AliESDHeader");
311  //task->SetConfigFileName(""); //Don't configure the analysis via configuration file.
312 
313  task->SetAnalysisMaker(maker);
314 
315  mgr->AddTask(task);
316 
317  //
318  // Select events trigger depending on trigger
319  //
320  maker->GetReader()->SwitchOnEventTriggerAtSE(); // on is default case
321  if(!simulation)
322  {
323  gROOT->LoadMacro("$ALICE_PHYSICS/PWGGA/CaloTrackCorrelations/macros/ConfigureAndGetEventTriggerMaskAndCaloTriggerString.C");
324  TString caloTriggerString = "";
325  UInt_t mask = ConfigureAndGetEventTriggerMaskAndCaloTriggerString(trigger, year, caloTriggerString);
326 
327  maker->GetReader()->SetFiredTriggerClassName(caloTriggerString);
328 
329  // For mixing with AliAnaParticleHadronCorrelation switch it off
330  if(mixOn)
331  {
332  maker->GetReader()->SwitchOffEventTriggerAtSE();
333  maker->GetReader()->SetEventTriggerMask(mask);
334  // what to do with caloTriggerString?
335 
336  // Careful, not all productions work with kMB, try kINT7, kINT1, kAnyINT
337  //reader->SetMixEventTriggerMask(AliVEvent::kMB);
338  maker->GetReader()->SetMixEventTriggerMask(AliVEvent::kINT7);
339 
340  printf("---Trigger selection done in AliCaloTrackReader!!!\n");
341  }
342  else
343  {
344  task ->SelectCollisionCandidates( mask );
345  }
346  }
347 
348  //
349  // Create containers
350  //
351  if(outputfile.Length()==0) outputfile = AliAnalysisManager::GetCommonFileName();
352 
353  AliAnalysisDataContainer *cout_pc = mgr->CreateContainer(kAnaGammaHadronCorr, TList::Class(),
354  AliAnalysisManager::kOutputContainer,
355  Form("%s",outputfile.Data()));
356 
357  AliAnalysisDataContainer *cout_cuts = mgr->CreateContainer(Form("Param_%s",kAnaGammaHadronCorr.Data()), TList::Class(),
358  AliAnalysisManager::kParamContainer,
359  "AnalysisParameters.root");
360 
361  // Create ONLY the output containers for the data produced by the task.
362  // Get and connect other common input/output containers via the manager as below
363  //==============================================================================
364  mgr->ConnectInput (task, 0, mgr->GetCommonInputContainer());
365  //if(!kData.Contains("delta") && outputAOD) mgr->ConnectOutput (task, 0, mgr->GetCommonOutputContainer());
366  mgr->ConnectOutput (task, 1, cout_pc);
367  mgr->ConnectOutput (task, 2, cout_cuts);
368 
369  return task;
370 }
371 
376  TString clustersArray, Bool_t tender,
377  TString calorimeter, Bool_t nonLinOn,
378  TString trigger, Bool_t rejectEMCTrig,
379  Int_t minCen, Int_t maxCen,
380  Bool_t printSettings, Int_t debug )
381 {
382  // Get the data type ESD or AOD
383  AliAnalysisManager * mgr = AliAnalysisManager::GetAnalysisManager();
384  TString inputDataType = mgr->GetInputEventHandler()->GetDataType();
385 
386  AliCaloTrackReader * reader = 0;
387  if (inputDataType == "AOD") reader = new AliCaloTrackAODReader();
388  else if(inputDataType == "ESD") reader = new AliCaloTrackESDReader();
389  else printf("AliCaloTrackReader::ConfigureReader() - Data not known InputData=%s\n",inputDataType.Data());
390 
391  reader->SetDebug(debug);//10 for lots of messages
392 
393  //
394  // MC settings
395  //
396  // Check if kine stack is available, independent of request of simulation
397 // Bool_t useKinematics = kFALSE;
398 // useKinematics = (mgr->GetMCtruthEventHandler())?kTRUE:kFALSE;
399 //
400 // if(simulation)
401 // {
402 // if (!useKinematics && inputDataType=="AOD") useKinematics = kTRUE; //AOD primary should be available ...
403 // }
404 
405  // In case of Pythia pt Hard bin simulations (jet-jet, gamma-jet)
406  // reject some special events that bother the cross section
407  if(simulation)
408  {
409  // Event rejection cuts for jet-jet simulations, do not use in other
410  reader->SetPtHardAndJetPtComparison(kTRUE);
411  reader->SetPtHardAndJetPtFactor(2);
412 
413  // Event rejection more suitable for gamma-jet simulations, do not use in other
414 // reader->SetPtHardAndClusterPtComparison(kTRUE);
415 // reader->SetPtHardAndClusterPtFactor(1.5);
416  }
417 
418  //---------------------------
419  // Detectors acceptance, open
420  //---------------------------
421 
422  reader->SwitchOnFiducialCut();
423 
424  // Tracks
425  reader->GetFiducialCut()->SetSimpleCTSFiducialCut(0.8, 0, 360) ;
426 
427  // EMCal/DCal
428  if ( calorimeter == "EMCAL" ) reader->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.70, 80, 187) ;
429  else if ( calorimeter == "DCAL" ) reader->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.70, 260, 327) ;
430 
431  //PHOS
432  reader->GetFiducialCut()->SetSimplePHOSFiducialCut (0.12, 250, 320) ;
433 
434  //------------------------
435  // Detector input filling
436  //------------------------
437 
438  //Min cluster/track E
439  reader->SetEMCALEMin(0.3);
440  reader->SetEMCALEMax(1000);
441  reader->SetPHOSEMin(0.3);
442  reader->SetPHOSEMax(1000);
443  reader->SetCTSPtMin(0.2);
444  reader->SetCTSPtMax(1000);
445 
448 
449  // Shower shape smearing
450  // Set it in the train configuration page not here for the moment
451  // if(simulation)
452  // {
453  // reader->SwitchOffShowerShapeSmearing(); // Active only on MC, off by default
454  // reader->SetShowerShapeSmearWidth(0.005);
455  // }
456 
457  //
458  // Tracks
459  //
460  reader->SwitchOnCTS();
461 
462  reader->SwitchOffUseTrackTimeCut();
463  reader->SetTrackTimeCut(0,50);
464 
465  reader->SwitchOffUseTrackDCACut();
466  //reader->SetTrackDCACut(0,0.0105);
467  //reader->SetTrackDCACut(1,0.035);
468  //reader->SetTrackDCACut(2,1.1);
469 
470  if(inputDataType=="ESD")
471  {
472  gROOT->LoadMacro("$ALICE_PHYSICS/PWGJE/macros/CreateTrackCutsPWGJE.C");
473 
474  //AliESDtrackCuts * esdTrackCuts = CreateTrackCutsPWGJE(10041004);
475  //reader->SetTrackCuts(esdTrackCuts);
476 
477  AliESDtrackCuts * esdTrackCuts = CreateTrackCutsPWGJE(10001008);
478  reader->SetTrackCuts(esdTrackCuts);
479  AliESDtrackCuts * esdTrackCuts2 = CreateTrackCutsPWGJE(10011008);
480  reader->SetTrackComplementaryCuts(esdTrackCuts2);
481 
483  }
484  else if(inputDataType=="AOD")
485  {
486  reader->SwitchOnAODHybridTrackSelection(); // Check that the AODs have Hybrids!!!!
488  reader->SetTrackStatus(AliVTrack::kITSrefit);
489 
490  //reader->SwitchOnAODPrimaryTrackSelection(); // Used in preliminary results of QM from Nicolas and Xiangrong?
491  //reader->SwitchOnTrackHitSPDSelection(); // Check that the track has at least a hit on the SPD, not much sense to use for hybrid or TPC only tracks
492  //reader->SetTrackFilterMask(128); // Filter bit, not mask, use if off hybrid, TPC only
493  }
494 
495  //
496  // Calorimeter
497  //
498  if(clustersArray == "" && !tender)
499  {
500  printf("**************** Standard EMCAL clusters branch analysis **************** \n");
502  // Check in ConfigureCaloUtils that the recalibration and bad map are ON
503  }
504  else
505  {
506  printf("**************** Input for analysis is Clusterizer %s **************** \n", clustersArray.Data());
507  reader->SetEMCALClusterListName(clustersArray);
509  }
510 
511  // Time cuts
513  if(simulation)
514  {
515  reader->SwitchOffUseEMCALTimeCut();
516  reader->SetEMCALTimeCut(-1e10,1e10); // Open time cut
517  }
518  else
519  {
520  reader->SwitchOnUseEMCALTimeCut();
521  reader->SetEMCALTimeCut(-25,20);
522  }
523 
524  // CAREFUL
525  if(nonLinOn) reader->SwitchOnClusterELinearityCorrection();
527 
528  if(calorimeter == "EMCAL" || calorimeter == "DCAL" )
529  {
530  reader->SwitchOnEMCALCells();
531  reader->SwitchOnEMCAL();
532  }
533 
534  if(calorimeter == "PHOS")
535  { // Should be on if QA is activated with correlation on
536  reader->SwitchOnPHOSCells();
537  reader->SwitchOnPHOS();
538  }
539 
540  //-----------------
541  // Event selection
542  //-----------------
543 
544  //if(!simulation) reader->SetFiredTriggerClassName("CEMC7EGA-B-NOPF-CENTNOTRD"); // L1 Gamma
545 
546  // Event triggered by EMCal selection settings
549 
550  if( rejectEMCTrig > 0 && !simulation && (trigger.Contains("EMC") || trigger.Contains("L")))
551  {
552  printf("=== Remove bad triggers === \n");
555 
556  // reader->SetTriggerPatchTimeWindow(8,9); // default values
557  // if (kRunNumber < 146861) reader->SetEventTriggerL0Threshold(3.);
558  // else if(kRunNumber < 154000) reader->SetEventTriggerL0Threshold(4.);
559  // else if(kRunNumber < 165000) reader->SetEventTriggerL0Threshold(5.5);
560  // //redefine for other periods, triggers
561  //
562  // if(kRunNumber < 172000)
563  // {
564  // reader->SetEventTriggerL1Bit(4,5); // current LHC11 data
565  // printf("\t Old L1 Trigger data format!\n");
566  // }
567  // else
568  // {
569  // reader->SetEventTriggerL1Bit(6,8); // LHC12-13 data
570  // printf("\t Current L1 Trigger data format!\n");
571  // }
572 
573  if(clustersArray != "" || tender)
574  {
575  printf("Trigger cluster calibration OFF\n");
577  }
578 
579  }
580 
581  //reader->RejectFastClusterEvents() ;
582 
583  reader->SetZvertexCut(10.); // Open cut
584  reader->SwitchOnPrimaryVertexSelection(); // and besides primary vertex
585  reader->SwitchOnRejectNoTrackEvents();
586 
587  reader->SwitchOffV0ANDSelection() ; // and besides v0 AND
588  reader->SwitchOffPileUpEventRejection(); // remove pileup by default off, apply it only for MB not for trigger
589 
590  if(col=="PbPb")
591  {
592  // Centrality
593  reader->SetCentralityClass("V0M");
594  reader->SetCentralityOpt(100); // 10 (c= 0-10, 10-20 ...), 20 (c= 0-5, 5-10 ...) or 100 (c= 1, 2, 3 ..)
595  reader->SetCentralityBin(minCen,maxCen); // Accept all events, if not select range
596  //reader->SwitchOnAcceptOnlyHIJINGLabels();
597  // Event plane (only used in Maker and mixing for AliAnaPi0/AliAnaHadronCorrelation for the moment)
598  reader->SetEventPlaneMethod("V0");
599  }
600 
601  if(printSettings) reader->Print("");
602 
603  return reader;
604 }
605 
610  Bool_t tender, TString calorimeter,
611  Bool_t nonLinOn,
612  Int_t year, Bool_t printSettings,
613  Int_t debug)
614 {
616 
617  cu->SetDebug(debug);
618 
619  // Remove clusters close to borders, at least max energy cell is 1 cell away
622 
624 
625  if ( calorimeter == "PHOS" )
626  {
627  if(year < 2014) cu->SetNumberOfSuperModulesUsed(3);
628  else cu->SetNumberOfSuperModulesUsed(4);
629  }
630  else
631  {
632  Int_t nSM = 20;
633  Int_t lastEMC = 11;
634  if (year == 2010) { nSM = 4; lastEMC = 3; }// EMCAL first year
635  else if (year < 2014) { nSM = 10; lastEMC = 9; }// EMCAL active 2011-2013
636 
638 
639  if (calorimeter.Contains("EMCAL"))
640  {
641  cu->SetFirstSuperModuleUsed( 0);
642  cu->SetLastSuperModuleUsed (lastEMC);
643  }
644  else if (calorimeter.Contains("DCAL"))
645  {
646  cu->SetFirstSuperModuleUsed(12);
647  cu->SetLastSuperModuleUsed (19);
648  }
649  else
650  {
653  }
654 
655  printf("AddTaskGammaHadronCorrelation::CalorimeterUtils() - nSM %d, first %d, last %d\n",
657  }
658 
659  // Search of local maxima in cluster
660  if(col=="pp")
661  {
662  cu->SetLocalMaximaCutE(0.1);
663  cu->SetLocalMaximaCutEDiff(0.03);
664  }
665  else
666  {
667  cu->SetLocalMaximaCutE(0.2);
668  cu->SetLocalMaximaCutEDiff(0.03);
669  }
670 
672 
674 
675  // EMCAL settings
676 
677  if(!simulation)
679 
680  AliEMCALRecoUtils * recou = cu->GetEMCALRecoUtils();
681 
682  // calibrations
683  Bool_t calibEner = kFALSE;
684  Bool_t calibTime = kFALSE;
687 
688  if( !tender )
689  {
690  cu->SwitchOnRecalibration(); // Check the reader if it is taken into account during filtering
692 
693  calibEner = kTRUE;
694  calibTime = kTRUE;
695  }
696 
697  if( simulation )
698  {
699  calibEner = kFALSE;
700  calibTime = kFALSE;
701 
702  cu->SwitchOffRecalibration(); // Check the reader if it is taken into account during filtering
704  }
705 
706 // gROOT->LoadMacro("$ALICE_PHYSICS/PWGPP/EMCAL/macros/ConfigureEMCALRecoUtils.C");
707 // ConfigureEMCALRecoUtils(recou,
708 // simulation,
709 // kTRUE, // exotic
710 // nonLinOn, // Non linearity
711 // calibEner, // E calib
712 // kTRUE, // bad map
713 // calibTime); // time calib
714 
715  cu->ConfigureEMCALRecoUtils(simulation,
716  kTRUE, // exotic
717  nonLinOn, // Non linearity
718  calibEner, // E calib
719  kTRUE, // bad map
720  calibTime); // time calib
721 
722  if( calibTime ) recou->SetExoticCellDiffTimeCut(50);
723 
724  if( nonLinOn ) cu->SwitchOnCorrectClusterLinearity();
725 
726  printf("ConfigureCaloUtils() - EMCAL Recalibration ON? %d %d\n",recou->IsRecalibrationOn(), cu->IsRecalibrationOn());
727  printf("ConfigureCaloUtils() - EMCAL BadMap ON? %d %d\n",recou->IsBadChannelsRemovalSwitchedOn(), cu->IsBadChannelsRemovalSwitchedOn());
728 
729  // PHOS
731 
732  if(printSettings) cu->Print("");
733 
734  return cu;
735 }
736 
742  TString calorimeter, Int_t year, Int_t tm,
743  Bool_t printSettings, Int_t debug)
744 {
745  AliAnaPhoton *ana = new AliAnaPhoton();
746 
747  // cluster selection cuts
748 
750 
751  ana->SwitchOffFiducialCut();
752 
753  ana->SetCalorimeter(calorimeter);
754  if(calorimeter == "DCAL")
755  {
756  TString calo = "EMCAL";
757  ana->SetCalorimeter(calo);
758  }
759 
760  ana->SwitchOnFillShowerShapeHistograms(); // Filled before photon shower shape selection
761 
762  //if(!simulation) ana->SwitchOnFillPileUpHistograms();
763 
764  if(tm) ana->SwitchOnTrackMatchRejection() ;
765  else ana->SwitchOffTrackMatchRejection() ;
766 
767  ana->SwitchOnTMHistoFill() ;
768 
769  if(calorimeter == "PHOS")
770  {
771  ana->SetNCellCut(2);// At least 3 cells
772  ana->SetMinPt(0.3);
773  ana->SetMinDistanceToBadChannel(2, 4, 5);
774  ana->SetTimeCut(-1e10,1e10); // open cut
775  }
776  else
777  {//EMCAL
778  ana->SetNCellCut(1);// At least 2 cells
779  ana->SetMinEnergy(0.3); // avoid mip peak at E = 260 MeV
780  ana->SetMaxEnergy(100);
781  ana->SetTimeCut(-1e10,1e10); // open cut, usual time window of [425-825] ns if time recalibration is off
782  // restrict to less than 100 ns when time calibration is on
783  ana->SetMinDistanceToBadChannel(2, 4, 6);
784 
785  // NLM cut, used in all, exclude clusters with more than 2 maxima
786  // Not needed if M02 cut is already strong or clusterizer V2
787  ana->SetNLMCut(1, 2) ;
788  }
789 
790  //PID cuts (shower shape)
791  ana->SwitchOnCaloPID(); // do PID selection, unless specified in GetCaloPID, selection not based on bayesian
792  AliCaloPID* caloPID = ana->GetCaloPID();
793  //Not used in bayesian
794 
795  // EMCAL
796 
797  //caloPID->SetEMCALLambda0CutMax(0.27);
798  caloPID->SetEMCALLambda0CutMax(10); // open, full shower shape needed for isolation studies
799  caloPID->SetEMCALLambda0CutMin(0.10);
800 
801  // Track matching
802  // tm = 1, fixed cuts
803  caloPID->SetEMCALDEtaCut(0.025);
804  caloPID->SetEMCALDPhiCut(0.030);
805 
806  // pT track dependent cuts
807  if(tm > 1) caloPID->SwitchOnEMCTrackPtDepResMatching();
808 
809  // PHOS
810  caloPID->SetPHOSDispersionCut(2.5);
811  caloPID->SetPHOSRCut(2.);
812  //if(kInputData=="AOD") caloPID->SetPHOSRCut(2000.); // Open cut since dX, dZ not stored
813 
814  // Branch AOD settings
815  ana->SetOutputAODName(Form("PhotonTrigger_%s",kAnaGammaHadronCorr.Data()));
816  ana->SetOutputAODClassName("AliCaloTrackParticleCorrelation");
817 
818  //Set Histograms name tag, bins and ranges
819  ana->AddToHistogramsName(Form("AnaPhoton_TM%d_",tm));
820 
821  SetAnalysisCommonParameters(ana,calorimeter,year,col,simulation,printSettings,debug) ; // see method below
822 
823  if(ana->GetFirstSMCoveredByTRD() > 0)
824  printf("AddTaskGammaHadronCorrelation() >>> Set first SM covered by TRD, SM=%d <<< year %d \n", ana->GetFirstSMCoveredByTRD(),year);
825 
826  // Number of particle type MC histograms
827  ana->FillNOriginHistograms (17); // 18 max
828  ana->FillNPrimaryHistograms(6); // 6 max
829 
830  return ana;
831 }
832 
838  Bool_t useSSIso, Bool_t useAsy,
839  TString col, Bool_t simulation,
840  TString calorimeter, Int_t year, Int_t tm,
841  Bool_t printSettings, Int_t debug )
842 {
843  // Configuration of pi0 event by event selection
844 
845  AliAnaPi0EbE *ana = new AliAnaPi0EbE();
846 
847  ana->SetDebug(debug);
848 
849  ana->SetAnalysisType((AliAnaPi0EbE::anaTypes)analysis);
850  TString opt = "";
851  if(analysis==AliAnaPi0EbE::kIMCaloTracks) opt = "Conv";
852  if(analysis==AliAnaPi0EbE::kSSCalo) opt = "SS";
853 
856 
858  //if(!simulation) ana->SwitchOnFillPileUpHistograms();
859  //if(!kTime && !simulation) ana->SwitchOnFillEMCALBCHistograms();
860 
861  if(tm) ana->SwitchOnTrackMatchRejection() ;
862  else ana->SwitchOffTrackMatchRejection() ;
863  ana->SwitchOffTMHistoFill() ;
864 
865  ana->SetCalorimeter(calorimeter);
866  if(calorimeter == "DCAL")
867  {
868  TString calo = "EMCAL";
869  ana->SetCalorimeter(calo);
870  }
871 
872  // Branch AOD settings
873  ana->SetOutputAODName(Form("%s%sTrigger_%s",particle.Data(), opt.Data(), kAnaGammaHadronCorr.Data()));
874  printf("***Out branch %s***\n",ana->GetOutputAODName().Data());
875  ana->SetOutputAODClassName("AliCaloTrackParticleCorrelation");
876 
877  if(analysis == AliAnaPi0EbE::kIMCaloTracks) ana->SetInputAODGammaConvName("PhotonsCTS");
878 
879  //Set Histograms name tag, bins and ranges
880 
881  ana->AddToHistogramsName(Form("Ana%s%sEbE_TM%d_",particle.Data(),opt.Data(),tm));
882 
883  SetAnalysisCommonParameters(ana,calorimeter,year,col,simulation,printSettings,debug); // see method below
884 
886  if(analysis!=AliAnaPi0EbE::kSSCalo)
887  {
888  ana->SetInputAODName(Form("PhotonTrigger_%s",kAnaGammaHadronCorr.Data()));
889 
890  ana->SetM02CutForInvMass(0.1,0.35); // Loose SS cut
891 
893  ana->SetR(0.4);
895 
896  if(useSSIso)
897  {
899  ana->AddToHistogramsName(Form("Ana%s%sEbEIsoDecay_TM%d_",particle.Data(),opt.Data(),tm));
900  ana->SetOutputAODName(Form("%s%sIsoDecayTrigger_%s",particle.Data(), opt.Data(), kAnaGammaHadronCorr.Data()));
901  }
902 
903  if ( calorimeter.Contains("CAL") && !simulation ) ana->SetPairTimeCut(100);
904 
906  nms->SetParticle(particle);
907 
908  //****
909  nms->SetInvMassCutMaxParameters(0,0,0); // Overrule the setting in SetParticle for Pi0 option
910  //****
911 
912  // Tighten a bit mass cut with respect to default window
913  if(particle=="Pi0") nms->SetInvMassCutRange(0.110,0.160);
914  if(particle=="Eta") nms->SetInvMassCutRange(0.520,0.580);
915 
916  //if(!particle.Contains("SideBand")) nms->SwitchOnAngleSelection();
917  //else nms->SwitchOnAngleSelection();
918 
920 
921  if(particle.Contains("Pi0SideBand")) // For pi0, do not consider left band
922  nms->SetSideBandCutRanges(-1,0,0.190,0.240);
923 
924  if(particle.Contains("EtaSideBand")) // For pi0, do not consider left band
925  nms->SetSideBandCutRanges(0.410,0.470,0.620,0.680);
926 
928  //nms->SetAngleMaxParam(2,0.2);
929  nms->SetHistoERangeAndNBins(0, 20, 80) ;
930  //nms->SetHistoIMRangeAndNBins(0, 1, 400);
931  }
932  else
933  {
934  // cluster splitting settings
935  ana->SetMinEnergy(6);
936  ana->SetMaxEnergy(100.);
937 
938  ana->SetNLMMinEnergy(0, 10);
939  ana->SetNLMMinEnergy(1, 6);
940  ana->SetNLMMinEnergy(2, 6);
941 
942  // NLM cut, used in all, exclude clusters with more than 2 maxima
943  ana->SetNLMCut(1, 2) ;
944 
945  //
946  ana->SetMinDistanceToBadChannel(2, 4, 6);
948  ana->SetTimeCut(-1e10,1e10); // Open time cut
949 
950  AliCaloPID* caloPID = ana->GetCaloPID();
951 
952  // Track matching
953  // tm = 1, fixed cuts
954  caloPID->SetEMCALDEtaCut(0.025);
955  caloPID->SetEMCALDPhiCut(0.030);
956 
957  // pT track dependent cuts
958  if(tm > 1) caloPID->SwitchOnEMCTrackPtDepResMatching();
959 
960  caloPID->SetSplitWidthSigma(3); // cut at 3 sigma of the mean pi0 peak.
961 
962  if(!useSSIso)
963  {
964  printf("Do not apply SS cut on merged pi0 analysis \n");
965  caloPID->SwitchOffSplitShowerShapeCut() ;
966  ana->AddToHistogramsName(Form("Ana%s%sEbE_OpenSS_TM%d_",particle.Data(),opt.Data(),tm));
967  ana->SetOutputAODName(Form("%s%sTrigger_%s_OpenSS",particle.Data(), opt.Data(), kAnaGammaHadronCorr.Data()));
968  caloPID->SetClusterSplittingM02Cut(0.1,10);
969  }
970  else
971  {
972  caloPID->SetClusterSplittingM02Cut(0.3,4); // Do the selection in the analysis class and not in the PID method to fill SS histograms
973  caloPID->SwitchOnSplitShowerShapeCut() ;
974  }
975 
976  if(useAsy)
977  {
978  caloPID->SwitchOnSplitAsymmetryCut() ;
980  ana->GetCaloPID()->SetSubClusterEnergyMinimum(1,0.5);
981  ana->GetCaloPID()->SetSubClusterEnergyMinimum(2,0.5);
982  }
983  else
984  {
985  caloPID->SwitchOffSplitAsymmetryCut() ;
986  if(!useSSIso)
987  {
988  ana->AddToHistogramsName(Form("Ana%s%sEbE_OpenSS_OpenAsy_TM%d_",particle.Data(),opt.Data(),tm));
989  ana->SetOutputAODName(Form("%s%sTrigger_%s_OpenSS_OpenAsy",particle.Data(), opt.Data(), kAnaGammaHadronCorr.Data()));
990  }
991  else
992  {
993  ana->AddToHistogramsName(Form("Ana%s%sEbE_OpenAsy_TM%d_",particle.Data(),opt.Data(),tm));
994  ana->SetOutputAODName(Form("%s%sTrigger_%s_OpenAsy",particle.Data(), opt.Data(), kAnaGammaHadronCorr.Data()));
995  }
996  }
997 
998  // For Pi0 only if SwitchOnSimpleSplitMassCut()
999  caloPID->SetPi0MassRange(0.10, 0.18);
1000  caloPID->SetEtaMassRange(0.50, 0.60);
1001  caloPID->SetPhotonMassRange(0.00, 0.08);
1002 
1003  caloPID->SetClusterSplittingMinNCells(6);
1004 
1005  //caloPID->SetSplitEnergyFractionMinimum(0, 0.95);
1006  //caloPID->SetSplitEnergyFractionMinimum(1, 0.95);
1007  //caloPID->SetSplitEnergyFractionMinimum(2, 0.8);
1008 
1009  if(col=="PbPb" || kAnaGammaHadronCorr.Contains("150"))
1010  {
1011  caloPID->SetClusterSplittingMinNCells(4);
1012  //caloPID->SetPi0MassShiftHighECell(0.005);
1013  }
1014  }
1016 
1017  return ana;
1018 }
1019 
1024  Int_t partInCone, Int_t thresType,
1025  Float_t cone, Float_t pth, Bool_t multi,
1026  TString col, Bool_t simulation,
1027  TString calorimeter, Int_t year, Int_t tm,
1028  Bool_t printSettings, Int_t debug )
1029 {
1031 
1032  ana->SetDebug(debug);
1033 
1034  ana->SetMinPt(5);
1035  ana->SetCalorimeter(calorimeter);
1036  if(calorimeter == "DCAL")
1037  {
1038  TString calo = "EMCAL";
1039  ana->SetCalorimeter(calo);
1040  }
1041 
1043  ana->SwitchOffCellHistoFill() ;
1044 
1045  ana->SwitchOffLeadingOnly();
1047  if( leading > 0 ) ana->SwitchOnLeadingOnly();
1048  if( leading == 2 ||
1049  leading == 4) ana->SwitchOnCheckNeutralClustersForLeading();
1050 
1051  // MC
1054 
1055  if(particle.Contains("Photon"))
1056  {
1058  ana->SetNDecayBits(5);
1059  ana->SwitchOnSSHistoFill();
1060  }
1061  else
1062  {
1063  ana->SwitchOffSSHistoFill();
1064  }
1065 
1067  ana->SetNPtTrigBins(6);
1068  //ana->SetPtTrigLimits(0,8); ana->SetPtTrigLimits(1,12); ana->SetPtTrigLimits(2,16); ana->SetPtTrigLimits(3,25);
1069 
1071  ana->SetNBackgroundBins(11);
1072  //ana->SetBackgroundLimits(0,0); ana->SetBackgroundLimits(1,0.2); ana->SetBackgroundLimits(2,3); ana->SetBackgroundLimits(3,0.4);
1073 
1074  if(!tm) ana->SwitchOnTMHistoFill();
1075  else ana->SwitchOffTMHistoFill();
1076 
1077  //if(!simulation) ana->SwitchOnFillPileUpHistograms();
1078 
1080  ana->SwitchOnFiducialCut();
1081 
1082  // Avoid borders of calorimeter
1083  if ( calorimeter == "EMCAL" ) ana->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.60, 86, 174) ;
1084  else if ( calorimeter == "DCAL" ) ana->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.60, 266, 314) ;
1085  else if ( calorimeter == "PHOS" ) ana->GetFiducialCut()->SetSimplePHOSFiducialCut (0.10, 266, 314) ;
1086 
1087  // Same Eta as EMCal, cut in phi if EMCAL was triggering
1088  if(particle=="Hadron" || particle.Contains("CTS"))
1089  {
1090  //if(trigger.Contains("EMC"))
1091  // ana->GetFiducialCut()->SetSimpleCTSFiducialCut (0.6, 260, 360) ;
1092  //else
1093  ana->GetFiducialCut()->SetSimpleCTSFiducialCut (0.6, 0, 360) ;
1094  }
1095 
1096  // Branch AOD settings
1097 
1098  ana->SetInputAODName(Form("%sTrigger_%s",particle.Data(),kAnaGammaHadronCorr.Data()));
1099  ana->SetAODObjArrayName(Form("IC%sTrigger_%s_R%1.1f_ThMin%1.1f",particle.Data(),kAnaGammaHadronCorr.Data(),cone,pth));
1100 
1101  //
1102  // Do settings for main isolation cut class
1103  //
1104  AliIsolationCut * ic = ana->GetIsolationCut();
1105  ic->SetDebug(debug);
1106  ic->SetParticleTypeInCone(partInCone);
1107  ic->SetICMethod(thresType);
1108  ic->SetPtFraction(0.1);
1109  ic->SetPtThreshold(0.5); // default, change in next lines
1110  ic->SetSumPtThreshold(1.0); // default, change in next lines
1111 
1112  if(cone > 0 && pth > 0)
1113  {
1114  ic->SetConeSize(cone);
1115  ic->SetPtThresholdMax(10000);
1116 
1117  if(thresType == AliIsolationCut::kPtThresIC)
1118  {
1119  printf("*** Iso *** PtThresMin = %1.1f GeV/c *** R = %1.1f ***\n",pth,cone);
1120  ic->SetPtThreshold(pth);
1121  }
1122 
1123  if(thresType == AliIsolationCut::kSumPtIC)
1124  {
1125  printf("*** Iso *** SumPtMin = %1.1f GeV/c *** R = %1.1f ***\n",pth,cone);
1126  ic->SetSumPtThreshold(pth);
1127  }
1128  }
1129  else
1130  {
1131  if(col=="pp")
1132  {
1133  ic->SetPtThreshold(0.5);
1134  ic->SetSumPtThreshold(1.0) ;
1135  ic->SetConeSize(0.4);
1136  }
1137  if(col=="PbPb")
1138  {
1139  ic->SetPtThreshold(3.);
1140  ic->SetSumPtThreshold(3.0) ;
1141  ic->SetConeSize(0.3);
1142  }
1143  }
1144 
1145 
1146  // Do or not do isolation with previously produced AODs.
1147  // No effect if use of SwitchOnSeveralIsolation()
1148  ana->SwitchOffReIsolation();
1149 
1150  // Multiple IC
1151  if(multi)
1152  {
1153  ic->SetConeSize(1.); // Take all for first iteration
1154  ic->SetPtThreshold(100);// Take all for first iteration
1155  ana->SwitchOnSeveralIsolation() ;
1156  ana->SetAODObjArrayName(Form("MultiIC%sTM%d",particle.Data(),tm));
1157 
1158  ana->SetNCones(3);
1159  ana->SetNPtThresFrac(2);
1160  ana->SetConeSizes(0,0.3); ana->SetConeSizes(1,0.4); ana->SetConeSizes(2,0.5);
1161  ana->SetPtThresholds(0, 0.5); ana->SetPtThresholds(1, 1); ana->SetPtThresholds(2, 1.5); ana->SetPtThresholds(3, 2);
1162  ana->SetPtFractions (0, 0.05) ; ana->SetPtFractions (1, 0.1); ana->SetPtFractions (2, 0.2) ; ana->SetPtFractions (3, 0.3) ;
1163  ana->SetSumPtThresholds(0, 0.5) ; ana->SetSumPtThresholds(1, 1) ; ana->SetSumPtThresholds(2, 1.5); ana->SetSumPtThresholds(3, 2) ;
1164  //ana->SetPtThresholds(0, 0.5);
1165 
1166  ana->SwitchOffTMHistoFill();
1167  ana->SwitchOffSSHistoFill();
1168  }
1169  else
1170  ana->SwitchOffSeveralIsolation() ;
1171 
1172  AliCaloPID* caloPID = ana->GetCaloPID();
1173 
1174  // Track matching
1175  // tm = 1, fixed cuts
1176  caloPID->SetEMCALDEtaCut(0.025);
1177  caloPID->SetEMCALDPhiCut(0.030);
1178 
1179  // pT track dependent cuts
1180  if(tm > 1) caloPID->SwitchOnEMCTrackPtDepResMatching();
1181 
1182  //Set Histograms name tag, bins and ranges
1183 
1184  if(!multi) ana->AddToHistogramsName(Form("AnaIsol%s_TM%d_",particle.Data(),tm));
1185  else ana->AddToHistogramsName(Form("AnaMultiIsol%s_TM%d_",particle.Data(),tm));
1186 
1187  SetAnalysisCommonParameters(ana,calorimeter,year,col,simulation,printSettings,debug); // see method below
1188 
1189  if(particle=="Hadron" || particle.Contains("CTS"))
1190  {
1191  ana->GetHistogramRanges()->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 200) ;
1192  ana->GetHistogramRanges()->SetHistoEtaRangeAndNBins(-1.5, 1.5, 300) ;
1193  }
1194 
1195  if(printSettings) ic ->Print("");
1196 
1197  return ana;
1198 }
1199 
1200 
1205  Bool_t bIsolated, Float_t shshMax,
1206  Int_t partInCone, Int_t thresType,
1207  Float_t cone, Float_t pth, Bool_t mixOn,
1208  TString col, Bool_t simulation,
1209  TString calorimeter, Int_t year, Int_t tm,
1210  Bool_t printSettings, Int_t debug )
1211 {
1213 
1214  ana->SetTriggerPtRange(5,100);
1215  ana->SetAssociatedPtRange(0.2,100);
1216  ana->SetDeltaPhiCutRange (TMath::DegToRad()*120.,TMath::DegToRad()*240.);
1217 
1218  // Underlying event
1219  ana->SetUeDeltaPhiCutRange(TMath::DegToRad()*60. ,TMath::DegToRad()*120.);
1221 
1222  ana->SwitchOffAbsoluteLeading(); // Select trigger leading particle of all the selected tracks
1223  ana->SwitchOffNearSideLeading(); // Select trigger leading particle of all the particles at +-90 degrees, default
1225 
1226  if(leading > 0 && leading < 3 ) ana->SwitchOnAbsoluteLeading();
1227  if(leading > 2 ) ana->SwitchOnNearSideLeading();
1228  if(leading == 2 || leading == 4 ) ana->SwitchOnCheckNeutralClustersForLeading();
1229 
1231  ana->SwitchOffCorrelationVzBin() ;
1233 
1235 
1237  if(particle.Contains("Photon"))
1238  {
1240  ana->SetNDecayBits(5);
1241  printf("**** SET M02 limits in correlation task *** \n");
1242  ana->SetM02Cut(0.10,shshMax);
1245  }
1246 
1247  ana->SetMCGenType(0,7);
1248 
1249  ana->SwitchOffLeadHadronSelection(); // Open cuts, just fill histograms
1253  ana->SetLeadHadronPhiCut(TMath::DegToRad()*130, TMath::DegToRad()*230.);
1254  ana->SetLeadHadronPtCut(0.5, 1000);
1255 
1256  // if triggering on PHOS and EMCAL is on
1257  ana->SwitchOffNeutralCorr(); // Do only correlation with TPC
1258  //ana->SetPi0AODBranchName("Pi0EMCAL_TrigEMC7_Cl_TM1");
1259 
1261 
1263 
1264  //if(!simulation) ana->SwitchOnFillPileUpHistograms();
1265 
1267 
1268  ana->SetNAssocPtBins(16); // set last bin [20,30] GeV/c
1269  // See AliAnaParticleCorrelation::InitParameters();
1270  // Default bins{0.2,0.5,1,2,3,4,5,6,7,8,9,10,12,14,16,20,30,40,50,100} GeV/c
1271  // If you want to change it:
1272  // ana->SetAssocPtBinLimit(0, 1) ;
1273  // ana->SetAssocPtBinLimit(1, 2) ;
1274  // ana->SetAssocPtBinLimit(2, 3) ;
1275  // ana->SetAssocPtBinLimit(3, 4) ;
1276  // ana->SetAssocPtBinLimit(4, 5) ;
1277  // ana->SetAssocPtBinLimit(5, 8) ;
1278  // ana->SetAssocPtBinLimit(6, 10) ;
1279  // ana->SetAssocPtBinLimit(7, 100);
1280 
1281  ana->SetNTriggerPtBins(5); // set last bin [25,30] GeV/c
1282  // See AliAnaParticleCorrelation::InitParameters();
1283  // Default bins{10,12,16,20,25,30,40,50,75,100} GeV/c
1284  // If you want to change it:
1285  // ana->SetTriggerPtBinLimit(0, 5) ;
1286  // ana->SetTriggerPtBinLimit(1, 10) ;
1287 
1288  ana->SelectIsolated(bIsolated); // do correlation with isolated photons
1289 
1290  // Mixing with own pool
1291  if(mixOn)
1292  {
1293  ana->SwitchOnOwnMix();
1295 
1296  if(bIsolated)
1297  {
1298  //Do settings for main isolation cut class
1299  AliIsolationCut * ic = ana->GetIsolationCut();
1300  ic->SetDebug(debug);
1301 
1302  if(cone >0 && pth > 0)
1303  {
1304  printf("*** Correl *** PtThres = %1.1f GeV/c *** R = %1.1f ***\n",pth,cone);
1305  ic->SetPtThreshold(pth);
1306  ic->SetConeSize(cone);
1307  }
1308  else
1309  {
1310  if(col=="pp")
1311  {
1312  ic->SetPtThreshold(0.5);
1313  ic->SetConeSize(0.4);
1314  }
1315  if(col=="PbPb")
1316  {
1317  ic->SetPtThreshold(3.);
1318  //ic->SetPtThreshold(1.);
1319  ic->SetConeSize(0.3);
1320  }
1321  }
1322 
1323  ic->SetPtFraction(0.1);
1324  ic->SetSumPtThreshold(1.0) ;
1325  ic->SetParticleTypeInCone(partInCone);
1326  ic->SetICMethod(thresType);
1327  }
1328  }
1329  else
1330  ana->SwitchOffOwnMix();
1331 
1332  ana->SetNZvertBin(20);
1333 
1334  if(col=="pp")
1335  {
1336  ana->SetNMaxEvMix(100);
1337  ana->SwitchOnTrackMultBins();
1338  ana->SetNTrackMultBin(10);
1339  ana->SetNRPBin(1);
1340  }
1341  else
1342  {
1343  ana->SetNMaxEvMix(10);
1344  ana->SwitchOffTrackMultBins(); // centrality bins
1345  ana->SetNCentrBin(12);
1346  ana->SetNRPBin(3);
1347  if(kAnaGammaHadronCorr.Contains("60_90"))
1348  {
1349  printf("*** Set mixing for peripheral\n");
1350  ana->SetNMaxEvMix(50);
1351  ana->SetNCentrBin(2);
1352  }
1353  }
1354 
1355  ana->SwitchOnFiducialCut();
1356 
1357  // Avoid borders of calorimeter, same as for isolation
1358  if ( calorimeter == "EMCAL" ) ana->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.60, 86, 174) ;
1359  else if ( calorimeter == "DCAL" ) ana->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.60, 266, 314) ;
1360  else if ( calorimeter == "PHOS" ) ana->GetFiducialCut()->SetSimplePHOSFiducialCut (0.10, 266, 314) ;
1361 
1362  // Input / output delta AOD settings
1363 
1364  ana->SetInputAODName(Form("%sTrigger_%s",particle.Data(),kAnaGammaHadronCorr.Data()));
1365  ana->SetAODObjArrayName(Form("%sHadronCorrIso%dTrigger_%s",particle.Data(),bIsolated,kAnaGammaHadronCorr.Data()));
1366  //ana->SetAODNamepTInConeHisto(Form("IC%s_%s_R%1.1f_ThMin%1.1f" ,particle.Data(),kAnaGammaHadronCorr.Data(),cone,pth));
1367 
1368  //Set Histograms name tag, bins and ranges
1369 
1370  ana->AddToHistogramsName(Form("Ana%sHadronCorr_Iso%d_TM%d_",particle.Data(),bIsolated,tm));
1371 
1372  SetAnalysisCommonParameters(ana,calorimeter,year,col,simulation,printSettings,debug); // see method below
1373 
1374  if(particle=="Hadron" || particle.Contains("CTS"))
1375  {
1376  ana->GetHistogramRanges()->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 200) ;
1377  ana->GetHistogramRanges()->SetHistoEtaRangeAndNBins(-1.5, 1.5, 300) ;
1378  }
1379 
1380  return ana;
1381 }
1382 
1387 {
1389 
1390  ana->SetDebug(debug);
1391 
1392  // selection cuts
1393 
1394  ana->SetMinPt(0.2);
1395  ana->SwitchOnFiducialCut();
1396  ana->GetFiducialCut()->SetSimpleCTSFiducialCut(0.8, 0, 360) ; //more restrictive cut in reader and after in isolation
1397 
1399 
1402 
1403  // Branch AOD settings
1404 
1405  ana->SetOutputAODName(Form("HadronTrigger_%s",kAnaGammaHadronCorr.Data()));
1406  ana->SetOutputAODClassName("AliCaloTrackParticle"); // use if no correlation done
1407 
1408  //Set Histograms name tag, bins and ranges
1409 
1410  ana->AddToHistogramsName("AnaHadrons_");
1411 
1412  SetAnalysisCommonParameters(ana,"CTS",2012,"pp",simulation,printSettings,debug); // see method below
1413 
1414  return ana;
1415 }
1416 
1421  TString calorimeter, Int_t year,
1422  Bool_t printSettings, Int_t debug )
1423 {
1425 
1426  ana->SetCalorimeter(calorimeter);
1427  if(calorimeter == "DCAL")
1428  {
1429  TString calo = "EMCAL";
1430  ana->SetCalorimeter(calo);
1431  }
1432 
1433  ana->SetTimeCut(-1e10,1e10); // Open time cut
1434 
1435  ana->SwitchOffCorrelation(); // make sure you switch in the reader PHOS and EMCAL cells and clusters if option is ON
1436 
1438 
1439  ana->SwitchOffFiducialCut();
1443  ana->SwitchOffStudyBadClusters() ;
1445 
1447 
1448  ana->AddToHistogramsName("QA_"); // Begining of histograms name
1449 
1450  SetAnalysisCommonParameters(ana,calorimeter,year,col,simulation,printSettings,debug); // see method below
1451 
1452  return ana;
1453 }
1454 
1459  Float_t pth,
1460  TString col, Bool_t simulation,
1461  TString calorimeter, Int_t year,
1462  Bool_t printSettings, Int_t debug )
1463 {
1465 
1466  // Trigger detector, acceptance and pT cut
1467  ana->SetTriggerDetector(calorimeter);
1468  if(calorimeter == "DCAL")
1469  {
1470  TString calo = "EMCAL";
1471  ana->SetTriggerDetector(calo);
1472  }
1473 
1474  ana->SetMinPt(2); // Trigger photon, pi0 minimum pT
1475  if ( calorimeter == "EMCAL" ) ana->GetFiducialCutForTrigger()->SetSimpleEMCALFiducialCut(0.60, 86, 174) ;
1476  else if ( calorimeter == "DCAL" ) ana->GetFiducialCutForTrigger()->SetSimpleEMCALFiducialCut(0.60, 266, 314) ;
1477  else if ( calorimeter == "PHOS" ) ana->GetFiducialCutForTrigger()->SetSimplePHOSFiducialCut (0.10, 266, 314) ;
1478 
1479  // Particles associated to trigger or isolation cone acceptance and pT cut
1480  ana->SetCalorimeter(calorimeter);
1481  if(calorimeter == "DCAL")
1482  {
1483  TString calo = "EMCAL";
1484  ana->SetCalorimeter(calo);
1485  }
1486 
1487  ana->SetMinChargedPt(0.2);
1488  ana->SetMinNeutralPt(0.3);
1489 
1490  if ( calorimeter == "EMCAL" ) ana->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.70, 80, 174) ;
1491  else if ( calorimeter == "DCAL" ) ana->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.70, 260, 327) ;
1492  else if ( calorimeter == "PHOS" ) ana->GetFiducialCut()->SetSimplePHOSFiducialCut (0.12, 250, 320) ;
1493  ana->GetFiducialCut()->SetSimpleCTSFiducialCut(0.9, 0, 360);
1494 
1495  // Isolation paramters
1496  AliIsolationCut * ic = ana->GetIsolationCut();
1497  ic->SetDebug(debug);
1498  ic->SetPtThreshold(pth);
1499  ic->SetConeSize(cone);
1500  ic->SetSumPtThreshold(1.0) ;
1501  ic->SetICMethod(thresType);
1502 
1503  ana->AddToHistogramsName("AnaGenKine_");
1504 
1505  SetAnalysisCommonParameters(ana,calorimeter,year,col,simulation,printSettings,debug); // see method below
1506 
1507  return ana;
1508 }
1509 
1510 
1518  TString calorimeter, Int_t year,
1519  TString col, Bool_t simulation,
1520  Bool_t printSettings, Int_t debug)
1521 {
1522  //
1523  // Histograms ranges
1524  //
1525  AliHistogramRanges* histoRanges = ana->GetHistogramRanges();
1526 
1527  histoRanges->SetHistoPtRangeAndNBins(0, 100, 200) ; // Energy and pt histograms
1528 
1529  if(calorimeter=="EMCAL")
1530  {
1531  ana->SetFirstSMCoveredByTRD(-1);
1532 
1533  if ( year == 2010 )
1534  {
1535  histoRanges->SetHistoPhiRangeAndNBins(78*TMath::DegToRad(), 122*TMath::DegToRad(), 78) ;
1536  histoRanges->SetHistoXRangeAndNBins(-230,90,120); // QA
1537  histoRanges->SetHistoYRangeAndNBins(370,450,40); // QA
1538  }
1539  else if ( year < 2014 )
1540  {
1541  histoRanges->SetHistoPhiRangeAndNBins(78*TMath::DegToRad(), 182*TMath::DegToRad(), 104) ;
1542  histoRanges->SetHistoXRangeAndNBins(-460,90,200); // QA
1543  histoRanges->SetHistoYRangeAndNBins(100,450,100); // QA
1544 
1545  if (year == 2011) ana->SetFirstSMCoveredByTRD( 6);
1546  else if(year == 2012 ||
1547  year == 2013) ana->SetFirstSMCoveredByTRD( 4);
1548  }
1549  else // Run2
1550  {
1551  histoRanges->SetHistoPhiRangeAndNBins(78*TMath::DegToRad(), 189*TMath::DegToRad(), 111) ;
1552  histoRanges->SetHistoXRangeAndNBins(-460,460,230); // QA, revise
1553  histoRanges->SetHistoYRangeAndNBins(-450,450,225); // QA, revise
1554  }
1555 
1556  histoRanges->SetHistoEtaRangeAndNBins(-0.72, 0.72, 144) ;
1557  }
1558  else if(calorimeter=="DCAL")
1559  {
1560  histoRanges->SetHistoPhiRangeAndNBins(260*TMath::DegToRad(), 327*TMath::DegToRad(), 67) ;
1561  histoRanges->SetHistoEtaRangeAndNBins(-0.72, 0.72, 144) ;
1562  histoRanges->SetHistoXRangeAndNBins(-460,460,230); // QA, revise
1563  histoRanges->SetHistoYRangeAndNBins(-450,450,225); // QA, revise
1564  ana->SetFirstSMCoveredByTRD( 0);
1565  }
1566  else if(calorimeter=="PHOS")
1567  {
1568  histoRanges->SetHistoPhiRangeAndNBins(250*TMath::DegToRad(), 320*TMath::DegToRad(), 70) ;
1569  histoRanges->SetHistoEtaRangeAndNBins(-0.13, 0.13, 130) ;
1570  }
1571  else if(calorimeter=="CTS")
1572  {
1573  ana->GetHistogramRanges()->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 200) ;
1574  ana->GetHistogramRanges()->SetHistoEtaRangeAndNBins(-1.5, 1.5, 300) ;
1575  }
1576 
1577  histoRanges->SetHistoShowerShapeRangeAndNBins(-0.1, 4.9, 500);
1578 
1579  // Invariant mass histo
1580  histoRanges->SetHistoMassRangeAndNBins(0., 1., 200) ;
1581  histoRanges->SetHistoAsymmetryRangeAndNBins(0., 1. , 100) ;
1582 
1583  // check if time calibration is on
1584  //histoRanges->SetHistoTimeRangeAndNBins(-1000.,1000,1000);
1585  histoRanges->SetHistoTimeRangeAndNBins(-400.,400,400);
1586  histoRanges->SetHistoDiffTimeRangeAndNBins(-200, 200, 800);
1587 
1588  // track-cluster residuals
1589  histoRanges->SetHistoTrackResidualEtaRangeAndNBins(-0.15,0.15,300);
1590  histoRanges->SetHistoTrackResidualPhiRangeAndNBins(-0.15,0.15,300);
1591  histoRanges->SetHistodRRangeAndNBins(0.,0.15,150);//QA
1592 
1593  // QA, electron, charged
1594  histoRanges->SetHistoPOverERangeAndNBins(0,2.,200);
1595  histoRanges->SetHistodEdxRangeAndNBins(0.,200.,200);
1596 
1597  // QA
1598  histoRanges->SetHistoFinePtRangeAndNBins(0, 10, 200) ; // bining for fhAmpId
1599  histoRanges->SetHistoVertexDistRangeAndNBins(0.,500.,500);
1600  histoRanges->SetHistoZRangeAndNBins(-350,350,175);
1601  histoRanges->SetHistoRRangeAndNBins(430,460,30);
1602  histoRanges->SetHistoV0SignalRangeAndNBins(0,5000,500);
1603  histoRanges->SetHistoV0MultiplicityRangeAndNBins(0,5000,500);
1604 
1605  // QA, correlation
1606  if(col=="PbPb")
1607  {
1608  histoRanges->SetHistoNClusterCellRangeAndNBins(0,100,100);
1609  histoRanges->SetHistoNClustersRangeAndNBins(0,500,50);
1610  histoRanges->SetHistoTrackMultiplicityRangeAndNBins(0,2000,200);
1611  }
1612  else
1613  {
1614  histoRanges->SetHistoNClusterCellRangeAndNBins(0,50,50);
1615  histoRanges->SetHistoNClustersRangeAndNBins(0,50,50);
1616  histoRanges->SetHistoTrackMultiplicityRangeAndNBins(0,200,200);
1617  }
1618 
1619  // xE, zT
1620  histoRanges->SetHistoRatioRangeAndNBins(0.,2.,200);
1621  histoRanges->SetHistoHBPRangeAndNBins (0.,10.,200);
1622 
1623  // Isolation
1624  histoRanges->SetHistoPtInConeRangeAndNBins(0, 50 , 250);
1625  histoRanges->SetHistoPtSumRangeAndNBins (0, 100, 250);
1626 
1627  //
1628  // MC histograms?
1629  //
1630  if(simulation) ana->SwitchOnDataMC() ;//Access MC stack and fill more histograms, AOD MC not implemented yet.
1631  else ana->SwitchOffDataMC() ;
1632 
1633  //
1634  // Debug
1635  //
1636  if(printSettings) ana->Print("");
1637 
1638  ana->SetDebug(debug); // 10 for lots of messages
1639 }
void SetSumPtThresholds(Int_t i, Float_t pt)
void SwitchOffSelectedClusterHistoFill()
Definition: AliAnaPi0EbE.h:112
void SetPtThreshold(Float_t pt)
AliAnalysisTaskCaloTrackCorrelation * AddTaskGammaHadronCorrelation(TString calorimeter="EMCAL", Bool_t simulation=kFALSE, Int_t year=2011, TString col="pp", Int_t rejectEMCTrig=0, TString clustersArray="", Bool_t tender=kFALSE, Bool_t nonLinOn=kFALSE, Float_t shshMax=0.27, Float_t isoCone=0.4, Float_t isoPtTh=0.5, Int_t isoMethod=AliIsolationCut::kPtThresIC, Int_t isoContent=AliIsolationCut::kNeutralAndCharged, Int_t leading=0, Int_t tm=1, Int_t minCen=-1, Int_t maxCen=-1, Bool_t mixOn=kTRUE, Bool_t qaAn=kFALSE, Bool_t chargedAn=kFALSE, TString outputfile="", Bool_t printSettings=kFALSE, Int_t debug=0, const char *trigSuffix="EMC7")
virtual void SwitchOnAODTrackSharedClusterSelection()
void SetHistoTrackMultiplicityRangeAndNBins(Int_t min, Int_t max, Int_t n)
virtual void SetDebug(Int_t d)
AliAnaParticleHadronCorrelation * ConfigureHadronCorrelationAnalysis(TString particle, Int_t leading, Bool_t bIsolated, Float_t shshMax, Int_t partInCone, Int_t thresType, Float_t cone, Float_t pth, Bool_t mixOn, TString col, Bool_t simulation, TString calorimeter, Int_t year, Int_t tm, Bool_t printSettings, Int_t debug)
void SetNumberOfCellsFromPHOSBorder(Int_t n)
void SetTimeCut(Double_t min, Double_t max)
void SetSimpleCTSFiducialCut(Float_t abseta, Float_t phimin, Float_t phimax)
void SetHistoMassRangeAndNBins(Float_t min, Float_t max, Int_t n)
virtual void AddToHistogramsName(TString add)
void SwitchOnFillAllTrackMatchingHistogram()
void SetFirstSuperModuleUsed(Int_t first)
void SetSplitWidthSigma(Float_t s)
Definition: AliCaloPID.h:283
virtual void SwitchOnFiducialCut()
void SetHistoXRangeAndNBins(Float_t min, Float_t max, Int_t n)
void SwitchOnTrackMatchRejection()
Definition: AliAnaPhoton.h:134
AliAnaChargedParticles * ConfigureChargedAnalysis(Bool_t simulation, Bool_t printSettings, Int_t debug)
void SetClusterSplittingM02Cut(Float_t min=0, Float_t max=100)
Definition: AliCaloPID.h:264
void SwitchOnLoadOwnEMCALGeometryMatrices()
void SwitchOnSplitShowerShapeCut()
Definition: AliCaloPID.h:260
void SwitchOnTMHistoFill()
Definition: AliAnaPhoton.h:96
void SetAssociatedPtRange(Float_t min, Float_t max)
Class with utils to perform Isolation Cuts.
virtual void SwitchOnAODHybridTrackSelection()
virtual void SetCentralityBin(Int_t min, Int_t max)
Int_t GetLastSuperModuleUsed() const
AliEMCALRecoUtils * GetEMCALRecoUtils() const
void SetLastSuperModuleUsed(Int_t last)
TString kAnaGammaHadronCorr
Global name to be composed of the settings, used to set the AOD branch name.
void SetConeSize(Float_t r)
void SetPtFraction(Float_t pt)
void SwitchOnSelectPairInIsolationCone()
Definition: AliAnaPi0EbE.h:127
void SetTimeCut(Double_t min, Double_t max)
Definition: AliAnaPi0EbE.h:142
void FillNPrimaryHistograms(Int_t n)
Definition: AliAnaPhoton.h:145
void SwitchOffClusterELinearityCorrection()
void SetCTSPtMax(Float_t pt)
void SwitchOffSplitShowerShapeCut()
Definition: AliCaloPID.h:261
void SetClusterSplittingMinNCells(Int_t c)
Definition: AliCaloPID.h:267
void SetHistoPtRangeAndNBins(Float_t min, Float_t max, Int_t n)
void SetDebug(Int_t d)
void SetMinChargedPt(Float_t pt)
AliAnaParticleIsolation * ConfigureIsolationAnalysis(TString particle, Int_t leading, Int_t partInCone, Int_t thresType, Float_t cone, Float_t pth, Bool_t multi, TString col, Bool_t simulation, TString calorimeter, Int_t year, Int_t tm, Bool_t printSettings, Int_t debug)
virtual void SetInputAODName(TString name)
void SetM02CutForInvMass(Float_t min=0, Float_t max=10)
Definition: AliAnaPi0EbE.h:117
virtual void SetTrackComplementaryCuts(AliESDtrackCuts *)
virtual AliIsolationCut * GetIsolationCut()
void SetEMCALEMin(Float_t en)
void SetMinDistanceToBadChannel(Float_t m1, Float_t m2, Float_t m3)
Definition: AliAnaPhoton.h:117
void SetInvMassCutRange(Double_t invmassmin, Double_t invmassmax)
void KeepNeutralMesonSelectionHistos(Bool_t keep)
void SetTriggerDetector(TString det)
Set the calorimeter for the analysis.
void SetSimpleEMCALFiducialCut(Float_t abseta, Float_t phimin, Float_t phimax)
AliCaloTrackReader * ConfigureReader(TString col, Bool_t simulation, TString clustersArray, Bool_t tender, TString calorimeter, Bool_t nonLinOn, TString trigger, Bool_t rejectEMCTrig, Int_t minCen, Int_t maxCen, Bool_t printSettings, Int_t debug)
virtual void SetOutputAODClassName(TString name)
virtual AliNeutralMesonSelection * GetNeutralMesonSelection()
void SetSumPtThreshold(Float_t s)
void SetMinDistanceToBadChannel(Float_t m1, Float_t m2, Float_t m3)
Definition: AliAnaPi0EbE.h:133
void ana(anaModes mode=mLocal)
Definition: ana.C:986
Class for the Calorimeter QA analysis.
void SetHistodRRangeAndNBins(Float_t min, Float_t max, Int_t n)
void SetLocalMaximaCutEDiff(Float_t c)
Get trigger particles/partons/jets and correlations at generator level.
void SwitchOffSplitAsymmetryCut()
Definition: AliCaloPID.h:257
void SetAnalysisType(anaTypes ana)
Definition: AliAnaPi0EbE.h:92
void SetNumberOfCellsFromEMCALBorder(Int_t n)
void SwitchOffTMHistoFill()
Definition: AliAnaPi0EbE.h:109
void SetNCellCut(Int_t n)
Definition: AliAnaPhoton.h:125
Select cluster pairs or single merged clusters with pi0 or eta invariant mass.
Definition: AliAnaPi0EbE.h:30
Bool_t IsBadChannelsRemovalSwitchedOn() const
void SetSimplePHOSFiducialCut(Float_t abseta, Float_t phimin, Float_t phimax)
void SetPtFractions(Int_t i, Float_t pt)
const TString calorimeter
Definition: anaM.C:36
virtual void SetPtHardAndJetPtComparison(Bool_t compare)
Bool_t IsRecalibrationOn() const
Bool_t IsRecalibrationOn() const
Some utilities for cluster and cell treatment.
void SetHistoTrackResidualEtaRangeAndNBins(Float_t min, Float_t max, Int_t n)
virtual void Print(const Option_t *opt) const
Print some relevant parameters set for the analysis.
void SetIsolationCandidateMinPt(Float_t min)
Definition: AliAnaPi0EbE.h:122
Class for event, clusters and tracks filtering and preparation for the AOD analysis.
void SetCTSPtMin(Float_t pt)
void SetPHOSRCut(Float_t rcut)
Definition: AliCaloPID.h:213
void SetEMCALDPhiCut(Float_t dcut)
Definition: AliCaloPID.h:221
void SetLeadHadronPtCut(Float_t min, Float_t max)
virtual void SetOutputAODName(TString name)
void SwitchOffRecalculateClusterTrackMatching()
void SetICMethod(Int_t i)
void SwitchOffVertexBCEventSelection()
void SetR(Float_t r)
Definition: AliAnaPi0EbE.h:121
virtual void SetCalorimeter(TString calo)
Set the calorimeter for the analysis. A string.
Correlate trigger particles (photon, pi0, tracks) and charged tracks: Azimuthal correlations, xE distributions.
virtual void SetNTrackMultBin(Int_t n=1)
Number of bins in centrality.
Base class for CaloTrackCorr analysis algorithms.
void SwitchOffTrackMatchRejection()
Definition: AliAnaPhoton.h:135
void SetDeltaPhiCutRange(Double_t phimin, Double_t phimax)
void SetUeDeltaPhiCutRange(Double_t uephimin, Double_t uephimax)
1 calorimeter cluster and 1 photon conversion pair invariant mass selection
Definition: AliAnaPi0EbE.h:88
void SetHistoPtInConeRangeAndNBins(Float_t min, Float_t max, Int_t n)
virtual void SetDebug(Int_t d)
void SetHistoPOverERangeAndNBins(Float_t min, Float_t max, Int_t n)
virtual AliFiducialCut * GetFiducialCut()
virtual AliFiducialCut * GetFiducialCut()
virtual void SetCentralityOpt(Int_t opt)
void SetHistoPhiRangeAndNBins(Float_t min, Float_t max, Int_t n)
int Int_t
Definition: External.C:63
void SetSideBandCutRanges(Double_t lmin, Double_t lmax, Double_t rmin, Double_t rmax)
AliAnaCalorimeterQA * ConfigureQAAnalysis(TString col, Bool_t simulation, TString calorimeter, Int_t year, Bool_t printSettings, Int_t debug)
virtual AliHistogramRanges * GetHistogramRanges()
unsigned int UInt_t
Definition: External.C:33
void SwitchOnSplitAsymmetryCut()
Definition: AliCaloPID.h:256
void SetHistoPtSumRangeAndNBins(Float_t min, Float_t max, Int_t n)
void SetAnalysisCommonParameters(AliAnaCaloTrackCorrBaseClass *ana, TString calorimeter, Int_t year, TString col, Bool_t simulation, Bool_t printSettings, Int_t debug)
float Float_t
Definition: External.C:68
void SetTrackStatus(ULong_t bit)
Bool_t IsBadChannelsRemovalSwitchedOn() const
Filter EMCal/PHOS clusters for photon analysis.
Definition: AliAnaPhoton.h:33
void SetHistoTrackResidualPhiRangeAndNBins(Float_t min, Float_t max, Int_t n)
void SetConeSizes(Int_t i, Float_t r)
void SwitchOnSelectIsolatedDecay()
Definition: AliAnaPi0EbE.h:124
void SetM02Cut(Float_t min=0, Float_t max=10)
void SetExoticCellDiffTimeCut(Float_t dt)
void SwitchOnSplitClusterDistToBad()
Definition: AliAnaPi0EbE.h:160
void GetAlienGlobalProductionVariables(Bool_t &simulation, TString &col, TString &period, Int_t &year, Bool_t print=kFALSE)
void SetEMCALClusterListName(TString name)
void SwitchOffTriggerClusterTimeRecal()
void SwitchOffTrackMatchRejection()
Definition: AliAnaPi0EbE.h:149
virtual void SetHistoERangeAndNBins(Float_t min, Float_t max, Int_t n)
Base class for event, clusters and tracks filtering and preparation for the analysis.
void SetInputAODGammaConvName(TString name)
Definition: AliAnaPi0EbE.h:98
void SetHistoNClusterCellRangeAndNBins(Int_t min, Int_t max, Int_t n)
void SwitchOffBadTriggerEventsRemoval()
void SetMinNeutralPt(Float_t pt)
anaTypes
Analysis types.
Definition: AliAnaPi0EbE.h:84
Int_t GetNumberOfSuperModulesUsed() const
void SetParticle(TString particleName)
Set some default parameters for selection of pi0 or eta.
void Print(const Option_t *opt) const
Print some relevant parameters set for the analysis.
void SetHistoYRangeAndNBins(Float_t min, Float_t max, Int_t n)
void SetTimeCut(Double_t min, Double_t max)
Definition: AliAnaPhoton.h:120
virtual void Print(const Option_t *opt) const
Print some relevant parameters set for the analysis.
void SetPtThresholdMax(Float_t pt)
virtual void SetZvertexCut(Float_t zcut=10.)
AliFiducialCut * GetFiducialCutForTrigger()
void SetEtaMassRange(Float_t min, Float_t max)
Definition: AliCaloPID.h:305
void SetNumberOfSuperModulesUsed(Int_t nSM)
UInt_t ConfigureAndGetEventTriggerMaskAndCaloTriggerString(TString trigger, Int_t year, TString &triggerString)
void SetPi0MassRange(Float_t min, Float_t max)
Definition: AliCaloPID.h:304
Class containing more common histogram axis types.
void SetEMCALLambda0CutMax(Float_t lcut)
Definition: AliCaloPID.h:202
void SwitchOnTrackMatchRejection()
Definition: AliAnaPi0EbE.h:148
virtual void SetNCentrBin(Int_t n=1)
Number of bins in reaction plain.
virtual void SetPtHardAndJetPtFactor(Float_t factor)
Class for event, clusters and tracks filtering and preparation for the ESD analysis.
virtual void SetEventPlaneMethod(TString m)
void SetPHOSEMin(Float_t en)
virtual void SetNRPBin(Int_t n=1)
Number of bins for vertex position.
void SetHistoNClustersRangeAndNBins(Int_t min, Int_t max, Int_t n)
void SetPHOSDispersionCut(Float_t dcut)
Definition: AliCaloPID.h:208
virtual void SetTrackCuts(AliESDtrackCuts *)
void SetNLMCut(Int_t min, Int_t max)
Definition: AliAnaPi0EbE.h:100
Track selection for correlation analysis.
AliAnaPhoton * ConfigurePhotonAnalysis(TString col, Bool_t simulation, TString calorimeter, Int_t year, Int_t tm, Bool_t printSettings, Int_t debug)
void SetMCGenType(Int_t min=0, Int_t max=6)
virtual void SetCentralityClass(TString name)
void SwitchOnFillShowerShapeHistograms()
Definition: AliAnaPhoton.h:73
Select clusters/tracks with low particle environment in their vecinity, isolated within a cone...
void SwitchOnBadTriggerEventsRemoval()
void SetHistoEtaRangeAndNBins(Float_t min, Float_t max, Int_t n)
void SwitchOnClusterELinearityCorrection()
void SwitchOffFillWeightHistograms()
Definition: AliAnaPi0EbE.h:152
void SetTrackTimeCut(Double_t a, Double_t b)
void SetEMCALTimeCut(Double_t a, Double_t b)
Configuration of the event trigger.
void SetHistoAsymmetryRangeAndNBins(Float_t min, Float_t max, Int_t n)
void SetHistoV0MultiplicityRangeAndNBins(Int_t min, Int_t max, Int_t n)
void SwitchOnEMCTrackPtDepResMatching()
Definition: AliCaloPID.h:228
virtual void SetAODObjArrayName(TString name)
Main class conecting the CaloTrackCorrelations package and Analysis Frame.
void SetNTriggerPtBins(Int_t n)
Set number of associated charged (neutral) hadrons pT bins.
void SwitchOffLoadOwnPHOSGeometryMatrices()
void SetHistoZRangeAndNBins(Float_t min, Float_t max, Int_t n)
void SetHistoV0SignalRangeAndNBins(Int_t min, Int_t max, Int_t n)
virtual void Print(const Option_t *) const
Print some relevant parameters set for the analysis.
Class that contains methods to select candidate cluster pairs to neutral meson.
void SetHistoTimeRangeAndNBins(Float_t min, Float_t max, Int_t n)
void SetHistoFinePtRangeAndNBins(Float_t min, Float_t max, Int_t n)
void SetEMCALEMax(Float_t en)
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)
void SetNAssocPtBins(Int_t n)
Set number of associated charged (neutral) hadrons pT bins.
Class for PID selection with calorimeters.
Definition: AliCaloPID.h:53
1 calorimeter cluster shower shape and split invariatn mass selection
Definition: AliAnaPi0EbE.h:87
void FillNOriginHistograms(Int_t n)
Definition: AliAnaPhoton.h:143
virtual void SetNMaxEvMix(Int_t n=20)
Number of bins in track multiplicity.
void SetNLMMinEnergy(Int_t i, Float_t min)
Definition: AliAnaPi0EbE.h:138
AliAnaGeneratorKine * ConfigureGenKineAnalysis(Int_t thresType, Float_t cone, Float_t pth, TString col, Bool_t simulation, TString calorimeter, Int_t year, Bool_t printSettings, Int_t debug)
void SetHistoShowerShapeRangeAndNBins(Float_t min, Float_t max, Int_t n)
bool Bool_t
Definition: External.C:53
void SetLeadHadronPhiCut(Float_t min, Float_t max)
void SetPHOSEMax(Float_t en)
void SetHistodEdxRangeAndNBins(Float_t min, Float_t max, Int_t n)
Class with utils specific to calorimeter clusters/cells.
virtual void SwitchOnConstrainTrackToVertex()
Int_t GetFirstSuperModuleUsed() const
void SetInvMassCutMaxParameters(Float_t a, Float_t b, Float_t c)
AliCalorimeterUtils * ConfigureCaloUtils(TString col, Bool_t simulation, Bool_t tender, TString calorimeter, Bool_t nonLinOn, Int_t year, Bool_t printSettings, Int_t debug)
void SetHistoHBPRangeAndNBins(Float_t min, Float_t max, Int_t n)
void SetEMCALLambda0CutMin(Float_t lcut)
Definition: AliCaloPID.h:205
void SetHistoRRangeAndNBins(Float_t min, Float_t max, Int_t n)
void SetParticleTypeInCone(Int_t i)
void SetHistoVertexDistRangeAndNBins(Float_t min, Float_t max, Int_t n)
void SetSubClusterEnergyMinimum(Int_t i, Float_t min)
Definition: AliCaloPID.h:273
void SetHistoDiffTimeRangeAndNBins(Float_t min, Float_t max, Int_t n)
void SetEMCALDEtaCut(Float_t dcut)
Definition: AliCaloPID.h:218
void SetPhotonMassRange(Float_t min, Float_t max)
Definition: AliCaloPID.h:306
void SwitchOffUseParametrizedTimeCut()
AliAnaPi0EbE * ConfigurePi0EbEAnalysis(TString particle, Int_t analysis, Bool_t useSSIso, Bool_t useAsy, TString col, Bool_t simulation, TString calorimeter, Int_t year, Int_t tm, Bool_t printSettings, Int_t debug)
void SetLocalMaximaCutE(Float_t cut)
void SetPtThresholds(Int_t i, Float_t pt)
Get year, collision type, mc/data type and period from alien global variables.
void SetHistoRatioRangeAndNBins(Float_t min, Float_t max, Int_t n)
void SwitchOffAllNLMHistoFill()
Definition: AliAnaPi0EbE.h:106
void SetNLMCut(Int_t min, Int_t max)
Definition: AliAnaPhoton.h:128
Int_t GetFirstSMCoveredByTRD() const
Time cut in ns.
Steering class of package CaloTrackCorrelartions.
void SetTriggerPtRange(Float_t min, Float_t max)