[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(AliQAv1::GetDetName(AliQAv1::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(AliQAv1::TASKINDEX_t task, TObjArray ** list)
{
  //Detector specific actions at end of cycle
  // do the QA checking
  AliQAChecker::Instance()->Run(AliQAv1::kPHOS, task, list) ;  
}

//____________________________________________________________________________ 
void AliPHOSQADataMakerSim::InitHits()
{
  // create Hits histograms in Hits subdir
  const Bool_t expert   = kTRUE ; 
  const Bool_t image    = kTRUE ; 
  TH1F * h0 = new TH1F("hPhosHits",    "Hits energy distribution in PHOS;Energy [MeV];Counts",       100, 0., 100.) ; 
  h0->Sumw2() ;
  Add2HitsList(h0, kHits, !expert, image) ;
  TH1I * h1 = new TH1I("hPhosHitsMul", "Hits multiplicity distribution in PHOS;# of Hits;Entries", 500, 0., 10000) ; 
  h1->Sumw2() ;
  Add2HitsList(h1, kHitsMul, !expert, image) ;
  
}

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

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

//____________________________________________________________________________
void AliPHOSQADataMakerSim::MakeHits()
{
  //make QA data from Hits
  
  // Check id histograms already created for this Event Specie
  if ( ! GetHitsData(kHits) )
    InitHits() ;
  
  
  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
  
  // Check id histograms already created for this Event Specie
  if ( ! GetDigitsData(kDigits) )
    InitDigits() ;
  
  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
  
  
  // Check id histograms already created for this Event Specie
  if ( ! GetSDigitsData(kSDigits) )
    InitSDigits() ;

  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() :
4e25ac79	47	AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kPHOS), "PHOS Quality Assurance Data Maker"),
8d8258f6	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	//____________________________________________________________________________
4e25ac79	75	void AliPHOSQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray ** list)
04236e67	76	{
	77	//Detector specific actions at end of cycle
	78	// do the QA checking
4e25ac79	79	AliQAChecker::Instance()->Run(AliQAv1::kPHOS, task, list) ;
04236e67	80	}
	81
	82	//____________________________________________________________________________
	83	void AliPHOSQADataMakerSim::InitHits()
	84	{
	85	// create Hits histograms in Hits subdir
7d297381	86	const Bool_t expert = kTRUE ;
7d297381	87	const Bool_t image = kTRUE ;
db72ff3b	88	TH1F * h0 = new TH1F("hPhosHits", "Hits energy distribution in PHOS;Energy [MeV];Counts", 100, 0., 100.) ;
04236e67	89	h0->Sumw2() ;
7d297381	90	Add2HitsList(h0, kHits, !expert, image) ;
db72ff3b	91	TH1I * h1 = new TH1I("hPhosHitsMul", "Hits multiplicity distribution in PHOS;# of Hits;Entries", 500, 0., 10000) ;
04236e67	92	h1->Sumw2() ;
7d297381	93	Add2HitsList(h1, kHitsMul, !expert, image) ;
57acd2d2	94
04236e67	95	}
	96
	97	//____________________________________________________________________________
	98	void AliPHOSQADataMakerSim::InitDigits()
	99	{
	100	// create Digits histograms in Digits subdir
7d297381	101	const Bool_t expert = kTRUE ;
7d297381	102	const Bool_t image = kTRUE ;
db72ff3b	103	TH1I * h0 = new TH1I("hPhosDigits", "Digits amplitude distribution in PHOS;Amplitude [ADC counts];Counts", 500, 0, 1000) ;
04236e67	104	h0->Sumw2() ;
7d297381	105	Add2DigitsList(h0, kDigits, !expert, image) ;
db72ff3b	106	TH1I * h1 = new TH1I("hPhosDigitsMul", "Digits multiplicity distribution in PHOS;# of Digits;Entries", 2000, 0, 10000) ;
04236e67	107	h1->Sumw2() ;
7d297381	108	Add2DigitsList(h1, kDigitsMul, !expert, image) ;
04236e67	109	}
	110
	111	//____________________________________________________________________________
	112	void AliPHOSQADataMakerSim::InitSDigits()
	113	{
	114	// create SDigits histograms in SDigits subdir
7d297381	115	const Bool_t expert = kTRUE ;
7d297381	116	const Bool_t image = kTRUE ;
db72ff3b	117	TH1F * h0 = new TH1F("hPhosSDigits", "SDigits energy distribution in PHOS; Energy [MeV];Counts", 500, 0., 1000.) ;
04236e67	118	h0->Sumw2() ;
7d297381	119	Add2SDigitsList(h0, kSDigits, !expert, image) ;
db72ff3b	120	TH1I * h1 = new TH1I("hPhosSDigitsMul", "SDigits multiplicity distribution in PHOS;# of SDigits;Entries", 500, 0, 1000) ;
04236e67	121	h1->Sumw2() ;
7d297381	122	Add2SDigitsList(h1, kSDigitsMul, !expert, image) ;
04236e67	123	}
	124
	125	//____________________________________________________________________________
8d8258f6	126	void AliPHOSQADataMakerSim::MakeHits()
04236e67	127	{
8d8258f6	128	//make QA data from Hits
8d8258f6	129
eca4fa66	130	// Check id histograms already created for this Event Specie
	131	if ( ! GetHitsData(kHits) )
	132	InitHits() ;
	133
	134
8d8258f6	135	TIter next(fHits) ;
	136	AliPHOSHit * hit ;
	137	while ( (hit = dynamic_cast<AliPHOSHit *>(next())) ) {
57acd2d2	138	GetHitsData(kHits)->Fill( hit->GetEnergy()) ;
8d8258f6	139	}
04236e67	140	}
	141
	142	//____________________________________________________________________________
	143	void AliPHOSQADataMakerSim::MakeHits(TTree * hitTree)
	144	{
8d8258f6	145	// make QA data from Hit Tree
	146
	147	TBranch * branch = hitTree->GetBranch("PHOS") ;
	148	if ( ! branch ) {
	149	AliWarning("PHOS branch in Hit Tree not found") ;
	150	} else {
	151	Int_t nHits = 0;
	152	branch->SetAddress(&fHits) ;
	153	for (Int_t ientry = 0 ; ientry < branch->GetEntries() ; ientry++) {
	154	branch->GetEntry(ientry) ;
	155	nHits += fHits->GetEntriesFast();
	156	MakeHits() ;
	157	fHits->Clear();
	158	}
	159	GetHitsData(1)->Fill(nHits) ;
	160	}
04236e67	161	}
	162
	163	//____________________________________________________________________________
	164	void AliPHOSQADataMakerSim::MakeDigits(TClonesArray * digits)
	165	{
	166	// makes data from Digits
eca4fa66	167
	168	// Check id histograms already created for this Event Specie
	169	if ( ! GetDigitsData(kDigits) )
	170	InitDigits() ;
	171
	172	GetDigitsData(1)->Fill(digits->GetEntriesFast()) ;
04236e67	173	TIter next(digits) ;
	174	AliPHOSDigit * digit ;
	175	while ( (digit = dynamic_cast<AliPHOSDigit *>(next())) ) {
57acd2d2	176	GetDigitsData(kDigits)->Fill( digit->GetEnergy()) ;
04236e67	177	}
	178	}
	179
	180	//____________________________________________________________________________
	181	void AliPHOSQADataMakerSim::MakeDigits(TTree * digitTree)
	182	{
	183	// makes data from Digit Tree
	184	TClonesArray * digits = new TClonesArray("AliPHOSDigit", 1000) ;
	185
	186	TBranch * branch = digitTree->GetBranch("PHOS") ;
	187	if ( ! branch ) {
	188	AliWarning("PHOS branch in Digit Tree not found") ;
	189	} else {
	190	branch->SetAddress(&digits) ;
	191	branch->GetEntry(0) ;
	192	MakeDigits(digits) ;
	193	}
	194	}
	195
	196	//____________________________________________________________________________
	197	void AliPHOSQADataMakerSim::MakeSDigits(TClonesArray * sdigits)
	198	{
	199	// makes data from SDigits
	200
eca4fa66	201
	202	// Check id histograms already created for this Event Specie
	203	if ( ! GetSDigitsData(kSDigits) )
	204	InitSDigits() ;
	205
	206	GetSDigitsData(1)->Fill(sdigits->GetEntriesFast()) ;
	207	TIter next(sdigits) ;
	208	AliPHOSDigit * sdigit ;
	209	while ( (sdigit = dynamic_cast<AliPHOSDigit *>(next())) ) {
	210	GetSDigitsData(kSDigits)->Fill( sdigit->GetEnergy()) ;
	211	}
04236e67	212	}
	213
	214	//____________________________________________________________________________
	215	void AliPHOSQADataMakerSim::MakeSDigits(TTree * sdigitTree)
	216	{
	217	// makes data from SDigit Tree
	218	TClonesArray * sdigits = new TClonesArray("AliPHOSDigit", 1000) ;
	219
	220	TBranch * branch = sdigitTree->GetBranch("PHOS") ;
	221	if ( ! branch ) {
	222	AliWarning("PHOS branch in SDigit Tree not found") ;
	223	} else {
	224	branch->SetAddress(&sdigits) ;
	225	branch->GetEntry(0) ;
	226	MakeSDigits(sdigits) ;
	227	}
	228	}
	229
	230	//____________________________________________________________________________
	231	void AliPHOSQADataMakerSim::StartOfDetectorCycle()
	232	{
	233	//Detector specific actions at start of cycle
	234
	235	}