HadMon.cxx

Go to the documentation of this file.

#include "HadMon.h"

#include <cmath>

const int HadMon::NROWCOLS = 7;
const double HadMon::ROWCOL_SPACING = 11.43; // cm

int HadMon::Index(int channel)
{
    // hard code channel map.
    // Channel numbers are 1-1 with buffer indices except:
    // - channels start with 1, indices start with 0
    // - buffer starts at index 104
    // - vert/horiz blocks are reversed.  Ie:
    //   channel = (01-48)(49-96)
    //   buffer  = (48-95)(00-47)

    if (channel < 1 || channel > 49) return -1;
    int offset = -1;
    if (channel == 49) offset = 104;
    else offset = channel + 151;
    return offset - 104;
}

int HadMon::Channel(int row, int col)
{
    return col + (row-1)*NROWCOLS;
}

double HadMon::PixelPosition(int rowcol)
{
    return (rowcol-4)*ROWCOL_SPACING;   // cm
}

double HadMon::GetStats(double& xmean, double &ymean,
                        double &xrms, double &yrms)
{
    double qy=0,qx=0,q=0,q2y2=0,q2x2=0,max=0;
    for (int col=1; col<=7; ++col) {
        for (int row=1; row <= 7; ++row) {
            double X = this->PixelPosition(col);
            double Y = this->PixelPosition(row);
            int index = Index(Channel(row,col));
            double Q = this->GetVoltage(index);

            q += Q;
            qx += Q*X;
            qy += Q*Y;
            q2x2 += Q*X*Q*X;
            q2y2 += Q*Y*Q*Y;
            if (max<q) max = q;
        }
    }
    if (q == 0.0) return q;

    xmean = qx/q;
    ymean = qy/q;
    xrms = sqrt(q2x2)/fabs(q);
    yrms = sqrt(q2y2)/fabs(q);
    return q;
}

---