[u/mrichter/AliRoot.git] / ITS / AliITSQASPDDataMakerSim.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 INFN Torino Feb 2008
//  M. Nicassio D. Elia INFN Bari April 2008
//  maria.nicassio@ba.infn.it

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

// --- AliRoot header files ---
#include "AliRun.h"
#include "AliITSQADataMakerSim.h"
#include "AliITSQASPDDataMakerSim.h"
#include "AliQA.h"
#include "AliQAChecker.h"
#include "AliITSdigit.h"   
#include "AliITSdigitSPD.h"
#include "AliITS.h"
#include "AliITSmodule.h"
#include "AliITShit.h"
#include "AliITSLoader.h"
#include "AliRunLoader.h"

ClassImp(AliITSQASPDDataMakerSim)

//____________________________________________________________________________ 
AliITSQASPDDataMakerSim::AliITSQASPDDataMakerSim(AliITSQADataMakerSim *aliITSQADataMakerSim) :
TObject(),
fAliITSQADataMakerSim(aliITSQADataMakerSim),
fSPDhHTask(0),
fSPDhSTask(0),
fSPDhDTask(0),
fGenOffsetH(0),
fGenOffsetS(0),
fGenOffsetD(0)
{
  //ctor used to discriminate OnLine-Offline analysis   
}

//____________________________________________________________________________ 
AliITSQASPDDataMakerSim::AliITSQASPDDataMakerSim(const AliITSQASPDDataMakerSim& qadm) :
TObject(),
fAliITSQADataMakerSim(qadm.fAliITSQADataMakerSim),
fSPDhHTask(qadm.fSPDhHTask),
fSPDhSTask(qadm.fSPDhSTask),
fSPDhDTask(qadm.fSPDhDTask),
fGenOffsetH(qadm.fGenOffsetH),
fGenOffsetS(qadm.fGenOffsetS),
fGenOffsetD(qadm.fGenOffsetD)
{
  //copy ctor 
  fAliITSQADataMakerSim->SetName((const char*)qadm.fAliITSQADataMakerSim->GetName()) ; 
  fAliITSQADataMakerSim->SetTitle((const char*)qadm.fAliITSQADataMakerSim->GetTitle());
  }

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

//____________________________________________________________________________ 
void AliITSQASPDDataMakerSim::StartOfDetectorCycle()
{
  //Detector specific actions at start of cycle
  AliDebug(1,"AliITSQADM::Start of SPD Cycle\n");
}

//____________________________________________________________________________ 
void AliITSQASPDDataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX_t /*task*/, TObjArray* /*list*/)
{
  // launch the QA checking
  AliDebug(1,"AliITSDM instantiates checker with Run(AliQA::kITS, task, list)\n"); 
  
  //AliQAChecker::Instance()->Run( AliQA::kITS , task, list);
}

//____________________________________________________________________________ 
void AliITSQASPDDataMakerSim::InitDigits()
{ 
  // Initialization for DIGIT data - SPD -
  fGenOffsetD = (fAliITSQADataMakerSim->fDigitsQAList[AliRecoParam::kDefault])->GetEntries();
  //fSPDhDTask must be incremented by one unit every time a histogram is ADDED to the QA List

  Char_t name[50];
  Char_t title[50];

  TH1F *hlayer = new TH1F("SPDLayPattern_SPD","Layer map - SPD",6,0.,6.);
  hlayer->GetXaxis()->SetTitle("Layer number");
  hlayer->GetYaxis()->SetTitle("Entries");
  fAliITSQADataMakerSim->Add2DigitsList(hlayer,fGenOffsetD,kTRUE);
  fSPDhDTask++;
  
  TH1F **hmod = new TH1F*[2];
  for (Int_t iLay=0; iLay<2; iLay++) {
    sprintf(name,"SPDModPattern_SPD%d",iLay+1);
    sprintf(title,"Module map - SPD Layer %d",iLay+1);
    hmod[iLay]=new TH1F(name,title,240,0,240);
    hmod[iLay]->GetXaxis()->SetTitle("Module number");
    hmod[iLay]->GetYaxis()->SetTitle("Entries");
    fAliITSQADataMakerSim->Add2DigitsList(hmod[iLay],1+iLay+fGenOffsetD);
    fSPDhDTask++;
  }
  
  TH1F *hcolumns = new TH1F("SPDColumns_SPD","Columns - SPD",160,0.,160.);
  hcolumns->GetXaxis()->SetTitle("Column number");
  hcolumns->GetYaxis()->SetTitle("Entries");
  fAliITSQADataMakerSim->Add2DigitsList(hcolumns,3+fGenOffsetD,kTRUE);
  fSPDhDTask++;

  TH1F *hrows = new TH1F("SPDRows_SPD","Rows - SPD",256,0.,256.);
  hrows->GetXaxis()->SetTitle("Row number");
  hrows->GetYaxis()->SetTitle("Entries");
  fAliITSQADataMakerSim->Add2DigitsList(hrows,4+fGenOffsetD,kTRUE);
  fSPDhDTask++;

  TH1F** hMultSPDdigits = new TH1F*[2];
  for (Int_t iLay=0; iLay<2; ++iLay) {
    sprintf(name,"SPDDigitMultiplicity_SPD%d",iLay+1);
    sprintf(title,"Digit multiplicity - SPD Layer %d",iLay+1);
    hMultSPDdigits[iLay]=new TH1F(name,title,200,0.,200.);
    hMultSPDdigits[iLay]->GetXaxis()->SetTitle("Digit multiplicity");
    hMultSPDdigits[iLay]->GetYaxis()->SetTitle("Entries");
    fAliITSQADataMakerSim->Add2DigitsList(hMultSPDdigits[iLay], 5+iLay+fGenOffsetD);
    fSPDhDTask++;
  }

  TH2F *hMultSPDdig2MultSPDdig1 
       = new TH2F("SPDDigitMultCorrelation_SPD","Digit multiplicity correlation - SPD",200,0.,200.,200,0.,200.);
  hMultSPDdig2MultSPDdig1->GetXaxis()->SetTitle("Digit multiplicity (Layer 1)");
  hMultSPDdig2MultSPDdig1->GetYaxis()->SetTitle("Digit multiplicity (Layer 2)");
  fAliITSQADataMakerSim->Add2DigitsList(hMultSPDdig2MultSPDdig1,7+fGenOffsetD);
  fSPDhDTask++;

  AliDebug(1,Form("%d SPD Digits histograms booked\n",fSPDhDTask));

}

