[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
  Bool_t expert   = kTRUE ; 
  TH1F * h0 = new TH1F("hPhosHits",    "Hits energy distribution in PHOS",       100, 0., 100.) ; 
  h0->Sumw2() ;
  Add2HitsList(h0, kHits, !expert) ;
  TH1I * h1 = new TH1I("hPhosHitsMul", "Hits multiplicity distribution in PHOS", 500, 0., 10000) ; 
  h1->Sumw2() ;
  Add2HitsList(h1, kHitsMul, !expert) ;
  
}

//____________________________________________________________________________ 
void AliPHOSQADataMakerSim::InitDigits()
{
  // create Digits histograms in Digits subdir
  Bool_t expert   = kTRUE ; 
  TH1I * h0 = new TH1I("hPhosDigits",    "Digits amplitude distribution in PHOS",    500, 0, 1000) ; 
  h0->Sumw2() ;
  Add2DigitsList(h0, kDigits, !expert) ;
  TH1I * h1 = new TH1I("hPhosDigitsMul", "Digits multiplicity distribution in PHOS", 2000, 0, 10000) ; 
  h1->Sumw2() ;
  Add2DigitsList(h1, kDigitsMul, !expert) ;
}

//____________________________________________________________________________ 
void AliPHOSQADataMakerSim::InitSDigits()
{
  // create SDigits histograms in SDigits subdir
  Bool_t expert   = kTRUE ; 
  TH1F * h0 = new TH1F("hPhosSDigits",    "SDigits energy distribution in PHOS",       500, 0., 1000.) ; 
  h0->Sumw2() ;
  Add2SDigitsList(h0, kSDigits, !expert) ;
  TH1I * h1 = new TH1I("hPhosSDigitsMul", "SDigits multiplicity distribution in PHOS", 500, 0,  1000) ; 
  h1->Sumw2() ;
  Add2SDigitsList(h1, kSDigitsMul, !expert) ;
}

//____________________________________________________________________________
void AliPHOSQADataMakerSim::MakeHits()
{
  //make QA data from Hits
  
  TIter next(fHits) ; 
  AliPHOSHit * hit ; 
  while ( (hit = dynamic_cast<AliPHOSHit *>(next())) ) {
    GetHitsData(kHits)->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(kDigits)->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(kSDigits)->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	//____________________________________________________________________________
57acd2d2	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
33055c0e	86	Bool_t expert = kTRUE ;
04236e67	87	TH1F * h0 = new TH1F("hPhosHits", "Hits energy distribution in PHOS", 100, 0., 100.) ;
04236e67	88	h0->Sumw2() ;
57acd2d2	89	Add2HitsList(h0, kHits, !expert) ;
8d8258f6	90	TH1I * h1 = new TH1I("hPhosHitsMul", "Hits multiplicity distribution in PHOS", 500, 0., 10000) ;
04236e67	91	h1->Sumw2() ;
57acd2d2	92	Add2HitsList(h1, kHitsMul, !expert) ;
57acd2d2	93
04236e67	94	}
	95
	96	//____________________________________________________________________________
	97	void AliPHOSQADataMakerSim::InitDigits()
	98	{
	99	// create Digits histograms in Digits subdir
33055c0e	100	Bool_t expert = kTRUE ;
5997244b	101	TH1I * h0 = new TH1I("hPhosDigits", "Digits amplitude distribution in PHOS", 500, 0, 1000) ;
04236e67	102	h0->Sumw2() ;
57acd2d2	103	Add2DigitsList(h0, kDigits, !expert) ;
5997244b	104	TH1I * h1 = new TH1I("hPhosDigitsMul", "Digits multiplicity distribution in PHOS", 2000, 0, 10000) ;
04236e67	105	h1->Sumw2() ;
57acd2d2	106	Add2DigitsList(h1, kDigitsMul, !expert) ;
04236e67	107	}
	108
	109	//____________________________________________________________________________
	110	void AliPHOSQADataMakerSim::InitSDigits()
	111	{
	112	// create SDigits histograms in SDigits subdir
33055c0e	113	Bool_t expert = kTRUE ;
5997244b	114	TH1F * h0 = new TH1F("hPhosSDigits", "SDigits energy distribution in PHOS", 500, 0., 1000.) ;
04236e67	115	h0->Sumw2() ;
57acd2d2	116	Add2SDigitsList(h0, kSDigits, !expert) ;
5997244b	117	TH1I * h1 = new TH1I("hPhosSDigitsMul", "SDigits multiplicity distribution in PHOS", 500, 0, 1000) ;
04236e67	118	h1->Sumw2() ;
57acd2d2	119	Add2SDigitsList(h1, kSDigitsMul, !expert) ;
04236e67	120	}
	121
	122	//____________________________________________________________________________
8d8258f6	123	void AliPHOSQADataMakerSim::MakeHits()
04236e67	124	{
8d8258f6	125	//make QA data from Hits
	126
	127	TIter next(fHits) ;
	128	AliPHOSHit * hit ;
	129	while ( (hit = dynamic_cast<AliPHOSHit *>(next())) ) {
57acd2d2	130	GetHitsData(kHits)->Fill( hit->GetEnergy()) ;
8d8258f6	131	}
04236e67	132	}
	133
	134	//____________________________________________________________________________
	135	void AliPHOSQADataMakerSim::MakeHits(TTree * hitTree)
	136	{
8d8258f6	137	// make QA data from Hit Tree
	138
	139	TBranch * branch = hitTree->GetBranch("PHOS") ;
	140	if ( ! branch ) {
	141	AliWarning("PHOS branch in Hit Tree not found") ;
	142	} else {
	143	Int_t nHits = 0;
	144	branch->SetAddress(&fHits) ;
	145	for (Int_t ientry = 0 ; ientry < branch->GetEntries() ; ientry++) {
	146	branch->GetEntry(ientry) ;
	147	nHits += fHits->GetEntriesFast();
	148	MakeHits() ;
	149	fHits->Clear();
	150	}
	151	GetHitsData(1)->Fill(nHits) ;
	152	}
04236e67	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())) ) {
57acd2d2	164	GetDigitsData(kDigits)->Fill( digit->GetEnergy()) ;
04236e67	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())) ) {
57acd2d2	193	GetSDigitsData(kSDigits)->Fill( sdigit->GetEnergy()) ;
04236e67	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	}