G4_SensitiveDetectors.C

Go to the documentation of this file.

#pragma once
#if ROOT_VERSION_CODE >= ROOT_VERSION(6,00,0)
#include <g4detectors/PHG4BlockSubsystem.h>
#include <fun4all/Fun4AllServer.h>
#include <TSQLServer.h>
#include <TSQLResult.h>
#include <TSQLRow.h>
class SubsysReco;
R__LOAD_LIBRARY(libg4detectors)
#endif

#define LogDebug(exp)   std::cout<<"DEBUG: "  <<__FILE__<<": "<<__LINE__<<": "<< exp << std::endl

using namespace std;

void SetupSensitiveDetectors(
                PHG4Reco *g4Reco,
                const int verbosity = 0
                ){

        GeomSvc *geom_svc = GeomSvc::instance();

        vector<string> sim_list = geom_svc->getDefaultSimList();

        for(int i=0; i<sim_list.size(); ++i){
                string name = sim_list[i];

                double size[3];
                double place[3];
                double rot[3];

                int id = geom_svc->getDetectorID(name);
                place[0] = geom_svc->getDetectorX0(name);
                place[1] = geom_svc->getDetectorY0(name);
                place[2] = geom_svc->getDetectorZ0(name);
                size[0] = geom_svc->getPlaneScaleX(id);
                size[1] = geom_svc->getPlaneScaleY(id);
                size[2] = geom_svc->getPlaneScaleZ(id);
                string material = geom_svc->getPlaneMaterial(id);

                if(verbosity > 2) {
                        LogDebug("");
                        cout
                                << "name: " << name
                                << ", id: " << id
                                << " {" << size[0] << ", " << size[1] << ", " << size[2] << "} "
                                << " {" << place[0] << ", " << place[1] << ", " << place[2] << "} "
                                << " {" << rot[0] << ", " << rot[1] << ", " << rot[2] << "} "
                                << material
                                << endl;
                }

                PHG4BlockSubsystem *box = new PHG4BlockSubsystem(name.c_str(), 0);
                box->SuperDetector(name.c_str());
                box->set_double_param("size_x", size[0]);
                box->set_double_param("size_y", size[1]);
                box->set_double_param("size_z", size[2]);
                box->set_double_param("place_x", place[0]);
                box->set_double_param("place_y", place[1]);
                box->set_double_param("place_z", place[2]);
                box->set_string_param("material", material.c_str());// G4_Si, G4_AIR, G4_Galactic
                //box->set_string_param("material", "G4_Galactic");// G4_Si, G4_AIR, G4_Galactic
                box->SetActive(1);
                g4Reco->registerSubsystem(box);
        }

        return;
}