generate_icarus_geometry_splitwire.pl File Reference

Go to the source code of this file.

Macros
#define	the   racetrack geometrical parameters

Functions
Math::BigFloat	precision (-16)
	GetOptions ("help|h"=>\$help,"suffix|s:s"=>\$suffix,"output|o:s"=>\$output,"concrete|c:s"=>\$thickness_over,"wires|w:s"=>\$wires,"vetocrt|v:s"=>\$crt)
	if (defined $help)
	if (!defined $suffix)
	gen_Define ()
	gen_Materials ()
	gen_RaceTracks ()
	gen_TPC ()
	gen_PMT ()
	gen_Mech_Structure ()
	gen_Cryostat ()
	if ($crt_on==1)
	gen_Enclosure ()
	gen_World ()
	write_fragments ()
sub	usage ()
sub	define_RaceTrack ()
	push (@gdmlFiles, $RACETRACK)
	open (RACETRACK) or die("Could not open file $RACETRACK for writing")
	for ($it=0;$it< $RaceTrack_number;$it++)
	open (PMT) or die("Could not open file $PMT for writing")
	open (MECH) or die("Could not open file $MECH for writing")
sub	gen_CRT ()

Variables
	else
	$inch = 2.54
	$YWirePitch = 0.3
	$UWirePitch = 0.3
	$VWirePitch = 0.3
	$UAngle = 60
	$VAngle = 60
	$SinUAngle = sin( deg2rad($UAngle) )
	$CosUAngle = cos( deg2rad($UAngle) )
	$TanUAngle = tan( deg2rad($UAngle) )
	$SinVAngle = sin( deg2rad($VAngle) )
	$CosVAngle = cos( deg2rad($VAngle) )
	$TanVAngle = tan( deg2rad($VAngle) )
	$UWire_ypitch = $UWirePitch / $CosUAngle
	$UWire_zpitch = $UWirePitch / $SinUAngle
	$VWire_ypitch = $VWirePitch / $CosVAngle
	$VWire_zpitch = $VWirePitch / $SinVAngle
	$TPCWireThickness = 0.015
	$TPCWireRadius = $TPCWireThickness /2
	$CommonWireLength = (1056 * $YWirePitch) / $SinUAngle
	$CPA_x = 0.15
	$WirePlaneSpacing = 0.3
	$MaxDrift = 148.2
	$SpaceWirePlToWall = 31.8
	$SpaceWirePlToWirePl = 85
	$SpaceTPCToFloor = 36
	$SpaceTPCToTopLAr = 29.5
	$UpstreamLArPadding = 90.0
	$DownstreamLArPadding = 90.0
	$CPA_epsi = 0.001
	$CPA_y = 320
	$CPA_z = 1800
	$CPAStrip_x = $CPA_x-$CPA_epsi
	$CPAStrip_y = 2.1
	$CPAStrip_pitch = 5
	$TPCActive_x = $MaxDrift
	$TPCActive_y = $CommonWireLength * $SinUAngle + 0.02
	$TPCActive_z = $CommonWireLength * $CosUAngle + (4640 - 1) * $UWire_zpitch - 0.000000000002
	$TPCActiveHalf_z = $TPCActive_z/2
	$TPCWirePlane_x = 2*$TPCWireThickness
	$TPCWirePlane_y = $TPCActive_y
	$TPCWirePlane_z = $TPCActive_z
	$TPCWirePlaneHalf_z = $TPCActiveHalf_z
	$TPC_x = $MaxDrift+ 6*$TPCWireThickness + 3*$WirePlaneSpacing + $CPA_x
	$DeltaLUCorner = $UWirePitch/($SinUAngle*$CosUAngle)
	$DeltaLVCorner = $VWirePitch/($SinVAngle*$CosVAngle)
	$LAr_x
	$LAr_y
	$LAr_z
	$LArHalf_z = $TPCActiveHalf_z + $UpstreamLArPadding
	$TPC_y = $LAr_y
	$TPC_z = $LAr_z
	$TPCHalf_z = $LArHalf_z
	$AlumThickness = 19
	$GaseousAr_y = 6.5
	$CryoDist = 20
	$Cryostat_x = $LAr_x + 2*$AlumThickness
	$Cryostat_y = $TPC_y + $GaseousAr_y + 2*$AlumThickness
	$Cryostat_z = $TPC_z + 2*$AlumThickness
	$LAr_x_orig = 2*$LAr_x + 2*$AlumThickness + $CryoDist
	$Cryostat_x_orig = $LAr_x_orig + 2*$AlumThickness
	$TPCinCryo_x [0] = - $TPC_x/2 - $CPA_x/2
	$posCat_x = 0
	$TPCinCryo_y = -$GaseousAr_y/2
	$TPCinCryo_z = 0.
	$TPCinCryo_zneg = (-$TPCHalf_z/2)
	$TPCinCryo_zpos = ( $TPCHalf_z/2)
	$posTPCActive0_x = 0.570000000000007
	$posTPCActive_y = 0
	$posTPCActive_z = $UpstreamLArPadding/2
	$struct_width = 10.
	$ext_struct_x = $struct_width
	$ext_struct_y = 341.8
	$ext_struct_z = 199.532
	$int_struct_x = $ext_struct_x
	$int_struct_y = $ext_struct_y - 2*$struct_width
	$int_struct_z = $ext_struct_z - 2*$struct_width
	$NumberPMT = 90
	$PMTthickness = 0.5
	$PMTradius = 4*$inch
	$PMTradiusOuter = 4*$inch
	$PMTradiusInner = $PMTradiusOuter - $PMTthickness
	$PMTGrid_offset = 2
	$PMTx = 0.557*$PMTradiusOuter
	$PMTtubeRmax = 4.225
	$PMTtubeRmin = 3.725
	$PMTtube_x = 2*$PMTtubeRmax
	$PMTPlane_x = 2*$PMTradius+$PMTGrid_offset+$PMTtube_x
	$PMTPlane_y = $ext_struct_y + 6.
	$PMTPlane_z = 9*$ext_struct_z + 2*$struct_width + 6.
	$PMTWiresOffset = 0.5
	$FoamPadding = 60
	$WarmVesselThickness = 27.4
	$WarmVesselThickness_x = 21.8
	$SpaceTop = 22.4
	$SpaceUpstream = 65.55
	$SpaceDownstream = 15.55
	$SpaceBottom = 18.0
	$ThermIns_x = 2*$Cryostat_x+2*$FoamPadding + 3*$CryoDist
	$ThermIns_y = $Cryostat_y+2*$FoamPadding+$SpaceTop+$SpaceBottom
	$ThermIns_z = $Cryostat_z+2*$FoamPadding+$SpaceUpstream+$SpaceDownstream
	$ThermInsInDetEncl_x = 0
	$ThermInsInDetEncl_z = 0
	$WarmVessel_x = $ThermIns_x+2*$WarmVesselThickness_x
	$WarmVessel_y = $ThermIns_y+2*$WarmVesselThickness
	$WarmVessel_z = $ThermIns_z+2*$WarmVesselThickness
	$WarmVesselInDetEncl_x = 0
	$WarmVesselInDetEncl_z = 0
	$Cryo1InWarmVessel_x = -$Cryostat_x/2 - $CryoDist
	$Cryo2InWarmVessel_x = $Cryostat_x/2 + $CryoDist
	$CryoInWarmVessel_z = 0
	$OriginXSet = 0
	$OriginYSet = 0
	$OriginZSet = 0
	$Building_y = 1040.0
	$ExpHall_y = 1170.0 + 80
	$Hall_x = 1890.0
	$Hall_y = $Building_y + $ExpHall_y
	$Hall_z = 3870.0
	$HallWallThicnekss = 60
	$Overburden_yDefault = 300
	$Overburden_x = 1830.0
	$Overburden_y = $concrete_on
	$Overburden_z = 3810.0
	$posOverburden_x = $OriginXSet
	$posOverburden_z = $OriginZSet
	$ExpHall_VertSpace = $ExpHall_y - $HallWallThicnekss - $Overburden_y/3
	$DetEnc_pad = 0.1
	$DetEnc_x = 1245
	$DetEnc_y = $ExpHall_VertSpace - 2*$DetEnc_pad
	$DetEnc_z = 3200
	$WarmVessel_FootHeight = 6.0
	$WarmVessel_FloorSpace = 10.16
	$WarmVessel_CenterToFloor = $WarmVessel_FootHeight + $WarmVessel_FloorSpace + $WarmVessel_y/2
	$World_x = 1e4
	$World_y = 1e4
	$World_z = 1e4
	$CRTSHELL_WV_OFFSET_Y = 15
	$CRTSHELL_WV_OFFSET_Z = 164.5865
	$Ground_y = 780.0
	$DetEncl_yOffset = $Ground_y - $Overburden_y/3 - $ExpHall_VertSpace/2 -$DetEnc_pad
	$ThermInsInDetEncl_y = -1*$DetEncl_yOffset
	$WarmVesselInDetEncl_y = $ThermInsInDetEncl_y
	$CryoInWarmVessel_y = $WarmVesselInDetEncl_y -$WarmVessel_y/2 + $FoamPadding + $WarmVesselThickness + $Cryostat_y/2
	exit
sub	gen_RaceTracks
	$RACETRACK = "icarus_racetracks" . $suffix . ".gdml"
print RACETRACK<< EOF;<?xml version='1.0'?>< gdml > EOF	$RTy_T = ($RaceTrack_ly - $RaceTrack_thickness )/2. + $RT_epsilon
	$RTy_B =-($RaceTrack_ly - $RaceTrack_thickness )/2. + $RT_epsilon
	$RTx_T = ($TPC_x - $RaceTrack_pitch -$RaceTrack_width )/2.
	$RTx_B = $RTx_T
	$RTx_U = $RTx_T
	$RTz_U = -($RaceTrackT_length - $RaceTrack_thickness)/2. - $RT_epsilon + 45 - ($RaceTrackT_length - $TPCActiveHalf_z)/2
	$RTx_D = $RTx_T
	$RTz_D = -$RTz_U
	$RTz_T = 45 - ($RaceTrackT_length - $TPCActiveHalf_z)/2
	$RTz_B = 45 - ($RaceTrackT_length - $TPCActiveHalf_z)/2
