AlgCaldetDigit.cxx

Go to the documentation of this file.

// $Id: AlgCaldetDigit.cxx,v 1.5 2005/04/08 12:26:29 tagg Exp $
//
// AlgCaldetDigit.cxx
//
// AlgCaldetDigit is a concrete Digit Algorithm class.
//
// Author:  G. Irwin 4/2000

#include <cassert>

#include "Calibrator/Calibrator.h"
#include "CandData/CandRecord.h"
#include "CandDigit/AlgCaldetDigit.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"
#include "CaldetRawDigitWrapper.h"
ClassImp(AlgCaldetDigit)

//______________________________________________________________________
//CVSID("$Id: AlgCaldetDigit.cxx,v 1.5 2005/04/08 12:26:29 tagg Exp $");

//______________________________________________________________________
AlgCaldetDigit::AlgCaldetDigit()
{
}

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


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

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

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

   assert(cx.GetDataIn());
   const CaldetRawDigitWrapper* wrapper = 
     dynamic_cast<const CaldetRawDigitWrapper *> (cx.GetDataIn());

   const RawDigit* digit = wrapper->fDigit;
   const RawDigit* fidDigit = wrapper->fFiducialDigit;
   
   const RawChannelId digitRCId = digit->GetChannel();

   assert(digit);

   // Get VldContext
   const VldContext &vldc = *(cx.GetCandRecord()->GetVldContext());
   
   // Use CandDigitHandle methods to set CandDigit member variables.
   cdh.SetChannelId(digitRCId);
   
   Calibrator& cal = Calibrator::Instance();

   double time =  cal.GetTimeFromTDC(digit->GetTDC(),
                                     digitRCId);

   double fidTime = 0;
   if(fidDigit)
     fidTime =cal.GetTimeFromTDC(fidDigit->GetTDC(),
                                 fidDigit->GetChannel());

   //   cout << digitRCId.AsString("e") << time << "\t" << fidTime << endl;

   // Add on the fiducial time if this is a ND digit.
   if(digitRCId.GetElecType() == ElecType::kQIE) 
     time = time + fidTime;


   Double_t adc_offset  = 0.;
   if (digitRCId.GetElecType() == ElecType::kQIE) {
     adc_offset = 50.;
     cdh.SetQieErrorBits(digit->GetErrorCode());
   }
   else
     cdh.SetVaErrorBits(digit->GetErrorCode());
   
   // Convert TDC counts to Double_t.  Don't subtract Trigger time offset.
   cdh.SetTime(time);
   
   // Set charge to offset-subtracted ADC value ( ~60. * p.e. )
   cdh.SetCharge(((Double_t) digit->GetADC()) - adc_offset);
   
   // 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 AlgCaldetDigit::Trace(const char * /* c */) const
{
}

---