G4_mRICH.C

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

#pragma once

#include "GlobalVariables.C"

#include <g4main/PHG4Reco.h>

#include <g4mrich/PHG4mRICHSubsystem.h>
#include <g4trackfastsim/PHG4TrackFastSim.h>
//#include <g4rich/PHG4mRICHSubsystem.h>
namespace Enable
{
  bool mRICH = false;
  bool mRICH_OVERLAPCHECK = false;
}  // namespace Enable

void mRICHInit()
{
}

//-1: single module
// 0: build h-side sectors and e-side wall
// 1: build h-side sectors
// 2: build e-side wall
// 3: build h-side wall
// 4: build h-side wall and e-side wall
// 5: build b-side sectors
// 6: build projected e-side wall

void mRICHSetup(PHG4Reco *g4Reco, const int detectorSetup = 1,  //1: full setup; 0:skeleton
                const int mRICHsystemSetup = 6                  //5//2//-1//
)
{
  bool OverlapCheck = Enable::OVERLAPCHECK || Enable::mRICH_OVERLAPCHECK;

  PHG4mRICHSubsystem *mRICH = new PHG4mRICHSubsystem("mRICH", 1);
  mRICH->set_int_param("detectorSetup", detectorSetup);
  mRICH->set_int_param("subsystemSetup", mRICHsystemSetup);
  mRICH->UseCalibFiles(PHG4DetectorSubsystem::xml);
  mRICH->SetCalibrationFileDir(string(getenv("CALIBRATIONROOT")) + string("/mRICH/Geometry/"));
  mRICH->OverlapCheck(OverlapCheck);
  //  mRICH->Verbosity(5);

  g4Reco->registerSubsystem(mRICH);

  /*if (TRACKING::FastKalmanFilter)
  {
    // project to an reference plane at z=170 cm
    TRACKING::FastKalmanFilter-> add_zplane_state("mRICH", -170);
    TRACKING::ProjectionNames.insert("mRICH");
  }*/
}

void mRICH_Eval(std::string outputfile, int verbosity = 0)
{
  gSystem->Load("libfun4all.so");
  gSystem->Load("libg4eval.so");
  Fun4AllServer *se = Fun4AllServer::instance();

  return;
}