[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 (keV)PRE(PMD)", 500, 0., 500.); 
    h0->Sumw2() ;
    Add2HitsList(h0, 0) ;

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

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

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

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

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

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

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

    TH1I *h3 = new TH1I("hCpvSDigitsMult","SDigits multiplicity distribution in CPV(PMD)", 500, 0., 1000.);
    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, 1000) ; 
    h2->Sumw2();
    Add2DigitsList(h2, 2);

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

}

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

    Int_t premul = 0, cpvmul = 0;
    Float_t edepkev = 0.;
    TIter next(hits); 
    AliPMDhit * hit; 
    
    while ( (hit = dynamic_cast<AliPMDhit *>(next())) )
      {
	if (hit->Z() > 361.5)
	  {
	    edepkev = (hit->GetEnergy())/1000.;
	    GetHitsData(0)->Fill(edepkev);
	    premul++;
	  }
	else if (hit->Z() < 361.5)
	  {
	    edepkev = (hit->GetEnergy())/1000.;
	    GetHitsData(1)->Fill(edepkev);
	    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;
    Float_t edepkev = 0.;

    TIter next(sdigits) ; 
    AliPMDsdigit * sdigit ; 
    while ( (sdigit = dynamic_cast<AliPMDsdigit *>(next())) )
    {
	if(sdigit->GetDetector() == 0)
	{
	  edepkev = (sdigit->GetCellEdep())/1000.;
	  GetSDigitsData(0)->Fill(edepkev);
	  premul++;
	}
	if(sdigit->GetDetector() == 1)
	{
	  edepkev = (sdigit->GetCellEdep())/1000.;
	  GetSDigitsData(1)->Fill(edepkev);
	  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
2c1131dd	74	TH1F *h0 = new TH1F("hPreHitsEdep","Hits energy distribution in (keV)PRE(PMD)", 500, 0., 500.);
a9523c3e	75	h0->Sumw2() ;
	76	Add2HitsList(h0, 0) ;
	77
2c1131dd	78	TH1F *h1 = new TH1F("hCpvHitsEdep","Hits energy distribution in (keV)CPV(PMD)", 500, 0., 500.);
a9523c3e	79	h1->Sumw2() ;
	80	Add2HitsList(h1, 1) ;
	81
2c1131dd	82	TH1I *h2 = new TH1I("hPreHitsMult","Hits multiplicity distribution in PRE(PMD)", 500, 0, 3000) ;
a9523c3e	83	h2->Sumw2() ;
	84	Add2HitsList(h2, 2) ;
	85
2c1131dd	86	TH1I *h3 = new TH1I("hCpvHitsMult","Hits multiplicity distribution in PRE(PMD)", 500, 0, 3000) ;
a9523c3e	87	h2->Sumw2() ;
	88	Add2HitsList(h3, 3) ;
	89	}
	90
	91	//____________________________________________________________________________
	92	void AliPMDQADataMakerSim::InitSDigits()
	93	{
	94	// create SDigits histograms in SDigits subdir
	95
2c1131dd	96	TH1F *h0 = new TH1F("hPreSDigitsEdep","SDigits energy distribution in(keV) PRE(PMD)", 500, 0., 500.);
a9523c3e	97	h0->Sumw2();
	98	Add2SDigitsList(h0, 0);
	99
2c1131dd	100	TH1F *h1 = new TH1F("hCpvSDigitsEdep","SDigits energy distribution in (keV)CPV(PMD)", 500, 0., 500.);
a9523c3e	101	h1->Sumw2();
	102	Add2SDigitsList(h1, 1);
	103
2c1131dd	104	TH1I *h2 = new TH1I("hPreSDigitsMult","SDigits multiplicity distribution in PRE(PMD)", 500, 0., 1000.);
a9523c3e	105	h2->Sumw2();
	106	Add2SDigitsList(h2, 2);
	107
2c1131dd	108	TH1I *h3 = new TH1I("hCpvSDigitsMult","SDigits multiplicity distribution in CPV(PMD)", 500, 0., 1000.);
a9523c3e	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
2c1131dd	127	TH1I *h2 = new TH1I("hPreDigitsMult","Digits multiplicity distribution in PRE(PMD)", 500, 0, 1000) ;
a9523c3e	128	h2->Sumw2();
	129	Add2DigitsList(h2, 2);
	130
2c1131dd	131	TH1I *h3 = new TH1I("hCpvDigitsMult","Digits multiplicity distribution in CPV(PMD)", 500, 0, 1000);
a9523c3e	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;
2c1131dd	143	Float_t edepkev = 0.;
a9523c3e	144	TIter next(hits);
a9523c3e	145	AliPMDhit * hit;
2c1131dd	146
a9523c3e	147	while ( (hit = dynamic_cast<AliPMDhit *>(next())) )
2c1131dd	148	{
a9523c3e	149	if (hit->Z() > 361.5)
2c1131dd	150	{
	151	edepkev = (hit->GetEnergy())/1000.;
	152	GetHitsData(0)->Fill(edepkev);
a9523c3e	153	premul++;
2c1131dd	154	}
a9523c3e	155	else if (hit->Z() < 361.5)
2c1131dd	156	{
	157	edepkev = (hit->GetEnergy())/1000.;
	158	GetHitsData(1)->Fill(edepkev);
a9523c3e	159	cpvmul++;
2c1131dd	160	}
a9523c3e	161	}
	162
	163	if(premul <= 0)
	164	{
	165	GetHitsData(2)->Fill(-1.);
	166	}
	167	else
	168	{
	169	GetHitsData(2)->Fill(premul);
	170	}
	171
	172	if(cpvmul <= 0)
	173	{
	174	GetHitsData(3)->Fill(-1.);
	175	}
	176	else
	177	{
	178	GetHitsData(3)->Fill(cpvmul);
	179	}
	180
	181	}
	182
	183	//____________________________________________________________________________
	184	void AliPMDQADataMakerSim::MakeHits(TTree * hitTree)
	185	{
	186	// make QA data from Hit Tree
	187
	188	TBranch * branch = hitTree->GetBranch("PMD") ;
	189	if ( ! branch )
	190	{
	191	AliWarning("PMD branch in Hit Tree not found") ;
	192	return;
	193	}
	194
	195	static TClonesArray statichits("AliPMDhit", 1000);
	196	statichits.Clear();
	197	TClonesArray *hits = &statichits;
	198	static TClonesArray staticdummy("AliPMDhit", 1000);
	199	staticdummy.Clear();
	200	TClonesArray *dummy = &staticdummy;
	201	branch->SetAddress(&dummy);
	202	Int_t index = 0 ;
	203
	204	for (Int_t ientry = 0 ; ientry < branch->GetEntries() ; ientry++) {
	205	branch->GetEntry(ientry) ;
	206	for (Int_t ihit = 0 ; ihit < dummy->GetEntries() ; ihit++) {
	207	AliPMDhit * hit = dynamic_cast<AliPMDhit *> (dummy->At(ihit)) ;
	208	new((hits)[index]) AliPMDhit(hit) ;
	209
	210	index++ ;
	211	}
	212	}
	213
	214	MakeHits(hits);
	215
	216	}
	217	//____________________________________________________________________________
	218	void AliPMDQADataMakerSim::MakeSDigits(TClonesArray * sdigits)
	219	{
	220	// makes data from SDigits
	221
	222	Int_t cpvmul = 0, premul = 0;
2c1131dd	223	Float_t edepkev = 0.;
a9523c3e	224
	225	TIter next(sdigits) ;
	226	AliPMDsdigit * sdigit ;
	227	while ( (sdigit = dynamic_cast<AliPMDsdigit *>(next())) )
	228	{
	229	if(sdigit->GetDetector() == 0)
	230	{
2c1131dd	231	edepkev = (sdigit->GetCellEdep())/1000.;
	232	GetSDigitsData(0)->Fill(edepkev);
	233	premul++;
a9523c3e	234	}
	235	if(sdigit->GetDetector() == 1)
	236	{
2c1131dd	237	edepkev = (sdigit->GetCellEdep())/1000.;
	238	GetSDigitsData(1)->Fill(edepkev);
	239	cpvmul++;
a9523c3e	240	}
	241
	242	}
	243	if (premul > 0) GetSDigitsData(2)->Fill(premul);
	244	if (cpvmul > 0) GetSDigitsData(3)->Fill(cpvmul);
	245
	246	}
	247
	248	//____________________________________________________________________________
	249	void AliPMDQADataMakerSim::MakeSDigits(TTree * sdigitTree)
	250	{
	251	// makes data from SDigit Tree
	252
	253	TClonesArray * sdigits = new TClonesArray("AliPMDsdigit", 1000) ;
	254
	255	TBranch * branch = sdigitTree->GetBranch("PMDSDigit") ;
	256	branch->SetAddress(&sdigits) ;
	257
	258	if ( ! branch )
	259	{
	260	AliWarning("PMD branch in SDigit Tree not found") ;
	261	}
	262	else
	263	{
	264	for (Int_t ient = 0; ient < branch->GetEntries(); ient++)
	265	{
	266	branch->GetEntry(ient) ;
	267	MakeSDigits(sdigits) ;
	268	}
	269	}
	270	}
	271
	272	//____________________________________________________________________________
	273	void AliPMDQADataMakerSim::MakeDigits(TClonesArray * digits)
	274	{
	275	// makes data from Digits
	276
	277	Int_t cpvmul = 0, premul = 0;
	278
	279	TIter next(digits) ;
	280	AliPMDdigit * digit ;
	281	while ( (digit = dynamic_cast<AliPMDdigit *>(next())) )
	282	{
	283	if(digit->GetDetector() == 0)
	284	{
	285	GetDigitsData(0)->Fill( digit->GetADC()) ;
	286	premul++;
	287	}
	288	if(digit->GetDetector() == 1)
	289	{
	290	GetDigitsData(1)->Fill( digit->GetADC());
	291	cpvmul++;
	292	}
	293	}
	294
	295	if (premul > 0) GetDigitsData(2)->Fill(premul);
	296	if (cpvmul > 0) GetDigitsData(3)->Fill(cpvmul);
	297
	298
	299	}
	300
	301	//____________________________________________________________________________
	302	void AliPMDQADataMakerSim::MakeDigits(TTree * digitTree)
	303	{
304	// makes data from Digit Tree
305
306	TClonesArray * digits = new TClonesArray("AliPMDdigit", 1000) ;
307
308	TBranch * branch = digitTree->GetBranch("PMDDigit") ;
309	branch->SetAddress(&digits) ;
310
311	if ( ! branch )
312	{
313	AliWarning("PMD branch in Digit Tree not found") ;
314	}
315	else
316	{
317	for (Int_t ient = 0; ient < branch->GetEntries(); ient++)
318	{
a9523c3e	319	branch->GetEntry(ient) ;
	320	MakeDigits(digits) ;
	321	}
	322
	323	}
	324	}
	325
	326
	327	//____________________________________________________________________________
	328
	329	void AliPMDQADataMakerSim::StartOfDetectorCycle()
	330	{
	331	//Detector specific actions at start of cycle
	332
	333	}
	334	//____________________________________________________________________________
	335
	336	void AliPMDQADataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX_t task, TObjArray * list)
	337	{
	338	//Detector specific actions at end of cycle
	339	// do the QA checking
	340	AliQAChecker::Instance()->Run(AliQA::kPMD, task, list) ;
	341	}