FindLineTrajectory.h

// -*- mode: c++; mode: auto-fill; mode: flyspell-prog; -*-


#ifndef FindLineTrajectory_h
#define FindLineTrajectory_h 1

#include "luiza/Processor.h"

using namespace gloria;

namespace luiza {

  typedef long long int type_fitcount;



  class FindLineTrajectory : public Processor {

  public:




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

    FindLineTrajectory() ;


    virtual void init() ;



    virtual void startAnalysis( GloriaDataContainer* ) ; 


    virtual void processData( GloriaDataContainer* data ) ;
    
    

    virtual void check( GloriaDataContainer* data ) ; 
    
    

    virtual void concludeAnalysis( GloriaDataContainer* ) ;


    virtual void end() ;

  protected:

    bool initFit(GloriaDataContainer*);


    void endFit();


    double MatrixFit();


    double SingleFit();


    double NominalFit();

    int DoAnalFit(double * pos, double *err, double velocity=0.);


    double GetFitChi2();

    int GaussjSolve(double * alfa, double * beta, int n);


    // Global processor parameters

    std::string _inputCollectionName ;

    std::string _outputTrackCollectionName ;

    std::string _outputPositionCollectionName ;


    // Parameters of the algorithm

    int _maxFrameObjects;

    double _chi2Max ;

    std::vector<std::string> _positionColumns;

    std::vector<float> _positionErrors;

    std::string _timeRecord;

    bool _searchMultipleTracks;

    bool _allowAmbiguousObjects;

    int  _maximumAmbiguousObjects;

    int _allowMissingObservations;

    double _missingObservationPenalty;

    bool   _useVelocityConstraint ;

    double _velocitySpread;

    double _raVelocity;

    double _decVelocity;

    bool _usePreselection;


    // Setup description

    int _allFrameCount;
    int _goodFrameCount;

    int * _frameSort;
    int * _frameID;
    double * _frameTime;
    double * _frameResolution;

    double _timeOffset;

    // Internal processor variables
    // ----------------------------

    // Arrays for selecting different hit combinations

    int * _frameHits;
    int * _frameChoice;
    type_fitcount * _frameMod;

    // Fitting algorithm arrays

    double * _frameX  ;
    double * _frameEx ;
    double * _frameY  ;
    double * _frameEy ;

    double * _fitX  ;
    double * _fitEx ;
    double * _fitY  ;
    double * _fitEy ;

    int _vecDim;
    int _arrayDim;

    double * _parVec;
    double * _fitArray;
    double * _parVecX;
    double * _parVecY;
    double * _parErrX;
    double * _parErrY;

    double * _nominalFitArrayX ;
    double * _nominalErrorX ;
    double * _nominalFitArrayY ;
    double * _nominalErrorY ;

    // Track table

    GloriaFitsTable * _trackTable;

    // few counter to show the final summary

    int _nLoops;

    int _nFrames;

    int _nTracks;



  } ;

}

#endif