print RACETRACK<< EOF;< define > < positionname="posRaceTrackTInTPC"unit="cm"x="$RTx_T"y="$RTy_T"z="$RTz_T"/> < positionname="posRaceTrackBInTPC"unit="cm"x="$RTx_B"y="$RTy_B"z="$RTz_B"/> < positionname="posRaceTrackUInTPC"unit="cm"x="$RTx_U"y="0"z="$RTz_U"/> < positionname="posRaceTrackDInTPC"unit="cm"x="$RTx_D"y="0"z="$RTz_D"/> </define > <!--+++++++++++++++++++RACETRACKSolids+++++++++++++++++++--> EOF print RACETRACK<< EOF; < solids >< boxname="RaceTrackTVolume"x="$RaceTrack_width"y="$RaceTrack_thickness"z="$RaceTrackT_length"lunit="cm"/> < boxname="RaceTrackBVolume"x="$RaceTrack_width"y="$RaceTrack_thickness"z="$RaceTrackB_length"lunit="cm"/> < boxname="RaceTrackUVolume"x="$RaceTrack_width"y="$RaceTrack_thickness"z="$RaceTrackU_length"lunit="cm"/> < tubename="RaceTrackTTubeVolume"rmax="$RaceTrack_ExternalRadius"rmin="$RaceTrack_InnerRadius"z="$RaceTrackTTube_length"deltaphi="360"aunit="deg"lunit="cm"/> < tubename="RaceTrackBTubeVolume"rmax="$RaceTrack_ExternalRadius"rmin="$RaceTrack_InnerRadius"z="$RaceTrackBTube_length"deltaphi="360"aunit="deg"lunit="cm"/> < tubename="RaceTrackUTubeVolume"rmax="$RaceTrack_ExternalRadius"rmin="$RaceTrack_InnerRadius"z="$RaceTrackUTube_length"deltaphi="360"aunit="deg"lunit="cm"/> </solids > EOF print RACETRACK << EOF;< structure > <!--+++++++++++++++++RACETRACKLogicalVolumes++++++++++++++++++++--> < volumename="volRaceTrackTTubeVolume"> < materialrefref="STEEL_STAINLESS_Fe7Cr2Ni"/> < solidrefref="RaceTrackTTubeVolume"/> </volume >< volumename="volRaceTrackBTubeVolume"> < materialrefref="STEEL_STAINLESS_Fe7Cr2Ni"/> < solidrefref="RaceTrackBTubeVolume"/> </volume >< volumename="volRaceTrackUTubeVolume"> < materialrefref="STEEL_STAINLESS_Fe7Cr2Ni"/> < solidrefref="RaceTrackUTubeVolume"/> </volume > EOF print RACETRACK << EOF; <!--+++++++++++++++++RACETRACK Top Volume++++++++++++++++++++--> < volumename="volRaceTrackTVolume"> < materialrefref="LAr"/> < solidrefref="RaceTrackTVolume"/> EOF	$xTpos =-($RaceTrack_width - $RaceTrack_pitch - $RT_epsilon )/2.
print RACETRACK<< EOF;</volume > EOF print RACETRACK<< EOF; <!--+++++++++++++++++RACETRACK end structure+++++++++++++++++++--> </structure > EOF print RACETRACK<< EOF;</gdml > EOF sub	gen_PMT
	$PMT = ">" . $PMT
print PMT<< EOF;<?xml version='1.0'?> < gdml > EOF print PMT<< EOF; <!--+++++++++++++++++++PMT Solids++++++++++++++++++++++--> EOF print PMT<< EOF;< solids > < spherename="PMTVolume"rmin="$PMTradiusInner"rmax="$PMTradiusOuter"deltaphi="360"deltatheta="90"aunit="deg"lunit="cm"/> < spherename="PMTPassSphere"rmin="$PMTradiusInner"rmax="$PMTradiusOuter"deltaphi="360"deltatheta="90"aunit="deg"lunit="cm"/> < tubename="PMTPassTube"rmax="$PMTtubeRmax"rmin="$PMTtubeRmin"z="$PMTtube_x"deltaphi="360"aunit="deg"lunit="cm"/> < unionname="PMTPassVolume"> < firstref="PMTPassSphere"/> < secondref="PMTPassTube"/> < positionrefref="PMTtube"/> < rotationrefref="rMinus90AboutZ"/> </union ></solids > EOF print PMT<< EOF;< structure > <!--+++++++++++++++++PMTLogicalVolumes++++++++++++++++++++--> < volumename="volOpDetSensitive"> < materialrefref="LAr"/> < solidrefref="PMTVolume"/> </volume >< volumename="volNotOpDetSensitive"> < materialrefref="Glass"/> < solidrefref="PMTPassVolume"/> </volume ></structure > EOF print PMT<< EOF;</gdml > EOF sub	gen_Mech_Structure
	$MECH = ">" . $MECH

Macro Definition Documentation

#define the   racetrack geometrical parameters

Function Documentation

sub define_RaceTrack

(

)

Definition at line 1511 of file generate_icarus_geometry_splitwire.pl.

{
# from "Design, construction and tests of the ICARUS T600 detector",
# ICARUS Collaboration (S. Amerio et al.). Jul 2004. 82 pp. 
# Nucl.Instrum.Meth. A527 (2004) 329-410
#
#  4 parts: Top (T) Bottom (B) Upstream (U) Downstream (D)
  $RT_epsilon = .0001 ;# extra safety space;
  $RaceTrack_d = 3.4;  # cm race track tube diameter
  $RaceTrack_TubeThick = 0.08; #cm race track tube thick
#  $RaceTrack_lz = 1810.0; #cm length of race track structure
  $RaceTrack_lz = 905; #cm length of race track structure (half to accommodate for split wire TPC)
#  $RaceTrack_ly = 320.0 - 2*$RT_epsilon;  #cm length of race track tubes
  $RaceTrack_ly = 323.6201 ;  #cm height of race track structure (modified to be external to TPCActive_y
  $RaceTrack_number = 29; #number of race tracks
  $RaceTrack_pitch = 4.96; #distance between each race track tube
  $RaceTrack_ExternalRadius = $RaceTrack_d * 0.5;  # cm  external race track tube radius
  $RaceTrack_InnerRadius = $RaceTrack_d * 0.5-$RaceTrack_TubeThick;  # cm internal race track tube radius
  $RaceTrack_width = ($RaceTrack_pitch +$RT_epsilon) * $RaceTrack_number;
  $RaceTrackT_length = $RaceTrack_lz + $RT_epsilon; # guess about Horizontal length
  $RaceTrackTTube_length = $RaceTrackT_length - $RT_epsilon; # guess about Horizontal length
  $RaceTrackB_length = $RaceTrack_lz + $RT_epsilon; # guess about Horizontal length
  $RaceTrackBTube_length = $RaceTrackB_length - $RT_epsilon; # guess about Horizontal length
  $RaceTrackU_length  = $RaceTrack_ly - 2* $RaceTrack_d; # guess of vertical length
  $RaceTrackUTube_length  = $RaceTrackU_length - $RT_epsilon; # guess of vertical length
  $RaceTrack_thickness=$RaceTrack_d+2*$RT_epsilon; # 
}

for

(

)

Definition at line 1671 of file generate_icarus_geometry_splitwire.pl.

{
  print RACETRACK <<EOF;
      <physvol>
        <volumeref ref="volRaceTrackTTubeVolume"/>
        <position name="posRTkTTube$it" unit="cm" x="$xTpos" y="0" z="0" />
        <rotationref ref="rIdentity" />
      </physvol>
EOF
  $xTpos+= $RaceTrack_pitch ;
}

sub gen_CRT

(

)

Definition at line 2153 of file generate_icarus_geometry_splitwire.pl.

{
# Create the CRT fragment file name,
# add file to list of output GDML fragments,
# and open it

# using a python script to generate the CRT geometry
    $CRTPYTHON="gen_crt_frags.py";
    my $ret=`python $CRTPYTHON `;
    

    $CRT = "icarus_crt" . $suffix . ".gdml";
    push (@gdmlFiles, $CRT);
# expecting the python script to generate a file "icarus_crt_test.gdml"
#    $CRT = ">" . $CRT;
#    open(CRT) or die("Could not open file $CRT for writing");


# The standard XML prefix and starting the gdml
#    print CRT <<EOF;
#<?xml version='1.0'?>
#<gdml>
#EOF


close(CRT);
}

sub gen_Cryostat

(

)

Definition at line 1139 of file generate_ICARUS.pl.

{

# Create the cryostat fragment file name,
# add file to list of output GDML fragments,
# and open it
    $CRYO = "dune_10kT_Cryostat" . $suffix . ".gdml";
    push (@gdmlFiles, $CRYO);
    $CRYO = ">" . $CRYO;
    open(CRYO) or die("Could not open file $CRYO for writing");


# The standard XML prefix and starting the gdml
    print CRYO <<EOF;
<?xml version='1.0'?>
<gdml>
EOF


# All the cryostat solids.
print CRYO <<EOF;
<solids>

    <box name="Cryostat" lunit="cm" 
      x="$Cryostat_x" 
      y="$Cryostat_y" 
      z="$Cryostat_z"/>
    <box name="ArgonInterior" lunit="cm" 
      x="$LAr_x"
      y="$LAr_y"
      z="$LAr_z"/>
    <box name="GaseousArgon" lunit="cm" 
      x="$LAr_x"
      y="$GaseousAr_y"
      z="$LAr_z"/>
    <subtraction name="SteelShell">
      <first ref="Cryostat"/>
      <second ref="ArgonInterior"/>
    </subtraction>

    <box name="Cathode" lunit="cm"
      x="$CPA_x"
      y="$TPC_y"
      z="$TPC_z"/>

</solids>
EOF

# Cryostat structure
print CRYO <<EOF;
<structure>
    <volume name="volSteelShell">
      <materialref ref="STEEL_STAINLESS_Fe7Cr2Ni" />
      <solidref ref="SteelShell" />
    </volume>
    <volume name="volGaseousArgon">
      <materialref ref="ArGas"/>
      <solidref ref="GaseousArgon"/>
    </volume>

    <volume name="volCathode">
      <materialref ref="STEEL_STAINLESS_Fe7Cr2Ni" />
      <solidref ref="Cathode" />
    </volume>


    <volume name="volCryostat">
      <materialref ref="LAr" />
      <solidref ref="Cryostat" />
      <physvol>
        <volumeref ref="volGaseousArgon"/>
        <position name="posGaseousArgon" unit="cm" x="0" y="$LAr_y/2-$GaseousAr_y/2" z="0"/>
      </physvol>
      <physvol>
        <volumeref ref="volSteelShell"/>
        <position name="posSteelShell" unit="cm" x="0" y="0" z="0"/>
      </physvol>

      <physvol>
        <volumeref ref="volTPC"/>
        <positionref ref="posTPC0inCryo"/>
        <rotationref ref="rIdentity"/>
      </physvol>
      <physvol>
        <volumeref ref="volCathode" />
        <positionref ref="posRightCathode"/>
      </physvol>
      <physvol>
        <volumeref ref="volTPC"/>
        <positionref ref="posTPC1inCryo"/>
        <rotationref ref="rPlus180AboutY"/>
      </physvol>

      <physvol>
        <volumeref ref="volTPC"/>
        <positionref ref="posTPC2inCryo"/>
        <rotationref ref="rIdentity"/>
      </physvol>
      <physvol>
        <volumeref ref="volCathode" />
        <positionref ref="posLeftCathode"/>
      </physvol>
      <physvol>
        <volumeref ref="volTPC"/>
        <positionref ref="posTPC3inCryo"/>
        <rotationref ref="rPlus180AboutY"/>
      </physvol>

EOF



print CRYO <<EOF;
    </volume>
</structure>
</gdml>
EOF

close(CRYO);
}

sub gen_Define

(

)

Definition at line 454 of file generate_icarus_geometry_splitwire.pl.

