[u/mrichter/AliRoot.git] / PHOS / AliPHOSQualAssDataMaker.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> 

// --- Standard library ---

// --- AliRoot header files ---
#include "AliESDCaloCluster.h"
#include "AliESDEvent.h"
#include "AliLog.h"
#include "AliPHOSDigit.h"
#include "AliPHOSHit.h"
#include "AliPHOSQualAssDataMaker.h"
#include "AliPHOSCpvRecPoint.h" 
#include "AliPHOSEmcRecPoint.h" 
#include "AliPHOSRecParticle.h" 
#include "AliPHOSTrackSegment.h" 

ClassImp(AliPHOSQualAssDataMaker)
           
//____________________________________________________________________________ 
  AliPHOSQualAssDataMaker::AliPHOSQualAssDataMaker() : 
  AliQualAssDataMaker(AliQualAss::GetDetName(AliQualAss::kPHOS), "PHOS Quality Assurance Data Maker"),
  fhHits(0x0), 
  fhHitsMul(0x0), 
  fhDigits(0x0),
  fhDigitsMul(0x0),
  fhSDigits(0x0),
  fhSDigitsMul(0x0),
  fhEmcRecPoints(0x0),
  fhEmcRecPointsMul(0x0),
  fhCpvRecPoints(0x0),
  fhCpvRecPointsMul(0x0),
  fhTrackSegments(0x0),
  fhTrackSegmentsMul(0x0),
  fhRecParticles(0x0),
  fhRecParticlesMul(0x0),
  fhESDs(0x0),
  fhESDsMul(0x0) 
{
  // ctor
  fDetectorDir = fOutput->GetDirectory(GetName()) ;  
  if (!fDetectorDir) 
    fDetectorDir = fOutput->mkdir(GetName()) ;  
}

//____________________________________________________________________________ 
AliPHOSQualAssDataMaker::AliPHOSQualAssDataMaker(const AliPHOSQualAssDataMaker& qadm) :
  AliQualAssDataMaker(), 
  fhHits(qadm.fhHits), 
  fhHitsMul(qadm.fhHitsMul), 
  fhDigits(qadm.fhDigits),
  fhDigitsMul(qadm.fhDigitsMul),
  fhSDigits(qadm.fhSDigits),
  fhSDigitsMul(qadm.fhSDigitsMul), 
  fhEmcRecPoints(qadm.fhEmcRecPoints),
  fhEmcRecPointsMul(qadm.fhEmcRecPointsMul), 
  fhCpvRecPoints(qadm.fhCpvRecPoints),
  fhCpvRecPointsMul(qadm.fhCpvRecPointsMul), 
  fhTrackSegments(qadm.fhTrackSegments),
  fhTrackSegmentsMul(qadm.fhTrackSegmentsMul), 
  fhRecParticles(qadm.fhRecParticles),
  fhRecParticlesMul(qadm.fhRecParticlesMul), 
  fhESDs(qadm.fhESDs), 
  fhESDsMul(qadm.fhESDsMul) 
{
  //copy ctor 
  SetName((const char*)qadm.GetName()) ; 
  SetTitle((const char*)qadm.GetTitle()); 
}

//__________________________________________________________________
AliPHOSQualAssDataMaker& AliPHOSQualAssDataMaker::operator = (const AliPHOSQualAssDataMaker& qadm )
{
  // Equal operator.
  this->~AliPHOSQualAssDataMaker();
  new(this) AliPHOSQualAssDataMaker(qadm);
  return *this;
}
 
//____________________________________________________________________________ 
void AliPHOSQualAssDataMaker::InitESDs()
{
  //create ESDs histograms in ESDs subdir
  fhESDs     = new TH1F("hPhosESDs",    "ESDs energy distribution in PHOS",       100, 0., 100.) ; 
  fhESDs->Sumw2() ; 
  fhESDsMul  = new TH1I("hPhosESDsMul", "ESDs multiplicity distribution in PHOS", 100, 0., 100) ; 
  fhESDsMul->Sumw2() ;
}

//____________________________________________________________________________ 
void AliPHOSQualAssDataMaker::InitHits()
{
  // create Hits histograms in Hits subdir
  fhHits     = new TH1F("hPhosHits",    "Hits energy distribution in PHOS",       100, 0., 100.) ; 
  fhHits->Sumw2() ;
  fhHitsMul  = new TH1I("hPhosHitsMul", "Hits multiplicity distribution in PHOS", 500, 0., 10000) ; 
  fhHitsMul->Sumw2() ;
}

