SetupGenfitExtrapolationModule.cc

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2010 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Tobias Schlüter, Tadeas Bilka                            *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
//
// Common setup for everything that uses genfit's extrapolation code.
//
 
#include <tracking/modules/genfitUtilities/SetupGenfitExtrapolationModule.h>
#include <tracking/modules/genfitUtilities/Geant4MaterialInterface.h>
 
#include <geometry/GeometryManager.h>
 
#include <tracking/gfbfield/GFGeant4Field.h>
#include <genfit/FieldManager.h>
#include <genfit/MaterialEffects.h>
#include <genfit/TGeoMaterialInterface.h>
#include <genfit/IO.h>
 
#include <boost/iostreams/stream_buffer.hpp>
#include <boost/iostreams/concepts.hpp>
 
#include <TGeoManager.h>
 
using namespace Belle2;
 
REG_MODULE(SetupGenfitExtrapolation)
 
namespace {
  template<size_t T_level>
  class genfitSink : public boost::iostreams::sink {
  public:
    std::streamsize write(const char* s, std::streamsize n)
    {
      B2DEBUG(T_level, s);
      return n;
    }
  };
 
  genfitSink<200> debugSink;
  boost::iostreams::stream_buffer<genfitSink<200> > debugStreamBuf(debugSink);
  genfitSink<100> errorSink;
  boost::iostreams::stream_buffer<genfitSink<100> > errorStreamBuf(errorSink);
  genfitSink<150> printSink;
  boost::iostreams::stream_buffer<genfitSink<150> > printStreamBuf(printSink);
 
  void setupGenfitStreams()
  {
    genfit::debugOut.rdbuf(&debugStreamBuf);
    genfit::errorOut.rdbuf(&errorStreamBuf);
    genfit::printOut.rdbuf(&printStreamBuf);
  }
}
 
SetupGenfitExtrapolationModule::SetupGenfitExtrapolationModule() :
  Module(), m_vxdAlignment()
{
 
  setDescription("Sets up material handling for genfit extrapolation.  Also setups up I/O streams for"
                 " genfit in order to integrate it into basf2 logging system. Also sets up update of VXDAlignment from DB (temporary).");
  setPropertyFlags(c_ParallelProcessingCertified);
 
  addParam("ignoreIfPresent", m_ignoreIfPresent,
           "If true this module will silently ignore if the geometry is already "
           "present and do nothing in that case. If false a B2FATAL will be "
           "if the geometry was already created before", m_ignoreIfPresent);
 
  //input
  addParam("whichGeometry", m_geometry,
           "Which geometry should be used, either 'TGeo' or 'Geant4'", m_geometry);
 
  // Energy loss, multiple scattering configuration.
  addParam("energyLossBetheBloch", m_energyLossBetheBloch,
           "activate the material effect: EnergyLossBetheBloch", m_energyLossBetheBloch);
  addParam("noiseBetheBloch", m_noiseBetheBloch,
           "activate the material effect: NoiseBetheBloch", m_noiseBetheBloch);
  addParam("noiseCoulomb", m_noiseCoulomb,
           "activate the material effect: NoiseCoulomb", m_noiseCoulomb);
  addParam("energyLossBrems", m_energyLossBrems,
           "activate the material effect: EnergyLossBrems", m_energyLossBrems);
  addParam("noiseBrems", m_noiseBrems,
           "activate the material effect: NoiseBrems", m_noiseBrems);
  addParam("noEffects", m_noEffects,
           "switch off all material effects in Genfit. This overwrites all "
           "individual material effects switches", m_noEffects);
  addParam("MSCModel", m_mscModel,
           "Multiple scattering model", m_mscModel);
  addParam("useVXDAlignment", m_useVXDAlignment,
           "Use VXD alignment from database?", m_useVXDAlignment);
}
 
void SetupGenfitExtrapolationModule::initialize()
{
  if (genfit::FieldManager::getInstance()->isInitialized() or genfit::MaterialEffects::getInstance()->isInitialized()) {
    if (m_ignoreIfPresent) {
      B2DEBUG(50, "Magnetic field or material handling already initialized.  Not touching settings.");
      return;
    } else {
      B2FATAL("Magnetic field or material handling already initialized.  Not touching settings.");
    }
  }
 
  setupGenfitStreams();
 
  genfit::FieldManager::getInstance()->init(new GFGeant4Field());
  genfit::FieldManager::getInstance()->useCache();
 
  if (!geometry::GeometryManager::getInstance().getTopVolume()) {
    B2FATAL("No geometry set up so far. Load the geometry module.");
  }
 
  if (m_geometry == "TGeo") {
    if (!gGeoManager) {
      B2INFO("Building TGeo representation.");
      geometry::GeometryManager& geoManager = geometry::GeometryManager::getInstance();
      geoManager.createTGeoRepresentation();
    }
    genfit::MaterialEffects::getInstance()->init(new genfit::TGeoMaterialInterface());
  } else if (m_geometry == "Geant4") {
    genfit::MaterialEffects::getInstance()->init(new Geant4MaterialInterface());
  } else {
    B2FATAL("Invalid choice of geometry interface.  Please use 'TGeo' or 'Geant4'.");
  }
 
  // activate / deactivate material effects in genfit
  if (m_noEffects) {
    genfit::MaterialEffects::getInstance()->setNoEffects(true);
  } else {
    genfit::MaterialEffects::getInstance()->setEnergyLossBetheBloch(m_energyLossBetheBloch);
    genfit::MaterialEffects::getInstance()->setNoiseBetheBloch(m_noiseBetheBloch);
    genfit::MaterialEffects::getInstance()->setNoiseCoulomb(m_noiseCoulomb);
    genfit::MaterialEffects::getInstance()->setEnergyLossBrems(m_energyLossBrems);
    genfit::MaterialEffects::getInstance()->setNoiseBrems(m_noiseBrems);
    genfit::MaterialEffects::getInstance()->setMscModel(m_mscModel);
  }
}