[u/mrichter/AliRoot.git] / T0 / AliT0CalibOffsetChannelsTask.cxx

#include "TChain.h"\r
#include "TTree.h"\r
#include "TH1D.h"\r
#include "TF1.h"\r
#include "TCanvas.h"\r
#include "TObjArray.h"\r
\r
#include "AliAnalysisTask.h"\r
#include "AliAnalysisManager.h"\r
\r
#include "AliESDEvent.h"\r
#include "AliESDInputHandler.h"\r
\r
#include "AliT0CalibOffsetChannelsTask.h"\r
\r
//#include "AliCDBMetaData.h"\r
//#include "AliCDBId.h"\r
//#include "AliCDBEntry.h"\r
//#include "AliCDBManager.h"\r
//#include "AliCDBStorage.h"\r
\r
// Task should calculate channels offset \r
// Authors: Alla \r
\r
ClassImp(AliT0CalibOffsetChannelsTask)\r
//________________________________________________________________________\r
AliT0CalibOffsetChannelsTask::AliT0CalibOffsetChannelsTask() \r
  : AliAnalysisTaskSE(),  fESD(0x0), fTzeroObject(0x0),\r
  fTzeroORA(0x0), fTzeroORC(0x0), fResolution(0x0), fTzeroORAplusORC(0x0),\r
    fRunNumber(0),fRefPMTA(12), fRefPMTC(0)\r
{\r
  // Constructor\r
\r
  for( int ip=0; ip < 24; ip++){\r
    fTimeDiff[ip] = 0;\r
    fCFD[ip]      = 0;\r
    fCDBdelays[ip]= 0;\r
    fCDBcfds[ip]= 0;\r
    if (ip<4 ) fCDBT0s[ip]= 0;\r
  }\r
\r
  // Define input and output slots here\r
  // Input slot #0 works with a TChain\r
  //  DefineInput(0,  TChain::Class());\r
  //  DefineOutput(1, TObjArray::Class());\r
}\r
\r
\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),fRefPMTA(12), fRefPMTC(0)\r
{\r
  // Constructor\r
 \r
  for( int ip=0; ip < 24; ip++){\r
    fTimeDiff[ip] = 0;\r
    fCFD[ip]      = 0;\r
    fCDBdelays[ip]= 0;\r
    fCDBcfds[ip]= 0;\r
    if (ip<4 ) fCDBT0s[ip]= 0;\r
\r
  }\r
 \r
  // Define input and output slots here\r
  // Input slot #0 works with a TChain\r
  DefineInput(0, TChain::Class());\r
  DefineOutput(1, TObjArray::Class());\r
  // Output slot #0 id reserved by the base class for AOD\r
  // Output slot #1 writes into a TH1 container\r
}\r
\r
//________________________________________________________________________\r
AliT0CalibOffsetChannelsTask::~AliT0CalibOffsetChannelsTask() \r
{\r
  // Destructor\r
  // printf("AliT0CalibOffsetChannels~AliT0CalibOffsetChannels() ");\r
  delete fTzeroORA;\r
  delete fTzeroORC;\r
  delete fResolution;\r
  delete fTzeroORAplusORC;\r
  for( Int_t  ip=0; ip < 24; ip++){\r
    delete fTimeDiff[ip];\r
    delete fCFD[ip];\r
  }\r
  \r
  delete fTzeroObject;\r
}\r
\r
//________________________________________________________________________\r
/*void AliT0CalibOffsetChannelsTaskX::ConnectInputData(Option_t *) {\r
  //\r
  //\r
  //\r
  TTree* tree=dynamic_cast<TTree*>(GetInputData(0));\r
  if (!tree) {\r
    printf("ERROR: Could not read chain from input slot 0");\r
  } \r
  else {\r
    AliESDInputHandler *esdH = dynamic_cast<AliESDInputHandler*> (AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler());\r
    if (!esdH) {\r
      printf ("ERROR: Could not get ESDInputHandler");\r
    } \r
    else {\r
      fESD = esdH->GetEvent();\r
      printf ("*** CONNECTED NEW EVENT ****");\r
    }\r
  }\r
}\r
*/\r
//________________________________________________________________________\r
void AliT0CalibOffsetChannelsTask::UserCreateOutputObjects()\r
{\r
  // Create histograms\r
  Float_t low = fCDBcfds[fRefPMTC] - 500;\r
  Float_t high = fCDBcfds[fRefPMTA] + 500;\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,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,-4000,4000);   //or A plus or C \r
  fResolution      = new TH1F("fResolution","fResolution",200,-2000,2000);// or A minus 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
  fTzeroObject->SetOwner(kTRUE);\r
  \r
  for (Int_t i=0; i<24; i++)\r
    fTzeroObject->AddAtAndExpand(fTimeDiff[i],i);\r