//____________________________________________________________________________ 
void AliPHOSQualAssDataMaker::InitDigits()
{
  // create Digits histograms in Digits subdir
  fhDigits     = new TH1I("hPhosDigits",    "Digits amplitude distribution in PHOS",    500, 0, 5000) ; 
  fhDigits->Sumw2() ;
  fhDigitsMul  = new TH1I("hPhosDigitsMul", "Digits multiplicity distribution in PHOS", 500, 0, 1000) ; 
  fhDigitsMul->Sumw2() ;
}

//____________________________________________________________________________ 
void AliPHOSQualAssDataMaker::InitRecParticles()
{
  // create Reconstructed particles histograms in RecParticles subdir
  fhRecParticles     = new TH1F("hPhosRecParticles",    "RecParticles energy distribution in PHOS",       100, 0., 100.) ; 
  fhRecParticles->Sumw2() ;
  fhRecParticlesMul  = new TH1I("hPhosRecParticlesMul", "RecParticles multiplicity distribution in PHOS", 100, 0,  100) ; 
  fhRecParticlesMul->Sumw2() ;
}

//____________________________________________________________________________ 
void AliPHOSQualAssDataMaker::InitRecPoints()
{
  // create Reconstructed Points histograms in RecPoints subdir
  fhEmcRecPoints     = new TH1F("hEmcPhosRecPoints",    "EMCA RecPoints energy distribution in PHOS",       100, 0., 100.) ; 
  fhEmcRecPoints->Sumw2() ;
  fhEmcRecPointsMul  = new TH1I("hEmcPhosRecPointsMul", "EMCA RecPoints multiplicity distribution in PHOS", 100, 0,  100) ; 
  fhEmcRecPointsMul->Sumw2() ;

  fhCpvRecPoints     = new TH1F("hCpvPhosRecPoints",    "CPV RecPoints energy distribution in PHOS",       100, 0., 100.) ; 
  fhCpvRecPoints->Sumw2() ;
  fhCpvRecPointsMul  = new TH1I("hCpvPhosRecPointsMul", "CPV RecPoints multiplicity distribution in PHOS", 100, 0,  100) ; 
  fhCpvRecPointsMul->Sumw2() ;
}

//____________________________________________________________________________ 
void AliPHOSQualAssDataMaker::InitSDigits()
{
  // create SDigits histograms in SDigits subdir
  fhSDigits     = new TH1F("hPhosSDigits",    "SDigits energy distribution in PHOS",       100, 0., 100.) ; 
  fhSDigits->Sumw2() ;
  fhSDigitsMul  = new TH1I("hPhosSDigitsMul", "SDigits multiplicity distribution in PHOS", 500, 0,  10000) ; 
  fhSDigitsMul->Sumw2() ;
}

//____________________________________________________________________________ 
void AliPHOSQualAssDataMaker::InitTrackSegments()
{
  // create Track Segments histograms in TrackSegments subdir
  fhTrackSegments     = new TH1F("hPhosTrackSegments",    "TrackSegments EMC-CPV distance in PHOS",       500, 0., 5000.) ; 
  fhTrackSegments->Sumw2() ;
  fhTrackSegmentsMul  = new TH1I("hPhosTrackSegmentsMul", "TrackSegments multiplicity distribution in PHOS", 100, 0,  100) ; 
  fhTrackSegmentsMul->Sumw2() ;
}

//____________________________________________________________________________
void AliPHOSQualAssDataMaker::MakeESDs()
{
  // make QA data from ESDs
  
  AliESDEvent * esd = dynamic_cast<AliESDEvent*>(fData) ; 
  Int_t maxClu = esd->GetNumberOfPHOSClusters() ; 
  Int_t index = 0, count = 0 ; 
  for ( index = 0 ; index < maxClu; index++ ) {
    AliESDCaloCluster * clu = esd->GetCaloCluster(index) ;
    fhESDs->Fill(clu->E()) ;
    count++ ; 
  }
  fhESDsMul->Fill(count) ;
}

//____________________________________________________________________________
void AliPHOSQualAssDataMaker::MakeHits()
{
  //make QA data from Hits

  TClonesArray * hits = dynamic_cast<TClonesArray*>(fData) ; 
  fhHitsMul->Fill(hits->GetEntriesFast()) ; 
  TIter next(hits) ; 
  AliPHOSHit * hit ; 
  while ( (hit = dynamic_cast<AliPHOSHit *>(next())) ) {
    fhHits->Fill( hit->GetEnergy()) ;
  }
} 
 
