[u/mrichter/AliRoot.git] / PHOS / AliPHOSQADataMakerSim.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$ */

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

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

// --- Standard library ---

// --- AliRoot header files ---
#include "AliESDCaloCluster.h"
#include "AliLog.h"
#include "AliPHOSDigit.h"
#include "AliPHOSHit.h"
#include "AliPHOSQADataMakerSim.h"
#include "AliQAChecker.h"

ClassImp(AliPHOSQADataMakerSim)
           
//____________________________________________________________________________ 
AliPHOSQADataMakerSim::AliPHOSQADataMakerSim() : 
  AliQADataMakerSim(AliQA::GetDetName(AliQA::kPHOS), "PHOS Quality Assurance Data Maker"),
  fHits(0x0)
{
  // ctor
  fHits = new TClonesArray("AliPHOSHit", 1000);
}

//____________________________________________________________________________ 
AliPHOSQADataMakerSim::AliPHOSQADataMakerSim(const AliPHOSQADataMakerSim& qadm) :
  AliQADataMakerSim(),
  fHits(0x0)
{
  //copy ctor 
  SetName((const char*)qadm.GetName()) ; 
  SetTitle((const char*)qadm.GetTitle()); 
  fHits = new TClonesArray("AliPHOSHit", 1000);
}

//__________________________________________________________________
AliPHOSQADataMakerSim& AliPHOSQADataMakerSim::operator = (const AliPHOSQADataMakerSim& qadm )
{
  // Assign operator.
  this->~AliPHOSQADataMakerSim();
  new(this) AliPHOSQADataMakerSim(qadm);
  return *this;
}
 
//____________________________________________________________________________ 
void AliPHOSQADataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX_t task, TObjArray * list)
{
  //Detector specific actions at end of cycle
  // do the QA checking
  AliQAChecker::Instance()->Run(AliQA::kPHOS, task, list) ;  
}

//____________________________________________________________________________ 
void AliPHOSQADataMakerSim::InitHits()
{
  // create Hits histograms in Hits subdir
  TH1F * h0 = new TH1F("hPhosHits",    "Hits energy distribution in PHOS",       100, 0., 100.) ; 
  h0->Sumw2() ;
  Add2HitsList(h0, 0) ;
  TH1I * h1 = new TH1I("hPhosHitsMul", "Hits multiplicity distribution in PHOS", 500, 0., 10000) ; 
  h1->Sumw2() ;
  Add2HitsList(h1, 1) ;
}

//____________________________________________________________________________ 
void AliPHOSQADataMakerSim::InitDigits()
{
  // create Digits histograms in Digits subdir
  TH1I * h0 = new TH1I("hPhosDigits",    "Digits amplitude distribution in PHOS",    500, 0, 5000) ; 
  h0->Sumw2() ;
  Add2DigitsList(h0, 0) ;
  TH1I * h1 = new TH1I("hPhosDigitsMul", "Digits multiplicity distribution in PHOS", 500, 0, 1000) ; 
  h1->Sumw2() ;
  Add2DigitsList(h1, 1) ;
}

//____________________________________________________________________________ 
void AliPHOSQADataMakerSim::InitSDigits()
{
  // create SDigits histograms in SDigits subdir
  TH1F * h0 = new TH1F("hPhosSDigits",    "SDigits energy distribution in PHOS",       100, 0., 100.) ; 
  h0->Sumw2() ;
  Add2SDigitsList(h0, 0) ;
  TH1I * h1 = new TH1I("hPhosSDigitsMul", "SDigits multiplicity distribution in PHOS", 500, 0,  10000) ; 
  h1->Sumw2() ;
  Add2SDigitsList(h1, 1) ;
}

//____________________________________________________________________________
void AliPHOSQADataMakerSim::MakeHits()
{
  //make QA data from Hits
  
  TIter next(fHits) ; 
  AliPHOSHit * hit ; 
  while ( (hit = dynamic_cast<AliPHOSHit *>(next())) ) {
    GetHitsData(0)->Fill( hit->GetEnergy()) ;
  }
}

