//
// $Id: SiMaterialGeometry_run2b.rcp,v 1.1 2006/02/09 18:46:47 shfu Exp $
//
// File: SiMaterialGeometry.rcp
// Purpose: Define the parameters needed to build the default SMT material
// Created: ??-Nov-1997   John Hobbs
//
// May 20, 2005 by Shaohua Fu
// added Layer0 material // SFL0_May_2005
// for now we only have the wafer material for Layer0
//
// June 28, 2005 by Mike Hildreth
// added segmentation to SMT Support trough to mock up SMT cable plant
//
// 12th July 2005 by Michiel Sanders
// added correct densities and fractional radiation lengths for 
// SMT support trough
//
//
// $Revision: 1.1 $
//

float thickness = 0.03     // cm
string material_name = "SMT_Wafer"
string material_symbol = "SmtSi"
float material_zeff = 14.0
float material_aeff = 28.085
float material_rho = 2.33  // gm/cm^3

// factors to multiply the default SMT material for better description of 
// multiple scattering
float MS_Factor_ladder = 2.5
float MS_Factor_Fdisk = 2.5 
float MS_Factor_Hdisk = 1.

// Dummy values now
float trough_thickness = 0.8 // cm
string trough_name = "SMT_support_trough"
string trough_symbol = "SMTst"
float trough_zeff = 6.0
float trough_aeff = 12.001
float trough_rho = 1.6 // gm/cm^3

// Dummy values now
float bulkhead_thickness = 0.47625 // cm
string bulkhead_name = "SMT_Bulkhead"
string bulkhead_symbol = "SMTbh"
float bulkhead_zeff = 13.0
float bulkhead_aeff = 26.9815
float bulkhead_rho = 2.7 // gm/cm^3

// The geometrical information for the bulkheads

float bulkhead_zoffs = -2.33 //cm 
float bulkhead_phioffs = 0.0
int bulkhead_layer_count = 4
int bulkhead_sides1 = 6    // Layer 1 parameters 
float bulkhead_inner_edge1 =  3.2 //3.4
float bulkhead_outer_edge1 =  3.6 //3.8
float bulkhead_phi_offset1 =  0.0
int bulkhead_sides2 = 6    // Layer 2 parameters 
float bulkhead_inner_edge2 =  5.2 //5.6
float bulkhead_outer_edge2 =  5.6 //6.0
float bulkhead_phi_offset2 =  0.523598775
int bulkhead_sides3 = 12   // Layer 3 parameters 
float bulkhead_inner_edge3 =  3.65 //3.8
float bulkhead_outer_edge3 =  3.85 //4.0
float bulkhead_phi_offset3 =  0.13089693
int bulkhead_sides4 = 12   // Layer 4 parameters 
float bulkhead_inner_edge4 =  4.95 //5.2
float bulkhead_outer_edge4 =  5.15 //5.4
float bulkhead_phi_offset4 =  0.39269981

// added SFL0_May_2005
bool use_Layer0_geometry = true;

float L0wafer_thickness = 0.03     // cm
string L0wafer_material_name = "Layer0_Wafer"
string L0wafer_material_symbol = "L0Si"
float L0wafer_material_zeff = 14.0
float L0wafer_material_aeff = 28.085
float L0wafer_material_rho = 2.33  // gm/cm^3

// factors to multiply the default Layer0 wafer material
// for better description of multiple scattering
float MS_Factor_L0_ladder = 2.5
// end of SFL0_May_2005

// The geometrical information for the trough
float trough_length = 50.0 // cm
float trough_radius = 15.0
float trough_dphi = 3.1415927

// segmentation of trough material is controlled by these parameters:

int n_trough_z_bins = 14
int n_trough_phi_bins = 37

float trough_z_bin_edges = ( -92.6, -52.7, -47.7, -42.7, -36.2, -23.0, -10.6, 0., 10.6, 23.0, 36.2, 42.7, 47.7, 52.7, 92.6 )
float trough_phi_bin_edges = ( 0., 13.0, 19.0, 29.0, 35.0, 45.0, 51.0, 61.0, 67.0, 77.0, 83.0, 97.0, 103.0, 113.0, 119.0, 129.0, 135.0, 145.0, 151.0, 161.0, 167.0, 211.0, 217.0, 227.0, 233.0, 243.0, 249.0, 259.0, 265.0, 275.0, 281.0, 291.0, 297.0, 307.0, 313.0, 323.0, 329.0, 360.)


