G4_Target.C

Go to the documentation of this file.

using namespace std;

void SetupTarget(
  PHG4Reco *g4Reco,
        const int use_g4steps,
        const double target_l,
        const double target_z
) {

  const double target_coil_pos_z = -130;
  const double collimator_pos_z = target_coil_pos_z - 302.36;

  //gSystem->Load("libfun4all.so");
        //gSystem->Load("libg4detectors");
        //gSystem->Load("libg4testbench");

  //Fun4AllServer *se = Fun4AllServer::instance();
  //PHG4Reco *g4Reco = (PHG4Reco *) se->getSubsysReco("PHG4RECO");

        PHG4CollimatorSubsystem* collimator = new PHG4CollimatorSubsystem("Collimator",0);
        collimator->SuperDetector("Collimator");
        collimator->set_double_param("place_z", collimator_pos_z);//-302.36 cm
        collimator->set_double_param("size_z",121.92);
        collimator->SetActive(1);
        g4Reco->registerSubsystem(collimator);

        PHG4TargetCoilSubsystem* coil_0 = new PHG4TargetCoilSubsystem("Coil", 0);
        coil_0->SuperDetector("Coil");
        coil_0->set_double_param("rot_x", 90.);
        coil_0->set_double_param("rot_y", 0.);
        coil_0->set_double_param("rot_z", 0.);
        coil_0->set_double_param("place_x", 0.);
        coil_0->set_double_param("place_y", (22.7+4.)/2);
        coil_0->set_double_param("place_z", target_coil_pos_z);
        coil_0->set_int_param("use_g4steps", use_g4steps);
        coil_0->SetActive(1);                                   // it is an active volume - save G4Hits
        g4Reco->registerSubsystem(coil_0);

        PHG4TargetCoilSubsystem* coil_1 = new PHG4TargetCoilSubsystem("Coil", 1);
        coil_1->SuperDetector("Coil");
        coil_1->set_double_param("rot_x", -90.);
        coil_1->set_double_param("rot_y", 0.);
        coil_1->set_double_param("rot_z", 0.);
        coil_1->set_double_param("place_x", 0.);
        coil_1->set_double_param("place_y", -(22.7+4.)/2);
        coil_1->set_double_param("place_z", target_coil_pos_z);
        coil_1->set_int_param("use_g4steps", use_g4steps);
        coil_1->SetActive(1);                                   // it is an active volume - save G4Hits
        g4Reco->registerSubsystem(coil_1);

        PHG4CylinderSubsystem* target = new PHG4CylinderSubsystem("Target", 0);
        target->SuperDetector("Target");
        target->set_double_param("length", target_l);
        target->set_double_param("rot_x", 0.);
        target->set_double_param("rot_y", 0.);
        target->set_double_param("rot_z", 0.);
        target->set_double_param("place_x", 0.);
        target->set_double_param("place_y", 0.);
        target->set_double_param("place_z", target_coil_pos_z + target_z);
        target->set_double_param("radius", 0.);
        target->set_double_param("thickness", (2.)/2);
        target->set_string_param("material", "Target");          // material of target
        target->set_int_param("lengthviarapidity", 0);
        target->set_int_param("use_g4steps", use_g4steps);
        target->SetActive(1);                                   // it is an active volume - save G4Hits
        g4Reco->registerSubsystem(target);

  return;
}