//____________________________________________________________________________
void AliPHOSQualAssDataMaker::MakeDigits()
{
  // makes data from Digits

  TClonesArray * digits = dynamic_cast<TClonesArray*>(fData) ; 
  fhDigitsMul->Fill(digits->GetEntriesFast()) ; 
  TIter next(digits) ; 
  AliPHOSDigit * digit ; 
  while ( (digit = dynamic_cast<AliPHOSDigit *>(next())) ) {
    fhDigits->Fill( digit->GetEnergy()) ;
  }  
}

//____________________________________________________________________________
void AliPHOSQualAssDataMaker::MakeRecParticles()
{
  // makes data from RecParticles

  TClonesArray * recparticles = dynamic_cast<TClonesArray*>(fData) ; 
  fhRecParticlesMul->Fill(recparticles->GetEntriesFast()) ; 
  TIter next(recparticles) ; 
  AliPHOSRecParticle * recparticle ; 
  while ( (recparticle = dynamic_cast<AliPHOSRecParticle *>(next())) ) {
    fhRecParticles->Fill( recparticle->Energy()) ;
  }
}

//____________________________________________________________________________
void AliPHOSQualAssDataMaker::MakeRecPoints()
{
  // makes data from RecPoints
  TObjArray * recpoints = dynamic_cast<TObjArray*>(fData) ;  
  TIter next(recpoints) ; 
  
  if ( strcmp(fData->GetName(), "EMCRECPOINTS") == 0 ) {
    fhEmcRecPointsMul->Fill(recpoints->GetEntriesFast()) ; 
    AliPHOSEmcRecPoint * rp ; 
    while ( (rp = dynamic_cast<AliPHOSEmcRecPoint *>(next())) ) {
      fhEmcRecPoints->Fill( rp->GetEnergy()) ;
    }
  } 
  else if  ( strcmp(fData->GetName(), "CPVRECPOINTS") == 0 ) {
    fhCpvRecPointsMul->Fill(recpoints->GetEntriesFast()) ; 
    AliPHOSCpvRecPoint * rp ; 
    while ( (rp = dynamic_cast<AliPHOSCpvRecPoint *>(next())) ) {
      fhCpvRecPoints->Fill( rp->GetEnergy()) ;
    }
  }  
}

//____________________________________________________________________________
void AliPHOSQualAssDataMaker::MakeSDigits()
{
  // makes data from SDigits

  TClonesArray * sdigits = dynamic_cast<TClonesArray*>(fData) ; 
  fhSDigitsMul->Fill(sdigits->GetEntriesFast()) ; 
  TIter next(sdigits) ; 
  AliPHOSDigit * sdigit ; 
  while ( (sdigit = dynamic_cast<AliPHOSDigit *>(next())) ) {
    fhSDigits->Fill( sdigit->GetEnergy()) ;
  } 
}

//____________________________________________________________________________
void AliPHOSQualAssDataMaker::MakeTrackSegments()
{
  // makes data from TrackSegments

  TClonesArray * tracksegments = dynamic_cast<TClonesArray*>(fData) ;

  fhTrackSegmentsMul->Fill(tracksegments->GetEntriesFast()) ; 
  TIter next(tracksegments) ; 
  AliPHOSTrackSegment * ts ; 
  while ( (ts = dynamic_cast<AliPHOSTrackSegment *>(next())) ) {
    fhTrackSegments->Fill( ts->GetCpvDistance()) ;
  } 
}
Commit	Line	Data
ddd1a39c	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
	31	// --- Standard library ---
	32
	33	// --- AliRoot header files ---
	34	#include "AliESDCaloCluster.h"
	35	#include "AliESDEvent.h"
	36	#include "AliLog.h"
	37	#include "AliPHOSDigit.h"
	38	#include "AliPHOSHit.h"
	39	#include "AliPHOSQualAssDataMaker.h"
	40	#include "AliPHOSCpvRecPoint.h"
	41	#include "AliPHOSEmcRecPoint.h"
	42	#include "AliPHOSRecParticle.h"
	43	#include "AliPHOSTrackSegment.h"
	44
	45	ClassImp(AliPHOSQualAssDataMaker)
	46
	47	//____________________________________________________________________________
	48	AliPHOSQualAssDataMaker::AliPHOSQualAssDataMaker() :
