In This Package:

CalibReadoutEventCnv.cc

Go to the documentation of this file.

#include "CalibReadoutEventCnv.h"
#include "PerBaseEvent/HeaderObjectCnv.h"
#include "Conventions/Electronics.h"
#include "Context/TimeStamp.h"

#include <vector>

using namespace DayaBay;
using namespace std;

CalibReadoutHeaderCnv::CalibReadoutHeaderCnv(ISvcLocator* svc)
      : RootIOTypedCnv<PerCalibReadoutHeader,CalibReadoutHeader>(
                                                       "PerCalibReadoutHeader",
                                                       classID(),svc)
{
}


CalibReadoutHeaderCnv::~CalibReadoutHeaderCnv()
{
}

StatusCode CalibReadoutHeaderCnv::PerToTran(const PerCalibReadoutHeader& perobj,
                                          DayaBay::CalibReadoutHeader& tranobj)
{
  MsgStream log(msgSvc(), "CalibReadoutHeaderCnv::PerToTran");
  StatusCode sc = HeaderObjectCnv::toTran(perobj,tranobj);
  if (sc.isFailure()) return sc;

  delete tranobj.calibReadout();
  if(perobj.site){
    CalibReadout* readout = 0;
    Detector detector(Site::Site_t(perobj.site),
                      DetectorId::DetectorId_t(perobj.detector));
    if(detector.isAD() || detector.isWaterShield()){ 
      CalibReadoutPmtCrate *out_crate = new CalibReadoutPmtCrate();
      readout = out_crate;
      out_crate->setDetector(detector);
      out_crate->setTriggerType( Trigger::TriggerType_t(perobj.triggerType) );
      out_crate->setTriggerNumber( perobj.eventNumber );
      timespec value = {perobj.triggerTimeSec,perobj.triggerTimeNanoSec};
      TimeStamp tstamp(value);
      out_crate->setTriggerTime(tstamp);
      
      CalibReadoutPmtCrate::PmtChannelReadouts out_channels;
      
      CalibReadoutPmtChannel* out_channel = 0;
      vector<double> timeArr;
      vector<double> chargeArr;
      if(detector.isAD()){
        // Fill AD hits
        for(unsigned int hitIdx=0; hitIdx<perobj.nHitsAD; hitIdx++){
          AdPmtSensor pmtId(perobj.ring[hitIdx], perobj.column[hitIdx],
                            detector.site(),detector.detectorId());
          if(!out_channel || (pmtId!=out_channel->pmtSensorId())){
            if(out_channel){
              out_channel->setTime(timeArr);
              out_channel->setCharge(chargeArr);
            }
            // create channel
            timeArr.clear();
            chargeArr.clear();
            CalibReadoutPmtChannel cur_channel(pmtId,out_crate);
            out_channels.push_back(cur_channel);
            out_channel=&(out_channels[out_channels.size()-1]);
          }
          timeArr.push_back(perobj.timeAD[hitIdx]);
          chargeArr.push_back(perobj.chargeAD[hitIdx]);
        }
        // Catch 2-inch calibration PMTs
        for(unsigned int hitIdx=0; hitIdx<perobj.nHitsAD_calib; hitIdx++){
          AdPmtSensor pmtId(0/*ring==0 for calib PMTs*/, 
                            perobj.topOrBottom[hitIdx] 
                            + 2*(perobj.acuColumn[hitIdx]-1),
                            detector.site(),detector.detectorId());
          if(!out_channel || (pmtId!=out_channel->pmtSensorId())){
            if(out_channel){
              out_channel->setTime(timeArr);
              out_channel->setCharge(chargeArr);
            }
            // create channel
            timeArr.clear();
            chargeArr.clear();
            CalibReadoutPmtChannel cur_channel(pmtId,out_crate);
            out_channels.push_back(cur_channel);
            out_channel=&(out_channels[out_channels.size()-1]);
          }
          timeArr.push_back(perobj.timeAD_calib[hitIdx]);
          chargeArr.push_back(perobj.chargeAD_calib[hitIdx]);
        }
        // Catch last channel
        if(out_channel){
          out_channel->setTime(timeArr);
          out_channel->setCharge(chargeArr);
        }
      }else if(detector.isWaterShield()){
        // Fill Water Pool hits
        for(unsigned int hitIdx=0; hitIdx<perobj.nHitsPool; hitIdx++){
          PoolPmtSensor pmtId(perobj.wallNumber[hitIdx], 
                              perobj.wallSpot[hitIdx],
                              perobj.inwardFacing[hitIdx],
                              detector.site(),detector.detectorId());
          if(!out_channel || (pmtId!=out_channel->pmtSensorId())){
            if(out_channel){
              out_channel->setTime(timeArr);
              out_channel->setCharge(chargeArr);
            }
            // create channel
            timeArr.clear();
            chargeArr.clear();
            CalibReadoutPmtChannel cur_channel(pmtId,out_crate);
            out_channels.push_back(cur_channel);
            out_channel=&(out_channels[out_channels.size()-1]);
          }
          timeArr.push_back(perobj.timePool[hitIdx]);
          chargeArr.push_back(perobj.chargePool[hitIdx]);
        }
        // Catch last channel
        if(out_channel){
          out_channel->setTime(timeArr);
          out_channel->setCharge(chargeArr);
        }
      }
      out_crate->setChannelReadout(out_channels);
      out_crate->setHeader( &tranobj );
    }else if(detector.detectorId() == DetectorId::kRPC){
      // FIXME: Add RPC here
      readout = 0;
    }

    tranobj.setCalibReadout( readout );
  }else{
    // Header has no readout
    tranobj.setCalibReadout(0);
  }

  return StatusCode::SUCCESS;
}                                       

