[u/mrichter/AliRoot.git] / EMCAL / AliEMCALQADataMakerSim.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.                  *
 **************************************************************************/

/*
  Produces the data needed to calculate the quality assurance. 
  All data must be mergeable objects.

  Based on PHOS code written by
  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 "AliEMCALDigit.h"
#include "AliEMCALHit.h"
#include "AliEMCALQADataMakerSim.h"
#include "AliQAChecker.h"

ClassImp(AliEMCALQADataMakerSim)
           
//____________________________________________________________________________ 
  AliEMCALQADataMakerSim::AliEMCALQADataMakerSim() : 
  AliQADataMakerSim(AliQA::GetDetName(AliQA::kEMCAL), "EMCAL Quality Assurance Data Maker")
{
  // ctor
}

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

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

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

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

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

//____________________________________________________________________________
void AliEMCALQADataMakerSim::MakeHits(TClonesArray * hits)
{
  //make QA data from Hits

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

}

//____________________________________________________________________________
void AliEMCALQADataMakerSim::MakeHits(TTree * hitTree)
{
  // make QA data from Hit Tree
  
  TClonesArray * hits = new TClonesArray("AliEMCALHit", 1000);
  
  TBranch * branch = hitTree->GetBranch("EMCAL") ;
  if ( ! branch ) {
    AliWarning("EMCAL branch in Hit Tree not found") ; 
  } else {
    TClonesArray * tmp =  new TClonesArray("AliEMCALHit", 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++) {
	AliEMCALHit * hit = dynamic_cast<AliEMCALHit *> (tmp->At(ihit)) ; 
	new((*hits)[index]) AliEMCALHit(*hit) ; 
	index++ ;
      } 
    } 	
    tmp->Delete() ; 
    delete tmp ; 
    MakeHits(hits) ; 
  }
}

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

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

}

//____________________________________________________________________________
void AliEMCALQADataMakerSim::MakeDigits(TTree * digitTree)
{
  // makes data from Digit Tree
  TClonesArray * digits = new TClonesArray("AliEMCALDigit", 1000) ; 
  
  TBranch * branch = digitTree->GetBranch("EMCAL") ;
  if ( ! branch ) {
    AliWarning("EMCAL branch in Digit Tree not found") ; 
  } else {
    branch->SetAddress(&digits) ;
    branch->GetEntry(0) ; 
    MakeDigits(digits) ; 
  }

}

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

}

//____________________________________________________________________________
void AliEMCALQADataMakerSim::MakeSDigits(TTree * sdigitTree)
{
  // makes data from SDigit Tree
  TClonesArray * sdigits = new TClonesArray("AliEMCALDigit", 1000) ; 
  
  TBranch * branch = sdigitTree->GetBranch("EMCAL") ;
  if ( ! branch ) {
    AliWarning("EMCAL branch in SDigit Tree not found") ; 
  } else {
    branch->SetAddress(&sdigits) ;
    branch->GetEntry(0) ;
    MakeSDigits(sdigits) ; 
  }

}

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