//____________________________________________________________________________
void AliITSQASPDDataMakerSim::MakeDigits(TTree *digits)
{ 
  // Fill QA for DIGIT - SPD -
  AliITS *fITS  = (AliITS*)gAlice->GetModule("ITS");
  fITS->SetTreeAddress();
  TClonesArray *iITSdigits  = fITS->DigitsAddress(0);  // 0->SPD

  Int_t nDigitsL1=0;
  Int_t nDigitsL2=0;

  for (Int_t imod=0; imod<240; ++imod){
    digits->GetEvent(imod);
    Int_t ndigits = iITSdigits->GetEntries();
    if (imod<80) {
      fAliITSQADataMakerSim->GetDigitsData(0+fGenOffsetD)->Fill(0.5,ndigits);
      fAliITSQADataMakerSim->GetDigitsData(1+fGenOffsetD)->Fill(imod,ndigits);
      nDigitsL1+=ndigits;
    }
    else {
      fAliITSQADataMakerSim->GetDigitsData(0+fGenOffsetD)->Fill(1,ndigits);
      fAliITSQADataMakerSim->GetDigitsData(2+fGenOffsetD)->Fill(imod,ndigits);
      nDigitsL2+=ndigits;
    }
    for (Int_t idig=0; idig<ndigits; ++idig) {
      AliITSdigit *dig=(AliITSdigit*)iITSdigits->UncheckedAt(idig);
      Int_t col=dig->GetCoord1();  // cell number z
      Int_t row=dig->GetCoord2();  // cell number x
      fAliITSQADataMakerSim->GetDigitsData(3+fGenOffsetD)->Fill(col);
      fAliITSQADataMakerSim->GetDigitsData(4+fGenOffsetD)->Fill(row);
    }
  }
  fAliITSQADataMakerSim->GetDigitsData(5+fGenOffsetD)->Fill(nDigitsL1);
  fAliITSQADataMakerSim->GetDigitsData(6+fGenOffsetD)->Fill(nDigitsL2);
  fAliITSQADataMakerSim->GetDigitsData(7+fGenOffsetD)->Fill(nDigitsL1,nDigitsL2);
}

//____________________________________________________________________________ 
void AliITSQASPDDataMakerSim::InitSDigits()
{ 
  // Initialization for SDIGIT data - SPD -
  fGenOffsetS = (fAliITSQADataMakerSim->fSDigitsQAList[AliRecoParam::kDefault])->GetEntries();
  //printf("--W-- AliITSQASPDDataMakerSim::InitSDigits()  fGenOffset= %d \n",fGenOffset);
  //fSPDhSTask must be incremented by one unit every time a histogram is ADDED to the QA List
  
  Char_t name[50];
  Char_t title[50];

  TH1F *hlayer = new TH1F("SPDLayPattern_SPD","Layer map - SPD",6,0.,6.);
  hlayer->GetXaxis()->SetTitle("Layer number");
  hlayer->GetYaxis()->SetTitle("Entries");
  fAliITSQADataMakerSim->Add2SDigitsList(hlayer,fGenOffsetS,kTRUE);
  fSPDhSTask++;

  TH1F **hmod = new TH1F*[2];
  for (Int_t iLay=0; iLay<2; ++iLay) {
    sprintf(name,"SPDModPattern_SPD%d",iLay+1);
    sprintf(title,"Module map - SPD Layer %d",iLay+1);
    hmod[iLay]=new TH1F(name,title,240,0,240);
    hmod[iLay]->GetXaxis()->SetTitle("Module number");
    hmod[iLay]->GetYaxis()->SetTitle("Entries");
    fAliITSQADataMakerSim->Add2SDigitsList(hmod[iLay],1+iLay+fGenOffsetS);
    fSPDhSTask++;
  }

  AliDebug(1,Form("%d SPD SDigits histograms booked\n",fSPDhSTask));

}

