Ge68Data.cc

Go to the documentation of this file.

#include "Ge68Data.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"

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

  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("Ge68Location", m_Ge68Location="/Event/Tag/Calib/Ge68",
                  "Location of the Ge68 Tag");
}

Ge68Data::~Ge68Data() {
}

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

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

StatusCode Ge68Data::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 Ge68, if not, skip.
  bool isGe68 = exist<DayaBay::HeaderObject>(evtSvc(), m_Ge68Location);
  if (!isGe68){
    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_ge68 = rectrigger.energy()/Gaudi::Units::MeV;
  float x_ge68 = rectrigger.position().x() / Gaudi::Units::mm;
  float y_ge68 = rectrigger.position().y() / Gaudi::Units::mm;
  float z_ge68 = 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* ge68Data = MakeHeaderObject();
  Context context = readoutHdr->context();
  context.SetDetId(DetectorId::kAll);
  ge68Data->setContext(context);
  ge68Data->setEarliest(readoutHdr->earliest());
  ge68Data->setLatest(readoutHdr->latest());

  // -- 3. Finally, put data into the user data header
  ge68Data->set("e_ge68",e_ge68); 
  ge68Data->set("x_ge68",x_ge68); 
  ge68Data->set("y_ge68",y_ge68); 
  ge68Data->set("z_ge68",z_ge68); 
  ge68Data->set("nPESUM_ge68",nPESUM); 
  ge68Data->set("nPERMS_ge68",nPERMS); 

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