// fractional radiation length per phi bin
float radlength_zbin0 = ( 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0119, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0143 )
float radlength_zbin1 = ( 0.0062, 0.1400, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1400, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1400, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1400, 0.0119, 0.1735, 0.0062, 0.1735, 0.0062, 0.1400, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1400, 0.0062, 0.1735, 0.0062, 0.1735, 0.0143 )
float radlength_zbin2 = ( 0.0062, 0.1400, 0.0062, 0.1567, 0.0062, 0.1567, 0.0062, 0.1400, 0.0062, 0.1567, 0.0062, 0.1567, 0.0062, 0.1400, 0.0062, 0.1567, 0.0062, 0.1567, 0.0062, 0.1400, 0.0119, 0.1567, 0.0062, 0.1567, 0.0062, 0.1400, 0.0062, 0.1567, 0.0062, 0.1567, 0.0062, 0.1400, 0.0062, 0.1567, 0.0062, 0.1567, 0.0143 )
float radlength_zbin3 = ( 0.0062, 0.1066, 0.0062, 0.1567, 0.0062, 0.1567, 0.0062, 0.1066, 0.0062, 0.1567, 0.0062, 0.1567, 0.0062, 0.1066, 0.0062, 0.1567, 0.0062, 0.1567, 0.0062, 0.1066, 0.0119, 0.1567, 0.0062, 0.1567, 0.0062, 0.1066, 0.0062, 0.1567, 0.0062, 0.1567, 0.0062, 0.1066, 0.0062, 0.1567, 0.0062, 0.1567, 0.0143 )
float radlength_zbin4 = ( 0.0062, 0.0731, 0.0062, 0.1066, 0.0062, 0.1066, 0.0062, 0.0731, 0.0062, 0.1066, 0.0062, 0.1066, 0.0062, 0.0731, 0.0062, 0.1066, 0.0062, 0.1066, 0.0062, 0.0731, 0.0119, 0.1066, 0.0062, 0.1066, 0.0062, 0.0731, 0.0062, 0.1066, 0.0062, 0.1066, 0.0062, 0.0731, 0.0062, 0.1066, 0.0062, 0.1066, 0.0143 )
float radlength_zbin5 = ( 0.0062, 0.0396, 0.0062, 0.0564, 0.0062, 0.0564, 0.0062, 0.0396, 0.0062, 0.0564, 0.0062, 0.0564, 0.0062, 0.0396, 0.0062, 0.0564, 0.0062, 0.0564, 0.0062, 0.0396, 0.0119, 0.0564, 0.0062, 0.0564, 0.0062, 0.0396, 0.0062, 0.0564, 0.0062, 0.0564, 0.0062, 0.0396, 0.0062, 0.0564, 0.0062, 0.0564, 0.0143 )
float radlength_zbin6 = ( 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0119, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0143 )
float radlength_zbin7 = ( 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0119, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0062, 0.0143 )
float radlength_zbin8 = ( 0.0062, 0.0396, 0.0062, 0.0564, 0.0062, 0.0564, 0.0062, 0.0396, 0.0062, 0.0564, 0.0062, 0.0564, 0.0062, 0.0396, 0.0062, 0.0564, 0.0062, 0.0564, 0.0062, 0.0396, 0.0119, 0.0564, 0.0062, 0.0564, 0.0062, 0.0396, 0.0062, 0.0564, 0.0062, 0.0564, 0.0062, 0.0396, 0.0062, 0.0564, 0.0062, 0.0564, 0.0143 )
float radlength_zbin9 = ( 0.0062, 0.0731, 0.0062, 0.1066, 0.0062, 0.1066, 0.0062, 0.0731, 0.0062, 0.1066, 0.0062, 0.1066, 0.0062, 0.0731, 0.0062, 0.1066, 0.0062, 0.1066, 0.0062, 0.0731, 0.0119, 0.1066, 0.0062, 0.1066, 0.0062, 0.0731, 0.0062, 0.1066, 0.0062, 0.1066, 0.0062, 0.0731, 0.0062, 0.1066, 0.0062, 0.1066, 0.0143 )
float radlength_zbin10 = ( 0.0062, 0.1066, 0.0062, 0.1567, 0.0062, 0.1567, 0.0062, 0.1066, 0.0062, 0.1567, 0.0062, 0.1567, 0.0062, 0.1066, 0.0062, 0.1567, 0.0062, 0.1567, 0.0062, 0.1066, 0.0119, 0.1567, 0.0062, 0.1567, 0.0062, 0.1066, 0.0062, 0.1567, 0.0062, 0.1567, 0.0062, 0.1066, 0.0062, 0.1567, 0.0062, 0.1567, 0.0143 )
float radlength_zbin11 = ( 0.0062, 0.1400, 0.0062, 0.1567, 0.0062, 0.1567, 0.0062, 0.1400, 0.0062, 0.1567, 0.0062, 0.1567, 0.0062, 0.1400, 0.0062, 0.1567, 0.0062, 0.1567, 0.0062, 0.1400, 0.0119, 0.1567, 0.0062, 0.1567, 0.0062, 0.1400, 0.0062, 0.1567, 0.0062, 0.1567, 0.0062, 0.1400, 0.0062, 0.1567, 0.0062, 0.1567, 0.0143 )
float radlength_zbin12 = ( 0.0062, 0.1400, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1400, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1400, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1400, 0.0119, 0.1735, 0.0062, 0.1735, 0.0062, 0.1400, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1400, 0.0062, 0.1735, 0.0062, 0.1735, 0.0143 )
float radlength_zbin13 = ( 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0119, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0062, 0.1735, 0.0143 )

