DriftProcessor.h

#ifndef DriftProcessor_h
#define DriftProcessor_h 1


//LCIO
#include <EVENT/SimTrackerHit.h>
#include <lcio.h>

//MARLIN
#include "marlin/Processor.h"

//C++
#include <string>

//GSL
#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>

#include "TPCGas.h"

namespace marlintpc{


/* The following comment is tranfered to doxygen documentation
 * which begins with a second asterix, or in case of on line comment an extra third slash 
 */

class DriftProcessor : public marlin::Processor {
  
 public:
  
  virtual Processor*  newProcessor() { return new DriftProcessor ; }
  
  DriftProcessor() ;

  virtual  ~DriftProcessor(){}  

  virtual void init() ;
  
  virtual void processRunHeader( lcio::LCRunHeader* run ) ;
  
  virtual void processEvent( lcio::LCEvent * evt ) ; 


  virtual void end();
  

protected:

  void MakeElectronDriftCoordinates();

  void GetElectronCoordinates(SimTrackerHit *);


  //Input/output collection names
  std::string _inputCollectionName;
  std::string _outputCollectionName;
  std::string _tpcConditionsColName;
  std::string _gas;
  
  float _EField;
  float _BField;

  int _transient;
 
  gsl_rng *r;
  // coordinates and variables for track
  double x_e,y_e,z_e,t_e; // e- coordinates
  // particle variables
  //int numberofelectrons;

  //TPC parameter
  double _TPCLength;
  double _TPCOuterRadius;
  double _TPCInnerRadius;
  double _distance;

  //Coordinates of the drifted electron.
  float xdrift,ydrift,zdrift,tdrift; 
  float _vdrift;  
  float _difflong; 
  float _difftrans; 
  double _meanAmp;

  TPCGas *_gasClass;

  double _ionisationEnergy;
  bool _fullCurlerBackwardCompatibility; 

  bool _isFirstEvent; 
  
  double _attachment; 
};

}//end of namespace marlintpc

#endif