[u/mrichter/AliRoot.git] / ITS / AliITSQADataMakerSim.cxx

/**************************************************************************
 * Copyright(c) 2007-2009, 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$ */

//  *************************************************************
//  Checks the quality assurance 
//  by comparing with reference data
//  contained in a DB
//  -------------------------------------------------------------
//  W. Ferrarese + P. Cerello Feb 2008
//  INFN Torino

// --- ROOT system ---
#include <TTree.h>
#include <TMath.h>
// --- Standard library ---

// --- AliRoot header files ---
#include "AliITSQADataMakerSim.h"
#include "AliITSQASPDDataMakerSim.h"
#include "AliITSQASDDDataMakerSim.h"
#include "AliITSQASSDDataMakerSim.h"
#include "AliLog.h"
#include "AliQAv1.h"
#include "AliQAChecker.h"
#include "AliITSQAChecker.h"
#include "AliRawReader.h"

ClassImp(AliITSQADataMakerSim)

//____________________________________________________________________________ 
AliITSQADataMakerSim::AliITSQADataMakerSim(Short_t subDet) :
AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kITS), "ITS Quality Assurance Data Maker"),
fSubDetector(subDet),
fSPDDataMaker(NULL),
fSDDDataMaker(NULL),
fSSDDataMaker(NULL)
{
  //ctor used to discriminate OnLine-Offline analysis
  if(fSubDetector < 0 || fSubDetector > 3) {
	AliError("Error: fSubDetector number out of range; return\n");
  }

  // Initialization for each subdetector 
  if(fSubDetector == 0 || fSubDetector == 1) {
    AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SPD DataMakerSim\n");
	fSPDDataMaker = new AliITSQASPDDataMakerSim(this);
  }
  if(fSubDetector == 0 || fSubDetector == 2) {
    AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SDD DataMakerSim\n");
	//printf("AliITSQADM::Create SDD DataMakerSim\n");		    
	fSDDDataMaker = new AliITSQASDDDataMakerSim(this);
  }
  if(fSubDetector == 0 || fSubDetector == 3) {
    AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SSD DataMakerSim\n");
	fSSDDataMaker = new AliITSQASSDDataMakerSim(this);
  }
}

//____________________________________________________________________________ 
AliITSQADataMakerSim::~AliITSQADataMakerSim(){
  // destructor
  if(fSPDDataMaker)delete fSPDDataMaker;
  if(fSDDDataMaker)delete fSDDDataMaker;
  if(fSSDDataMaker)delete fSSDDataMaker;
}

//____________________________________________________________________________ 
AliITSQADataMakerSim::AliITSQADataMakerSim(const AliITSQADataMakerSim& qadm) :
AliQADataMakerSim(),
fSubDetector(qadm.fSubDetector),
fSPDDataMaker(NULL),
fSDDDataMaker(NULL),
fSSDDataMaker(NULL)
{
  //copy ctor 
  SetName((const char*)qadm.GetName()) ; 
  SetTitle((const char*)qadm.GetTitle());
}

//__________________________________________________________________
AliITSQADataMakerSim& AliITSQADataMakerSim::operator = (const AliITSQADataMakerSim& qac )
{
  // Equal operator.
  this->~AliITSQADataMakerSim();
  new(this) AliITSQADataMakerSim(qac);
  return *this;
}

//____________________________________________________________________________ 
void AliITSQADataMakerSim::StartOfDetectorCycle()
{
  //Detector specific actions at start of cycle
  AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Start of ITS Cycle\n");
  ResetEventTrigClasses(); // reset triggers list to select all histos
  ResetEvCountCycle();
  //  
  if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->StartOfDetectorCycle();
  if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->StartOfDetectorCycle();
  if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->StartOfDetectorCycle();
}