// averate density per phi bin, assuming R = 15 cm and thickness = 0.8 cm
float rho_zbin0 = ( 0.3601, 3.2879, 0.3601, 3.2877, 0.3601, 3.2877, 0.3601, 3.2879, 0.3601, 3.2877, 0.3601, 3.2877, 0.3601, 3.2879, 0.3601, 3.2877, 0.3601, 3.2877, 0.3601, 3.2879, 0.7162, 3.2877, 0.3601, 3.2877, 0.3601, 3.2879, 0.3601, 3.2877, 0.3601, 3.2877, 0.3601, 3.2879, 0.3601, 3.2877, 0.3601, 3.2877, 0.8656 )
float rho_zbin1 = ( 0.3601, 2.6736, 0.3601, 3.2877, 0.3601, 3.2877, 0.3601, 2.6736, 0.3601, 3.2877, 0.3601, 3.2877, 0.3601, 2.6736, 0.3601, 3.2877, 0.3601, 3.2877, 0.3601, 2.6736, 0.7162, 3.2877, 0.3601, 3.2877, 0.3601, 2.6736, 0.3601, 3.2877, 0.3601, 3.2877, 0.3601, 2.6736, 0.3601, 3.2877, 0.3601, 3.2877, 0.8656 )
float rho_zbin2 = ( 0.3601, 2.6736, 0.3601, 2.9788, 0.3601, 2.9788, 0.3601, 2.6736, 0.3601, 2.9788, 0.3601, 2.9788, 0.3601, 2.6736, 0.3601, 2.9788, 0.3601, 2.9788, 0.3601, 2.6736, 0.7162, 2.9788, 0.3601, 2.9788, 0.3601, 2.6736, 0.3601, 2.9788, 0.3601, 2.9788, 0.3601, 2.6736, 0.3601, 2.9788, 0.3601, 2.9788, 0.8656 )
float rho_zbin3 = ( 0.3601, 2.0736, 0.3601, 2.9788, 0.3601, 2.9788, 0.3601, 2.0736, 0.3601, 2.9788, 0.3601, 2.9788, 0.3601, 2.0736, 0.3601, 2.9788, 0.3601, 2.9788, 0.3601, 2.0736, 0.7162, 2.9788, 0.3601, 2.9788, 0.3601, 2.0736, 0.3601, 2.9788, 0.3601, 2.9788, 0.3601, 2.0736, 0.3601, 2.9788, 0.3601, 2.9788, 0.8656 )
float rho_zbin4 = ( 0.3601, 1.4880, 0.3601, 2.0736, 0.3601, 2.0736, 0.3601, 1.4880, 0.3601, 2.0736, 0.3601, 2.0736, 0.3601, 1.4880, 0.3601, 2.0736, 0.3601, 2.0736, 0.3601, 1.4880, 0.7162, 2.0736, 0.3601, 2.0736, 0.3601, 1.4880, 0.3601, 2.0736, 0.3601, 2.0736, 0.3601, 1.4880, 0.3601, 2.0736, 0.3601, 2.0736, 0.8656 )
float rho_zbin5 = ( 0.3601, 0.9169, 0.3601, 1.2006, 0.3601, 1.2006, 0.3601, 0.9169, 0.3601, 1.2006, 0.3601, 1.2006, 0.3601, 0.9169, 0.3601, 1.2006, 0.3601, 1.2006, 0.3601, 0.9169, 0.7162, 1.2006, 0.3601, 1.2006, 0.3601, 0.9169, 0.3601, 1.2006, 0.3601, 1.2006, 0.3601, 0.9169, 0.3601, 1.2006, 0.3601, 1.2006, 0.8656 )
float rho_zbin6 = ( 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.7162, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.8656 )
float rho_zbin7 = ( 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.7162, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.3601, 0.8656 )
float rho_zbin8 = ( 0.3601, 0.9169, 0.3601, 1.2006, 0.3601, 1.2006, 0.3601, 0.9169, 0.3601, 1.2006, 0.3601, 1.2006, 0.3601, 0.9169, 0.3601, 1.2006, 0.3601, 1.2006, 0.3601, 0.9169, 0.7162, 1.2006, 0.3601, 1.2006, 0.3601, 0.9169, 0.3601, 1.2006, 0.3601, 1.2006, 0.3601, 0.9169, 0.3601, 1.2006, 0.3601, 1.2006, 0.8656 )
float rho_zbin9 = ( 0.3601, 1.4880, 0.3601, 2.0736, 0.3601, 2.0736, 0.3601, 1.4880, 0.3601, 2.0736, 0.3601, 2.0736, 0.3601, 1.4880, 0.3601, 2.0736, 0.3601, 2.0736, 0.3601, 1.4880, 0.7162, 2.0736, 0.3601, 2.0736, 0.3601, 1.4880, 0.3601, 2.0736, 0.3601, 2.0736, 0.3601, 1.4880, 0.3601, 2.0736, 0.3601, 2.0736, 0.8656 )
float rho_zbin10 = ( 0.3601, 2.0736, 0.3601, 2.9788, 0.3601, 2.9788, 0.3601, 2.0736, 0.3601, 2.9788, 0.3601, 2.9788, 0.3601, 2.0736, 0.3601, 2.9788, 0.3601, 2.9788, 0.3601, 2.0736, 0.7162, 2.9788, 0.3601, 2.9788, 0.3601, 2.0736, 0.3601, 2.9788, 0.3601, 2.9788, 0.3601, 2.0736, 0.3601, 2.9788, 0.3601, 2.9788, 0.8656 )
float rho_zbin11 = ( 0.3601, 2.6736, 0.3601, 2.9788, 0.3601, 2.9788, 0.3601, 2.6736, 0.3601, 2.9788, 0.3601, 2.9788, 0.3601, 2.6736, 0.3601, 2.9788, 0.3601, 2.9788, 0.3601, 2.6736, 0.7162, 2.9788, 0.3601, 2.9788, 0.3601, 2.6736, 0.3601, 2.9788, 0.3601, 2.9788, 0.3601, 2.6736, 0.3601, 2.9788, 0.3601, 2.9788, 0.8656 )
float rho_zbin12 = ( 0.3601, 2.6736, 0.3601, 3.2877, 0.3601, 3.2877, 0.3601, 2.6736, 0.3601, 3.2877, 0.3601, 3.2877, 0.3601, 2.6736, 0.3601, 3.2877, 0.3601, 3.2877, 0.3601, 2.6736, 0.7162, 3.2877, 0.3601, 3.2877, 0.3601, 2.6736, 0.3601, 3.2877, 0.3601, 3.2877, 0.3601, 2.6736, 0.3601, 3.2877, 0.3601, 3.2877, 0.8656 )
float rho_zbin13 = ( 0.3601, 3.2879, 0.3601, 3.2877, 0.3601, 3.2877, 0.3601, 3.2879, 0.3601, 3.2877, 0.3601, 3.2877, 0.3601, 3.2879, 0.3601, 3.2877, 0.3601, 3.2877, 0.3601, 3.2879, 0.7162, 3.2877, 0.3601, 3.2877, 0.3601, 3.2879, 0.3601, 3.2877, 0.3601, 3.2877, 0.3601, 3.2879, 0.3601, 3.2877, 0.3601, 3.2877, 0.8656 )