---

RecoCsIData.cxx

Go to the documentation of this file.

// $Header: /nfs/slac/g/glast/ground/cvs/reconstruction/src/data/RecoCsIData.cxx,v 1.2 2000/12/11 16:38:45 burnett Exp $


// Implementation of the RecoCsIData class for testing reconstruction

#include "reconstruction/data/RecoCsIData.h"
#include "instrument/Calorimeter.h"


RecoCsIData::RecoCsIData (int numLayers)
{
    for(int i=0; i<numLayers; i++)  {    
        calorList.push_back((new std::vector<Xtal>));
    }
}

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

int RecoCsIData::nHits (unsigned int layer) const
{
    return calorList[layer]->size();
}

float RecoCsIData::energy (unsigned int layer, unsigned int n) const
{
    return (*calorList[layer])[n].energy;
}

Point RecoCsIData::xtalPos (unsigned int layer, unsigned int n) const
{
    return  (*calorList[layer])[n].pos;
}


idents::ModuleId RecoCsIData::moduleId (unsigned int layer, unsigned int n) const
{
    return (*calorList[layer])[n].module;
}

idents::XtalId RecoCsIData::xtalId (unsigned int layer, unsigned int n) const
{
    return (*calorList[layer])[n].id;
}

float RecoCsIData::Lresp (unsigned int layer, unsigned int n) const
{
    return (*calorList[layer])[n].Lresp;
}

float RecoCsIData::Rresp (unsigned int layer, unsigned int n) const
{
    return (*calorList[layer])[n].Rresp;
}

const std::vector<double>& RecoCsIData::Diodes_Energy (unsigned int layer, unsigned int n) const
{
    return (*calorList[layer])[n].Diodes_Energy;
}

void RecoCsIData::load (const CsIDetector& xtal, idents::ModuleId moduleId)
{
    if (xtal.hit()) {
        calorList[xtal.layer()]->push_back(Xtal(xtal.xtalPos(), xtal.energy(),moduleId ,
            xtal.xtalId(), xtal.Lresp(), xtal.Rresp(),
                                xtal.getDiodesEnergy()));

        // put the "raw" information into the basic list
        idents::XtalId id( moduleId, xtal.layer(), xtal.xtalId());
        int     left =  static_cast<int>(xtal.Lresp()*1e6),
            right = static_cast<int>(xtal.Rresp()*1e6);
        m_xtals[id]=std::pair<int,int>(left,right);
    }
}




void RecoCsIData::readData (std::istream& in)
{
   float Left, Right;
   int X, Y, Z, E;

   int numLayers;
   in>> numLayers;

   for(int i=0; i<numLayers; i++) {
       int numX;
       in>>numX;
       unsigned st, mod;
       for(int j=0; j<numX; j++) {
         in>>mod>>st>>E>>X>>Y>>Z>>Left>>Right;
         calorList[i]->push_back(Xtal(Point(X/1e3, Y/1e3, Z/1e3), E/1e6, idents::ModuleId(mod), st,
           Left/1e6, Right/1e6));
       }
   }
}

void RecoCsIData::writeData (std::ostream& out) const
{
   int numLayers = calorList.size();

   out<<numLayers<<'\n';
   if(out.eof() ) return;  // happens to be the first
   for(int i=0; i<numLayers; i++) {
       int numX = nHits(i);
       out<<numX<<'\n';
       for(int j=0; j<numX; j++) {
            out<<moduleId(i, j)<<' '<<xtalId(i, j)<<' '
               <<int(1e6*energy(i,j))       <<' '
               <<int(1e3*xtalPos(i, j).x())<<' '
               <<int(1e3*xtalPos(i, j).y())<<' '
               <<int(1e3*xtalPos(i, j).z())<<' '
                           <<int(1e6*Lresp(i, j))<<' '
                           <<int(1e6*Rresp(i, j))<<'\n';
       }
   }
}

void RecoCsIData::clear ()
{
    m_xtals.clear();
    for(unsigned i=0; i<calorList.size(); i++) {
        calorList[i]->clear();
    }
}

void RecoCsIData::printOn(std::ostream& out) const
{
    out << "RecoCsIData" << std::endl;
    writeData(out);
}

---