[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::operator = (const AliQADataMakerSim& qadm )
{
  // Assignment operator.
  this->~AliQADataMakerSim();
  new(this) AliQADataMakerSim(qadm);
  return *this;
}

//____________________________________________________________________________
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 ; 
	
	EndOfDetectorCycle(task, list) ; 
	TDirectory * subDir = fDetectorDir->GetDirectory(AliQA::GetTaskName(task)) ; 
	if (subDir) { 
		subDir->cd() ; 
		list->Write() ; 
	}
}
 
//____________________________________________________________________________
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 run, Int_t cycles)
{
  // general intialisation
	
	fRun = run ;
	if (cycles > 0)
		SetCycle(cycles) ;  
	TObjArray * rv = NULL ; 
	if ( task == AliQA::kHITS ) {
		fHitsQAList = new TObjArray(100) ;	 
		InitHits() ;
		rv = fHitsQAList ;
	} else if ( task == AliQA::kSDIGITS ) {
		fSDigitsQAList = new TObjArray(100) ; 
		InitSDigits() ;
		rv = fSDigitsQAList ;
   } else if ( task == AliQA::kDIGITS ) {
	   fDigitsQAList = new TObjArray(100); 
	   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(AliQA::TASKINDEX_t task, const Bool_t sameCycle) 
{ 
  // Finishes a cycle of QA data acquistion
	if ( !sameCycle || fCurrentCycle == -1) {
		ResetCycle() ;
	if (fOutput) 
		fOutput->Close() ; 
	fOutput = AliQA::GetQADataFile(GetName(), fRun, fCurrentCycle) ; 	
	}	

	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)) ;  
	subDir->cd() ; 
	
	TObjArray * list = 0x0 ; 
  
	if ( task == AliQA::kHITS )  
		list = fHitsQAList ;
	else if ( task == AliQA::kSDIGITS )  
		list = fSDigitsQAList ;
	else  if ( task == AliQA::kDIGITS ) 
		list = fDigitsQAList ;
	
// Should be the choice of detectors
//	TIter next(list) ;
//	TH1 * h ; 
//	while ( (h = dynamic_cast<TH1 *>(next())) )
//		h->Reset() ;  
//
	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
	61	//__________________________________________________________________
	62	AliQADataMakerSim& AliQADataMakerSim::operator = (const AliQADataMakerSim& qadm )
	63	{
	64	// Assignment operator.
	65	this->~AliQADataMakerSim();
	66	new(this) AliQADataMakerSim(qadm);
	67	return *this;
	68	}
	69
	70	//____________________________________________________________________________
96d67a8d	71	void AliQADataMakerSim::EndOfCycle(AliQA::TASKINDEX_t task)
04236e67	72	{
	73	// Finishes a cycle of QA data acquistion
	74	TObjArray * list = 0x0 ;
	75
	76	if ( task == AliQA::kHITS )
	77	list = fHitsQAList ;
	78	else if ( task == AliQA::kSDIGITS )
	79	list = fSDigitsQAList ;
	80	else if ( task == AliQA::kDIGITS )
	81	list = fDigitsQAList ;
	82
	83	EndOfDetectorCycle(task, list) ;
	84	TDirectory * subDir = fDetectorDir->GetDirectory(AliQA::GetTaskName(task)) ;
ccbf0759	85	if (subDir) {
	86	subDir->cd() ;
	87	list->Write() ;
	88	}
04236e67	89	}
	90
	91	//____________________________________________________________________________
96d67a8d	92	void AliQADataMakerSim::Exec(AliQA::TASKINDEX_t task, TObject * data)
04236e67	93	{
	94	// creates the quality assurance data for the various tasks (Hits, SDigits, Digits, ESDs)
	95
	96	if ( task == AliQA::kHITS ) {
	97	AliDebug(1, "Processing Hits QA") ;
	98	TClonesArray * arr = dynamic_cast<TClonesArray *>(data) ;
	99	if (arr) {
	100	MakeHits(arr) ;
	101	} else {
	102	TTree * tree = dynamic_cast<TTree *>(data) ;
	103	if (tree) {
	104	MakeHits(tree) ;
	105	} else {
	106	AliWarning("data are neither a TClonesArray nor a TTree") ;
	107	}
	108	}
	109	} else if ( task == AliQA::kSDIGITS ) {
	110	AliDebug(1, "Processing SDigits QA") ;
	111	TClonesArray * arr = dynamic_cast<TClonesArray *>(data) ;
	112	if (arr) {
	113	MakeSDigits(arr) ;
	114	} else {
	115	TTree * tree = dynamic_cast<TTree *>(data) ;
	116	if (tree) {
	117	MakeSDigits(tree) ;
	118	} else {
	119	AliWarning("data are neither a TClonesArray nor a TTree") ;
	120	}
	121	}
	122	} else if ( task == AliQA::kDIGITS ) {
	123	AliDebug(1, "Processing Digits QA") ;
	124	TClonesArray * arr = dynamic_cast<TClonesArray *>(data) ;
	125	if (arr) {
	126	MakeDigits(arr) ;
	127	} else {
	128	TTree * tree = dynamic_cast<TTree *>(data) ;
	129	if (tree) {
	130	MakeDigits(tree) ;
	131	} else {
	132	AliWarning("data are neither a TClonesArray nor a TTree") ;
	133	}
	134	}
	135	}
	136	}
	137
	138	//____________________________________________________________________________
96d67a8d	139	TObjArray * AliQADataMakerSim::Init(AliQA::TASKINDEX_t task, Int_t run, Int_t cycles)
04236e67	140	{
	141	// general intialisation
	142
	143	fRun = run ;
	144	if (cycles > 0)
	145	SetCycle(cycles) ;
	146	TObjArray * rv = NULL ;
	147	if ( task == AliQA::kHITS ) {
	148	fHitsQAList = new TObjArray(100) ;
	149	InitHits() ;
	150	rv = fHitsQAList ;
	151	} else if ( task == AliQA::kSDIGITS ) {
	152	fSDigitsQAList = new TObjArray(100) ;
	153	InitSDigits() ;
	154	rv = fSDigitsQAList ;
	155	} else if ( task == AliQA::kDIGITS ) {
	156	fDigitsQAList = new TObjArray(100);
	157	InitDigits() ;
	158	rv = fDigitsQAList ;
	159	}
	160
	161	return rv ;
	162	}
	163
	164	//____________________________________________________________________________
96d67a8d	165	void AliQADataMakerSim::Init(AliQA::TASKINDEX_t task, TObjArray * list, Int_t run, Int_t cycles)
04236e67	166	{
	167	// Intialisation by passing the list of QA data booked elsewhere
	168
	169	fRun = run ;
	170	if (cycles > 0)
	171	SetCycle(cycles) ;
	172
	173	if ( task == AliQA::kHITS ) {
	174	fHitsQAList = list ;
	175	} else if ( task == AliQA::kSDIGITS) {
	176	fSDigitsQAList = list ;
	177	} else if ( task == AliQA::kDIGITS ) {
	178	fDigitsQAList = list ;
	179	}
	180	}
	181
	182	//____________________________________________________________________________
96d67a8d	183	void AliQADataMakerSim::StartOfCycle(AliQA::TASKINDEX_t task, const Bool_t sameCycle)
04236e67	184	{
	185	// Finishes a cycle of QA data acquistion
	186	if ( !sameCycle \|\| fCurrentCycle == -1) {
	187	ResetCycle() ;
	188	if (fOutput)
	189	fOutput->Close() ;
	190	fOutput = AliQA::GetQADataFile(GetName(), fRun, fCurrentCycle) ;
	191	}
	192
	193	AliInfo(Form(" Run %d Cycle %d task %s file %s",
	194	fRun, fCurrentCycle, AliQA::GetTaskName(task).Data(), fOutput->GetName() )) ;
	195
	196	fDetectorDir = fOutput->GetDirectory(GetDetectorDirName()) ;
	197	if (!fDetectorDir)
	198	fDetectorDir = fOutput->mkdir(GetDetectorDirName()) ;
	199
	200	TDirectory * subDir = fDetectorDir->GetDirectory(AliQA::GetTaskName(task)) ;
	201	if (!subDir)
	202	subDir = fDetectorDir->mkdir(AliQA::GetTaskName(task)) ;
	203	subDir->cd() ;
	204
	205	TObjArray * list = 0x0 ;
	206
	207	if ( task == AliQA::kHITS )
	208	list = fHitsQAList ;
	209	else if ( task == AliQA::kSDIGITS )
	210	list = fSDigitsQAList ;
	211	else if ( task == AliQA::kDIGITS )
	212	list = fDigitsQAList ;
075f5ecf	213
	214	// Should be the choice of detectors
	215	// TIter next(list) ;
	216	// TH1 * h ;
	217	// while ( (h = dynamic_cast<TH1 *>(next())) )
	218	// h->Reset() ;
	219	//
04236e67	220	StartOfDetectorCycle() ;
04236e67	221	}