{

# Create the <define> fragment file name, 
# add file to list of fragments,
# and open it
    $DEF = "icarus_Def" . $suffix . ".gdml";
    push (@gdmlFiles, $DEF);
    $DEF = ">" . $DEF;
    open(DEF) or die("Could not open file $DEF for writing");


print DEF <<EOF;
<?xml version='1.0'?>
<gdml>
<define>

<!--
-->
   <position name="PMTtube" unit="cm" x="0" y="0" z="@{[1.5*$PMTtube_x]}"/>
   <position name="posActiveInTPC0"   unit="cm" x="$posTPCActive0_x" y="$posTPCActive_y" z="$posTPCActive_z"/>
   <position name="posTPC00inCryo"    unit="cm" x="$TPCinCryo_x[0]" y="$TPCinCryo_y"    z="$TPCinCryo_zneg" />
   <position name="posTPC01inCryo"    unit="cm" x="$TPCinCryo_x[0]" y="$TPCinCryo_y"    z="$TPCinCryo_zpos" />
   <position name="posTPC10inCryo"    unit="cm" x="$TPCinCryo_x[1]" y="$TPCinCryo_y"    z="$TPCinCryo_zneg" />
   <position name="posTPC11inCryo"    unit="cm" x="$TPCinCryo_x[1]" y="$TPCinCryo_y"    z="$TPCinCryo_zpos" />
   <position name="posCathode"  unit="cm" x="$posCat_x"  y="$TPCinCryo_y"    z="$TPCinCryo_z" />
   <position name="posCryo1InWarmVessel"  unit="cm" x="$Cryo1InWarmVessel_x" y="$CryoInWarmVessel_y" z="$CryoInWarmVessel_z" />
   <position name="posCryo2InWarmVessel"  unit="cm" x="$Cryo2InWarmVessel_x" y="$CryoInWarmVessel_y" z="$CryoInWarmVessel_z" />
   <position name="posDetEncInWorld" unit="cm" x="$OriginXSet"     y="$DetEncl_yOffset"     z="$OriginZSet"/>
   <position name="posCenter"           unit="cm" x="0" y="0" z="0"/>
   <position name="posThermInsInDetEncl" unit="cm" x="$ThermInsInDetEncl_x" y="$ThermInsInDetEncl_y" z="$ThermInsInDetEncl_z"/>
   <position name="posWarmVesselInDetEncl" unit="cm" x="$WarmVesselInDetEncl_x" y="$WarmVesselInDetEncl_y" z="$WarmVesselInDetEncl_z"/>
   <position name="posCRTShellInDetEncl" unit="cm" x="0" y="$CRTSHELL_WV_OFFSET_Y" z="$CRTSHELL_WV_OFFSET_Z"/>
   <position name="posBuildingInWorld" unit="cm" x="0" y="@{[$Ground_y + $Building_y/2]}" z="0"/>
   <position name="posExpHallInWorld" unit="cm" x="0" y="@{[$Ground_y - $ExpHall_y/2 ]}" z="0"/>
   <rotation name="rPlus90AboutZPlus90AboutY"  unit="deg" x="0" y="90" z="90"/>
   <rotation name="rPlus90AboutX"       unit="deg" x="90" y="0" z="0"/>
   <rotation name="rMinus90AboutX"      unit="deg" x="-90" y="0" z="0"/>
   <rotation name="rPlus90AboutY"       unit="deg" x="0" y="90"   z="0"/>
   <rotation name="rPlus90AboutZ"       unit="deg" x="0" y="0"   z="90"/>
   <rotation name="rMinus90AboutZ"      unit="deg" x="0" y="0"   z="-90"/>
   <rotation name="rMinus90AboutY"      unit="deg" x="0" y="270" z="0"/>
   <rotation name="rMinus90AboutYMinus90AboutX"       unit="deg" x="270" y="270" z="0"/>
   <rotation name="rPlus90VAngleAboutX" unit="deg" x="@{[90-$VAngle]}" y="0"   z="0"/>  
   <rotation name="rPlus90UAngleAboutX" unit="deg" x="@{[90+$UAngle]}" y="0"   z="0"/>
   <rotation name="rPlus180AboutX"      unit="deg" x="180" y="0"   z="0"/>
   <rotation name="rPlus180AboutY"      unit="deg" x="0" y="180"   z="0"/>
   <rotation name="rPlus180AboutZ"      unit="deg" x="0" y="0"   z="180"/>
   <rotation name="rPlus180AboutXZ"     unit="deg" x="180" y="0"   z="180"/>
   <rotation name="rIdentity"           unit="deg" x="0" y="0"   z="0"/>
   <rotation name="rPlusUAngleAboutX"   unit="deg" x="$UAngle" y="0" z="0"/>
   <rotation name="rMinusVAngleAboutX"  unit="deg" x="300" y="0" z="0"/>
</define>
</gdml>
EOF
    close (DEF);
}

sub gen_Enclosure

(

)

Definition at line 2192 of file generate_icarus_geometry_splitwire.pl.

{

# Create the detector enclosure fragment file name,
# add file to list of output GDML fragments,
# and open it
    $ENCL = "icarus_DetEnclosure" . $suffix . ".gdml";
    push (@gdmlFiles, $ENCL);
    $ENCL = ">" . $ENCL;
    open(ENCL) or die("Could not open file $ENCL for writing");


# The standard XML prefix and starting the gdml
    print ENCL <<EOF;
<?xml version='1.0'?>
<gdml>
EOF

# All the detector enclosure solids.
print ENCL <<EOF;
<solids>

    <box name="DetEnclosure" lunit="cm" 
      x="$DetEnc_x"
      y="$DetEnc_y"
      z="$DetEnc_z"/>
     
    <box name="ThermIns" lunit="cm" 
      x="$ThermIns_x"
      y="$ThermIns_y"
      z="$ThermIns_z"/>

    <box name="WarmVessel" lunit="cm" 
      x="$WarmVessel_x"
      y="$WarmVessel_y"
      z="$WarmVessel_z"/>

    <box name="ThermInsInterior" lunit="cm" 
      x="@{[$ThermIns_x - 2*$FoamPadding]}"
      y="@{[$ThermIns_y - 2*$FoamPadding]}"
      z="@{[$ThermIns_z - 2*$FoamPadding]}"/>

    <box name="WarmVesselInterior" lunit="cm" 
      x="@{[$WarmVessel_x - 2*$WarmVesselThickness_x]}"
      y="@{[$WarmVessel_y - 2*$WarmVesselThickness]}"
      z="@{[$WarmVessel_z - 2*$WarmVesselThickness]}"/>

    <subtraction name="ThermInsShell">
      <first ref="ThermIns"/>
      <second ref="ThermInsInterior"/>
    </subtraction>

    <subtraction name="WarmVesselShell">
      <first ref="WarmVessel"/>
      <second ref="WarmVesselInterior"/>
    </subtraction>

</solids>
EOF


# Detector enclosure structure    
    print ENCL <<EOF;
<structure>

    <volume name="volThermIns">
      <materialref ref="Polyurethane"/> 
      <solidref ref="ThermInsShell"/>
    </volume>

    <volume name="volWarmVessel">
      <materialref ref="STEEL_STAINLESS_Fe7Cr2Ni_WV"/>
      <solidref ref="WarmVesselShell"/>
    </volume>

    <volume name="volDetEnclosure">
      <materialref ref="Air"/>
      <solidref ref="DetEnclosure"/>

    <physvol>
      <volumeref ref="volThermIns"/>
      <positionref ref="posThermInsInDetEncl"/>
    </physvol>

    <physvol>
      <volumeref ref="volCRT_Shell"/>
      <positionref ref="posCRTShellInDetEncl"/>
    </physvol>

    <physvol>
      <volumeref ref="volWarmVessel"/>
      <positionref ref="posWarmVesselInDetEncl"/>
    </physvol>

    <physvol>
      <volumeref ref="volCryostat"/>
      <positionref ref="posCryo1InWarmVessel"/>
    </physvol>

    <physvol>
      <volumeref ref="volCryostat"/>
      <positionref ref="posCryo2InWarmVessel"/>
    </physvol>
        
EOF
if ($crt_on==1) 
{
#+++++++++++++++++++++++++ Position CRT modules ++++++++++++++++++++++++++
    #$YPosMINOSSide = $Cryostat_y/2+$MINOSModWidth/2;


} #end of crt_on 

 else { 

print ENCL <<EOF;

        <!-- This GDML version has no crt and uses much less memory -->

EOF

}

print ENCL <<EOF;

</volume>

</structure>
</gdml>
EOF

close(ENCL);
}

sub gen_Materials

(

)

Definition at line 517 of file generate_icarus_geometry_splitwire.pl.

