new reconstruction with pass0 calibration

[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)\r
{\r
  // Constructor\r
\r
  for( int ip=0; ip < 24; ip++){\r
    fTimeDiff[ip] = 0;\r
    fCFD[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)\r
{\r
  // Constructor\r
 \r
  for( int ip=0; ip < 24; ip++){\r
    fTimeDiff[ip] = 0;\r
    fCFD[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
  // 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
  delete [] fTimeDiff;\r
  delete [] fCFD;\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
  for (Int_t i=0; i<24; i++) {\r
    fTimeDiff[i]   = new TH1F (Form("CFD1minCFD%d",i+1),"fTimeDiff",300, -300, 300);\r
    //   fCFD[i]        = new TH1F(Form("CFD%d",i+1),"CFD",500, 2000, 3000);//6000, 7000);\r
    fCFD[i]        = new TH1F(Form("CFD%d",i+1),"CFD",500, -1000, 1000);//6000, 7000);\r
  }\r
\r
  fTzeroORAplusORC = new TH1F("fTzeroORAplusORC","ORA+ORC /2",400,-2000,2000);   //or A plus or C \r
  fResolution      = new TH1F("fResolution","fResolution",400,-2000,2000);// or A minus or C spectrum\r
  fTzeroORA        = new TH1F("fTzeroORA","fTzeroORA",400,-2000,2000);// or A spectrum\r
  fTzeroORC        = new TH1F("fTzeroORC","fTzeroORC",400,-2000,2000);// 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
  fESD = dynamic_cast<AliESDEvent*>(InputEvent());\r
  if (!fESD) {\r
    printf("ERROR: fESD not available\n");\r
    return;\r
  }\r
\r
  fRunNumber =  fESD->GetRunNumber() ; \r
\r
  const Double32_t* time = fESD->GetT0time();\r
  for (Int_t i=0; i<12; i++) {\r
     if( time[i] != 0 ){\r
     printf(" %i time %f \n", i, time[i] );\r
    fCFD[i]->Fill( time[i]);\r
      if(  time[0] != 0 ) \r
	fTimeDiff[i]->Fill( time[i]-time[0]);\r
    }\r
  }\r
  for (Int_t i=12; i<24; i++) {\r
    if( time[i] != 0) {\r
      fCFD[i]->Fill( time[i]);\r
      if( time[12] != 0 ) \r
	fTimeDiff[i]->Fill( time[i]-time[12]);\r
    }\r
  }\r
  const Double32_t* mean = fESD->GetT0TOF();\r
  Double32_t meanTOF = mean[0] ;\r
  Double32_t orA = mean[1] ;\r
  Double32_t orC = mean[2];\r
 \r
  if(orA<9999) fTzeroORA->Fill(orA);\r
  if(orC<9999) fTzeroORC->Fill(orC);\r
  if(orA<9999 && orC<9999) fResolution->Fill((orA-orC)/2.);\r
  if(orA<9999 && orC<9999) fTzeroORAplusORC->Fill(meanTOF); \r
\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