BDHadMuMon.cxx

Go to the documentation of this file.

#include "BDHadMuMon.h"


#include <Conventions/Munits.h>
#include <MessageService/MsgService.h>
CVSID("$Id: BDHadMuMon.cxx,v 1.6 2005/05/16 16:23:02 bv Exp $");


#include <cmath>
#include <string>
using namespace std;

BDHadMuMon::BDHadMuMon()
    : BDSwicDevice(), fNrowcol(0),fNchannels(0),fRowColSpacing(0.0)
{
}
    
BDHadMuMon::BDHadMuMon(const RawBeamData& hadmu_data)
    : BDSwicDevice(hadmu_data), fNrowcol(0),fNchannels(0),fRowColSpacing(0.0)
{
    this->UpdateCache();
}
BDHadMuMon::~BDHadMuMon()
{
}

void BDHadMuMon::SetData(const RawBeamData& hadmu_data)
{
    this->BDSwicDevice::SetData(hadmu_data);
    this->UpdateCache();
}

void BDHadMuMon::UpdateCache()
{
    if (!this->IsValid()) {
        MSG("BD",Msg::kDebug)
            << "UpdateCache called with invalid device data\n";
        return;
    }
    RawBeamData& rbd = this->GetData();
    string name = rbd.GetName();
    if (name == "E:HADMDS") {   // hadron monitor
        fNrowcol = 7;
        fNchannels = 7*7;
        fRowColSpacing = 11.43*Munits::cm;
    }
    else {                      // muon monitor
        fNrowcol = 9;
        fNchannels = 9*9;
        fRowColSpacing = 25.4*Munits::cm;
    }

}

int BDHadMuMon::Index(int channel)
{
    --channel;                  // internally we use 0 based counting

    //   offset  = (00-47)(48-95)
    //   channel = (48-95)(00-47)
    if (channel < 0 || channel >= fNchannels) return -1;
    if (channel < 48) return channel + 48;
    return channel - 48;        // array offset is by definition 0 based
}
int BDHadMuMon::Channel(int row, int col)
{
    --row;                      // internally we use
    --col;                      // 0 based counting

    int chan = -2;

    if (row<0||row>=fNrowcol||col<0||col>=fNrowcol) {
        MSG("BD",Msg::kWarning)
            << "BDHadMuMon::Channel: row/col out of bounds, r,c = "
            << row << "," << col << endl;
        return -2;
    }

    // The PTB decided to transpose the channel numbering between the two....
    if (fNrowcol == 7)          // hadron monitor
        chan = (6-row)*fNrowcol + col;
    else if (fNrowcol == 9)             // muon monitors
        chan = col*fNrowcol + (8-row);
    else
        MSG("BD",Msg::kWarning)
            << "BDHadMuMon::Channel: Unknown number or rows/columns: " << fNrowcol << endl;

    return chan+1;              // externally Had/Mu channes are 1 based counting
}
double BDHadMuMon::PixelPosition(int rowcol)
{
    --rowcol;                   // internally we use 0 based counting
    return (rowcol - (fNrowcol-1)/2) * fRowColSpacing;
}
double BDHadMuMon::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=0; col<fNrowcol; ++col) {
        double X = (col-(fNrowcol-1)/2)* fRowColSpacing;

        for (int row=0; row <fNrowcol; ++row) {
            double Y = (row-(fNrowcol-1)/2)*fRowColSpacing;

            int offset = this->Index(Channel(row+1,col+1));
            double Q = this->GetCharge(offset);

            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;

}

---