ADLEDData.cc

Go to the documentation of this file.

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

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

  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("ADLEDLocation", m_ADLEDLocation="/Event/Tag/Calib/ADLED",
                  "Location of the ADLED Tag");
}

ADLEDData::~ADLEDData() {
}

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

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

StatusCode ADLEDData::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 ADLED, if not, skip.
  bool isADLED = exist<DayaBay::HeaderObject>(evtSvc(), m_ADLEDLocation);
  if (!isADLED){
    return StatusCode::SUCCESS;
  }

//  // Skip SAB Data 
//  Site::Site_t siteId = readoutHdr->context().GetSite();
//  
//  if (siteId == Site::kSAB) {
//    return StatusCode::SUCCESS;
//  }

  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 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_adled = rectrigger.energy()/Gaudi::Units::MeV;
  float x_adled = rectrigger.position().x() / Gaudi::Units::mm;
  float y_adled = rectrigger.position().y() / Gaudi::Units::mm;
  float z_adled = 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* adledData = MakeHeaderObject();
  Context context = readoutHdr->context();
  context.SetDetId(DetectorId::kAll);
  adledData->setContext(context);
  adledData->setEarliest(readoutHdr->earliest());
  adledData->setLatest(readoutHdr->latest());

  // -- 3. Finally, put data into the user data header
  adledData->set("e_adled",e_adled); 
  adledData->set("x_adled",x_adled); 
  adledData->set("y_adled",y_adled); 
  adledData->set("z_adled",z_adled); 
  adledData->set("nPESUM_adled",nPESUM); 
  adledData->set("nPERMS_adled",nPERMS); 

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