BeamMonSpill.cxx

Go to the documentation of this file.

#include "BeamMonSpill.h"

#include <cmath>

#include <DatabaseInterface/DbiResultSet.h>
#include <DatabaseInterface/DbiOutRowStream.h>

#include <Validity/VldRange.h>

#include <MessageService/MsgService.h>
CVSID("$Id: BeamMonSpill.cxx,v 1.17 2007/05/07 19:06:16 mdier Exp $");

BeamMonSpill::BeamMonSpill()
    : fDaeTime(VldTimeStamp(2038,1,1,0,0,0)) // end of time
    , fVmeTime(VldTimeStamp(2038,1,1,0,0,0)) // end of time
    , fTor101(0.0), fTr101d(0.0), fTortgt(0.0), fTrtgtd(0.0)
    , fHornCur(0.0)
    , fTargProfX(0.0), fTargProfY(0.0)
    , fProfWidX(0.0), fProfWidY(0.0)
    , fHadInt(0.0), fMuInt1(0.0), fMuInt2(0.0), fMuInt3(0.0)
{
    fStatus.integer = 0;
    for (int ind=0; ind<6; ++ind) 
        fTargBpmX[ind] = fTargBpmY[ind] = fBpmInt[ind] = 0.0;
}

BeamMonSpill::~BeamMonSpill()
{
}

DbiTableRow* BeamMonSpill::CreateTableRow() const
{
    return new BeamMonSpill;
}

void BeamMonSpill::Fill(DbiResultSet& rs, const DbiValidityRec* /*vrec*/)
{
    int ns=0;

    rs >> fDaeTime;
    rs >> ns;
    fDaeTime.Add(VldTimeStamp(0,ns));

    rs >> fVmeTime;
    rs >> ns;
    fVmeTime.Add(VldTimeStamp(0,ns));

    rs >> fTor101 >> fTr101d >> fTortgt >> fTrtgtd;
    rs >> fHornCur;

    for (int ind=0; ind<6; ++ind) 
        rs >> fTargBpmX[ind] >> fTargBpmY[ind] >> fBpmInt[ind];

    rs >> fTargProfX >> fTargProfY >> fProfWidX >> fProfWidY;
    rs >> fHadInt >> fMuInt1 >> fMuInt2 >> fMuInt3;    
    rs >> fStatus.integer;
}

void BeamMonSpill::Store(DbiOutRowStream& ors, const DbiValidityRec* /*vrec*/) const
{
    ors << fDaeTime << fDaeTime.GetNanoSec()
        << fVmeTime << fVmeTime.GetNanoSec();
    ors << fTor101 << fTr101d << fTortgt << fTrtgtd;
    ors << fHornCur;
    for (int ind=0; ind<6; ++ind) 
        ors << fTargBpmX[ind] << fTargBpmY[ind] << fBpmInt[ind];
    ors << fTargProfX << fTargProfY << fProfWidX << fProfWidY;
    ors << fHadInt << fMuInt1 << fMuInt2 << fMuInt3;    
    ors << fStatus.integer;
}

double BeamMonSpill::BpmAtTarget(double &xmean, double &ymean, double &xrms, double &yrms) const
{
    xmean=ymean=xrms=yrms=0.0;

    double Itot = 0, ix=0,iy=0,ix2=0,iy2=0;
    for (int ind=0; ind<6; ++ind) {
        if (fBpmInt[ind] == 0.0) continue;
        
        double X = fTargBpmX[ind];
        double Y = fTargBpmY[ind];
        double I = fBpmInt[ind];
        I = I > 0 ? I : 0;
        Itot += I;

        ix += I*X;
        iy += I*Y;

        ix2 += I*X*X;
        iy2 += I*Y*Y;
    }

    if (Itot <= 0.0) return 0.0;

    xmean = ix/Itot;
    double xtmp = ix2/Itot-xmean*xmean;
    if (xtmp > 0.0) xrms = sqrt(xtmp);

    ymean = iy/Itot;
    double ytmp = iy2/Itot-ymean*ymean;
    if (ytmp > 0.0) yrms = sqrt(ytmp);

    return Itot;
}


BeamType::BeamType_t BeamMonSpill::BeamType() const
{
    BeamType::BeamType_t btype = BeamType::kUnknown;
    //
    Int_t db_btyp = this->GetStatusBits().beam_type;    

    // Overwrite the beamtype that is in the database for spills
    // after May 2006
    // Change triggered by problems observed for Summer 2006 spills
    // Some spills during pME & pHE running had a beamtype in the
    // database of 1, which means le running. And a large fraction of
    // spills in September had a beamtype of 4, i.e. pME.
    // The most likely casue for this is that the timing used to
    // determine the target position when filling the database was screwed up. 

    double spilltime = this->SpillTime().GetSec();
    if (spilltime > 1149202720 && spilltime <= 1150047812) db_btyp=4;  // pME
    else if (spilltime > 1150047812 && spilltime <= 1155664632 ) db_btyp=5;  // pHE
    else if (spilltime > 1155664632 ) db_btyp=1;  // LE

    
    if (db_btyp==1){
        if (this->SpillTime().GetSec() <= 1112010000){ // early real LE data
            if (fHornCur < -190)                     // nominal at -200kA
                btype = BeamType::kL000z200i;        // 
        }
        else {                                       // later on LE-10 data
            if (fHornCur < -190)                     //  -200 kA
                btype = BeamType::kL010z200i;
            else if (fHornCur < -175){               //  -185kA 
                btype = BeamType::kL010z185i;        //  nominal
            }
            else if (fHornCur < -160)
                btype = BeamType::kL010z170i;
            else if (fHornCur > -5)
                btype = BeamType::kL010z000i;
        }        
    }
    else if (db_btyp==2 || db_btyp==4){                 // (p)ME
        if (this->SpillTime().GetSec() > 1149202720 &&
            this->SpillTime().GetSec() <= 1150047812){  // pME-150
            if (fHornCur < -190)                        // 
                btype = BeamType::kL150z200i;           //  -200 kA
        }
        else {
            if (fHornCur < -190)                        // pME-100 
                btype = BeamType::kL100z200i;           //  -200 kA
        }
    }
    else if (db_btyp==3 || db_btyp==5){             // (p)HE
        if (fHornCur < -190)                        // pHE-250 
            btype = BeamType::kL250z200i;           //  -200 kA
    }

    return btype;
}

VldTimeStamp BeamMonSpill::SpillTime() const
{
  //    if (fStatus.bits.time_source) return fDaeTime;
  if (fStatus.integer&0x100) return fDaeTime;
    return fVmeTime;
}


void BeamMonSpill::SetBPM(const vector<double> &xpos, const vector<double> &ypos,
                          const vector<double> &intensities)
{
    for (size_t ind=0;ind<6 && ind<xpos.size(); ++ind)
        fTargBpmX[ind] = xpos[ind];
    for (size_t ind=0;ind<6 && ind<ypos.size(); ++ind)
        fTargBpmY[ind] = ypos[ind];
    for (size_t ind=0;ind<6 && ind<intensities.size(); ++ind)
        fBpmInt[ind] = intensities[ind];
}


ClassImp(BeamMonSpill)

---