BeamSpotMonitorModule.cc

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#include <reconstruction/modules/ipMonitor/BeamSpotMonitorModule.h>
#include <framework/datastore/StoreObjPtr.h>
#include <framework/dataobjects/EventMetaData.h>
 
#include <TMatrixDSym.h>
 
using namespace Belle2;
 
//-----------------------------------------------------------------
//                 Register the Module
//-----------------------------------------------------------------
REG_MODULE(BeamSpotMonitor);
 
//-----------------------------------------------------------------
//                 Implementation
//-----------------------------------------------------------------
 
BeamSpotMonitorModule::BeamSpotMonitorModule() : Module()
{
  // Set module properties
  setDescription("Module for the monitoring of the BeamSpot position and size");
 
  // Parameter definitions
  addParam("outputFileName", m_rootFileName, "Name of output root file.", std::string("BeamSpotMonitor.root"));
}
 
void BeamSpotMonitorModule::initialize()
{
  m_EventMetaData.isRequired();
 
  TDirectory* olddir = gDirectory;
  m_rootFilePtr = new TFile(m_rootFileName.c_str(), "RECREATE");
 
  //tree initialization
  m_tree = new TTree("bspt", "RECREATE");
  b_exp = m_tree->Branch("exp", &m_exp, "exp/i");
  b_run = m_tree->Branch("run", &m_run, "run/i");
  b_x = m_tree->Branch("x", &m_x, "x/d");
  b_y = m_tree->Branch("y", &m_y, "y/d");
  b_z = m_tree->Branch("z", &m_z, "z/d");
  b_xErr = m_tree->Branch("xErr", &m_xErr, "xErr/d");
  b_yErr = m_tree->Branch("yErr", &m_yErr, "yErr/d");
  b_zErr = m_tree->Branch("zErr", &m_zErr, "zErr/d");
  b_xSize = m_tree->Branch("xSize", &m_xSize, "xSize/d");
  b_ySize = m_tree->Branch("ySize", &m_ySize, "ySize/d");
  b_zSize = m_tree->Branch("zSize", &m_zSize, "zSize/d");
 
  olddir->cd();
 
}
 
void BeamSpotMonitorModule::beginRun()
{
  if (! m_BeamSpotDB.isValid()) {
    B2WARNING("No valid BeamSpot for the requested IoV");
  } else {
    m_BeamSpot = *m_BeamSpotDB;
  }
 
}
 
void BeamSpotMonitorModule::event()
{
  m_exp = m_EventMetaData->getExperiment();
  m_run = m_EventMetaData->getRun();
  B2DEBUG(25, "monitoring beam spot for experiment = " << m_exp << ", run = " << m_run);
 
  if (! m_BeamSpotDB.isValid())
    return;
 
  //retrieve vertex position
  m_x = m_BeamSpot.getIPPosition().X();
  m_y = m_BeamSpot.getIPPosition().Y();
  m_z = m_BeamSpot.getIPPosition().Z();
 
  //retrieve vertex position error
  TMatrixDSym posErr = m_BeamSpot.getIPPositionCovMatrix();
  m_xErr = sqrt(posErr[0][0]);
  m_yErr = sqrt(posErr[1][1]);
  m_zErr = sqrt(posErr[2][2]);
 
  //retrieve beam spot size
  TMatrixDSym size = m_BeamSpot.getSizeCovMatrix();
  m_xSize = sqrt(size[0][0]);
  m_ySize = sqrt(size[1][1]);
  m_zSize = sqrt(size[2][2]);
 
  m_tree->Fill();
 
}
 
void BeamSpotMonitorModule::terminate()
{
  TDirectory* olddir = gDirectory;
  if (m_rootFilePtr != nullptr) {
 
    m_rootFilePtr->cd();
 
    //write the tree
    m_tree->Write();
 
    m_rootFilePtr->Close();
 
    olddir->cd();
  }
}