[u/mrichter/AliRoot.git] / STEER / AliQADataMakerSim.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/


/* $Id$ */

//
//  Base Class
//  Produces the data needed to calculate the quality assurance. 
//  All data must be mergeable objects.
//  Y. Schutz CERN July 2007
//

// --- ROOT system ---
#include <TFile.h>
#include <TTree.h>
#include <TClonesArray.h>

// --- Standard library ---

// --- AliRoot header files ---
#include "AliLog.h"
#include "AliQADataMakerSim.h"

ClassImp(AliQADataMakerSim)
             
//____________________________________________________________________________ 
AliQADataMakerSim::AliQADataMakerSim(const char * name, const char * title) : 
  AliQADataMaker(name, title), 
  fDigitsQAList(0x0), 
  fHitsQAList(0x0),
  fSDigitsQAList(0x0)
{
	// ctor
	fDetectorDirName = GetName() ; 
}

//____________________________________________________________________________ 
AliQADataMakerSim::AliQADataMakerSim(const AliQADataMakerSim& qadm) :
  AliQADataMaker(qadm.GetName(), qadm.GetTitle()), 
  fDigitsQAList(qadm.fDigitsQAList),
  fHitsQAList(qadm.fHitsQAList),
  fSDigitsQAList(qadm.fSDigitsQAList) 
{
  //copy ctor
  fDetectorDirName = GetName() ; 
}

//____________________________________________________________________________ 
AliQADataMakerSim::~AliQADataMakerSim()
{
	//dtor: delete the TObjArray and thei content
	if ( fDigitsQAList ) { 
		if ( fDigitsQAList->IsOwner() )
			fDigitsQAList->Delete() ;     
		delete fDigitsQAList ;     
	}
	if ( fHitsQAList ) {
		if ( fHitsQAList->IsOwner() ) 
			fHitsQAList->Delete() ;
		delete fHitsQAList ;
	}
	if ( fSDigitsQAList ) { 
		if ( fSDigitsQAList->IsOwner() ) 
			fSDigitsQAList->Delete() ; 
		delete fSDigitsQAList ; 
	}
}

//__________________________________________________________________
AliQADataMakerSim& AliQADataMakerSim::operator = (const AliQADataMakerSim& qadm )
{
  // Assignment operator.
  this->~AliQADataMakerSim();
  new(this) AliQADataMakerSim(qadm);
  return *this;
}

//____________________________________________________________________________
void AliQADataMakerSim::EndOfCycle() 
{ 
  // Finishes a cycle of QA for all tasks
  EndOfCycle(AliQA::kHITS) ; 
  EndOfCycle(AliQA::kSDIGITS) ; 
  EndOfCycle(AliQA::kDIGITS) ;
  ResetCycle() ; 
}

//____________________________________________________________________________
void AliQADataMakerSim::EndOfCycle(AliQA::TASKINDEX_t task) 
{ 
  // Finishes a cycle of QA data acquistion
	TObjArray * list = 0x0 ; 
	
	if ( task == AliQA::kHITS ) 
		list = fHitsQAList ; 
	else if ( task == AliQA::kSDIGITS )
		list = fSDigitsQAList ; 
	else if ( task == AliQA::kDIGITS ) 
		list = fDigitsQAList ; 
  
  if ( ! list ) 
    return ; 
	EndOfDetectorCycle(task, list) ; 
	TDirectory * subDir = fDetectorDir->GetDirectory(AliQA::GetTaskName(task)) ; 
	if (subDir) { 
		subDir->cd() ; 
    TIter next(list) ; 
    TH1 * obj ; 
    while ( (obj = dynamic_cast<TH1 *>(next())) ) {
      TString name(obj->GetTitle()) ;
      if (name.Contains(AliQA::GetExpert())) {
        name.ReplaceAll(AliQA::GetExpert(), "") ;
        obj->SetTitle(name) ; 
        obj->Write() ;
      }
    }
    if (WriteExpert()) {
      TDirectory * expertDir = subDir->GetDirectory(AliQA::GetExpert()) ; 
      if ( expertDir ) {
        expertDir->cd() ;
        next.Reset() ; 
        while ( (obj = dynamic_cast<TH1 *>(next())) ) {
          TString name(obj->GetTitle()) ;
          if (name.Contains(AliQA::GetExpert())) 
            continue ; 
          name.ReplaceAll(AliQA::GetExpert(), "") ;
          obj->SetTitle(name) ; 
          obj->Write() ;
        }      
      }
    }
  }
  fOutput->Save() ; 
  ResetCycle() ; 
}
 
