THILFragment.cxx

Go to the documentation of this file.

#include "THILFragment.h"
#include "Globals.h"
 
THILFragment::THILFragment() : TFragment()
{
   /// Default constructor
   Clear();
}
 
THILFragment::THILFragment(const THILFragment& rhs) : TFragment(rhs)
{
   /// Copy constructor
   fGermaniumId     = rhs.fGermaniumId;
   fGermaniumEnergy = rhs.fGermaniumEnergy;
   fGermaniumTime   = rhs.fGermaniumTime;
   fSiliconId       = rhs.fSiliconId;
   fSiliconEnergy   = rhs.fSiliconEnergy;
   fEventNumber     = rhs.fEventNumber;
   fEventTime       = rhs.fEventTime;
   fUsTime          = rhs.fUsTime;
}
 
THILFragment::~THILFragment()
{
   /// Default destructor
}
 
void THILFragment::Clear(Option_t* opt)
{
   /// Clears all fields of the THILFragment and TFragment
   TFragment::Clear(opt);
   fGermaniumId.clear();
   fGermaniumEnergy.clear();
   fGermaniumTime.clear();
   fSiliconId.clear();
   fSiliconEnergy.clear();
   fEventNumber = 0;
   fEventTime   = 0;
   fUsTime      = 0;
}
 
TObject* THILFragment::Clone(const char*) const
{
   auto result = new THILFragment;
   *result     = *this;
   result->ClearTransients();
   return result;
}
 
bool THILFragment::Good(bool verbose) const
{
   /// This function checks whether this event is good, i.e. the length of all germanium
   /// and all silicon vectors matches and there is at least one silicon and one germanium hit.
 
   if(fGermaniumId.empty() || fSiliconId.empty()) {
      if(verbose) std::cout << this << ": Either no germanium IDs (" << fGermaniumId.size() << ") or no silicon IDs (" << fSiliconId.size() << ")" << std::endl;
      return false;
   }
   if(fGermaniumId.size() != fGermaniumEnergy.size() || fGermaniumId.size() != fGermaniumTime.size()) {
      if(verbose) std::cout << this << ": Size of germanium IDs (" << fGermaniumId.size() << ") doesn't match size of germanium energies (" << fGermaniumEnergy.size() << ") or size of germanium times (" << fGermaniumTime.size() << ")" << std::endl;
      return false;
   }
   if(fSiliconId.size() != fSiliconEnergy.size()) {
      if(verbose) std::cout << this << ": Size of silicon IDs (" << fSiliconId.size() << ") doesn't match size of silicon energies (" << fSiliconEnergy.size() << ")" << std::endl;
      return false;
   }
 
   return true;
}
 
void THILFragment::Print(std::ostream& out) const
{
   std::stringstream str;
   str << "THILFragment " << this << ", #" << fEventNumber << ": " << fGermaniumId.size() << " germaniums and " << fSiliconId.size() << " silicons at event time " << fEventTime << ", and microsecond time " << fUsTime << std::endl;
   str << "#   Ge# GeEn GeTime Si# SiEn" << std::endl;
   for(size_t h = 0; h < fGermaniumId.size(); ++h) {
      if(h < fSiliconId.size()) {
         str << std::setw(2) << h << "  " << std::setw(3) << (int)fGermaniumId.at(h) << " " << std::setw(4) << fGermaniumEnergy.at(h) << " " << std::setw(6) << fGermaniumTime.at(h) << " " << std::setw(3) << (int)fSiliconId.at(h) << " " << std::setw(4) << fSiliconEnergy.at(h) << std::endl;
      } else {
         str << std::setw(2) << h << "  " << std::setw(3) << (int)fGermaniumId.at(h) << " " << std::setw(4) << fGermaniumEnergy.at(h) << " " << std::setw(6) << fGermaniumTime.at(h) << std::endl;
      }
   }
   str << "============================" << std::endl;
   out << str.str();
}