[u/mrichter/AliRoot.git] / FMD / AliFMDQADataMakerRec.cxx

/**************************************************************************
 * Copyright(c) 2004, 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.                  *
 **************************************************************************/
// --- ROOT system ---
#include <iostream>
#include <TClonesArray.h>
#include <TFile.h> 
#include <TH1F.h> 
#include <TH1I.h> 

// --- AliRoot header files ---
#include "AliESDEvent.h"
#include "AliLog.h"
#include "AliFMDQADataMakerRec.h"
#include "AliFMDDigit.h"
#include "AliFMDRecPoint.h"
#include "AliQAChecker.h"
#include "AliESDFMD.h"
#include "AliFMDParameters.h"
#include "AliFMDRawReader.h"
#include "AliRawReader.h"

//_____________________________________________________________________
// This is the class that collects the QA data for the FMD during
// reconstruction.  
//
// The following data types are picked up:
// - rec points
// - esd data
// - raws
// Author : Hans Hjersing Dalsgaard, hans.dalsgaard@cern.ch
//_____________________________________________________________________

ClassImp(AliFMDQADataMakerRec)
#if 0
; // For Emacs - do not delete!
#endif
           
//_____________________________________________________________________
AliFMDQADataMakerRec::AliFMDQADataMakerRec() : 
  AliQADataMakerRec(AliQA::GetDetName(AliQA::kFMD), 
		    "FMD Quality Assurance Data Maker"),
  fDigitsArray("AliFMDDigit", 0),
  fRecPointsArray("AliFMDRecPoint", 1000)
{
  // ctor
 
}

//_____________________________________________________________________
AliFMDQADataMakerRec::AliFMDQADataMakerRec(const AliFMDQADataMakerRec& qadm) 
  : AliQADataMakerRec(AliQA::GetDetName(AliQA::kFMD), 
		      "FMD Quality Assurance Data Maker"),
    fDigitsArray(qadm.fDigitsArray),
    fRecPointsArray(qadm.fRecPointsArray)
{
  // copy ctor 
  // Parameters: 
  //    qadm    Object to copy from
  
}
//_____________________________________________________________________
AliFMDQADataMakerRec& AliFMDQADataMakerRec::operator = (const AliFMDQADataMakerRec& qadm ) 
{
  fDigitsArray = qadm.fDigitsArray;
  fRecPointsArray = qadm.fRecPointsArray;
  
  return *this;
}
//_____________________________________________________________________
AliFMDQADataMakerRec::~AliFMDQADataMakerRec()
{
 
}


//_____________________________________________________________________ 

void 
AliFMDQADataMakerRec::EndOfDetectorCycle(AliQA::TASKINDEX_t task, 
					 TObjArray * list)
{
  // Detector specific actions at end of cycle
  // do the QA checking
  AliLog::Message(5,"FMD: end of detector cycle",
		  "AliFMDQADataMakerRec","AliFMDQADataMakerRec",
		  "AliFMDQADataMakerRec::EndOfDetectorCycle",
		  "AliFMDQADataMakerRec.cxx",95);
  AliQAChecker::Instance()->Run(AliQA::kFMD, task, list);
}

//_____________________________________________________________________ 
void AliFMDQADataMakerRec::InitESDs()
{
  // create Digits histograms in Digits subdir
  TH1F* hEnergyOfESD = new TH1F("hEnergyOfESD","Energy distribution",100,0,3);
  hEnergyOfESD->SetXTitle("Edep/Emip");
  hEnergyOfESD->SetYTitle("Counts");
  Add2ESDsList(hEnergyOfESD, 0);
    
}

