//! for visualization purpose only
double
CEmc_Vis(PHG4Reco* g4Reco, double radius, const int crossings, const int absorberactive = 0)
{
double emc_inner_radius = 95.; // emc inner radius from engineering drawing
double cemcthickness = 12.7;
double emc_outer_radius = emc_inner_radius + cemcthickness; // outer radius
if (radius > emc_inner_radius) {
cout << "inconsistency: preshower radius+thickness: " << radius
<< " larger than emc inner radius: " << emc_inner_radius
<< endl;
gSystem->Exit(-1);
}
//---------------
// Load libraries
//---------------
gSystem->Load("libg4detectors.so");
gSystem->Load("libg4testbench.so");
// the radii are only to determined the thickness of the cemc
radius = emc_inner_radius;
int ilayer = Min_cemc_layer;
PHG4SpacalSubsystem *cemc;
cemc = new PHG4SpacalSubsystem("CEMC", ilayer);
cemc ->get_geom().set_radius(emc_inner_radius);
cemc ->get_geom().set_thickness(cemcthickness);
cemc ->get_geom().set_virualize_fiber(true);
cemc ->get_geom().set_construction_verbose(2);
cemc->SetActive();
cemc->SuperDetector("CEMC");
if (absorberactive) cemc->SetAbsorberActive();
cemc->OverlapCheck(overlapcheck);
g4Reco->registerSubsystem( cemc );
if (ilayer > Max_cemc_layer)
{
cout << "layer discrepancy, current layer " << ilayer
<< " max cemc layer: " << Max_cemc_layer << endl;
}
radius += cemcthickness;
radius += no_overlapp;
PHG4CylinderSubsystem *cyl = new PHG4CylinderSubsystem("EMCELECTRONICS", 0);
cyl->SetRadius(radius);
cyl->SetMaterial("G4_TEFLON"); // plastic
cyl->SetThickness(0.5);
if (absorberactive) cyl->SetActive();
g4Reco->registerSubsystem( cyl );
radius += 0.5;
radius += no_overlapp;
return radius;
}
//! EMCal setup macro - 1D azimuthal projective SPACAL
double
CEmc_1DProjectiveSpacal(PHG4Reco* g4Reco, double radius, const int crossings, const int absorberactive = 0)
{
double emc_inner_radius = 95.; // emc inner radius from engineering drawing
double cemcthickness = 12.7;
double emc_outer_radius = emc_inner_radius + cemcthickness; // outer radius
if (radius > emc_inner_radius) {
cout << "inconsistency: preshower radius+thickness: " << radius
<< " larger than emc inner radius: " << emc_inner_radius
<< endl;
gSystem->Exit(-1);
}
//---------------
// Load libraries
//---------------
gSystem->Load("libg4detectors.so");
gSystem->Load("libg4testbench.so");
// boundary check
if (radius> emc_inner_radius - 1.5 - no_overlapp)
{
cout <<"G4_CEmc_Spacal.C::CEmc() - expect radius < "<<emc_inner_radius - 1.5 - no_overlapp<<" to install SPACAL"<<endl;
exit(1);
}
radius = emc_inner_radius - 1.5 - no_overlapp;
// 1.5cm thick teflon as an approximation for EMCAl light collection + electronics (10% X0 total estimated)
PHG4CylinderSubsystem *cyl = new PHG4CylinderSubsystem("CEMC_ELECTRONICS", 0);
cyl->SuperDetector("CEMC_ELECTRONICS");
cyl->set_double_param("radius",radius);
cyl->set_string_param("material","G4_TEFLON");
cyl->set_double_param("thickness",1.5);
if (absorberactive) cyl->SetActive();
g4Reco->registerSubsystem( cyl );
radius += 1.5;
radius += no_overlapp;
int ilayer = Min_cemc_layer;
PHG4SpacalSubsystem *cemc;
cemc = new PHG4SpacalSubsystem("CEMC", ilayer);
cemc ->get_geom().set_radius(emc_inner_radius);
cemc ->get_geom().set_thickness(cemcthickness);
cemc ->get_geom().set_construction_verbose(1);
cemc->SetActive();
cemc->SuperDetector("CEMC");
if (absorberactive) cemc->SetAbsorberActive();
cemc->OverlapCheck(overlapcheck);
g4Reco->registerSubsystem( cemc );
if (ilayer > Max_cemc_layer)
{
cout << "layer discrepancy, current layer " << ilayer
<< " max cemc layer: " << Max_cemc_layer << endl;
}
radius += cemcthickness;
radius += no_overlapp;
// 0.5cm thick Stainless Steel as an approximation for EMCAl support system
cyl = new PHG4CylinderSubsystem("CEMC_SPT", 0);
cyl->SuperDetector("CEMC_SPT");
cyl->set_double_param("radius",radius);
cyl->set_string_param("material","SS310"); // SS310 Stainless Steel
cyl->set_double_param("thickness",0.5);
if (absorberactive)
cyl->SetActive();
g4Reco->registerSubsystem(cyl);
radius += 0.5;
radius += no_overlapp;
return radius;
}
//! 2D full projective SPACAL
double
CEmc_2DProjectiveSpacal(PHG4Reco* g4Reco, double radius, const int crossings,
const int absorberactive = 0)
{
double emc_inner_radius = 95.; // emc inner radius from engineering drawing
double cemcthickness = 21.00000 - no_overlapp;
double emc_outer_radius = emc_inner_radius + cemcthickness; // outer radius
if (radius > emc_inner_radius)
{
cout << "inconsistency: preshower radius+thickness: " << radius
<< " larger than emc inner radius: " << emc_inner_radius << endl;
gSystem->Exit(-1);
}
//---------------
// Load libraries
//---------------
gSystem->Load("libg4detectors.so");
// the radii are only to determined the thickness of the cemc
radius = emc_inner_radius;
//---------------
// Load libraries
//---------------
// 1.5cm thick teflon as an approximation for EMCAl light collection + electronics (10% X0 total estimated)
PHG4CylinderSubsystem *cyl = new PHG4CylinderSubsystem("CEMC_ELECTRONICS", 0);
cyl->set_double_param("radius",radius);
cyl->set_string_param("material","G4_TEFLON");
cyl->set_double_param("thickness",1.5- no_overlapp);
cyl->SuperDetector("CEMC_ELECTRONICS");
cyl->OverlapCheck(overlapcheck);
if (absorberactive) cyl->SetActive();
g4Reco->registerSubsystem( cyl );
radius += 1.5;
cemcthickness -= 1.5+no_overlapp;
// 0.5cm thick Stainless Steel as an approximation for EMCAl support system
cyl = new PHG4CylinderSubsystem("CEMC_SPT", 0);
cyl->SuperDetector("CEMC_SPT");
cyl->set_double_param("radius",radius +cemcthickness - 0.5 );
cyl->set_string_param("material","SS310"); // SS310 Stainless Steel
cyl->set_double_param("thickness",0.5 - no_overlapp);
cyl->OverlapCheck(overlapcheck);
if (absorberactive)
cyl->SetActive();
g4Reco->registerSubsystem(cyl);
cemcthickness -= 0.5+no_overlapp;
//---------------
// Load libraries
//---------------
int ilayer = Min_cemc_layer;
PHG4SpacalSubsystem *cemc;
cemc = new PHG4SpacalSubsystem("CEMC", ilayer);
cemc->get_geom().set_config(
PHG4CylinderGeom_Spacalv1::kFullProjective_2DTaper_SameLengthFiberPerTower);
// load from hard code map for now. Need to move to calibration file or database.
cemc->get_geom().load_demo_sector_tower_map3();
cemc->get_geom().set_radius(radius);
cemc->get_geom().set_thickness(cemcthickness);
cemc->get_geom().set_construction_verbose(1);
cemc->SetActive();
cemc->SuperDetector("CEMC");
if (absorberactive)
cemc->SetAbsorberActive();
cemc->OverlapCheck(overlapcheck);
g4Reco->registerSubsystem(cemc);
if (ilayer > Max_cemc_layer)
{
cout << "layer discrepancy, current layer " << ilayer
<< " max cemc layer: " << Max_cemc_layer << endl;
}
radius += cemcthickness;
radius += no_overlapp;
return radius;
}
//! test for 2D projective SPACAL using 1D projective modules. For evaluation only
double
CEmc_Proj(PHG4Reco* g4Reco, double radius, const int crossings, const int absorberactive = 0)
{
double emc_inner_radius = 95.; // emc inner radius from engineering drawing
double cemcthickness = 12.9+1.5;
double emc_outer_radius = emc_inner_radius + cemcthickness; // outer radius
if (radius > emc_inner_radius) {
cout << "inconsistency: preshower radius+thickness: " << radius
<< " larger than emc inner radius: " << emc_inner_radius
<< endl;
gSystem->Exit(-1);
}
//---------------
// Load libraries
//---------------
gSystem->Load("libg4detectors.so");
gSystem->Load("libg4testbench.so");
cout <<"Get_Spacal_Tilt() = "<<Get_Spacal_Tilt()<<endl;
// the radii are only to determined the thickness of the cemc
//double emc_inner_radius = radius; // emc inner radius from engineering spreadsheet
// double cemcthickness = 12.9+1.5;
// double emc_outer_radius = emc_inner_radius + cemcthickness; // outer radius
radius = emc_outer_radius;
int ilayer = Min_cemc_layer;
PHG4SpacalSubsystem *cemc;
cemc = new PHG4SpacalSubsystem("CEMC", ilayer);
cemc ->get_geom().set_radius(emc_inner_radius);
cemc ->get_geom().set_thickness(cemcthickness);
cemc ->get_geom().set_construction_verbose(2);
cemc ->get_geom().set_config(PHG4CylinderGeom_Spacalv1::kProjective_PolarTaper);
cemc ->get_geom().set_azimuthal_tilt(Get_Spacal_Tilt());
// cemc ->get_geom().set_azimuthal_tilt(-2.3315/2./95*1.9);
// cemc ->get_geom().set_azimuthal_tilt(6.28/256*2);
// cemc ->get_geom().set_azimuthal_tilt(-6.28/256*4);
// cemc ->get_geom().set_azimuthal_seg_visible(true);
// cemc ->get_geom().set_virualize_fiber(false);
// cemc ->get_geom().set_assembly_spacing(0.001);
cemc ->get_geom().set_polar_taper_ratio(1.128);
// cemc ->get_geom().set_polar_taper_ratio(1.123);
// cemc ->get_geom().set_polar_taper_ratio(1.117);
// cemc ->get_geom().set_absorber_mat("G4_AIR");
// cemc ->get_geom().set_azimuthal_n_sec(256/2);
cemc->SetActive();
cemc->SuperDetector("CEMC");
if (absorberactive) cemc->SetAbsorberActive();
cemc->OverlapCheck(overlapcheck);
g4Reco->registerSubsystem( cemc );
if (ilayer > Max_cemc_layer)
{
cout << "layer discrepancy, current layer " << ilayer
<< " max cemc layer: " << Max_cemc_layer << endl;
}
radius += cemcthickness;
radius += no_overlapp;
PHG4CylinderSubsystem *cyl = new PHG4CylinderSubsystem("EMCELECTRONICS", 0);
cyl->set_double_param("radius",radius);
cyl->set_string_param("material","G4_TEFLON");
cyl->set_double_param("thickness",0.5);
if (absorberactive) cyl->SetActive();
g4Reco->registerSubsystem( cyl );
radius += 0.5;
radius += no_overlapp;
return radius;
}
//! 2D full projective SPACAL
double
CEmc_2DProjectiveSpacal(PHG4Reco *g4Reco, double radius, const int crossings,
const int absorberactive = 0)
{
double emc_inner_radius = 92; // emc inner radius from engineering drawing
double cemcthickness = 24.00000 - no_overlapp;
//max radius is 116 cm;
double emc_outer_radius = emc_inner_radius + cemcthickness; // outer radius
assert(emc_outer_radius < 116);
if (radius > emc_inner_radius)
{
cout << "inconsistency: preshower radius+thickness: " << radius
<< " larger than emc inner radius: " << emc_inner_radius << endl;
gSystem->Exit(-1);
}
//---------------
// Load libraries
//---------------
gSystem->Load("libg4detectors.so");
// the radii are only to determined the thickness of the cemc
radius = emc_inner_radius;
//---------------
// Load libraries
//---------------
// 1.5cm thick teflon as an approximation for EMCAl light collection + electronics (10% X0 total estimated)
PHG4CylinderSubsystem *cyl = new PHG4CylinderSubsystem("CEMC_ELECTRONICS", 0);
cyl->set_double_param("radius", radius);
cyl->set_string_param("material", "G4_TEFLON");
cyl->set_double_param("thickness", 1.5 - no_overlapp);
cyl->SuperDetector("CEMC_ELECTRONICS");
cyl->OverlapCheck(overlapcheck);
if (absorberactive) cyl->SetActive();
g4Reco->registerSubsystem(cyl);
radius += 1.5;
cemcthickness -= 1.5 + no_overlapp;
// 0.5cm thick Stainless Steel as an approximation for EMCAl support system
cyl = new PHG4CylinderSubsystem("CEMC_SPT", 0);
cyl->SuperDetector("CEMC_SPT");
cyl->set_double_param("radius", radius + cemcthickness - 0.5);
cyl->set_string_param("material", "SS310"); // SS310 Stainless Steel
cyl->set_double_param("thickness", 0.5 - no_overlapp);
cyl->OverlapCheck(overlapcheck);
if (absorberactive)
cyl->SetActive();
g4Reco->registerSubsystem(cyl);
cemcthickness -= 0.5 + no_overlapp;
//---------------
// Load libraries
//---------------
int ilayer = 0;
PHG4SpacalSubsystem *cemc;
const bool use_2015_design = false;
if (use_2015_design)
{
cemc = new PHG4SpacalSubsystem("CEMC", ilayer);
cemc->set_int_param("config", PHG4CylinderGeom_Spacalv1::kFullProjective_2DTaper_SameLengthFiberPerTower);
cemc->set_double_param("radius", radius); // overwrite minimal radius
cemc->set_double_param("thickness", cemcthickness); // overwrite thickness
cemc->set_int_param("azimuthal_n_sec", 32);
// cemc->set_int_param("construction_verbose", 2);
cemc->SetActive();
cemc->SuperDetector("CEMC");
if (absorberactive)
cemc->SetAbsorberActive();
cemc->OverlapCheck(overlapcheck);
}
else
{
cemc = new PHG4SpacalSubsystem("CEMC", ilayer);
cemc->set_int_param("virualize_fiber", 0);
cemc->set_int_param("azimuthal_seg_visible", 1);
cemc->set_int_param("construction_verbose", 0);
cemc->Verbosity(0);
cemc->UseCalibFiles(PHG4DetectorSubsystem::xml);
cemc->SetCalibrationFileDir(string(getenv("CALIBRATIONROOT")) + string("/CEMC/Geometry_2017ProjTilted/"));
cemc->set_double_param("radius", radius); // overwrite minimal radius
cemc->set_double_param("thickness", cemcthickness); // overwrite thickness
cemc->SetActive();
cemc->SuperDetector("CEMC");
if (absorberactive)
cemc->SetAbsorberActive();
cemc->OverlapCheck(overlapcheck);
}
g4Reco->registerSubsystem(cemc);
if (ilayer > Max_cemc_layer)
{
cout << "layer discrepancy, current layer " << ilayer
<< " max cemc layer: " << Max_cemc_layer << endl;
}
radius += cemcthickness;
radius += no_overlapp;
return radius;
}
