PWG4/macros/ConfigAnalysisCalorimeterQA.C

   1 /* $Id: $ */
   2
   3 //------------------------------------
   4 // Configuration macro example:
   5 //
   6 // Calorimeters QA: Validation of data (MC)
   7 // Valid for ESDs, comment/uncomment below  in the reader part for AODs
   8 //
   9 // Author : Gustavo Conesa Balbastre (INFN-LNF)
  10 //------------------------------------
  11
  12 AliAnaPartCorrMaker*  ConfigAnalysis()
  13 {
  14         //
  15         // Configuration goes here
  16         //
  17         printf("======================== \n");
  18         printf("ConfigAnalysis() \n");
  19         printf("======================== \n");
  20
  21
  22
  23         //-----------------------------------------------------------
  24         // Reader
  25         //-----------------------------------------------------------
  26         //For this particular analysis few things done by the reader.
  27         //Nothing else needs to be set.
  28         AliCaloTrackESDReader *reader = new AliCaloTrackESDReader();
  29         reader->SetDebug(-1);
  30         reader->SwitchOnStack();
  31         reader->SwitchOffAODMCParticles();
  32         reader->SetDeltaAODFileName(""); //Do not create deltaAOD file, this analysis do not create branches.
  33         reader->Print("");
  34
  35         //For AODs:
  36 //      AliCaloTrackAODReader *reader = new AliCaloTrackAODReader();
  37 //      reader->SetDebug(-1);
  38 //      reader->SwitchOffStack();
  39 //      reader->SwitchOnAODMCParticles();
  40 //      reader->SetDeltaAODFileName(""); //Do not create deltaAOD file, this analysis do not create branches.
  41 //      reader->Print("");
  42
  43
  44         //---------------------------------------------------------------------
  45         // Analysis algorithm
  46         //---------------------------------------------------------------------
  47
  48         AliFiducialCut * fidCut = new AliFiducialCut();
  49         fidCut->DoCTSFiducialCut(kFALSE) ;
  50         fidCut->DoEMCALFiducialCut(kTRUE) ;
  51         fidCut->DoPHOSFiducialCut(kTRUE) ;
  52
  53
  54         AliAnaCalorimeterQA *anaEMCAL = new AliAnaCalorimeterQA();
  55         anaEMCAL->SetDebug(-1); //10 for lots of messages
  56         anaEMCAL->SetCalorimeter("EMCAL");
  57         anaEMCAL->SwitchOnDataMC() ;//Access MC stack and fill more histograms
  58         //anaEMCAL->SetStyleMacro("$MACROS/style.C") ;
  59         anaEMCAL->AddToHistogramsName("AnaCaloQA_EMCAL_");
  60         anaEMCAL->SetFiducialCut(fidCut);
  61         anaEMCAL->SwitchOnFiducialCut();
  62         anaEMCAL->Print("");
  63
  64         AliAnaCalorimeterQA *anaPHOS = new AliAnaCalorimeterQA();
  65         anaPHOS->SetDebug(-1); //10 for lots of messages
  66         anaPHOS->SetCalorimeter("PHOS");
  67         anaPHOS->SwitchOnDataMC() ;//Access MC stack and fill more histograms
  68         //anaPHOS->SetStyleMacro("$MACROS/style.C") ;
  69         anaPHOS->AddToHistogramsName("AnaCaloQA_PHOS_");
  70         anaPHOS->SetFiducialCut(fidCut);
  71         anaPHOS->SwitchOnFiducialCut();
  72         anaPHOS->Print("");
  73
  74
  75         //---------------------------------------------------------------------
  76         // Set  analysis algorithm and reader
  77         //---------------------------------------------------------------------
  78         maker = new AliAnaPartCorrMaker();
  79         maker->SetReader(reader);//pointer to reader
  80         maker->AddAnalysis(anaEMCAL,1);
  81         maker->AddAnalysis(anaPHOS,0);
  82         //maker->SetAnaDebug(0)  ;
  83         maker->SwitchOnHistogramsMaker()  ;
  84         maker->SwitchOffAODsMaker()  ;//No AODs created in this task.
  85
  86         maker->Print("");
  87         //
  88         printf("======================== \n");
  89         printf("END ConfigAnalysis() \n");
  90         printf("======================== \n");
  91         return maker ;
  92 }
  93