VTXNoiseClusters.h

#ifndef VTXNoiseClusters_h
#define VTXNoiseClusters_h 1


//#ifdef USE_ROOT 
// we need some root histograms with cluster size distribution for this processor

#include "marlin/Processor.h"
#include "marlin/EventModifier.h"
#include "lcio.h"
#include <string>
#include <vector>
#include <gsl/gsl_rng.h>

#include "TH2F.h"
#include "TFile.h"

#ifdef MARLIN_USE_AIDA
#include <AIDA/AIDA.h>
typedef std::vector< AIDA::IHistogram1D* > Hist1DVec ;
typedef std::vector< AIDA::IHistogram2D* > Hist2DVec ;
#endif

using namespace lcio ;
using namespace marlin ;

class VXDGeometry ;


class VTXNoiseClusters : public  Processor, public EventModifier{
  
 public:
  
  virtual Processor*  newProcessor() { return new VTXNoiseClusters ; }
  
  
  VTXNoiseClusters() ;
  
  virtual const std::string & name() const { return Processor::name() ; }

  virtual void init() ;
  
  virtual void processRunHeader( LCRunHeader* run ) ;
  
  void modifyEvent( LCEvent * evt ) ; 
  
  virtual void processEvent( LCEvent * evt ) ; 
  
  
  virtual void check( LCEvent * evt ) ; 
  
  
  virtual void end() ;


 protected:

  std::string _colNameVTX ;
  FloatVec _densities ;
  StringVec _rootNames ;  
//   float _pointResoRPhiVTX ;
//   float _pointResoZVTX ;
  int   _ranSeed  ;

  int _nRun ;
  int _nEvt ;

  gsl_rng* _rng ;
  VXDGeometry* _vxdGeo ;

  std::vector<TH2F*> _hist ;
  TFile* _hfile ;

#ifdef MARLIN_USE_AIDA
   Hist1DVec _hist1DVec ;
   Hist2DVec _hist2DVec ;
#endif

} ;

#endif
//#endif // USE_ROOT