// $Id$
//
//    File: JEventProcessor_CalCal.h
// Created: Mon Feb  2 08:42:18 AM EST 2026
// Creator: staylor (on Linux ifarm2402.jlab.org 5.14.0-611.20.1.el9_7.x86_64 x86_64)
//

/// For more information on the syntax changes between JANA1 and JANA2, visit: https://jeffersonlab.github.io/JANA2/#/jana1to2/jana1-to-jana2

#ifndef _JEventProcessor_CalCal_
#define _JEventProcessor_CalCal_

#include <JANA/JEventProcessor.h>
#include <JANA/Services/JLockService.h> // Required for accessing services
#include <TH1F.h>
#include <TH2F.h>
#include <PID/DChargedTrackHypothesis.h>
#include <ECAL/DECALShower.h>
#include <ECAL/DECALGeometry.h>
#include <DVector3.h>

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

    private:
        void Init() override;                       ///< Called once at program start.
        void BeginRun(const std::shared_ptr<const JEvent>& event) override; ///< Called everytime a new run number is detected.
        void Process(const std::shared_ptr<const JEvent>& event) override;  ///< Called every event.
        void EndRun() override;                     ///< Called everytime run number changes, provided BeginRun has been called.
        void Finish() override;                     ///< Called after last event of last event source has been processed.
  void LeptonAnalysis(const DChargedTrackHypothesis *hyp,TH2F *histo,
		      TH2F *histo2,TH1F *thisto);
  void GammaAnalysis(const DChargedTrackHypothesis *hyp,
		     vector<const DECALShower *>&ecalshowers);

  std::shared_ptr<JLockService> lockService; //Used to access all the services, its value should be set inside Init()

  TH2F *hEoverP_ECAL_Electron,*hEoverP_ECAL_Positron,*hdEdx_CDC,*hdEdx_CDC_cut;
  TH1F *hProtonGammaTimeDiff,*h2GammaMass;
  TH2F *h2GammaMass_vs_E,*hGainScaleFactor,*hChannelDeltaT;
  TH1F *hDeltaT_Electron,*hDeltaT_Positron;
  TH2F *hEfitOverP_Electron,*hEfitOverP_Positron;
  TH1F *hFitProb1,*h2GammaMassImproved;

  int MIN_NUM_ECAL_BLOCKS,MAX_NUM_ECAL_BLOCKS;
  double TRACK_CL_CUT,VAR_E_PAR0,VAR_E_PAR1,VAR_E_PAR2;

  const DECALGeometry *dECALGeom=NULL;
  vector<double>scale_factors,old_gains;
};

#endif // _JEventProcessor_CalCal_