---

CalRecLogs.cpp

Go to the documentation of this file.

#include "CalRecon/CalRecLogs.h"
#define LOGWIDTH 3.05

//----------------- CalRecLog ------------------
//################################################
CalRecLog::CalRecLog(int ilayer, CalDetGeo::axis v, int icolumn,idents::ModuleId mod) : 
CalADCLog(ilayer,v,icolumn,mod)
//################################################
{
        clear();
}

//################################################
CalRecLog::CalRecLog(int ilayer, CalDetGeo::axis v, int icolumn) : 
CalADCLog(ilayer,v,icolumn)
//################################################
{
        clear();
}
//################################################
double CalRecLog::asymmetry() const
//################################################
{
        double asymmetry = 0.;
        if (m_pEnergy > 0 && m_nEnergy > 0)
                asymmetry = (m_pEnergy-m_nEnergy)/(m_pEnergy+m_nEnergy);
        return asymmetry;
}
//################################################
void CalRecLog::writeOut() const
//################################################
{
        if(energy()>0){
                std::cout << " ID " << logID() <<  " - "<< modId()<<" " << layer()<<" "<< view()
                                  <<" " << column() << " ";
                std::cout << " E-+ " << energy() << " " << negEnergy() << " " << posEnergy() ;
                std::cout << " Asy " << asymmetry();
                std::cout << " Pos " << position().x() << " " << position().y() << " " << position().z();
                std::cout << "\n";
        }
}
//######################################################
void CalRecLog::draw(gui::DisplayRep& v) const
//######################################################
{
        double MINENE = 10.;
        double FACTOR = 0.1;
//      double delta = 0.2*calorimeterGeo::logWidth();
        double delta = 0.2*LOGWIDTH;
        double x = position().x();
        double y = position().y();
        double z = position().z();
        double ene = energy();
        if (ene <= MINENE) return;

        v.setColor("black");
        v.moveTo(Point(x-delta,y,z));
        v.lineTo(Point(x+delta,y,z));
        v.moveTo(Point(x,y-delta,z));
        v.lineTo(Point(x,y+delta,z));

        v.moveTo(Point(x,y,z));
        v.lineTo(Point(x,y,z-log(FACTOR*ene)));
}
//------------- private --------------------------
//################################################
void CalRecLog::clear()
//################################################
{
        CalADCLog::clear();

        for (int irange = 0; irange < CALNRANGES; irange++) {
                m_negEnergy[irange] = 0;
                m_posEnergy[irange] = 0;
        }

        m_nEnergy = 0.;
        m_pEnergy = 0.;
        m_position = Point(0.,0.,0.);
}

//------------ private ---------------------------
//######################################
void CalRecLogs::ini(int nModX, int nModY, int nLogs, int nLayers, int nViews)
//######################################
{
        for (int iy = 1; iy <= nModY; iy++){            
                for (int ix = 1; ix <= nModX; ix++){
                        idents::ModuleId mod(ix,iy);
                        
                        for (int ilayer = 0; ilayer < nLayers; ilayer++) {
                                for (int v=0; v < nViews; v++){
                                        CalDetGeo::axis view = CalDetGeo::makeAxis(v);
                                        for (int ilog = 0; ilog < nLogs; ilog++) {
                                                m_List.push_back(new CalRecLog(ilayer,view,ilog,mod));
                                        }
                                }
                        }
                }
        }
}
//######################################
CalRecLogs::~CalRecLogs()
//######################################
{
    if(m_calDisp) m_calDisp->clearLogsDisp();
        for (int ilog = 0; ilog < m_List.size(); ilog++) delete m_List[ilog];
}
//######################################
void CalRecLogs::clear()
//######################################
{
        for (int ilog = 0; ilog < m_List.size(); ilog++) {
                m_List[ilog]->clear();
        }
}
//######################################
CalRecLog* CalRecLogs::getLogID(int logID) const
//######################################
{
        int jlog = -1;
        for (int ilog = 0; ilog < m_List.size(); ilog++) {
                if (m_List[ilog]->logID() == logID) {
                        jlog  = ilog;
                        break;
                }
        }
        return m_List[jlog];
}

//######################################
void CalRecLogs::writeOut() const
//######################################
{
        if (m_List.size()<=0) return;

        std::cout << " CalRecLogs " << m_List.size() <<"\n";
        for (int i = 0; i < m_List.size();i++) {
                m_List[i]->writeOut();
        }
}
//######################################
void CalRecLogs::draw(gui::DisplayRep& v) const
//######################################
{
        if (m_List.size() <=0) return;
        for (int i = 0; i < m_List.size();i++) m_List[i]->draw(v);
}

---