TimePixClusterProjectionSeparatorProcessor.h

#ifndef TimePixClusterProjectionSeparatorProcessor_h
#define TimePixClusterProjectionSeparatorProcessor_h 1

// C++
#include <iostream>
#include <queue>
#include <string>
#include <sstream>
#include <cmath>
#include <map>
#include <vector>
#include <iomanip>

// LCIO
#include "lcio.h"
#include <EVENT/LCCollection.h>
#include <EVENT/TrackerData.h>
#include <IMPL/LCCollectionVec.h>
#include <IMPL/TrackerHitImpl.h>
#include <IMPL/TrackerDataImpl.h>
#include <IMPL/LCRelationImpl.h>
#include "UTIL/LCTOOLS.h"
#include <IMPL/LCFlagImpl.h>

// Marlin
#include "marlin/Processor.h"
#include "TimePixPixelMode.h"
#include "TimePixMapListener.h"
#include "lccd/ConditionsMap.hh"
#include <marlin/ConditionsProcessor.h>

// Aida
#ifdef MARLIN_USE_AIDA
#include <AIDA/AIDA.h>
#include <marlin/AIDAProcessor.h>
#include <AIDA/IPlotter.h>
#include <AIDA/IPlotterFactory.h>
#include <RAIDA/RAIDAUtil.h>
#include <RAIDA/IBrushStyleROOT.h>
#include <RAIDA/IPlotterFactoryROOT.h>
#include <AIDA/IMarkerStyle.h>

//Root
#include "TROOT.h"
#include "TCanvas.h"
#include "TFile.h"
#include "TH2F.h"
#include "THStack.h"
#include "TLine.h"
#include "TH1F.h"

#endif


namespace AIDA{

  class IHistogram1D;
  class IHistogram2D;

}

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 TimePixClusterProjectionSeparatorProcessor : public marlin::Processor {
        
    public:
        
        /* This method will be called by the marlin package
         * It returns a processor of the currend type (here TimePixClusterProjectionSeparatorProcessor)
         */
        virtual Processor*  newProcessor() { return new TimePixClusterProjectionSeparatorProcessor ; }
        
        /* the default constructor
         * here the processor parameters are registered to the marlin package
         * other initialisation should be placed in the init method
         */
        TimePixClusterProjectionSeparatorProcessor() ;
        
        /* Called at the beginning of the job before anything is read.
         * Use to initialize the processor, e.g. book histograms
         */
        virtual void init() ;
        
        /* Called for every run.
         * in this processor it is not used
         */
        virtual void processRunHeader( lcio::LCRunHeader* run ) ;
        
        /* Called for every event - the working horse.
         * Here the real conversion take place
         */
        virtual void processEvent( lcio::LCEvent * evt ) ; 
        
        /* This method is only called if the check flag is set (default) in the main processor
         * (where the input files are given)
         * It should be used to create check plots which are not needed in a complete analysis
         * (where normaly the check flag is deactivated)
         */
        virtual void check( lcio::LCEvent * evt ) ; 
        
        
        /* Called after data processing for clean up.
         * e.g. saving control histogramms, delete created objects
         * (which are not stored in a collection/event), etc.
         */
        virtual void end() ;

        /* Finds the next local maximum in an AIDA::IHistogram1D,
         * starting at the binnumber given by int startbin
         * and returns the bin number of the maximum 
         * if there is no maximum to be found it returns 0
        */
        virtual int locMaxFinder( int startbin, AIDA::IHistogram1D * h1) ;

        /* Finds the next local minimum in an AIDA::IHistogram1D,
         * starting at the binnumber given by int startbin
         * and returns the bin number of the minimum 
         * if there is no minimum to be found it returns 0
        */
        virtual int locMinFinder( int startbin,AIDA::IHistogram1D * h1) ;

        /* Finds the bin at which the cluster should be separated
         * starting at the local minimum given by the binnumber locmin
         * if the cluster should not be separated it returns 0
        */
        virtual int sepBinFinder( int locmin, AIDA::IHistogram1D * h1) ;

    protected:
        
        /* normally the doxygen documentation is palced before the method/variable
         * but for one line comments in case of variables it may be more convenied for you
         * to place it after the variable. In this case you use an arrow to tell doxygen
         * to which variable the commet belongs to.
         */
        
        std::string _inputColName ;  
        std::string _outputColName ; 
        std::string _outputTrackerDataColName ; 
        std::string _relationsColName ; 

        int _nColumn; 
        int _nRow; 
        int _clusterSizeToSep; 
        int _minResultingClusterSize; 
        bool _keep_hs_sep; 
        
        int _outputIsTransient;

        bool isFirstCollection;

        int _nRun ; 

        double a; 
        double b; 
        
        tpcconddata::TimePixMap * _timePixMap;

    private:
        std::vector<AIDA::IHistogram1D *> hProjectXVec;
        TH1F *_pixelCountHisto;
//      AIDA::IHistogram1D * highTOT;
        AIDA::IHistogram2D * aidahistocluster;
        AIDA::IHistogram2D ** aidaonehit;


    } ;
    
}

#endif