a5fa6165	49	AliQualAssDataMaker(AliQualAss::GetDetName(AliQualAss::kPHOS), "PHOS Quality Assurance Data Maker"),
ddd1a39c	50	fhHits(0x0),
	51	fhHitsMul(0x0),
	52	fhDigits(0x0),
	53	fhDigitsMul(0x0),
	54	fhSDigits(0x0),
	55	fhSDigitsMul(0x0),
	56	fhEmcRecPoints(0x0),
	57	fhEmcRecPointsMul(0x0),
	58	fhCpvRecPoints(0x0),
	59	fhCpvRecPointsMul(0x0),
	60	fhTrackSegments(0x0),
	61	fhTrackSegmentsMul(0x0),
	62	fhRecParticles(0x0),
	63	fhRecParticlesMul(0x0),
	64	fhESDs(0x0),
	65	fhESDsMul(0x0)
	66	{
	67	// ctor
	68	fDetectorDir = fOutput->GetDirectory(GetName()) ;
	69	if (!fDetectorDir)
	70	fDetectorDir = fOutput->mkdir(GetName()) ;
	71	}
	72
	73	//____________________________________________________________________________
	74	AliPHOSQualAssDataMaker::AliPHOSQualAssDataMaker(const AliPHOSQualAssDataMaker& qadm) :
	75	AliQualAssDataMaker(),
	76	fhHits(qadm.fhHits),
	77	fhHitsMul(qadm.fhHitsMul),
	78	fhDigits(qadm.fhDigits),
	79	fhDigitsMul(qadm.fhDigitsMul),
	80	fhSDigits(qadm.fhSDigits),
	81	fhSDigitsMul(qadm.fhSDigitsMul),
	82	fhEmcRecPoints(qadm.fhEmcRecPoints),
	83	fhEmcRecPointsMul(qadm.fhEmcRecPointsMul),
	84	fhCpvRecPoints(qadm.fhCpvRecPoints),
	85	fhCpvRecPointsMul(qadm.fhCpvRecPointsMul),
	86	fhTrackSegments(qadm.fhTrackSegments),
	87	fhTrackSegmentsMul(qadm.fhTrackSegmentsMul),
	88	fhRecParticles(qadm.fhRecParticles),
	89	fhRecParticlesMul(qadm.fhRecParticlesMul),
	90	fhESDs(qadm.fhESDs),
	91	fhESDsMul(qadm.fhESDsMul)
	92	{
	93	//copy ctor
	94	SetName((const char*)qadm.GetName()) ;
	95	SetTitle((const char*)qadm.GetTitle());
	96	}
	97
	98	//__________________________________________________________________
	99	AliPHOSQualAssDataMaker& AliPHOSQualAssDataMaker::operator = (const AliPHOSQualAssDataMaker& qadm )
	100	{
	101	// Equal operator.
	102	this->~AliPHOSQualAssDataMaker();
	103	new(this) AliPHOSQualAssDataMaker(qadm);
	104	return *this;
	105	}
	106
	107	//____________________________________________________________________________
	108	void AliPHOSQualAssDataMaker::InitESDs()
	109	{
	110	//create ESDs histograms in ESDs subdir
	111	fhESDs = new TH1F("hPhosESDs", "ESDs energy distribution in PHOS", 100, 0., 100.) ;
a5fa6165	112	fhESDs->Sumw2() ;
ddd1a39c	113	fhESDsMul = new TH1I("hPhosESDsMul", "ESDs multiplicity distribution in PHOS", 100, 0., 100) ;
a5fa6165	114	fhESDsMul->Sumw2() ;
ddd1a39c	115	}
	116
	117	//____________________________________________________________________________
	118	void AliPHOSQualAssDataMaker::InitHits()
	119	{
	120	// create Hits histograms in Hits subdir
	121	fhHits = new TH1F("hPhosHits", "Hits energy distribution in PHOS", 100, 0., 100.) ;
a5fa6165	122	fhHits->Sumw2() ;
ddd1a39c	123	fhHitsMul = new TH1I("hPhosHitsMul", "Hits multiplicity distribution in PHOS", 500, 0., 10000) ;
a5fa6165	124	fhHitsMul->Sumw2() ;
ddd1a39c	125	}
	126
	127	//____________________________________________________________________________
	128	void AliPHOSQualAssDataMaker::InitDigits()
	129	{
	130	// create Digits histograms in Digits subdir
	131	fhDigits = new TH1I("hPhosDigits", "Digits amplitude distribution in PHOS", 500, 0, 5000) ;
a5fa6165	132	fhDigits->Sumw2() ;
ddd1a39c	133	fhDigitsMul = new TH1I("hPhosDigitsMul", "Digits multiplicity distribution in PHOS", 500, 0, 1000) ;
a5fa6165	134	fhDigitsMul->Sumw2() ;
ddd1a39c	135	}
	136
	137	//____________________________________________________________________________
	138	void AliPHOSQualAssDataMaker::InitRecParticles()
	139	{
	140	// create Reconstructed particles histograms in RecParticles subdir
	141	fhRecParticles = new TH1F("hPhosRecParticles", "RecParticles energy distribution in PHOS", 100, 0., 100.) ;
a5fa6165	142	fhRecParticles->Sumw2() ;
ddd1a39c	143	fhRecParticlesMul = new TH1I("hPhosRecParticlesMul", "RecParticles multiplicity distribution in PHOS", 100, 0, 100) ;
a5fa6165	144	fhRecParticlesMul->Sumw2() ;
ddd1a39c	145	}
	146
	147	//____________________________________________________________________________
	148	void AliPHOSQualAssDataMaker::InitRecPoints()
	149	{
	150	// create Reconstructed Points histograms in RecPoints subdir
	151	fhEmcRecPoints = new TH1F("hEmcPhosRecPoints", "EMCA RecPoints energy distribution in PHOS", 100, 0., 100.) ;
a5fa6165	152	fhEmcRecPoints->Sumw2() ;
ddd1a39c	153	fhEmcRecPointsMul = new TH1I("hEmcPhosRecPointsMul", "EMCA RecPoints multiplicity distribution in PHOS", 100, 0, 100) ;
a5fa6165	154	fhEmcRecPointsMul->Sumw2() ;
ddd1a39c	155
ddd1a39c	156	fhCpvRecPoints = new TH1F("hCpvPhosRecPoints", "CPV RecPoints energy distribution in PHOS", 100, 0., 100.) ;
a5fa6165	157	fhCpvRecPoints->Sumw2() ;
ddd1a39c	158	fhCpvRecPointsMul = new TH1I("hCpvPhosRecPointsMul", "CPV RecPoints multiplicity distribution in PHOS", 100, 0, 100) ;
a5fa6165	159	fhCpvRecPointsMul->Sumw2() ;
ddd1a39c	160	}
	161
	162	//____________________________________________________________________________
	163	void AliPHOSQualAssDataMaker::InitSDigits()
	164	{
	165	// create SDigits histograms in SDigits subdir
	166	fhSDigits = new TH1F("hPhosSDigits", "SDigits energy distribution in PHOS", 100, 0., 100.) ;
a5fa6165	167	fhSDigits->Sumw2() ;
ddd1a39c	168	fhSDigitsMul = new TH1I("hPhosSDigitsMul", "SDigits multiplicity distribution in PHOS", 500, 0, 10000) ;
a5fa6165	169	fhSDigitsMul->Sumw2() ;
ddd1a39c	170	}
	171
	172	//____________________________________________________________________________
	173	void AliPHOSQualAssDataMaker::InitTrackSegments()
	174	{
	175	// create Track Segments histograms in TrackSegments subdir
	176	fhTrackSegments = new TH1F("hPhosTrackSegments", "TrackSegments EMC-CPV distance in PHOS", 500, 0., 5000.) ;
a5fa6165	177	fhTrackSegments->Sumw2() ;
ddd1a39c	178	fhTrackSegmentsMul = new TH1I("hPhosTrackSegmentsMul", "TrackSegments multiplicity distribution in PHOS", 100, 0, 100) ;
a5fa6165	179	fhTrackSegmentsMul->Sumw2() ;
ddd1a39c	180	}
	181
	182	//____________________________________________________________________________
	183	void AliPHOSQualAssDataMaker::MakeESDs()
	184	{
	185	// make QA data from ESDs
	186
	187	AliESDEvent * esd = dynamic_cast<AliESDEvent*>(fData) ;
	188	Int_t maxClu = esd->GetNumberOfPHOSClusters() ;
	189	Int_t index = 0, count = 0 ;
	190	for ( index = 0 ; index < maxClu; index++ ) {
	191	AliESDCaloCluster * clu = esd->GetCaloCluster(index) ;
	192	fhESDs->Fill(clu->E()) ;
	193	count++ ;
	194	}
	195	fhESDsMul->Fill(count) ;
	196	}
	197
	198	//____________________________________________________________________________
	199	void AliPHOSQualAssDataMaker::MakeHits()
	200	{
	201	//make QA data from Hits
	202
	203	TClonesArray * hits = dynamic_cast<TClonesArray*>(fData) ;
	204	fhHitsMul->Fill(hits->GetEntriesFast()) ;
	205	TIter next(hits) ;
	206	AliPHOSHit * hit ;
	207	while ( (hit = dynamic_cast<AliPHOSHit *>(next())) ) {
	208	fhHits->Fill( hit->GetEnergy()) ;
	209	}
	210	}
	211
	212	//____________________________________________________________________________
	213	void AliPHOSQualAssDataMaker::MakeDigits()
	214	{
	215	// makes data from Digits
	216
	217	TClonesArray * digits = dynamic_cast<TClonesArray*>(fData) ;
	218	fhDigitsMul->Fill(digits->GetEntriesFast()) ;
	219	TIter next(digits) ;
	220	AliPHOSDigit * digit ;
	221	while ( (digit = dynamic_cast<AliPHOSDigit *>(next())) ) {
	222	fhDigits->Fill( digit->GetEnergy()) ;
	223	}
	224	}
	225
	226	//____________________________________________________________________________
	227	void AliPHOSQualAssDataMaker::MakeRecParticles()
	228	{
	229	// makes data from RecParticles
	230
	231	TClonesArray * recparticles = dynamic_cast<TClonesArray*>(fData) ;
	232	fhRecParticlesMul->Fill(recparticles->GetEntriesFast()) ;
	233	TIter next(recparticles) ;
	234	AliPHOSRecParticle * recparticle ;
	235	while ( (recparticle = dynamic_cast<AliPHOSRecParticle *>(next())) ) {
	236	fhRecParticles->Fill( recparticle->Energy()) ;
	237	}
	238	}
	239
	240	//____________________________________________________________________________
	241	void AliPHOSQualAssDataMaker::MakeRecPoints()
	242	{
	243	// makes data from RecPoints
244	TObjArray * recpoints = dynamic_cast<TObjArray*>(fData) ;
245	TIter next(recpoints) ;
246
247	if ( strcmp(fData->GetName(), "EMCRECPOINTS") == 0 ) {
248	fhEmcRecPointsMul->Fill(recpoints->GetEntriesFast()) ;
249	AliPHOSEmcRecPoint * rp ;
250	while ( (rp = dynamic_cast<AliPHOSEmcRecPoint *>(next())) ) {
251	fhEmcRecPoints->Fill( rp->GetEnergy()) ;
252	}
253	}
254	else if ( strcmp(fData->GetName(), "CPVRECPOINTS") == 0 ) {
255	fhCpvRecPointsMul->Fill(recpoints->GetEntriesFast()) ;
256	AliPHOSCpvRecPoint * rp ;
257	while ( (rp = dynamic_cast<AliPHOSCpvRecPoint *>(next())) ) {
258	fhCpvRecPoints->Fill( rp->GetEnergy()) ;
259	}
260	}
261	}
262
263	//____________________________________________________________________________
264	void AliPHOSQualAssDataMaker::MakeSDigits()
265	{
266	// makes data from SDigits
267
268	TClonesArray * sdigits = dynamic_cast<TClonesArray*>(fData) ;
269	fhSDigitsMul->Fill(sdigits->GetEntriesFast()) ;
270	TIter next(sdigits) ;
271	AliPHOSDigit * sdigit ;
272	while ( (sdigit = dynamic_cast<AliPHOSDigit *>(next())) ) {
273	fhSDigits->Fill( sdigit->GetEnergy()) ;
274	}
275	}
276
277	//____________________________________________________________________________
278	void AliPHOSQualAssDataMaker::MakeTrackSegments()
279	{
280	// makes data from TrackSegments
281
282	TClonesArray * tracksegments = dynamic_cast<TClonesArray*>(fData) ;
283
284	fhTrackSegmentsMul->Fill(tracksegments->GetEntriesFast()) ;
285	TIter next(tracksegments) ;
286	AliPHOSTrackSegment * ts ;
287	while ( (ts = dynamic_cast<AliPHOSTrackSegment *>(next())) ) {
288	fhTrackSegments->Fill( ts->GetCpvDistance()) ;
289	}
290	}