---

TrackerVeto.cxx

Go to the documentation of this file.

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

#include "reconstruction/TrackerVeto.h"    // include class declaration
#include "reconstruction/GlastFit.h"
#include "reconstruction/GlastRecon.h"
#include "reconstruction/LbldData.h"
#include "reconstruction/TrackerRecon.h"

// implementation of the TrackerVeto class

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

TrackerVeto::~TrackerVeto ()
{
}

void    TrackerVeto::setUpData ()
{
    i_trkVeto_code = data()->addElement("trkVeto_code");
    i_trkVeto_Xsigma = data()->addElement("trkVeto_Xsigma");
    i_trkVeto_Ysigma = data()->addElement("trkVeto_Ysigma");
}

void    TrackerVeto::compute ()
{
    // get the TrackerRecon object which has reconstructed the
    // data within the silicon tracker
    TrackerRecon*   tr = recon()->getTracker();
    
    // find the best x and y fits to the gamma ray
    GFtrack* _Xbest = tr->gammaFit()->getXpair()->getBest();
    GFtrack* _Ybest = tr->gammaFit()->getYpair()->getBest();

    // local variables
    int indexhit = 0;
    double  xsigma = 0, ysigma = 0;

    // determine incidence w. ACD system
    bool    Xveto = _Xbest->veto(indexhit, xsigma);
    bool    Yveto = _Ybest->veto(indexhit, ysigma);

    // determine the id code based upon the results
    float   idcode = 0;
    if (Xveto)  idcode = 1;
    if (Yveto)  idcode += 2;

    // set the state of this object to reflect the results of the
    // computation
    ReconAnalysis::result (new Result(idcode, xsigma, ysigma));

    // add the data to the list of archival data elements.
    tr->data.addData(i_trkVeto_code, idcode);
    tr->data.addData(i_trkVeto_Xsigma, xsigma);
    tr->data.addData(i_trkVeto_Ysigma, ysigma);
}

---