// Author: David Lawrence  June 25, 2004
//
//
// MyProcessor.h
//
/// Example program for a Hall-D analyzer which uses DANA
///

#ifndef _MyProcessor_
#define _MyProcessor_

#include <fstream>
#include <map>
#include <string>
#include <vector>

#include "JANA/JEventProcessor.h"
#include "DVector3.h"

class DKinematicData;
class DGeometry;
class DSCHit;

class MyProcessor:public jana::JEventProcessor
{
 public:
  MyProcessor(std::vector<std::string>& args, 
	      int max_events,
	      string &output_file, 
	      bool batch_mode); 
  
  jerror_t init(void);							
  jerror_t brun(jana::JEventLoop *eventLoop, int runnumber);
  jerror_t evnt(jana::JEventLoop *eventLoop, int eventnumber);		
  jerror_t erun(void);	
  jerror_t fini(void);				       

  bool GetSCPlanesBarrel(std::vector<DVector3>& norms, 
			 std::vector<DVector3>& pts, 
			 std::vector<double>& ang_edges_low);

  jerror_t MatchToSC(const DKinematicData* trk, 
		     const std::vector<const DSCHit*>& sc_hits, 
		     const DSCHit*& assoc_hit) ;

  static std::string NtupleHeader;

 private:
  std::vector<std::string> _args;
  int _max_events;
  std::string _output_file;
  bool _batch_mode;

  std::ofstream* _ofs;
  int _runno;
  std::map<std::string, double> _geomap;

  // Geometry stuff
  DGeometry* _dgeo;
  std::vector<double> _sc_phi_low_edges;
  std::vector<DVector3> _sc_barrel_norms;
  std::vector<DVector3> _sc_barrel_pts;
  std::vector<DVector3> _sc_nose_norms;
  std::vector<DVector3> _sc_nose_pts;
  
  
};

#endif