// $Id$
//
//    File: JEventProcessor_pi0.h
// Created: Mon Feb 20 11:37:43 EST 2017
// Creator: zihlmann (on Linux gluon47.jlab.org 2.6.32-431.20.3.el6.x86_64 x86_64)
//

#ifndef _JEventProcessor_pi0_
#define _JEventProcessor_pi0_

#include <JANA/JEventProcessor.h>
using namespace jana;
using namespace std;

#include "TRIGGER/DL1Trigger.h"
#include <TRACKING/DTrackCandidate.h>
#include <TRACKING/DTrackTimeBased.h>
#include <PID/DNeutralParticle.h>
#include <PID/DNeutralShower.h>
#include <START_COUNTER/DSCHit.h>
#include <PID/DBeamPhoton.h>
#include <PID/DKinematicData.h>
#include <PID/DEventRFBunch.h>
#include <DAQ/DCAEN1290TDCHit.h>
#include <DAQ/DCODAROCInfo.h>
#include <RF/DRFTDCDigiTime.h>

#include <TTree.h>
#include <TFile.h>
#include <TH1.h>
#include <TH2.h>
#include <TF1.h>
#include <TMath.h>


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

  double sc_angle_corr;
  vector<vector<DVector3> >sc_pos;   // SC geometry vector
  vector<vector<DVector3> >sc_norm;  
  int TOF_TDC_SHIFT;

  TH1D *eventtype;
  TH1D *BeamPhotonE;
  TH1D *BeamPhotonT;
  TH1D *FoundBeamPhotonE;
  TH1D *FoundBeamPhotonT;
  TH1D *FoundBeamPhotons;
  TH1D *rfTOF;
  TH1D *rfTOFAfter;
  TH1D *SCE[30];
  TH1D *SCT[30];
  TH1D *twophotonm;
  TH1D *twophotonmextended;
  TH1D *rfbunch;
  TH1D *rfbeam;
  TH1D *rfbeamH;
  TH1D *rfbeamM;
  TH1D *rfbeamHNOBCAL;
  TH1D *rfbeamMNOBCAL;
  TH1D *rfbeamNoOffset;
  TH1D *PhotonTBcal;
  TH1D *PhotonTFcal;
  TH1D *FoundPhotonTBcal;
  TH1D *FoundPhotonTFcal;
  TH1D *FoundPhotonEBcal;
  TH1D *FoundPhotonEFcal;
  TH1D *MissingP;
  TH1D *MissingM;
  TH1D *MissingE;
  TH1D *MissingPperp;
  TH1D *PHI;
  TH1D *THETA;
  TH1D *DPHI;
  TH2D *PE1vsPE2;
  TH2D *mictime;
  TH2D *hodtime;
  TH2D *tSC_vs_tgamma;
  TH2D *ESC_vs_Mgamma;
  TH1D *sctimeraw;
  TH1D *scenergyraw;

  int ECounter;
  char RootFile[128];
  TFile *F1;
  TTree *t3;
  //ULong64_t EventNUM;
  int EventNUM;
  float RFT;
  int Nbeam;
  float BeamE[500];
  float BeamT[500];
  int Nneut;
  float NeutralX[40];
  float NeutralY[40];
  float NeutralZ[40];
  float NeutralT[40];
  float NeutralE[40];
  int   DetSys[40];
  int Nstart;
  float SC_E[60];
  float SC_T[60];
  int SC_S[60];

 private:
  jerror_t WriteHistograms(void);
  jerror_t CreateHistogramsAndTree(void);
  jerror_t init(void);
  jerror_t brun(jana::JEventLoop *eventLoop, int32_t runnumber);
  jerror_t evnt(jana::JEventLoop *eventLoop, uint64_t eventnumber);
  jerror_t erun(void);
  jerror_t fini(void);
};

#endif // _JEventProcessor_pi0_