AliT0CalibOffsetChannelsTask::AliT0CalibOffsetChannelsTask() \r
: AliAnalysisTaskSE(), fESD(0x0), fTzeroObject(0x0),\r
fTzeroORA(0x0), fTzeroORC(0x0), fResolution(0x0), fTzeroORAplusORC(0x0),\r
- fRunNumber(0)\r
+ fRunNumber(0),fRefPMTA(12), fRefPMTC(0)\r
{\r
// Constructor\r
\r
AliT0CalibOffsetChannelsTask::AliT0CalibOffsetChannelsTask(const char *name) \r
: AliAnalysisTaskSE(name), fESD(0), fTzeroObject(0),\r
fTzeroORA(0x0), fTzeroORC(0x0), fResolution(0x0), fTzeroORAplusORC(0x0),\r
- fRunNumber(0)\r
+ fRunNumber(0),fRefPMTA(12), fRefPMTC(0)\r
{\r
// Constructor\r
\r
void AliT0CalibOffsetChannelsTask::UserCreateOutputObjects()\r
{\r
// Create histograms\r
+ Float_t low = fCDBcfds[fRefPMTC] - 500;\r
+ Float_t high = fCDBcfds[fRefPMTA] + 500;\r
+ printf(" AliT0CalibOffsetChannelsTask::UserCreateOutputObjects ::low %f high %f \n", low, high);\r
for (Int_t i=0; i<24; i++) {\r
fTimeDiff[i] = new TH1F (Form("CFD1minCFD%d",i+1),"fTimeDiff",150, -300, 300);\r
- fCFD[i] = new TH1F(Form("CFD%d",i+1),"CFD",250, 2000, 3000);//6000, 7000);\r
+ fCFD[i] = new TH1F(Form("CFD%d",i+1),"CFD",250,low, high);//6000, 7000);\r
// fCFD[i] = new TH1F(Form("CFD%d",i+1),"CFD",250, -1000, 1000);//6000, 7000);\r
}\r
\r
- fTzeroORAplusORC = new TH1F("fTzeroORAplusORC","ORA+ORC /2",200,-2500,2500); //or A plus or C \r
+ fTzeroORAplusORC = new TH1F("fTzeroORAplusORC","ORA+ORC /2",400,-4000,4000); //or A plus or C \r
fResolution = new TH1F("fResolution","fResolution",400,-2500,2500);// or A minus or C spectrum\r
- fTzeroORA = new TH1F("fTzeroORA","fTzeroORA",200,-2500,2500);// or A spectrum\r
- fTzeroORC = new TH1F("fTzeroORC","fTzeroORC",200,-2500,2500);// or C spectrum\r
+ fTzeroORA = new TH1F("fTzeroORA","fTzeroORA",200,-4000,4000);// or A spectrum\r
+ fTzeroORC = new TH1F("fTzeroORC","fTzeroORC",200,-4000,4000);// or C spectrum\r
\r
\r
fTzeroObject = new TObjArray(0);\r
fCFD[i]->Fill( time[i] );\r
// printf(" time %f ocdb %f \n", time[i],fCDBcfds[i]); \r
\r
- if( time[0] != 0 ) {\r
- diff = time[i]-time[0] + fCDBdelays[i];\r
+ if( time[fRefPMTC] != 0 ) {\r
+ diff = time[i]-time[fRefPMTC] + fCDBdelays[i];\r
fTimeDiff[i]->Fill( diff);\r
}\r
}\r
if( time[i] != 0 && amp[i]>0.1) {\r
fCFD[i]->Fill( time[i]);\r
// printf(" time %f ocdb %f \n", time[i],fCDBcfds[i]); \r
- if( time[12] != 0 ) {\r
- diff = time[i]-time[12] + fCDBdelays[i];\r
+ if( time[fRefPMTA] != 0 ) {\r
+ diff = time[i]-time[fRefPMTA] + fCDBdelays[i];\r
fTimeDiff[i]->Fill( diff);\r
}\r
}\r