//____________________________________________________________________________
void AliQADataMakerSim::Exec(AliQA::TASKINDEX_t task, TObject * data) 
{ 
  // creates the quality assurance data for the various tasks (Hits, SDigits, Digits, ESDs)
    
	if ( task == AliQA::kHITS ) {  
		AliDebug(1, "Processing Hits QA") ; 
		TClonesArray * arr = dynamic_cast<TClonesArray *>(data) ; 
		if (arr) { 
			MakeHits(arr) ;
		} else {
			TTree * tree = dynamic_cast<TTree *>(data) ; 
			if (tree) {
				MakeHits(tree) ; 
			} else {
				AliWarning("data are neither a TClonesArray nor a TTree") ; 
			}
		}
	} else if ( task == AliQA::kSDIGITS ) {
		AliDebug(1, "Processing SDigits QA") ; 
		TClonesArray * arr = dynamic_cast<TClonesArray *>(data) ; 
		if (arr) { 
			MakeSDigits(arr) ;
		} else {
			TTree * tree = dynamic_cast<TTree *>(data) ; 
			if (tree) {
				MakeSDigits(tree) ; 
			} else {
				AliWarning("data are neither a TClonesArray nor a TTree") ; 
			}
		}
	} else if ( task == AliQA::kDIGITS ) {
		AliDebug(1, "Processing Digits QA") ; 
		TClonesArray * arr = dynamic_cast<TClonesArray *>(data) ; 
		if (arr) { 
			MakeDigits(arr) ;
		} else {
			TTree * tree = dynamic_cast<TTree *>(data) ; 
			if (tree) {
				MakeDigits(tree) ; 
			} else {
				AliWarning("data are neither a TClonesArray nor a TTree") ; 
			}
		}
	}
}

//____________________________________________________________________________ 
TObjArray *  AliQADataMakerSim::Init(AliQA::TASKINDEX_t task, Int_t cycles)
{
  // general intialisation
	
	if (cycles > 0)
		SetCycle(cycles) ;  
	TObjArray * rv = NULL ; 
	if ( task == AliQA::kHITS ) {
		if ( ! fHitsQAList ) {
			fHitsQAList = new TObjArray(100) ;	 
   fHitsQAList->SetName(Form("%s/%s", GetName(), AliQA::GetTaskName(task).Data())) ; 
			InitHits() ;
		}
		rv = fHitsQAList ;
	} else if ( task == AliQA::kSDIGITS ) {
		if ( ! fSDigitsQAList ) {
			fSDigitsQAList = new TObjArray(100) ; 
   fSDigitsQAList->SetName(Form("%s/%s", GetName(), AliQA::GetTaskName(task).Data())) ; 
			InitSDigits() ;
		}
		rv = fSDigitsQAList ;
   } else if ( task == AliQA::kDIGITS ) {
	   if ( ! fDigitsQAList ) {
		   fDigitsQAList = new TObjArray(100) ;
     fDigitsQAList->SetName(Form("%s/%s", GetName(), AliQA::GetTaskName(task).Data())) ; 
		   InitDigits() ;
	   }
	   rv =  fDigitsQAList ;
   }
  
	return rv ; 
}

//____________________________________________________________________________ 
void AliQADataMakerSim::Init(AliQA::TASKINDEX_t task, TObjArray * list, Int_t run, Int_t cycles)
{
  // Intialisation by passing the list of QA data booked elsewhere
  
	fRun = run ;
	if (cycles > 0)
		SetCycle(cycles) ;  
	
	if ( task == AliQA::kHITS ) {
		fHitsQAList = list ;	 
	} else if ( task == AliQA::kSDIGITS) {
		fSDigitsQAList = list ; 
	} else if ( task == AliQA::kDIGITS ) {
		fDigitsQAList = list ; 
	} 
}

//____________________________________________________________________________
void AliQADataMakerSim::StartOfCycle(Int_t run) 
{ 
  // Finishes a cycle of QA for all tasks
  Bool_t samecycle = kFALSE ; 
  StartOfCycle(AliQA::kHITS,    run, samecycle) ;
  samecycle = kTRUE ; 
  StartOfCycle(AliQA::kSDIGITS, run, samecycle) ;
  StartOfCycle(AliQA::kDIGITS,  run, samecycle) ;
}

