In This Package:

AdRecHistoCheck.cc

Go to the documentation of this file.

#include "AdRecHistoCheck.h"

#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "GaudiKernel/DataObject.h"
#include "StatisticsSvc/IStatisticsSvc.h"

#include "Event/RecHeader.h"
#include "Event/RecTrigger.h"

#include "TH1F.h"

using namespace DayaBay;
using namespace CLHEP;

AdRecHistoCheck::AdRecHistoCheck(const std::string& name,
                                 ISvcLocator* pSvcLocator) : 
  GaudiAlgorithm(name,pSvcLocator)
{
}

StatusCode AdRecHistoCheck::initialize()
{
  if ( service("StatisticsSvc", p_svc).isFailure() ) {
    error()<<" No StatisticsSvc available ! "<<endreq;
    return StatusCode::FAILURE;
  } 
  
  TH1F*  h1 = new TH1F("recX","recX", 50, -1000., 1000. );
  h1->GetXaxis()->SetTitle("mm");
  if ( p_svc->put("/File1/AdRec/recX", h1).isFailure() ) {
    error()<<"Could't register recX "<<endreq; 
  }
  
  TH1F*  h2 = new TH1F("recY","recY", 50, -1000., 1000. );
  h2->GetXaxis()->SetTitle("mm");
  if ( p_svc->put("/File1/AdRec/recY", h2).isFailure() ) {
    error()<<"Could't register recY "<<endreq; 
  }

  TH1F*  h3 = new TH1F("recZ","recZ", 50, -1000., 1000. );
  h3->GetXaxis()->SetTitle("mm");
  if ( p_svc->put("/File1/AdRec/recZ", h3).isFailure() ) {
    error()<<"Could't register recZ "<<endreq; 
  }

  info() << "Finished booking Histograms" << endreq;
  return StatusCode::SUCCESS;
}

StatusCode AdRecHistoCheck::execute()
{
  
  TH1* h1 = p_svc->getTH1F("/File1/AdRec/recX");
  TH1* h2 = p_svc->getTH1F("/File1/AdRec/recX");
  TH1* h3 = p_svc->getTH1F("/File1/AdRec/recX");
  if (h1 && h2 && h3 ) {
    
    DayaBay::RecHeader* recHeader 
      = get<DayaBay::RecHeader>(DayaBay::RecHeaderLocation::Default);
    
    if(!recHeader) {
      error()<< "Could not find recHeader" << endreq;
      return StatusCode::FAILURE; 
    }else {

      const DayaBay::RecTrigger* recTrigger = &(recHeader->recTrigger());
      if(!recTrigger) {
        info() << "No recTrigger in recHeader this cycle." << endreq;
        return StatusCode::SUCCESS;
      }
      
        // Get reconstructed pos 
      CLHEP::HepLorentzVector pos = recTrigger->position();

      h1->Fill(pos.x());
      h2->Fill(pos.y());
      h3->Fill(pos.z());

    }
  } else {
    error() <<" Could not retrieve Histograms "<<endreq;
    return StatusCode::FAILURE;
  } 
  
  return StatusCode::SUCCESS;
}

StatusCode AdRecHistoCheck::finalize()
{
  debug() << "Finalizing..." << endreq;
  
  return StatusCode::SUCCESS;
}

---