1 /**************************************************************************
2 * Copyright(c) 1998-2008, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
16 /////////////////////////////////////////////////////////////
18 // Class to analyze ZDC data //
20 /////////////////////////////////////////////////////////////
28 #include "AliAnalysisManager.h"
29 #include "AliInputEventHandler.h"
30 #include "AliVEvent.h"
32 #include "AliESDEvent.h"
33 #include "AliESDHeader.h"
34 #include "AliESDInputHandler.h"
35 #include "AliESDZDC.h"
36 #include "AliMultiplicity.h"
37 #include "AliAODHandler.h"
38 #include "AliAODEvent.h"
39 #include "AliAODVertex.h"
40 #include "AliAODMCHeader.h"
41 #include "AliMCEventHandler.h"
42 #include "AliMCEvent.h"
43 #include "AliHeader.h"
44 #include "AliAODMCParticle.h"
45 #include "AliAnalysisTaskSE.h"
46 #include "AliGenEventHeader.h"
47 #include "AliGenHijingEventHeader.h"
48 #include "AliPhysicsSelectionTask.h"
49 #include "AliPhysicsSelection.h"
50 #include "AliBackgroundSelection.h"
51 #include "AliTriggerAnalysis.h"
52 #include "AliCentrality.h"
53 #include "AliAnalysisTaskZDCpA.h"
55 ClassImp(AliAnalysisTaskZDCpA)
58 //________________________________________________________________________
59 AliAnalysisTaskZDCpA::AliAnalysisTaskZDCpA():
84 // Default constructor
87 //________________________________________________________________________
88 AliAnalysisTaskZDCpA::AliAnalysisTaskZDCpA(const char *name):
89 AliAnalysisTaskSE(name),
113 // Output slot #1 writes into a TList container
114 DefineOutput(1, TList::Class());
118 //________________________________________________________________________
119 AliAnalysisTaskZDCpA& AliAnalysisTaskZDCpA::operator=(const AliAnalysisTaskZDCpA& c)
122 // Assignment operator
125 AliAnalysisTaskSE::operator=(c);
130 //________________________________________________________________________
131 AliAnalysisTaskZDCpA::AliAnalysisTaskZDCpA(const AliAnalysisTaskZDCpA& ana):
132 AliAnalysisTaskSE(ana),
134 fIsMCInput(ana.fIsMCInput),
135 fOutput(ana.fOutput),
136 fhTDCZNSum(ana.fhTDCZNSum),
137 fhTDCZNDiff(ana.fhTDCZNDiff),
138 fhZNCSpectrum(ana.fhZNCSpectrum),
139 fhZNASpectrum(ana.fhZNASpectrum),
140 fhZPCSpectrum(ana.fhZPCSpectrum),
141 fhZPASpectrum(ana.fhZPASpectrum),
142 fhZEM1Spectrum(ana.fhZEM1Spectrum),
143 fhZEM2Spectrum(ana.fhZEM2Spectrum),
144 fhZNCpmc(ana.fhZNCpmc),
145 fhZNApmc(ana.fhZNApmc),
146 fhZPCpmc(ana.fhZPCpmc),
147 fhZPApmc(ana.fhZPApmc),
148 fhZNCCentroid(ana.fhZNCCentroid),
149 fhZNACentroid(ana.fhZNACentroid),
150 fhZNCemd(ana.fhZNCemd),
151 fhZNAemd(ana.fhZNAemd),
152 fhPMCZNCemd(ana.fhPMCZNCemd),
153 fhPMCZNAemd(ana.fhPMCZNAemd),
154 fDebunch(ana.fDebunch)
161 //________________________________________________________________________
162 AliAnalysisTaskZDCpA::~AliAnalysisTaskZDCpA()
165 if(fOutput && !AliAnalysisManager::GetAnalysisManager()->IsProofMode()){
166 delete fOutput; fOutput=0;
171 //________________________________________________________________________
172 void AliAnalysisTaskZDCpA::UserCreateOutputObjects()
174 // Create the output containers
178 //fOutput->SetName("output");
180 fhTDCZNSum = new TH1F("fhTDCZNSum","TDC_{ZNC}+TDC_{ZNA}",60,-30.,30.);
181 fhTDCZNSum->GetXaxis()->SetTitle("TDC_{ZNC}+TDC_{ZNA} (ns)");
182 fOutput->Add(fhTDCZNSum);
184 fhTDCZNDiff = new TH1F("fhTDCZNDiff","TDC_{ZNC}-TDC_{ZNA}",60,-30.,30.);
185 fhTDCZNDiff->GetXaxis()->SetTitle("TDC_{ZNC}-TDC_{ZNA} (ns)");
186 fOutput->Add(fhTDCZNDiff);
188 fhZNCSpectrum = new TH1F("fhZNCSpectrum", "ZNC signal", 250,0., 2500.);
189 fOutput->Add(fhZNCSpectrum);
190 fhZNASpectrum = new TH1F("fhZNASpectrum", "ZNA signal", 250,0., 2500.) ;
191 fOutput->Add(fhZNASpectrum);
192 fhZPCSpectrum = new TH1F("fhZPCSpectrum", "ZPC signal", 200,0., 2000.) ;
193 fOutput->Add(fhZPCSpectrum);
194 fhZPASpectrum = new TH1F("fhZPASpectrum", "ZPA signal", 200,0., 2000.) ;
195 fOutput->Add(fhZPASpectrum);
196 fhZEM1Spectrum = new TH1F("fhZEM1Spectrum", "ZEM1 signal", 100,0., 2500.);
197 fOutput->Add(fhZEM1Spectrum);
198 fhZEM2Spectrum = new TH1F("fhZEM2Spectrum", "ZEM2 signal", 100,0., 2500.);
199 fOutput->Add(fhZEM2Spectrum);
201 fhZNCpmc = new TH1F("fhZNCpmc","ZNC PMC",200, 0., 2000.);
202 fOutput->Add(fhZNCpmc);
203 fhZNApmc = new TH1F("fhZNApmc","ZNA PMC",200, 0., 12000.);
204 fOutput->Add(fhZNApmc);
205 fhZPCpmc = new TH1F("fhZPCpmc","ZPC PMC",200, 0., 2000.);
206 fOutput->Add(fhZPCpmc);
207 fhZPApmc = new TH1F("fhZPApmc","ZPA PMC",200, 0., 2000.);
208 fOutput->Add(fhZPApmc);
210 fhZNCCentroid = new TH2F("fhZNCCentroid","Centroid over ZNC",70,-3.5,3.5,70,-3.5,3.5);
211 fOutput->Add(fhZNCCentroid);
212 fhZNACentroid = new TH2F("fhZNACentroid","Centroid over ZNA",70,-3.5,3.5,70,-3.5,3.5);
213 fOutput->Add(fhZNACentroid);
215 fhZNCemd = new TH1F("fhZNCemd","ZNC signal lg",200,0.,6000.);
216 fOutput->Add(fhZNCemd);
217 fhZNAemd = new TH1F("fhZNAemd","ZNA signal lg",200,0.,6000.);
218 fOutput->Add(fhZNAemd);
219 fhPMCZNCemd = new TH1F("fhPMCZNCemd","ZNC PMC lg",200, 10., 6000.);
220 fOutput->Add(fhPMCZNCemd);
221 fhPMCZNAemd = new TH1F("fhPMCZNAemd","ZNA PMC lg",200, 10., 6000.);
222 fOutput->Add(fhPMCZNAemd);
224 fDebunch = new TH2F("fDebunch","ZN TDC sum vs. diff", 120,-30,30,120,-30,-30);
225 fOutput->Add(fDebunch);
227 PostData(1, fOutput);
230 //________________________________________________________________________
231 void AliAnalysisTaskZDCpA::UserExec(Option_t */*option*/)
233 // Execute analysis for current event:
234 if(fDebug>1) printf(" **** AliAnalysisTaskZDCpA::UserExec() \n");
237 Printf("ERROR: InputEvent not available");
242 AliESDEvent* esd = dynamic_cast<AliESDEvent*> (InputEvent());
244 // Select PHYSICS events (type=7, for data)
245 //if(!fIsMCInput && esd->GetEventType()!=7) return;
247 // ********* MC INFO *********************************
250 AliMCEventHandler* eventHandler = dynamic_cast<AliMCEventHandler*> (AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler());
252 Printf("ERROR: Could not retrieve MC event handler");
256 AliMCEvent* mcEvent = eventHandler->MCEvent();
258 Printf("ERROR: Could not retrieve MC event");
263 // ****************************************************
265 AliAnalysisManager *am = AliAnalysisManager::GetAnalysisManager();
267 AliESDZDC *esdZDC = esd->GetESDZDC();
269 // if((((AliInputEventHandler*)(am->GetInputEventHandler()))->IsEventSelected())){
271 fhZNCSpectrum->Fill(esdZDC->GetZDCN1Energy());
272 fhZNASpectrum->Fill(esdZDC->GetZDCN2Energy());
273 fhZPCSpectrum->Fill(esdZDC->GetZDCP1Energy());
274 fhZPASpectrum->Fill(esdZDC->GetZDCP2Energy());
275 fhZEM1Spectrum->Fill(esdZDC->GetZDCEMEnergy(0));
276 fhZEM2Spectrum->Fill(esdZDC->GetZDCEMEnergy(1));
278 const Double_t * towZNC = esdZDC->GetZN1TowerEnergy();
279 const Double_t * towZPC = esdZDC->GetZP1TowerEnergy();
280 const Double_t * towZNA = esdZDC->GetZN2TowerEnergy();
281 const Double_t * towZPA = esdZDC->GetZP2TowerEnergy();
283 fhZNCpmc->Fill(towZNC[0]);
284 fhZNApmc->Fill(towZNA[0]);
285 fhZPCpmc->Fill(towZPC[0]);
286 fhZPApmc->Fill(towZPA[0]);
288 Double_t xyZNC[2]={-99.,-99.}, xyZNA[2]={-99.,-99.};
289 esdZDC->GetZNCentroidInpp(xyZNC, xyZNA);
291 fhZNCCentroid->Fill(xyZNC[0], xyZNC[1]);
292 fhZNACentroid->Fill(xyZNA[0], xyZNA[1]);
294 const Double_t * towZNCLG = esdZDC->GetZN1TowerEnergyLR();
295 const Double_t * towZNALG = esdZDC->GetZN2TowerEnergyLR();
296 Double_t znclg=0., znalg=0.;
297 for(Int_t iq=0; iq<5; iq++){
298 znclg += towZNCLG[iq];
299 znalg += towZNALG[iq];
301 fhZNCemd->Fill(znclg/2.);
302 fhZNAemd->Fill(znalg/2.);
303 fhPMCZNCemd->Fill(towZNCLG[0]);
304 fhPMCZNAemd->Fill(towZNALG[0]);
308 Float_t tdcC=999., tdcA=999;
309 Float_t tdcSum=999., tdcDiff=999;
310 for(int i=0; i<4; i++){
311 if(esdZDC->GetZDCTDCData(10,i) != 0.){
312 tdcC = esdZDC->GetZDCTDCCorrected(10,i);
313 if(esdZDC->GetZDCTDCData(12,i) != 0.){
314 tdcA = esdZDC->GetZDCTDCCorrected(12,i);
319 if(tdcSum!=999.) fhTDCZNSum->Fill(tdcSum);
320 if(tdcDiff!=999.)fhTDCZNDiff->Fill(tdcDiff);
321 if(tdcSum!=999. && tdcDiff!=999.) fDebunch->Fill(tdcDiff, tdcSum);
324 PostData(1, fOutput);
330 //________________________________________________________________________
331 void AliAnalysisTaskZDCpA::Terminate(Option_t */*option*/)
333 // Terminate analysis