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:
  VTXDigiProcessor(const VTXDigiProcessor&) = delete;
  VTXDigiProcessor& operator=(const VTXDigiProcessor&) = delete;

  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