#ifndef _myAlg_H
#define _myAlg_H 1

#include "IAlgBase.h"
#include "TF1.h"

#include <vector>
#include <utility>
#include <fstream>
//#include "Track.h"
#include "TRandom.h" 
//#include "TrackStore.h"
#include "TObjArray.h"
#include "math.h"

class TFile;
class TF1;
class TH1F;
class TH2F;
class IMilleTool;
using namespace std;
class myAlg: public IAlgBase {

 public:

  /** constructer */
  myAlg();

  /** destructer */
  virtual ~myAlg();

  virtual bool initialize(const jobParams& tparams);

  virtual bool execute();

  virtual bool finalize();


protected:

  /** parameters for this job */	 
   jobParams params() const;

  /* /\** root file *\/ */
  /* TFile* rootFile() const;  */
  std::ofstream output;
  IMilleTool* m_align;

 private:

  //  TrackStore *my_tracks;
  //  Tracks input_tracks;

  double offset(const double sizeOff) const;
 
  jobParams m_params;
  // TFile* m_rootFile;
  //  TF1* m_GaussGenerator;

  /* // histos */
  /* TObjArray *m_mis_aft; */
  /* TObjArray *m_pullHistos; */
};

inline jobParams myAlg::params() const{
  return m_params;
} 


/* inline TFile* myAlg::rootFile() const{ */
/*   return m_rootFile; */
/* } */


/* inline double myAlg::offset(const double sizeOff) const{ */
/*   return (m_GaussGenerator->GetRandom()*sizeOff); */
/* } */

#endif