[u/mrichter/AliRoot.git] / PMD / AliPMDQADataMakerSim.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.
  B.K. Nandi
*/

// --- 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 "AliLog.h"
#include "AliPMDhit.h"
#include "AliPMDsdigit.h"
#include "AliPMDdigit.h"
#include "AliPMDQADataMakerSim.h"
#include "AliQAChecker.h"

ClassImp(AliPMDQADataMakerSim)
           
//____________________________________________________________________________ 
AliPMDQADataMakerSim::AliPMDQADataMakerSim() : 
    AliQADataMakerSim(AliQA::GetDetName(AliQA::kPMD), "PMD Quality Assurance Data Maker")
{
    // ctor
}

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

//__________________________________________________________________
AliPMDQADataMakerSim& AliPMDQADataMakerSim::operator = (const AliPMDQADataMakerSim& qadm )
{
    // Assign operator.
    this->~AliPMDQADataMakerSim();
    new(this) AliPMDQADataMakerSim(qadm);
    return *this;
}
 
//____________________________________________________________________________ 
void AliPMDQADataMakerSim::InitHits()
{
    // create Hits histograms in Hits subdir

    TH1F *h0 = new TH1F("hPreHitsEdep","Hits energy distribution in PRE(PMD)", 100, 0., 2000.); 
    h0->Sumw2() ;
    Add2HitsList(h0, 0) ;

    TH1F *h1 = new TH1F("hCpvHitsEdep","Hits energy distribution in CPV(PMD)", 100, 0., 2000.); 
    h1->Sumw2() ;
    Add2HitsList(h1, 1) ;

    TH1I *h2 = new TH1I("hPreHitsMult","Hits multiplicity distribution in PRE(PMD)", 500, 0, 10000) ; 
    h2->Sumw2() ;
    Add2HitsList(h2, 2) ;

    TH1I *h3 = new TH1I("hCpvHitsMult","Hits multiplicity distribution in PRE(PMD)", 500, 0, 10000) ; 
    h2->Sumw2() ;
    Add2HitsList(h3, 3) ;
}

//____________________________________________________________________________ 
void AliPMDQADataMakerSim::InitSDigits()
{
    // create SDigits histograms in SDigits subdir

    TH1F *h0 = new TH1F("hPreSDigitsEdep","SDigits energy distribution in PRE(PMD)", 100, 0., 2000.);
    h0->Sumw2();
    Add2SDigitsList(h0, 0);

    TH1F *h1 = new TH1F("hCpvSDigitsEdep","SDigits energy distribution in CPV(PMD)", 100, 0., 2000.);
    h1->Sumw2();
    Add2SDigitsList(h1, 1);

    TH1I *h2 = new TH1I("hPreSDigitsMult","SDigits multiplicity distribution in PRE(PMD)", 500, 0., 5000.);
    h2->Sumw2();
    Add2SDigitsList(h2, 2);

    TH1I *h3 = new TH1I("hCpvSDigitsMult","SDigits multiplicity distribution in CPV(PMD)", 500, 0., 5000.);
    h3->Sumw2();
    Add2SDigitsList(h3, 3);

}

//____________________________________________________________________________
void AliPMDQADataMakerSim::InitDigits()
{
    // create Digits histograms in Digits subdir

    TH1F *h0 = new TH1F("hPreDigitsEdep","Digits energy distribution in PRE(PMD)", 100, 0., 2000.);
    h0->Sumw2();
    Add2DigitsList(h0, 0);

    TH1F *h1 = new TH1F("hCpvDigitsEdep","Digits energy distribution in CPV(PMD)", 100, 0., 2000.); 
    h1->Sumw2();
    Add2DigitsList(h1, 1);

    TH1I *h2 = new TH1I("hPreDigitsMult","Digits multiplicity distribution in PRE(PMD)", 500, 0, 5000) ; 
    h2->Sumw2();
    Add2DigitsList(h2, 2);

    TH1I *h3 = new TH1I("hCpvDigitsMult","Digits multiplicity distribution in CPV(PMD)", 500, 0, 5000);
    h3->Sumw2();
    Add2DigitsList(h3, 3);

}

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

    Int_t premul = 0, cpvmul = 0;
  
    TIter next(hits); 
    AliPMDhit * hit; 

    while ( (hit = dynamic_cast<AliPMDhit *>(next())) )
    {
	if (hit->Z() > 361.5)
	{
	    GetHitsData(0)->Fill(hit->GetEnergy());
	    premul++;
	}
	else if (hit->Z() < 361.5)
	{
	    GetHitsData(1)->Fill(hit->GetEnergy());
	    cpvmul++;
	}
    }

    if(premul <= 0)
    {
	GetHitsData(2)->Fill(-1.); 
    }
    else
    {
	GetHitsData(2)->Fill(premul); 
    }

    if(cpvmul <= 0)
    {
	GetHitsData(3)->Fill(-1.); 
    }
    else
    {
	GetHitsData(3)->Fill(cpvmul); 
    }

}