//____________________________________________________________________________ 
void AliITSQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray** list)
{
  //
  AliInfo(Form("End of Dedetctor Cycle called for %s\n",AliQAv1::GetTaskName(task).Data() ));
  ResetEventTrigClasses(); // reset triggers list to select all histos
  //
  // launch the QA checking
  for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
    if (! IsValidEventSpecie(specie, list) ) continue; 
    SetEventSpecie(AliRecoParam::ConvertIndex(specie));
    //
    Int_t idnumber=list[specie]->GetUniqueID();
    //printf("specie %s \t id number == %d\n",AliRecoParam::GetEventSpecieName(specie),idnumber);
    if(idnumber==0) {
      //AliInfo(Form("No check for %s\n",AliQAv1::GetTaskName(task).Data() ))
      continue;
    } //skip kDigitsR and not filled TobjArray specie
    else{
      AliDebug(AliQAv1::GetQADebugLevel(),"AliITSDM instantiates checker with Run(AliQAv1::kITS, task, list)\n"); 
      if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->EndOfDetectorCycle(task, list); //[specie]);
      if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->EndOfDetectorCycle(task, list); //[specie]);
      if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->EndOfDetectorCycle(task, list); //[specie]);
      //
      AliQAChecker *qac = AliQAChecker::Instance();
      AliITSQAChecker *qacb = (AliITSQAChecker *) qac->GetDetQAChecker(0);
      Int_t subdet=GetSubDet();
      qacb->SetSubDet(subdet);
      
      if(subdet== 0 ){
	qacb->SetTaskOffset(fSPDDataMaker->GetOffset(task,specie),fSDDDataMaker->GetOffset(task,specie),fSSDDataMaker->GetOffset(task,specie)); //Setting the offset for the QAChecker list
	qacb->SetHisto(fSPDDataMaker->GetTaskHisto(task), fSDDDataMaker->GetTaskHisto(task), fSSDDataMaker->GetTaskHisto(task));	
      }
      else
	if(subdet!=0){
	  Int_t offset=GetDetTaskOffset(subdet, task);
	  qacb->SetDetTaskOffset(subdet,offset);
	  Int_t histo=GetDetTaskHisto(subdet, task);
	  qacb->SetDetHisto(subdet,histo);
	}
      qac->Run( AliQAv1::kITS , task, list);
    }//end else unique id 
  }//end for
}

//____________________________________________________________________________ 
void AliITSQADataMakerSim::InitDigits()
{  

 // Initialization for Digits data 
  fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(60);
  if(fSubDetector == 0 || fSubDetector == 1) {
    AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitDigits\n");
    
    fSPDDataMaker->InitDigits();
  }
  if(fSubDetector == 0 || fSubDetector == 2) {
    AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitDigits\n");
    
    fSDDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
    fSDDDataMaker->InitDigits();
  }
  if(fSubDetector == 0 || fSubDetector == 3) {
    AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitDigits\n");
    
    fSSDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
    fSSDDataMaker->InitDigits();
  }
  //
  ClonePerTrigClass(AliQAv1::kDIGITS); // this should be the last line
}

//____________________________________________________________________________
void AliITSQADataMakerSim::MakeDigits()
{ 
  // Fill QA for digits   
  if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeDigits(); 
  if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeDigits(); 
  if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeDigits();
  //
  IncEvCountCycleDigits();
  IncEvCountTotalDigits();
  //
}

//____________________________________________________________________________
void AliITSQADataMakerSim::MakeDigits(TTree * digits)
{ 
  // Fill QA for digits   
  if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeDigits(digits); 
  if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeDigits(digits); 
  if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeDigits(digits);
  //
  IncEvCountCycleDigits();
  IncEvCountTotalDigits();
  //
}

//____________________________________________________________________________ 
void AliITSQADataMakerSim::InitSDigits()
{
  // Initialization for SDigits
  fSDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(70);
  if(fSubDetector == 0 || fSubDetector == 1) {
    AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitSDigits\n");

    fSPDDataMaker->InitSDigits();
  }
  if(fSubDetector == 0 || fSubDetector == 2){ 
    AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitSDigits\n");

    fSDDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
    fSDDDataMaker->InitSDigits();
  }
  if(fSubDetector == 0 || fSubDetector == 3) {
    AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitSDigits\n");

    fSSDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
    fSSDDataMaker->InitSDigits();
  }
  //
  ClonePerTrigClass(AliQAv1::kSDIGITS); // this should be the last line
}

