// $Id$
//
//    File: JEventProcessor_AccPi0.h
// Created: Ср. янв. 30 16:59:40 CET 2013
// Creator: alexander (on Linux altro.site 3.4.11-2.16-desktop x86_64)
//

#ifndef _JEventProcessor_AccPi0_
#define _JEventProcessor_AccPi0_

#include <JANA/JEventProcessor.h>

#include <TH2D.h>
#include <TTree.h>
#include <TLorentzVector.h>

class JEventProcessor_AccPi0:public jana::JEventProcessor{
	public:
		JEventProcessor_AccPi0();
		~JEventProcessor_AccPi0();
		const char* className(void){return "JEventProcessor_AccPi0";}

	private:
		jerror_t init(void);						///< Called once at program start.
		jerror_t brun(jana::JEventLoop *eventLoop, int runnumber);	///< Called everytime a new run number is detected.
		jerror_t evnt(jana::JEventLoop *eventLoop, int eventnumber);	///< Called every event.
		jerror_t erun(void);						///< Called everytime run number changes, provided brun has been called.
		jerror_t fini(void);						///< Called after last event of last event source has been processed.
		
		// Output tree
// 		TTree *fTreeAccPi0;
		//
		Int_t fNtrack;
		TH1D *fHistNtrack;
		TH1D *fHistTypeGen;

		
		// Gammas
		Int_t fNGammaGen;
		Int_t fNGamma;
		TH1D *fHistNGamma;
		TH1D *fHistNGammaGen;
		TH1D *fHistThetaGammaGen;
		TH1D *fHistThetaGamma;
		TH1D *fHistThetaGammaGenDet;
		TH1D *fHistPhiGammaGen;
		TH1D *fHistPhiGamma;
		TH1D *fHistPhiGammaGenDet;
		TH1D *fHistEGammaGen;
		TH1D *fHistEGamma;
		TH1D *fHistEGammaGenDet;
		TH1D *fHistFOMgamma;
		TH2D *fHistEGammaVSThetaGamma;
		TH2D *fHistPhiGammaVSThetaGamma;
		TH2D *fHistPhiGammaGenVSThetaGammaGen;
		
		// Pi0
		Int_t fNPi0Gen;
		Int_t fNPi0;
		TH1D *fHistNPi0;
		TH1D *fHistNPi0Gen;
		TH1D *fHistThetaPi0Gen;
		TH1D *fHistThetaPi0;
		TH1D *fHistThetaPi0GenDet;
		TH1D *fHistPhiPi0Gen;
		TH1D *fHistPhiPi0;
		TH1D *fHistPhiPi0GenDet;
		TH2D *fHistPhiPi0VSThetaPi0;
		TH2D *fHistPhiPi0ShiftVSThetaPi0; // +90 degree shift
		TH2D *fHistPhiPi0GenVSThetaPi0Gen;
		TH2D *fHistPhiPi0GenDetVSThetaPi0GenDet;
		TH1D *fHistEPi0Gen;
		TH1D *fHistEPi0;
		TH1D *fHistEPi0GenDet;
		TH1D *fHistMgg;
		TH1D *fHistMggNoCut;
		TH1D *fHistMggNoCutAll;
		TH2D *fHistMggNoCutVsThetaPi0;
		
		// Other particles
		Int_t fNOtherGen;
		Int_t fNOther;
		TH1D *fHistNOther;
		TH1D *fHistNOtherGen;
		TH1D *fHistThetaOtherGen;
		TH1D *fHistPhiOtherGen;
		TH1D *fHistThetaOther;
		TH1D *fHistPhiOther;
		TH1D *fHistEOtherGen;
		TH1D *fHistEOther;
		
		//
		void ResetEvnt();
};

#endif // _JEventProcessor_AccPi0_