[u/mrichter/AliRoot.git] / VZERO / AliVZEROQADataMakerSim.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: AliVZEROQADataMakerSim.cxx 23123 2007-12-18 09:08:18Z hristov $ */

//---
//  Produces the data needed to calculate the quality assurance. 
//  All data must be mergeable objects.
//  Author : BC
//---

// --- ROOT system ---
#include <TClonesArray.h>
#include <TFile.h> 
#include <TH1F.h> 
#include <TDirectory.h>
// --- Standard library ---

// --- AliRoot header files ---
#include "AliESDEvent.h"
#include "AliLog.h"
#include "AliVZEROdigit.h" 
#include "AliVZEROhit.h"
#include "AliVZEROQADataMakerSim.h"
#include "AliQAChecker.h"

ClassImp(AliVZEROQADataMakerSim)
           
//____________________________________________________________________________ 
  AliVZEROQADataMakerSim::AliVZEROQADataMakerSim() : 
  AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kVZERO), "VZERO Quality Assurance Data Maker")

{
  // constructor

  
}

//____________________________________________________________________________ 
AliVZEROQADataMakerSim::AliVZEROQADataMakerSim(const AliVZEROQADataMakerSim& qadm) :
  AliQADataMakerSim() 
{
  //copy constructor 

  SetName((const char*)qadm.GetName()) ; 
  SetTitle((const char*)qadm.GetTitle()); 
}

//__________________________________________________________________
AliVZEROQADataMakerSim& AliVZEROQADataMakerSim::operator = (const AliVZEROQADataMakerSim& qadm )
{
  // Assign operator.
  this->~AliVZEROQADataMakerSim();
  new(this) AliVZEROQADataMakerSim(qadm);
  return *this;
}
//____________________________________________________________________________
void AliVZEROQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray ** list)
{
  //Detector specific actions at end of cycle
  // do the QA checking
  AliQAChecker::Instance()->Run(AliQAv1::kVZERO, task, list) ;
}

 
//____________________________________________________________________________ 
void AliVZEROQADataMakerSim::InitHits()
{
 
  // create Hits histograms in Hits subdir
  
  TH1I * h0 = new TH1I("hHitMultiplicity", "Hit multiplicity distribution in VZERO", 300, 0, 299) ; 
  h0->Sumw2() ;
  Add2HitsList(h0, 0) ;  
  
  TH1I * h1 = new TH1I("hHitCellNumber", "Hit cell distribution in VZERO", 80, 0, 79) ; 
  h1->Sumw2() ;
  Add2HitsList(h1, 1) ;  
    
}

//____________________________________________________________________________ 
void AliVZEROQADataMakerSim::InitDigits()
{
  // create Digits histograms in Digits subdir

  
  char TDCname[100];
  char ADCname[100];
  TH1I *fhDigTDC[64]; 
  TH1I *fhDigADC[64]; 
  char texte[100];

  // create Digits histograms in Digits subdir
  TH1I * h0 = new TH1I("hDigitMultiplicity", "Digits multiplicity distribution in VZERO", 100, 0, 99) ; 
  h0->Sumw2() ;
  Add2DigitsList(h0, 0) ;
  
  for (Int_t i=0; i<64; i++)
    {
       sprintf(TDCname, "hDigitTDC%d", i);
       sprintf(texte,"Digit TDC in cell %d",i);    
       fhDigTDC[i] = new TH1I(TDCname,texte,300,0.,149.);
       
       sprintf(ADCname,"hDigitADC%d",i);
       sprintf(texte,"Digit ADC in cell %d",i);
       fhDigADC[i]= new TH1I(ADCname,texte,1024,0.,1023.);
       
       Add2DigitsList(fhDigTDC[i],i+1);
       Add2DigitsList(fhDigADC[i],i+1+64);  
     }  
}


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

    GetHitsData(0)->Fill(hits->GetEntriesFast()) ;    // fills Hit multiplicity
    Int_t nhits = hits->GetEntriesFast();
    for (Int_t ihit=0;ihit<nhits;ihit++) 
	{
	   AliVZEROhit  * VZEROHit   = (AliVZEROhit*) hits->UncheckedAt(ihit);
	   if (!VZEROHit) {
 	      AliError("The unchecked hit doesn't exist");
	      break;
	   }
	   GetHitsData(1)->Fill(VZEROHit->Cell());
	}
}


