VTXDigiProcessor.h

#ifndef VTXDigiProcessor_h
#define VTXDigiProcessor_h 1

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

// STUFF needed for GEAR
#include <marlin/Global.h>
#include <gear/GEAR.h>
#include <gear/VXDParameters.h>
#include <gear/VXDLayerLayout.h>


using namespace lcio ;
using namespace marlin ;


class VTXDigiProcessor : public Processor {
  
 public:
  
  virtual Processor*  newProcessor() { return new VTXDigiProcessor ; }
  
  
  VTXDigiProcessor() ;
  
  virtual void init() ;
  
  virtual void processRunHeader( LCRunHeader* run ) ;
  
  virtual void processEvent( LCEvent * evt ) ; 
  
  
  virtual void check( LCEvent * evt ) ; 
  
  
  virtual void end() ;

  // find phi in correct range, taken from gear::VXDParameters
  double correctPhiRange( double Phi ) const ;  
  
  
 protected:

  std::string _colNameVTX ;
  std::string _colNameSIT ;
  std::string _colNameSET ;
  std::string _outColNameVTX ;
  std::string _outColNameSIT ;
  std::string _outColNameSET ;

  std::vector<int> _activeSETLayers ;

  int _nRun ;
  int _nEvt ;
  int _smearAlongLadders;
  int _debug;
  int _removeDRays;
  float _pointResoRPhi,_pointResoRPhi_VTX,_pointResoRPhi_SIT,_pointResoRPhi_SET;
  float _pointResoZ,_pointResoZ_VTX,_pointResoZ_SIT,_pointResoZ_SET;
  float _momCut;

  FloatVec _vxdEff ;
  std::vector< std::pair<long, long> > _vxdCount ;

  gsl_rng * _rng ;


} ;

#endif