//____________________________________________________________________________
void AliPHOSQADataMakerSim::MakeHits(TTree * hitTree)
{
  // make QA data from Hit Tree
  
  TBranch * branch = hitTree->GetBranch("PHOS") ;
  if ( ! branch ) {
    AliWarning("PHOS branch in Hit Tree not found") ; 
  } else {
    Int_t nHits = 0;
    branch->SetAddress(&fHits) ;
    for (Int_t ientry = 0 ; ientry < branch->GetEntries() ; ientry++) {
      branch->GetEntry(ientry) ;
      nHits += fHits->GetEntriesFast();
      MakeHits() ; 
      fHits->Clear();
    } 	
    GetHitsData(1)->Fill(nHits) ;
  }
}

//____________________________________________________________________________
void AliPHOSQADataMakerSim::MakeDigits(TClonesArray * digits)
{
  // makes data from Digits

    GetDigitsData(1)->Fill(digits->GetEntriesFast()) ; 
    TIter next(digits) ; 
    AliPHOSDigit * digit ; 
    while ( (digit = dynamic_cast<AliPHOSDigit *>(next())) ) {
      GetDigitsData(0)->Fill( digit->GetEnergy()) ;
    }  
}

//____________________________________________________________________________
void AliPHOSQADataMakerSim::MakeDigits(TTree * digitTree)
{
	// makes data from Digit Tree
	TClonesArray * digits = new TClonesArray("AliPHOSDigit", 1000) ; 

	TBranch * branch = digitTree->GetBranch("PHOS") ;
	if ( ! branch ) {
		AliWarning("PHOS branch in Digit Tree not found") ; 
	} else {
		branch->SetAddress(&digits) ;
		branch->GetEntry(0) ; 
		MakeDigits(digits) ; 
	}
}

//____________________________________________________________________________
void AliPHOSQADataMakerSim::MakeSDigits(TClonesArray * sdigits)
{
  // makes data from SDigits
  
	GetSDigitsData(1)->Fill(sdigits->GetEntriesFast()) ; 
    TIter next(sdigits) ; 
    AliPHOSDigit * sdigit ; 
    while ( (sdigit = dynamic_cast<AliPHOSDigit *>(next())) ) {
      GetSDigitsData(0)->Fill( sdigit->GetEnergy()) ;
    } 
}

//____________________________________________________________________________
void AliPHOSQADataMakerSim::MakeSDigits(TTree * sdigitTree)
{
	// makes data from SDigit Tree
	TClonesArray * sdigits = new TClonesArray("AliPHOSDigit", 1000) ; 

	TBranch * branch = sdigitTree->GetBranch("PHOS") ;
	if ( ! branch ) {
		AliWarning("PHOS branch in SDigit Tree not found") ; 
	} else {
		branch->SetAddress(&sdigits) ;
		branch->GetEntry(0) ;
		MakeSDigits(sdigits) ; 
	}
}

//____________________________________________________________________________ 
void AliPHOSQADataMakerSim::StartOfDetectorCycle()
{
  //Detector specific actions at start of cycle
  
}
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
	19	/*
	20	Produces the data needed to calculate the quality assurance.
	21	All data must be mergeable objects.
	22	Y. Schutz CERN July 2007
	23	*/
	24
	25	// --- ROOT system ---
	26	#include <TClonesArray.h>
	27	#include <TFile.h>
	28	#include <TH1F.h>
	29	#include <TH1I.h>
	30	#include <TH2F.h>
	31	#include <TTree.h>
	32
	33	// --- Standard library ---
	34
	35	// --- AliRoot header files ---
	36	#include "AliESDCaloCluster.h"
	37	#include "AliLog.h"
	38	#include "AliPHOSDigit.h"
	39	#include "AliPHOSHit.h"
	40	#include "AliPHOSQADataMakerSim.h"
	41	#include "AliQAChecker.h"
	42
	43	ClassImp(AliPHOSQADataMakerSim)
	44
	45	//____________________________________________________________________________
