//
// $Id: MaterialGeometry_run2b.rcp,v 1.2 2006/02/14 17:47:20 shfu Exp $
//
// File: MaterialGeometry.rcp
// Purpose: Parameters defining material geometry (incld: beampipe, solnd)
//          with appended "old" code from (Hobbs et al.)
// Created: 25-MAR-1999, Abid Patwa, John Hobbs
//
// Modified: 01-APR-1999, Abid Patwa
//
// Modified: 20-May-2005, Shaohua Fu
// Added a rcp option to swich from RunIIa to RunIIb beampipe
// RunIIb beampipe is smaller than RunIIa one because of Layer0 detector
//
// $Revision: 1.2 $
//
// --------------------------------------------------------------------
//  Beampipe
// --------------------------------------------------------------------
// Beryllium (Be) beampipe  
// Parameters provided by B. Cooper -- private comm.;  Apr. 1, 1999
// see Eng. draw: \\d0server1\users\cooper\Dwgs_ps\Beamtube\358897b.ps
// bp_r of Be pipe = OD(1.5") - wall thickness(0.02"); bp_z is
// nominal pipe length (93.62" upto stainless steel transitions);
// bp_MatThick = wall thickness (0.02");
// Xo(Be) = 65.19 g/cm^2 and 35.38 cm.  (may contain 2% BeO)
// (all dimensions in cm, density=Rho in g-cm^-3)

string bp_MatName = "Beampipe"
string bp_MatSymbol = "bp"
float  bp_MatA = 9.012182
float  bp_MatZ = 4.0
float  bp_MatRho = 1.8480       // Now in g/cm^3
float  bp_MatX0 = 65.19         // Now in g/cm^2
float  bp_MatThick = 0.0508     // 0.02"*2.54
float  bp_Radius = 1.905        // cm
float  bp_Halflength = 118.8974 // half length in cm


// added SFL0_May_2005

bool   UseRunIIbBeampipe = true

string run2b_bp_MatName = "Beampipe"
string run2b_bp_MatSymbol = "bp"
float  run2b_bp_MatA = 9.012182
float  run2b_bp_MatZ = 4.0
float  run2b_bp_MatRho = 1.8480       // Now in g/cm^3
float  run2b_bp_MatX0 = 65.19         // Now in g/cm^2
float  run2b_bp_MatThick = 0.0508     // 0.02"*2.54
float  run2b_bp_Radius = 1.474        // cm
float  run2b_bp_Halflength = 87.7062  // half length in cm

// end of SFL0_May_2005


// --------------------------------------------------------------------
//  Solenoid
// --------------------------------------------------------------------
// reference document: Design of 2 T Super.Solenoid for Fermilab D0 Upgrade
//                     CEC July 12, 1993, Fermilab-Pub-098/047-E
// solenoid composed: Aluminum (Al) Shells/Shield, Cu+NbTi Conductor, and
//                    G-10/Kapton Coil Insulation;  
//                    (dimensions below in cm, density=Rho in g/cm^3)             
//                    There are 5.97 cm of Al, 0.28 cm of Cu+NbTi and 
//                    0.44 cm of G10 (at X0 = 22 g/cm^2 )
//
// All material is represented as a single, effective material located at
// the inner radius of the coil.  The A, Z and density is assumed to be Al,
// since it dominates.  However, the actual radiation length percentage is taken
// from the above reference document (0.88 X0).  Thus, an effective thickness
// is calculated as Thickness = %X0 * X0 / Rho = 0.88 * 24.6 / 2.7 = 8.02 cm.  (HLM)

string solnd_MatName = "Solenoid"
string solnd_MatSymbol = "solnd"
float  solnd_MatA = 26.981538             
float  solnd_MatZ = 13.0
float  solnd_MatRho = 2.7000       // Now in g/cm^3
float  solnd_MatX0 = 24.6          // Weighted average with thicknesses above
float  solnd_MatThick = 8.02       // Effective thickness (cm) -- see above.
float  solnd_Radius = 61.21        // Inner coil radius (not cryostat)
float  solnd_Halflength = 137.5    // Full length = 275.0 cm

// ----------------------------------------------------------------------

bool useCft = true
RCP CftMaterialGeometry = <cft_geometry CftMaterialGeometry>

bool useSmt = true
RCP SmtMaterialGeometry = <silicon_geometry SiMaterialGeometry_run2b>

bool useCal = false
RCP CalMaterialGeometry = <calorimeter_geometry CalMaterials>

bool useCps = true
RCP CpsMaterialGeometry = <cps_geometry CpsMaterialGeometry>

bool useFPS = true
RCP FpsMaterialGeometry = <fps_geometry fps_material>

bool useMuon = false  // Muon has no material geometry now...