{

# Create the <materials> fragment file name,
# add file to list of output GDML fragments,
# and open it
    $MAT = "icarus_Materials" . $suffix . ".gdml";
    push (@gdmlFiles, $MAT);
    $MAT = ">" . $MAT;
    open(MAT) or die("Could not open file $MAT for writing");


  print MAT <<EOF;
<materials>
  <element name="videRef" formula="VACUUM" Z="1">  <atom value="1"/> </element>
  <element name="bromine" formula="Br" Z="35"> <atom value="79.904"/> </element>
  <element name="hydrogen" formula="H" Z="1">  <atom value="1.0079"/> </element>
  <element name="nitrogen" formula="N" Z="7">  <atom value="14.0067"/> </element>
  <element name="oxygen" formula="O" Z="8">  <atom value="15.999"/> </element>
  <element name="aluminum" formula="Al" Z="13"> <atom value="26.9815"/>  </element>
  <element name="silicon" formula="Si" Z="14"> <atom value="28.0855"/>  </element>
  <element name="carbon" formula="C" Z="6">  <atom value="12.0107"/>  </element>
  <element name="potassium" formula="K" Z="19"> <atom value="39.0983"/>  </element>
  <element name="chromium" formula="Cr" Z="24"> <atom value="51.9961"/>  </element>
  <element name="iron" formula="Fe" Z="26"> <atom value="55.8450"/>  </element>
  <element name="nickel" formula="Ni" Z="28"> <atom value="58.6934"/>  </element>
  <element name="calcium" formula="Ca" Z="20"> <atom value="40.078"/>   </element>
  <element name="magnesium" formula="Mg" Z="12"> <atom value="24.305"/>   </element>
  <element name="sodium" formula="Na" Z="11"> <atom value="22.99"/>    </element>
  <element name="titanium" formula="Ti" Z="22"> <atom value="47.867"/>   </element>
  <element name="argon" formula="Ar" Z="18"> <atom value="39.9480"/>  </element>
  <element name="sulphur" formula="S" Z="16"> <atom value="32.065"/>  </element>
  <element name="phosphorus" formula="P" Z="15"> <atom value="30.973"/>  </element>
  <element name="copper" formula="Cu" Z="29"> <atom value="63.5463"/>  </element>
  <element name="manganese" formula="Mn" Z="25"> <atom value="54.938043"/>  </element>
  <element name="vanadium" formula="V" Z="23"> <atom value="50.94151"/>  </element>

  <material name="Vacuum" formula="Vacuum">
   <D value="1.e-25" unit="g/cm3"/>
   <fraction n="1.0" ref="videRef"/>
  </material>

  <material name="ALUMINUM_Al" formula="ALUMINUM_Al">
   <D value="2.6990" unit="g/cm3"/>
   <fraction n="1.0000" ref="aluminum"/>
  </material>

  <material name="ALUMINUM_CRYO" formula="ALUMINUM_CRYO">
   <D value="0.5991" unit="g/cm3"/>
   <fraction n="1.0000" ref="aluminum"/>
  </material>

  <material name="SILICON_Si" formula="SILICON_Si">
   <D value="2.3300" unit="g/cm3"/>
   <fraction n="1.0000" ref="silicon"/>
  </material>

  <material name="epoxy_resin" formula="C38H40O6Br4">
   <D value="1.1250" unit="g/cm3"/>
   <composite n="38" ref="carbon"/>
   <composite n="40" ref="hydrogen"/>
   <composite n="6" ref="oxygen"/>
   <composite n="4" ref="bromine"/>
  </material>

  <material name="SiO2" formula="SiO2">
   <D value="2.2" unit="g/cm3"/>
   <composite n="1" ref="silicon"/>
   <composite n="2" ref="oxygen"/>
  </material>

  <material name="Al2O3" formula="Al2O3">
   <D value="3.97" unit="g/cm3"/>
   <composite n="2" ref="aluminum"/>
   <composite n="3" ref="oxygen"/>
  </material>

  <material name="Fe2O3" formula="Fe2O3">
   <D value="5.24" unit="g/cm3"/>
   <composite n="2" ref="iron"/>
   <composite n="3" ref="oxygen"/>
  </material>

  <material name="CaO" formula="CaO">
   <D value="3.35" unit="g/cm3"/>
   <composite n="1" ref="calcium"/>
   <composite n="1" ref="oxygen"/>
  </material>

  <material name="MgO" formula="MgO">
   <D value="3.58" unit="g/cm3"/>
   <composite n="1" ref="magnesium"/>
   <composite n="1" ref="oxygen"/>
  </material>

  <material name="Na2O" formula="Na2O">
   <D value="2.27" unit="g/cm3"/>
   <composite n="2" ref="sodium"/>
   <composite n="1" ref="oxygen"/>
  </material>

  <material name="TiO2" formula="TiO2">
   <D value="4.23" unit="g/cm3"/>
   <composite n="1" ref="titanium"/>
   <composite n="2" ref="oxygen"/>
  </material>

  <material name="FeO" formula="FeO">
   <D value="5.745" unit="g/cm3"/>
   <composite n="1" ref="iron"/>
   <composite n="1" ref="oxygen"/>
  </material>

  <material name="CO2" formula="CO2">
   <D value="1.562" unit="g/cm3"/>
   <composite n="1" ref="iron"/>
   <composite n="2" ref="oxygen"/>
  </material>

  <material name="P2O5" formula="P2O5">
   <D value="1.562" unit="g/cm3"/>
   <composite n="2" ref="phosphorus"/>
   <composite n="5" ref="oxygen"/>
  </material>

  <material formula=" " name="DUSEL_Rock">
    <D value="2.82" unit="g/cm3"/>
    <fraction n="0.5267" ref="SiO2"/>
    <fraction n="0.1174" ref="FeO"/>
    <fraction n="0.1025" ref="Al2O3"/>
    <fraction n="0.0473" ref="MgO"/>
    <fraction n="0.0422" ref="CO2"/>
    <fraction n="0.0382" ref="CaO"/>
    <fraction n="0.0240" ref="carbon"/>
    <fraction n="0.0186" ref="sulphur"/>
    <fraction n="0.0053" ref="Na2O"/>
    <fraction n="0.00070" ref="P2O5"/>
    <fraction n="0.0771" ref="oxygen"/>
  </material> 

  <material name="fibrous_glass">
   <D value="2.74351" unit="g/cm3"/>
   <fraction n="0.600" ref="SiO2"/>
   <fraction n="0.118" ref="Al2O3"/>
   <fraction n="0.001" ref="Fe2O3"/>
   <fraction n="0.224" ref="CaO"/>
   <fraction n="0.034" ref="MgO"/>
   <fraction n="0.010" ref="Na2O"/>
   <fraction n="0.013" ref="TiO2"/>
  </material>

  <material name="FR4">
   <D value="1.98281" unit="g/cm3"/>
   <fraction n="0.47" ref="epoxy_resin"/>
   <fraction n="0.53" ref="fibrous_glass"/>
  </material>

  <material name="STEEL_STAINLESS_Fe7Cr2Ni" formula="STEEL_STAINLESS_Fe7Cr2Ni">
   <D value="7.9300" unit="g/cm3"/>
   <fraction n="0.0010" ref="carbon"/>
   <fraction n="0.1792" ref="chromium"/>
   <fraction n="0.7298" ref="iron"/>
   <fraction n="0.0900" ref="nickel"/>
  </material>

  <material name="STEEL_STAINLESS_Fe7Cr2Ni_WV" formula="STEEL_STAINLESS_Fe7Cr2Ni_WV">
   <!-- <D value="0.3844" unit="g/cm3"/> -->
   <D value="0.38897" unit="g/cm3"/>
   <fraction n="0.0010" ref="carbon"/>
   <fraction n="0.1792" ref="chromium"/>
   <fraction n="0.7298" ref="iron"/>
   <fraction n="0.0900" ref="nickel"/>
  </material>

  <material name="LAr" formula="LAr">
   <D value="1.40" unit="g/cm3"/>
   <fraction n="1.0000" ref="argon"/>
  </material>

  <material name="ArGas" formula="ArGas">
   <D value="0.00166" unit="g/cm3"/>
   <fraction n="1.0" ref="argon"/>
  </material>

  <material formula=" " name="Air">
   <D value="0.001205" unit="g/cm3"/>
   <fraction n="0.781154" ref="nitrogen"/>
   <fraction n="0.209476" ref="oxygen"/>
   <fraction n="0.00934" ref="argon"/>
  </material>

  <material formula=" " name="G10">
   <D value="1.7" unit="g/cm3"/>
   <fraction n="0.2805" ref="silicon"/>
   <fraction n="0.3954" ref="oxygen"/>
   <fraction n="0.2990" ref="carbon"/>
   <fraction n="0.0251" ref="hydrogen"/>
  </material>

  <material formula=" " name="Granite">
   <D value="2.7" unit="g/cm3"/>
   <fraction n="0.438" ref="oxygen"/>
   <fraction n="0.257" ref="silicon"/>
   <fraction n="0.222" ref="sodium"/>
   <fraction n="0.049" ref="aluminum"/>
   <fraction n="0.019" ref="iron"/>
   <fraction n="0.015" ref="potassium"/>
  </material>

  <material formula=" " name="ShotRock">
   <D value="@{[2.7*0.6]}" unit="g/cm3"/>
   <fraction n="0.438" ref="oxygen"/>
   <fraction n="0.257" ref="silicon"/>
   <fraction n="0.222" ref="sodium"/>
   <fraction n="0.049" ref="aluminum"/>
   <fraction n="0.019" ref="iron"/>
   <fraction n="0.015" ref="potassium"/>
  </material>

  <material formula=" " name="Dirt">
   <D value="1.7" unit="g/cm3"/>
   <fraction n="0.438" ref="oxygen"/>
   <fraction n="0.257" ref="silicon"/>
   <fraction n="0.222" ref="sodium"/>
   <fraction n="0.049" ref="aluminum"/>
   <fraction n="0.019" ref="iron"/>
   <fraction n="0.015" ref="potassium"/>
  </material>

  <material formula=" " name="Concrete">
   <D value="2.3" unit="g/cm3"/>
   <fraction n="0.530" ref="oxygen"/>
   <fraction n="0.335" ref="silicon"/>
   <fraction n="0.060" ref="calcium"/>
   <fraction n="0.015" ref="sodium"/>
   <fraction n="0.020" ref="iron"/>
   <fraction n="0.040" ref="aluminum"/>
  </material>

  <material formula="H2O" name="Water">
   <D value="1.0" unit="g/cm3"/>
   <fraction n="0.1119" ref="hydrogen"/>
   <fraction n="0.8881" ref="oxygen"/>
  </material>

  <material formula="Ti" name="Titanium">
   <D value="4.506" unit="g/cm3"/>
   <fraction n="1." ref="titanium"/>
  </material>

  <material name="TPB" formula="TPB">
   <D value="1.40" unit="g/cm3"/>
   <fraction n="1.0000" ref="argon"/>
  </material>

  <material name="Glass">
   <D value="2.74351" unit="g/cm3"/>
   <fraction n="0.600" ref="SiO2"/>
   <fraction n="0.118" ref="Al2O3"/>
   <fraction n="0.001" ref="Fe2O3"/>
   <fraction n="0.224" ref="CaO"/>
   <fraction n="0.034" ref="MgO"/>
   <fraction n="0.010" ref="Na2O"/>
   <fraction n="0.013" ref="TiO2"/>
  </material>

  <material name="Acrylic">
   <D value="1.19" unit="g/cm3"/>
   <fraction n="0.600" ref="carbon"/>
   <fraction n="0.320" ref="oxygen"/>
   <fraction n="0.080" ref="hydrogen"/>
  </material>

  <material name="Polystyrene">
   <D unit="g/cm3" value="1.06"/>
   <fraction n="0.077418" ref="hydrogen"/>
   <fraction n="0.922582" ref="carbon"/>
  </material>

  <material name="Polyurethane" formula="C27H36N2O10" >
   <D value="0.9" unit="g/cm3"/>
   <composite n="27" ref="carbon"/>
   <composite n="36" ref="hydrogen"/>
   <composite n="2" ref="nitrogen"/>
   <composite n="10" ref="oxygen"/>
  </material>

  <material name="STEEL_A992">
   <D unit="g/cm3" value="7.85"/>
   <fraction n="0.0022"  ref="carbon"/>
   <fraction n="0.005"   ref="copper"/>
   <fraction n="0.01"    ref="manganese"/>
   <fraction n="0.0044"  ref="nickel"/>
   <fraction n="0.00034" ref="phosphorus"/>
   <fraction n="0.0039"  ref="silicon"/>
   <fraction n="0.00044" ref="sulphur"/>
   <fraction n="0.001"   ref="vanadium"/>
   <fraction n="0.97272" ref="iron" />
  </material>


</materials>
EOF

close(MAT);
}

gen_Mech_Structure

(

)

gen_PMT

(

)

gen_RaceTracks

(

)

sub gen_TPC

(

)

Definition at line 830 of file generate_icarus_geometry_splitwire.pl.