//_____________________________________________________________________ 
/*void AliFMDQADataMakerRec::InitDigits()
{
  // create Digits histograms in Digits subdir
  TH1I* hADCCounts      = new TH1I("hADCCounts","Dist of ADC counts",
				   1024,0,1024);
  hADCCounts->SetXTitle("ADC counts");
  hADCCounts->SetYTitle("");
  Add2DigitsList(hADCCounts, 0);
  
}
*/
//_____________________________________________________________________ 
void AliFMDQADataMakerRec::InitRecPoints()
{
  TH1F* hEnergyOfRecpoints = new TH1F("hEnergyOfRecpoints",
				      "Energy Distribution",100,0,3);
  hEnergyOfRecpoints->SetXTitle("Edep/Emip");
  hEnergyOfRecpoints->SetYTitle("");
  Add2RecPointsList(hEnergyOfRecpoints,0);
}

//_____________________________________________________________________ 
void AliFMDQADataMakerRec::InitRaws()
{
  
  TH1I* hADCCounts;
  for(Int_t det = 1; det<=3; det++) {
    Int_t firstring = (det==1 ? 1 : 0);
    for(Int_t iring = firstring;iring<=1;iring++) {
      for(Int_t b = 0; b<=1;b++) {
	
	//Hexadecimal board numbers 0x0, 0x1, 0x10, 0x10;
	UInt_t board = (iring == 1 ? 0 : 1);
	board = board + b*16;
	Char_t ring = (iring == 1 ? 'I' : 'O');
	
	hADCCounts      = new TH1I(Form("hADCCounts_FMD%d%c_board%d",
					det, ring, board), "ADC counts",
				   1024,0,1023);
	hADCCounts->SetXTitle("ADC counts");
	hADCCounts->SetYTitle("");
	Int_t index = GetHalfringIndex(det, ring, board/16);
	Add2RawsList(hADCCounts, index);
	
      }
    }
  }
}

//_____________________________________________________________________
void AliFMDQADataMakerRec::MakeESDs(AliESDEvent * esd)
{
  if(!esd) {
    AliError("FMD ESD object not found!!") ; 
    return;
  }
  AliESDFMD* fmd = esd->GetFMDData();
  if (!fmd) return;
  
  for(UShort_t det=1;det<=3;det++) {
    for (UShort_t ir = 0; ir < 2; ir++) {
      Char_t   ring = (ir == 0 ? 'I' : 'O');
      UShort_t nsec = (ir == 0 ? 20  : 40);
      UShort_t nstr = (ir == 0 ? 512 : 256);
      for(UShort_t sec =0; sec < nsec;  sec++)  {
	for(UShort_t strip = 0; strip < nstr; strip++) {
	  Float_t mult = fmd->Multiplicity(det,ring,sec,strip);
	  if(mult == AliESDFMD::kInvalidMult) continue;
	  
	  GetESDsData(0)->Fill(mult);
	}
      }
    }
  }
}

/*
//_____________________________________________________________________
void AliFMDQADataMakerRec::MakeDigits(TClonesArray * digits)
{
  // makes data from Digits  
  if(!digits)  {
    AliError("FMD Digit object not found!!") ;
    return;
  }
  for(Int_t i=0;i<digits->GetEntriesFast();i++) {
    //Raw ADC counts
    AliFMDDigit* digit = static_cast<AliFMDDigit*>(digits->At(i));
    GetDigitsData(0)->Fill(digit->Counts());
  }
}

//_____________________________________________________________________
void AliFMDQADataMakerRec::MakeDigits(TTree * digitTree)
{
  
  fDigitsArray.Clear();
  TBranch*      branch = digitTree->GetBranch("FMD");
  if (!branch) {
    AliWarning("FMD branch in Digit Tree not found") ; 
    return;
  } 
  TClonesArray* digitsAddress = &fDigitsArray;
  branch->SetAddress(&digitsAddress);
  branch->GetEntry(0); 
  MakeDigits(digitsAddress);
}
*/
//_____________________________________________________________________
void AliFMDQADataMakerRec::MakeRaws(AliRawReader* rawReader)
{
 
  AliFMDRawReader fmdReader(rawReader,0);
  TClonesArray* digitsAddress = &fDigitsArray;
  
  rawReader->Reset();
		
	digitsAddress->Clear();
	fmdReader.ReadAdcs(digitsAddress);
	for(Int_t i=0;i<digitsAddress->GetEntriesFast();i++) {
	  //Raw ADC counts
	  AliFMDDigit* digit = static_cast<AliFMDDigit*>(digitsAddress->At(i));
	  UShort_t det = digit->Detector();
	  Char_t ring  = digit->Ring();
	  UShort_t sec = digit->Sector();
	  UShort_t strip = digit->Strip();
	  UInt_t ddl, board, chip, channel;
	  AliFMDParameters* pars = AliFMDParameters::Instance();
	  pars->Detector2Hardware(det,ring,sec,strip,ddl,board,chip,channel);
	  Int_t index = GetHalfringIndex(det, ring, board/16);
	  
	  GetRawsData(index)->Fill(digit->Counts());
	  
	}
}