//____________________________________________________________________________
void AliITSQASPDDataMakerSim::MakeSDigits(TTree *sdigits)
{ 
  // Fill QA for SDIGIT - SPD -
  TBranch *brchSDigits = sdigits->GetBranch("ITS");
  for (Int_t imod=0; imod<240; ++imod){
    TClonesArray * sdig = new TClonesArray( "AliITSpListItem",1000 );
    brchSDigits->SetAddress( &sdig );
    brchSDigits->GetEvent(imod);
    Int_t nsdig=sdig->GetEntries();
    if (imod<80) {
      fAliITSQADataMakerSim->GetSDigitsData(0+fGenOffsetS)->Fill(0.5,nsdig);
      fAliITSQADataMakerSim->GetSDigitsData(1+fGenOffsetS)->Fill(imod,nsdig);
    }
    else {
      fAliITSQADataMakerSim->GetSDigitsData(0+fGenOffsetS)->Fill(1,nsdig);
      fAliITSQADataMakerSim->GetSDigitsData(2+fGenOffsetS)->Fill(imod,nsdig);
    }
    delete sdig;
  }

}

//____________________________________________________________________________ 
void AliITSQASPDDataMakerSim::InitHits()
{ 
  // Initialization for HITS data - SPD -
  fGenOffsetH = (fAliITSQADataMakerSim->fHitsQAList[AliRecoParam::kDefault])->GetEntries();
  //printf("--W-- AliITSQASPDDataMakerSim::InitHits()  fGenOffset= %d \n",fGenOffset);
  //fSPDhHTask must be incremented by one unit every time a histogram is ADDED to the QA List
  Char_t name[50];
  Char_t title[50];

  TH1F *hlayer = new TH1F("SPDLayPattern_SPD","Layer map - SPD",6,0.,6.);
  hlayer->GetXaxis()->SetTitle("Layer number");
  hlayer->GetYaxis()->SetTitle("Entries");
  fAliITSQADataMakerSim->Add2HitsList(hlayer,fGenOffsetH,kTRUE);
  fSPDhHTask++;

  TH1F **hmod = new TH1F*[2];
  for (Int_t iLay=0; iLay<2; ++iLay) {
    sprintf(name,"SPDModPattern_SPD%d",iLay+1);
    sprintf(title,"Module map - SPD Layer %d",iLay+1);
    hmod[iLay]=new TH1F(name,title,240,0,240);
    hmod[iLay]->GetXaxis()->SetTitle("Module number");
    hmod[iLay]->GetYaxis()->SetTitle("Entries");
    fAliITSQADataMakerSim->Add2HitsList(hmod[iLay],1+iLay+fGenOffsetH);
    fSPDhHTask++;
  }

  TH1F *hhitlenght = new TH1F("SPDLenght_SPD","Hit lenght along y_{loc} coord",210,0.,210.);
  hhitlenght->GetXaxis()->SetTitle("Hit lenght [#mum]");
  hhitlenght->GetYaxis()->SetTitle("# hits");
  fAliITSQADataMakerSim->Add2HitsList(hhitlenght,3+fGenOffsetH);
  fSPDhHTask++;

  TH1F *hEdepos = new TH1F("SPDEnergyDeposit_SPD","Deposited energy distribution (y_{loc}>180 #mum)",150,0.,300.); 
  hEdepos->GetXaxis()->SetTitle("Deposited energy [keV]"); 
  hEdepos->GetYaxis()->SetTitle("# hits");
  fAliITSQADataMakerSim->Add2HitsList(hEdepos,4+fGenOffsetH);
  fSPDhHTask++;

  AliDebug(1,Form("%d SPD Hits histograms booked\n",fSPDhHTask));

}

