IonsInVoxelsProcessor.h

#ifndef IonsInVoxelsProcessor_h
#define IonsInVoxelsProcessor_h 1

//MARLINTPC
#include "VoxelTPC.h"

//LCIO
#include <lcio.h>

//MARLIN
#include <marlin/Processor.h>

//C++
#include <string>

//GEAR
#include <gear/PadRowLayout2D.h>

#include <TTree.h>
#include <TFile.h>


namespace marlintpc{

class VoxelTPC;

/* 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 IonsInVoxelsProcessor : public marlin::Processor {

 
 public:
  
  virtual Processor*  newProcessor() { return new IonsInVoxelsProcessor ; }  
  
  IonsInVoxelsProcessor() ;
  virtual ~IonsInVoxelsProcessor();
  
  virtual void init() ;
  
  virtual void processRunHeader(lcio::LCRunHeader* run ) ;
  
  virtual void processEvent(lcio::LCEvent * evt ) ; 

  virtual void end();
    

  int getPad(gear::PadRowLayout2D const* padLayout, SimTrackerHit* hit);
  
  
protected:

  //input and output collections
  std::string _inputIonsColName;
  std::string _outputVoxelColName;

  //padLayout pointer
  gear::PadRowLayout2D const * _padLayout;

  //the pad layout type (cartesian or polar)
  int _padLayoutType;

  //pointer to the VoxelTPC which stores the ions
  marlintpc::VoxelTPC* _voxelTPC;

  //input parameters from steering file
  bool _outputIsPersistent; 
  int _outputType;
  int _nEvt;
  double _binLength;
  double _ionisationEnergy;
  int _writeAfterNEvents;
  std::string _rootFileName;

  //TPC parameters
  double _maxDriftLength;

  // Parameters for the root tree
  std::vector <double> _leaf_chargePerVolume; 
  std::vector <int> _leaf_charge;
  std::vector <int> _leaf_padIndex; 
  std::vector <int> _leaf_zIndex; 
  std::vector <double> _leaf_x,_leaf_y,_leaf_z; 

  TTree *_chargeInVoxelsTree; // the root tree
  TFile * _rootFile; // the root file

  void fillTree(LCCollectionVec* outputVoxelCollection);

  bool _useNegativeHalf;
};

}//end of namespace marlintpc

#endif