#ifndef ALIPHOSANALYZE_H
#define ALIPHOSANALYZE_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

/* $Id$ */

//_________________________________________________________________________
// Algorythm class to analyze PHOSv1 events:
// Construct histograms and displays them.
// Use the macro EditorBar.C for best access to the functionnalities
//*--
//*-- Author : Yves Schutz (SUBATECH)

// --- ROOT system ---

class TFile ;
class TH1F ;
class TH2F ;

// --- Standard library ---

// --- AliRoot header files ---

class AliPHOSv1 ;
class AliPHOSGeometry ;
class AliPHOSGetter ;

class AliPHOSAnalyze : public TObject {

public:

  AliPHOSAnalyze() ;              // ctor
  AliPHOSAnalyze(Text_t * name) ; // ctor
  AliPHOSAnalyze(const AliPHOSAnalyze & ana) ; // cpy ctor                   
  virtual ~AliPHOSAnalyze() ;     // dtor

  void DrawRecon(Int_t Nevent= 0,Int_t Nmod = 1,
		 const char* branchName = "PHOSRP",
		 const char* branchTitle = "Default") ; 
  // draws positions of entering of primaries and reconstructed objects in PHOS

  void InvariantMass(const char* RecPartTitle = "Default") ;      // Photons invariant mass distributions

  void EnergyResolution (const char* RecPartTitle = "Default") ;  // analyzes Energy resolution ;

  void PositionResolution(const char* RecPartTitle = "Default") ; // analyzes Position resolution ;

  void Contamination(const char* RecPartTitle = "Default") ;      // Counts contamination of photon spectrum

  void Ls() ; //Prints PHOS-related contents of TreeS, TreeD and TreeR

  void SetEnergyCorrection(const Float_t ecor){fCorrection = ecor ;} 

  AliPHOSAnalyze & operator = (const AliPHOSAnalyze & rvalue)  {
    // assignement operator requested by coding convention but not needed
    abort() ;
    return *this ; 
  }
 
private:

  Float_t CorrectedEnergy(const Float_t ReconstEnergy)const
    {return ReconstEnergy * fCorrection;} 
  //Converts reconstructed energy (energy of the EMCRecPoint) to the energy of primary
  //The coeficient shoud be (and was) calculated usin Erec vs. Eprim plot 
  //(see Energy Resolution function). However, if one change parameters of reconstruction 
  //or geometry, one have to recalculate coefficient!

private:

  Float_t fCorrection ;               //! Conversion coefficient between True and Reconstructed energies
  Int_t   fEvt ;                      //! the evt number being processed 
  TString ffileName ;                 //! the root file that contains the data


ClassDef(AliPHOSAnalyze,1)  // PHOSv1 event analyzis algorithm

};

#endif // AliPHOSANALYZE_H