adding LED Ref FEE channels to PED DA and execute LED DA for calibration events only

[u/mrichter/AliRoot.git] / prod / acrcaf / qa_pp / AddTaskCalorimeterQA.C

AliAnalysisTaskParticleCorrelation *AddTaskCalorimeterQA(TString data, Bool_t kUseKinematics = kFALSE, Bool_t kPrintSettings = kFALSE)\r
{\r
  // Creates a PartCorr task for calorimeters performance studies, configures it and adds it to the analysis manager.\r
  \r
  // Get the pointer to the existing analysis manager via the static access method.\r
  //==============================================================================\r
  AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();\r
  if (!mgr) {\r
    ::Error("AddTaskPartCorr", "No analysis manager to connect to.");\r
    return NULL;\r
  }  \r
  \r
  // Check the analysis type using the event handlers connected to the analysis manager.\r
  //==============================================================================\r
  if (!mgr->GetInputEventHandler()) {\r
    ::Error("AddTaskPartCorr", "This task requires an input event handler");\r
    return NULL;\r
  }\r
   TString inputDataType = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"\r
  \r
   // Configure analysis\r
   //===========================================================================\r
  \r
  //Reader\r
  //For this particular analysis few things done by the reader.\r
  //Nothing else needs to be set.\r
  AliCaloTrackReader * reader = 0x0;\r
  if(data=="AOD")      reader = new AliCaloTrackAODReader();\r
  else if(data=="ESD") reader = new AliCaloTrackESDReader();\r
  //reader->SetDebug(10);//10 for lots of messages\r
   if(kPrintSettings) reader->Print("");\r
  \r
  if(kUseKinematics){\r
		if(inputDataType == "ESD"){\r
			reader->SwitchOnStack();          \r
			reader->SwitchOffAODMCParticles(); \r
		}\r
		else if(inputDataType == "AOD"){\r
			reader->SwitchOffStack();          \r
			reader->SwitchOnAODMCParticles(); \r
		}\r
   }\r
	\r
   reader->SetDeltaAODFileName(""); //Do not create deltaAOD file, this analysis do not create branches.\r
\r
	\r
  // ##### Analysis algorithm settings ####\r
   //Only needed now for MC data\r
   //AliFiducialCut * fidCut = new AliFiducialCut();\r
   //fidCut->DoCTSFiducialCut(kFALSE) ;\r
   //fidCut->DoEMCALFiducialCut(kTRUE) ;\r
   //fidCut->DoPHOSFiducialCut(kTRUE) ;\r
		\r
  AliAnaCalorimeterQA *emcalQA = new AliAnaCalorimeterQA();\r
  //emcalQA->SetDebug(2); //10 for lots of messages\r
  emcalQA->SetCalorimeter("EMCAL");\r
  if(kUseKinematics) emcalQA->SwitchOnDataMC() ;//Access MC stack and fill more histograms, AOD MC not implemented yet.\r
  else  emcalQA->SwitchOffDataMC() ;\r
  emcalQA->AddToHistogramsName("EMCAL_"); //Begining of histograms name\r
  //emcalQA->SetFiducialCut(fidCut);\r
  emcalQA->SwitchOffFidutialCut();\r
  if(kPrintSettings) emcalQA->Print("");	\r
  //emcalQA->SwithOnPlotsMaking();\r
  //emcalQA->GetMCAnalysisUtils()->SetDebug(10);\r
	\r
   AliAnaCalorimeterQA *phosQA = new AliAnaCalorimeterQA();\r
   //phosQA->SetDebug(2); //10 for lots of messages\r
   phosQA->SetCalorimeter("PHOS");\r
   if(kUseKinematics) phosQA->SwitchOnDataMC() ;//Access MC stack and fill more histograms, AOD MC not implemented yet.\r
   else  phosQA->SwitchOffDataMC() ;  \r
   phosQA->AddToHistogramsName("PHOS_");//Begining of histograms name\r
   //phosQA->SetFiducialCut(fidCut);\r
   phosQA->SwitchOffFidutialCut();\r
   if(kPrintSettings)phosQA->Print("");	\r
   //phosQA->GetMCAnalysisUtils()->SetDebug(10);\r
   //phosQA->SwithOnPlotsMaking();\r
  \r
  // #### Configure Maker ####\r
  AliAnaPartCorrMaker * maker = new AliAnaPartCorrMaker();\r
  maker->SetReader(reader);//pointer to reader\r
  maker->AddAnalysis(emcalQA,0);\r
  maker->AddAnalysis(phosQA,1);\r
  maker->SetAnaDebug(0)  ;\r
  maker->SwitchOnHistogramsMaker()  ;\r
  maker->SwitchOffAODsMaker()  ;\r
  if(kPrintSettings) maker->Print("");\r
 \r
\r
  printf("======================== \n");\r
  printf(" End Configuration of Calorimeter QA \n");\r
  printf("======================== \n");\r
  \r
   // Create task\r
   //===========================================================================\r
  AliAnalysisTaskParticleCorrelation * task = new AliAnalysisTaskParticleCorrelation ("CalorimeterPerformance");\r
  task->SetConfigFileName(""); //Don't configure the analysis via configuration file.\r
  task->SetDebugLevel(0);\r
  task->SetAnalysisMaker(maker);				\r
  mgr->AddTask(task);\r
  \r
  AliAnalysisDataContainer *cout_pc = mgr->CreateContainer("Calo.Performance",TList::Class(),\r
							   AliAnalysisManager::kOutputContainer, "Calo.Performance.root");\r
  \r
  // Create ONLY the output containers for the data produced by the task.\r
  // Get and connect other common input/output containers via the manager as below\r
  //==============================================================================\r
  mgr->ConnectInput  (task, 0, mgr->GetCommonInputContainer());\r
  mgr->ConnectOutput (task, 1, cout_pc);\r
  \r
  return task;\r
}\r
\r
\r