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

 
//____________________________________________________________________________ 
void AliVZEROQADataMakerSim::InitHits()
{
  // create Hits histograms in Hits subdir
  char  timename[10];
  char  chargename[12];
  TH1F  *fhHitTime[72];
  TH1F  *fhHitCharge[72];
  char  texte[25];
 
  // create Hits histograms in Hits subdir
  TH1F * h0 = new TH1F("hHitMultiplicity", "Hits multiplicity distribution in VZERO", 300, 0., 299.) ; 
  h0->Sumw2() ;
  Add2HitsList(h0, 0) ;  
    
  for (Int_t i=0; i<72; i++)
    {
       sprintf(timename,"hHitTime%d",i);
       sprintf(texte,"Hit Time in VZERO cell %d",i);    
       fhHitTime[i] = new TH1F(timename,texte,300,0.,149.);
       
       sprintf(chargename,"hHitCharge%d",i);
       sprintf(texte,"Hit Charge in VZERO cell %d",i);
       fhHitCharge[i]= new TH1F(chargename,texte,1024,0.,1023.);
       
       Add2HitsList( fhHitTime[i],i+1);
       Add2HitsList( fhHitCharge[i],i+1+72);
    }

}

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

  
  char TDCname[10];
  char ADCname[12];
  TH1F *fhDigTDC[64]; 
  TH1F *fhDigADC[64]; 
  char texte[30];

  // 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 TH1F(TDCname,texte,300,0.,149.);
       
       sprintf(ADCname,"hDigitADC%d",i);
       sprintf(texte,"Digit ADC in cell %d",i);
       fhDigADC[i]= new TH1F(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;
	   }
	   Int_t cell = VZEROHit->Cell();
	   GetHitsData(cell+1)->Fill(VZEROHit->Tof());
	   GetHitsData(cell+1+72)->Fill(VZEROHit->Charge());
	}
}


//____________________________________________________________________________

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;
	  }
	  Int_t cell = VZEROHit->Cell();
	  GetHitsData(cell+1)->Fill(VZEROHit->Tof());
	  GetHitsData(cell+1+72)->Fill(VZEROHit->Charge());
	}
    }
  }
}


//____________________________________________________________________________
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("VZERO") ;
    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

}