//____________________________________________________________________________
void AliITSQASPDDataMakerSim::MakeHits(TTree *hits)
{ 
  // Fill QA for HITS - SPD -
  AliITS *fITS  = (AliITS*)gAlice->GetModule("ITS");
  fITS->SetTreeAddress();
  Int_t nmodules;
  fITS->InitModules(-1,nmodules); //-1->number of modules taken from AliITSgeom class kept in fITSgeom
                                  //nmodules is set

  fITS->FillModules(hits,0);

  for (Int_t imod=0; imod<240; ++imod){
    AliITSmodule *module = fITS->GetModule(imod);
    TObjArray *arrHits = module->GetHits();
    Int_t nhits = arrHits->GetEntriesFast();
    if (imod<80) {
      fAliITSQADataMakerSim->GetHitsData(fGenOffsetH)->Fill(0.5,nhits);
      fAliITSQADataMakerSim->GetHitsData(1+fGenOffsetH)->Fill(imod,nhits);
    } else {
      fAliITSQADataMakerSim->GetHitsData(fGenOffsetH)->Fill(1,nhits);
      fAliITSQADataMakerSim->GetHitsData(2+fGenOffsetH)->Fill(imod,nhits);
    }
    for (Int_t iHit=0; iHit<nhits; ++iHit) {
      AliITShit *hit = (AliITShit*) arrHits->At(iHit);
      Double_t xl,yl,zl,xl0,yl0,zl0;
      Double_t tof,tof0;
      hit->GetPositionL(xl,yl,zl,tof);
      hit->GetPositionL0(xl0,yl0,zl0,tof0);
      Float_t dyloc=TMath::Abs(yl-yl0)*10000.;
      fAliITSQADataMakerSim->GetHitsData(3+fGenOffsetH)->Fill(dyloc);
      Float_t edep=hit->GetIonization()*1000000;
      if(dyloc>180.){
        fAliITSQADataMakerSim->GetHitsData(4+fGenOffsetH)->Fill(edep);
      }
    }
  }
}


//_______________________________________________________________

Int_t AliITSQASPDDataMakerSim::GetOffset(AliQA::TASKINDEX_t task){
  // Returns histogram offset according to the specified task
  Int_t offset=0;
  if( task == AliQA::kHITS){
    offset=fGenOffsetH;
  }
  else if( task == AliQA::kSDIGITS) {
    offset=fGenOffsetS;
  }
  else if( task == AliQA::kDIGITS) {
    offset=fGenOffsetD;
  }
  else {
    AliInfo("No task has been selected. TaskHisto set to zero.\n");
  }

  return offset;
}


//_______________________________________________________________

Int_t AliITSQASPDDataMakerSim::GetTaskHisto(AliQA::TASKINDEX_t task) {
  // Returns the number of booked histograms for the selected task
  Int_t histotot=0;
  if( task == AliQA::kHITS) {
    histotot=fSPDhHTask ;
  }
  else if( task == AliQA::kSDIGITS) {
    histotot=fSPDhSTask;
  }
  else if( task == AliQA::kDIGITS) {
    histotot=fSPDhDTask ;
  }
  else {
    AliInfo("No task has been selected. TaskHisto set to zero.\n");
  }
  return histotot;

}
Commit	Line	Data
379510c2	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	// Checks the quality assurance
	19	// by comparing with reference data
	20	// contained in a DB
	21	// -------------------------------------------------------------
	22	// W. Ferrarese + P. Cerello INFN Torino Feb 2008
	23	// M. Nicassio D. Elia INFN Bari April 2008
	24	// maria.nicassio@ba.infn.it
	25
	26	// --- ROOT system ---
	27	#include <TTree.h>
	28	#include <TH2.h>
	29	#include <TH1.h>
	30	// --- Standard library ---
	31
	32	// --- AliRoot header files ---
	33	#include "AliRun.h"
	34	#include "AliITSQADataMakerSim.h"
	35	#include "AliITSQASPDDataMakerSim.h"
	36	#include "AliQA.h"
	37	#include "AliQAChecker.h"
	38	#include "AliITSdigit.h"
	39	#include "AliITSdigitSPD.h"
	40	#include "AliITS.h"
	41	#include "AliITSmodule.h"
	42	#include "AliITShit.h"
	43	#include "AliITSLoader.h"
	44	#include "AliRunLoader.h"
	45
	46	ClassImp(AliITSQASPDDataMakerSim)
	47
	48	//____________________________________________________________________________
	49	AliITSQASPDDataMakerSim::AliITSQASPDDataMakerSim(AliITSQADataMakerSim *aliITSQADataMakerSim) :
	50	TObject(),
	51	fAliITSQADataMakerSim(aliITSQADataMakerSim),
7a0e5776	52	fSPDhHTask(0),
	53	fSPDhSTask(0),
	54	fSPDhDTask(0),
3f905799	55	fGenOffsetH(0),
	56	fGenOffsetS(0),
	57	fGenOffsetD(0)
379510c2	58	{
	59	//ctor used to discriminate OnLine-Offline analysis
	60	}
	61
	62	//____________________________________________________________________________
	63	AliITSQASPDDataMakerSim::AliITSQASPDDataMakerSim(const AliITSQASPDDataMakerSim& qadm) :
	64	TObject(),
	65	fAliITSQADataMakerSim(qadm.fAliITSQADataMakerSim),
