SimpleDiscDigiProcessor.h

/* -*- Mode: C++; tab-width: 2; indent-tabs-mode: nil; c-basic-offset: 2 -*- */
#ifndef SimpleDiscDigiProcessor_h
#define SimpleDiscDigiProcessor_h 1

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

using namespace lcio ;
using namespace marlin ;

class FTDLayer ;

namespace EVENT{
  class SimTrackerHit;
}

class SimpleDiscDigiProcessor : public Processor {
  
public:
  
  virtual Processor*  newProcessor() { return new SimpleDiscDigiProcessor ; }
  
  
  SimpleDiscDigiProcessor() ;
  
  virtual void init() ;
  
  virtual void processRunHeader( LCRunHeader* run ) ;
  
  virtual void processEvent( LCEvent * evt ) ; 
  
  
  virtual void check( LCEvent * evt ) ; 
  
  
  virtual void end() ;
  
  
protected:
  
  void process_hits_loi( LCEvent* evt, LCCollection* STHcol );

  void process_hits_new( LCEvent* evt, LCCollection* STHcol );

  std::string _inColName ;
  std::string _outColName ;
  std::string _outRelColName ;
  
  int _sub_det_id ;
  
  int _nRun ;
  int _nEvt ;
  
  bool _SimHits_encoded_with_cellID;

  bool _use_FTDLayerLayout_from_GEAR;

  float _pointReso;
  bool _keepHitsFromDeltas;

  std::vector< double > _FTDZCoordinate;
  std::vector< double > _diskInnerRadius;
  std::vector< double > _diskOuterRadius;
  
  bool hasCorrectZPos ( SimTrackerHit* hit );
  
  int getPetalNumber ( int layer , double x , double y );
  int getSensorNumber ( int layer , double x , double y );
  
  int _petalsPerDisk;
  int _sensorsPerPetal;
  
  // gsl random number generator
  gsl_rng * r ;
  
  
} ;

#endif