\r
  for (Int_t i=0; i<24; i++)\r
    fTzeroObject->AddAtAndExpand(fCFD[i],i+24); //24 - 48\r
\r
  fTzeroObject->AddAtAndExpand(fTzeroORAplusORC, 48);\r
  fTzeroObject->AddAtAndExpand(fResolution, 49);\r
  fTzeroObject->AddAtAndExpand(fTzeroORA, 50);\r
  fTzeroObject->AddAtAndExpand(fTzeroORC, 51);\r
\r
  PostData(1, fTzeroObject);\r
  // Called once\r
}\r
\r
//________________________________________________________________________\r
void AliT0CalibOffsetChannelsTask::UserExec(Option_t *) \r
{\r
  // Main loop\r
  // Called for each event\r
\r
  // Post output data.\r
\r
  fESD = dynamic_cast<AliESDEvent*>(InputEvent());\r
  if (!fESD) {\r
    printf("ERROR: fESD not available\n");\r
    return;\r
  }\r
  Int_t trigT0 = fESD->GetT0Trig();\r
  Bool_t eq = kTRUE;\r
  fRunNumber =  fESD->GetRunNumber() ;\r
  if( fRunNumber<165747) eq = kFALSE;\r
    \r
  const Double32_t* time = fESD->GetT0time();\r
  const Double32_t* amp = fESD->GetT0amplitude();\r
  \r
  Double32_t diff;\r
  for (Int_t i=0; i<24; i++) {\r
    if( time[i] > 0  && amp[i]>0.1 ){\r
      if (eq)	{\r
	fCFD[i]->Fill( time[i] );//////!!!!!\r
	if(  time[fRefPMTC] > 0 && i<12)   {\r
	  diff =  time[i]-time[fRefPMTC];\r
	  fTimeDiff[i]->Fill( diff);\r
	}\r
	if(  time[fRefPMTA] >0  && i>11)  {\r
	  diff =  time[i]-time[fRefPMTA] ;\r
	  fTimeDiff[i]->Fill( diff);\r
	}\r
      } //eq=1\r
      else  {\r
	fCFD[i]->Fill( time[i] + fCDBdelays[i] );\r
	if(  time[fRefPMTC] > 0 && i<12) {\r
	  diff =  time[i]-time[fRefPMTC] + fCDBdelays[i];\r
	  fTimeDiff[i]->Fill( diff);\r
	} //C\r
	if(  time[fRefPMTA] >0  && i>11) {\r
	  diff =  time[i]-time[fRefPMTA] + fCDBdelays[i];\r
	  fTimeDiff[i]->Fill( diff);\r
	} //A\r
      } //eq=0\r
    }\r
    \r
  }\r
  if (trigT0>5) {\r
    const Double32_t* mean = fESD->GetT0TOF();\r
    Double32_t meanTOF = mean[0]  +  fCDBT0s[0] ;\r
    Double32_t orA = mean[1]  +  fCDBT0s[1] ;\r
    Double32_t orC = mean[2] + fCDBT0s[2] ;\r
    \r
    if(orA<99999) fTzeroORA->Fill(orA);\r
    if(orC<99999) fTzeroORC->Fill(orC);\r
    if(orA<99999 && orC<99999) fResolution->Fill((orA-orC)/2.);\r
    if(orA<99999 && orC<99999) fTzeroORAplusORC->Fill(meanTOF); \r
  } //if TVDC on\r
    //  printf("%f   %f  %f\n",orA,orC,meanTOF);\r
  PostData(1, fTzeroObject);\r
  }      \r
 //________________________________________________________________________\r
  void AliT0CalibOffsetChannelsTask::Terminate(Option_t *) \r
{\r
  \r
   // Called once at the end of the query\r
}\r
 \r
Commit	Line	Data
b5e0299c	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
4f0fdf74	30	fRunNumber(0),fRefPMTA(12), fRefPMTC(0)\r
b5e0299c	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
b9dd81a7	37	fCDBdelays[ip]= 0;\r
b9dd81a7	38	fCDBcfds[ip]= 0;\r
e36e79d8	39	if (ip<4 ) fCDBT0s[ip]= 0;\r
b5e0299c	40	}\r
	41	\r
	42	// Define input and output slots here\r
	43	// Input slot #0 works with a TChain\r
	44	// DefineInput(0, TChain::Class());\r
	45	// DefineOutput(1, TObjArray::Class());\r
	46	}\r
	47	\r
	48	\r
	49	//________________________________________________________________________\r
	50	AliT0CalibOffsetChannelsTask::AliT0CalibOffsetChannelsTask(const char *name) \r
	51	: AliAnalysisTaskSE(name), fESD(0), fTzeroObject(0),\r
	52	fTzeroORA(0x0), fTzeroORC(0x0), fResolution(0x0), fTzeroORAplusORC(0x0),\r