7a0e5776	66	fSPDhHTask(qadm.fSPDhHTask),
	67	fSPDhSTask(qadm.fSPDhSTask),
	68	fSPDhDTask(qadm.fSPDhDTask),
3f905799	69	fGenOffsetH(qadm.fGenOffsetH),
	70	fGenOffsetS(qadm.fGenOffsetS),
	71	fGenOffsetD(qadm.fGenOffsetD)
379510c2	72	{
	73	//copy ctor
	74	fAliITSQADataMakerSim->SetName((const char*)qadm.fAliITSQADataMakerSim->GetName()) ;
	75	fAliITSQADataMakerSim->SetTitle((const char*)qadm.fAliITSQADataMakerSim->GetTitle());
	76	}
	77
	78	//__________________________________________________________________
	79	AliITSQASPDDataMakerSim& AliITSQASPDDataMakerSim::operator = (const AliITSQASPDDataMakerSim& qac )
	80	{
	81	// Equal operator.
	82	this->~AliITSQASPDDataMakerSim();
	83	new(this) AliITSQASPDDataMakerSim(qac);
	84	return *this;
	85	}
	86
	87	//____________________________________________________________________________
	88	void AliITSQASPDDataMakerSim::StartOfDetectorCycle()
	89	{
	90	//Detector specific actions at start of cycle
	91	AliDebug(1,"AliITSQADM::Start of SPD Cycle\n");
	92	}
	93
	94	//____________________________________________________________________________
	95	void AliITSQASPDDataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX_t /task/, TObjArray* /list/)
	96	{
	97	// launch the QA checking
	98	AliDebug(1,"AliITSDM instantiates checker with Run(AliQA::kITS, task, list)\n");
	99
	100	//AliQAChecker::Instance()->Run( AliQA::kITS , task, list);
	101	}
	102
	103	//____________________________________________________________________________
	104	void AliITSQASPDDataMakerSim::InitDigits()
	105	{
	106	// Initialization for DIGIT data - SPD -
57acd2d2	107	fGenOffsetD = (fAliITSQADataMakerSim->fDigitsQAList[AliRecoParam::kDefault])->GetEntries();
7a0e5776	108	//fSPDhDTask must be incremented by one unit every time a histogram is ADDED to the QA List
379510c2	109
	110	Char_t name[50];
	111	Char_t title[50];
	112
26ee9565	113	TH1F *hlayer = new TH1F("SPDLayPattern_SPD","Layer map - SPD",6,0.,6.);
379510c2	114	hlayer->GetXaxis()->SetTitle("Layer number");
379510c2	115	hlayer->GetYaxis()->SetTitle("Entries");
954ef57a	116	fAliITSQADataMakerSim->Add2DigitsList(hlayer,fGenOffsetD,kTRUE);
7a0e5776	117	fSPDhDTask++;
379510c2	118
	119	TH1F *hmod = new TH1F[2];
	120	for (Int_t iLay=0; iLay<2; iLay++) {
26ee9565	121	sprintf(name,"SPDModPattern_SPD%d",iLay+1);
379510c2	122	sprintf(title,"Module map - SPD Layer %d",iLay+1);
	123	hmod[iLay]=new TH1F(name,title,240,0,240);
	124	hmod[iLay]->GetXaxis()->SetTitle("Module number");
	125	hmod[iLay]->GetYaxis()->SetTitle("Entries");
3f905799	126	fAliITSQADataMakerSim->Add2DigitsList(hmod[iLay],1+iLay+fGenOffsetD);
7a0e5776	127	fSPDhDTask++;
379510c2	128	}
379510c2	129
26ee9565	130	TH1F *hcolumns = new TH1F("SPDColumns_SPD","Columns - SPD",160,0.,160.);
379510c2	131	hcolumns->GetXaxis()->SetTitle("Column number");
379510c2	132	hcolumns->GetYaxis()->SetTitle("Entries");
954ef57a	133	fAliITSQADataMakerSim->Add2DigitsList(hcolumns,3+fGenOffsetD,kTRUE);
7a0e5776	134	fSPDhDTask++;
379510c2	135
26ee9565	136	TH1F *hrows = new TH1F("SPDRows_SPD","Rows - SPD",256,0.,256.);
379510c2	137	hrows->GetXaxis()->SetTitle("Row number");
379510c2	138	hrows->GetYaxis()->SetTitle("Entries");
954ef57a	139	fAliITSQADataMakerSim->Add2DigitsList(hrows,4+fGenOffsetD,kTRUE);
7a0e5776	140	fSPDhDTask++;
379510c2	141
	142	TH1F** hMultSPDdigits = new TH1F*[2];
	143	for (Int_t iLay=0; iLay<2; ++iLay) {
26ee9565	144	sprintf(name,"SPDDigitMultiplicity_SPD%d",iLay+1);
379510c2	145	sprintf(title,"Digit multiplicity - SPD Layer %d",iLay+1);
	146	hMultSPDdigits[iLay]=new TH1F(name,title,200,0.,200.);
	147	hMultSPDdigits[iLay]->GetXaxis()->SetTitle("Digit multiplicity");
	148	hMultSPDdigits[iLay]->GetYaxis()->SetTitle("Entries");
3f905799	149	fAliITSQADataMakerSim->Add2DigitsList(hMultSPDdigits[iLay], 5+iLay+fGenOffsetD);
7a0e5776	150	fSPDhDTask++;
379510c2	151	}
379510c2	152
954ef57a	153	TH2F *hMultSPDdig2MultSPDdig1
26ee9565	154	= new TH2F("SPDDigitMultCorrelation_SPD","Digit multiplicity correlation - SPD",200,0.,200.,200,0.,200.);
379510c2	155	hMultSPDdig2MultSPDdig1->GetXaxis()->SetTitle("Digit multiplicity (Layer 1)");
379510c2	156	hMultSPDdig2MultSPDdig1->GetYaxis()->SetTitle("Digit multiplicity (Layer 2)");
3f905799	157	fAliITSQADataMakerSim->Add2DigitsList(hMultSPDdig2MultSPDdig1,7+fGenOffsetD);
7a0e5776	158	fSPDhDTask++;
379510c2	159
7a0e5776	160	AliDebug(1,Form("%d SPD Digits histograms booked\n",fSPDhDTask));
379510c2	161
	162	}
	163
	164	//____________________________________________________________________________
	165	void AliITSQASPDDataMakerSim::MakeDigits(TTree *digits)
	166	{
	167	// Fill QA for DIGIT - SPD -
	168	AliITS fITS = (AliITS)gAlice->GetModule("ITS");
	169	fITS->SetTreeAddress();
	170	TClonesArray *iITSdigits = fITS->DigitsAddress(0); // 0->SPD
	171
	172	Int_t nDigitsL1=0;
	173	Int_t nDigitsL2=0;
	174
	175	for (Int_t imod=0; imod<240; ++imod){
	176	digits->GetEvent(imod);
	177	Int_t ndigits = iITSdigits->GetEntries();
	178	if (imod<80) {
3f905799	179	fAliITSQADataMakerSim->GetDigitsData(0+fGenOffsetD)->Fill(0.5,ndigits);
3f905799	180	fAliITSQADataMakerSim->GetDigitsData(1+fGenOffsetD)->Fill(imod,ndigits);
379510c2	181	nDigitsL1+=ndigits;
	182	}
	183	else {
3f905799	184	fAliITSQADataMakerSim->GetDigitsData(0+fGenOffsetD)->Fill(1,ndigits);
3f905799	185	fAliITSQADataMakerSim->GetDigitsData(2+fGenOffsetD)->Fill(imod,ndigits);
379510c2	186	nDigitsL2+=ndigits;
	187	}
	188	for (Int_t idig=0; idig<ndigits; ++idig) {
	189	AliITSdigit dig=(AliITSdigit)iITSdigits->UncheckedAt(idig);
	190	Int_t col=dig->GetCoord1(); // cell number z
	191	Int_t row=dig->GetCoord2(); // cell number x
3f905799	192	fAliITSQADataMakerSim->GetDigitsData(3+fGenOffsetD)->Fill(col);
3f905799	193	fAliITSQADataMakerSim->GetDigitsData(4+fGenOffsetD)->Fill(row);
379510c2	194	}
379510c2	195	}
3f905799	196	fAliITSQADataMakerSim->GetDigitsData(5+fGenOffsetD)->Fill(nDigitsL1);
	197	fAliITSQADataMakerSim->GetDigitsData(6+fGenOffsetD)->Fill(nDigitsL2);
	198	fAliITSQADataMakerSim->GetDigitsData(7+fGenOffsetD)->Fill(nDigitsL1,nDigitsL2);
379510c2	199	}
	200
	201	//____________________________________________________________________________
	202	void AliITSQASPDDataMakerSim::InitSDigits()
	203	{
	204	// Initialization for SDIGIT data - SPD -
57acd2d2	205	fGenOffsetS = (fAliITSQADataMakerSim->fSDigitsQAList[AliRecoParam::kDefault])->GetEntries();
c71529b0	206	//printf("--W-- AliITSQASPDDataMakerSim::InitSDigits() fGenOffset= %d \n",fGenOffset);
7a0e5776	207	//fSPDhSTask must be incremented by one unit every time a histogram is ADDED to the QA List
c71529b0	208
379510c2	209	Char_t name[50];
	210	Char_t title[50];
	211
26ee9565	212	TH1F *hlayer = new TH1F("SPDLayPattern_SPD","Layer map - SPD",6,0.,6.);
379510c2	213	hlayer->GetXaxis()->SetTitle("Layer number");
379510c2	214	hlayer->GetYaxis()->SetTitle("Entries");
954ef57a	215	fAliITSQADataMakerSim->Add2SDigitsList(hlayer,fGenOffsetS,kTRUE);
7a0e5776	216	fSPDhSTask++;
379510c2	217
	218	TH1F *hmod = new TH1F[2];
	219	for (Int_t iLay=0; iLay<2; ++iLay) {
26ee9565	220	sprintf(name,"SPDModPattern_SPD%d",iLay+1);
379510c2	221	sprintf(title,"Module map - SPD Layer %d",iLay+1);
	222	hmod[iLay]=new TH1F(name,title,240,0,240);
	223	hmod[iLay]->GetXaxis()->SetTitle("Module number");
	224	hmod[iLay]->GetYaxis()->SetTitle("Entries");
3f905799	225	fAliITSQADataMakerSim->Add2SDigitsList(hmod[iLay],1+iLay+fGenOffsetS);
7a0e5776	226	fSPDhSTask++;
379510c2	227	}
379510c2	228
7a0e5776	229	AliDebug(1,Form("%d SPD SDigits histograms booked\n",fSPDhSTask));
379510c2	230
	231	}
	232
	233	//____________________________________________________________________________
	234	void AliITSQASPDDataMakerSim::MakeSDigits(TTree *sdigits)
	235	{
	236	// Fill QA for SDIGIT - SPD -
	237	TBranch *brchSDigits = sdigits->GetBranch("ITS");
	238	for (Int_t imod=0; imod<240; ++imod){
	239	TClonesArray * sdig = new TClonesArray( "AliITSpListItem",1000 );
	240	brchSDigits->SetAddress( &sdig );
	241	brchSDigits->GetEvent(imod);
	242	Int_t nsdig=sdig->GetEntries();
	243	if (imod<80) {
3f905799	244	fAliITSQADataMakerSim->GetSDigitsData(0+fGenOffsetS)->Fill(0.5,nsdig);
3f905799	245	fAliITSQADataMakerSim->GetSDigitsData(1+fGenOffsetS)->Fill(imod,nsdig);
379510c2	246	}
379510c2	247	else {
3f905799	248	fAliITSQADataMakerSim->GetSDigitsData(0+fGenOffsetS)->Fill(1,nsdig);
3f905799	249	fAliITSQADataMakerSim->GetSDigitsData(2+fGenOffsetS)->Fill(imod,nsdig);
379510c2	250	}
	251	delete sdig;
	252	}
	253
	254	}
	255
	256	//____________________________________________________________________________
	257	void AliITSQASPDDataMakerSim::InitHits()
	258	{
	259	// Initialization for HITS data - SPD -
57acd2d2	260	fGenOffsetH = (fAliITSQADataMakerSim->fHitsQAList[AliRecoParam::kDefault])->GetEntries();
c71529b0	261	//printf("--W-- AliITSQASPDDataMakerSim::InitHits() fGenOffset= %d \n",fGenOffset);
7a0e5776	262	//fSPDhHTask must be incremented by one unit every time a histogram is ADDED to the QA List
379510c2	263	Char_t name[50];
	264	Char_t title[50];
	265
26ee9565	266	TH1F *hlayer = new TH1F("SPDLayPattern_SPD","Layer map - SPD",6,0.,6.);
379510c2	267	hlayer->GetXaxis()->SetTitle("Layer number");
379510c2	268	hlayer->GetYaxis()->SetTitle("Entries");
954ef57a	269	fAliITSQADataMakerSim->Add2HitsList(hlayer,fGenOffsetH,kTRUE);
7a0e5776	270	fSPDhHTask++;
379510c2	271
	272	TH1F *hmod = new TH1F[2];
	273	for (Int_t iLay=0; iLay<2; ++iLay) {
26ee9565	274	sprintf(name,"SPDModPattern_SPD%d",iLay+1);
379510c2	275	sprintf(title,"Module map - SPD Layer %d",iLay+1);
	276	hmod[iLay]=new TH1F(name,title,240,0,240);
	277	hmod[iLay]->GetXaxis()->SetTitle("Module number");
	278	hmod[iLay]->GetYaxis()->SetTitle("Entries");
3f905799	279	fAliITSQADataMakerSim->Add2HitsList(hmod[iLay],1+iLay+fGenOffsetH);
7a0e5776	280	fSPDhHTask++;
379510c2	281	}
379510c2	282
26ee9565	283	TH1F *hhitlenght = new TH1F("SPDLenght_SPD","Hit lenght along y_{loc} coord",210,0.,210.);
379510c2	284	hhitlenght->GetXaxis()->SetTitle("Hit lenght [#mum]");
379510c2	285	hhitlenght->GetYaxis()->SetTitle("# hits");
3f905799	286	fAliITSQADataMakerSim->Add2HitsList(hhitlenght,3+fGenOffsetH);
7a0e5776	287	fSPDhHTask++;
379510c2	288
26ee9565	289	TH1F *hEdepos = new TH1F("SPDEnergyDeposit_SPD","Deposited energy distribution (y_{loc}>180 #mum)",150,0.,300.);
379510c2	290	hEdepos->GetXaxis()->SetTitle("Deposited energy [keV]");
379510c2	291	hEdepos->GetYaxis()->SetTitle("# hits");
3f905799	292	fAliITSQADataMakerSim->Add2HitsList(hEdepos,4+fGenOffsetH);
7a0e5776	293	fSPDhHTask++;
379510c2	294
7a0e5776	295	AliDebug(1,Form("%d SPD Hits histograms booked\n",fSPDhHTask));
379510c2	296
	297	}
	298
	299	//____________________________________________________________________________
	300	void AliITSQASPDDataMakerSim::MakeHits(TTree *hits)
	301	{
	302	// Fill QA for HITS - SPD -
	303	AliITS fITS = (AliITS)gAlice->GetModule("ITS");
	304	fITS->SetTreeAddress();
	305	Int_t nmodules;
	306	fITS->InitModules(-1,nmodules); //-1->number of modules taken from AliITSgeom class kept in fITSgeom
	307	//nmodules is set
	308
	309	fITS->FillModules(hits,0);
	310
	311	for (Int_t imod=0; imod<240; ++imod){
	312	AliITSmodule *module = fITS->GetModule(imod);
	313	TObjArray *arrHits = module->GetHits();
	314	Int_t nhits = arrHits->GetEntriesFast();
	315	if (imod<80) {
3f905799	316	fAliITSQADataMakerSim->GetHitsData(fGenOffsetH)->Fill(0.5,nhits);
3f905799	317	fAliITSQADataMakerSim->GetHitsData(1+fGenOffsetH)->Fill(imod,nhits);
379510c2	318	} else {
3f905799	319	fAliITSQADataMakerSim->GetHitsData(fGenOffsetH)->Fill(1,nhits);
3f905799	320	fAliITSQADataMakerSim->GetHitsData(2+fGenOffsetH)->Fill(imod,nhits);
379510c2	321	}
379510c2	322	for (Int_t iHit=0; iHit<nhits; ++iHit) {
	323	AliITShit hit = (AliITShit) arrHits->At(iHit);
	324	Double_t xl,yl,zl,xl0,yl0,zl0;
	325	Double_t tof,tof0;
	326	hit->GetPositionL(xl,yl,zl,tof);
	327	hit->GetPositionL0(xl0,yl0,zl0,tof0);
	328	Float_t dyloc=TMath::Abs(yl-yl0)*10000.;
3f905799	329	fAliITSQADataMakerSim->GetHitsData(3+fGenOffsetH)->Fill(dyloc);
379510c2	330	Float_t edep=hit->GetIonization()*1000000;
379510c2	331	if(dyloc>180.){
3f905799	332	fAliITSQADataMakerSim->GetHitsData(4+fGenOffsetH)->Fill(edep);
379510c2	333	}
	334	}
	335	}
	336	}
