---

TrackerRecon.h

Go to the documentation of this file.

//-------------------------------------------------------------------
//
//     TrackerRecon:
//
//          Steers the Silicon-Tracker Reconstruction    
//
//                    Bill Atwood
//                    B. Atwood, JA Hernando, Santa Cruz, 02/05/99
//
//-------------------------------------------------------------------

#ifndef __TRACKERRECON_H
#define __TRACKERRECON_H 1

#include "reconstruction/Recon.h"
#include "data/SiData.h"
#include "reconstruction/CalRecon.h"
#include "reconstruction/MCRecon.h"
#include "geometry/Ray.h"
#include "reconstruction/GlastFit.h"
#include "reconstruction/KalFit.h"

// forward declarations
class Point;
class Vector;
//class CalRecon;
// class GlastFit;
class SiClusters;
class ReconVisitor;
class GlastFit;
class GlastFit;
class GlastRecon;
class GFpair;
class GFgamma; 
class GFbase;
class ActiveDistance;
class EnergyCorrection;
class ExtraHits;
class TowerBoundaries;
class TrackerVeto;

class TrackerRecon : public Recon
{
public:

    TrackerRecon(); // test constructor, tracker only
    TrackerRecon (GlastRecon* c);    
    ~TrackerRecon ();
    void initialize();
   
        // access functions
    Point center () {return(x1);}
    Vector direction () {return(t1);}
    const inline CalRecon* getCalRecon() {return cal;}
        
        // Kalman Track list
    std::vector<KalTrack> xtrackList;
    std::vector<KalTrack> ytrackList;

        // operations
    void accept (ReconVisitor& d);
    void clear ();
    void draw (gui::DisplayRep& v);
    // Do the Tracker Reconstruction
    void reconstruct(const SiData* stripData);

    // access methods

    // access the silicon cluster finding data
    SiClusters* getClusters ();

    // access the best CsICorrEnergy estimate
    double  getCsICorrEnergy () const;

    // provide access to the GlastFit gamma reconstruction tracks
    GFgamma*    gammaFit ();

    // private structure - defines instrument parameters
    // HSG: 7/10/99
    // this is sketchy, but it won't change all of the calls
    // in TrackerRecon.cxx...(change it later)
    class   Parameters {
    public:
        int num_convs;          // number of converters
        float conv_gap;         // converter gap ?
        float conv_rad_len;     // radiation length of converters
        float si_thickness;     // thickness of the silicon
        float tower_width;      // width of a tower module

    protected:
        // initialize all of the parameters
        void    initialize ();

        // make TrackerRecon a friend so that only TrackerRecon
        // can do the initialization
        friend class TrackerRecon;
    };

    // provide access to the standard tracker parameters
    static const   Parameters& params ();

    // need for the Csi_corr_energy
    enum configuration {BASELINE, SUPERGLAST}; // need for the Csi_corr_energy

private:
    //  Original selection of the best/pair GlastFit
    //----------------------------------------------
    GlastFit* bestTrkHypo (enum SiData::Axis axis, float numSigma,
                             float energy, int firstLayer);
    
    GlastFit* calTrkHypo  (enum SiData::Axis axis, float numSigma,
                             float energy, int firstLayer, const Point& calHit);
    
    GlastFit* pairTrkHypo (enum SiData::Axis axis, float numSigma,
                             GlastFit& bestFit);
    
    GlastFit* nextTrkHypo (enum SiData::Axis axis, float numSigma,
                             float energy, int firstLayer);

    // Kalman Filter Private Ntuple for debuging purposes (JAH 011399)
    // void KalPrivateTup();
    
    // Tracker Analysis functions
    //----------------------------
    double TAna_trackerVeto(double& xsigma, double& ysigma);
    double TAna_kink(double& pairKink,double& fitKinkN, double& pairKinkN);
    double TAna_densityHits(double& previousHits, double& nextHits);
    double TAna_towerBoundary(double& xboundary, double& yboundary);
    double TAna_activeDistance();
    double TAna_extraHits(double& nlayers, double& istHitCount, double& outHits, double& shwHits1, double& shwHits2);
    double TAna_corrEnergy(double shwHits1, double shwHits2, double boundary);
    