8d8258f6	46	AliPHOSQADataMakerSim::AliPHOSQADataMakerSim() :
	47	AliQADataMakerSim(AliQA::GetDetName(AliQA::kPHOS), "PHOS Quality Assurance Data Maker"),
	48	fHits(0x0)
04236e67	49	{
04236e67	50	// ctor
8d8258f6	51	fHits = new TClonesArray("AliPHOSHit", 1000);
04236e67	52	}
	53
	54	//____________________________________________________________________________
	55	AliPHOSQADataMakerSim::AliPHOSQADataMakerSim(const AliPHOSQADataMakerSim& qadm) :
8d8258f6	56	AliQADataMakerSim(),
8d8258f6	57	fHits(0x0)
04236e67	58	{
	59	//copy ctor
	60	SetName((const char*)qadm.GetName()) ;
	61	SetTitle((const char*)qadm.GetTitle());
8d8258f6	62	fHits = new TClonesArray("AliPHOSHit", 1000);
04236e67	63	}
	64
	65	//__________________________________________________________________
	66	AliPHOSQADataMakerSim& AliPHOSQADataMakerSim::operator = (const AliPHOSQADataMakerSim& qadm )
	67	{
	68	// Assign operator.
	69	this->~AliPHOSQADataMakerSim();
	70	new(this) AliPHOSQADataMakerSim(qadm);
	71	return *this;
	72	}
	73
	74	//____________________________________________________________________________
92a357bf	75	void AliPHOSQADataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX_t task, TObjArray * list)
04236e67	76	{
	77	//Detector specific actions at end of cycle
	78	// do the QA checking
	79	AliQAChecker::Instance()->Run(AliQA::kPHOS, task, list) ;
	80	}
	81
	82	//____________________________________________________________________________
	83	void AliPHOSQADataMakerSim::InitHits()
	84	{
	85	// create Hits histograms in Hits subdir
	86	TH1F * h0 = new TH1F("hPhosHits", "Hits energy distribution in PHOS", 100, 0., 100.) ;
	87	h0->Sumw2() ;
	88	Add2HitsList(h0, 0) ;
8d8258f6	89	TH1I * h1 = new TH1I("hPhosHitsMul", "Hits multiplicity distribution in PHOS", 500, 0., 10000) ;
04236e67	90	h1->Sumw2() ;
	91	Add2HitsList(h1, 1) ;
	92	}
	93
	94	//____________________________________________________________________________
	95	void AliPHOSQADataMakerSim::InitDigits()
	96	{
	97	// create Digits histograms in Digits subdir
	98	TH1I * h0 = new TH1I("hPhosDigits", "Digits amplitude distribution in PHOS", 500, 0, 5000) ;
	99	h0->Sumw2() ;
	100	Add2DigitsList(h0, 0) ;
	101	TH1I * h1 = new TH1I("hPhosDigitsMul", "Digits multiplicity distribution in PHOS", 500, 0, 1000) ;
	102	h1->Sumw2() ;
	103	Add2DigitsList(h1, 1) ;
	104	}
	105
	106	//____________________________________________________________________________
	107	void AliPHOSQADataMakerSim::InitSDigits()
	108	{
	109	// create SDigits histograms in SDigits subdir
	110	TH1F * h0 = new TH1F("hPhosSDigits", "SDigits energy distribution in PHOS", 100, 0., 100.) ;
	111	h0->Sumw2() ;
	112	Add2SDigitsList(h0, 0) ;
	113	TH1I * h1 = new TH1I("hPhosSDigitsMul", "SDigits multiplicity distribution in PHOS", 500, 0, 10000) ;
	114	h1->Sumw2() ;
	115	Add2SDigitsList(h1, 1) ;
	116	}
	117
	118	//____________________________________________________________________________
