PrimaryIonisationProcessor.h

#ifndef PrimaryIonisationProcessor_h
#define PrimaryIonisationProcessor_h 1

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

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

//C++
#include <string>

//GSL
//#include <gsl/gsl_rng.h>
//#include <gsl/gsl_randist.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 PrimaryIonisationProcessor : public marlin::Processor {
  
 public:
  
  virtual Processor*  newProcessor() { return new PrimaryIonisationProcessor ; }
  
  
  PrimaryIonisationProcessor() ;
  virtual  ~PrimaryIonisationProcessor(){ }
  
  virtual void init() ;
  
  virtual void processRunHeader( lcio::LCRunHeader* run ) ;
  
  virtual void processEvent( lcio::LCEvent * evt ) ; 
  
   
  virtual void SetParticleParameters(lcio::MCParticle* mcparticle);

  virtual void InitializeTrackPoint();

  virtual  void CalculateNumberOfElectrons();
  virtual void CalculateNumberOfClusters();
  virtual void ReadEPerClusterInfoFile();
  virtual void ReadClusterDensityFile();

  virtual void CalculateTrackPointwithBfield();

  virtual void CalculateTrackPointwithoutBfield();
  int GetTotalNumberOfElectrons(){return InTPC;}

  virtual  void CalculateClusterDistance();

  virtual  void CalculateEnergyLoss();

  virtual  void RotateToTrack(double x, double y, double z);

  virtual  void RotateToDelta(double x, double y, double z);

  virtual  void MultipleScattering(double phi, double theta);
  virtual  void SetNextStartPoint();
  
  


 protected:

  // Input/output collection names.
  std::string _inputCollectionName;
  std::string _outputCollectionName;


  std::string clusterdensityfile;  
  std::string eperclusterfile;  

  //For internal bookkeeping
  int _nRun ;
  int _nEvt ;
 
  //TPC parameters
  float _TPCRadius; 
  float _TPCInnerRadius; 
  float _BField;

  float _distance;
  int _transient;

  bool _groupElectrons; // flag whether to group electrons at the same position to one hit

  int Curler;
  //  gsl_rng *r; 

  double aheed[8];
  int numberofelectrons;
  double numberofclusters;
  double clusterdistance;

  std::vector<double> clusterelectrons[2];
  double _clusterprobcutoff; // slightly smaller than where the integrated probability is cut off
  
  int InTPC; 
  // coordinates and variables for track
  double xa,ya,za; 
  double xv,yv,zv; 
  double xend,yend,zend, tend; 
  double x_e,y_e,z_e,t_e; 
  double dx,dy,dz,dr; 

  double dp;

  // particle variables
  double px,py,pz,pt,p,p_const; 
  double energy;
  double mass;
  double charge;
  int pdg;
  double delta_x,delta_y,delta_z;
} ;

}//end of namespace marlintpc

#endif