G4_DSTReader_JLeic.C

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4_DSTReader_JLeic.C

#ifndef MACRO_G4DSTREADERJLEIC_C
#define MACRO_G4DSTREADERJLEIC_C

#include <GlobalVariables.C>

#include <G4_Barrel_Hcal_JLeic.C>
#include <G4_BeamLine_JLeic.C>
#include <G4_BlackHole.C>
#include <G4_CTD_JLeic.C>
#include <G4_DIRC_JLeic.C>
#include <G4_DRich_JLeic.C>
#include <G4_EndCap_Electron_JLeic.C>
#include <G4_EndCap_Hadron_JLeic.C>
#include <G4_Gem_JLeic.C>
#include <G4_Magnet_Cleo.C>
#include <G4_Pipe_EIC.C>
#include <G4_VTX_JLeic.C>

#include <g4eval/PHG4DSTReader.h>

#include <fun4all/Fun4AllServer.h>

R__LOAD_LIBRARY(libg4eval.so)




namespace Enable
{
  bool DSTREADER = false;
  int DSTREADER_VERBOSITY = 0;
}  // namespace Enable

namespace G4DSTREADER
{
  bool save_g4_raw = true;
}  // namespace G4DSTREADER

void G4DSTreader(const string &outputFile = "G4JLeic.root")
{
  int verbosity = max(Enable::VERBOSITY, Enable::DSTREADER_VERBOSITY);

  // save a comprehensive  evaluation file
  PHG4DSTReader *ana = new PHG4DSTReader(outputFile);
  ana->set_save_particle(true);
  ana->set_load_all_particle(false);
  ana->set_load_active_particle(true);
  ana->set_save_vertex(true);

  ana->Verbosity(verbosity);

  if (G4DSTREADER::save_g4_raw)
  {
    if (Enable::PIPE)
    {
      if (Enable::ABSORBER || Enable::PIPE_ABSORBER)
        ana->AddNode("PIPE");
    }
    if (Enable::VTX)
    {
      ana->AddNode("JLVTX");
    }
    if (Enable::CTD)
    {
      ana->AddNode("JLCTD");
    }
    if (Enable::DIRC)
    {
      ana->AddNode("JLDIRC");
    }

    if (Enable::MAGNET)
    {
      if (Enable::ABSORBER || Enable::MAGNET_ABSORBER)
        ana->AddNode("MAGNET");
    }

    if (Enable::BARREL_HCAL)
    {
      ana->AddNode("BARRELHCAL");
    }

    if (Enable::GEM)
    {
      ana->AddNode("GEMHADRON");
      ana->AddNode("GEMELECTRON");
    }

    if (Enable::DRICH)
    {
      ana->AddNode("DRICH");
    }

    if (Enable::ENDCAP_ELECTRON)
    {
      ana->AddNode("ECELECTRON");
    }

    if (Enable::ENDCAP_HADRON)
    {
      ana->AddNode("ECHADRON");
    }
    if (Enable::BEAMLINE)
    {
      ana->AddNode("BEAMLINE");
      if (Enable::BEAMLINE_ABSORBER)
      {
        ana->AddNode("ABSORBER_BEAMLINE");
      }
    }
    if (Enable::BLACKHOLE)
    {
      ana->AddNode("BH_1");
      ana->AddNode("BH_FORWARD_PLUS");
      ana->AddNode("BH_FORWARD_NEG");
    }
  }
  Fun4AllServer *se = Fun4AllServer::instance();
  se->registerSubsystem(ana);
}
#endif  // MACRO_G4DSTREADERJLEIC_C