//_____________________________________________________________________
void AliFMDQADataMakerRec::MakeRecPoints(TTree* clustersTree)
{
  // makes data from RecPoints
  AliFMDParameters* pars = AliFMDParameters::Instance();
  fRecPointsArray.Clear();
  TBranch *fmdbranch = clustersTree->GetBranch("FMD");
  if (!fmdbranch) { 
    AliError("can't get the branch with the FMD recpoints !");
    return;
  }
  
  TClonesArray* RecPointsAddress = &fRecPointsArray;
  
  fmdbranch->SetAddress(&RecPointsAddress);
  fmdbranch->GetEntry(0);
  TIter next(RecPointsAddress) ; 
  AliFMDRecPoint * rp ; 
  while ((rp = static_cast<AliFMDRecPoint*>(next()))) {
    
    GetRecPointsData(0)->Fill(rp->Edep()/pars->GetEdepMip()) ;
  
  }

}

//_____________________________________________________________________ 
void AliFMDQADataMakerRec::StartOfDetectorCycle()
{
  // What 
  // to 
  // do?
}
//_____________________________________________________________________ 
Int_t AliFMDQADataMakerRec::GetHalfringIndex(UShort_t det, 
					     Char_t ring, 
					     UShort_t board) {
  
  UShort_t iring  =  (ring == 'I' ? 1 : 0);
  
  Int_t index = (((det-1) << 2) | (iring << 1) | (board << 0));
  
  return index-2;
  
}

//_____________________________________________________________________ 


