MarlinKalTest.h

#ifndef MarlinKalTest_h
#define MarlinKalTest_h

#include "MarlinTrk/IMarlinTrkSystem.h"

#ifdef MARLINTRK_DIAGNOSTICS_ON
#include "MarlinTrk/DiagnosticsController.h"
#endif

#include "gear/GearMgr.h"

//LCIO:
#include "lcio.h"
#include "UTIL/BitField64.h" 
#include "UTIL/LCTOOLS.h"
#include <LCRTRelations.h>

#include "streamlog/streamlog.h"

#include "TObjArray.h"
#include "TVector3.h"

#include <cmath>
#include <vector>


class TKalDetCradle ;
class TVKalDetector ;
class ILDVMeasLayer ;
class THelicalTrack ;

class ILDCylinderMeasLayer;

namespace EVENT{
  class TrackerHit ;
}

namespace MarlinTrk{
  
  class MarlinKalTest : public MarlinTrk::IMarlinTrkSystem {
    
  public:
    
    friend class MarlinKalTestTrack;
    
    // define some configuration constants
    static const bool FitBackward   = kIterBackward ;
    static const bool FitForward    = kIterForward ;
    static const bool OrderOutgoing  = true ;
    static const bool OrderIncoming  = false ;
    
    
    MarlinKalTest( const gear::GearMgr& gearMgr) ;
    
    ~MarlinKalTest() ;
    
    void init() ; 
    
    virtual std::string name() { return "KalTest" ; }

    MarlinTrk::IMarlinTrack* createTrack()  ;
    
    
  protected:
    
    void includeMultipleScattering( bool on )  ;
    
    void includeEnergyLoss( bool on )  ;
    
    void storeActiveMeasurementModuleIDs(TVKalDetector* detector);  
    
    void getSensitiveMeasurementModules( int detElementID, std::vector< const ILDVMeasLayer *>& measmodules) const; 
    
    void getSensitiveMeasurementModulesForLayer( int layerID, std::vector<const ILDVMeasLayer *>& measmodules) const;
    
    //  void init(bool MSOn, bool EnergyLossOn) ;
    bool is_initialised ;
    
    //** find the measurment layer for a given hit 
    const ILDVMeasLayer* findMeasLayer( EVENT::TrackerHit * trkhit) const ; 
    //** find the measurment layer for a given det element ID and point in space 
    const ILDVMeasLayer* findMeasLayer( int detElementID, const TVector3& point) const ;
    
    // get the last layer crossed by the helix when extrapolating from the present position to the pca to point
    const ILDVMeasLayer* getLastMeasLayer(THelicalTrack const& helix, TVector3 const& point) const ;
    
    const ILDCylinderMeasLayer* getIPLayer() const { return _ipLayer; }
    
    const ILDCylinderMeasLayer* _ipLayer ;
    
    const gear::GearMgr* _gearMgr ;
    
    TKalDetCradle* _det ;            // the detector cradle
    
    std::multimap< int,const ILDVMeasLayer *> _active_measurement_modules;
    
    std::multimap< int,const ILDVMeasLayer *> _active_measurement_modules_by_layer;
    
    
#ifdef MARLINTRK_DIAGNOSTICS_ON

  private:    
    MarlinTrk::DiagnosticsController _diagnostics;

  public:    

    virtual void * getDiagnositicsPointer() { return &_diagnostics ; }
            
#endif

    
  } ;
  
} // end of namespace MarlinTrk

#endif