---

RecoSiData.cxx

Go to the documentation of this file.

// $Id: RecoSiData.cxx,v 1.2 2000/12/11 16:38:45 burnett Exp $
//
// Implement access to Silicon strip data

#include "reconstruction/data/RecoSiData.h"

#include "instrument/SiTracker.h"
#include <algorithm>
using namespace std;

class RecoSiData::Strip {
private:
    friend class RecoSiData;
    
public:
    Strip () {}
    
private:
    Strip (Point p, Id m, unsigned int n, unsigned int t = 0)
        : stripIndex(n), stripType(t), pos(p), module(m)
    {}
    
    unsigned int stripIndex;
    unsigned int stripType; 
    Point pos;
    Id module;
};


RecoSiData::RecoSiData ()
{
    for(int i = 0; i < SiTracker::numTrays(); i++) {
        xhitList.push_back(new vector<Strip>);
        yhitList.push_back(new vector<Strip>);
    }
    m_total_hits=0;
}

RecoSiData::RecoSiData (unsigned int n)
{
    for(unsigned i=0; i<n; i++) {       
        xhitList.push_back(new vector<Strip>);
        yhitList.push_back(new vector<Strip>);
    }
}


RecoSiData::~RecoSiData ()
{
    clear();
    for(unsigned i=0; i<xhitList.size(); i++) {
        delete xhitList[i];
        delete yhitList[i];
    }
    
}


const SiData::Id& RecoSiData::moduleId (enum RecoSiData::Axis a, 
                                  unsigned int tray, 
                                  unsigned int n      ) const
{
    return a==X? (*xhitList[tray])[n].module
        : (*yhitList[tray])[n].module;
}

void RecoSiData::load (const SiDetector& plane, idents::ModuleId moduleId)
{
    static double tray_max=48.10, tray_spacing=3.20;
    
    CoordTransform T 
        = plane.localToGlobal(); // for converting the coords following
    
    // figure out which tray this is from the transform
    Vector a(1,0,0),b(1,0,0); a.transform(T);
    double rot = a*b; 
    Axis axis = (rot>0.5)? X : Y;
    Point p(0,0,0); p.transform(T); 
    double z = p.z();
    int tray = static_cast<int>((tray_max-z)/tray_spacing);
    tray = tray< 0?  0: tray;
    tray = tray>17? 17: tray;
    
    // final loop over hits
    for( SiDetector::const_iterator it3=plane.begin(); 
    it3!=plane.end(); ++it3)                      {
        int wire_num = (*it3).index();
        int wire_typ = (*it3).noise();
        //************** remove noise  *********************
        if( wire_typ !=0 ) continue; 
        // *************************************************
        Point p = Point(SiDetector::localX(wire_num),0,0);
        p.transform(T);
        if( axis == X) { 
            xhitList[tray]->push_back(
                Strip(p, moduleId, wire_num, wire_typ));
        } else {
            yhitList[tray]->push_back(
                Strip(p, moduleId, wire_num, wire_typ));
        }
        m_total_hits++; 
    }
}

void RecoSiData::clear () {
    for(unsigned i=0; i<xhitList.size(); i++) {
        xhitList[i]->clear();
        yhitList[i]->clear();
    }
    m_total_hits=0;
}

int RecoSiData::nHits (enum RecoSiData::Axis a, int tray) const
{
    return a==X?  xhitList[tray]->size()
        :  yhitList[tray]->size() ;
    
}

Point RecoSiData::hit (enum RecoSiData::Axis a, 
                   unsigned int tray, 
                   unsigned int n       ) const
{
    return a==X? (*xhitList[tray])[n].pos
        : (*yhitList[tray])[n].pos;
}

unsigned int RecoSiData::hitId (enum RecoSiData::Axis a, 
                            unsigned int tray, 
                            unsigned int n       ) const
{
    return a==X? (*xhitList[tray])[n].stripIndex
        : (*yhitList[tray])[n].stripIndex;
}
#if 0
unsigned int RecoSiData::hitType (enum RecoSiData::Axis a, 
                              unsigned int tray, 
                              unsigned int n       ) const
{
    return a==X? (*xhitList[tray])[n].stripType
        : (*yhitList[tray])[n].stripType;
}
#endif
int RecoSiData::totalHits () const
{
    return m_total_hits;
}

void RecoSiData::readData (istream& in)
{
    int x, y, Z;
    
    int numLayers;
    in>> numLayers;
    if( in.eof() )  // this is the first one
        return;
    
    for(int i=0; i<numLayers; i++) {
        int numX;
        in>>numX;
        unsigned int mod;
        unsigned st, type;
        int j;
        for( j=0; j<numX; j++) {
            in>>mod>>st>>type>>x>>y>>Z;
            xhitList[i]->push_back(
                Strip(Point(x/1e3, y/1e3, Z/1e3), idents::ModuleId(mod), st, type));
        }
        int numY;
        in>>numY;
        for(j=0; j<numY; j++) {
            in>>mod>>st>>type>>x>>y>>Z;
            yhitList[i]->push_back(
                Strip(Point(x/1e3, y/1e3, Z/1e3), idents::ModuleId(mod), st, type));
        }
    }
}

void RecoSiData::writeData (ostream& out)
{
    int numLayers = xhitList.size();
    
    out<<numLayers<<'\n';
    
    for(int i=0; i<numLayers; i++) {
        int numX = nHits(X, i);
        out<<numX<<'\n';
        int j;
        for(j=0; j<numX; j++) {
            out<<moduleId(X, i, j)
                <<' '<<hitId(X, i, j) <<' '
//                <<hitType(X, i, j)<<' '
                <<int(1e3*hit(X, i, j).x())<<' '
                <<int(1e3*hit(X, i, j).y())<<' '
                <<int(1e3*hit(X, i, j).z())<<'\n';
        }
        int numY = nHits(Y, i);
        out<<numY<<'\n';
        for(j=0; j<numY; j++) {
            out<<moduleId(Y, i, j) 
                <<' '<<hitId(Y, i, j)<<' '
//                <<hitType(Y, i, j)<<' '
                <<int(1e3*hit(Y, i, j).x())<<' '
                <<int(1e3*hit(Y, i, j).y())<<' '
                <<int(1e3*hit(Y, i, j).z())<<'\n';
        }
    }
}

void RecoSiData::printOn (ostream& cout) const
{
    unsigned tray, i;
    cout << "\nRecoSiData:\n";
    for(tray=0; tray < xhitList.size() && tray < yhitList.size(); tray++)
    {
        unsigned nx= nHits(X,tray),
            ny= nHits(Y,tray);
        if( nx==0 && ny==0 ) continue;
        cout <<  tray << "X";
        for( i =0; i<nx; i++)
            cout << '\t'<< hit(X,tray,i) 
//            <<"\t " << hitType(X,tray,i) 
            << '\n';
        cout <<  tray << "Y";
        for( i =0; i<ny; i++)
            cout << '\t'<< hit(Y,tray,i)  
//            <<"\t " << hitType(Y,tray,i) 
            << '\n';
    }
    return;
}

// Class RecoSiData_Hit

---