EMShowerFinder.h

#ifndef EMShowerFinder_h
#define EMShowerFinder_h 1

#include <iostream>
#include <vector>
#include <map>

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

#include <EVENT/LCCollection.h>
#include <IMPL/LCCollectionVec.h>
#include <IMPL/ClusterImpl.h>

#include <EVENT/CalorimeterHit.h>
#include "UTIL/CellIDDecoder.h"
#include <LCRTRelations.h>


#include "ClusterShapes.h"
#include "KITutil.h"
#include "Phys_Geom_Database.h" 


// MarlinCED is only used for debugging
#include <MarlinCED.h>


using namespace lcio ;
using namespace marlin ;


typedef struct {

  CalorimeterHit* ECALHit;
  std::vector<PROTSEED2*> relatedCores;
  std::vector<double> probabilitiesForThisECALHit;
  std::vector<double> distancesToCoresForThisECALHit;
  std::vector<double> estimatedEnergyPerCore;

} ECALHitWithAttributes;



// integer runtime extension which flags CalorimeterHits in EMShowers
struct isPartOfEMShowerCandidate : LCIntExtension<isPartOfEMShowerCandidate> {};



class EMShowerFinder : public Processor {

 public:

  virtual Processor* newProcessor() { return new EMShowerFinder ; }

  EMShowerFinder() ;

  virtual void init() ;
  virtual void processRunHeader( LCRunHeader* run ) ;
  virtual void processEvent( LCEvent * evt ) ;   
  virtual void check( LCEvent * evt ) ; 
  virtual void end() ;

  
 private:

  
 protected:

  std::string _colNameECAL;
  std::string _collectionNameOfEMShowerCandidates;
  std::string _ToClean;
  int _CleanCut;
  int _N;
  vector<float> _miipstep;
  int _MinHit0;
  int _MinHitSplit;
  double _Rcut;
  double _Distcut; 
  double _Coscut;
  double _energyDeviationCut;
  double _probabilityDensityCut;

  int _debugLevel;  
  int _drawOnCED;

  int _nRun;
  int _nEvt;

} ;

#endif