//____________________________________________________________________________

void AliVZEROQADataMakerSim::MakeHits(TTree *hitTree)
{
  //fills QA histos for Hits
  TClonesArray * hits = new TClonesArray("AliVZEROhit", 1000);
  
  TBranch * branch = hitTree->GetBranch("VZERO") ;
  if ( ! branch ) {
    AliWarning("VZERO branch in Hit Tree not found") ;
  } else {

   if (branch) {
      branch->SetAddress(&hits);
    }else{
      AliError("Branch VZERO hit not found");
      exit(111);
    } 
    Int_t ntracks    = (Int_t) hitTree->GetEntries();
    
    if (ntracks<=0) return;
    // Start loop on tracks in the hits containers
    for (Int_t track=0; track<ntracks;track++) {
      branch->GetEntry(track);
      GetHitsData(0)->Fill(hits->GetEntriesFast()) ;    // fills Hit multiplicity
      Int_t nhits = hits->GetEntriesFast();
      for (Int_t ihit=0;ihit<nhits;ihit++) 
	{
	  AliVZEROhit  * VZEROHit   = (AliVZEROhit*) hits->UncheckedAt(ihit);
	  if (!VZEROHit) {
 	    AliError("The unchecked hit doesn't exist");
	    break;
	  }
	  GetHitsData(1)->Fill(VZEROHit->Cell());	 
	}
    }
  }
}


//____________________________________________________________________________
void AliVZEROQADataMakerSim::MakeDigits(TClonesArray * digits)
{
  // makes data from Digits

    GetDigitsData(0)->Fill(digits->GetEntriesFast()) ; 
    TIter next(digits) ; 
    AliVZEROdigit *VZERODigit ; 
    while ( (VZERODigit = dynamic_cast<AliVZEROdigit *>(next())) ) {
         Int_t   PMNumber  = VZERODigit->PMNumber();         
         GetDigitsData(PMNumber +1)->Fill( VZERODigit->Time()) ;    // in 100 of picoseconds
	 GetDigitsData(PMNumber +1+64)->Fill( VZERODigit->ADC()) ;
    }  
}


//____________________________________________________________________________
void AliVZEROQADataMakerSim::MakeDigits(TTree *digitTree)
{
    // makes data from Digit Tree
	
    TClonesArray * digits = new TClonesArray("AliVZEROdigit", 1000) ; 

    TBranch * branch = digitTree->GetBranch("VZERODigit") ;
    if ( ! branch ) {
         AliWarning("VZERO branch in Digit Tree not found") ; 
    } else {
         branch->SetAddress(&digits) ;
         branch->GetEntry(0) ; 
         MakeDigits(digits) ; 
    }  
}


//____________________________________________________________________________
void AliVZEROQADataMakerSim::StartOfDetectorCycle()
{
  //Detector specific actions at start of cycle

}
Commit	Line	Data
508b9fc0	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: AliVZEROQADataMakerSim.cxx 23123 2007-12-18 09:08:18Z hristov $ */
	18
	19	//---
	20	// Produces the data needed to calculate the quality assurance.
	21	// All data must be mergeable objects.
	22	// Author : BC
	23	//---
	24
	25	// --- ROOT system ---
	26	#include <TClonesArray.h>
	27	#include <TFile.h>
	28	#include <TH1F.h>
	29	#include <TDirectory.h>
	30	// --- Standard library ---
	31
	32	// --- AliRoot header files ---
	33	#include "AliESDEvent.h"
	34	#include "AliLog.h"
	35	#include "AliVZEROdigit.h"
	36	#include "AliVZEROhit.h"
	37	#include "AliVZEROQADataMakerSim.h"
	38	#include "AliQAChecker.h"
	39
	40	ClassImp(AliVZEROQADataMakerSim)
	41
	42	//____________________________________________________________________________
	43	AliVZEROQADataMakerSim::AliVZEROQADataMakerSim() :
