AlgDigit.cxx

Go to the documentation of this file.

// $Id: AlgDigit.cxx,v 1.26 2009/10/12 18:39:08 tagg Exp $
//
// AlgDigit.cxx
//
// AlgDigit is a concrete Digit Algorithm class.
//
// Author:  G. Irwin 4/2000

#include <cassert>

#include "Calibrator/Calibrator.h"
#include "CandData/CandRecord.h"
#include "CandDigit/AlgDigit.h"
#include "CandDigit/CandDigitHandle.h"
#include "Candidate/CandContext.h"
#include "Conventions/Munits.h"
#include "MessageService/MsgService.h"
#include "Plex/PlexHandle.h"
#include "Plex/PlexSEIdAltL.h"
#include "RawData/RawChannelId.h"
#include "RawData/RawDigit.h"
#include "RawData/RawHeader.h"
#include "RawData/RawRecord.h"
#include "Validity/VldContext.h"

ClassImp(AlgDigit)

//______________________________________________________________________
//CVSID("$Id: AlgDigit.cxx,v 1.26 2009/10/12 18:39:08 tagg Exp $");

//______________________________________________________________________
AlgDigit::AlgDigit()
{
}

//______________________________________________________________________
AlgDigit::~AlgDigit()
{
}

//______________________________________________________________________
void AlgDigit::RunAlg(AlgConfig & /* ac */, CandHandle &ch,
                                                        CandContext &cx)
{

#ifdef ENABLE_DEBUG_MSG
   MSG("Alg", Msg::kDebug) << "Starting AlgDigit::RunAlg()" << endl;
#endif

   assert(ch.InheritsFrom("CandDigitHandle"));
   CandDigitHandle &cdh = (CandDigitHandle &) ch;

   const RawDigit *digit =
                        dynamic_cast<const RawDigit *> (cx.GetDataIn());
   assert(digit);

   const Calibrator& cal = Calibrator::Instance();
   
   const RawChannelId digitRCId = digit->GetChannel();

// Get VldContext
   const VldContext &vldc = *(cx.GetCandRecord()->GetVldContext());

// Use CandDigitHandle methods to set CandDigit member variables.
   cdh.SetChannelId(digitRCId);

// Take off a 50 ADC offset from QIE channels.
   Double_t adc_offset  = 0.;
   if (digitRCId.GetElecType()==ElecType::kQIE) {
     adc_offset = 50.;
     cdh.SetQieErrorBits(digit->GetErrorCode());
   }
   else
     cdh.SetVaErrorBits(digit->GetErrorCode());

// Set charge to offset-subtracted ADC value ( ~60. * p.e. )
   cdh.SetCharge(((Double_t) digit->GetADC()) - adc_offset);

// Convert TDC counts to Double_t.  Don't subtract Trigger time offset.
   Double_t t = cal.
                 GetTimeFromTDC(digit->GetTDC(), digitRCId);

// Add CrateT0 (A.Weber)
// this is needed to account for the counter->reset hold-off during the spill
// Removed again my N.Tagg. This does the wrong thing in 99.9% of all cases,
// as documented in Doc-DB 1870.  
//   t += (digit->GetCrateT0()).GetNanoSec() * Munits::ns;

// And account for snarls overlapping 1 sec boundary. (A.Weber)
// This assumes that the crate t0 is the t0 coresponding to the hit.
// t now be <0 or >1
   t += ( (digit->GetCrateT0()).GetSec() -
          (vldc.GetTimeStamp()).GetSec() ) * Munits::s;

#ifdef ENABLE_DEBUG_MSG
   MSG("Alg", Msg::kDebug) << "rdid = " << digitRCId << endl
                        << "    tdc = " << digit->GetTDC() << " nsec = "
   << cal.GetTimeFromTDC(digit->GetTDC(),
                  digitRCId) << " t = " << t << endl << "    CrT0_ns = "
                        << (digit->GetCrateT0()).GetNanoSec()*Munits::ns
                       << " CrT0_s = " << (digit->GetCrateT0()).GetSec()
          << " TrigTime_s = " << (vldc.GetTimeStamp()).GetSec() << endl;
#endif

   cdh.SetTime(t);

// Get a PlexHandle for PlexStripEndId to RawChannelId conversions
   PlexHandle ph(vldc);

// Fill and calibrate PlexSEIdAltL
   PlexSEIdAltL altlist =
                      ph.GetSEIdAltL(digitRCId, &cal,
                      (int)(digit->GetADC()-adc_offset), cdh.GetTime());

   cdh.SetPlexSEIdAltL(altlist);
}

//______________________________________________________________________
void AlgDigit::Trace(const char * /* c */) const
{
}

---