4f0fdf74	53	fRunNumber(0),fRefPMTA(12), fRefPMTC(0)\r
b5e0299c	54	{\r
	55	// Constructor\r
	56	\r
	57	for( int ip=0; ip < 24; ip++){\r
	58	fTimeDiff[ip] = 0;\r
	59	fCFD[ip] = 0;\r
b9dd81a7	60	fCDBdelays[ip]= 0;\r
b9dd81a7	61	fCDBcfds[ip]= 0;\r
e36e79d8	62	if (ip<4 ) fCDBT0s[ip]= 0;\r
b9dd81a7	63	\r
b5e0299c	64	}\r
	65	\r
	66	// Define input and output slots here\r
	67	// Input slot #0 works with a TChain\r
	68	DefineInput(0, TChain::Class());\r
	69	DefineOutput(1, TObjArray::Class());\r
	70	// Output slot #0 id reserved by the base class for AOD\r
	71	// Output slot #1 writes into a TH1 container\r
	72	}\r
	73	\r
	74	//________________________________________________________________________\r
	75	AliT0CalibOffsetChannelsTask::~AliT0CalibOffsetChannelsTask() \r
	76	{\r
	77	// Destructor\r
	78	// printf("AliT0CalibOffsetChannels~AliT0CalibOffsetChannels() ");\r
	79	delete fTzeroORA;\r
	80	delete fTzeroORC;\r
	81	delete fResolution;\r
	82	delete fTzeroORAplusORC;\r
b0ab3f59	83	for( Int_t ip=0; ip < 24; ip++){\r
	84	delete fTimeDiff[ip];\r
	85	delete fCFD[ip];\r
	86	}\r
	87	\r
b5e0299c	88	delete fTzeroObject;\r
	89	}\r
	90	\r
	91	//________________________________________________________________________\r
	92	/void AliT0CalibOffsetChannelsTaskX::ConnectInputData(Option_t ) {\r
	93	//\r
	94	//\r
	95	//\r
	96	TTree* tree=dynamic_cast<TTree*>(GetInputData(0));\r
	97	if (!tree) {\r
	98	printf("ERROR: Could not read chain from input slot 0");\r
	99	} \r
	100	else {\r
	101	AliESDInputHandler esdH = dynamic_cast<AliESDInputHandler> (AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler());\r
	102	if (!esdH) {\r
	103	printf ("ERROR: Could not get ESDInputHandler");\r
	104	} \r
	105	else {\r
	106	fESD = esdH->GetEvent();\r
	107	printf ("* CONNECTED NEW EVENT **");\r
	108	}\r
	109	}\r
	110	}\r
	111	*/\r
	112	//________________________________________________________________________\r
	113	void AliT0CalibOffsetChannelsTask::UserCreateOutputObjects()\r
	114	{\r
	115	// Create histograms\r
4f0fdf74	116	Float_t low = fCDBcfds[fRefPMTC] - 500;\r
4f0fdf74	117	Float_t high = fCDBcfds[fRefPMTA] + 500;\r
b5e0299c	118	for (Int_t i=0; i<24; i++) {\r
b9dd81a7	119	fTimeDiff[i] = new TH1F (Form("CFD1minCFD%d",i+1),"fTimeDiff",150, -300, 300);\r
4f0fdf74	120	fCFD[i] = new TH1F(Form("CFD%d",i+1),"CFD",250,low, high);//6000, 7000);\r
b9dd81a7	121	// fCFD[i] = new TH1F(Form("CFD%d",i+1),"CFD",250, -1000, 1000);//6000, 7000);\r
b5e0299c	122	}\r
b5e0299c	123	\r
a4d36291	124	fTzeroORAplusORC = new TH1F("fTzeroORAplusORC","ORA+ORC /2",200,-4000,4000); //or A plus or C \r
13e2fbbd	125	fResolution = new TH1F("fResolution","fResolution",200,-2000,2000);// or A minus or C spectrum\r
e5835843	126	fTzeroORA = new TH1F("fTzeroORA","fTzeroORA",200,-4000,4000);// or A spectrum\r
e5835843	127	fTzeroORC = new TH1F("fTzeroORC","fTzeroORC",200,-4000,4000);// or C spectrum\r
b5e0299c	128	\r
	129	\r
	130	fTzeroObject = new TObjArray(0);\r
	131	fTzeroObject->SetOwner(kTRUE);\r
	132	\r
	133	for (Int_t i=0; i<24; i++)\r
	134	fTzeroObject->AddAtAndExpand(fTimeDiff[i],i);\r
	135	\r
	136	for (Int_t i=0; i<24; i++)\r
	137	fTzeroObject->AddAtAndExpand(fCFD[i],i+24); //24 - 48\r
	138	\r
	139	fTzeroObject->AddAtAndExpand(fTzeroORAplusORC, 48);\r
	140	fTzeroObject->AddAtAndExpand(fResolution, 49);\r
	141	fTzeroObject->AddAtAndExpand(fTzeroORA, 50);\r
	142	fTzeroObject->AddAtAndExpand(fTzeroORC, 51);\r
	143	\r
	144	PostData(1, fTzeroObject);\r
	145	// Called once\r
	146	}\r
	147	\r
	148	//________________________________________________________________________\r
	149	void AliT0CalibOffsetChannelsTask::UserExec(Option_t *) \r
	150	{\r
	151	// Main loop\r
	152	// Called for each event\r
	153	\r
	154	// Post output data.\r
a4d36291	155	\r
b5e0299c	156	fESD = dynamic_cast<AliESDEvent*>(InputEvent());\r
	157	if (!fESD) {\r
	158	printf("ERROR: fESD not available\n");\r
	159	return;\r
	160	}\r
a4d36291	161	Int_t trigT0 = fESD->GetT0Trig();\r
	162	Bool_t eq = kTRUE;\r
	163	fRunNumber = fESD->GetRunNumber() ;\r
	164	if( fRunNumber<165747) eq = kFALSE;\r
	165	\r
b5e0299c	166	const Double32_t* time = fESD->GetT0time();\r
b9dd81a7	167	const Double32_t* amp = fESD->GetT0amplitude();\r
a4d36291	168	\r
b9dd81a7	169	Double32_t diff;\r
a4d36291	170	for (Int_t i=0; i<24; i++) {\r
	171	if( time[i] > 0 && amp[i]>0.1 ){\r
	172	if (eq) {\r
13e2fbbd	173	fCFD[i]->Fill( time[i] );//////!!!!!\r
a4d36291	174	if( time[fRefPMTC] > 0 && i<12) {\r
	175	diff = time[i]-time[fRefPMTC];\r
	176	fTimeDiff[i]->Fill( diff);\r
	177	}\r
	178	if( time[fRefPMTA] >0 && i>11) {\r
	179	diff = time[i]-time[fRefPMTA] ;\r
	180	fTimeDiff[i]->Fill( diff);\r
	181	}\r
	182	} //eq=1\r
	183	else {\r
	184	fCFD[i]->Fill( time[i] + fCDBdelays[i] );\r
	185	if( time[fRefPMTC] > 0 && i<12) {\r
	186	diff = time[i]-time[fRefPMTC] + fCDBdelays[i];\r
	187	fTimeDiff[i]->Fill( diff);\r
	188	} //C\r
	189	if( time[fRefPMTA] >0 && i>11) {\r
	190	diff = time[i]-time[fRefPMTA] + fCDBdelays[i];\r
	191	fTimeDiff[i]->Fill( diff);\r
	192	} //A\r
	193	} //eq=0\r
b5e0299c	194	}\r
a4d36291	195	\r
b5e0299c	196	}\r
a4d36291	197	if (trigT0>5) {\r
	198	const Double32_t* mean = fESD->GetT0TOF();\r
	199	Double32_t meanTOF = mean[0] + fCDBT0s[0] ;\r
	200	Double32_t orA = mean[1] + fCDBT0s[1] ;\r
	201	Double32_t orC = mean[2] + fCDBT0s[2] ;\r
	202	\r
13e2fbbd	203	if(orA<99999) fTzeroORA->Fill(orA);\r
	204	if(orC<99999) fTzeroORC->Fill(orC);\r
	205	if(orA<99999 && orC<99999) fResolution->Fill((orA-orC)/2.);\r
	206	if(orA<99999 && orC<99999) fTzeroORAplusORC->Fill(meanTOF); \r
a4d36291	207	} //if TVDC on\r
a4d36291	208	// printf("%f %f %f\n",orA,orC,meanTOF);\r
b5e0299c	209	PostData(1, fTzeroObject);\r
a4d36291	210	} \r
b5e0299c	211	//________________________________________________________________________\r
a4d36291	212	void AliT0CalibOffsetChannelsTask::Terminate(Option_t *) \r
b5e0299c	213	{\r
	214	\r
	215	// Called once at the end of the query\r
	216	}\r
	217	\r