//____________________________________________________________________________ 
void AliITSQADataMakerSim::MakeSDigits()
{
  // Fill QA for sdigits
  if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeSDigits(); 
  if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeSDigits(); 
  if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeSDigits();
  //
  IncEvCountCycleSDigits();
  IncEvCountTotalSDigits();
  //
}

//____________________________________________________________________________ 
void AliITSQADataMakerSim::MakeSDigits(TTree * sdigits)
{
  // Fill QA for recpoints
  if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeSDigits(sdigits); 
  if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeSDigits(sdigits); 
  if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeSDigits(sdigits);
  //
  IncEvCountCycleSDigits();
  IncEvCountTotalSDigits();
  //
}

//____________________________________________________________________________ 
void AliITSQADataMakerSim::InitHits()
{
  // Initialization for hits
  fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(50);
  if(fSubDetector == 0 || fSubDetector == 1) {
    AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitHits\n");
    fSPDDataMaker->InitHits();
  }
  if(fSubDetector == 0 || fSubDetector == 2) {
    AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitHits\n");

    fSDDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
    fSDDDataMaker->InitHits();
  }
  if(fSubDetector == 0 || fSubDetector == 3) {
    AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitHits\n");

    fSSDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
    fSSDDataMaker->InitHits();
  }
  //
  ClonePerTrigClass(AliQAv1::kHITS); // this should be the last line
}

//____________________________________________________________________________ 
void AliITSQADataMakerSim::MakeHits()
{
  // Fill QA for hits
  if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeHits(); 
  if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeHits(); 
  if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeHits();
  //
  IncEvCountCycleHits();
  IncEvCountTotalHits();
  //
}

//____________________________________________________________________________ 
void AliITSQADataMakerSim::MakeHits(TTree * hits)
{
  // Fill QA for hits
  if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeHits(hits);
  if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeHits(hits);
  if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeHits(hits);
  //
  IncEvCountCycleHits();
  IncEvCountTotalHits();
  //
}

//_________________________________________________________________
Int_t AliITSQADataMakerSim::GetDetTaskOffset(Int_t subdet,AliQAv1::TASKINDEX_t task)
{

  //return the offset for each subdetector
  switch(subdet)
    {
      Int_t offset;
    case 1:
      offset=fSPDDataMaker->GetOffset(task);
      return offset;
      break;
    case 2:
      offset=fSDDDataMaker->GetOffset(task);
      return offset;
      break;
    case 3:
      offset=fSSDDataMaker->GetOffset(task);
      return offset;
      break;
    default:
      AliWarning("No specific subdetector (SPD, SDD, SSD) selected!! Offset set to zero \n");
      offset=0;
      return offset;
      break;
    }
  //return offset;
}

