---

ActiveDistance.cxx

Go to the documentation of this file.

// $Id: ActiveDistance.cxx,v 1.6 2000/12/11 16:38:44 burnett Exp $
// 
//  Original author: Sawyer Gillespie
//                   hgillesp@u.washington.edu
//

#include "reconstruction/ActiveDistance.h"    // include class declaration

#include "reconstruction/GlastFit.h"
#include "reconstruction/GlastRecon.h"
#include "reconstruction/LbldData.h"
#include "reconstruction/TrackerRecon.h"

// implementation of the ActiveDistance class

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

ActiveDistance::~ActiveDistance ()
{
}

void    ActiveDistance::setUpData ()
{
    i_activeDist = data()->addElement("Active_Dist");
}

void    ActiveDistance::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 vertex of the gamma conversion & direction of gamma
    Vector  t0 = gamma->direction();
    Point   x0 = gamma->vertex();

    // native code - this code has been brought in from TrackerRecon
    // -------------------------------------------------------------
    float activeDist = -20.;
    Point firstHit = gamma->getFirstHit();
    Ray prjRay(firstHit, -t0);
    Point x_prj;
    
#if 0 // need to figure out how to  do this
    GlastMed *glast = GlastMed::instance();
    
    int istLayer = gamma->firstLayer();
    
    float firstProp = 0.;
    if (istLayer > 0) { //go to next gap up + espsilon
        firstProp = SiTracker::traySpacing()/fabs(t0.z());
        x_prj = prjRay.position(firstProp);
    }
    else {
        x_prj = firstHit;
        firstProp = 0.;
    }
    
    const Medium *med0 = glast->inside(x_prj);
    prjRay = Ray(x_prj, -t0);
    const Medium *detMed = 0;
    if(med0){
        firstProp = med0->distanceToLeave(prjRay, detMed, 100.);
        activeDist = -2. - firstProp;
    }
    
    Detector *det =  0;
    if(detMed) {
        med0 = detMed;
        x_prj = prjRay.position(firstProp);
        prjRay = Ray(x_prj, -t0);
        firstProp = med0->distanceToLeave(prjRay, detMed, 100.);
        activeDist = -4. - firstProp;
        if(detMed) det = detMed->detector();
        if(det) {
            if(!strcmp(det->nameOf(), "SiDetector")) {
                x_prj = prjRay.position(firstProp);
                activeDist = static_cast<const SiDetectorMed*>(det)->data().insideActiveArea(x_prj);
            }
        }
    }
    //-----------------------------------------------------------
    // end native code
  #endif    
  
    // set the state of this object to reflect the results of the
    // computation
    ReconAnalysis::result (new Result(activeDist));

    // add the data to the list of archival data elements.
    tr->data.addData(i_activeDist, activeDist);
}

---