--- Unused/pams/geometry/calbgeo/calbgeo1.g	2004-02-05 20:55:13.000000000 -0500
+++ New/pams/geometry/calbgeo/calbgeo1.g	2006-10-01 16:41:58.000000000 -0400
@@ -34,7 +34,7 @@
 external etsphit
 *
       Content CALB,CHLV,CPHI,CSUP,CPBP,CSCI,CSMD,CSMG,CSDA,CSMC,CSMB,CSME,
-              CSHI,CBTW
+              CSHI,CBTW,CPHT
 
       Structure CALG { version,  Rmin,     Etacut,   CrackWd,
                        FrontThk, CompThk,  AirThk,   BackThk,  SpaceThk, 
@@ -102,8 +102,8 @@
       Nmodule  = {60,60}     ! number of modules
       Shift    = {75,105}    ! starting azimuth of the first module   
       ModMap   = {1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
-                  0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
-                  0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
+                  1,1,1,1,1,1,1,1,1,1,1,1,0,1,1,
+                  1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
                   1,1,1,1,1,1,1,1,1,1,1,1,1,1,1}  ! Populated modules map
    Endfill
 *
@@ -185,10 +185,33 @@
           if(ii.eq.2.and.calg_ModMap(imod).gt.0) {
             Position CPHI AlphaZ=angular_offset Ncopy=imod;
           }
+          if(ii.eq.2.and.imod.eq.28.and.calg_ModMap(imod).eq.0) {
+!         position Crystal Detector in place of BEMC module 88
+            Create and Position CPHT AlphaZ=angular_offset Ncopy=imod;
+          }
 
       enddo
 EndBlock
 * -----------------------------------------------------------------------------
+Block CPHT
+      Material  CSI A=129.974  Z=54.023  Dens=4.51  RadL=1.85  AbsL=37.0 isvol=1
+      attribute CPHT  seen=1   colo=7
+      Shape  PCON Phi1=-3.0 DPhi=6.0 Nz=3,
+                      zi  = { 0,          cut_length,  Hleng},
+                      rmn = { Calg_rmin,  Calg_Rmin,   cut_radius},
+                      rmx = { Rmax,       Rmax,        Rmax };
+                      
+      Call CALBPAR(ag_imed,'ABSORBER')
+
+* define Birks law parameters
+      Call GSTPAR (ag_imed,'BIRK1',1.)
+      Call GSTPAR (ag_imed,'BIRK2',0.013)
+      Call GSTPAR (ag_imed,'BIRK3',9.6E-6)
+
+      HITS CPHT eta:0.05:(-1,0) phi:0.66:(-3.0,3.0) Eloss:0:(0,50) Birk:0:C(0,10)
+
+EndBlock
+* -----------------------------------------------------------------------------
 Block CPHI corresponds to a single module
       attribute CPHI  seen=1   colo=5
       Shape  PCON Phi1=-3.0 DPhi=6.0 Nz=3,