//
// EOF
//
Commit	Line	Data
c9dd1c4d	1	/**************************************************************************
	2	* Copyright(c) 2004, 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	// --- ROOT system ---
	16	#include <iostream>
	17	#include <TClonesArray.h>
	18	#include <TFile.h>
	19	#include <TH1F.h>
	20	#include <TH1I.h>
	21
	22	// --- AliRoot header files ---
	23	#include "AliESDEvent.h"
	24	#include "AliLog.h"
	25	#include "AliFMDQADataMakerRec.h"
	26	#include "AliFMDDigit.h"
	27	#include "AliFMDRecPoint.h"
	28	#include "AliQAChecker.h"
	29	#include "AliESDFMD.h"
	30	#include "AliFMDParameters.h"
86c6cec8	31	#include "AliFMDRawReader.h"
3ceaa9ad	32	#include "AliRawReader.h"
c9dd1c4d	33
	34	//_____________________________________________________________________
	35	// This is the class that collects the QA data for the FMD during
	36	// reconstruction.
	37	//
	38	// The following data types are picked up:
c9dd1c4d	39	// - rec points
c9dd1c4d	40	// - esd data
c9dd1c4d	41	// - raws
	42	// Author : Hans Hjersing Dalsgaard, hans.dalsgaard@cern.ch
	43	//_____________________________________________________________________
	44
	45	ClassImp(AliFMDQADataMakerRec)
	46	#if 0
	47	; // For Emacs - do not delete!
	48	#endif
	49
	50	//_____________________________________________________________________
	51	AliFMDQADataMakerRec::AliFMDQADataMakerRec() :
	52	AliQADataMakerRec(AliQA::GetDetName(AliQA::kFMD),
a7e41e8d	53	"FMD Quality Assurance Data Maker"),
198942c7	54	fDigitsArray("AliFMDDigit", 0),
198942c7	55	fRecPointsArray("AliFMDRecPoint", 1000)
c9dd1c4d	56	{
c9dd1c4d	57	// ctor
56236ce9	58
c9dd1c4d	59	}
	60
	61	//_____________________________________________________________________
a7e41e8d	62	AliFMDQADataMakerRec::AliFMDQADataMakerRec(const AliFMDQADataMakerRec& qadm)
	63	: AliQADataMakerRec(AliQA::GetDetName(AliQA::kFMD),
	64	"FMD Quality Assurance Data Maker"),
	65	fDigitsArray(qadm.fDigitsArray),
	66	fRecPointsArray(qadm.fRecPointsArray)
c9dd1c4d	67	{
ffa78f64	68	// copy ctor
c9dd1c4d	69	// Parameters:
	70	// qadm Object to copy from
	71
	72	}
9bd2ccc2	73	//_____________________________________________________________________
a7e41e8d	74	AliFMDQADataMakerRec& AliFMDQADataMakerRec::operator = (const AliFMDQADataMakerRec& qadm )
	75	{
	76	fDigitsArray = qadm.fDigitsArray;
	77	fRecPointsArray = qadm.fRecPointsArray;
	78
	79	return *this;
	80	}
	81	//_____________________________________________________________________
9bd2ccc2	82	AliFMDQADataMakerRec::~AliFMDQADataMakerRec()
9bd2ccc2	83	{
56236ce9	84
9bd2ccc2	85	}
c9dd1c4d	86
	87
	88	//_____________________________________________________________________
	89
	90	void
92a357bf	91	AliFMDQADataMakerRec::EndOfDetectorCycle(AliQA::TASKINDEX_t task,
c9dd1c4d	92	TObjArray * list)
	93	{
	94	// Detector specific actions at end of cycle
	95	// do the QA checking
ffa78f64	96	AliLog::Message(5,"FMD: end of detector cycle",
	97	"AliFMDQADataMakerRec","AliFMDQADataMakerRec",
	98	"AliFMDQADataMakerRec::EndOfDetectorCycle",
	99	"AliFMDQADataMakerRec.cxx",95);
c9dd1c4d	100	AliQAChecker::Instance()->Run(AliQA::kFMD, task, list);
	101	}
	102
	103	//_____________________________________________________________________
	104	void AliFMDQADataMakerRec::InitESDs()
	105	{
	106	// create Digits histograms in Digits subdir
	107	TH1F* hEnergyOfESD = new TH1F("hEnergyOfESD","Energy distribution",100,0,3);
	108	hEnergyOfESD->SetXTitle("Edep/Emip");
	109	hEnergyOfESD->SetYTitle("Counts");
	110	Add2ESDsList(hEnergyOfESD, 0);
	111
	112	}
	113
	114	//_____________________________________________________________________
86c6cec8	115	/*void AliFMDQADataMakerRec::InitDigits()
c9dd1c4d	116	{
	117	// create Digits histograms in Digits subdir
	118	TH1I* hADCCounts = new TH1I("hADCCounts","Dist of ADC counts",
	119	1024,0,1024);
c9dd1c4d	120	hADCCounts->SetXTitle("ADC counts");
c9dd1c4d	121	hADCCounts->SetYTitle("");
c9dd1c4d	122	Add2DigitsList(hADCCounts, 0);
c9dd1c4d	123
c9dd1c4d	124	}
86c6cec8	125	*/
c9dd1c4d	126	//_____________________________________________________________________
	127	void AliFMDQADataMakerRec::InitRecPoints()
	128	{
	129	TH1F* hEnergyOfRecpoints = new TH1F("hEnergyOfRecpoints",
	130	"Energy Distribution",100,0,3);
	131	hEnergyOfRecpoints->SetXTitle("Edep/Emip");
	132	hEnergyOfRecpoints->SetYTitle("");
	133	Add2RecPointsList(hEnergyOfRecpoints,0);
	134	}
	135
	136	//_____________________________________________________________________
	137	void AliFMDQADataMakerRec::InitRaws()
	138	{
d01e669e	139
	140	TH1I* hADCCounts;
	141	for(Int_t det = 1; det<=3; det++) {
	142	Int_t firstring = (det==1 ? 1 : 0);
	143	for(Int_t iring = firstring;iring<=1;iring++) {
	144	for(Int_t b = 0; b<=1;b++) {
	145
	146	//Hexadecimal board numbers 0x0, 0x1, 0x10, 0x10;
	147	UInt_t board = (iring == 1 ? 0 : 1);
	148	board = board + b*16;
	149	Char_t ring = (iring == 1 ? 'I' : 'O');
	150
	151	hADCCounts = new TH1I(Form("hADCCounts_FMD%d%c_board%d",
	152	det, ring, board), "ADC counts",
198942c7	153	1024,0,1023);
d01e669e	154	hADCCounts->SetXTitle("ADC counts");
	155	hADCCounts->SetYTitle("");
	156	Int_t index = GetHalfringIndex(det, ring, board/16);
	157	Add2RawsList(hADCCounts, index);
	158
	159	}
	160	}
	161	}
c9dd1c4d	162	}
	163
	164	//_____________________________________________________________________
	165	void AliFMDQADataMakerRec::MakeESDs(AliESDEvent * esd)
	166	{
	167	if(!esd) {
	168	AliError("FMD ESD object not found!!") ;
	169	return;
	170	}
	171	AliESDFMD* fmd = esd->GetFMDData();
	172	if (!fmd) return;
	173
	174	for(UShort_t det=1;det<=3;det++) {
	175	for (UShort_t ir = 0; ir < 2; ir++) {
	176	Char_t ring = (ir == 0 ? 'I' : 'O');
	177	UShort_t nsec = (ir == 0 ? 20 : 40);
	178	UShort_t nstr = (ir == 0 ? 512 : 256);
	179	for(UShort_t sec =0; sec < nsec; sec++) {
	180	for(UShort_t strip = 0; strip < nstr; strip++) {
	181	Float_t mult = fmd->Multiplicity(det,ring,sec,strip);
	182	if(mult == AliESDFMD::kInvalidMult) continue;
	183
	184	GetESDsData(0)->Fill(mult);
	185	}
	186	}
	187	}
	188	}
	189	}
	190