{
# Create the TPC fragment file name,
# add file to list of output GDML fragments,
# and open it
    $TPC = "icarus_TPC" . $suffix . ".gdml";
    push (@gdmlFiles, $TPC);
    $TPC = ">" . $TPC;
    open(TPC) or die("Could not open file $TPC for writing");


# The standard XML prefix and starting the gdml
    print TPC <<EOF;
<?xml version='1.0'?>
<gdml>
EOF


# All the TPC solids save the wires.
print TPC <<EOF;
<solids>
    <box name="TPC" lunit="cm" 
      x="$TPC_x" 
      y="$TPC_y" 
      z="$TPCHalf_z"/>
    <box name="TPCPlane" lunit="cm" 
      x="$TPCWirePlane_x" 
      y="$TPCWirePlane_y" 
      z="$TPCWirePlaneHalf_z"/>
    <box name="TPCActive" lunit="cm"
      x="$TPCActive_x"
      y="$TPCActive_y"
      z="$TPCActiveHalf_z"/>
EOF


#++++++++++++++++++++++++++++ Wire Solids ++++++++++++++++++++++++++++++

# Set number of wires to default to zero, when $wires_on = 0, for a low memory 
# version. But if $wires_on = 1, calculate the number of wires on each side of each
# plane to be used in the for loops

my $NumberHorizontalWires = 0;
my $NumberCornerUWires = 0;
my $NumberCommonUWires = 0;
my $NumberCornerVWires = 0;
my $NumberCommonVWires = 0;


if ($wires_on == 1)
{
   # Number of wires in one corner
   #$NumberCornerVWires = int( $TPCWirePlane_y/$VWire_ypitch );
   #$NumberCornerWWires = int( $TPCWirePlane_y/$WWire_ypitch );
    $NumberCornerUWires = 480;
    $NumberCornerVWires = 480;

   # Number of wires in one corner + the number to cover the whole corner.
    $NumberCornerExtUWires = 528;
    $NumberCornerExtVWires = 528;

   # Total number of wires touching one vertical (longer) side
   # Note that the total number of wires per plane is this + another set of corner wires
   # $NumberSideUWires = int( $TPCWirePlane_z/$UWire_zpitch );
   # $NumberSideVWires = int( $TPCWirePlane_z/$VWire_zpitch );
   # $NumberSideWWires = int( $TPCWirePlane_z/$WWire_zpitch );

   # Number of wires per side that aren't cut off by the corner
   # $NumberCommonUWires = $NumberSideUWires - $NumberCornerUWires;
   # $NumberCommonVWires = $NumberSideVWires - $NumberCornerVWires;
    #$NumberCommonWWires = $NumberSideWWires - $NumberCornerWWires;
    #$NumberCommonUWires  = 4640;
    #$NumberCommonVWires  = 4640;
   $NumberCommonUWires  = 2056;
   $NumberCommonVWires  = 2056;

   # number of wires on the vertical plane
   #$NumberHorizontalWires = int( ($TPCWirePlane_y-$TPCWireThickness)/$UWirePitch );
   #Number of wires inthe Y plane-->Horizontal plane Induction I
    $NumberHorizontalWires = 1056;

}

# These XML comments throughout make the GDML file easier to navigate
print TPC <<EOF;

<!--+++++++++++++++++++ Y Wire Solids ++++++++++++++++++++++-->

EOF

if ($wires_on==1) 
{

#CommonWire = wires with same length

   print TPC <<EOF;
    <tube name="TPCWireYCommon"
      rmax="$TPCWireRadius"
      z="$TPCWirePlaneHalf_z"
      deltaphi="360"
      aunit="deg"
      lunit="cm"/>
EOF

} else { 

print TPC <<EOF;

        <!-- This GDML version has no wires and uses much less memory -->

EOF

}

print TPC <<EOF;


<!--+++++++++++++++++++ U Wire Solids ++++++++++++++++++++++-->


EOF

# The corner wires for the U plane
if ($wires_on==1) 
{
#FAKE CORNER 
    $length = $CommonWireLength + $DeltaLUCorner/2;

   for ($i = 0; $i < $NumberCornerExtUWires; ++$i)
    {
         $length -= $DeltaLUCorner;

print TPC <<EOF;
    <tube name="TPCWireU$i"
      rmax="$TPCWireRadius"
      z="$length"
      deltaphi="360"
      aunit="deg"
      lunit="cm"/>
EOF

#print(" $i $length \n");

   } #ends FAKE CORNER

#REAL CORNER 
    $length = $CommonWireLength;

   for ($i = 0; $i < $NumberCornerUWires; ++$i)
    {
         $length -= $DeltaLUCorner;

print TPC <<EOF;
    <tube name="TPCWireCornerU$i"
      rmax="$TPCWireRadius"
      z="$length"
      deltaphi="360"
      aunit="deg"
      lunit="cm"/>
EOF

#print(" $i $length \n");

    } #ends REAL CORNER
    
print TPC <<EOF;
    <tube name="TPCWireUCommon"
      rmax="$TPCWireRadius"
      z="$CommonWireLength"
      deltaphi="360"
      aunit="deg"
      lunit="cm"/>
EOF

} else { 

print TPC <<EOF;

                   <!-- no wires in this GDML -->

EOF

}


print TPC <<EOF;


<!--+++++++++++++++++++ V Wire Solids ++++++++++++++++++++++-->


EOF

# The corner wires for the V plane
if ($wires_on==1) 
{
# FAKE CORNER
   $length = $CommonWireLength + $DeltaLVCorner/2;

    for ($i = 0; $i < $NumberCornerExtVWires; ++$i)
   {

        $length -= $DeltaLVCorner;

   print TPC <<EOF;
    <tube name="TPCWireV$i"
      rmax="$TPCWireRadius"
      z="$length"
      deltaphi="360"
      aunit="deg"
      lunit="cm"/>
EOF

#print(" $i $length \n");

    } #ends FAKE CORNER

# REAL CORNER
   $length = $CommonWireLength;

    for ($i = 0; $i < $NumberCornerVWires; ++$i)
   {

        $length -= $DeltaLVCorner;

   print TPC <<EOF;
    <tube name="TPCWireCornerV$i"
      rmax="$TPCWireRadius"
      z="$length"
      deltaphi="360"
      aunit="deg"
      lunit="cm"/>
EOF

#print(" $i $length \n");

    } #ends REAL CORNER


   print TPC <<EOF;
    <tube name="TPCWireVCommon"
      rmax="$TPCWireRadius"
      z="$CommonWireLength"
      deltaphi="360"
      aunit="deg"
      lunit="cm"/>
EOF

} else { 

print TPC <<EOF;

                   <!-- no wires in this GDML -->

EOF

}

# Begin structure and create the vertical wire logical volume
print TPC <<EOF;
</solids>
<structure>
    <volume name="volTPCActive">
      <materialref ref="LAr"/>
      <solidref ref="TPCActive"/>
    </volume>

<!--+++++++++++++++++ Wire Logical Volumes ++++++++++++++++++++-->

EOF


if ($wires_on==1) 
{ 

  # Common Y wire logical volume, referenced many times
  print TPC <<EOF;
    <volume name="volTPCWireYCommon">
      <materialref ref="STEEL_STAINLESS_Fe7Cr2Ni" />
      <solidref ref="TPCWireYCommon" />
    </volume>
EOF

  # Fake Corner U wires logical volumes
  for ($i = 0; $i < $NumberCornerExtUWires; ++$i)
  {
  print TPC <<EOF;
    <volume name="volTPCWireU$i">
      <materialref ref="STEEL_STAINLESS_Fe7Cr2Ni" />
      <solidref ref="TPCWireU$i" />
    </volume>
EOF

  }
  
  # Real Corner U wires logical volumes
  for ($i = 0; $i < $NumberCornerUWires; ++$i)
  {
  print TPC <<EOF;
    <volume name="volTPCWireCornerU$i">
      <materialref ref="STEEL_STAINLESS_Fe7Cr2Ni" />
      <solidref ref="TPCWireCornerU$i" />
    </volume>
EOF

  }
  
  # Common U wire logical volume, referenced many times
  print TPC <<EOF;
    <volume name="volTPCWireUCommon">
      <materialref ref="STEEL_STAINLESS_Fe7Cr2Ni" />
      <solidref ref="TPCWireUCommon" />
    </volume>
EOF

  # Fake Corner V wires logical volumes
  for ($i = 0; $i < $NumberCornerExtVWires; ++$i)
  {
  print TPC <<EOF;
    <volume name="volTPCWireV$i">
      <materialref ref="STEEL_STAINLESS_Fe7Cr2Ni" />
      <solidref ref="TPCWireV$i" />
    </volume>
EOF

  }
  
  # Real Corner V wires logical volumes
  for ($i = 0; $i < $NumberCornerVWires; ++$i)
  {
  print TPC <<EOF;
    <volume name="volTPCWireCornerV$i">
      <materialref ref="STEEL_STAINLESS_Fe7Cr2Ni" />
      <solidref ref="TPCWireCornerV$i" />
    </volume>
EOF

  }
  # Common V wire logical volume, referenced many times
  print TPC <<EOF;
    <volume name="volTPCWireVCommon">
      <materialref ref="STEEL_STAINLESS_Fe7Cr2Ni" />
      <solidref ref="TPCWireVCommon" />
    </volume>
EOF

} else { 


print TPC <<EOF;

        <!-- This GDML version has no wires and uses much less memory -->

EOF

}


#+++++++++++++++++++++++++ Position physical wires ++++++++++++++++++++++++++

#            ++++++++++++++++++++++  Y Plane  +++++++++++++++++++++++

# Create Y plane logical volume
print TPC <<EOF;


<!--+++++++++++++++++++++ Y Plane ++++++++++++++++++++++++-->


    <volume name="volTPCPlaneY">
      <materialref ref="LAr"/>
      <solidref ref="TPCPlane"/>
EOF

if ($wires_on==0)
{
print TPC <<EOF;

           <!-- no wires -->

EOF

} else {

    $ypos = ($NumberHorizontalWires+1)*$YWirePitch/2;

    for ($i = 0; $i < $NumberHorizontalWires; ++$i)
    {
        $ypos -= $YWirePitch;


print TPC <<EOF;
      <physvol>
        <volumeref ref="volTPCWireYCommon"/>
        <position name="posTPCWireY$i" unit="cm" x="0" y="$ypos " z="0"/>
        <rotationref ref="rIdentity"/>
      </physvol>
EOF

        #print("0 $ypos \n");

        #$ypos -= $YWirePitch;

    }

} #ends else


#            ++++++++++++++++++++++  U Plane  +++++++++++++++++++++++

# End U plane and create U plane logical volume
print TPC <<EOF;
    </volume>


<!--+++++++++++++++++++++ U Plane ++++++++++++++++++++++++-->


    <volume name="volTPCPlaneU">
      <materialref ref="LAr"/>
      <solidref ref="TPCPlane"/>
EOF

if ($wires_on==0)
{
print TPC <<EOF;

           <!-- no wires -->

EOF

} else {

#CORNERS

    $ypos1 = $TPCActive_y/2;
    $zpos1 = -$TPCActiveHalf_z/2 + $CommonWireLength * $CosUAngle + $UWire_zpitch - $UWire_zpitch/2 - $UWire_zpitch/2;

    $ypos2 = -$TPCActive_y/2;
    $zpos2 = -$TPCActiveHalf_z/2;

   #Positioning of 480 real U corner wires.
   for ($i = 0; $i < $NumberCornerUWires; ++$i)
    {
         $ypos1 += $UWire_ypitch  ;
         $zpos2 -= $UWire_zpitch  ;

         $ypos = ($ypos1+$ypos2)/2;
         $zpos = ($zpos1+$zpos2)/2;

        #print("U Corner wires: $i $zpos $ypos , ");
        
print TPC <<EOF;
      <physvol>
        <volumeref ref="volTPCWireCornerU$i"/>
        <position name="posTPCWireU$i" unit="cm" x="0" y="$ypos " z="$zpos"/>
        <rotationref ref="rPlusUAngleAboutX"/>
      </physvol>
EOF
    }

    $ypos1 = $TPCActive_y/2 - $UWire_ypitch + $UWire_ypitch/2;
    $zpos1 = -$TPCActiveHalf_z/2 + $CommonWireLength * $CosUAngle + $UWire_zpitch - $UWire_zpitch/2;

    $ypos2 = - $TPCActive_y/2;
    $zpos2 = -$TPCActiveHalf_z/2;
    
    #Positioning of 528 fake U corners wires.
    for ($i = 0; $i < $NumberCornerExtUWires; ++$i)
    {
      $ypos1 += $UWire_ypitch;
      $zpos2 -= $UWire_zpitch;

      $ypos = ($ypos1+$ypos2)/2;
      $zpos = ($zpos1+$zpos2)/2;

      $ypos = - $ypos;
      $zpos = - $zpos;

print TPC <<EOF;
      <physvol>
        <volumeref ref="volTPCWireU$i"/>
        <position name="posTPCWireU@{[$i+$NumberCommonUWires+$NumberCornerUWires]}" unit="cm" x="0" y="$ypos " z="$zpos"/>
        <rotationref ref="rPlusUAngleAboutX"/>
      </physvol>
EOF
    } #ends CORNER


#Common Wires
    
   $zpos = (-$TPCActiveHalf_z +  $CommonWireLength * $CosUAngle + $UWire_zpitch)/2. - $UWire_zpitch/2;
   #print("common wires $zpos \n");

    for ($i = 0; $i < $NumberCommonUWires; ++$i)
    {

print TPC <<EOF;
      <physvol>
        <volumeref ref="volTPCWireUCommon"/>
        <position name="posTPCWireU@{[$i+$NumberCornerUWires]}" unit="cm" x="0" y="0 " z="$zpos"/>
        <rotationref ref="rPlusUAngleAboutX"/>
      </physvol>
EOF

    #print("U wires $i $zpos 0 \n");
        $zpos += $UWire_zpitch;

    }


} #ends else

#            ++++++++++++++++++++++  V Plane  +++++++++++++++++++++++

# End V plane and create V plane logical volume
print TPC <<EOF;
    </volume>

<!--+++++++++++++++++++++ V Plane ++++++++++++++++++++++++-->


    <volume name="volTPCPlaneV">
      <materialref ref="LAr"/>
      <solidref ref="TPCPlane"/>
EOF

if ($wires_on==0)
{
print TPC <<EOF;

           <!-- no wires -->

EOF

} else {

#CORNERS

    $ypos1 = $TPCActive_y/2;
    $zpos1 = -$TPCActiveHalf_z/2;

    $ypos2 = -$TPCActive_y/2;
    $zpos2 = -$TPCActiveHalf_z/2 + $CommonWireLength * $CosVAngle + $VWire_zpitch - $VWire_zpitch/2 - $VWire_zpitch/2 ;

   #Positioning of 480 real V corners wires.

   for ($i = 0; $i < $NumberCornerVWires; ++$i)
    {
         $ypos1 -= $VWire_ypitch  ;
         $zpos2 -= $VWire_zpitch  ;

         $ypos = ($ypos1+$ypos2)/2;
         $zpos = ($zpos1+$zpos2)/2;

        #print("V Corner wires: $i $zpos $ypos , ");


print TPC <<EOF;
      <physvol>
        <volumeref ref="volTPCWireCornerV$i"/>
        <position name="posTPCWireV$i" unit="cm" x="0" y="$ypos " z="$zpos"/>
        <rotationref ref="rMinusVAngleAboutX"/>
      </physvol>
EOF
    }
    
    $ypos1 = $TPCActive_y/2 + $VWire_ypitch - $VWire_ypitch/2;
    $zpos1 = -$TPCActiveHalf_z/2 + $VWire_zpitch - $VWire_zpitch/2;
    
    $ypos2 = -$TPCActive_y/2;
    $zpos2 = -$TPCActiveHalf_z/2 + $CommonWireLength * $CosVAngle;

    #Positioning of 528 fake V corners wires.
    for ($i = 0; $i < $NumberCornerExtVWires; ++$i)
    {
       $ypos1 -= $VWire_ypitch;
       $zpos2 -= $VWire_zpitch;

       $ypos = ($ypos1+$ypos2)/2;
       $zpos = ($zpos1+$zpos2)/2;

       $ypos = - $ypos;
       $zpos = - $zpos;

print TPC <<EOF;
      <physvol>
        <volumeref ref="volTPCWireV$i"/>
        <position name="posTPCWireV@{[$i+$NumberCommonVWires+$NumberCornerVWires]}" unit="cm" x="0" y="$ypos " z="$zpos"/>
        <rotationref ref="rMinusVAngleAboutX"/>
      </physvol>
EOF

        #print(" $zpos $ypos \n");

    } #ends CORNERS

   #Common Wires

   $zpos = (-$TPCActiveHalf_z +  $CommonWireLength * $CosVAngle + $VWire_zpitch) / 2. - $VWire_zpitch/2;

    for ($i = 0; $i < $NumberCommonVWires; ++$i)
    {

print TPC <<EOF;
      <physvol>
        <volumeref ref="volTPCWireVCommon"/>
        <position name="posTPCWireV@{[$i+$NumberCornerVWires]}" unit="cm" x="0" y="0 " z="$zpos"/>
        <rotationref ref="rMinusVAngleAboutX"/>
      </physvol>
EOF

        #print("V wires: $i $zpos  0 \n");
        $zpos += $VWire_zpitch;

    }


} #ends else

print TPC <<EOF;
    </volume>
EOF

#+++++++++++++++++++++ Position physical wires Above +++++++++++++++++++++

my $VolY_x = (-$TPC_x/2) + 3*$WirePlaneSpacing; #+ $TPCWirePlane_x/2;    
my $VolU_x = (-$TPC_x/2) + 2*$WirePlaneSpacing; #+ $TPCWirePlane_x/2;    
my $VolV_x = (-$TPC_x/2) + 1*$WirePlaneSpacing; #+ $TPCWirePlane_x/2;  


#wrap up the TPC file
print TPC <<EOF;

    <volume name="volTPC0">
      <materialref ref="LAr" />
      <solidref ref="TPC" />
      <physvol>
       <volumeref ref="volTPCPlaneY" />
       <position name="posTPCPlaneY" unit="cm" x="$VolY_x" y="0" z="@{[$UpstreamLArPadding/2]}" />
     </physvol>
     <physvol>
       <volumeref ref="volTPCPlaneU" />
       <position name="posTPCPlaneU" unit="cm" x="$VolU_x" y="0" z="@{[$UpstreamLArPadding/2]}" />
     </physvol>
     <physvol>
       <volumeref ref="volTPCPlaneV" />
       <position name="posTPCPlaneV" unit="cm" x="$VolV_x" y="0" z="@{[$UpstreamLArPadding/2]}" />
     </physvol>
     <physvol>
     <volumeref ref="volTPCActive"/>
     <positionref ref="posActiveInTPC0"/>
     </physvol> 
     <physvol> 
       <volumeref ref="volRaceTrackTVolume"/>
       <positionref ref="posRaceTrackTInTPC"/>
       <rotationref ref="rIdentity"/>
     </physvol>
     <physvol> 
       <volumeref ref="volRaceTrackBVolume"/>
       <positionref ref="posRaceTrackBInTPC"/>
       <rotationref ref="rIdentity"/>
     </physvol>
     <physvol> 
       <volumeref ref="volRaceTrackUVolume"/>
       <positionref ref="posRaceTrackUInTPC"/>
       <rotationref ref="rPlus90AboutX"/>
     </physvol>
    </volume>

</structure>
</gdml>
EOF

    close(GDML);

} #end of sub gen_TPC