    // Original TrackerRecon
    //--------------------------------------
    void iniTrackerRecon();
    void iniClusters(const SiData*);
    void iniCalorimeterInput();
    
    
    void moreGlastTracks();
    void loadResultsIntoData(const SiData*);

    // Update TrackerRecon
    //-----------------------------------
    void RCreconstruct(const SiData*);
    
    bool seedCandidates(GFdata::type type);
    bool candidates(GFdata::type , GFdata::type seeds);
    void selectGamma(); 
    void selectParticle();
    
    bool findSeedCandidates(std::vector<GFdata>& candidates, GFdata::type type, SiData::Axis); 
    bool findSeedCandidates(std::vector<GFdata>& candidates, GFdata::type type, SiData::Axis, int iplane, int itower = 0);
    bool findCandidates(std::vector<GFdata>& candidates, const GFdata& Xcandidate, const GFdata& Ycandidate, GFdata::type type);
    GFdata GFconstructor(GFdata::type type, int ilayer, const Ray testRay, SiData::Axis axis = SiData::X);
    void incorporate(std::vector<GFdata>& candidates, const GFdata);
    
    double corrEnergy(int ilevel, enum configuration, int nShowerHits, double slope);
    Point createPRef(SiData::Axis axis, int ilayer, const Point& PIni);
    
    void loadTrackerRecon();
    void loadTuple(const SiData*);
        
    std::vector<GFdata> m_Xcandidates;
    std::vector<GFdata> m_Ycandidates;
    
    std::vector<GFdata> m_candidates;
    
    GFgamma* _mGFgamma;
    

    // Silicon Hits
    LbldData::iterator i_Conv_Layer_Hits;
    LbldData::iterator i_max_controller_hits;
    LbldData::iterator i_First_Conv_Layer;
    LbldData::iterator i_nConv_Layers_Hit;
    LbldData::iterator i_Conv_Layer_Hits_Gamma;
    LbldData::iterator i_nConv_Layers_Gamma;
    LbldData::iterator i_First_Conv_Layer_Gamma;
    
    // Tracker Reconstruction
    LbldData::iterator i_Kal_Energy;
    LbldData::iterator i_Kal_Energy_X;
    LbldData::iterator i_Kal_Energy_Y;
    
    // Event Reconstruction
    LbldData::iterator i_No_XTrks;
    LbldData::iterator i_No_YTrks;
    LbldData::iterator i_No_Trks;
    
    // Gamma quality
    LbldData::iterator i_fit_type;
    LbldData::iterator i_fit_topo;
    LbldData::iterator i_Chisq_1st;
    LbldData::iterator i_Chisq;
    LbldData::iterator i_No_Hits;
    LbldData::iterator i_No_SHits;
    LbldData::iterator i_No_Gaps;
    LbldData::iterator i_No_Gaps_1st;
    LbldData::iterator i_No_Noise;
    LbldData::iterator i_No_Noise_1st;
    LbldData::iterator i_qual_X;
    LbldData::iterator i_qual_Y;
    LbldData::iterator i_qual;
    
    // Pair reconstruction
    LbldData::iterator i_Fit_XNhits;
    LbldData::iterator i_Fit_YNhits;
    LbldData::iterator i_Fit_XChisq;
    LbldData::iterator i_Fit_YChisq;
    LbldData::iterator i_Fit_XChisq_1st;
    LbldData::iterator i_Fit_YChisq_1st;
    LbldData::iterator i_Fit_XKalThetaMS;
    LbldData::iterator i_Fit_YKalThetaMS;
    LbldData::iterator i_Fit_XKalEne;
    LbldData::iterator i_Fit_YKalEne;
    LbldData::iterator i_Fit_x0;
    LbldData::iterator i_Fit_y0;
    LbldData::iterator i_Fit_z0;
    LbldData::iterator i_Fit_xdir;
    LbldData::iterator i_Fit_ydir;
    LbldData::iterator i_Fit_zdir;
    
    LbldData::iterator i_Pair_XNhits;
    LbldData::iterator i_Pair_YNhits;
    LbldData::iterator i_Pair_XChisq;
    LbldData::iterator i_Pair_YChisq;
    LbldData::iterator i_Pair_XChisq_1st;
    LbldData::iterator i_Pair_YChisq_1st;
    LbldData::iterator i_Pair_XKalThetaMS;
    LbldData::iterator i_Pair_YKalThetaMS;
    LbldData::iterator i_Pair_XKalEne;
    LbldData::iterator i_Pair_YKalEne;
    LbldData::iterator i_Pair_x0;
    LbldData::iterator i_Pair_y0;
    LbldData::iterator i_Pair_z0;
    LbldData::iterator i_Pair_xdir;
    LbldData::iterator i_Pair_ydir;
    LbldData::iterator i_Pair_zdir;
    
