---

TowerBoundaries.cxx

Go to the documentation of this file.

// $Id: TowerBoundaries.cxx,v 1.1.1.1 1999/12/28 01:36:03 burnett Exp $
// 
//  Original author: Sawyer Gillespie
//                   hgillesp@u.washington.edu
//

#include "reconstruction/TowerBoundaries.h"    // include class declaration
#include "reconstruction/GlastFit.h"
#include "reconstruction/GlastRecon.h"
#include "reconstruction/LbldData.h"
#include "reconstruction/TrackerRecon.h"
#include "instrument/Glast.h"

// implementation of the TowerBoundaries class

TowerBoundaries::TowerBoundaries (GlastRecon* r)
: ReconAnalysis (r)
{
}

TowerBoundaries::~TowerBoundaries ()
{
}

void    TowerBoundaries::setUpData ()
{
    i_Tower_Bnd = data()->addElement("Conv_Twr_Dist");
    i_Fit_xv = data()->addElement("Fit_xV");
    i_Fit_yv = data()->addElement("Fit_yV");
}

void    TowerBoundaries::compute ()
{
    // get the TrackerRecon object which has reconstructed the
    // data within the silicon tracker
    TrackerRecon*   tr = recon()->getTracker();
    
    // declare local variables

    // get the gamma-ray fit from TrackerRecon
    GFgamma*    gamma = tr->gammaFit ();

    // get the instrument parameters from TrackerRecon
    const TrackerRecon::Parameters& s_parms = tr->params ();

    // get the vertex of the gamma conversion?
    Point vertex = gamma->getFirstHit();
    float x_ist = vertex.x();
    float y_ist = vertex.y();

    // compute the results    
    float num_2 = Glast::xNum/2.; 
    float x_twr_bnd = fmod((x_ist + s_parms.tower_width*num_2), s_parms.tower_width);
    x_twr_bnd = x_twr_bnd - s_parms.tower_width*.5;
    float y_twr_bnd = fmod((y_ist + s_parms.tower_width*num_2), s_parms.tower_width);
    y_twr_bnd = y_twr_bnd - s_parms.tower_width*.5;
    
    x_twr_bnd = fabs(x_twr_bnd);
    y_twr_bnd = fabs(y_twr_bnd);
    float twr_bnd   = (x_twr_bnd > y_twr_bnd) ? x_twr_bnd : y_twr_bnd;
    
    // set the state of this object to reflect the results of the
    // computation
    ReconAnalysis::result (new Result(twr_bnd, x_twr_bnd, y_twr_bnd));

    // add the data to the list of archival data elements.
    tr->data.addData(i_Tower_Bnd, twr_bnd);
    tr->data.addData(i_Fit_xv, x_twr_bnd);
    tr->data.addData(i_Fit_yv, y_twr_bnd);
}

---