sub gen_World

(

)

Definition at line 2332 of file generate_icarus_geometry_splitwire.pl.

GetOptions	(	"help|h"	,
		\	$help,
		"suffix|s:s"	,
		\	$suffix,
		"output|o:s"	,
		\	$output,
		"concrete|c:s"	,
		\	$thickness_over,
		"wires|w:s"	,
		\	$wires,
		"vetocrt|v:s"	,
		\	$crt
	)

if

(

defined

$help

)

Definition at line 25 of file generate_icarus_geometry_splitwire.pl.

{
    # If the user requested help, print the usage notes and exit.
    usage();
    exit;
}

if

(

!defined

$suffix

)

Definition at line 32 of file generate_icarus_geometry_splitwire.pl.

{
    # The user didn't supply a suffix, so append nothing to the file
    # names.
    $suffix = "";
}

if

(

$crt_on = =1

)

Definition at line 422 of file generate_icarus_geometry_splitwire.pl.

{gen_CRT();}     # places CRT: CERN modules top, eves; MINOS modules sides, DC bottom

open

(

RACETRACK

)

open

(

PMT

)

open

(

MECH

)

Math::BigFloat precision

(

-

16

)

push	(	@	gdmlFiles,
			$RACETRACK
	)

sub usage

(

)

Definition at line 436 of file generate_icarus_geometry_splitwire.pl.

{
    print "Usage: $0 [-h|--help] [-o|--output <fragments-file>] [-s|--suffix <string>]";
    print " [-c|--concrete <double>] [-w|--wires <wire or no wire geometry>] [-v <crt or no crt geometry>] \n";
    print "       if -o is omitted, output goes to STDOUT; <fragments-file.xml> is input to make_gdml.pl\n";
    print "       -s <string> appends the string to the file names; useful for multiple detector versions\n";
    print "       -c <double> set the thickness in cm of the concrete overburden\n";
    print "          (default is 300 and 0 means no overburden)\n";
    print "       -w <1> geometry with wires, <0> geometry with no wires\n";
    print "       -v <1> geometry with CRT, <0> geometry with no CRT \n";
    print "       -h prints this message, then quits\n";
    print "Remind: the file without wires has to be <filename_nowires.gdml> \n";
}

sub write_fragments

(

)

Definition at line 1168 of file generate_gdml_simple.pl.

{
    # The output file is a list of the GDML sub-files created by this
    # script.

    if ( ! defined $output )
    {
        $output = "-"; # write to STDOUT 
    }

    # Set up the output file.
    $OUTPUT = ">" . $output;
    open(OUTPUT) or die("Could not open file $OUTPUT");

    print OUTPUT <<EOF;
<?xml version='1.0'?>

<!-- Input to Geometry/gdml/make_gdml.pl; define the GDML fragments
     that will be zipped together to create a detector description. 
     -->

<config>

   <constantfiles>

      <!-- These files contain GDML <constant></constant>
           blocks. They are read in separately, so they can be
           interpreted into the remaining GDML. See make_gdml.pl for
           more information. 
           -->
           
EOF

    foreach $filename (@defFiles)
    {
        print OUTPUT <<EOF;
      <filename> $filename </filename>
EOF
    }

    print OUTPUT <<EOF;

   </constantfiles>

   <gdmlfiles>

      <!-- The GDML file fragments to be zipped together. -->

EOF

    foreach $filename (@gdmlFiles)
    {
        print OUTPUT <<EOF;
      <filename> $filename </filename>
EOF
    }

    print OUTPUT <<EOF;

   </gdmlfiles>

</config>
EOF

    close(OUTPUT);
}

Variable Documentation

$AlumThickness = 19

Definition at line 196 of file generate_icarus_geometry_splitwire.pl.

$Building_y = 1040.0

Definition at line 348 of file generate_icarus_geometry_splitwire.pl.

$CommonWireLength = (1056 * $YWirePitch) / $SinUAngle

Definition at line 112 of file generate_icarus_geometry_splitwire.pl.

$CosUAngle = cos( deg2rad($UAngle) )

Definition at line 91 of file generate_icarus_geometry_splitwire.pl.

$CosVAngle = cos( deg2rad($VAngle) )

Definition at line 95 of file generate_icarus_geometry_splitwire.pl.

$CPA_epsi = 0.001

Definition at line 138 of file generate_icarus_geometry_splitwire.pl.

$CPA_x = 0.15

Definition at line 119 of file generate_icarus_geometry_splitwire.pl.