    LbldData::iterator i_e_frac;
    LbldData::iterator i_errSlopeX;
    LbldData::iterator i_errSlopeY;
    LbldData::iterator i_weightXSlope;
    LbldData::iterator i_weightYSlope; 
    LbldData::iterator i_xeneXSlope;
    LbldData::iterator i_xeneYSlope;    
    
    // gamma reconstruction
    LbldData::iterator i_First_XHit;
    LbldData::iterator i_First_YHit;
    LbldData::iterator i_Diff_1st_Hit;
    LbldData::iterator i_Gamma_x0;
    LbldData::iterator i_Gamma_y0;
    LbldData::iterator i_Gamma_z0;
    LbldData::iterator i_Gamma_xdir;
    LbldData::iterator i_Gamma_ydir;
    LbldData::iterator i_Gamma_zdir;
    
    // analysis
    LbldData::iterator i_CsI_corr_Energy;
    LbldData::iterator i_1st_count;
    LbldData::iterator i_srp_hit_rat;
    LbldData::iterator i_out_hit_rat;
    LbldData::iterator i_activeDist;
    LbldData::iterator i_showerHits1;
    LbldData::iterator i_showerHits2;

    LbldData::iterator i_Gamma_DLT;
    LbldData::iterator i_Tower_Bnd;
    LbldData::iterator i_Fit_xv;
    LbldData::iterator i_Fit_yv;
    LbldData::iterator i_Zdiff_XY;
    LbldData::iterator i_Zdiff_plane;
    LbldData::iterator i_trkVeto_code;
    LbldData::iterator i_trkVeto_Xsigma;
    LbldData::iterator i_trkVeto_Ysigma;

    LbldData::iterator i_fit_kink;
    LbldData::iterator i_pair_kink;
    LbldData::iterator i_fit_kinkN;
    LbldData::iterator i_pair_kinkN;
    
    LbldData::iterator i_hits_previous;
    LbldData::iterator i_hits_first;
    LbldData::iterator i_hits_next;

    // HMA (GSFC) silicon clusters
    LbldData::iterator i_FitSide_nXused, i_FitSide_nYused;

    //  Common Data Members 
    Point x1;
    Vector t1;
    Vector t0;
    Vector t2;
    Point x0;
    Point x2;
    
    //  Si variables
    SiClusters* _mclsData;
    int m_totalHits;
    int m_firstConvLayer;
    int m_numLayersHit;
    
    //  Cal variables
    const CalRecon* cal;
    double m_CsIEnergy;
    double m_CsICorrEnergy;
    Point m_aveCalHit;
    
    // Internal Additional members 
    int m_firstLayer;
    int m_fitType;
    double m_xFactor;
    double m_yFactor; 
    
    // HMA (GSFC) variables used for clusters of silicon hits
    Point clusterFitPos;
    Vector clusterFitDir;
    int nHitsUsedCluster;
    
    // Analysis Objects 
    // ------------------------------------------------------------------------------
    ActiveDistance*     m_activeDist;   // object which handles Active_Dist comp.
    EnergyCorrection*   m_Ecorr;        // object which handles the energy correction
    ExtraHits*          m_extraHits;    // object which examines hit patterns
    TowerBoundaries*    m_twrBounds;    // object which examines twr boundaries
    TrackerVeto*        m_tkrVeto;      // object which examines tkr veto handling
    // ------------------------------------------------------------------------------

    // instrument-wide parameters
    static  Parameters  s_parms;
};

// inline declarations

// TrackerRecon::Parameters class

inline  const TrackerRecon::Parameters& TrackerRecon::params ()
{
    return s_parms;
}

// TrackerRecon class

inline  double  TrackerRecon::getCsICorrEnergy () const
{
    return m_CsICorrEnergy;
}

inline  SiClusters* TrackerRecon::getClusters ()
{
    return _mclsData;
}

#endif

---