GEMProcessor.h

#ifndef GEMProcessor_h
#define GEMProcessor_h 1

//C++
#include <string>
#include <vector>
#include "TF1.h"

//LCIO
#include <IMPL/LCCollectionVec.h>
#include <marlin/Processor.h>
#include <lcio.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 GEMProcessor : public marlin::Processor {
  
 public:
  
  virtual Processor *newProcessor(){ return new GEMProcessor ; }
   
  GEMProcessor() ;

  virtual ~GEMProcessor(){}
  
  virtual void init() ;
  
  virtual void processRunHeader(lcio::LCRunHeader* run ) ;
  
  virtual void processEvent(lcio::LCEvent * evt ) ; 
    
  void LoadGEMSettings();

  int GEM_Collection(int GEMNr, int enumber);

  double GEM_Gain(int GEMNr, int enumber);

  int GEM_Extraction(int GEMNr, int enumber);

  double CalculateAttachment(double E);

  double AttachmentParam(double E, double a0, double a1, double a2, 
                         double a3, double a4, double a5);

  int CalculateNrElectrons(double attachment1, double attachment2, 
                           double attachment3, int charge);


  virtual void end() ;

 
protected:

  std::string _inputCollectionName; 
  std::string _outputCollectionName;
  std::string _tpcConditionsColName;
  std::string _altTPCConditionsName;

  std::string paramfile;
  std::string gas; 
  
  float UGEM1; 
  float UGEM2; 
  float UGEM3; 

  float ETrans1; 
  float ETrans2; 
  float EInd;    

  float _gainFactor;

  int _transient;

  float _driftField;

  double Extr[3];
  double Coll[3];
  double Gain[3];
  double gainmean[3];

  /* Attachment 
   */
  double A1, A2, A3;

  //For internal bookkeeping
  int _nRun;
  int _nEvt;

  double _meanAmp;

  double _ionisationEnergy;
  bool _fullCurlerBackwardCompatibility; 
  bool _GEMGainPolya;//< Flag to swicth between polya and exponential gain in GEMs
  //  TF1 *_gainPolya;///<Polya function for gain in GEMs, if option selected
  // store one for each GEM, to avoid all the recalculations in TF1::GetRandom after changing parameters
  std::map<int, TF1*> _gainPolya;
  double _polyaThetaP1;

}; //end of class definition

} //end of namespace marlintpc

#endif