removed separate histograms for Calib and physics data
[u/mrichter/AliRoot.git] / T0 / AliT0CalibOffsetChannelsTask.cxx
1 #include "TChain.h"\r
2 #include "TTree.h"\r
3 #include "TH1D.h"\r
4 #include "TF1.h"\r
5 #include "TCanvas.h"\r
6 #include "TObjArray.h"\r
7 \r
8 #include "AliAnalysisTask.h"\r
9 #include "AliAnalysisManager.h"\r
10 \r
11 #include "AliESDEvent.h"\r
12 #include "AliESDInputHandler.h"\r
13 \r
14 #include "AliT0CalibOffsetChannelsTask.h"\r
15 \r
16 //#include "AliCDBMetaData.h"\r
17 //#include "AliCDBId.h"\r
18 //#include "AliCDBEntry.h"\r
19 //#include "AliCDBManager.h"\r
20 //#include "AliCDBStorage.h"\r
21 \r
22 // Task should calculate channels offset \r
23 // Authors: Alla \r
24 \r
25 ClassImp(AliT0CalibOffsetChannelsTask)\r
26 //________________________________________________________________________\r
27 AliT0CalibOffsetChannelsTask::AliT0CalibOffsetChannelsTask() \r
28   : AliAnalysisTaskSE(),  fESD(0x0), fTzeroObject(0x0),\r
29   fTzeroORA(0x0), fTzeroORC(0x0), fResolution(0x0), fTzeroORAplusORC(0x0),\r
30   fRunNumber(0)\r
31 {\r
32   // Constructor\r
33 \r
34   for( int ip=0; ip < 24; ip++){\r
35     fTimeDiff[ip] = 0;\r
36     fCFD[ip]      = 0;\r
37   }\r
38 \r
39   // Define input and output slots here\r
40   // Input slot #0 works with a TChain\r
41   //  DefineInput(0,  TChain::Class());\r
42   //  DefineOutput(1, TObjArray::Class());\r
43 }\r
44 \r
45 \r
46 //________________________________________________________________________\r
47 AliT0CalibOffsetChannelsTask::AliT0CalibOffsetChannelsTask(const char *name) \r
48   : AliAnalysisTaskSE(name), fESD(0), fTzeroObject(0),\r
49   fTzeroORA(0x0), fTzeroORC(0x0), fResolution(0x0), fTzeroORAplusORC(0x0),\r
50   fRunNumber(0)\r
51 {\r
52   // Constructor\r
53  \r
54   for( int ip=0; ip < 24; ip++){\r
55     fTimeDiff[ip] = 0;\r
56     fCFD[ip]      = 0;\r
57   }\r
58  \r
59   // Define input and output slots here\r
60   // Input slot #0 works with a TChain\r
61   DefineInput(0, TChain::Class());\r
62   DefineOutput(1, TObjArray::Class());\r
63   // Output slot #0 id reserved by the base class for AOD\r
64   // Output slot #1 writes into a TH1 container\r
65 }\r
66 \r
67 //________________________________________________________________________\r
68 AliT0CalibOffsetChannelsTask::~AliT0CalibOffsetChannelsTask() \r
69 {\r
70   // Destructor\r
71   // printf("AliT0CalibOffsetChannels~AliT0CalibOffsetChannels() ");\r
72   delete fTzeroORA;\r
73   delete fTzeroORC;\r
74   delete fResolution;\r
75   delete fTzeroORAplusORC;\r
76   for( Int_t  ip=0; ip < 24; ip++){\r
77     delete fTimeDiff[ip];\r
78     delete fCFD[ip];\r
79   }\r
80   \r
81   delete fTzeroObject;\r
82 }\r
83 \r
84 //________________________________________________________________________\r
85 /*void AliT0CalibOffsetChannelsTaskX::ConnectInputData(Option_t *) {\r
86   //\r
87   //\r
88   //\r
89   TTree* tree=dynamic_cast<TTree*>(GetInputData(0));\r
90   if (!tree) {\r
91     printf("ERROR: Could not read chain from input slot 0");\r
92   } \r
93   else {\r
94     AliESDInputHandler *esdH = dynamic_cast<AliESDInputHandler*> (AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler());\r
95     if (!esdH) {\r
96       printf ("ERROR: Could not get ESDInputHandler");\r
97     } \r
98     else {\r
99       fESD = esdH->GetEvent();\r
100       printf ("*** CONNECTED NEW EVENT ****");\r
101     }\r
102   }\r
103 }\r
104 */\r
105 //________________________________________________________________________\r
106 void AliT0CalibOffsetChannelsTask::UserCreateOutputObjects()\r
107 {\r
108   // Create histograms\r
109   for (Int_t i=0; i<24; i++) {\r
110     fTimeDiff[i]   = new TH1F (Form("CFD1minCFD%d",i+1),"fTimeDiff",300, -300, 300);\r
111     fCFD[i]        = new TH1F(Form("CFD%d",i+1),"CFD",500, 2000, 3000);//6000, 7000);\r
112     // fCFD[i]        = new TH1F(Form("CFD%d",i+1),"CFD",500, -1000, 1000);//6000, 7000);\r
113   }\r
114 \r
115   fTzeroORAplusORC = new TH1F("fTzeroORAplusORC","ORA+ORC /2",400,-2000,2000);   //or A plus or C \r
116   fResolution      = new TH1F("fResolution","fResolution",400,-2000,2000);// or A minus or C spectrum\r
117   fTzeroORA        = new TH1F("fTzeroORA","fTzeroORA",400,-2000,2000);// or A spectrum\r
118   fTzeroORC        = new TH1F("fTzeroORC","fTzeroORC",400,-2000,2000);// or C spectrum\r
119 \r
120   \r
121   fTzeroObject     = new TObjArray(0);\r
122   fTzeroObject->SetOwner(kTRUE);\r
123   \r
124   for (Int_t i=0; i<24; i++)\r
125     fTzeroObject->AddAtAndExpand(fTimeDiff[i],i);\r
126 \r
127   for (Int_t i=0; i<24; i++)\r
128     fTzeroObject->AddAtAndExpand(fCFD[i],i+24); //24 - 48\r
129 \r
130   fTzeroObject->AddAtAndExpand(fTzeroORAplusORC, 48);\r
131   fTzeroObject->AddAtAndExpand(fResolution, 49);\r
132   fTzeroObject->AddAtAndExpand(fTzeroORA, 50);\r
133   fTzeroObject->AddAtAndExpand(fTzeroORC, 51);\r
134 \r
135   PostData(1, fTzeroObject);\r
136   // Called once\r
137 }\r
138 \r
139 //________________________________________________________________________\r
140 void AliT0CalibOffsetChannelsTask::UserExec(Option_t *) \r
141 {\r
142   // Main loop\r
143   // Called for each event\r
144 \r
145   // Post output data.\r
146   fESD = dynamic_cast<AliESDEvent*>(InputEvent());\r
147   if (!fESD) {\r
148     printf("ERROR: fESD not available\n");\r
149     return;\r
150   }\r
151 \r
152   fRunNumber =  fESD->GetRunNumber() ; \r
153 \r
154   const Double32_t* time = fESD->GetT0time();\r
155   for (Int_t i=0; i<12; i++) {\r
156     if( time[i] != 0 ){\r
157       fCFD[i]->Fill( time[i]);\r
158       if(  time[0] != 0 ) \r
159         fTimeDiff[i]->Fill( time[i]-time[0]);\r
160     }\r
161   }\r
162   for (Int_t i=12; i<24; i++) {\r
163     if( time[i] != 0) {\r
164       fCFD[i]->Fill( time[i]);\r
165       if( time[12] != 0 ) \r
166         fTimeDiff[i]->Fill( time[i]-time[12]);\r
167     }\r
168   }\r
169   const Double32_t* mean = fESD->GetT0TOF();\r
170   Double32_t meanTOF = mean[0] ;\r
171   Double32_t orA = mean[1] ;\r
172   Double32_t orC = mean[2];\r
173  \r
174   if(orA<99999) fTzeroORA->Fill(orA);\r
175   if(orC<99999) fTzeroORC->Fill(orC);\r
176   if(orA<99999 && orC<99999) fResolution->Fill((orA-orC)/2.);\r
177   if(orA<99999 && orC<99999) fTzeroORAplusORC->Fill(meanTOF); \r
178 \r
179   //  printf("%f   %f  %f\n",orA,orC,meanTOF);\r
180   PostData(1, fTzeroObject);\r
181 }      \r
182  //________________________________________________________________________\r
183 void AliT0CalibOffsetChannelsTask::Terminate(Option_t *) \r
184 {\r
185   \r
186    // Called once at the end of the query\r
187 }\r
188  \r