//____________________________________________________________________________
void AliQADataMakerSim::StartOfCycle(AliQA::TASKINDEX_t task, Int_t run, const Bool_t sameCycle) 
{ 
  // Finishes a cycle of QA data acquistion
  if ( run > 0 ) 
    fRun = run ; 
	if ( !sameCycle || fCurrentCycle == -1) {
		ResetCycle() ;
	if (fOutput) 
		fOutput->Close() ; 
	fOutput = AliQA::GetQADataFile(GetName(), fRun) ; 	
	}	

	AliInfo(Form(" Run %d Cycle %d task %s file %s", 
				 fRun, fCurrentCycle, AliQA::GetTaskName(task).Data(), fOutput->GetName() )) ;

	fDetectorDir = fOutput->GetDirectory(GetDetectorDirName()) ; 
	if (!fDetectorDir)
		fDetectorDir = fOutput->mkdir(GetDetectorDirName()) ; 

	TDirectory * subDir = fDetectorDir->GetDirectory(AliQA::GetTaskName(task)) ; 
	if (!subDir)
		subDir = fDetectorDir->mkdir(AliQA::GetTaskName(task)) ;  
  
  TDirectory * expertDir = subDir->GetDirectory(AliQA::GetExpert()) ; 
  if (!expertDir)
    expertDir = subDir->mkdir(AliQA::GetExpert()) ; 

	subDir->cd() ; 
	  
	StartOfDetectorCycle() ; 
}
Commit	Line	Data
04236e67	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
	16
	17	/* $Id$ */
	18
202374b1	19	//
	20	// Base Class
	21	// Produces the data needed to calculate the quality assurance.
	22	// All data must be mergeable objects.
	23	// Y. Schutz CERN July 2007
	24	//
04236e67	25
04236e67	26	// --- ROOT system ---
04236e67	27	#include <TFile.h>
04236e67	28	#include <TTree.h>
	29	#include <TClonesArray.h>
	30
	31	// --- Standard library ---
	32
	33	// --- AliRoot header files ---
	34	#include "AliLog.h"
	35	#include "AliQADataMakerSim.h"
04236e67	36
	37	ClassImp(AliQADataMakerSim)
	38
	39	//____________________________________________________________________________
	40	AliQADataMakerSim::AliQADataMakerSim(const char * name, const char * title) :
	41	AliQADataMaker(name, title),
	42	fDigitsQAList(0x0),
	43	fHitsQAList(0x0),
	44	fSDigitsQAList(0x0)
	45	{
	46	// ctor
	47	fDetectorDirName = GetName() ;
	48	}
	49
	50	//____________________________________________________________________________
	51	AliQADataMakerSim::AliQADataMakerSim(const AliQADataMakerSim& qadm) :
	52	AliQADataMaker(qadm.GetName(), qadm.GetTitle()),
	53	fDigitsQAList(qadm.fDigitsQAList),
	54	fHitsQAList(qadm.fHitsQAList),
	55	fSDigitsQAList(qadm.fSDigitsQAList)
	56	{
	57	//copy ctor
	58	fDetectorDirName = GetName() ;
	59	}
	60
63c6f8ae	61	//____________________________________________________________________________
	62	AliQADataMakerSim::~AliQADataMakerSim()
	63	{
	64	//dtor: delete the TObjArray and thei content
7ff8385d	65	if ( fDigitsQAList ) {
b3cc11cb	66	if ( fDigitsQAList->IsOwner() )
b3cc11cb	67	fDigitsQAList->Delete() ;
7ff8385d	68	delete fDigitsQAList ;
	69	}
	70	if ( fHitsQAList ) {
b3cc11cb	71	if ( fHitsQAList->IsOwner() )
b3cc11cb	72	fHitsQAList->Delete() ;
7ff8385d	73	delete fHitsQAList ;
	74	}
	75	if ( fSDigitsQAList ) {
b3cc11cb	76	if ( fSDigitsQAList->IsOwner() )
b3cc11cb	77	fSDigitsQAList->Delete() ;
7ff8385d	78	delete fSDigitsQAList ;
7ff8385d	79	}
63c6f8ae	80	}
63c6f8ae	81
04236e67	82	//__________________________________________________________________
	83	AliQADataMakerSim& AliQADataMakerSim::operator = (const AliQADataMakerSim& qadm )
	84	{
	85	// Assignment operator.
	86	this->~AliQADataMakerSim();
	87	new(this) AliQADataMakerSim(qadm);
	88	return *this;
	89	}
	90
930e6e3e	91	//____________________________________________________________________________
	92	void AliQADataMakerSim::EndOfCycle()
	93	{
	94	// Finishes a cycle of QA for all tasks
	95	EndOfCycle(AliQA::kHITS) ;
	96	EndOfCycle(AliQA::kSDIGITS) ;
	97	EndOfCycle(AliQA::kDIGITS) ;
	98	ResetCycle() ;
	99	}
	100
