StZdcTriggerDetector.cxx

/***************************************************************************
 *
 * $Id: StZdcTriggerDetector.cxx,v 2.14 2006/09/14 00:02:53 ullrich Exp $
 *
 * Author: Thomas Ullrich, Sep 1999
 ***************************************************************************
 *
 * Description:
 *
 ***************************************************************************
 *
 * $Log: StZdcTriggerDetector.cxx,v $
 * Revision 2.14  2006/09/14 00:02:53  ullrich
 * Removed argument (run) in constructor. Not needed anymore.
 *
 * Revision 2.13  2006/08/21 19:43:35  ullrich
 * Run number becomes constructor argument. Needed for ZDC. (Akio)
 *
 * Revision 2.12  2006/02/06 18:52:59  ullrich
 * Fixed bug in constructor. mAdc[3] and mAdc[7] were wrong.
 *
 * Revision 2.11  2004/04/14 22:47:06  ullrich
 * SMD problem fixed.
 *
 * Revision 2.10  2004/04/06 19:39:44  ullrich
 * Added ZDC SMD support.
 *
 * Revision 2.9  2004/02/11 01:42:09  ullrich
 * Added new constructor to load data from StTriggerData.
 *
 * Revision 2.8  2002/03/05 17:11:52  ullrich
 * Corrected bug in constructor. Wrong table entry was used.
 *
 * Revision 2.7  2001/09/14 19:11:11  ullrich
 * Load corrected vertex-Z from trigger table.
 *
 * Revision 2.6  2001/07/12 22:58:33  ullrich
 * Added variable to store the vertex_z from timing info.
 *
 * Revision 2.5  2001/04/05 04:00:59  ullrich
 * Replaced all (U)Long_t by (U)Int_t and all redundant ROOT typedefs.
 *
 * Revision 2.4  2000/07/13 12:51:09  ullrich
 * Added new method numberOfZdcWords to replace old one with wrong name.
 *
 * Revision 2.3  1999/12/21 15:09:28  ullrich
 * Modified to cope with new compiler version on Sun (CC5.0).
 *
 * Revision 2.2  1999/10/28 22:28:21  ullrich
 * Adapted new StArray version. First version to compile on Linux and Sun.
 *
 * Revision 2.1  1999/10/13 19:45:50  ullrich
 * Initial Revision
 *
 **************************************************************************/
#include <algorithm>
#include "StZdcTriggerDetector.h"
#include "tables/St_dst_TrgDet_Table.h"
#include "StTriggerData.h"
#if !defined(ST_NO_NAMESPACES)
using std::fill_n;
using std::copy;
#endif

static const char rcsid[] = "$Id: StZdcTriggerDetector.cxx,v 2.14 2006/09/14 00:02:53 ullrich Exp $";

ClassImp(StZdcTriggerDetector)

StZdcTriggerDetector::StZdcTriggerDetector()
{
    fill_n(mAdc, static_cast<int>(mMaxZdcWords), 0);
    fill_n(mTdc, static_cast<int>(mMaxZdcWords), 0);
    fill_n(mZdcSmdEast, static_cast<int>(mMaxZdcWords), 0);
    fill_n(mZdcSmdWest, static_cast<int>(mMaxZdcWords), 0);
    fill_n(mSumAdc, 2, 0);
    mSum = 0;
    mVertexZ = 0;
}

StZdcTriggerDetector::StZdcTriggerDetector(const dst_TrgDet_st& t)
{
    copy(t.adcZDC+0, t.adcZDC+mMaxZdcWords, mAdc);
    copy(t.tdcZDC+0, t.tdcZDC+mMaxZdcWords, mTdc);
    fill_n(mZdcSmdEast, static_cast<int>(mMaxZdcWords), 0);
    fill_n(mZdcSmdWest, static_cast<int>(mMaxZdcWords), 0);
    mSumAdc[east] = t.adcZDCEast;
    mSumAdc[west] = t.adcZDCWest;
    mSum          = t.adcZDCsum;
    mVertexZ      = t.ZDCvertexZ;
}

