#ifndef TFgdStoppingCosmicControlSample_h
#define TFgdStoppingCosmicControlSample_h
#include "IGeomInfo.hxx"
#include "IFGDGeom.hxx"
#include "IFGDTimeBinner.hxx"
#include "TH1F.h"

#include <IControlSampleBase.hxx>

// Class to select events based on FGD tracks.  Selection criteria is
// that there must be one (and only one) iso-recon reconstructed track 
// in each FGD (using all FGD hits).
class IFgdStoppingCosmicControlSample : public IControlSampleBase {

public:
  IFgdStoppingCosmicControlSample();
  ~IFgdStoppingCosmicControlSample(){}

private:

  IFGDTimeBinner* fgdTimeBinner;

  bool IsEventSelectedInternal(COMET::ICOMETEvent& event) ;
  void InitHistograms();

  Bool_t fInitialized;

  // variables for upstream/downstream planes
  Int_t usFGD1XModule, usFGD1YModule, dsFGD1XModule, dsFGD1YModule;
  Int_t usFGD2XModule, usFGD2YModule, dsFGD2XModule, dsFGD2YModule;

  // edge cut: do not accept events with hits within edgeCut 
  // of the periphery of the FGD in which the muon stops
  Double_t fEdgeCut;


  // histograms
  TH1F *hNTimeBins;
  TH1F *hNFgdHits, *hNFgd1Hits, *hNFgd2Hits;
  TH1F *hT0, *hT1, *hDT10;

};

#endif