VTXBgClusters.h

#ifndef VTXBgClusters_h
#define VTXBgClusters_h 1

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

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

//#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 VXDGeometry ;

class VTXBgClusters : public Processor {
  
 public:

  VTXBgClusters(const VTXBgClusters&) = delete;
  VTXBgClusters& operator=(const VTXBgClusters&) = delete;
  
  virtual Processor*  newProcessor() { return new VTXBgClusters ; }
  
  
  VTXBgClusters() ;
  
  virtual void init() ;
  
  virtual void processRunHeader( LCRunHeader* run ) ;
  
  virtual void processEvent( LCEvent * evt ) ; 
  
  
  virtual void check( LCEvent * evt ) ; 
  
  
  virtual void end() ;


//   // find phi of point in projection 2D, taken from gear::VXDParameters
//   double getPhiPoint( gear::Vector3D p ) const ;
//   // find phi in correct range, taken from gear::VXDParameters
//   double correctPhiRange( double Phi ) const ;  
  
  
 protected:

  std::string _colNameVTX{};
//   std::string _outColNameVTX ;

  int _nRun{};
  int _nEvt{};
  int _debug{};
  int _mod{};
  int _removeDRays{};
//   float _pointResoRPhi,_pointResoRPhi_VTX;
//   float _pointResoZ,_pointResoZ_VTX;
  float _momCut{};
  float _epi{}; //epitaxial thickness in unit of 50 um (standard thickness)
  float _pitch[6]{}; //pitch in um
  float _it[6]{}; //integration time in ns

  //  gsl_rng * r ;
  VXDGeometry* _vxdGeo{};

#ifdef MARLIN_USE_AIDA
   Hist1DVec _hist1DVec{};
   Hist2DVec _hist2DVec{};
#endif

} ;

#endif