/** AnaCHI
 * @author P. Russotto
 * @since 13.04.12
 **
 ** Base class for tasks in the cbmroot framework.
 ** Derived classes should implement the Exec method.
 **/

#ifndef AnaCHI_H
#define AnaCHI_H

#include "FairLogger.h"
#include "FairRootManager.h"
#include "FairTask.h"
#include "TRootEvent.h"
#include "TEvtBuffer.h"
#include <TRandom.h>
#include <TRandom2.h>
#include <TH1F.h>
#include <TH2F.h>
#include <TProfile.h>
#include <TF1.h>
#include <TMath.h>
# define UMA 931.5

class TGlobalVar : public TNamed
{
 public:
  int CHImulti;
  int MBallmulti;
  double CHIET12;
  double CHIRP;
  double MBallRP;

 TGlobalVar(){
  CHImulti=-1;
  MBallmulti=-1;
  CHIET12=-1;
  CHIRP=-1000.;
  MBallRP=-1000.  ;

 }
 virtual ~TGlobalVar(){
 }
 ClassDef(TGlobalVar,1);
};

class AnaCHI : public FairTask
{

  public:

    /** Default constructor **/
    AnaCHI();


    /** Standard constructor
    *@param name        Name of task
    *@param iVerbose    Verbosity level
    **/
    AnaCHI(const char* name, Int_t iVerbose = 1);


    /** Destructor **/
    virtual ~AnaCHI();


    /** Virtual method Init **/
    InitStatus Init();
    
    /** Virtual method Exec **/
    void Exec(Option_t* opt);
    
    void Finish();
    void Reset();
    void SetParTask();



  private:
    AnaCHI(const AnaCHI&){};
    AnaCHI& operator=(const AnaCHI&);

    int nevt;

    TClonesArray *fRootMBallEventCopy;
    TClonesArray *fRootLANDEventCopy;
    TClonesArray* fRootCHIEventCopy;
    TClonesArray* fRootCHITSCopy;
    TClonesArray* fRootMBSTSCopy;

    TRootTS *fCHITS;
    TRootTS *fMBSTS;
    TRootCHIEvent *fCHIEvent;
    TRootMBallEvent *fMBallEvent;
    TRootLANDEvent *fLANDEvent;
    unsigned long int *diffTS;
    TH1F *h1_ZTot;
 
    TH2F *h2_MultCHI_MultMB;
    TH2F *h2_PhiCHI_PhiMB;
    TH1F *h1_DPhiCHIMB;
    TH1F *h1_ThetaMB;
    TH1F *h1_ThetaMBMean;
    TH1F *h1_MBMultRatio;

    TH1F *h1_PHIMBR5;
    TH1F *h1_PHIMBR6;
    TH1F *h1_PHIMBR7;
    TH1F *h1_PHIMBR8;

    TH1F *h1_NDETMBR5;
    TH1F *h1_NDETMBR6;
    TH1F *h1_NDETMBR7;
    TH1F *h1_NDETMBR8;
    TH1F *h1_NDETMBR;

    TH2F *h2_Etrans12_ZR7;
    TProfile *hp_Etrans12_ZR7;
    TH2F *h2_Etrans12_Mult;
    TProfile *hp_Etrans12_Mult;
    TH2F *h2_YZ;
    TH2F *h2_Ygb;
    TH2F *h2_Ygbcut;
    TH2F *h2_YptH;
    TH1F *h1_Y;
    TH1F *h1_EMB;   
    TRandom2 *rr;
    TRandom2 *er;
    TRandom2 *qq;

    TH2F *h2_YZ_ETlt200;
    TH2F *h2_Ygb_ETlt200;
    TH2F *h2_YptH_ETlt200;

    TH2F *h2_YZ_ETgt600;
    TH2F *h2_Ygb_ETgt600;
    TH2F *h2_YptH_ETgt600;

    TH2F *h2_YZ_ETbt200400;
    TH2F *h2_Ygb_ETbt200400;
    TH2F *h2_YptH_ETbt200400;

    TH2F *h2_YZ_ETbt400600;
    TH2F *h2_Ygb_ETbt400600;
    TH2F *h2_YptH_ETbt400600;

    TH1F *h1_rp;
    TH1F *h1_rpETgt600;
    TH1F *h1_rpETbt200400;
    TH1F *h1_rpETbt400600;
    TH1F *h1_rpETlt200;

    TH1F *h1_rpres180;
    TH1F *h1_rpres180ETgt600;
    TH1F *h1_rpres180ETbt200400;
    TH1F *h1_rpres180ETbt400600;
    TH1F *h1_rpres180ETlt200;
    TH1F *h1_rpres90;
    TH1F *h1_Etrans12;

    TH2F *h2_Ygb_1HS;
    TH2F *h2_Ygb_2HS;
    TH2F *h2_Ygb_3HS;
    TH2F *h2_Ygb_1HPT;
 
    TH2F *h2_Ygb_4HeS;
    TH2F *h2_Ygb_3HeS;
    TH2F *h2_Ygb_4HePT;
 
    TH2F *h2_Ygb_LiS;
    TH2F *h2_Ygb_LiPT;

    TH2F *h2_Ygb_BeS;
    TH2F *h2_Ygb_BePT;


    int Ntot[5];
    int Nback[5];
    double ratiomm;
    TGlobalVar *fGlobalVar;


    ClassDef(AnaCHI,1);
    
};


#endif