StFpsRawHitMaker.cxx

// \class StFpsRawHitMaker
// \author Akio Ogawa
//
//  $Id: StFpsRawHitMaker.cxx,v 1.4 2017/01/30 18:10:16 akio Exp $
//  $Log: StFpsRawHitMaker.cxx,v $
//  Revision 1.4  2017/01/30 18:10:16  akio
//  remove LOG_INFO
//
//  Revision 1.3  2017/01/30 17:49:28  akio
//  adding FPost
//
//  Revision 1.2  2015/09/02 14:55:18  akio
//  Modified to work with StFmsDbMaker
//

#include "StFpsRawHitMaker.h"

#include "St_base/StMessMgr.h"
#include "St_base/Stypes.h"
#include "RTS/src/DAQ_FPS/daq_fps.h"
#include "RTS/src/DAQ_READER/daq_dta.h"
#include "StChain/StRtsTable.h"
#include "StEvent/StEvent.h"
#include "StEvent/StFmsCollection.h"
#include "StEvent/StFmsHit.h"
#include "StRoot/StFmsDbMaker/StFmsDbMaker.h"

StFpsRawHitMaker::StFpsRawHitMaker(const Char_t* name) :
  StRTSBaseMaker( "adc", name ), mPrePost(0), mFmsCollection(0), mFmsDbMaker(0) {}

StFpsRawHitMaker::~StFpsRawHitMaker(){}

Int_t StFpsRawHitMaker::InitRun(Int_t runNumber){
  LOG_DEBUG << "StFpsRawHitMaker::InitRun with run = "  << runNumber << endm;
  mFmsDbMaker = static_cast<StFmsDbMaker*>(GetMaker("fmsDb")); 
  if(!mFmsDbMaker){
    LOG_ERROR  << "StFpsRawHitMaker::InitRun Failed to get StFmsDbMaker" << endm;
    return kStFatal;
  }
  return kStOK;
}

Int_t StFpsRawHitMaker::prepareEnvironment(){
  StEvent* eventPtr=0;
  eventPtr= (StEvent*)StRTSBaseMaker::GetInputDS("StEvent");
  mFmsCollection=0;
  if(eventPtr){
    LOG_DEBUG <<"StFpsRawHitMaker::prepareEnvironment found StEvent" <<endm;
    mFmsCollection=eventPtr->fmsCollection();
  }else{
    LOG_ERROR <<"StFpsRawHitMaker::prepareEnvironment found no StEvent"<<endm;
    return kStFatal;
  }
  if(!mFmsCollection){
    LOG_DEBUG <<"StFpsRawHitMaker::prepareEnvironment creating new StFmsCollection" <<endm;
    mFmsCollection=new StFmsCollection;
    if(!mFmsCollection){
      LOG_ERROR <<"StFpsRawHitMaker::prepareEnvironment could not create StFmsCollection" <<endm;
      return kStFatal;
    }
    LOG_DEBUG <<"StFpsRawHitMaker::prepareEnvironment setting new StFmsCollection to StEvent" <<endm;
    eventPtr->setFmsCollection(mFmsCollection);
  }  
  return kStOK;
}

Int_t StFpsRawHitMaker::Make(){
  StRtsTable* rts_tbl=0;
  int n[2]={0,0}, ngood[2]={0,0};
  if(int ret=prepareEnvironment()!=kStOK ) return ret;

  // loop over DaqElements of FPS data in daq file
  // A DaqElement contains data from one QT baord and one bunch crossings
  while((rts_tbl = GetNextDaqElement("fps/adc"))){    

    // 2015  
    //int xing=rts_tbl->Sector(); 
    //if(xing>=128) xing-=256;
    //int qt=rts_tbl->Rdo();
 
    // 2017
    int fpsfpost=rts_tbl->Sector();
    int xing=rts_tbl->Pad();  if(xing>=128) xing-=256;
    int qt=rts_tbl->Row(); 

    int ndata=rts_tbl->GetNRows();
    LOG_DEBUG << Form("FPS: fps/fpost=%1d xing=%4d QT=%02d NData=%d",fpsfpost,xing,qt,ndata)<<endm;
    if(mPrePost==0 && xing!=0) continue;  
    for(StRtsTable::iterator it=rts_tbl->begin(); it!=rts_tbl->end(); it++){
      fps_adc_t *a=(fps_adc_t *)*it;
      //fps_adc_t *a=static_cast<fps_adc_t *>(*it);  //would not work! use "reinterpret_cast".
      int ch=a->ch;
      int adc=a->adc;
      int tdc=a->tdc;  
      int slatid,q,l,s,det,crate;
      if(fpsfpost==1){
          slatid = mFmsDbMaker->fpsSlatidFromQT(qt,ch); //Get SlatId from QT address and channel
          mFmsDbMaker->fpsQLSfromSlatId(slatid,&q,&l,&s); //Get Quad/Layer/Slat#s from SlatId
          LOG_DEBUG << Form("FPS: xing=%4d QT%02d ch%02d Slaiid=%3d Q%1dL%1dS%02d ADC=%4d TDC=%2d",xing,qt,ch,slatid,q,l,s,adc,tdc)<<endm;
          det=kFpsDetId;
          crate=kFpsQtCrate;
      }else{
          slatid = mFmsDbMaker->fpostSlatidFromQT(qt,ch); //Get SlatId from QT address and channel
          mFmsDbMaker->fpostQLSfromSlatId(slatid,&q,&l,&s); //Get Quad/Layer/Slat#s from SlatId
          LOG_DEBUG << Form("FPO: xing=%4d QT%02d ch%02d Slaiid=%3d Q%1dL%1dS%02d ADC=%4d TDC=%2d",xing,qt,ch,slatid,q,l,s,adc,tdc)<<endm;
          det=kFpostDetId;
          crate=kFpostQtCrate;
      }
      n[fpsfpost-1]++;
      int flag=0;       
      if(slatid<0)          { /* LOG_WARN << "Invalid SlatId = "<<slatid<<endm;*/      flag=1; }
      if(q<1 || l<1 || s<1) { /* LOG_WARN << Form("Invalid Q/L/S = %d/%d/%d",q,l,s);*/ flag=1; }
      if(flag==0){
        StFmsHit* hit = new StFmsHit();
        hit->setDetectorId(det);         
        hit->setChannel(slatid);
        hit->setQtCrate(crate);
        hit->setQtSlot(qt);
        hit->setQtChannel(ch);
        hit->setAdc(adc);
        hit->setTdc(tdc);
        float gain=.0;
        if(fpsfpost==0) {mFmsDbMaker->fpsGain(slatid);}
        else            {mFmsDbMaker->fpostGain(slatid);}
        float nmip=0.0;
        if(gain>0.0) nmip=adc/gain;
        hit->setEnergy(nmip);
        mFmsCollection->addHit(hit);
        if(Debug()) hit->print();        
        ngood[fpsfpost-1]++;
      }
    }
  }
  LOG_DEBUG << Form("StFpsRawHitMaker:: FPS   found total %d hits and %d good",n[0],ngood[0]) << endm;
  LOG_DEBUG<< Form("StFpsRawHitMaker:: FPOST found total %d hits and %d good",n[1],ngood[1]) << endm;
  return kStOK;
}

ClassImp(StFpsRawHitMaker)