8 #include "AliAnalysisTask.h"
9 #include "AliAnalysisManager.h"
11 #include "AliESDEvent.h"
12 #include "AliESDInputHandler.h"
14 #include "AliT0CalibOffsetChannelsTask.h"
16 //#include "AliCDBMetaData.h"
17 //#include "AliCDBId.h"
18 //#include "AliCDBEntry.h"
19 //#include "AliCDBManager.h"
20 //#include "AliCDBStorage.h"
22 // Task should calculate channels offset
25 ClassImp(AliT0CalibOffsetChannelsTask)
26 //________________________________________________________________________
27 AliT0CalibOffsetChannelsTask::AliT0CalibOffsetChannelsTask()
28 : AliAnalysisTaskSE(), fESD(0x0), fTzeroObject(0x0),
29 fTzeroORA(0x0), fTzeroORC(0x0), fResolution(0x0), fTzeroORAplusORC(0x0),
30 fRunNumber(0),fRefPMTA(12), fRefPMTC(0)
34 for( int ip=0; ip < 24; ip++){
39 if (ip<4 ) fCDBT0s[ip]= 0;
42 // Define input and output slots here
43 // Input slot #0 works with a TChain
44 // DefineInput(0, TChain::Class());
45 // DefineOutput(1, TObjArray::Class());
49 //________________________________________________________________________
50 AliT0CalibOffsetChannelsTask::AliT0CalibOffsetChannelsTask(const char *name)
51 : AliAnalysisTaskSE(name), fESD(0), fTzeroObject(0),
52 fTzeroORA(0x0), fTzeroORC(0x0), fResolution(0x0), fTzeroORAplusORC(0x0),
53 fRunNumber(0),fRefPMTA(12), fRefPMTC(0)
57 for( int ip=0; ip < 24; ip++){
62 if (ip<4 ) fCDBT0s[ip]= 0;
66 // Define input and output slots here
67 // Input slot #0 works with a TChain
68 DefineInput(0, TChain::Class());
69 DefineOutput(1, TObjArray::Class());
70 // Output slot #0 id reserved by the base class for AOD
71 // Output slot #1 writes into a TH1 container
74 //________________________________________________________________________
75 AliT0CalibOffsetChannelsTask::~AliT0CalibOffsetChannelsTask()
78 // printf("AliT0CalibOffsetChannels~AliT0CalibOffsetChannels() ");
82 delete fTzeroORAplusORC;
83 for( Int_t ip=0; ip < 24; ip++){
91 //________________________________________________________________________
92 /*void AliT0CalibOffsetChannelsTaskX::ConnectInputData(Option_t *) {
96 TTree* tree=dynamic_cast<TTree*>(GetInputData(0));
98 printf("ERROR: Could not read chain from input slot 0");
101 AliESDInputHandler *esdH = dynamic_cast<AliESDInputHandler*> (AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler());
103 printf ("ERROR: Could not get ESDInputHandler");
106 fESD = esdH->GetEvent();
107 printf ("*** CONNECTED NEW EVENT ****");
112 //________________________________________________________________________
113 void AliT0CalibOffsetChannelsTask::UserCreateOutputObjects()
116 Float_t low = fCDBcfds[fRefPMTC] - 500;
117 Float_t high = fCDBcfds[fRefPMTA] + 500;
118 for (Int_t i=0; i<24; i++) {
119 fTimeDiff[i] = new TH1F (Form("CFD1minCFD%d",i+1),"fTimeDiff",150, -300, 300);
120 fCFD[i] = new TH1F(Form("CFD%d",i+1),"CFD",250,low, high);//6000, 7000);
121 // fCFD[i] = new TH1F(Form("CFD%d",i+1),"CFD",250, -1000, 1000);//6000, 7000);
124 fTzeroORAplusORC = new TH1F("fTzeroORAplusORC","ORA+ORC /2",200,-4000,4000); //or A plus or C
125 fResolution = new TH1F("fResolution","fResolution",200,-2000,2000);// or A minus or C spectrum
126 fTzeroORA = new TH1F("fTzeroORA","fTzeroORA",200,-4000,4000);// or A spectrum
127 fTzeroORC = new TH1F("fTzeroORC","fTzeroORC",200,-4000,4000);// or C spectrum
130 fTzeroObject = new TObjArray(0);
131 fTzeroObject->SetOwner(kTRUE);
133 for (Int_t i=0; i<24; i++)
134 fTzeroObject->AddAtAndExpand(fTimeDiff[i],i);
136 for (Int_t i=0; i<24; i++)
137 fTzeroObject->AddAtAndExpand(fCFD[i],i+24); //24 - 48
139 fTzeroObject->AddAtAndExpand(fTzeroORAplusORC, 48);
140 fTzeroObject->AddAtAndExpand(fResolution, 49);
141 fTzeroObject->AddAtAndExpand(fTzeroORA, 50);
142 fTzeroObject->AddAtAndExpand(fTzeroORC, 51);
144 PostData(1, fTzeroObject);
148 //________________________________________________________________________
149 void AliT0CalibOffsetChannelsTask::UserExec(Option_t *)
152 // Called for each event
156 fESD = dynamic_cast<AliESDEvent*>(InputEvent());
158 printf("ERROR: fESD not available\n");
161 Int_t trigT0 = fESD->GetT0Trig();
163 fRunNumber = fESD->GetRunNumber() ;
164 if( fRunNumber<165747) eq = kFALSE;
166 const Double32_t* time = fESD->GetT0time();
167 const Double32_t* amp = fESD->GetT0amplitude();
170 for (Int_t i=0; i<24; i++) {
171 if( time[i] > 0 && amp[i]>0.1 ){
173 fCFD[i]->Fill( time[i] );//////!!!!!
174 if( time[fRefPMTC] > 0 && i<12) {
175 diff = time[i]-time[fRefPMTC];
176 fTimeDiff[i]->Fill( diff);
178 if( time[fRefPMTA] >0 && i>11) {
179 diff = time[i]-time[fRefPMTA] ;
180 fTimeDiff[i]->Fill( diff);
184 fCFD[i]->Fill( time[i] + fCDBdelays[i] );
185 if( time[fRefPMTC] > 0 && i<12) {
186 diff = time[i]-time[fRefPMTC] + fCDBdelays[i];
187 fTimeDiff[i]->Fill( diff);
189 if( time[fRefPMTA] >0 && i>11) {
190 diff = time[i]-time[fRefPMTA] + fCDBdelays[i];
191 fTimeDiff[i]->Fill( diff);
198 const Double32_t* mean = fESD->GetT0TOF();
199 Double32_t meanTOF = mean[0] + fCDBT0s[0] ;
200 Double32_t orA = mean[1] + fCDBT0s[1] ;
201 Double32_t orC = mean[2] + fCDBT0s[2] ;
203 if(orA<99999) fTzeroORA->Fill(orA);
204 if(orC<99999) fTzeroORC->Fill(orC);
205 if(orA<99999 && orC<99999) fResolution->Fill((orA-orC)/2.);
206 if(orA<99999 && orC<99999) fTzeroORAplusORC->Fill(meanTOF);
208 // printf("%f %f %f\n",orA,orC,meanTOF);
209 PostData(1, fTzeroObject);
211 //________________________________________________________________________
212 void AliT0CalibOffsetChannelsTask::Terminate(Option_t *)
215 // Called once at the end of the query