$CPA_y = 320

Definition at line 139 of file generate_icarus_geometry_splitwire.pl.

$CPA_z = 1800

Definition at line 140 of file generate_icarus_geometry_splitwire.pl.

$CPAStrip_pitch = 5

Definition at line 143 of file generate_icarus_geometry_splitwire.pl.

$CPAStrip_x = $CPA_x-$CPA_epsi

Definition at line 141 of file generate_icarus_geometry_splitwire.pl.

$CPAStrip_y = 2.1

Definition at line 142 of file generate_icarus_geometry_splitwire.pl.

$CRTSHELL_WV_OFFSET_Y = 15

Definition at line 391 of file generate_icarus_geometry_splitwire.pl.

$CRTSHELL_WV_OFFSET_Z = 164.5865

Definition at line 392 of file generate_icarus_geometry_splitwire.pl.

$Cryo1InWarmVessel_x = -$Cryostat_x/2 - $CryoDist

Definition at line 323 of file generate_icarus_geometry_splitwire.pl.

$Cryo2InWarmVessel_x = $Cryostat_x/2 + $CryoDist

Definition at line 324 of file generate_icarus_geometry_splitwire.pl.

$CryoDist = 20

Definition at line 198 of file generate_icarus_geometry_splitwire.pl.

$CryoInWarmVessel_y = $WarmVesselInDetEncl_y -$WarmVessel_y/2 + $FoamPadding + $WarmVesselThickness + $Cryostat_y/2

Definition at line 401 of file generate_icarus_geometry_splitwire.pl.

$CryoInWarmVessel_z = 0

Definition at line 327 of file generate_icarus_geometry_splitwire.pl.

$Cryostat_x = $LAr_x + 2*$AlumThickness

Definition at line 200 of file generate_icarus_geometry_splitwire.pl.

$Cryostat_x_orig = $LAr_x_orig + 2*$AlumThickness

Definition at line 211 of file generate_icarus_geometry_splitwire.pl.

$Cryostat_y = $TPC_y + $GaseousAr_y + 2*$AlumThickness

Definition at line 204 of file generate_icarus_geometry_splitwire.pl.

$Cryostat_z = $TPC_z + 2*$AlumThickness

Definition at line 205 of file generate_icarus_geometry_splitwire.pl.

$DeltaLUCorner = $UWirePitch/($SinUAngle*$CosUAngle)

Definition at line 174 of file generate_icarus_geometry_splitwire.pl.

$DeltaLVCorner = $VWirePitch/($SinVAngle*$CosVAngle)

Definition at line 175 of file generate_icarus_geometry_splitwire.pl.

$DetEnc_pad = 0.1

Definition at line 376 of file generate_icarus_geometry_splitwire.pl.

$DetEnc_x = 1245

Definition at line 377 of file generate_icarus_geometry_splitwire.pl.

$DetEnc_y = $ExpHall_VertSpace - 2*$DetEnc_pad

Definition at line 378 of file generate_icarus_geometry_splitwire.pl.

$DetEnc_z = 3200

Definition at line 379 of file generate_icarus_geometry_splitwire.pl.

$DetEncl_yOffset = $Ground_y - $Overburden_y/3 - $ExpHall_VertSpace/2 -$DetEnc_pad

Definition at line 398 of file generate_icarus_geometry_splitwire.pl.

$DownstreamLArPadding = 90.0

Definition at line 133 of file generate_icarus_geometry_splitwire.pl.

$ExpHall_VertSpace = $ExpHall_y - $HallWallThicnekss - $Overburden_y/3

Definition at line 375 of file generate_icarus_geometry_splitwire.pl.

$ExpHall_y = 1170.0 + 80

Definition at line 349 of file generate_icarus_geometry_splitwire.pl.

$ext_struct_x = $struct_width

Definition at line 245 of file generate_icarus_geometry_splitwire.pl.

$ext_struct_y = 341.8

Definition at line 246 of file generate_icarus_geometry_splitwire.pl.

$ext_struct_z = 199.532

Definition at line 247 of file generate_icarus_geometry_splitwire.pl.

$FoamPadding = 60

Definition at line 285 of file generate_icarus_geometry_splitwire.pl.

$GaseousAr_y = 6.5

Definition at line 197 of file generate_icarus_geometry_splitwire.pl.

$Ground_y = 780.0

Definition at line 397 of file generate_icarus_geometry_splitwire.pl.

$Hall_x = 1890.0

Definition at line 351 of file generate_icarus_geometry_splitwire.pl.

$Hall_y = $Building_y + $ExpHall_y

Definition at line 352 of file generate_icarus_geometry_splitwire.pl.

$Hall_z = 3870.0

Definition at line 353 of file generate_icarus_geometry_splitwire.pl.

$HallWallThicnekss = 60

Definition at line 355 of file generate_icarus_geometry_splitwire.pl.

$inch = 2.54

Definition at line 76 of file generate_icarus_geometry_splitwire.pl.

$int_struct_x = $ext_struct_x

Definition at line 248 of file generate_icarus_geometry_splitwire.pl.

$int_struct_y = $ext_struct_y - 2*$struct_width

Definition at line 249 of file generate_icarus_geometry_splitwire.pl.

$int_struct_z = $ext_struct_z - 2*$struct_width

Definition at line 250 of file generate_icarus_geometry_splitwire.pl.

$LAr_x

Initial value:

=     $CPA_x 
              + 2*($TPC_x + $SpaceWirePlToWall)

Definition at line 180 of file generate_icarus_geometry_splitwire.pl.

$LAr_x_orig = 2*$LAr_x + 2*$AlumThickness + $CryoDist

Definition at line 209 of file generate_icarus_geometry_splitwire.pl.

$LAr_y

Initial value:

=     $TPCActive_y 
              + $SpaceTPCToFloor 
              + $SpaceTPCToTopLAr

Definition at line 182 of file generate_icarus_geometry_splitwire.pl.

$LAr_z

Initial value:

=     $TPCActive_z
              + $UpstreamLArPadding 
              + $DownstreamLArPadding

Definition at line 185 of file generate_icarus_geometry_splitwire.pl.

$LArHalf_z = $TPCActiveHalf_z + $UpstreamLArPadding

Definition at line 188 of file generate_icarus_geometry_splitwire.pl.

$MaxDrift = 148.2

Definition at line 121 of file generate_icarus_geometry_splitwire.pl.

$MECH = ">" . $MECH

Definition at line 1854 of file generate_icarus_geometry_splitwire.pl.

$NumberPMT = 90

Definition at line 255 of file generate_icarus_geometry_splitwire.pl.

$OriginXSet = 0

Definition at line 342 of file generate_icarus_geometry_splitwire.pl.

$OriginYSet = 0

Definition at line 343 of file generate_icarus_geometry_splitwire.pl.

$OriginZSet = 0

Definition at line 344 of file generate_icarus_geometry_splitwire.pl.

$Overburden_x = 1830.0

Definition at line 364 of file generate_icarus_geometry_splitwire.pl.

$Overburden_y = $concrete_on

Definition at line 365 of file generate_icarus_geometry_splitwire.pl.

$Overburden_yDefault = 300

Definition at line 356 of file generate_icarus_geometry_splitwire.pl.

$Overburden_z = 3810.0

Definition at line 366 of file generate_icarus_geometry_splitwire.pl.

$PMT = ">" . $PMT

Definition at line 1769 of file generate_icarus_geometry_splitwire.pl.

$PMTGrid_offset = 2

Definition at line 260 of file generate_icarus_geometry_splitwire.pl.

$PMTPlane_x = 2*$PMTradius+$PMTGrid_offset+$PMTtube_x

Definition at line 267 of file generate_icarus_geometry_splitwire.pl.

$PMTPlane_y = $ext_struct_y + 6.

Definition at line 269 of file generate_icarus_geometry_splitwire.pl.

$PMTPlane_z = 9*$ext_struct_z + 2*$struct_width + 6.

Definition at line 271 of file generate_icarus_geometry_splitwire.pl.

$PMTradius = 4*$inch

Definition at line 257 of file generate_icarus_geometry_splitwire.pl.

$PMTradiusInner = $PMTradiusOuter - $PMTthickness

Definition at line 259 of file generate_icarus_geometry_splitwire.pl.

$PMTradiusOuter = 4*$inch

Definition at line 258 of file generate_icarus_geometry_splitwire.pl.

$PMTthickness = 0.5

Definition at line 256 of file generate_icarus_geometry_splitwire.pl.

$PMTtube_x = 2*$PMTtubeRmax

Definition at line 265 of file generate_icarus_geometry_splitwire.pl.

$PMTtubeRmax = 4.225

Definition at line 263 of file generate_icarus_geometry_splitwire.pl.

$PMTtubeRmin = 3.725

Definition at line 264 of file generate_icarus_geometry_splitwire.pl.

$PMTWiresOffset = 0.5

Definition at line 272 of file generate_icarus_geometry_splitwire.pl.

$PMTx = 0.557*$PMTradiusOuter

Definition at line 261 of file generate_icarus_geometry_splitwire.pl.

$posCat_x = 0

Definition at line 219 of file generate_icarus_geometry_splitwire.pl.

$posOverburden_x = $OriginXSet

Definition at line 370 of file generate_icarus_geometry_splitwire.pl.

$posOverburden_z = $OriginZSet

Definition at line 372 of file generate_icarus_geometry_splitwire.pl.

$posTPCActive0_x = 0.570000000000007

Definition at line 235 of file generate_icarus_geometry_splitwire.pl.

$posTPCActive_y = 0

Definition at line 237 of file generate_icarus_geometry_splitwire.pl.

$posTPCActive_z = $UpstreamLArPadding/2

Definition at line 238 of file generate_icarus_geometry_splitwire.pl.

$RACETRACK = "icarus_racetracks" . $suffix . ".gdml"

Definition at line 1549 of file generate_icarus_geometry_splitwire.pl.

$RTx_B = $RTx_T

Definition at line 1566 of file generate_icarus_geometry_splitwire.pl.

$RTx_D = $RTx_T

Definition at line 1569 of file generate_icarus_geometry_splitwire.pl.

$RTx_T = ($TPC_x - $RaceTrack_pitch -$RaceTrack_width )/2.

Definition at line 1565 of file generate_icarus_geometry_splitwire.pl.

$RTx_U = $RTx_T

Definition at line 1567 of file generate_icarus_geometry_splitwire.pl.

$RTy_B =-($RaceTrack_ly - $RaceTrack_thickness )/2. + $RT_epsilon

Definition at line 1564 of file generate_icarus_geometry_splitwire.pl.

print RACETRACK<<EOF; <?xml version='1.0'?><gdml> EOF $RTy_T = ($RaceTrack_ly - $RaceTrack_thickness )/2. + $RT_epsilon

Definition at line 1563 of file generate_icarus_geometry_splitwire.pl.

$RTz_B = 45 - ($RaceTrackT_length - $TPCActiveHalf_z)/2

Definition at line 1572 of file generate_icarus_geometry_splitwire.pl.

$RTz_D = -$RTz_U

Definition at line 1570 of file generate_icarus_geometry_splitwire.pl.

$RTz_T = 45 - ($RaceTrackT_length - $TPCActiveHalf_z)/2

Definition at line 1571 of file generate_icarus_geometry_splitwire.pl.

