PWG4/macros/electrons/ConfigAnalysisElectron.C

   1 /* $Id: $ */\r
   2 /* $Log$ */\r
   3 \r
   4 //------------------------------------\r
   5 // Configuration macro example:\r
   6 //\r
   7 // Do EMCal electron analysis with ESDs\r
   8 //\r
   9 //------------------------------------\r
  10 \r
  11 AliAnaPartCorrMaker*  ConfigAnalysis()\r
  12 {\r
  13         //\r
  14         // Configuration goes here\r
  15         // \r
  16         printf("======================== \n");\r
  17         printf("ConfigAnalysisElectron() \n");\r
  18         printf("======================== \n");\r
  19         \r
  20         \r
  21         //Detector Fidutial Cuts\r
  22         //AliFidutialCut * fidCut = new AliFidutialCut();\r
  23         //fidCut->DoEMCALFidutialCut(kTRUE) ;\r
  24         \r
  25         //fidCut->SetSimpleEMCALFidutialCut(0.7,80.,190.);\r
  26         \r
  27         //fidCut->Print("");\r
  28         \r
  29         \r
  30         //-----------------------------------------------------------  \r
  31         // Reader\r
  32         //-----------------------------------------------------------\r
  33         AliCaloTrackESDReader *reader = new AliCaloTrackESDReader();\r
  34         reader->SetDebug(10);//10 for lots of messages\r
  35         \r
  36         //Switch on or off the detectors information that you want\r
  37         reader->SwitchOffCTS();\r
  38         reader->SwitchOffEMCAL();\r
  39         reader->SwitchOffEMCALCells();  \r
  40         reader->SwitchOffPHOS();\r
  41         reader->SwitchOffPHOSCells();   \r
  42         //Min particle pT\r
  43         //reader->SetEMCALPtMin(1.); \r
  44         \r
  45         //reader->SetFidutialCut(fidCut);\r
  46         reader->Print("");\r
  47         \r
  48         \r
  49         //---------------------------------------------------------------------\r
  50         // Analysis algorithm\r
  51         //---------------------------------------------------------------------\r
  52         \r
  53         AliAnaElectron *anaelectron = new AliAnaElectron();\r
  54         anaelectron->SetDebug(10); //10 for lots of messages\r
  55         anaelectron->SetCalorimeter("EMCAL");\r
  56         anaelectron->SetpOverEmin(0.8);\r
  57         anaelectron->SetpOverEmax(1.2);\r
  58         anaelectron->SetResidualCut(0.05);\r
  59         anaelectron->SwitchOnDataMC();\r
  60         anaelectron->SetMinPt(1.);\r
  61         //Set Histrograms bins and ranges\r
  62         anaelectron->SetHistoPtRangeAndNBins(0, 50, 100) ;\r
  63         anaelectron->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 100) ;\r
  64         anaelectron->SetHistoEtaRangeAndNBins(-0.7, 0.7, 100) ; anaelectron->Print("");\r
  65         \r
  66         //Detector Fidutial Cuts\r
  67         AliFidutialCut * fidCut2 = new AliFidutialCut();\r
  68         fidCut2->DoEMCALFidutialCut(kTRUE) ;\r
  69         fidCut2->SetSimpleEMCALFidutialCut(0.7,80.,190.);\r
  70         fidCut2->Print("");\r
  71                 \r
  72         \r
  73         //---------------------------------------------------------------------\r
  74         // Set  analysis algorithm and reader\r
  75         //---------------------------------------------------------------------\r
  76         maker = new AliAnaPartCorrMaker();\r
  77         maker->SetReader(reader);//pointer to reader\r
  78         maker->AddAnalysis(anaelectron,0);\r
  79         maker->SetAnaDebug(10)  ;\r
  80         maker->SwitchOnHistogramsMaker()  ;\r
  81         //maker->SwitchOnAODsMaker()  ;\r
  82         \r
  83         maker->Print("");\r
  84         //\r
  85         printf("======================== \n");\r
  86         printf("END ConfigAnalysisElectron() \n");\r
  87         printf("======================== \n");\r
  88         return maker ;\r
  89 }\r