DigiSignal.cxx

Go to the documentation of this file.

#include "DigiSignal.h"
#include "DigiPE.h"
#include "Conventions/Munits.h"

// TString is needed for Form()
#include "TString.h"

// A global variable to hold
// the id of the next created Signal.
static UInt_t sSignalIdCounter = 0;

DigiSignal::DigiSignal() :
  fCharge(0),
  fPes(0),
  fNumPe(0),
  fHits(0),
  fHitWeights(0),
  fTruth(kUnknown)
{
  fSignalId = sSignalIdCounter++;
}

DigiSignal::DigiSignal(Float_t charge, 
                       PEList_t& pes,
                       HitList_t& hits,
                       HitWeightList_t& weights,
                       DigiPmtTruth_t truth) :
  fCharge(charge),
  fPes(pes),
  fNumPe((int)pes.size()),
  fHits(hits),
  fHitWeights(weights),
  fTruth(truth)
{
  fSignalId = sSignalIdCounter++;
}

// Copy constructor:
// copies the ID, too.
DigiSignal::DigiSignal( const DigiSignal& rhs ) :
  TObject(rhs),
  fSignalId(rhs.fSignalId),
  fCharge(rhs.fCharge),
  fPes(rhs.fPes),
  fNumPe(rhs.GetNumPe()),
  fHits(rhs.fHits),
  fHitWeights(rhs.fHitWeights),
  fTruth(rhs.fTruth)
{
}

DigiSignal::~DigiSignal()
{
  // Note: we don't own the Hits.
}

void DigiSignal::AddDigiPE(const DigiPE* digipe) 
{
  if(digipe) {
    fPes.push_back(digipe);
    const DigiScintHit* hit = digipe->GetHitPointer();
    fNumPe+=1;
    if(hit) AddHit(hit, 1.0);
  }
}

void DigiSignal::AddHit(const DigiScintHit* hitptr, Double_t weight )
{
  //Make sure we don't already have this hit.
  UInt_t n = fHits.size();
  for(UInt_t i=0; i<n; i++) {
    const DigiScintHit* oldhit = GetHit(i);
    if (oldhit->GetUniqueID() == hitptr->GetUniqueID()) {
      // Add this weight.
      fHitWeights[i] += weight;
      return;
    }
  }

  // Add it onto the list
  TRef ref((TObject*)hitptr);
  fHits.push_back(ref);
  fHitWeights.push_back(weight);
}

void DigiSignal::Merge(DigiSignal& other)
{
  fCharge += other.fCharge;
  fPes.insert(fPes.end(),other.fPes.begin(),other.fPes.end());
  fNumPe+=other.GetNumPe();

  UInt_t n = other.GetNumberOfHits();
  for(UInt_t i=0; i<n; i++) {
    AddHit(other.GetHit(i), other.GetHitWeight(i));
  }

  fTruth |= other.fTruth;
}

//_____________________________________________________________________________
const char* AsString(DigiSignal::DigiPmtTruth_t truth)
{
  // Global function to get a string.

  if(truth == DigiSignal::kUnknown) return "|Unknown|";
   
   std::string result("|");
   if(truth & DigiSignal::kGenuine)            result = result + "Genuine|";
   if(truth & DigiSignal::kDarkNoise)          result = result + "DarkNoise|";
   if(truth & DigiSignal::kFibreLight)         result = result + "FibreLight|";
   if(truth & DigiSignal::kCrosstalk)          result = result + "CrossTalk|";
   if(truth & DigiSignal::kCrosstalkOptical)   result = result + "CrossTalkOptical|";
   if(truth & DigiSignal::kLeakFromNextBucket) result = result + "LeakFromNextBucket|";
   if(truth & DigiSignal::kLeakFromPrevBucket) result = result + "LeakFromPrevBucket|";
   // Hack: use ROOT's circular buffer of strings. neat,huh?
   return Form("%s",result.c_str()); 
}

//......................................................................

ostream& operator<<(ostream& os, const DigiSignal& s) 
{ return s.FormatToOStream(os); }

//_____________________________________________________________________________
void DigiSignal::Print(Option_t *option) const
{
   std::cout << this->AsString(option) << std::endl;
}

//_____________________________________________________________________________
std::ostream& DigiSignal::FormatToOStream(std::ostream& os,
                                          Option_t *option) const
{
   os << AsString(option);
   return os;
}

//_____________________________________________________________________________
const char* DigiSignal::AsString(Option_t *option) const
{
  const char* fmt_default = 
    "DigiSignal [%d] %s Q=%.1ffC Npe=%d Nhits=%d";
  const char* fmt_extended =
    "DigiSignal ID=%d Truth=%s Charge=%.1ffC Npe=%d Nhits=%d";

  const char* fmt = fmt_default;

  switch(option[0]) {
  case 'e': // (E)xtended printout.
  case 'E':
    fmt = fmt_extended;
    break;
  default:
    fmt = fmt_default;
  }
   
  DigiSignal::DigiPmtTruth_t truth = GetTruth();
  const char* truthstr = ::AsString(truth);
  return Form(fmt,
              GetId(),
              truthstr,
              GetCharge()/Munits::fC,
              GetNumPe(),
              GetNumberOfHits()
              );

}

---