$RTz_U = -($RaceTrackT_length - $RaceTrack_thickness)/2. - $RT_epsilon + 45 - ($RaceTrackT_length - $TPCActiveHalf_z)/2

Definition at line 1568 of file generate_icarus_geometry_splitwire.pl.

$SinUAngle = sin( deg2rad($UAngle) )

Definition at line 90 of file generate_icarus_geometry_splitwire.pl.

$SinVAngle = sin( deg2rad($VAngle) )

Definition at line 94 of file generate_icarus_geometry_splitwire.pl.

$SpaceBottom = 18.0

Definition at line 291 of file generate_icarus_geometry_splitwire.pl.

$SpaceDownstream = 15.55

Definition at line 290 of file generate_icarus_geometry_splitwire.pl.

$SpaceTop = 22.4

Definition at line 288 of file generate_icarus_geometry_splitwire.pl.

$SpaceTPCToFloor = 36

Definition at line 128 of file generate_icarus_geometry_splitwire.pl.

$SpaceTPCToTopLAr = 29.5

Definition at line 129 of file generate_icarus_geometry_splitwire.pl.

$SpaceUpstream = 65.55

Definition at line 289 of file generate_icarus_geometry_splitwire.pl.

$SpaceWirePlToWall = 31.8

Definition at line 124 of file generate_icarus_geometry_splitwire.pl.

$SpaceWirePlToWirePl = 85

Definition at line 125 of file generate_icarus_geometry_splitwire.pl.

$struct_width = 10.

Definition at line 244 of file generate_icarus_geometry_splitwire.pl.

$TanUAngle = tan( deg2rad($UAngle) )

Definition at line 92 of file generate_icarus_geometry_splitwire.pl.

$TanVAngle = tan( deg2rad($VAngle) )

Definition at line 96 of file generate_icarus_geometry_splitwire.pl.

$ThermIns_x = 2*$Cryostat_x+2*$FoamPadding + 3*$CryoDist

Definition at line 294 of file generate_icarus_geometry_splitwire.pl.

$ThermIns_y = $Cryostat_y+2*$FoamPadding+$SpaceTop+$SpaceBottom

Definition at line 295 of file generate_icarus_geometry_splitwire.pl.

$ThermIns_z = $Cryostat_z+2*$FoamPadding+$SpaceUpstream+$SpaceDownstream

Definition at line 296 of file generate_icarus_geometry_splitwire.pl.

$ThermInsInDetEncl_x = 0

Definition at line 301 of file generate_icarus_geometry_splitwire.pl.

$ThermInsInDetEncl_y = -1*$DetEncl_yOffset

Definition at line 399 of file generate_icarus_geometry_splitwire.pl.

$ThermInsInDetEncl_z = 0

Definition at line 303 of file generate_icarus_geometry_splitwire.pl.

$TPC_x = $MaxDrift+ 6*$TPCWireThickness + 3*$WirePlaneSpacing + $CPA_x

Definition at line 171 of file generate_icarus_geometry_splitwire.pl.

$TPC_y = $LAr_y

Definition at line 192 of file generate_icarus_geometry_splitwire.pl.

$TPC_z = $LAr_z

Definition at line 193 of file generate_icarus_geometry_splitwire.pl.

$TPCActive_x = $MaxDrift

Definition at line 149 of file generate_icarus_geometry_splitwire.pl.

$TPCActive_y = $CommonWireLength * $SinUAngle + 0.02

Definition at line 150 of file generate_icarus_geometry_splitwire.pl.

$TPCActive_z = $CommonWireLength * $CosUAngle + (4640 - 1) * $UWire_zpitch - 0.000000000002

Definition at line 156 of file generate_icarus_geometry_splitwire.pl.

$TPCActiveHalf_z = $TPCActive_z/2

Definition at line 157 of file generate_icarus_geometry_splitwire.pl.

$TPCHalf_z = $LArHalf_z

Definition at line 194 of file generate_icarus_geometry_splitwire.pl.

$TPCinCryo_x[1] = - $TPC_x/2 - $CPA_x/2

Definition at line 214 of file generate_icarus_geometry_splitwire.pl.

$TPCinCryo_y = -$GaseousAr_y/2

Definition at line 227 of file generate_icarus_geometry_splitwire.pl.

$TPCinCryo_z = 0.

Definition at line 228 of file generate_icarus_geometry_splitwire.pl.

$TPCinCryo_zneg = (-$TPCHalf_z/2)

Definition at line 232 of file generate_icarus_geometry_splitwire.pl.

$TPCinCryo_zpos = ( $TPCHalf_z/2)

Definition at line 233 of file generate_icarus_geometry_splitwire.pl.

$TPCWirePlane_x = 2*$TPCWireThickness

Definition at line 162 of file generate_icarus_geometry_splitwire.pl.

$TPCWirePlane_y = $TPCActive_y

Definition at line 163 of file generate_icarus_geometry_splitwire.pl.

$TPCWirePlane_z = $TPCActive_z

Definition at line 164 of file generate_icarus_geometry_splitwire.pl.

$TPCWirePlaneHalf_z = $TPCActiveHalf_z

Definition at line 165 of file generate_icarus_geometry_splitwire.pl.

$TPCWireRadius = $TPCWireThickness /2

Definition at line 110 of file generate_icarus_geometry_splitwire.pl.

$TPCWireThickness = 0.015

Definition at line 109 of file generate_icarus_geometry_splitwire.pl.

$UAngle = 60

Definition at line 87 of file generate_icarus_geometry_splitwire.pl.

$UpstreamLArPadding = 90.0

Definition at line 132 of file generate_icarus_geometry_splitwire.pl.

$UWire_ypitch = $UWirePitch / $CosUAngle

Definition at line 100 of file generate_icarus_geometry_splitwire.pl.

$UWire_zpitch = $UWirePitch / $SinUAngle

Definition at line 101 of file generate_icarus_geometry_splitwire.pl.

$UWirePitch = 0.3

Definition at line 82 of file generate_icarus_geometry_splitwire.pl.

$VAngle = 60

Definition at line 88 of file generate_icarus_geometry_splitwire.pl.

$VWire_ypitch = $VWirePitch / $CosVAngle

Definition at line 105 of file generate_icarus_geometry_splitwire.pl.

$VWire_zpitch = $VWirePitch / $SinVAngle

Definition at line 106 of file generate_icarus_geometry_splitwire.pl.

$VWirePitch = 0.3

Definition at line 83 of file generate_icarus_geometry_splitwire.pl.

$WarmVessel_CenterToFloor = $WarmVessel_FootHeight + $WarmVessel_FloorSpace + $WarmVessel_y/2

Definition at line 382 of file generate_icarus_geometry_splitwire.pl.

$WarmVessel_FloorSpace = 10.16

Definition at line 381 of file generate_icarus_geometry_splitwire.pl.

$WarmVessel_FootHeight = 6.0

Definition at line 380 of file generate_icarus_geometry_splitwire.pl.

$WarmVessel_x = $ThermIns_x+2*$WarmVesselThickness_x

Definition at line 307 of file generate_icarus_geometry_splitwire.pl.

$WarmVessel_y = $ThermIns_y+2*$WarmVesselThickness

Definition at line 308 of file generate_icarus_geometry_splitwire.pl.

$WarmVessel_z = $ThermIns_z+2*$WarmVesselThickness

Definition at line 309 of file generate_icarus_geometry_splitwire.pl.

$WarmVesselInDetEncl_x = 0

Definition at line 311 of file generate_icarus_geometry_splitwire.pl.

$WarmVesselInDetEncl_y = $ThermInsInDetEncl_y

Definition at line 400 of file generate_icarus_geometry_splitwire.pl.

$WarmVesselInDetEncl_z = 0

Definition at line 313 of file generate_icarus_geometry_splitwire.pl.

$WarmVesselThickness = 27.4

Definition at line 286 of file generate_icarus_geometry_splitwire.pl.

$WarmVesselThickness_x = 21.8

Definition at line 287 of file generate_icarus_geometry_splitwire.pl.

$WirePlaneSpacing = 0.3

Definition at line 120 of file generate_icarus_geometry_splitwire.pl.

$World_x = 1e4

Definition at line 387 of file generate_icarus_geometry_splitwire.pl.

$World_y = 1e4

Definition at line 388 of file generate_icarus_geometry_splitwire.pl.

$World_z = 1e4

Definition at line 389 of file generate_icarus_geometry_splitwire.pl.

print RACETRACK<<EOF; </volume> EOF print RACETRACK<<EOF;<!--+++++++++++++++++ RACETRACK Upstream Volume ++++++++++++++++++++--><volumename="volRaceTrackUVolume"><materialrefref="LAr"/><solidrefref="RaceTrackUVolume"/> EOF $xTpos =-($RaceTrack_width - $RaceTrack_pitch - $RT_epsilon )/2.

Definition at line 1670 of file generate_icarus_geometry_splitwire.pl.

$YWirePitch = 0.3

Definition at line 81 of file generate_icarus_geometry_splitwire.pl.

else

Initial value:

{


    $suffix = "-" . $suffix

Definition at line 39 of file generate_icarus_geometry_splitwire.pl.

*echo Unrecognized option dohelp exit

Definition at line 430 of file generate_icarus_geometry_splitwire.pl.

print PMT<<EOF;<?xml version='1.0'?><gdml> EOF print PMT<<EOF;<!--+++++++++++++++++++ PMT Solids ++++++++++++++++++++++--> EOF print PMT<<EOF;<solids><spherename="PMTVolume"rmin="$PMTradiusInner"rmax="$PMTradiusOuter"deltaphi="360"deltatheta="90"aunit="deg"lunit="cm"/><spherename="PMTPassSphere"rmin="$PMTradiusInner"rmax="$PMTradiusOuter"deltaphi="360"deltatheta="90"aunit="deg"lunit="cm"/><tubename="PMTPassTube"rmax="$PMTtubeRmax"rmin="$PMTtubeRmin"z="$PMTtube_x"deltaphi="360"aunit="deg"lunit="cm"/><unionname="PMTPassVolume"><firstref="PMTPassSphere"/><secondref="PMTPassTube"/><positionrefref="PMTtube"/><rotationrefref="rMinus90AboutZ"/></union></solids> EOF print PMT<<EOF;<structure><!--+++++++++++++++++PMTLogicalVolumes++++++++++++++++++++--><volumename="volOpDetSensitive"><materialrefref="LAr"/><solidrefref="PMTVolume"/></volume><volumename="volNotOpDetSensitive"><materialrefref="Glass"/><solidrefref="PMTPassVolume"/></volume></structure> EOF print PMT<<EOF;</gdml> EOF sub gen_Mech_Structure

Initial value:

{

    $MECH = "icarus_Mech_Structure" . $suffix . ".gdml"

Definition at line 1850 of file generate_icarus_geometry_splitwire.pl.

print RACETRACK<<EOF; </volume> EOF print RACETRACK<<EOF;<!--+++++++++++++++++ RACETRACK end structure +++++++++++++++++++--></structure> EOF print RACETRACK<<EOF;</gdml> EOF sub gen_PMT

Initial value:

{

    $PMT = "icarus_pmt" . $suffix . ".gdml"

Definition at line 1765 of file generate_icarus_geometry_splitwire.pl.

sub gen_RaceTracks

Initial value:

{


#define the
    define_RaceTrack()

Definition at line 1543 of file generate_icarus_geometry_splitwire.pl.