Corrected analysis configuration and execution files after changes in PartCorr frame...
[u/mrichter/AliRoot.git] / PWG4 / macros / ConfigAnalysisGammaHadronCorrelation.C
CommitLineData
d92b41ad 1/* $Id: $ */
2/* $Log$ */
3
4//------------------------------------
5// Configuration macro example:
6//
7// Do prompt photon - hadron correlation analysis with ESDs
3e0577a2 8// First find photons with AliAnaPhoton, then
9// isolate them with AliAnaParticleIsolation and finally correlate
10// them with AliAnaParticleHadronCorrelation
d92b41ad 11//
12// Author : Gustavo Conesa Balbastre (INFN-LNF)
13//------------------------------------
14
c90ac396 15AliAnaPartCorrMaker* ConfigAnalysis()
d92b41ad 16{
3e0577a2 17 //
18 // Configuration goes here
19 //
20 printf("======================== \n");
21 printf("ConfigAnalysis() \n");
22 printf("======================== \n");
23
24
25 //Detector Fidutial Cuts
26 AliFidutialCut * fidCut = new AliFidutialCut();
27 fidCut->DoCTSFidutialCut(kTRUE) ;
28 fidCut->DoEMCALFidutialCut(kTRUE) ;
29 fidCut->DoPHOSFidutialCut(kTRUE) ;
30
31 fidCut->SetSimpleCTSFidutialCut(0.9,0.,360.);
32 fidCut->SetSimpleEMCALFidutialCut(0.7,80.,190.);
33 fidCut->SetSimplePHOSFidutialCut(0.13,220.,320.);
34
35 fidCut->Print("");
36
37 //-----------------------------------------------------------
38 // Reader
39 //-----------------------------------------------------------
40 AliCaloTrackESDReader *reader = new AliCaloTrackESDReader();
41 reader->SetDebug(-1);
42
43 //Switch on or off the detectors information that you want
44 reader->SwitchOnEMCAL();
45 reader->SwitchOnCTS();
46 reader->SwitchOnPHOS();
47 reader->SwitchOffEMCALCells();
48 reader->SwitchOffPHOSCells();
49
50
51 //Min particle pT
52 reader->SetEMCALPtMin(0.5);
53 reader->SetPHOSPtMin(0.5);
54 reader->SetCTSPtMin(0.2);
55
56 reader->SetFidutialCut(fidCut);
57 reader->Print("");
58
59
60 //---------------------------------------------------------------------
61 // Analysis algorithm
62 //---------------------------------------------------------------------
63 //<<<< first analysis >>> select the photons
64 //Detector Fidutial Cuts for analysis part
65 AliFidutialCut * fidCut2 = new AliFidutialCut();
66 fidCut2->DoCTSFidutialCut(kFALSE) ;
67 fidCut2->DoEMCALFidutialCut(kTRUE) ;
68 fidCut2->DoPHOSFidutialCut(kFALSE) ;
69
70 fidCut2->SetSimpleCTSFidutialCut(0.9,0.,360.);
71 fidCut2->SetSimpleEMCALFidutialCut(0.7,80.,190.);
72 fidCut2->SetSimplePHOSFidutialCut(0.13,220.,320.);
73 fidCut2->Print("");
74
75 AliCaloPID * pid = new AliCaloPID();
76 // use selection with simple weights
77 pid->SetPHOSPhotonWeight(0.7); pid->SetPHOSPi0Weight(0.7);
78 pid->SetEMCALPhotonWeight(0.7); pid->SetEMCALPi0Weight(0.7);
79
80 pid->Print("");
81
82
83 AliAnaPhoton *anaphoton = new AliAnaPhoton();
84 anaphoton->SetDebug(-1); //10 for lots of messages
85 anaphoton->SetMinPt(0.2);
86 anaphoton->SetMinDistanceToBadChannel(2, 4, 5);
87 anaphoton->SetCaloPID(pid);
88 anaphoton->SetFidutialCut(fidCut2); //More acceptance selections if needed at this level
89 anaphoton->SetCalorimeter("PHOS");
90 anaphoton->SwitchOffDataMC() ;//Access MC stack and fill more histograms
91 anaphoton->SwitchOnCaloPID();
92 anaphoton->SwitchOffCaloPIDRecalculation(); //recommended for EMCAL
93 anaphoton->SwitchOnFidutialCut();
94 anaphoton->SetOutputAODName("Photons");
95 anaphoton->SetOutputAODClassName("AliAODPWG4ParticleCorrelation");
96 //Set Histrograms bins and ranges
97 // anaphoton->SetHistoPtRangeAndNBins(0, 50, 100) ;
98 // anaphoton->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 100) ;
99 // anaphoton->SetHistoEtaRangeAndNBins(-0.7, 0.7, 100) ;
100 anaphoton->Print("");
101
102 // >>>> Second Analysis <<<< Isolate the photons
103 AliIsolationCut * ic = new AliIsolationCut();
104 ic->SetConeSize(0.5);
105 ic->SetPtThreshold(1.);
106 ic->SetICMethod(AliIsolationCut::kPtThresIC);
107 ic->Print("");
108
109 AliAnaParticleIsolation *anaisol = new AliAnaParticleIsolation();
110 anaisol->SetDebug(-1);
111 anaisol->SetMinPt(5);
112 anaisol->SetInputAODName("Photons");
113 anaisol->SetCalorimeter("PHOS");
114 anaisol->SwitchOffDataMC() ;//Access MC stack and fill more histograms
115 //Select clusters with no pair, if both clusters with pi0 mass
116 anaisol->SwitchOffInvariantMass();
117 //anaisol->SetNeutralMesonSelection(nms);
118 //Do isolation cut
119 anaisol->SetIsolationCut(ic);
120 //Do or not do isolation with previously produced AODs.
121 //No effect if use of SwitchOnSeveralIsolation()
122 anaisol->SwitchOffReIsolation();
123 //Multiple IC
124 anaisol->SwitchOffSeveralIsolation() ;
125
126 anaisol->Print("");
127
128 //<<<Third analysis>>> Isolated Photon- hadron correlation
129 AliAnaParticleHadronCorrelation *anacorr = new AliAnaParticleHadronCorrelation();
130 anacorr->SetInputAODName("Photons");
131 anacorr->SetDebug(-1);
132 anacorr->SetCaloPID(pid);
133 anacorr->SwitchOnCaloPID();
134 anacorr->SwitchOnCaloPIDRecalculation(); //recommended for EMCAL
135 anacorr->SwitchOffFidutialCut();
136 anacorr->SetPtCutRange(1,100);
137 anacorr->SetDeltaPhiCutRange(1.5,4.5);
138 //Set Histrograms bins and ranges
139 // anaphoton->SetHistoPtRangeAndNBins(0, 50, 100) ;
140 // anaphoton->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 100) ;
141 // anaphoton->SetHistoEtaRangeAndNBins(-0.7, 0.7, 100) ;
142 anacorr->Print("");
143
144 //---------------------------------------------------------------------
145 // Set analysis algorithm and reader
146 //---------------------------------------------------------------------
147 maker = new AliAnaPartCorrMaker();
148 maker->SetReader(reader);//pointer to reader
149 maker->AddAnalysis(anaphoton,0);
150 maker->AddAnalysis(anaisol,1);
151 maker->AddAnalysis(anacorr,2);
152 maker->SetAnaDebug(-1) ;
153 maker->SwitchOnHistogramsMaker() ;
154 maker->SwitchOnAODsMaker() ;
155
156 maker->Print("");
157 //
158 printf("======================== \n");
159 printf("END ConfigAnalysis() \n");
160 printf("======================== \n");
161 return maker ;
d92b41ad 162}
163