---

TkrClustersRep.cxx

Go to the documentation of this file.

#include "TkrRecon/TkrClustersRep.h"


//------------------------------------------------------------------------------
//  Code to implement the TkrClustersRep class. 
//  This class does the display of clusters in the Silicon
//  Tracker. 
//
//#############################################################################
//  Constructor for the class
TkrClustersRep::TkrClustersRep(SiClusters** ppClus)
//#############################################################################
{
        //Store pointer to the pointer to the SiClusters data
        ppClusters = ppClus;
}

//##############################################
//  Display routine for the class
void TkrClustersRep::update()
//##############################################
{
        //Recover pointer to the data
    SiClusters* pClusters = *ppClusters;

    //Zero out the pointer so we don't accidentally try to redraw the event
    *ppClusters = 0;

        //Make sure pointer is valid before trying to draw anything
        if (pClusters)
        {
                int    nHits      = pClusters->nHits();
                double stripPitch = pClusters->stripPitch();
                double towerPitch  = pClusters->towerPitch();
                
                setColor("green");

                //Loop over all cluster hits in the SiClusters vector
                while(nHits--)
                {
                        SiCluster* pCluster = pClusters->getHit(nHits);
                        Point      clusPos  = pCluster->position();
        
                        double     x        = clusPos.x();
                        double     y        = clusPos.y();
                        double     z        = clusPos.z();
        
                        double     delta    =  0.6;
                        double     Offset   = -0.5 * towerPitch;

                        //Draw a cross at the position of the cluster center
                        if (pCluster->v() == SiCluster::X)
                        {
                Offset += y;
                moveTo(Point(x - delta, Offset, z - delta));
                                lineTo(Point(x + delta, Offset, z + delta));
                                moveTo(Point(x - delta, Offset, z + delta));
                                lineTo(Point(x + delta, Offset, z - delta));
                        }
                        else
                        {
                Offset += x;
                                moveTo(Point(Offset, y - delta, z - delta));
                                lineTo(Point(Offset, y + delta, z + delta));
                                moveTo(Point(Offset, y - delta, z + delta));
                                lineTo(Point(Offset, y + delta, z - delta));
                        }

                        //Now draw the width of the cluster (if more than one)
                        if (pCluster->size() > 1)
                        {
                                double stripWid = 0.5 * pCluster->size() * stripPitch;

                                if (pCluster->v() == SiCluster::X)
                                {
                                        moveTo(Point(x - stripWid, Offset, z));
                                        lineTo(Point(x + stripWid, Offset, z));
                                }
                                else
                                {
                                        moveTo(Point(Offset, y - stripWid, z));
                                        lineTo(Point(Offset, y + stripWid, z));
                                }
                        }
                }
        }

    return;
}

---