86c6cec8	191	/*
c9dd1c4d	192	//_____________________________________________________________________
	193	void AliFMDQADataMakerRec::MakeDigits(TClonesArray * digits)
	194	{
	195	// makes data from Digits
	196	if(!digits) {
ffa78f64	197	AliError("FMD Digit object not found!!") ;
ffa78f64	198	return;
c9dd1c4d	199	}
56236ce9	200	for(Int_t i=0;i<digits->GetEntriesFast();i++) {
c9dd1c4d	201	//Raw ADC counts
56236ce9	202	AliFMDDigit* digit = static_cast<AliFMDDigit*>(digits->At(i));
c9dd1c4d	203	GetDigitsData(0)->Fill(digit->Counts());
	204	}
	205	}
	206
	207	//_____________________________________________________________________
	208	void AliFMDQADataMakerRec::MakeDigits(TTree * digitTree)
	209	{
	210
56236ce9	211	fDigitsArray.Clear();
ffa78f64	212	TBranch* branch = digitTree->GetBranch("FMD");
ffa78f64	213	if (!branch) {
c9dd1c4d	214	AliWarning("FMD branch in Digit Tree not found") ;
	215	return;
	216	}
56236ce9	217	TClonesArray* digitsAddress = &fDigitsArray;
56236ce9	218	branch->SetAddress(&digitsAddress);
ffa78f64	219	branch->GetEntry(0);
56236ce9	220	MakeDigits(digitsAddress);
c9dd1c4d	221	}
86c6cec8	222	*/
c9dd1c4d	223	//_____________________________________________________________________
86c6cec8	224	void AliFMDQADataMakerRec::MakeRaws(AliRawReader* rawReader)
c9dd1c4d	225	{
198942c7	226
86c6cec8	227	AliFMDRawReader fmdReader(rawReader,0);
86c6cec8	228	TClonesArray* digitsAddress = &fDigitsArray;
3ceaa9ad	229
3ceaa9ad	230	rawReader->Reset();
78328afd	231
	232	digitsAddress->Clear();
	233	fmdReader.ReadAdcs(digitsAddress);
	234	for(Int_t i=0;i<digitsAddress->GetEntriesFast();i++) {
d01e669e	235	//Raw ADC counts
	236	AliFMDDigit* digit = static_cast<AliFMDDigit*>(digitsAddress->At(i));
	237	UShort_t det = digit->Detector();
	238	Char_t ring = digit->Ring();
	239	UShort_t sec = digit->Sector();
	240	UShort_t strip = digit->Strip();
	241	UInt_t ddl, board, chip, channel;
	242	AliFMDParameters* pars = AliFMDParameters::Instance();
	243	pars->Detector2Hardware(det,ring,sec,strip,ddl,board,chip,channel);
	244	Int_t index = GetHalfringIndex(det, ring, board/16);
	245
	246	GetRawsData(index)->Fill(digit->Counts());
	247
78328afd	248	}
c9dd1c4d	249	}
	250
	251	//_____________________________________________________________________
	252	void AliFMDQADataMakerRec::MakeRecPoints(TTree* clustersTree)
	253	{
	254	// makes data from RecPoints
198942c7	255	AliFMDParameters* pars = AliFMDParameters::Instance();
56236ce9	256	fRecPointsArray.Clear();
c9dd1c4d	257	TBranch *fmdbranch = clustersTree->GetBranch("FMD");
	258	if (!fmdbranch) {
	259	AliError("can't get the branch with the FMD recpoints !");
	260	return;
	261	}
	262
56236ce9	263	TClonesArray* RecPointsAddress = &fRecPointsArray;
9bd2ccc2	264
56236ce9	265	fmdbranch->SetAddress(&RecPointsAddress);
c9dd1c4d	266	fmdbranch->GetEntry(0);
56236ce9	267	TIter next(RecPointsAddress) ;
c9dd1c4d	268	AliFMDRecPoint * rp ;
c9dd1c4d	269	while ((rp = static_cast<AliFMDRecPoint*>(next()))) {
198942c7	270
	271	GetRecPointsData(0)->Fill(rp->Edep()/pars->GetEdepMip()) ;
	272
c9dd1c4d	273	}
c9dd1c4d	274
	275	}
	276
	277	//_____________________________________________________________________
	278	void AliFMDQADataMakerRec::StartOfDetectorCycle()
	279	{
	280	// What
	281	// to
	282	// do?
	283	}
	284	//_____________________________________________________________________
d01e669e	285	Int_t AliFMDQADataMakerRec::GetHalfringIndex(UShort_t det,
	286	Char_t ring,
	287	UShort_t board) {
	288
	289	UShort_t iring = (ring == 'I' ? 1 : 0);
	290
	291	Int_t index = (((det-1) << 2) \| (iring << 1) \| (board << 0));
	292
	293	return index-2;
	294
	295	}
	296
	297	//_____________________________________________________________________
	298
	299
c9dd1c4d	300	//
	301	// EOF
	302	//