macros/diskinematics_fun4all/G4_FastSim_Reco_EIC.C

0001 using namespace std;
0002
0003 void FastSim_Reco_EIC(int verbosity = 0) {
0004
0005   //---------------
0006   // Load libraries
0007   //---------------
0008
0009   gSystem->Load("libfun4all.so");
0010   gSystem->Load("libg4hough.so");
0011
0012   //---------------
0013   // Fun4All server
0014   //---------------
0015
0016   Fun4AllServer *se = Fun4AllServer::instance();
0017
0018   PHG4TrackFastSim* kalman = new PHG4TrackFastSim("PHG4TrackFastSim");
0019   kalman->Verbosity(0);
0020
0021   kalman->set_use_vertex_in_fitting(true);
0022   kalman->set_vertex_xy_resolution(50E-4);
0023   kalman->set_vertex_z_resolution(50E-4);
0024
0025   kalman->set_detector_type(PHG4TrackFastSim::Vertical_Plane); // Vertical_Plane, Cylinder
0026   kalman->set_phi_resolution(50E-4);
0027   kalman->set_r_resolution(1.);
0028
0029   kalman->set_pat_rec_hit_finding_eff(1.);
0030   kalman->set_pat_rec_noise_prob(0.);
0031
0032   std::string phg4hits_names[] = {"G4HIT_SVTX", "G4HIT_SILICON_TRACKER", "G4HIT_MAPS", "G4HIT_FGEM_0", "G4HIT_FGEM_1", "G4HIT_FGEM_2", "G4HIT_FGEM_3", "G4HIT_FGEM_4", "G4HIT_EGEM_0", "G4HIT_EGEM_1", "G4HIT_EGEM_2"};
0033   kalman->set_phg4hits_names(phg4hits_names, 11);
0034   kalman->set_sub_top_node_name("SVTX");
0035   kalman->set_trackmap_out_name("SvtxTrackMap_FastSim");
0036
0037   // Saved track states (projections)
0038   //std::string state_names[] = {"CEMC", "HCALIN", "HCALOUT", "FEMC", "FHCAL", "EEMC"};
0039   std::string state_names[] = {"CEMC", "FEMC", "FHCAL", "EEMC"};
0040   kalman->set_state_names(state_names, 4);
0041
0042   kalman->set_fit_alg_name("KalmanFitterRefTrack");//
0043   kalman->set_primary_assumption_pid(-211); // -211 = pi-
0044   kalman->set_do_evt_display(false);
0045
0046   se->registerSubsystem(kalman);
0047
0048 }
0049