PulseShapeProcessor.h

#ifndef PULSESHAPEPROCESSOR_H
#define PULSESHAPEPROCESSOR_H

#include <marlin/Processor.h>
#include <lcio.h>
#include <string>
#include <utility>
#include <cstdlib>
#include <map>

//GEAR
#include "gear/TPCParameters.h"

// Aida
#ifdef MARLIN_USE_AIDA
//AIDA
#include <marlin/AIDAProcessor.h>
#endif

//ROOT
#include "TH2.h"

namespace marlintpc
{
class PulseShapeProcessor : public marlin::Processor
{
 public:
  
  virtual Processor*  newProcessor() { return new PulseShapeProcessor; }
  
  PulseShapeProcessor();
  
  virtual void init();
  
  virtual void processRunHeader(lcio::LCRunHeader* run);
  
  virtual void processEvent(lcio::LCEvent* evt);
  
  virtual void check(lcio::LCEvent* evt);
  
  virtual void end();
  
  
 protected:
  /* the place for protected and private member data and functions */
  std::string _inputColName; 

  float _normPulseHistoMax;
  bool _usePolynomialInFit;

  int _nNeighborsUsed;
  bool _useOnlyHighQualityPulses;//< flag to set use of only high quality pulses in average pulse shapes

  TH2* _pulseHighestHisto2D;
  TH2* _pulseNeighborsHisto2D;
  TH2* _normPulseHighestHisto2D;
  TH2* _normPulseNeighborsHisto2D;

  std::map<int, TH2*> _pulseHighestHistoVSCharge;

};
} // namespace marlintpc
#endif // PULSESHAPEPROCESSOR_H