Co60Data.cc

Go to the documentation of this file.

#include "Co60Data.h"

#include "Event/ReadoutHeader.h"
#include "Event/CalibReadoutHeader.h"
#include "Event/CalibReadout.h"
#include "Event/UserDataHeader.h"
#include "Event/RecHeader.h"
#include "Event/RecTrigger.h"
#include "Event/HeaderObject.h"

#include "GaudiKernel/SystemOfUnits.h"

#include "TMath.h"
#include "TString.h"

Co60Data::Co60Data(const std::string& name, ISvcLocator* pSvcLocator)
  : DybAlgorithm<DayaBay::UserDataHeader>(name,pSvcLocator) {
  m_location = "/Event/Data/Calib/Co60";

  declareProperty("ReadoutHeaderLocation", m_readoutHeaderLocation="/Event/Readout/ReadoutHeader",
                  "Location of the ReadoutHeader");
  declareProperty("CalibStatsLocation", m_calibStatsLocation="/Event/Data/CalibStats",
                  "Location of the CalibStats User data");
  declareProperty("Co60Location", m_Co60Location="/Event/Tag/Calib/Co60",
                  "Location of the Co60 Tag");
}

Co60Data::~Co60Data() {
}

StatusCode Co60Data::initialize() {
  return this->GaudiAlgorithm::initialize();
}

StatusCode Co60Data::finalize() {
  return this->GaudiAlgorithm::finalize();
}

StatusCode Co60Data::execute() {
  
  StatusCode sc = StatusCode::SUCCESS;

  const DayaBay::ReadoutHeader *readoutHdr = 0; 
  const DayaBay::RecHeader *recHdr = 0; 
  DayaBay::UserDataHeader *calibStatsHdr = 0; 

  // Check if it is tagged as Co60, if not, skip.
  bool isCo60 = exist<DayaBay::HeaderObject>(evtSvc(), m_Co60Location);
  if (!isCo60){
    return StatusCode::SUCCESS;
  }

  // Get Calibstats
  if ( !exist<DayaBay::UserDataHeader>(evtSvc(), m_calibStatsLocation) ) {
    warning() << "Cannot find header at " << m_calibStatsLocation << endreq;
    return StatusCode::FAILURE;
  } else {
    calibStatsHdr = get<DayaBay::UserDataHeader>(m_calibStatsLocation);
  }
  
  // Get readout
  if (!exist<DayaBay::ReadoutHeader>(evtSvc(), "/Event/Readout/ReadoutHeader")) {
    warning() << "Failed to get readout" << endreq;
  } else {
    readoutHdr = get<DayaBay::ReadoutHeader>("/Event/Readout/ReadoutHeader");
  }

  // Get recon header
  if (!exist<DayaBay::RecHeader>(evtSvc(), "/Event/Rec/AdSimple")) {
    warning() << "Failed to get recheader" << endreq;
    return StatusCode::FAILURE;
  } else {
    recHdr = get<DayaBay::RecHeader>("/Event/Rec/AdSimple");
  }
  const DayaBay::RecTrigger& rectrigger = recHdr->recTrigger();
  
  // Skip SAB Data 
  Site::Site_t siteId = readoutHdr->context().GetSite();
  if (siteId == Site::kSAB) {
    return StatusCode::SUCCESS;
  }

  // Skip if recon has failed
  float recES = rectrigger.energyStatus();
  float recPS = rectrigger.positionStatus();
  if (recES != 1  || recPS != 1){
    debug() << "Recon failed. skipping" << endreq;
    return StatusCode::SUCCESS;
  }

  // Record data
  // -- 1a.  Get data from recheader
  float e_co60 = rectrigger.energy()/Gaudi::Units::MeV;
  float x_co60 = rectrigger.position().x() / Gaudi::Units::mm;
  float y_co60 = rectrigger.position().y() / Gaudi::Units::mm;
  float z_co60 = rectrigger.position().z() / Gaudi::Units::mm;
  // -- 1b.  Get data from calib stats
  float nPESUM = calibStatsHdr->getFloat("nPESUM");
  float nPERMS = calibStatsHdr->getFloat("nPERMS");

  // -- 2. Create and prepare user data header
  DayaBay::UserDataHeader* co60Data = MakeHeaderObject();
  Context context = readoutHdr->context();
  context.SetDetId(DetectorId::kAll);
  co60Data->setContext(context);
  co60Data->setEarliest(readoutHdr->earliest());
  co60Data->setLatest(readoutHdr->latest());

  // -- 3. Finally, put data into the user data header
  co60Data->set("e_co60",e_co60); 
  co60Data->set("x_co60",x_co60); 
  co60Data->set("y_co60",y_co60); 
  co60Data->set("z_co60",z_co60); 
  co60Data->set("nPESUM_co60",nPESUM); 
  co60Data->set("nPERMS_co60",nPERMS); 

//==============================
  return StatusCode::SUCCESS;
}