1 - Rename AliAnaMaker and AliAnaBaseClass to AliAnaPartCorrMaker and AliAnaPartCorrB...
[u/mrichter/AliRoot.git] / PWG4 / macros / ConfigAnalysisGammaHadronCorrelation.C
1 /* $Id: $ */
2 /* $Log$ */
3
4 //------------------------------------
5 // Configuration macro example:
6 //
7 // Do prompt photon - hadron correlation analysis with ESDs
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   //Detector Fidutial Cuts
23   AliFidutialCut * fidCut = new AliFidutialCut();
24   fidCut->DoCTSFidutialCut(kTRUE) ;
25   fidCut->DoEMCALFidutialCut(kTRUE) ;
26   fidCut->DoPHOSFidutialCut(kTRUE) ;
27   
28   fidCut->SetSimpleCTSFidutialCut(0.9,0.,360.);
29   fidCut->SetSimpleEMCALFidutialCut(0.7,80.,190.);
30   fidCut->SetSimplePHOSFidutialCut(0.13,220.,320.);
31  
32   fidCut->Print("");
33   
34   //-----------------------------------------------------------  
35   // Reader
36   //-----------------------------------------------------------
37   AliCaloTrackReader *reader = new AliCaloTrackESDReader();
38   reader->SetDebug(-1);
39
40   //Switch on or off the detectors information that you want
41   reader->SwitchOnEMCAL();
42   reader->SwitchOnCTS();
43   reader->SwitchOnPHOS();
44   reader->SwitchOffEMCALCells();
45   reader->SwitchOffPHOSCells();
46  
47
48   //Min particle pT
49   reader->SetEMCALPtMin(0.5); 
50   reader->SetPHOSPtMin(0.5);
51   reader->SetCTSPtMin(0.2);
52   
53   reader->SetFidutialCut(fidCut);
54   reader->Print("");
55   
56   
57   //---------------------------------------------------------------------
58   // Analysis algorithm
59   //---------------------------------------------------------------------
60   
61   //Detector Fidutial Cuts for analysis part
62   AliFidutialCut * fidCut2 = new AliFidutialCut();
63   fidCut2->DoCTSFidutialCut(kFALSE) ;
64   fidCut2->DoEMCALFidutialCut(kTRUE) ;
65   fidCut2->DoPHOSFidutialCut(kFALSE) ;
66   
67   fidCut2->SetSimpleCTSFidutialCut(0.9,0.,360.);
68   fidCut2->SetSimpleEMCALFidutialCut(0.7,80.,190.);
69   fidCut2->SetSimplePHOSFidutialCut(0.13,220.,320.);
70   fidCut2->Print("");
71
72   AliCaloPID * pid = new AliCaloPID();
73   // use selection with simple weights
74   pid->SetPHOSPhotonWeight(0.7);    pid->SetPHOSPi0Weight(0.7); 
75   pid->SetEMCALPhotonWeight(0.7);    pid->SetEMCALPi0Weight(0.7);
76   
77   pid->Print("");
78
79   AliIsolationCut * ic = new AliIsolationCut();
80   ic->SetConeSize(0.4);
81   ic->SetPtThreshold(1.);
82   ic->SetICMethod(AliIsolationCut::kPtThresIC);
83   ic->Print("");
84
85   //Photon Analysis
86   AliAnaGammaDirect *ana = new AliAnaGammaDirect();
87   ana->SetDebug(-1);
88   ana->SetMinPt(5.);
89   ana->SetCaloPID(pid);
90   ana->SetFidutialCut(fidCut2);
91   ana->SetIsolationCut(ic) ;
92   ana->SetDetector("PHOS");
93   ana->SwitchOnIsolation();
94   ana->SwitchOnCaloPID();
95   ana->SwitchOffCaloPIDRecalculation(); //recommended for EMCAL
96   ana->SwitchOffFidutialCut();
97   ana->SwitchOffInvariantMass();
98   ana->SwitchOffDataMC() ;
99
100   ana->Print("");
101
102   //Photon- hadron correlation
103   AliAnaParticleHadronCorrelation *ana2 = new AliAnaParticleHadronCorrelation();
104   ana2->SetDebug(-1);
105   ana2->SetCaloPID(pid);
106   ana2->SwitchOnCaloPID();
107   ana2->SwitchOnCaloPIDRecalculation(); //recommended for EMCAL
108   ana2->SwitchOffFidutialCut();
109   ana2->SetPtCutRange(1,100);
110   ana2->SetDeltaPhiCutRange(1.5,4.5);
111   
112   ana2->Print("");
113   //---------------------------------------------------------------------
114   // Set  analysis algorithm and reader
115   //---------------------------------------------------------------------
116   maker = new AliAnaPartCorrMaker();
117   maker->SetReader(reader);//pointer to reader
118   maker->AddAnalysis(ana,0);
119   maker->AddAnalysis(ana2,1);
120   maker->SetAODBranchName("PhotonHadronCorrelation");
121   maker->SetAnaDebug(-1)  ;
122   maker->SwitchOnHistogramsMaker()  ;
123   //maker->SwitchOffHistogramsMaker() ;  
124   maker->SwitchOnAODsMaker()  ;
125   //maker->SwitchOffAODsMaker() ; 
126   
127   maker->Print("");
128   //
129   printf("======================== \n");
130   printf("END ConfigAnalysis() \n");
131   printf("======================== \n");
132   return maker ;
133 }
134