8d8258f6	119	void AliPHOSQADataMakerSim::MakeHits()
04236e67	120	{
8d8258f6	121	//make QA data from Hits
	122
	123	TIter next(fHits) ;
	124	AliPHOSHit * hit ;
	125	while ( (hit = dynamic_cast<AliPHOSHit *>(next())) ) {
	126	GetHitsData(0)->Fill( hit->GetEnergy()) ;
	127	}
04236e67	128	}
	129
	130	//____________________________________________________________________________
	131	void AliPHOSQADataMakerSim::MakeHits(TTree * hitTree)
	132	{
8d8258f6	133	// make QA data from Hit Tree
	134
	135	TBranch * branch = hitTree->GetBranch("PHOS") ;
	136	if ( ! branch ) {
	137	AliWarning("PHOS branch in Hit Tree not found") ;
	138	} else {
	139	Int_t nHits = 0;
	140	branch->SetAddress(&fHits) ;
	141	for (Int_t ientry = 0 ; ientry < branch->GetEntries() ; ientry++) {
	142	branch->GetEntry(ientry) ;
	143	nHits += fHits->GetEntriesFast();
	144	MakeHits() ;
	145	fHits->Clear();
	146	}
	147	GetHitsData(1)->Fill(nHits) ;
	148	}
04236e67	149	}
	150
	151	//____________________________________________________________________________
	152	void AliPHOSQADataMakerSim::MakeDigits(TClonesArray * digits)
	153	{
	154	// makes data from Digits
	155
	156	GetDigitsData(1)->Fill(digits->GetEntriesFast()) ;
	157	TIter next(digits) ;
	158	AliPHOSDigit * digit ;
	159	while ( (digit = dynamic_cast<AliPHOSDigit *>(next())) ) {
	160	GetDigitsData(0)->Fill( digit->GetEnergy()) ;
	161	}
	162	}
	163
	164	//____________________________________________________________________________
	165	void AliPHOSQADataMakerSim::MakeDigits(TTree * digitTree)
	166	{
	167	// makes data from Digit Tree
	168	TClonesArray * digits = new TClonesArray("AliPHOSDigit", 1000) ;
	169
	170	TBranch * branch = digitTree->GetBranch("PHOS") ;
	171	if ( ! branch ) {
	172	AliWarning("PHOS branch in Digit Tree not found") ;
	173	} else {
	174	branch->SetAddress(&digits) ;
	175	branch->GetEntry(0) ;
	176	MakeDigits(digits) ;
	177	}
	178	}
	179
	180	//____________________________________________________________________________
	181	void AliPHOSQADataMakerSim::MakeSDigits(TClonesArray * sdigits)
	182	{
	183	// makes data from SDigits
	184
	185	GetSDigitsData(1)->Fill(sdigits->GetEntriesFast()) ;
	186	TIter next(sdigits) ;
	187	AliPHOSDigit * sdigit ;
	188	while ( (sdigit = dynamic_cast<AliPHOSDigit *>(next())) ) {
	189	GetSDigitsData(0)->Fill( sdigit->GetEnergy()) ;
	190	}
	191	}
	192
	193	//____________________________________________________________________________
	194	void AliPHOSQADataMakerSim::MakeSDigits(TTree * sdigitTree)
	195	{
	196	// makes data from SDigit Tree
	197	TClonesArray * sdigits = new TClonesArray("AliPHOSDigit", 1000) ;
	198
	199	TBranch * branch = sdigitTree->GetBranch("PHOS") ;
	200	if ( ! branch ) {
	201	AliWarning("PHOS branch in SDigit Tree not found") ;
	202	} else {
	203	branch->SetAddress(&sdigits) ;
	204	branch->GetEntry(0) ;
	205	MakeSDigits(sdigits) ;
	206	}
	207	}
	208
	209	//____________________________________________________________________________
	210	void AliPHOSQADataMakerSim::StartOfDetectorCycle()
	211	{
	212	//Detector specific actions at start of cycle
213
214	}