//_________________________________________________________________
Int_t AliITSQADataMakerSim::GetDetTaskHisto(Int_t subdet,AliQAv1::TASKINDEX_t task)
{
  //return of the number of histograms for each task and for each sub detector
  switch(subdet)
    {

      Int_t histo;
    case 1:
      histo=fSPDDataMaker->GetOffset(task);
      return histo;
      break;
    case 2:
      histo=fSDDDataMaker->GetOffset(task);
      return histo;
      break;
    case 3:
      histo=fSSDDataMaker->GetOffset(task);
      return histo;
      break;
    default:
      AliWarning("No specific subdetector (SPD, SDD, SSD) selected!! Offset set to zero \n");
      histo=0;
      return histo;
      break;
    }
  //return offset;
}
Commit	Line	Data
	1	/**************************************************************************
	2	* Copyright(c) 2007-2009, 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	/* $Id$ */
	17
	18	// *************************************************************
	19	// Checks the quality assurance
	20	// by comparing with reference data
	21	// contained in a DB
	22	// -------------------------------------------------------------
	23	// W. Ferrarese + P. Cerello Feb 2008
	24	// INFN Torino
	25
	26	// --- ROOT system ---
	27	#include <TTree.h>
	28	#include <TMath.h>
	29	// --- Standard library ---
	30
	31	// --- AliRoot header files ---
	32	#include "AliITSQADataMakerSim.h"
	33	#include "AliITSQASPDDataMakerSim.h"
	34	#include "AliITSQASDDDataMakerSim.h"
	35	#include "AliITSQASSDDataMakerSim.h"
	36	#include "AliLog.h"
	37	#include "AliQAv1.h"
	38	#include "AliQAChecker.h"
	39	#include "AliITSQAChecker.h"
	40	#include "AliRawReader.h"
	41
	42	ClassImp(AliITSQADataMakerSim)
	43
	44	//____________________________________________________________________________
	45	AliITSQADataMakerSim::AliITSQADataMakerSim(Short_t subDet) :
	46	AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kITS), "ITS Quality Assurance Data Maker"),
	47	fSubDetector(subDet),
	48	fSPDDataMaker(NULL),
	49	fSDDDataMaker(NULL),
	50	fSSDDataMaker(NULL)
	51	{
	52	//ctor used to discriminate OnLine-Offline analysis
	53	if(fSubDetector < 0 \|\| fSubDetector > 3) {
	54	AliError("Error: fSubDetector number out of range; return\n");
	55	}
	56
	57	// Initialization for each subdetector
	58	if(fSubDetector == 0 \|\| fSubDetector == 1) {
	59	AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SPD DataMakerSim\n");
	60	fSPDDataMaker = new AliITSQASPDDataMakerSim(this);
	61	}
	62	if(fSubDetector == 0 \|\| fSubDetector == 2) {
	63	AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SDD DataMakerSim\n");
	64	//printf("AliITSQADM::Create SDD DataMakerSim\n");
	65	fSDDDataMaker = new AliITSQASDDDataMakerSim(this);
	66	}
	67	if(fSubDetector == 0 \|\| fSubDetector == 3) {
	68	AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SSD DataMakerSim\n");
	69	fSSDDataMaker = new AliITSQASSDDataMakerSim(this);
	70	}
	71	}
	72
	73	//____________________________________________________________________________
	74	AliITSQADataMakerSim::~AliITSQADataMakerSim(){
	75	// destructor
	76	if(fSPDDataMaker)delete fSPDDataMaker;
	77	if(fSDDDataMaker)delete fSDDDataMaker;
	78	if(fSSDDataMaker)delete fSSDDataMaker;
	79	}
	80
	81	//____________________________________________________________________________
	82	AliITSQADataMakerSim::AliITSQADataMakerSim(const AliITSQADataMakerSim& qadm) :
	83	AliQADataMakerSim(),
	84	fSubDetector(qadm.fSubDetector),
	85	fSPDDataMaker(NULL),
	86	fSDDDataMaker(NULL),
	87	fSSDDataMaker(NULL)
	88	{
	89	//copy ctor
	90	SetName((const char*)qadm.GetName()) ;
	91	SetTitle((const char*)qadm.GetTitle());
	92	}
	93
	94	//__________________________________________________________________
	95	AliITSQADataMakerSim& AliITSQADataMakerSim::operator = (const AliITSQADataMakerSim& qac )
	96	{
	97	// Equal operator.
	98	this->~AliITSQADataMakerSim();
	99	new(this) AliITSQADataMakerSim(qac);
	100	return *this;
	101	}
	102
	103	//____________________________________________________________________________
	104	void AliITSQADataMakerSim::StartOfDetectorCycle()
	105	{
	106	//Detector specific actions at start of cycle
	107	AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Start of ITS Cycle\n");
	108	ResetEventTrigClasses(); // reset triggers list to select all histos
	109	ResetEvCountCycle();
	110	//
	111	if(fSubDetector == 0 \|\| fSubDetector == 1) fSPDDataMaker->StartOfDetectorCycle();
	112	if(fSubDetector == 0 \|\| fSubDetector == 2) fSDDDataMaker->StartOfDetectorCycle();
	113	if(fSubDetector == 0 \|\| fSubDetector == 3) fSSDDataMaker->StartOfDetectorCycle();
	114	}
	115
	116	//____________________________________________________________________________
	117	void AliITSQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray** list)
	118	{
	119	//
	120	AliInfo(Form("End of Dedetctor Cycle called for %s\n",AliQAv1::GetTaskName(task).Data() ));
	121	ResetEventTrigClasses(); // reset triggers list to select all histos
	122	//
	123	// launch the QA checking
	124	for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
	125	if (! IsValidEventSpecie(specie, list) ) continue;
	126	SetEventSpecie(AliRecoParam::ConvertIndex(specie));
	127	//
	128	Int_t idnumber=list[specie]->GetUniqueID();
	129	//printf("specie %s \t id number == %d\n",AliRecoParam::GetEventSpecieName(specie),idnumber);
	130	if(idnumber==0) {
	131	//AliInfo(Form("No check for %s\n",AliQAv1::GetTaskName(task).Data() ))
	132	continue;
	133	} //skip kDigitsR and not filled TobjArray specie
	134	else{
	135	AliDebug(AliQAv1::GetQADebugLevel(),"AliITSDM instantiates checker with Run(AliQAv1::kITS, task, list)\n");
	136	if(fSubDetector == 0 \|\| fSubDetector == 1) fSPDDataMaker->EndOfDetectorCycle(task, list); //[specie]);
	137	if(fSubDetector == 0 \|\| fSubDetector == 2) fSDDDataMaker->EndOfDetectorCycle(task, list); //[specie]);
	138	if(fSubDetector == 0 \|\| fSubDetector == 3) fSSDDataMaker->EndOfDetectorCycle(task, list); //[specie]);
	139	//
	140	AliQAChecker *qac = AliQAChecker::Instance();
	141	AliITSQAChecker qacb = (AliITSQAChecker ) qac->GetDetQAChecker(0);
	142	Int_t subdet=GetSubDet();
	143	qacb->SetSubDet(subdet);
	144
	145	if(subdet== 0 ){
	146	qacb->SetTaskOffset(fSPDDataMaker->GetOffset(task,specie),fSDDDataMaker->GetOffset(task,specie),fSSDDataMaker->GetOffset(task,specie)); //Setting the offset for the QAChecker list
	147	qacb->SetHisto(fSPDDataMaker->GetTaskHisto(task), fSDDDataMaker->GetTaskHisto(task), fSSDDataMaker->GetTaskHisto(task));
	148	}
	149	else
	150	if(subdet!=0){
	151	Int_t offset=GetDetTaskOffset(subdet, task);
	152	qacb->SetDetTaskOffset(subdet,offset);
	153	Int_t histo=GetDetTaskHisto(subdet, task);
	154	qacb->SetDetHisto(subdet,histo);
	155	}
	156	qac->Run( AliQAv1::kITS , task, list);
	157	}//end else unique id
	158	}//end for
	159	}
	160
	161	//____________________________________________________________________________
	162	void AliITSQADataMakerSim::InitDigits()
	163	{
	164
	165	// Initialization for Digits data
	166	fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(60);
	167	if(fSubDetector == 0 \|\| fSubDetector == 1) {
	168	AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitDigits\n");
	169
	170	fSPDDataMaker->InitDigits();
	171	}
	172	if(fSubDetector == 0 \|\| fSubDetector == 2) {
	173	AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitDigits\n");
	174
	175	fSDDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
	176	fSDDDataMaker->InitDigits();
	177	}
	178	if(fSubDetector == 0 \|\| fSubDetector == 3) {
	179	AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitDigits\n");
	180
	181	fSSDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
	182	fSSDDataMaker->InitDigits();
	183	}
	184	//
	185	ClonePerTrigClass(AliQAv1::kDIGITS); // this should be the last line
	186	}
	187
	188	//____________________________________________________________________________
	189	void AliITSQADataMakerSim::MakeDigits()
	190	{
	191	// Fill QA for digits
	192	if(fSubDetector == 0 \|\| fSubDetector == 1) fSPDDataMaker->MakeDigits();
	193	if(fSubDetector == 0 \|\| fSubDetector == 2) fSDDDataMaker->MakeDigits();
	194	if(fSubDetector == 0 \|\| fSubDetector == 3) fSSDDataMaker->MakeDigits();
	195	//
	196	IncEvCountCycleDigits();
	197	IncEvCountTotalDigits();
	198	//
	199	}
	200
	201	//____________________________________________________________________________
	202	void AliITSQADataMakerSim::MakeDigits(TTree * digits)
	203	{
	204	// Fill QA for digits
	205	if(fSubDetector == 0 \|\| fSubDetector == 1) fSPDDataMaker->MakeDigits(digits);
	206	if(fSubDetector == 0 \|\| fSubDetector == 2) fSDDDataMaker->MakeDigits(digits);
	207	if(fSubDetector == 0 \|\| fSubDetector == 3) fSSDDataMaker->MakeDigits(digits);
	208	//
	209	IncEvCountCycleDigits();
	210	IncEvCountTotalDigits();
	211	//
	212	}
	213
	214	//____________________________________________________________________________
	215	void AliITSQADataMakerSim::InitSDigits()
	216	{
	217	// Initialization for SDigits
	218	fSDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(70);
	219	if(fSubDetector == 0 \|\| fSubDetector == 1) {
	220	AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitSDigits\n");
	221
	222	fSPDDataMaker->InitSDigits();
	223	}
	224	if(fSubDetector == 0 \|\| fSubDetector == 2){
	225	AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitSDigits\n");
	226
	227	fSDDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
	228	fSDDDataMaker->InitSDigits();
	229	}
	230	if(fSubDetector == 0 \|\| fSubDetector == 3) {
	231	AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitSDigits\n");
	232
	233	fSSDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
	234	fSSDDataMaker->InitSDigits();
	235	}
	236	//
	237	ClonePerTrigClass(AliQAv1::kSDIGITS); // this should be the last line
	238	}
	239
	240	//____________________________________________________________________________
	241	void AliITSQADataMakerSim::MakeSDigits()
	242	{
	243	// Fill QA for sdigits
	244	if(fSubDetector == 0 \|\| fSubDetector == 1) fSPDDataMaker->MakeSDigits();
	245	if(fSubDetector == 0 \|\| fSubDetector == 2) fSDDDataMaker->MakeSDigits();
	246	if(fSubDetector == 0 \|\| fSubDetector == 3) fSSDDataMaker->MakeSDigits();
	247	//
	248	IncEvCountCycleSDigits();
	249	IncEvCountTotalSDigits();
	250	//
	251	}
	252
	253	//____________________________________________________________________________
	254	void AliITSQADataMakerSim::MakeSDigits(TTree * sdigits)
	255	{
	256	// Fill QA for recpoints
	257	if(fSubDetector == 0 \|\| fSubDetector == 1) fSPDDataMaker->MakeSDigits(sdigits);
	258	if(fSubDetector == 0 \|\| fSubDetector == 2) fSDDDataMaker->MakeSDigits(sdigits);
	259	if(fSubDetector == 0 \|\| fSubDetector == 3) fSSDDataMaker->MakeSDigits(sdigits);
	260	//
	261	IncEvCountCycleSDigits();
	262	IncEvCountTotalSDigits();
	263	//
	264	}
	265
	266	//____________________________________________________________________________
	267	void AliITSQADataMakerSim::InitHits()
	268	{
	269	// Initialization for hits
	270	fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(50);
	271	if(fSubDetector == 0 \|\| fSubDetector == 1) {
	272	AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitHits\n");
	273	fSPDDataMaker->InitHits();
	274	}
	275	if(fSubDetector == 0 \|\| fSubDetector == 2) {
	276	AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitHits\n");
	277
	278	fSDDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
	279	fSDDDataMaker->InitHits();
	280	}
	281	if(fSubDetector == 0 \|\| fSubDetector == 3) {
	282	AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitHits\n");
	283
	284	fSSDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
	285	fSSDDataMaker->InitHits();
	286	}
	287	//
	288	ClonePerTrigClass(AliQAv1::kHITS); // this should be the last line
	289	}
	290
	291	//____________________________________________________________________________
	292	void AliITSQADataMakerSim::MakeHits()
	293	{
	294	// Fill QA for hits
	295	if(fSubDetector == 0 \|\| fSubDetector == 1) fSPDDataMaker->MakeHits();
	296	if(fSubDetector == 0 \|\| fSubDetector == 2) fSDDDataMaker->MakeHits();
	297	if(fSubDetector == 0 \|\| fSubDetector == 3) fSSDDataMaker->MakeHits();
	298	//
	299	IncEvCountCycleHits();
	300	IncEvCountTotalHits();
	301	//
	302	}
	303
	304	//____________________________________________________________________________
	305	void AliITSQADataMakerSim::MakeHits(TTree * hits)
	306	{
	307	// Fill QA for hits
	308	if(fSubDetector == 0 \|\| fSubDetector == 1) fSPDDataMaker->MakeHits(hits);
	309	if(fSubDetector == 0 \|\| fSubDetector == 2) fSDDDataMaker->MakeHits(hits);
	310	if(fSubDetector == 0 \|\| fSubDetector == 3) fSSDDataMaker->MakeHits(hits);
	311	//
	312	IncEvCountCycleHits();
	313	IncEvCountTotalHits();
	314	//
	315	}
	316
	317	//_________________________________________________________________
	318	Int_t AliITSQADataMakerSim::GetDetTaskOffset(Int_t subdet,AliQAv1::TASKINDEX_t task)
	319	{
	320
	321	//return the offset for each subdetector
	322	switch(subdet)
	323	{
	324	Int_t offset;
	325	case 1:
	326	offset=fSPDDataMaker->GetOffset(task);
	327	return offset;
	328	break;
	329	case 2:
	330	offset=fSDDDataMaker->GetOffset(task);
	331	return offset;
	332	break;
	333	case 3:
	334	offset=fSSDDataMaker->GetOffset(task);
	335	return offset;
	336	break;
	337	default:
	338	AliWarning("No specific subdetector (SPD, SDD, SSD) selected!! Offset set to zero \n");
	339	offset=0;
	340	return offset;
	341	break;
	342	}
	343	//return offset;
	344	}
	345
	346	//_________________________________________________________________
	347	Int_t AliITSQADataMakerSim::GetDetTaskHisto(Int_t subdet,AliQAv1::TASKINDEX_t task)
	348	{
	349	//return of the number of histograms for each task and for each sub detector
	350	switch(subdet)
	351	{
	352
	353	Int_t histo;
	354	case 1:
	355	histo=fSPDDataMaker->GetOffset(task);
	356	return histo;
	357	break;
	358	case 2:
	359	histo=fSDDDataMaker->GetOffset(task);
	360	return histo;
	361	break;
	362	case 3:
	363	histo=fSSDDataMaker->GetOffset(task);
	364	return histo;
	365	break;
	366	default:
	367	AliWarning("No specific subdetector (SPD, SDD, SSD) selected!! Offset set to zero \n");
	368	histo=0;
	369	return histo;
	370	break;
	371	}
	372	//return offset;
	373	}