9 #include "AliAnalysisTask.h"
10 #include "AliAnalysisManager.h"
12 #include "AliESDEvent.h"
13 #include "AliESDInputHandler.h"
15 #include "AliT0AnalysisTaskQA.h"
17 //#include "AliCDBMetaData.h"
18 //#include "AliCDBId.h"
19 //#include "AliCDBEntry.h"
20 //#include "AliCDBManager.h"
21 //#include "AliCDBStorage.h"
23 // Task should calculate channels offset
26 ClassImp(AliT0AnalysisTaskQA)
27 //________________________________________________________________________
28 AliT0AnalysisTaskQA::AliT0AnalysisTaskQA()
29 : AliAnalysisTaskSE(), fESD(0x0), fTzeroObject(0x0),
30 fTzeroORA(0x0), fTzeroORC(0x0), fResolution(0x0), fTzeroORAplusORC(0x0),
31 fRunNumber(0),fTimeVSAmplitude(0x0),fCFDVSPmtId(0x0),fSPDVertexVST0Vertex(0x0),
32 fOrAvsNtracks(0), fOrCvsNtracks(0), fT0vsNtracks(0),
33 fEffAC(0), fEffA(0), fEffC(0), ftracksEffSPD(0)
38 // Define input and output slots here
39 // Input slot #0 works with a TChain
41 // ########### NEVER define slots in the IO constructor
42 // DefineInput(0, TChain::Class());
43 // DefineOutput(1, TObjArray::Class());
47 //________________________________________________________________________
48 AliT0AnalysisTaskQA::AliT0AnalysisTaskQA(const char *name)
49 : AliAnalysisTaskSE(name), fESD(0x0), fTzeroObject(0x0),
50 fTzeroORA(0x0), fTzeroORC(0x0), fResolution(0x0), fTzeroORAplusORC(0x0),
51 fRunNumber(0),fTimeVSAmplitude(0x0),fCFDVSPmtId(0x0),fSPDVertexVST0Vertex(0x0),
52 fOrAvsNtracks(0), fOrCvsNtracks(0), fT0vsNtracks(0),
53 fEffAC(0), fEffA(0), fEffC(0), ftracksEffSPD(0)
57 // Define input and output slots here
58 // Input slot #0 works with a TChain
59 DefineInput(0, TChain::Class());
60 DefineOutput(1, TObjArray::Class());
61 // Output slot #0 id reserved by the base class for AOD
62 // Output slot #1 writes into a TH1 container
65 //________________________________________________________________________
66 AliT0AnalysisTaskQA::~AliT0AnalysisTaskQA()
69 // printf("AliT0CalibOffsetChannels~AliT0CalibOffsetChannels() ");
70 delete fTzeroObject; // This deletes everything ...
73 //------------------------------------------------------------------
74 void AliT0AnalysisTaskQA::UserCreateOutputObjects()
77 fTimeVSAmplitude = new TH2F*[NPMT0];
79 for (Int_t i=0; i<NPMT0; i++) {
80 fTimeVSAmplitude[i]= new TH2F (Form("fTimeVSAmplitude%d",i+1),"fTimeVsAmplitude",500, 0, 50,1500,1000,7000);
83 fTzeroORAplusORC = new TH1F("fTzeroORAplusORC","ORA+ORC /2",100,-2000,2000); //or A plus or C
84 fResolution = new TH1F("fResolution","fResolution",100,-2000,2000);// or A minus or C spectrum
85 fTzeroORA = new TH1F("fTzeroORA","fTzeroORA",100,-2000,2000);// or A spectrum
86 fTzeroORC = new TH1F("fTzeroORC","fTzeroORC",100,-2000,2000);// or C spectrum
87 fCFDVSPmtId = new TH2F("fCFDVSPmtId","fCFDVSPmtId",24,0,24,1500,1000,7000); //
88 fSPDVertexVST0Vertex = new TH2F("fSPDVertexVST0Vertex","fSPDVertexVST0Vertex",30,-30,30,30,-30,30);
89 fOrAvsNtracks = new TH2F("fAvstracks", "Avstracks",200, 0, 1000, 500, -1000, 1000);
90 fOrCvsNtracks = new TH2F("fCvstracks", "Cvstracks",200, 0, 1000, 500, -1000, 1000);
91 fT0vsNtracks = new TH2F("fT0ACvstrackles", "T0ACvstracks",200, 0, 1000, 500, -1000, 1000);
92 fEffAC = new TH1F("EffAC","T0Eff ", 200, 0, 1000);
93 fEffA = new TH1F("EffA","T0AEff ", 200, 0, 1000);
94 fEffC = new TH1F("EffC","T0CEff ", 200, 0, 1000);
95 ftracksEffSPD= new TH1F("ftracksEffSPD","SPDeff", 200, 0, 1000);
97 fTzeroObject = new TObjArray();
98 fTzeroObject->SetOwner(kTRUE);
100 for (Int_t i=0; i<24; i++)
101 fTzeroObject->AddAtAndExpand(fTimeVSAmplitude[i],i);
103 fTzeroObject->AddAtAndExpand(fCFDVSPmtId,24);
104 fTzeroObject->AddAtAndExpand(fSPDVertexVST0Vertex,25);
105 fTzeroObject->AddAtAndExpand(fTzeroORAplusORC, 26);
106 fTzeroObject->AddAtAndExpand(fResolution, 27);
107 fTzeroObject->AddAtAndExpand(fTzeroORA, 28);
108 fTzeroObject->AddAtAndExpand(fTzeroORC, 29);
109 fTzeroObject->AddAtAndExpand(fT0vsNtracks, 30);
110 fTzeroObject->AddAtAndExpand(fOrAvsNtracks,31);
111 fTzeroObject->AddAtAndExpand(fOrCvsNtracks, 32);
112 fTzeroObject->AddAtAndExpand(fEffC, 33);
113 fTzeroObject->AddAtAndExpand(fEffA, 34);
114 fTzeroObject->AddAtAndExpand(fEffAC, 35);
115 fTzeroObject->AddAtAndExpand(ftracksEffSPD, 36);
117 PostData(1, fTzeroObject);
121 //________________________________________________________________________
122 void AliT0AnalysisTaskQA::UserExec(Option_t *)
125 // Called for each event
128 fESD = dynamic_cast<AliESDEvent*>(InputEvent());
130 printf("ERROR: fESD not available\n");
134 fRunNumber = fESD->GetRunNumber() ;
136 const Double32_t* time = fESD->GetT0time();
137 const Double32_t* amplitude = fESD->GetT0amplitude();
140 for (Int_t i=0; i<24; i++) {
141 if(time[i]<9999 &&abs(time[i])>1e-8 && amplitude[i]<9999&&abs(amplitude[i])>1e-8 )
143 // cout<<"time "<<time[i]<<" amplitude "<<amplitude[i]<<endl;
144 fTimeVSAmplitude[i]->Fill(amplitude[i],time[i]);
145 fCFDVSPmtId->Fill(i,time[i]);
149 const Double32_t* mean = fESD->GetT0TOF();
150 Double32_t orA = mean[1];
151 Double32_t orC = mean[2];
152 Int_t ntracks = fESD->GetNumberOfTracks();
155 fTzeroORA->Fill(orA);
156 fOrAvsNtracks->Fill(ntracks, orA);
157 fEffA->Fill(ntracks);
160 fTzeroORC->Fill(orC);
161 fOrCvsNtracks->Fill(ntracks, orC);
162 fEffC->Fill(ntracks);
164 if(orA<99999 && orC<99999) {
165 fResolution->Fill((orA-orC)/2.);
166 fTzeroORAplusORC->Fill(mean[0]);
167 fEffAC->Fill(ntracks);
168 fT0vsNtracks->Fill(ntracks, mean[0]);
171 Double32_t t0vertex = fESD->GetT0zVertex();
172 Double32_t esdzvertex;
173 const AliESDVertex * esdvertex = fESD->GetPrimaryVertex();
176 nofcontrib=esdvertex->GetNContributors();
177 if(nofcontrib>0) ftracksEffSPD->Fill(ntracks);
178 if(esdvertex && t0vertex<999)
182 esdzvertex=esdvertex->GetZv();
183 fSPDVertexVST0Vertex->Fill(t0vertex,esdzvertex);
187 PostData(1, fTzeroObject);
189 //________________________________________________________________________
190 void AliT0AnalysisTaskQA::Terminate(Option_t *)
193 // Called once at the end of the query