cevtStationRecData_ROOT.cc

Go to the documentation of this file.

#include <cevt/StationRecData.h>
#include <io/cevtStationRecData_ROOT.h>

using namespace io;


ClassImp(cevtStationRecData_ROOT)


cevtStationRecData_ROOT::cevtStationRecData_ROOT(const cevt::StationRecData& d) :
  fTotalMuonHitsPileUpCorrected(d.GetTotalMuonHitsPileUpCorrected()),
  fHasFiducialAreaInfo(d.HasFiducialAreaInfo()),
  fFiducialArea(d.GetFiducialArea()),
  fNumberPadsInFiducialArea(d.GetNumberPadsInFiducialArea()),
  fNumberPadsOutFiducialArea(d.GetNumberPadsOutFiducialArea()),
  fHasSlantMassInfo(d.HasSlantMassInfo()),
  fAverageSlantMass(d.GetAverageSlantMass()),
  fAverageSlantMassInFiducialArea(d.GetAverageSlantMassInFiducialArea()),
  fAverageSlantMassOutFiducialArea(d.GetAverageSlantMassOutFiducialArea())
{ 

  for (int ii = 0; ii != static_cast<int>(cevt::StationConstants::eLastSource); ++ii) {
    fTotalHits[ii] = 
      d.GetTotalHits(static_cast<cevt::StationConstants::SignalComponent>(ii));

    fTotalHitsPileUpCorrected[ii] =
      d.GetTotalHitsPileUpCorrected(static_cast<cevt::StationConstants::SignalComponent>(ii));
      
    fTotalHitsInFiducialArea[ii] = 
      d.GetTotalHitsInFiducialArea(static_cast<cevt::StationConstants::SignalComponent>(ii));
    
    fTotalHitsOutFiducialArea[ii] =
      d.GetTotalHitsOutFiducialArea(static_cast<cevt::StationConstants::SignalComponent>(ii));

    fTotalHitsInFiducialAreaPileUpCorrected[ii] = 
      d.GetTotalHitsInFiducialAreaPileUpCorrected(static_cast<cevt::StationConstants::SignalComponent>(ii));

    fTotalHitsOutFiducialAreaPileUpCorrected[ii] =  
      d.GetTotalHitsOutFiducialAreaPileUpCorrected(static_cast<cevt::StationConstants::SignalComponent>(ii));
  }

  if (d.HasHitsInFiducialArea())
    fHitsInFiducialArea = 
      MultiTimeDistributionI_ROOT(
                           d.HitsInFiducialAreaBegin(), 
                           d.HitsInFiducialAreaEnd()
                           );

  if (d.HasHitsOutFiducialArea())
    fHitsOutFiducialArea = 
      MultiTimeDistributionI_ROOT(
                           d.HitsOutFiducialAreaBegin(), 
                           d.HitsOutFiducialAreaEnd() 
                           );

}


void
cevtStationRecData_ROOT::operator>>(cevt::StationRecData& d)
{
  d.SetTotalMuonHitsPileUpCorrected(fTotalMuonHitsPileUpCorrected);

  for (int ii = 0; ii != static_cast<int>(cevt::StationConstants::eLastSource); ++ii) {
    d.SetTotalHits(fTotalHits[ii],static_cast<cevt::StationConstants::SignalComponent>(ii));
    d.SetTotalHitsPileUpCorrected(fTotalHitsPileUpCorrected[ii],static_cast<cevt::StationConstants::SignalComponent>(ii));
    d.SetTotalHitsInFiducialArea(fTotalHitsInFiducialArea[ii],static_cast<cevt::StationConstants::SignalComponent>(ii)); 
    d.SetTotalHitsOutFiducialArea(fTotalHitsOutFiducialArea[ii],static_cast<cevt::StationConstants::SignalComponent>(ii)); 
    d.SetTotalHitsInFiducialAreaPileUpCorrected(fTotalHitsInFiducialAreaPileUpCorrected[ii],static_cast<cevt::StationConstants::SignalComponent>(ii)); 
    d.SetTotalHitsOutFiducialAreaPileUpCorrected(fTotalHitsOutFiducialAreaPileUpCorrected[ii],static_cast<cevt::StationConstants::SignalComponent>(ii));
  }

  
  d.SetHasFiducialAreaInfo(fHasFiducialAreaInfo);
  d.SetFiducialArea(fFiducialArea);
  d.SetNumberPadsInFiducialArea(fNumberPadsInFiducialArea);
  d.SetNumberPadsOutFiducialArea(fNumberPadsOutFiducialArea);

  d.SetHasSlantMassInfo(fHasSlantMassInfo);
  d.SetAverageSlantMass(fAverageSlantMass);
  d.SetAverageSlantMassInFiducialArea(fAverageSlantMassInFiducialArea);
  d.SetAverageSlantMassOutFiducialArea(fAverageSlantMassOutFiducialArea);


  //d.SetHitsInFiducialArea(fHitsInFiducialArea);
  //d.SetHitsOutFiducialArea(fHitsOutFiducialArea);

  if (fHitsInFiducialArea) {
    for (const auto& st : fHitsInFiducialArea.fTimeDistributions) {
      //utl::TimeDistributionI t(st.second.fSlotSize);
      //st.second >> t;

      const cevt::StationConstants::SignalComponent source = 
        static_cast<cevt::StationConstants::SignalComponent>(st.first);

      if (!d.HasHitsInFiducialArea(source))
         d.MakeHitsInFiducialArea(source);
      st.second >> d.GetHitsInFiducialArea(source);
    }
  }

  if (fHitsOutFiducialArea) {
    for (const auto& st : fHitsOutFiducialArea.fTimeDistributions) {
      //utl::TimeDistributionI t(st.second.fSlotSize);
      //st.second >> t;

      const cevt::StationConstants::SignalComponent source = 
        static_cast<cevt::StationConstants::SignalComponent>(st.first);

      if (!d.HasHitsOutFiducialArea(source))
         d.MakeHitsOutFiducialArea(source);
      st.second >> d.GetHitsOutFiducialArea(source);
    }
  }

}