04236e67	101	//____________________________________________________________________________
96d67a8d	102	void AliQADataMakerSim::EndOfCycle(AliQA::TASKINDEX_t task)
04236e67	103	{
	104	// Finishes a cycle of QA data acquistion
	105	TObjArray * list = 0x0 ;
	106
	107	if ( task == AliQA::kHITS )
	108	list = fHitsQAList ;
	109	else if ( task == AliQA::kSDIGITS )
	110	list = fSDigitsQAList ;
	111	else if ( task == AliQA::kDIGITS )
	112	list = fDigitsQAList ;
930e6e3e	113
	114	if ( ! list )
	115	return ;
04236e67	116	EndOfDetectorCycle(task, list) ;
04236e67	117	TDirectory * subDir = fDetectorDir->GetDirectory(AliQA::GetTaskName(task)) ;
ccbf0759	118	if (subDir) {
ccbf0759	119	subDir->cd() ;
6e65319c	120	TIter next(list) ;
	121	TH1 * obj ;
	122	while ( (obj = dynamic_cast<TH1 *>(next())) ) {
	123	TString name(obj->GetTitle()) ;
	124	if (name.Contains(AliQA::GetExpert())) {
	125	name.ReplaceAll(AliQA::GetExpert(), "") ;
	126	obj->SetTitle(name) ;
	127	obj->Write() ;
	128	}
	129	}
b1af1125	130	if (WriteExpert()) {
	131	TDirectory * expertDir = subDir->GetDirectory(AliQA::GetExpert()) ;
	132	if ( expertDir ) {
	133	expertDir->cd() ;
	134	next.Reset() ;
	135	while ( (obj = dynamic_cast<TH1 *>(next())) ) {
	136	TString name(obj->GetTitle()) ;
	137	if (name.Contains(AliQA::GetExpert()))
	138	continue ;
	139	name.ReplaceAll(AliQA::GetExpert(), "") ;
	140	obj->SetTitle(name) ;
	141	obj->Write() ;
	142	}
	143	}
6e65319c	144	}
6e65319c	145	}
a595ac34	146	fOutput->Save() ;
930e6e3e	147	ResetCycle() ;
04236e67	148	}
	149
	150	//____________________________________________________________________________
96d67a8d	151	void AliQADataMakerSim::Exec(AliQA::TASKINDEX_t task, TObject * data)
04236e67	152	{
	153	// creates the quality assurance data for the various tasks (Hits, SDigits, Digits, ESDs)
	154
	155	if ( task == AliQA::kHITS ) {
	156	AliDebug(1, "Processing Hits QA") ;
	157	TClonesArray * arr = dynamic_cast<TClonesArray *>(data) ;
	158	if (arr) {
	159	MakeHits(arr) ;
	160	} else {
	161	TTree * tree = dynamic_cast<TTree *>(data) ;
	162	if (tree) {
	163	MakeHits(tree) ;
	164	} else {
	165	AliWarning("data are neither a TClonesArray nor a TTree") ;
	166	}
	167	}
	168	} else if ( task == AliQA::kSDIGITS ) {
	169	AliDebug(1, "Processing SDigits QA") ;
	170	TClonesArray * arr = dynamic_cast<TClonesArray *>(data) ;
	171	if (arr) {
	172	MakeSDigits(arr) ;
	173	} else {
	174	TTree * tree = dynamic_cast<TTree *>(data) ;
	175	if (tree) {
	176	MakeSDigits(tree) ;
	177	} else {
	178	AliWarning("data are neither a TClonesArray nor a TTree") ;
	179	}
	180	}
	181	} else if ( task == AliQA::kDIGITS ) {
	182	AliDebug(1, "Processing Digits QA") ;
	183	TClonesArray * arr = dynamic_cast<TClonesArray *>(data) ;
	184	if (arr) {
	185	MakeDigits(arr) ;
	186	} else {
	187	TTree * tree = dynamic_cast<TTree *>(data) ;
	188	if (tree) {
	189	MakeDigits(tree) ;
	190	} else {
	191	AliWarning("data are neither a TClonesArray nor a TTree") ;
	192	}
	193	}
	194	}
	195	}
	196
	197	//____________________________________________________________________________