StZdcTriggerDetector::StZdcTriggerDetector(const StTriggerData& t)
{
    //
    //  This is a temporary fix only. In future this
    //  class will become obsolete and users should
    //  get this info from StTriggerData.
    //  This info is only a subset of what is available
    //  in StTriggerData.
    //  tu 2/10/2004
    //
    
    // unattenuated
    mAdc[0] = t.zdcUnAttenuated(west); 
    mAdc[1] = t.zdcADC(west, 3);
    mAdc[2] = t.zdcADC(west, 2);
    mAdc[3] = t.zdcADC(west, 1);  // was 2 bug ?
    mAdc[4] = t.zdcUnAttenuated(east); 
    mAdc[5] = t.zdcADC(east, 3);
    mAdc[6] = t.zdcADC(east, 2);
    mAdc[7] = t.zdcADC(east, 1);  // was 2 bug ?
    // TDC
    mAdc[8] = t.zdcTDC(east);
    mAdc[9] = t.zdcTDC(west);
    // attenuated
    mAdc[10] = t.zdcAttenuated(west);
    mAdc[11] = t.zdcAtAddress(11);
    mAdc[12] = t.zdcAtAddress(12);
    mAdc[13] = t.zdcAttenuated(east);
    mAdc[14] = t.zdcAtAddress(14);
    mAdc[15] = t.zdcAtAddress(15);
    mVertexZ = t.zdcVertexZ();

    mSumAdc[east] = t.zdcAttenuated(east);
    mSumAdc[west] = t.zdcAttenuated(west);
    mSum          = t.zdcAtAddress(14);

    for (int i=0; i<mMaxZdcWords; i++) 
        mZdcSmdEast[i] = t.zdcSMD(east, ((i>7) ? 1 : 0), ((i>7) ? (i-7) : (i+1)));
    for (int i=0; i<mMaxZdcWords; i++) 
        mZdcSmdWest[i] = t.zdcSMD(west, ((i>7) ? 1 : 0), ((i>7) ? (i-7) : (i+1)));
    
    fill_n(mTdc, static_cast<int>(mMaxZdcWords), 0); // always 0
}

StZdcTriggerDetector::~StZdcTriggerDetector() {/* noop */}

unsigned int
StZdcTriggerDetector::numberOfZdcCounters() const {return mMaxZdcWords;}

unsigned int
StZdcTriggerDetector::numberOfZdcWords() const {return mMaxZdcWords;}

float
StZdcTriggerDetector::adc(unsigned int i) const
{
    if (i < mMaxZdcWords)
        return mAdc[i];
    else
        return 0;
}

float
StZdcTriggerDetector::tdc(unsigned int i) const
{
    if (i < mMaxZdcWords)
        return mTdc[i];
    else
        return 0;
}

float
StZdcTriggerDetector::adcSum(StBeamDirection dir) const
{
    return mSumAdc[dir];
}


float
StZdcTriggerDetector::zdcSmd(StBeamDirection eastwest, int verthori, int strip) const
{
    
    if (eastwest<0 || eastwest>1) return 0.;
    if (verthori<0 || verthori>1) return 0.;
    if (strip > 8) return 0.;
    
    return ((eastwest<1) ? (mZdcSmdEast[strip -1 + verthori*8]) : (mZdcSmdWest[strip-1 + verthori*8])); 
}


float
StZdcTriggerDetector::adcSum() const {return mSum;}

float
StZdcTriggerDetector::vertexZ() const {return mVertexZ;}

void
StZdcTriggerDetector::setAdc(unsigned int i, float val)
{
    if (i < mMaxZdcWords)
        mAdc[i] = val;
}

void
StZdcTriggerDetector::setTdc(unsigned int i, float val)
{
    if (i < mMaxZdcWords)
        mTdc[i] = val;
}

void
StZdcTriggerDetector::setAdcSum(StBeamDirection dir, float val)
{
    mSumAdc[dir] = val;
}

void
StZdcTriggerDetector::setAdcSum(float val)
{
    mSum = val;
}

void
StZdcTriggerDetector::setVertexZ(float val)
{
    mVertexZ = val;
}

void 
StZdcTriggerDetector::setZdcSmd(StBeamDirection eastwest, int verthori, int strip, float val)
{
    if (eastwest<0 || eastwest>1) return ;
    if (verthori<0 || verthori>1) return ;
    if (strip > 8) return ;
    
    if (eastwest==0) mZdcSmdEast[strip -1 + verthori*8] = val;
    else if (eastwest==1) mZdcSmdWest[strip -1 + verthori*8] = val;
    
}