//____________________________________________________________________________
void AliPMDQADataMakerSim::MakeHits(TTree * hitTree)
{
    // make QA data from Hit Tree

    TBranch * branch = hitTree->GetBranch("PMD") ;
    if ( ! branch )
    {
	AliWarning("PMD branch in Hit Tree not found") ;
	return;
    }

    static TClonesArray statichits("AliPMDhit", 1000);
    statichits.Clear();
    TClonesArray *hits = &statichits;
    static TClonesArray staticdummy("AliPMDhit", 1000);
    staticdummy.Clear();
    TClonesArray *dummy = &staticdummy;
    branch->SetAddress(&dummy);
    Int_t index = 0 ;  

    for (Int_t ientry = 0 ; ientry < branch->GetEntries() ; ientry++) {
	branch->GetEntry(ientry) ; 
	for (Int_t ihit = 0 ; ihit < dummy->GetEntries() ; ihit++) {
	    AliPMDhit * hit = dynamic_cast<AliPMDhit *> (dummy->At(ihit)) ; 
	    new((*hits)[index]) AliPMDhit(*hit) ; 

	    index++ ;
	} 
    } 	

    MakeHits(hits);

}
//____________________________________________________________________________
void AliPMDQADataMakerSim::MakeSDigits(TClonesArray * sdigits)
{
    // makes data from SDigits

    Int_t cpvmul = 0, premul = 0;

    TIter next(sdigits) ; 
    AliPMDsdigit * sdigit ; 
    while ( (sdigit = dynamic_cast<AliPMDsdigit *>(next())) )
    {
	if(sdigit->GetDetector() == 0)
	{
	    GetSDigitsData(0)->Fill( sdigit->GetCellEdep());
	    premul++;
	}
	if(sdigit->GetDetector() == 1)
	{
	    GetSDigitsData(1)->Fill( sdigit->GetCellEdep());
	    cpvmul++;
	}
	
    } 
    if (premul > 0) GetSDigitsData(2)->Fill(premul);
    if (cpvmul > 0) GetSDigitsData(3)->Fill(cpvmul);
    
}

//____________________________________________________________________________
void AliPMDQADataMakerSim::MakeSDigits(TTree * sdigitTree)
{
    // makes data from SDigit Tree

    TClonesArray * sdigits = new TClonesArray("AliPMDsdigit", 1000) ; 
    
    TBranch * branch = sdigitTree->GetBranch("PMDSDigit") ;
    branch->SetAddress(&sdigits) ;

    if ( ! branch )
    {
	AliWarning("PMD branch in SDigit Tree not found") ; 
    }
    else
    {
	for (Int_t ient = 0; ient < branch->GetEntries(); ient++)
	{
	    branch->GetEntry(ient) ;
	    MakeSDigits(sdigits) ; 
	}
    }
}

//____________________________________________________________________________
void AliPMDQADataMakerSim::MakeDigits(TClonesArray * digits)
{
    // makes data from Digits
    
    Int_t cpvmul = 0, premul = 0;

    TIter next(digits) ; 
    AliPMDdigit * digit ; 
    while ( (digit = dynamic_cast<AliPMDdigit *>(next())) )
    {
	if(digit->GetDetector() == 0)
	{
	    GetDigitsData(0)->Fill( digit->GetADC()) ;
	    premul++;
	}
	if(digit->GetDetector() == 1)
	{
	    GetDigitsData(1)->Fill( digit->GetADC());
	    cpvmul++;
	}
    }  

    if (premul > 0) GetDigitsData(2)->Fill(premul);
    if (cpvmul > 0) GetDigitsData(3)->Fill(cpvmul);


}