930e6e3e	198	TObjArray * AliQADataMakerSim::Init(AliQA::TASKINDEX_t task, Int_t cycles)
04236e67	199	{
	200	// general intialisation
	201
04236e67	202	if (cycles > 0)
	203	SetCycle(cycles) ;
	204	TObjArray * rv = NULL ;
	205	if ( task == AliQA::kHITS ) {
63c6f8ae	206	if ( ! fHitsQAList ) {
63c6f8ae	207	fHitsQAList = new TObjArray(100) ;
a2b64fbd	208	fHitsQAList->SetName(Form("%s/%s", GetName(), AliQA::GetTaskName(task).Data())) ;
63c6f8ae	209	InitHits() ;
63c6f8ae	210	}
04236e67	211	rv = fHitsQAList ;
04236e67	212	} else if ( task == AliQA::kSDIGITS ) {
63c6f8ae	213	if ( ! fSDigitsQAList ) {
63c6f8ae	214	fSDigitsQAList = new TObjArray(100) ;
a2b64fbd	215	fSDigitsQAList->SetName(Form("%s/%s", GetName(), AliQA::GetTaskName(task).Data())) ;
63c6f8ae	216	InitSDigits() ;
63c6f8ae	217	}
04236e67	218	rv = fSDigitsQAList ;
04236e67	219	} else if ( task == AliQA::kDIGITS ) {
63c6f8ae	220	if ( ! fDigitsQAList ) {
63c6f8ae	221	fDigitsQAList = new TObjArray(100) ;
a2b64fbd	222	fDigitsQAList->SetName(Form("%s/%s", GetName(), AliQA::GetTaskName(task).Data())) ;
63c6f8ae	223	InitDigits() ;
63c6f8ae	224	}
04236e67	225	rv = fDigitsQAList ;
	226	}
	227
	228	return rv ;
	229	}
	230
	231	//____________________________________________________________________________
96d67a8d	232	void AliQADataMakerSim::Init(AliQA::TASKINDEX_t task, TObjArray * list, Int_t run, Int_t cycles)
04236e67	233	{
	234	// Intialisation by passing the list of QA data booked elsewhere
	235
	236	fRun = run ;
	237	if (cycles > 0)
	238	SetCycle(cycles) ;
	239
	240	if ( task == AliQA::kHITS ) {
	241	fHitsQAList = list ;
	242	} else if ( task == AliQA::kSDIGITS) {
	243	fSDigitsQAList = list ;
	244	} else if ( task == AliQA::kDIGITS ) {
	245	fDigitsQAList = list ;
	246	}
	247	}
	248
	249	//____________________________________________________________________________
930e6e3e	250	void AliQADataMakerSim::StartOfCycle(Int_t run)
	251	{
	252	// Finishes a cycle of QA for all tasks
	253	Bool_t samecycle = kFALSE ;
	254	StartOfCycle(AliQA::kHITS, run, samecycle) ;
	255	samecycle = kTRUE ;
	256	StartOfCycle(AliQA::kSDIGITS, run, samecycle) ;
	257	StartOfCycle(AliQA::kDIGITS, run, samecycle) ;
	258	}
	259
	260	//____________________________________________________________________________
	261	void AliQADataMakerSim::StartOfCycle(AliQA::TASKINDEX_t task, Int_t run, const Bool_t sameCycle)
04236e67	262	{
04236e67	263	// Finishes a cycle of QA data acquistion
930e6e3e	264	if ( run > 0 )
930e6e3e	265	fRun = run ;
04236e67	266	if ( !sameCycle \|\| fCurrentCycle == -1) {
	267	ResetCycle() ;
	268	if (fOutput)
	269	fOutput->Close() ;
930e6e3e	270	fOutput = AliQA::GetQADataFile(GetName(), fRun) ;
04236e67	271	}
	272
	273	AliInfo(Form(" Run %d Cycle %d task %s file %s",
	274	fRun, fCurrentCycle, AliQA::GetTaskName(task).Data(), fOutput->GetName() )) ;
	275
	276	fDetectorDir = fOutput->GetDirectory(GetDetectorDirName()) ;
	277	if (!fDetectorDir)
	278	fDetectorDir = fOutput->mkdir(GetDetectorDirName()) ;
	279
	280	TDirectory * subDir = fDetectorDir->GetDirectory(AliQA::GetTaskName(task)) ;
	281	if (!subDir)
	282	subDir = fDetectorDir->mkdir(AliQA::GetTaskName(task)) ;
6e65319c	283
	284	TDirectory * expertDir = subDir->GetDirectory(AliQA::GetExpert()) ;
	285	if (!expertDir)
	286	expertDir = subDir->mkdir(AliQA::GetExpert()) ;
	287
04236e67	288	subDir->cd() ;
930e6e3e	289
04236e67	290	StartOfDetectorCycle() ;
04236e67	291	}