7a0e5776	337
	338
	339	//_______________________________________________________________
	340
	341	Int_t AliITSQASPDDataMakerSim::GetOffset(AliQA::TASKINDEX_t task){
	342	// Returns histogram offset according to the specified task
	343	Int_t offset=0;
	344	if( task == AliQA::kHITS){
	345	offset=fGenOffsetH;
	346	}
	347	else if( task == AliQA::kSDIGITS) {
	348	offset=fGenOffsetS;
	349	}
	350	else if( task == AliQA::kDIGITS) {
	351	offset=fGenOffsetD;
	352	}
	353	else {
	354	AliInfo("No task has been selected. TaskHisto set to zero.\n");
	355	}
	356
	357	return offset;
	358	}
	359
	360
	361	//_______________________________________________________________
	362
	363	Int_t AliITSQASPDDataMakerSim::GetTaskHisto(AliQA::TASKINDEX_t task) {
	364	// Returns the number of booked histograms for the selected task
	365	Int_t histotot=0;
	366	if( task == AliQA::kHITS) {
	367	histotot=fSPDhHTask ;
	368	}
	369	else if( task == AliQA::kSDIGITS) {
	370	histotot=fSPDhSTask;
	371	}
	372	else if( task == AliQA::kDIGITS) {
	373	histotot=fSPDhDTask ;
	374	}
	375	else {
	376	AliInfo("No task has been selected. TaskHisto set to zero.\n");
	377	}
	378	return histotot;
	379
	380	}