AlgMCDigit.cxx

Go to the documentation of this file.

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

#include <cassert>

#include "Calibrator/Calibrator.h"
#include "CandData/CandRecord.h"
#include "CandDigit/AlgMCDigit.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/RawHeader.h"
#include "RawData/RawRecord.h"
#include "RawData/RawQieMCDigit.h"
#include "RawData/RawVaMCDigit.h"
#include "Validity/VldContext.h"

ClassImp(AlgMCDigit)

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

//______________________________________________________________________
AlgMCDigit::AlgMCDigit()
{
}

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

//______________________________________________________________________
void AlgMCDigit::RunAlg(AlgConfig & /* ac */, CandHandle &ch,
                                                        CandContext &cx)
{
   MSG("Alg", Msg::kDebug) << "Starting AlgMCDigit::RunAlg()" << endl;
   const Calibrator& cal = Calibrator::Instance();

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

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

   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
   t += (digit->GetCrateT0()).GetNanoSec() * Munits::ns;

// And account for snarls overlapping 1 sec boundary.
// 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());
   altlist.ClearWeights();
   altlist.DropZeroWeights();

   const RawMCDigitMixIn *truthDigit =
                            dynamic_cast<const RawMCDigitMixIn*>(digit);
   if (truthDigit) {
     PlexStripEndId seid(truthDigit->GetTrueSEIdEncoded());
     PlexPixelSpotId spotid = ph.GetPixelSpotId(seid);
     altlist.AddStripEndId(seid,spotid,1.0);
   } 

   cdh.SetPlexSEIdAltL(altlist);
}

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

---