[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")
{
  // ctor
}

//____________________________________________________________________________ 
AliPHOSQADataMakerSim::AliPHOSQADataMakerSim(const AliPHOSQADataMakerSim& qadm) :
  AliQADataMakerSim()
{
  //copy ctor 
  SetName((const char*)qadm.GetName()) ; 
  SetTitle((const char*)qadm.GetTitle()); 
}

//__________________________________________________________________
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(TClonesArray * hits)
{
	//make QA data from Hits

    GetHitsData(1)->Fill(hits->GetEntriesFast()) ; 
    TIter next(hits) ; 
    AliPHOSHit * hit ; 
    while ( (hit = dynamic_cast<AliPHOSHit *>(next())) ) {
      GetHitsData(0)->Fill( hit->GetEnergy()) ;
    }
}

//____________________________________________________________________________
void AliPHOSQADataMakerSim::MakeHits(TTree * hitTree)
{
	// make QA data from Hit Tree
	
	TClonesArray * hits = new TClonesArray("AliPHOSHit", 1000);

	TBranch * branch = hitTree->GetBranch("PHOS") ;
	if ( ! branch ) {
		AliWarning("PHOS branch in Hit Tree not found") ; 
	} else {
		TClonesArray * tmp =  new TClonesArray("AliPHOSHit", 1000) ;
		branch->SetAddress(&tmp) ;
		Int_t index = 0 ;  
		for (Int_t ientry = 0 ; ientry < branch->GetEntries() ; ientry++) {
			branch->GetEntry(ientry) ; 
			for (Int_t ihit = 0 ; ihit < tmp->GetEntries() ; ihit++) {
				AliPHOSHit * hit = dynamic_cast<AliPHOSHit *> (tmp->At(ihit)) ; 
				new((*hits)[index]) AliPHOSHit(*hit) ; 
				index++ ;
			} 
		} 	
		tmp->Delete() ; 
		delete tmp ; 
		MakeHits(hits) ; 
	}
}

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