4e25ac79	44	AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kVZERO), "VZERO Quality Assurance Data Maker")
508b9fc0	45
	46	{
	47	// constructor
	48
	49
	50	}
	51
	52	//____________________________________________________________________________
	53	AliVZEROQADataMakerSim::AliVZEROQADataMakerSim(const AliVZEROQADataMakerSim& qadm) :
	54	AliQADataMakerSim()
	55	{
	56	//copy constructor
	57
	58	SetName((const char*)qadm.GetName()) ;
	59	SetTitle((const char*)qadm.GetTitle());
	60	}
	61
	62	//__________________________________________________________________
	63	AliVZEROQADataMakerSim& AliVZEROQADataMakerSim::operator = (const AliVZEROQADataMakerSim& qadm )
	64	{
	65	// Assign operator.
	66	this->~AliVZEROQADataMakerSim();
	67	new(this) AliVZEROQADataMakerSim(qadm);
	68	return *this;
	69	}
	70	//____________________________________________________________________________
4e25ac79	71	void AliVZEROQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray ** list)
508b9fc0	72	{
	73	//Detector specific actions at end of cycle
	74	// do the QA checking
4e25ac79	75	AliQAChecker::Instance()->Run(AliQAv1::kVZERO, task, list) ;
508b9fc0	76	}
	77
	78
	79	//____________________________________________________________________________
	80	void AliVZEROQADataMakerSim::InitHits()
	81	{
508b9fc0	82
508b9fc0	83	// create Hits histograms in Hits subdir
095fab41	84
095fab41	85	TH1I * h0 = new TH1I("hHitMultiplicity", "Hit multiplicity distribution in VZERO", 300, 0, 299) ;
508b9fc0	86	h0->Sumw2() ;
508b9fc0	87	Add2HitsList(h0, 0) ;
095fab41	88
	89	TH1I * h1 = new TH1I("hHitCellNumber", "Hit cell distribution in VZERO", 80, 0, 79) ;
	90	h1->Sumw2() ;
	91	Add2HitsList(h1, 1) ;
508b9fc0	92
508b9fc0	93	}
	94
	95	//____________________________________________________________________________
	96	void AliVZEROQADataMakerSim::InitDigits()
	97	{
	98	// create Digits histograms in Digits subdir
	99
	100
a3f8715e	101	char TDCname[100];
a3f8715e	102	char ADCname[100];
095fab41	103	TH1I *fhDigTDC[64];
095fab41	104	TH1I *fhDigADC[64];
a3f8715e	105	char texte[100];
508b9fc0	106
508b9fc0	107	// create Digits histograms in Digits subdir
095fab41	108	TH1I * h0 = new TH1I("hDigitMultiplicity", "Digits multiplicity distribution in VZERO", 100, 0, 99) ;
508b9fc0	109	h0->Sumw2() ;
	110	Add2DigitsList(h0, 0) ;
	111
	112	for (Int_t i=0; i<64; i++)
	113	{
	114	sprintf(TDCname, "hDigitTDC%d", i);
	115	sprintf(texte,"Digit TDC in cell %d",i);
095fab41	116	fhDigTDC[i] = new TH1I(TDCname,texte,300,0.,149.);
508b9fc0	117
	118	sprintf(ADCname,"hDigitADC%d",i);
	119	sprintf(texte,"Digit ADC in cell %d",i);
095fab41	120	fhDigADC[i]= new TH1I(ADCname,texte,1024,0.,1023.);
508b9fc0	121
	122	Add2DigitsList(fhDigTDC[i],i+1);
	123	Add2DigitsList(fhDigADC[i],i+1+64);
	124	}
	125	}
	126
	127
	128	//____________________________________________________________________________
	129	void AliVZEROQADataMakerSim::MakeHits(TClonesArray * hits)
	130	{
	131	//make QA data from Hits
	132
	133	GetHitsData(0)->Fill(hits->GetEntriesFast()) ; // fills Hit multiplicity
	134	Int_t nhits = hits->GetEntriesFast();
	135	for (Int_t ihit=0;ihit<nhits;ihit++)
	136	{
	137	AliVZEROhit * VZEROHit = (AliVZEROhit*) hits->UncheckedAt(ihit);
	138	if (!VZEROHit) {
	139	AliError("The unchecked hit doesn't exist");
	140	break;
	141	}
095fab41	142	GetHitsData(1)->Fill(VZEROHit->Cell());
508b9fc0	143	}
	144	}
	145
	146
	147	//____________________________________________________________________________
	148
	149	void AliVZEROQADataMakerSim::MakeHits(TTree *hitTree)
	150	{
	151	//fills QA histos for Hits
	152	TClonesArray * hits = new TClonesArray("AliVZEROhit", 1000);
	153
	154	TBranch * branch = hitTree->GetBranch("VZERO") ;
	155	if ( ! branch ) {
	156	AliWarning("VZERO branch in Hit Tree not found") ;
	157	} else {
	158
	159	if (branch) {
	160	branch->SetAddress(&hits);
	161	}else{
	162	AliError("Branch VZERO hit not found");
	163	exit(111);
	164	}
	165	Int_t ntracks = (Int_t) hitTree->GetEntries();
	166
	167	if (ntracks<=0) return;
	168	// Start loop on tracks in the hits containers
	169	for (Int_t track=0; track<ntracks;track++) {
	170	branch->GetEntry(track);
	171	GetHitsData(0)->Fill(hits->GetEntriesFast()) ; // fills Hit multiplicity
	172	Int_t nhits = hits->GetEntriesFast();
	173	for (Int_t ihit=0;ihit<nhits;ihit++)
	174	{
	175	AliVZEROhit * VZEROHit = (AliVZEROhit*) hits->UncheckedAt(ihit);
	176	if (!VZEROHit) {
	177	AliError("The unchecked hit doesn't exist");
	178	break;
	179	}
095fab41	180	GetHitsData(1)->Fill(VZEROHit->Cell());
508b9fc0	181	}
	182	}
	183	}
	184	}
	185
	186
	187	//____________________________________________________________________________
	188	void AliVZEROQADataMakerSim::MakeDigits(TClonesArray * digits)
	189	{
	190	// makes data from Digits
	191
	192	GetDigitsData(0)->Fill(digits->GetEntriesFast()) ;
	193	TIter next(digits) ;
	194	AliVZEROdigit *VZERODigit ;
	195	while ( (VZERODigit = dynamic_cast<AliVZEROdigit *>(next())) ) {
	196	Int_t PMNumber = VZERODigit->PMNumber();
	197	GetDigitsData(PMNumber +1)->Fill( VZERODigit->Time()) ; // in 100 of picoseconds
	198	GetDigitsData(PMNumber +1+64)->Fill( VZERODigit->ADC()) ;
	199	}
	200	}
	201
	202
	203	//____________________________________________________________________________
	204	void AliVZEROQADataMakerSim::MakeDigits(TTree *digitTree)
	205	{
	206	// makes data from Digit Tree
	207
095fab41	208	TClonesArray * digits = new TClonesArray("AliVZEROdigit", 1000) ;
508b9fc0	209
095fab41	210	TBranch * branch = digitTree->GetBranch("VZERODigit") ;
508b9fc0	211	if ( ! branch ) {
	212	AliWarning("VZERO branch in Digit Tree not found") ;
	213	} else {
	214	branch->SetAddress(&digits) ;
	215	branch->GetEntry(0) ;
	216	MakeDigits(digits) ;
	217	}
	218	}
	219
	220
	221	//____________________________________________________________________________
	222	void AliVZEROQADataMakerSim::StartOfDetectorCycle()
	223	{
	224	//Detector specific actions at start of cycle
	225
	226	}