//____________________________________________________________________________
void AliPMDQADataMakerSim::MakeDigits(TTree * digitTree)
{
    // makes data from Digit Tree

    TClonesArray * digits = new TClonesArray("AliPMDdigit", 1000) ; 
    
    TBranch * branch = digitTree->GetBranch("PMDDigit") ;
    branch->SetAddress(&digits) ;

    if ( ! branch )
    {
	AliWarning("PMD branch in Digit Tree not found") ; 
    }
    else
    {
	for (Int_t ient = 0; ient < branch->GetEntries(); ient++)
	{
	    
	    branch->GetEntry(ient) ; 
	    MakeDigits(digits) ; 
	}
	
    }
}


//____________________________________________________________________________ 

void AliPMDQADataMakerSim::StartOfDetectorCycle()
{
    //Detector specific actions at start of cycle
    
}
//____________________________________________________________________________ 

void AliPMDQADataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX_t task, TObjArray * list)
{
    //Detector specific actions at end of cycle
    // do the QA checking
    AliQAChecker::Instance()->Run(AliQA::kPMD, task, list) ;  
}
Commit	Line	Data
a9523c3e	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	/*
	18	Produces the data needed to calculate the quality assurance.
	19	All data must be mergeable objects.
	20	B.K. Nandi
	21	*/
	22
	23	// --- ROOT system ---
	24	#include <TClonesArray.h>
	25	#include <TFile.h>
	26	#include <TH1F.h>
	27	#include <TH1I.h>
	28	#include <TH2F.h>
	29	#include <TTree.h>
	30
	31	// --- Standard library ---
	32
	33	// --- AliRoot header files ---
	34
	35	#include "AliLog.h"
	36	#include "AliPMDhit.h"
	37	#include "AliPMDsdigit.h"
	38	#include "AliPMDdigit.h"
	39	#include "AliPMDQADataMakerSim.h"
	40	#include "AliQAChecker.h"
	41
	42	ClassImp(AliPMDQADataMakerSim)
	43
	44	//____________________________________________________________________________
	45	AliPMDQADataMakerSim::AliPMDQADataMakerSim() :
	46	AliQADataMakerSim(AliQA::GetDetName(AliQA::kPMD), "PMD Quality Assurance Data Maker")
	47	{
	48	// ctor
	49	}
	50
	51	//____________________________________________________________________________
	52	AliPMDQADataMakerSim::AliPMDQADataMakerSim(const AliPMDQADataMakerSim& qadm) :
	53	AliQADataMakerSim()
	54	{
	55	//copy ctor
	56	SetName((const char*)qadm.GetName()) ;
	57	SetTitle((const char*)qadm.GetTitle());
	58	}
	59
	60	//__________________________________________________________________
	61	AliPMDQADataMakerSim& AliPMDQADataMakerSim::operator = (const AliPMDQADataMakerSim& qadm )
	62	{
	63	// Assign operator.
	64	this->~AliPMDQADataMakerSim();
65	new(this) AliPMDQADataMakerSim(qadm);
66	return *this;
67	}
68
69	//____________________________________________________________________________
70	void AliPMDQADataMakerSim::InitHits()
71	{
72	// create Hits histograms in Hits subdir
73
74	TH1F *h0 = new TH1F("hPreHitsEdep","Hits energy distribution in PRE(PMD)", 100, 0., 2000.);
75	h0->Sumw2() ;
76	Add2HitsList(h0, 0) ;
77
78	TH1F *h1 = new TH1F("hCpvHitsEdep","Hits energy distribution in CPV(PMD)", 100, 0., 2000.);
79	h1->Sumw2() ;
80	Add2HitsList(h1, 1) ;
81
82	TH1I *h2 = new TH1I("hPreHitsMult","Hits multiplicity distribution in PRE(PMD)", 500, 0, 10000) ;
83	h2->Sumw2() ;
84	Add2HitsList(h2, 2) ;
85
86	TH1I *h3 = new TH1I("hCpvHitsMult","Hits multiplicity distribution in PRE(PMD)", 500, 0, 10000) ;
87	h2->Sumw2() ;
88	Add2HitsList(h3, 3) ;
89	}
90
91	//____________________________________________________________________________
92	void AliPMDQADataMakerSim::InitSDigits()
93	{
94	// create SDigits histograms in SDigits subdir
95
96	TH1F *h0 = new TH1F("hPreSDigitsEdep","SDigits energy distribution in PRE(PMD)", 100, 0., 2000.);
97	h0->Sumw2();
98	Add2SDigitsList(h0, 0);
99
100	TH1F *h1 = new TH1F("hCpvSDigitsEdep","SDigits energy distribution in CPV(PMD)", 100, 0., 2000.);
101	h1->Sumw2();
102	Add2SDigitsList(h1, 1);
103
104	TH1I *h2 = new TH1I("hPreSDigitsMult","SDigits multiplicity distribution in PRE(PMD)", 500, 0., 5000.);
105	h2->Sumw2();
106	Add2SDigitsList(h2, 2);
107
108	TH1I *h3 = new TH1I("hCpvSDigitsMult","SDigits multiplicity distribution in CPV(PMD)", 500, 0., 5000.);
109	h3->Sumw2();
110	Add2SDigitsList(h3, 3);
111
112	}
113
114	//____________________________________________________________________________
115	void AliPMDQADataMakerSim::InitDigits()
116	{
117	// create Digits histograms in Digits subdir
118
119	TH1F *h0 = new TH1F("hPreDigitsEdep","Digits energy distribution in PRE(PMD)", 100, 0., 2000.);
120	h0->Sumw2();
121	Add2DigitsList(h0, 0);
122
123	TH1F *h1 = new TH1F("hCpvDigitsEdep","Digits energy distribution in CPV(PMD)", 100, 0., 2000.);
124	h1->Sumw2();
125	Add2DigitsList(h1, 1);
126
127	TH1I *h2 = new TH1I("hPreDigitsMult","Digits multiplicity distribution in PRE(PMD)", 500, 0, 5000) ;
128	h2->Sumw2();
129	Add2DigitsList(h2, 2);
130
131	TH1I *h3 = new TH1I("hCpvDigitsMult","Digits multiplicity distribution in CPV(PMD)", 500, 0, 5000);
132	h3->Sumw2();
133	Add2DigitsList(h3, 3);
134
135	}
136
137	//____________________________________________________________________________
138	void AliPMDQADataMakerSim::MakeHits(TClonesArray *hits)
139	{
140	//make QA data from Hits
141
142	Int_t premul = 0, cpvmul = 0;
143
144	TIter next(hits);
145	AliPMDhit * hit;
146
147	while ( (hit = dynamic_cast<AliPMDhit *>(next())) )
148	{
149	if (hit->Z() > 361.5)
150	{
151	GetHitsData(0)->Fill(hit->GetEnergy());
152	premul++;
153	}
154	else if (hit->Z() < 361.5)
155	{
156	GetHitsData(1)->Fill(hit->GetEnergy());
157	cpvmul++;
158	}
159	}
160
161	if(premul <= 0)
162	{
163	GetHitsData(2)->Fill(-1.);
164	}
165	else
166	{
167	GetHitsData(2)->Fill(premul);
168	}
169
170	if(cpvmul <= 0)
171	{
172	GetHitsData(3)->Fill(-1.);
173	}
174	else
175	{
176	GetHitsData(3)->Fill(cpvmul);
177	}
178
179	}
180
181	//____________________________________________________________________________
182	void AliPMDQADataMakerSim::MakeHits(TTree * hitTree)
183	{
184	// make QA data from Hit Tree
185
186	TBranch * branch = hitTree->GetBranch("PMD") ;
187	if ( ! branch )
188	{
189	AliWarning("PMD branch in Hit Tree not found") ;
190	return;
191	}
192
193	static TClonesArray statichits("AliPMDhit", 1000);
194	statichits.Clear();
195	TClonesArray *hits = &statichits;
196	static TClonesArray staticdummy("AliPMDhit", 1000);
197	staticdummy.Clear();
198	TClonesArray *dummy = &staticdummy;
199	branch->SetAddress(&dummy);
200	Int_t index = 0 ;
201
202	for (Int_t ientry = 0 ; ientry < branch->GetEntries() ; ientry++) {
203	branch->GetEntry(ientry) ;
204	for (Int_t ihit = 0 ; ihit < dummy->GetEntries() ; ihit++) {
205	AliPMDhit * hit = dynamic_cast<AliPMDhit *> (dummy->At(ihit)) ;
206	new((hits)[index]) AliPMDhit(hit) ;
207
208	index++ ;
209	}
210	}
211
212	MakeHits(hits);
213
214	}
215	//____________________________________________________________________________
216	void AliPMDQADataMakerSim::MakeSDigits(TClonesArray * sdigits)
217	{
218	// makes data from SDigits
219
220	Int_t cpvmul = 0, premul = 0;
221
222	TIter next(sdigits) ;
223	AliPMDsdigit * sdigit ;
224	while ( (sdigit = dynamic_cast<AliPMDsdigit *>(next())) )
225	{
226	if(sdigit->GetDetector() == 0)
227	{
228	GetSDigitsData(0)->Fill( sdigit->GetCellEdep());
229	premul++;
230	}
231	if(sdigit->GetDetector() == 1)
232	{
233	GetSDigitsData(1)->Fill( sdigit->GetCellEdep());
234	cpvmul++;
235	}
236
237	}
238	if (premul > 0) GetSDigitsData(2)->Fill(premul);
239	if (cpvmul > 0) GetSDigitsData(3)->Fill(cpvmul);
240
241	}
242
243	//____________________________________________________________________________
244	void AliPMDQADataMakerSim::MakeSDigits(TTree * sdigitTree)
245	{
246	// makes data from SDigit Tree
247
248	TClonesArray * sdigits = new TClonesArray("AliPMDsdigit", 1000) ;
249
250	TBranch * branch = sdigitTree->GetBranch("PMDSDigit") ;
251	branch->SetAddress(&sdigits) ;
252
253	if ( ! branch )
254	{
255	AliWarning("PMD branch in SDigit Tree not found") ;
256	}
257	else
258	{
259	for (Int_t ient = 0; ient < branch->GetEntries(); ient++)
260	{
261	branch->GetEntry(ient) ;
262	MakeSDigits(sdigits) ;
263	}
264	}
265	}
266
267	//____________________________________________________________________________
268	void AliPMDQADataMakerSim::MakeDigits(TClonesArray * digits)
269	{
270	// makes data from Digits
271
272	Int_t cpvmul = 0, premul = 0;
273
274	TIter next(digits) ;
275	AliPMDdigit * digit ;
276	while ( (digit = dynamic_cast<AliPMDdigit *>(next())) )
277	{
278	if(digit->GetDetector() == 0)
279	{
280	GetDigitsData(0)->Fill( digit->GetADC()) ;
281	premul++;
282	}
283	if(digit->GetDetector() == 1)
284	{
285	GetDigitsData(1)->Fill( digit->GetADC());
286	cpvmul++;
287	}
288	}
289
290	if (premul > 0) GetDigitsData(2)->Fill(premul);
291	if (cpvmul > 0) GetDigitsData(3)->Fill(cpvmul);
292
293
294	}
295
296	//____________________________________________________________________________
297	void AliPMDQADataMakerSim::MakeDigits(TTree * digitTree)
298	{
299	// makes data from Digit Tree
300
301	TClonesArray * digits = new TClonesArray("AliPMDdigit", 1000) ;
302
303	TBranch * branch = digitTree->GetBranch("PMDDigit") ;
304	branch->SetAddress(&digits) ;
305
306	if ( ! branch )
307	{
308	AliWarning("PMD branch in Digit Tree not found") ;
309	}
310	else
311	{
312	for (Int_t ient = 0; ient < branch->GetEntries(); ient++)
313	{
314
315	branch->GetEntry(ient) ;
316	MakeDigits(digits) ;
317	}
318
319	}
320	}
321
322
323	//____________________________________________________________________________
324
325	void AliPMDQADataMakerSim::StartOfDetectorCycle()
326	{
327	//Detector specific actions at start of cycle
328
329	}
330	//____________________________________________________________________________
331
332	void AliPMDQADataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX_t task, TObjArray * list)
333	{
334	//Detector specific actions at end of cycle
335	// do the QA checking
336	AliQAChecker::Instance()->Run(AliQA::kPMD, task, list) ;
337	}