TTigressHit.cxx

Go to the documentation of this file.

#include "TTigressHit.h"
 
#include "TClass.h"
 
#include "TTigress.h"
 
TTigressHit::TTigressHit()
{
   Clear();
}
 
TTigressHit::TTigressHit(const TTigressHit& rhs) : TDetectorHit(rhs)
{
   rhs.Copy(*this);
}
 
TTigressHit::TTigressHit(const TFragment& frag)
{
   frag.Copy(*this);
}
 
void TTigressHit::CopyFragment(const TFragment& frag)
{
   frag.Copy(*this);
}
 
void TTigressHit::Clear(Option_t* opt)
{
   TDetectorHit::Clear(opt);
   fTimeFit   = 0.0;
   fSig2Noise = 0.0;
 
   fCoreSet  = false;
   fBgoFired = false;
   fSegments.clear();
}
 
void TTigressHit::Copy(TObject& rhs) const
{
   TDetectorHit::Copy(rhs);
   static_cast<TTigressHit&>(rhs).fTimeFit  = fTimeFit;
   static_cast<TTigressHit&>(rhs).fSegments = fSegments;
   static_cast<TTigressHit&>(rhs).fCoreSet  = fCoreSet;
   static_cast<TTigressHit&>(rhs).fBgoFired = fBgoFired;
}
 
void TTigressHit::Copy(TObject& rhs, bool waveform) const
{
   Copy(rhs);
   if(waveform) {
      CopyWave(rhs);
   }
}
 
TVector3 TTigressHit::GetPosition(Double_t dist) const
{
   return TTigress::GetPosition(*this, dist);
}
 
TVector3 TTigressHit::GetLastPosition(Double_t dist) const
{
   const TDetectorHit* seg = nullptr;
   if(GetNSegments() > 0) {
      seg = &GetSegmentHit(GetNSegments() - 1);   // returns the last segment in the segment vector.
   } else {
      seg = this;   // if no segments, use the core. pcb.
   }
 
   return TTigress::GetPosition(seg->GetDetector(), seg->GetCrystal(), seg->GetSegment(), dist);
}
 
void TTigressHit::Print(Option_t*) const
{
   Print(std::cout);
}
 
void TTigressHit::Print(std::ostream& out) const
{
   std::ostringstream str;
   //TString sopt(opt);
   str << "==== TigressHit @ " << this << " ====" << std::endl;
   str << "\t" << GetName() << std::endl;
   str << "\tCharge:    " << GetCharge() << std::endl;
   str << "\tEnergy:    " << GetEnergy() << std::endl;
   str << "\tTime:      " << GetTime() << std::endl;
   str << "\tCore set:  " << (CoreSet() ? "true" : "false") << std::endl;
   str << "\tBGO Fired: " << (BGOFired() ? "true" : "false") << std::endl;
   str << "\tTime:      " << GetTimeStamp() << std::endl;
   str << "\thit contains " << GetNSegments() << " segments." << std::endl;
   //if(sopt.Contains("all")) {
   str << "Name           Charge" << std::endl;
   for(int x = 0; x < GetNSegments(); x++) {
      str << "\t\t" << GetSegmentHit(x).GetName() << "  |   " << GetSegmentHit(x).GetCharge() << std::endl;
   }
   auto p = GetPosition();
   str << p.GetName() << " " << p.GetTitle() << " (x,y,z)=(" << p.X() << "," << p.Y() << "," << p.Z() << ") (rho,theta,phi)=(" << p.Mag() << "," << p.Theta() << "," << p.Phi() << ")" << std::endl;
   //}
   str << "============================" << std::endl;
   out << str.str();
}
 
bool TTigressHit::Compare(const TTigressHit& lhs, const TTigressHit& rhs)
{
   if(lhs.GetDetector() == rhs.GetDetector()) {
      return (lhs.GetCrystal() < rhs.GetCrystal());
   }
   return (lhs.GetDetector() < rhs.GetDetector());
}
 
bool TTigressHit::CompareEnergy(const TTigressHit& lhs, const TTigressHit& rhs)
{
   return (lhs.GetEnergy()) > rhs.GetEnergy();
}
 
void TTigressHit::SumHit(TTigressHit* hit)
{
   if(this != hit) {
 
      // Should always be true when called by addback construction due to energy ordering during detector construction
      if(GetEnergy() > hit->GetEnergy()) {
         SetTime(GetTime());   // Needs to be call before energy sum to ensure and kIsTimeSet using original energy
                               // for any adjustment
         for(int x = 0; x < hit->GetNSegments(); x++) {
            AddSegment((hit->fSegments[x]));
         }
      } else {
         SetTime(hit->GetTime());
         SetAddress(hit->GetAddress());
         SetCfd(hit->GetCfd());
 
         // Maybe overkill, but consistent
         std::vector<TDetectorHit> fSegmentHold = hit->fSegments;
         for(int x = 0; x < GetNSegments(); x++) {
            fSegmentHold.push_back(fSegments[x]);
         }
         fSegments = fSegmentHold;
      }
 
      SetEnergy(GetEnergy() + hit->GetEnergy());
 
      if(hit->BGOFired()) {
         SetBGOFired(true);
      }
   }
}
 
void TTigressHit::SetWavefit(const TFragment& frag)
{
   TPulseAnalyzer pulse(frag);
   if(pulse.IsSet()) {
      fTimeFit   = static_cast<Float_t>(pulse.fit_newT0());
      fSig2Noise = static_cast<Float_t>(pulse.get_sig2noise());
   }
}
 
void TTigressHit::SetWavefit()
{
   TPulseAnalyzer pulse(*GetWaveform(), 0, GetName());
   if(pulse.IsSet()) {
      fTimeFit   = static_cast<Float_t>(pulse.fit_newT0());
      fSig2Noise = static_cast<Float_t>(pulse.get_sig2noise());
   }
}