Added new background scheeme, did some cleanup. added new bethe block parameters...
[u/mrichter/AliRoot.git] / 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