void HCALInner_Towers(int verbosity = 0) {
gSystem->Load("libg4calo.so");
gSystem->Load("libcalo_reco.so");
Fun4AllServer *se = Fun4AllServer::instance();
HcalRawTowerBuilder *TowerBuilder = new HcalRawTowerBuilder("HcalInRawTowerBuilder");
TowerBuilder->Detector("HCALIN");
TowerBuilder->set_sim_tower_node_prefix("SIM");
TowerBuilder->Verbosity(verbosity);
se->registerSubsystem( TowerBuilder );
// From 2016 Test beam sim
RawTowerDigitizer *TowerDigitizer = new RawTowerDigitizer("HcalInRawTowerDigitizer");
TowerDigitizer->Detector("HCALIN");
// TowerDigitizer->set_raw_tower_node_prefix("RAW_LG");
TowerDigitizer->set_digi_algorithm(
RawTowerDigitizer::kSimple_photon_digitalization);
TowerDigitizer->set_pedstal_central_ADC(0);
TowerDigitizer->set_pedstal_width_ADC(1); // From Jin's guess. No EMCal High Gain data yet! TODO: update
TowerDigitizer->set_photonelec_ADC(32. / 5.);
TowerDigitizer->set_photonelec_yield_visible_GeV(32. / 5 / (0.4e-3));
TowerDigitizer->set_zero_suppression_ADC(-0); // no-zero suppression
se->registerSubsystem(TowerDigitizer);
//Default sampling fraction for SS310
double visible_sample_fraction_HCALIN = 0.0631283 ; //, /gpfs/mnt/gpfs04/sphenix/user/jinhuang/prod_analysis/hadron_shower_res_nightly/./G4Hits_sPHENIX_pi-_eta0_16GeV-0000.root_qa.rootQA_Draw_HCALIN_G4Hit.pdf
if(inner_hcal_material_Al) visible_sample_fraction_HCALIN = 0.162166; //for "G4_Al", Abhisek Sen <*****@*****.**>
RawTowerCalibration *TowerCalibration = new RawTowerCalibration("HcalInRawTowerCalibration");
TowerCalibration->Detector("HCALIN");
// TowerCalibration->set_raw_tower_node_prefix("RAW_LG");
// TowerCalibration->set_calib_tower_node_prefix("CALIB_LG");
TowerCalibration->set_calib_algorithm(RawTowerCalibration::kSimple_linear_calibration);
TowerCalibration->set_calib_const_GeV_ADC(0.4e-3 / visible_sample_fraction_HCALIN);
TowerCalibration->set_pedstal_ADC(0);
TowerCalibration->set_zero_suppression_GeV(-1); // no-zero suppression
se->registerSubsystem(TowerCalibration);
return;
}
void HCALOuter_Towers(int verbosity = 0) {
gSystem->Load("libfun4all.so");
gSystem->Load("libg4detectors.so");
Fun4AllServer *se = Fun4AllServer::instance();
HcalRawTowerBuilder* TowerBuilder = new HcalRawTowerBuilder("HcalOutRawTowerBuilder");
TowerBuilder->Detector("HCALOUT");
TowerBuilder->set_sim_tower_node_prefix("SIM");
TowerBuilder->Verbosity(verbosity);
se->registerSubsystem( TowerBuilder );
// From 2016 Test beam sim
RawTowerDigitizer *TowerDigitizer = new RawTowerDigitizer("HcalOutRawTowerDigitizer");
TowerDigitizer->Detector("HCALOUT");
// TowerDigitizer->set_raw_tower_node_prefix("RAW_LG");
TowerDigitizer->set_digi_algorithm(
RawTowerDigitizer::kSimple_photon_digitalization);
TowerDigitizer->set_pedstal_central_ADC(0);
TowerDigitizer->set_pedstal_width_ADC(1); // From Jin's guess. No EMCal High Gain data yet! TODO: update
TowerDigitizer->set_photonelec_ADC(16. / 5.);
TowerDigitizer->set_photonelec_yield_visible_GeV(16. / 5 / (0.2e-3));
TowerDigitizer->set_zero_suppression_ADC(-0); // no-zero suppression
se->registerSubsystem(TowerDigitizer);
const double visible_sample_fraction_HCALOUT = 3.38021e-02; // /gpfs/mnt/gpfs04/sphenix/user/jinhuang/prod_analysis/hadron_shower_res_nightly/./G4Hits_sPHENIX_pi-_eta0_16GeV.root_qa.rootQA_Draw_HCALOUT_G4Hit.pdf
RawTowerCalibration *TowerCalibration = new RawTowerCalibration("HcalOutRawTowerCalibration");
TowerCalibration->Detector("HCALOUT");
// TowerCalibration->set_raw_tower_node_prefix("RAW_LG");
// TowerCalibration->set_calib_tower_node_prefix("CALIB_LG");
TowerCalibration->set_calib_algorithm(RawTowerCalibration::kSimple_linear_calibration);
TowerCalibration->set_calib_const_GeV_ADC(0.2e-3 / visible_sample_fraction_HCALOUT);
TowerCalibration->set_pedstal_ADC(0);
TowerCalibration->set_zero_suppression_GeV(-1); // no-zero suppression
se->registerSubsystem(TowerCalibration);
return;
}