StatusCode CalibReadoutHeaderCnv::TranToPer(
                                    const DayaBay::CalibReadoutHeader& tranobj,
                                    PerCalibReadoutHeader& perobj)
{
  MsgStream log(msgSvc(), "CalibReadoutHeaderCnv::TranToPer");
  StatusCode sc = HeaderObjectCnv::toPer(tranobj,perobj);
  if (sc.isFailure()) return sc;

  // Clear persistent object
  perobj.clear();

  const CalibReadout* readout = tranobj.calibReadout();
  if(!readout){
    return StatusCode::SUCCESS;
  }
  
  // Convert trigger information
  perobj.site = (unsigned int)(readout->detector().site());
  perobj.detector = (unsigned int)(readout->detector().detectorId());
  perobj.eventNumber = readout->triggerNumber();
  perobj.triggerType = readout->triggerType();
  perobj.triggerTimeSec =readout->triggerTime().GetSec();
  perobj.triggerTimeNanoSec = readout->triggerTime().GetNanoSec();
  
  // Convert hit information
  if(readout->detector().detectorId() == DetectorId::kUnknown){
    // Unknown detector, skip
    ; 
  }else if( readout->detector().isAD() || readout->detector().isWaterShield() ){
    // PMT data conversion
    const CalibReadoutPmtCrate& pmtCrate = 
      dynamic_cast<const CalibReadoutPmtCrate&>(*readout);
    const CalibReadoutPmtCrate::PmtChannelReadouts& chReadouts 
      = pmtCrate.channelReadout();
    CalibReadoutPmtCrate::PmtChannelReadouts::const_iterator channelIter,
      channelDone=chReadouts.end();
    if(readout->detector().isAD()){
      // Loop over AD hits
      for(channelIter=chReadouts.begin();channelIter!=channelDone;
          ++channelIter){
        const CalibReadoutPmtChannel& channel = *channelIter;
        AdPmtSensor pmtId(channel.pmtSensorId().fullPackedData());
        unsigned int nHits = channel.size();
        if(pmtId.ring()==0){
          // Catch 2inch calibration PMTs, save separately
          for(unsigned int hitIdx=0; hitIdx<nHits; hitIdx++){
            perobj.timeAD_calib.push_back( channel.time(hitIdx) );
            perobj.chargeAD_calib.push_back( channel.charge(hitIdx) );
            perobj.hitCountAD_calib.push_back( hitIdx );
            perobj.topOrBottom.push_back( ((pmtId.column()+1)%2) + 1 );
            perobj.acuColumn.push_back( (pmtId.column()+1)/2 );
            perobj.nHitsAD_calib++;
          }
        }else{
          // 8inch AD PMT
          for(unsigned int hitIdx=0; hitIdx<nHits; hitIdx++){
            perobj.timeAD.push_back( channel.time(hitIdx) );
            perobj.chargeAD.push_back( channel.charge(hitIdx) );
            perobj.hitCountAD.push_back( hitIdx );
            perobj.ring.push_back( pmtId.ring() );
            perobj.column.push_back( pmtId.column() );
            perobj.nHitsAD++;
          }       
        }
      }
    }else if(readout->detector().isWaterShield()){
      // Loop over Water Pool hits
      for(channelIter=chReadouts.begin();channelIter!=channelDone;
          ++channelIter){
        const CalibReadoutPmtChannel& channel = *channelIter;
        PoolPmtSensor pmtId(channel.pmtSensorId().fullPackedData());
        unsigned int nHits = channel.size();
        for(unsigned int hitIdx=0; hitIdx<nHits; hitIdx++){
          perobj.timePool.push_back( channel.time(hitIdx) );
          perobj.chargePool.push_back( channel.charge(hitIdx) );
          perobj.hitCountPool.push_back( hitIdx );
          perobj.wallNumber.push_back( pmtId.wallNumber() );
          perobj.wallSpot.push_back( pmtId.wallSpot() );
          perobj.inwardFacing.push_back( pmtId.inwardFacing() );
          perobj.nHitsPool++;
        }
      }
    }
  }else{
    // RPC data conversion
  }

  return StatusCode::SUCCESS;
}

StatusCode CalibReadoutHeaderCnv::fillRepRefs(IOpaqueAddress*, DataObject* dobj)
{
    MsgStream log(msgSvc(), "CalibReadoutHeaderCnv::fillRepRefs");
    CalibReadoutHeader* rh = dynamic_cast<CalibReadoutHeader*>(dobj);

    log << MSG::DEBUG
        << "Saving links to " << rh->inputHeaders().size() 
        << " input headers" << endreq;

    StatusCode sc = HeaderObjectCnv::fillPer(m_rioSvc,*rh,*m_perOutObj);
    if (sc.isFailure()) {
        log << MSG::ERROR << "Failed to fill HeaderObject part" << endreq;
        return sc;
    }

    // ... fill CalibReadoutHeader part...
    return sc;
}

StatusCode CalibReadoutHeaderCnv::fillObjRefs(IOpaqueAddress*, DataObject* dobj)
{
    MsgStream log(msgSvc(), "CalibReadoutHeaderCnv::fillObjRefs");
    HeaderObject* hobj = dynamic_cast<HeaderObject*>(dobj);
    StatusCode sc = HeaderObjectCnv::fillTran(m_rioSvc,*m_perInObj,*hobj);
    if (sc.isFailure()) {
        log << MSG::ERROR << "Failed to fill HeaderObject part" << endreq;
        return sc;
    }

    log << MSG::DEBUG
        << "Restored links to " << hobj->inputHeaders().size() 
        << " input headers" << endreq;

    // ... fill CalibReadoutHeader part...
    return sc;
}

---