Macros
#define	the solids

Functions
	GetOptions ("input\|i:s"=>\$input,"help\|h"=>\$help,"suffix\|s:s"=>\$suffix,"output\|o:s"=>\$output)

	if (defined $help)

	if (!defined $suffix)

foreach	$parameter ($xmldata->findnodes('/parameters/geometry/parameter'))

	gen_defs ()

	gen_rotations ()

	gen_materials ()

	gen_tpcplanevert ()

	gen_tpcplane ()

	gen_groundplate ()

	gen_cathode ()

	gen_fieldcage ()

	gen_tpc ()

	if ($pmt_switch eq"on")

	gen_cryostat ()

	gen_enclosure ()

	gen_world ()

	write_fragments ()

sub	usage ()

	push (@gdmlFiles, $PMT)

	open (PMT) or die("Could not open file $PMT for writing")

Variables
	else

	$parser = new XML::LibXML

	$xmldata = $parser->parse_file($input)

	$SinUVAngle = sin( deg2rad($UVAngle) )

	$CosUVAngle = cos( deg2rad($UVAngle) )

	$TanUVAngle = tan( deg2rad($UVAngle) )

	$wires_on =1

	$wire_int =20

	$NumberOfTPCPlanes =3

	$tpc_neg_length =1000

	$tpc_neg_width =200

	$tpc_neg_height =200

	$wire_frame_width =9.5

	$wire_plane_height =240

	$wire_plane_length =1042

	$wires_plength =($wire_plane_length - 2*$wire_frame_width)

	$wires_pwidth =($wire_plane_height - 2*$wire_frame_width)

	$field_cage_width = 200

	$field_cage_height = 180

	$field_cage_cross_length = sqrt(($field_cage_width)2+($field_cage_height-50)2)

	$field_cage_length = 1000

	$field_cage_loop_interval = 1

	$spacers_on_off = "off"

	$electronics_height =12

	$pmt_switch ="off"

	exit

sub	gen_pmt

	$PMT = ">" . $PMT

Macro Definition Documentation

#define the solids

Function Documentation

foreach $parameter ( $xmldata->'/parameters/geometry/parameter' )

Definition at line 48 of file generate_gdml_simple.pl.

 {
     # Get the name and value attributes for that parameter:
     $name = $parameter->getAttribute("name");
     $value = $parameter->getAttribute("value");
 
     # Here's the clever part: The following eval creates a variable
     # with the same name as $name. For example, if $name eq "TPCDepth",
     # then the following statement assigns the value to $TPCDepth. The
     # value is in quotes, because some of the parameters have text
     # strings in them (like "kInch").
 
     eval "\$$name = '$value'";
 }

sub gen_cathode ( )

Definition at line 586 of file generate_gdml_simple.pl.

                   {
 
     $CATHODE = "microboone/micro-cathode" . $suffix . ".gdml";
     push (@gdmlFiles, $CATHODE); # Add file to list of GDML fragments
     $CATHODE = ">" . $CATHODE;
     open(CATHODE) or die("Could not open file $CATHODE for writing");
 
     print CATHODE <<EOF;
 
 <?xml version='1.0'?>
 <gdml>
 <solids>
  <box name="CathodePlate"
   lunit="cm"
   x="kCathodePlateDepth"
   y="kCathodeHeight"
   z="kCathodeLength"/>
 </solids>
 <structure>
  <volume name="volCathodePlate">
   <materialref ref="STEEL_STAINLESS_Fe7Cr2Ni"/>
   <solidref ref="CathodePlate"/>
  </volume>
 </structure>
 </gdml>
 EOF
 
 }

gen_cryostat ( )

sub gen_defs ( )

Definition at line 129 of file generate_gdml_simple.pl.

 {
     # Set up the output file.
     $CONSTANTS = "microboone/micro-defs" . $suffix . ".gdml";
     push (@gdmlFiles, $CONSTANTS); # Add file to list of constant files
     $CONSTANTS = ">" . $CONSTANTS;
     open(CONSTANTS) or die("Could not open file $CONSTANTS for writing");
 
     # Create some math constants.
     my $pi = pi;
 
     # Though it's not strictly necessary, make each sub-file a valid
     # XML (if not GDML) document; that makes it accessible to an XML
     # parser if needed.
 
     # Here is a neat way to print out a block of text without getting
     # involved in a lot of messy quoting and formatting with print
     # statements.
 
     print CONSTANTS <<EOF;
 <?xml version='1.0'?>
 <define>
 <constant name="kInch"  value="2.54" />
 <constant name="kPi"    value="$pi" />
 <constant name="kDetEnclosureWidth"       value="$DetEnclosureWidth" />
 <constant name="kDetEnclosureHeight"      value="$DetEnclosureHeight" />
 <constant name="kDetEnclosureLength"      value="$DetEnclosureLength" />
 <constant name="kDirtThickness"           value="$DirtThickness" />
 <constant name="kWorldW"                  value="100.0*kDetEnclosureWidth"/>
 <constant name="kWorldH"                  value="100.0*kDetEnclosureHeight"/>
 <constant name="kWorldL"                  value="100.0*kDetEnclosureLength"/>
 
 <constant name="kTPCWidth"                value="$TPCWidth" />
 <constant name="kTPCLength"               value="$TPCLength" />
 <constant name="kTPCDepth"                value="$TPCDepth" />
 <constant name="kTPCWallThickness"        value="$TPCWallThickness" />
 
 <constant name="kTPCWirePlaneThickness"   value="$TPCWirePlaneThickness" />
 <constant name="kTPCWireThickness"        value="$TPCWireThickness" />
 <constant name="kTPCWirePlaneWidth"       value="$wires_pwidth" />
 <constant name="kTPCWirePlaneLength"      value="$wires_plength" />
 <constant name="kWireFrameDepth"          value="9" /> 
 <constant name="kWireFrameWidth"          value="$wire_frame_width" /> 
 <constant name="kWirePlaneHeight"         value="$wire_plane_height" /> 
 <constant name="kWirePlaneLength"         value="$wire_plane_length" /> 
 <constant name="kWireFrameVInHeight"      value="0.5*(kWirePlaneHeight-3*kWireFrameWidth)" /> 
 
 <constant name="kTPCWirePitch"            value="$TPCWirePitch"/>
 <constant name="kSinUVAngle"              value="$SinUVAngle"/>
 <constant name="kCosUVAngle"              value="$CosUVAngle"/>
 <constant name="kTanUVAngle"              value="$TanUVAngle"/>
 <constant name="kTPCWireXPitch"           value="kTPCWirePitch/kCosUVAngle"/>
 
 <constant name="kCathodeFrameWidth"     value="9" /> 
 <constant name="kCathodeFrameDepth"     value="5" /> 
 <constant name="kCathodePlateDepth"     value="0.1" /> 
 <constant name="kCathodeWidth"          value="5.1" /> 
 <constant name="kCathodeHeight"         value="240" /> 
 <constant name="kCathodeLength"         value="1042" /> 
 <constant name="kCathodeFrameVInHeight"         value="0.5*(kCathodeHeight-3*kCathodeFrameWidth)" /> 
 
 <constant name="kGroundPlateWidth"      value="224" /> 
 <constant name="kGroundPlateHeight"     value="0.1" /> 
 <constant name="kGroundPlateLength"             value="1100" /> 
 <constant name="kGroundBeamWidth"       value="2.5" /> 
 <constant name="kGroundBeamHeight"      value="2.5" /> 
 <constant name="kGroundBeamLength"              value="kGroundPlateLength" />
 <constant name="kGroundBeamThickness"   value=".15" /> 
 </define>
 EOF
 
    close(CONSTANTS);
 }

sub gen_enclosure ( )

Definition at line 998 of file generate_gdml_simple.pl.

 {
     # Set up the output file.
     $GDML = "micro-enclosure" . $suffix . ".gdml";
     push (@gdmlFiles, $GDML); # Add file to list of GDML fragments
     $GDML = ">" . $GDML;
     open(GDML) or die("Could not open file $GDML for writing");
 
     print GDML <<EOF;
 <?xml version='1.0'?>
 <gdml>
 <solids>
  <box name="DetEnclosure" lunit="cm"
    x="kDetEnclosureWidth" y="$CryostatRadius" z="kDetEnclosureLength"
  />
 </solids>
 
 <structure>
  <volume name="volDetEnclosure">
   <materialref ref="Air"/>
   <solidref ref="DetEnclosure"/>
   <physvol>
    <volumeref ref="volCryostat"/>
    <position name="posCryostat" unit="cm" x="0" y="0" z="0"/>
   </physvol>
  </volume>
 </structure>
 </gdml>
 EOF
 
    close(GDML);
 }

sub gen_fieldcage ( )

Definition at line 491 of file generate_gdml_simple.pl.

                     {
 
     # Set up the output file.
     $FIELDCAGE = "microboone/micro-fieldcage.gdml";
     push (@gdmlFiles, $FIELDCAGE); # Add file to list of constant files
     $FIELDCAGE = ">" . $FIELDCAGE;
     open(FIELDCAGE) or die("Could not open file $FIELDCAGE for writing");
 
     # Print the Field Cage constants
     print FIELDCAGE <<EOF;
 <define>
  <constant name="kFieldCageTPCClearance"        value="5*kInch" />
 
  <constant name="kFieldCageTubeRadius"          value="0.5*kInch" />
  <constant name="kFieldCageTubeThickness"       value="0.25*kInch" />
  <constant name="kFieldCageBeamDepth"           value="12.5"/>
  <constant name="kFieldCageBeamWidth"           value="2.5"/>
  <constant name="kFieldCageCrossDepth"          value="2.5"/>
  <constant name="kFieldCageCrossWidth"          value="4"/>
  <constant name="kFieldCageCrossLength"         value="$field_cage_cross_length"/>
 
  <constant name="kTPCTotalLength"               value="$field_cage_length"/>
  <constant name="kTPCTotalWidth"                value="$field_cage_width"/>
  <constant name="kTPCTotalHeight"               value="$field_cage_height"/>
 
  <constant name="kFieldCageLoopLength"          value="kTPCTotalLength+2*(kFieldCageTPCClearance+2*kFieldCageTubeRadius)"/>
  <constant name="kFieldCageLoopWidth"           value="$field_cage_width"/>
  <constant name="kFieldCageLoopHeight"          value="kTPCTotalHeight+2*(kFieldCageTPCClearance+2*kFieldCageTubeRadius)"/>
 
  <constant name="kFieldCageCornerRadius"        value="0.5*kFieldCageTPCClearance"/>
  <constant name="kFieldCageCornerY"             value="(0.5*kFieldCageLoopHeight)-kFieldCageCornerRadius-kFieldCageTubeRadius"/>
  <constant name="kFieldCageCornerZ"             value="(0.5*kFieldCageLoopLength)-kFieldCageCornerRadius-kFieldCageTubeRadius"/>
 
  <constant name="kFieldCageHeight"              value="kFieldCageLoopHeight+2*(kFieldCageBeamDepth+kFieldCageCrossDepth)"/>
  <constant name="kFieldCageLength"              value="kFieldCageLoopLength+2*(kFieldCageBeamDepth+kFieldCageCrossDepth)"/>
  <constant name="kFieldCageWidth"               value="$field_cage_width"/>
 
  <constant name="kFieldCageBeamYInt"            value="0.5*(kFieldCageLoopHeight-50)"/>
  <constant name="kFieldCageBeamZPos"            value="0.5*(kFieldCageLoopLength)"/>
  <constant name="kFieldCageBeamYPos"            value="0.5*(kFieldCageLoopHeight)"/>
  <constant name="kFieldCageBeamZInt"            value="0.5*(kFieldCageLoopLength-40)"/>
 
  <constant name="kFieldCageCrossYPos"           value="0.5*(kFieldCageLoopHeight+kFieldCageCrossDepth)+kFieldCageBeamDepth"/>
  <constant name="kFieldCageCrossZPos"           value="0.5*(kFieldCageLoopLength+kFieldCageCrossDepth)+kFieldCageBeamDepth"/>
 </define>
 EOF
 
     # Prints Field Cage solids
     print FIELDCAGE <<EOF;
 <solids>
  <tube name="FieldCageTubeZ"
   rmin="kFieldCageTubeRadius - kFieldCageTubeThickness"
   rmax="kFieldCageTubeRadius"  
   z="kFieldCageLoopLength-2*(kFieldCageCornerRadius+kFieldCageTubeRadius)" 
   deltaphi="2*kPi" 
   aunit="rad" 
   lunit="cm"/> 
  <tube name="FieldCageTubeY" 
   rmin="kFieldCageTubeRadius - kFieldCageTubeThickness"
   rmax="kFieldCageTubeRadius"  
   z="kFieldCageLoopHeight-2*(kFieldCageCornerRadius+kFieldCageTubeRadius)"  
   deltaphi="2*kPi" 
   aunit="rad" 
   lunit="cm"/> 
 </solids> 
 EOF
 
     # Prints Field Cage tube loop sub-structure
     print FIELDCAGE <<EOF;
 <structure> 
  <volume name="volFieldCageTubeTop"> 
   <materialref ref="STEEL_STAINLESS_Fe7Cr2Ni"/> 
   <solidref ref="FieldCageTubeZ"/> 
  </volume> 
  <volume name="volFieldCageTubeBot"> 
   <materialref ref="STEEL_STAINLESS_Fe7Cr2Ni"/> 
   <solidref ref="FieldCageTubeZ"/> 
  </volume> 
  <volume name="volFieldCageTubeFront"> 
   <materialref ref="STEEL_STAINLESS_Fe7Cr2Ni"/> 
   <solidref ref="FieldCageTubeY"/> 
  </volume> 
  <volume name="volFieldCageTubeBack"> 
   <materialref ref="STEEL_STAINLESS_Fe7Cr2Ni"/> 
   <solidref ref="FieldCageTubeY"/> 
  </volume> 
 </structure>
 EOF
 
 close(FIELDCAGE); 
 
 }

sub gen_groundplate ( )

Definition at line 615 of file generate_gdml_simple.pl.

                       {
 #this subroutine will produce the gdml fragment for ground plate
 
     $GROUNDPLATE = "microboone/micro-groundplate" . $suffix . ".gdml";
     push (@gdmlFiles, $GROUNDPLATE); # Add file to list of GDML fragments
     $GROUNDPLATE = ">" . $GROUNDPLATE;
     open(GROUNDPLATE) or die("Could not open file $GROUNDPLATE for writing");
 
     print GROUNDPLATE <<EOF;
 <?xml version='1.0'?>
 <gdml>
 <solids>
  <box name="GroundPlate"
   lunit="cm"
   x="kGroundPlateWidth"
   y="kGroundPlateHeight"
   z="kGroundPlateLength"/>
 </solids>
 <structure>
  <volume name="volGroundPlate">
   <materialref ref="STEEL_STAINLESS_Fe7Cr2Ni"/>
   <solidref ref="GroundPlate"/>
  </volume>
 </structure>
 </gdml>
 EOF
 
     close(GROUNDPLATE);
 }

sub gen_materials ( )

Definition at line 251 of file generate_gdml_simple.pl.

 {
     # Create the materials file name and open it.
     $MATERIALS = "materials" . $suffix . ".gdml";
     push (@gdmlFiles, $MATERIALS); # Add file to list of GDML fragments
     $MATERIALS = ">" . $MATERIALS;
     open(MATERIALS) or die("Could not open file $MATERIALS for writing");
 
     # Write the standard XML prefix.
     print MATERIALS <<EOF;
 <?xml version='1.0'?>
 EOF
 
     # Go back the DOM structure read in near the beginning of the
     # program. For each <materials /> element (and there'll probably
     # be only one):
     foreach $materials ( $xmldata->findnodes('/parameters/materials') )
     {
         # Convert that element back to text, and write it out.
         print MATERIALS $materials->toString;
     }
 
     close (MATERIALS);
 }

sub gen_rotations ( )

Definition at line 204 of file generate_gdml_simple.pl.

 {
     my $WirePlusRotation = $UVAngle + 90;
     my $WireMinusRotation = $UVAngle - 90;
 
     $ROTATIONS = "microboone/micro-rotations" . $suffix . ".gdml";
     push (@gdmlFiles, $ROTATIONS); # Add file to list of GDML fragments
     $ROTATIONS = ">" . $ROTATIONS;
     open(ROTATIONS) or die("Could not open file $ROTATIONS for writing");
 
     print ROTATIONS <<EOF;
 <?xml version='1.0'?>
 <define>
    <rotation name="rPlus30AboutX"  unit="deg" x="30"  y="0"   z="0"/>
    <rotation name="rPlus60AboutX"  unit="deg" x="60"  y="0"   z="0"/>
    <rotation name="rPlus90AboutX"  unit="deg" x="90"  y="0"   z="0"/>
    <rotation name="rMinus90AboutX"  unit="deg" x="-90"  y="0"   z="0"/>
    <rotation name="rPlusUVAngleAboutX"  unit="deg" x="150" y="0"   z="0"/>
    <rotation name="rPlus150AboutX"      unit="deg" x="150" y="0"   z="0"/>
    <rotation name="rPlus180AboutX"      unit="deg" x="180" y="0"   z="0"/>
    <rotation name="rMinusUVAngleAboutX" unit="deg" x="-30" y="0"   z="0"/>
    <rotation name="rPlus30AboutY"  unit="deg" x="0"   y="30"  z="0"/>
    <rotation name="rPlus60AboutY"  unit="deg" x="0"   y="60"  z="0"/>
    <rotation name="rPlus90AboutY"  unit="deg" x="0"   y="90"  z="0"/>
    <rotation name="rPlus180AboutY" unit="deg" x="0"   y="180" z="0"/>
    <rotation name="rMinus90AboutY" 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="rPlus180AboutZ"      unit="deg" x="0"   y="0"   z="180"/>
    <rotation name="rMinus180AboutZ"     unit="deg" x="0"   y="0"   z="-180"/>
    <rotation name="rMinus90AboutYPlus180AboutZ" unit="deg" x="0" y="-90" z="180"/>
    <rotation name="rMinus90AboutYMinus90AboutZ" unit="deg" x="0" y="-90" z="-90"/>
    <rotation name="rPlus90AboutYPlus180AboutZ" unit="deg" x="0" y="90" z="180"/>
    <rotation name="rMinus90AboutYPlus90AboutZ" unit="deg" x="0" y="-90" z="90"/>
    <rotation name="rPlus90AboutYMinus90AboutZ" unit="deg" x="0" y="90" z="-90"/>
    <rotation name="rPlus90AboutXPlus90AboutZ"  unit="deg" x="90" y="0"   z="90"/>
    <rotation name="rPlus90AboutXPlus180AboutZ" unit="deg" x="90" y="0"   z="180"/>
    <rotation name="rPlus90AboutXMinus90AboutY" unit="deg" x="90" y="-90" z="0"/>
    <rotation name="rPlus90AboutXMinus90AboutZ" unit="deg" x="90" y="0"   z="-90"/>
    <rotation name="rPlus90AboutXPlus90AboutY"  unit="deg"  x="90" y="90" z="0"/>
    <rotation name="rPMTRotation1"  unit="deg" x="90"  y="270"   z="0"/>
 </define>
 EOF
     close (ROTATIONS);
 }

sub gen_tpc ( )

Definition at line 647 of file generate_gdml_simple.pl.

 {
     # Set up the output file.
     $GDML = "microboone/micro-tpc" . $suffix . ".gdml";
     push (@gdmlFiles, $GDML); # Add file to list of GDML fragments
     $GDML = ">" . $GDML;
     open(GDML) or die("Could not open file $GDML for writing");
 
     print GDML <<EOF;
 <?xml version='1.0'?>
 <gdml>
 <solids>
  <box name="TPC"
   lunit="cm"
   x="kTPCDepth+kTPCWallThickness+(2*kTPCWirePlaneThickness)"
   y="kTPCWidth+(2*kTPCWallThickness)+($electronics_height)"
   z="kTPCLength+(2*kTPCWallThickness)"/>
 </solids>
 <structure>
  <volume name="volTPC">
   <materialref ref="LAr"/>
  <solidref ref="TPC"/>
 EOF
 
     # Ground Plane physical volumes
     # Center = (9,124,0)
     $ground_plate_X=9;
     $ground_plate_Y=124;
     $ground_plate_Z=0;
 
     print GDML <<EOF;
   <physvol>
    <volumeref ref="volGroundPlate"/>
    <position name="posGroundPlate" unit="cm" x="$ground_plate_X+0.25" y="$ground_plate_Y-0.5*(kGroundBeamHeight)" z="0"/>
   </physvol>
 EOF
 
 
     # Cathode Plane physical volumes
     # Center = (0.5*(kTPCWidth-kCathodeWidth),0,0)
     $Cathode_X="0.5*(kTPCWidth-kCathodeWidth)";
     $Cathode_plate_offset="0.5*kCathodePlateDepth";
     $Cathode_frame_position="$Cathode_X+($Cathode_plate_offset)";
 
     print GDML <<EOF;
   <physvol>
    <volumeref ref="volCathodePlate"/>
    <position name="posCathodePlate" unit="cm" x="$Cathode_X-0.5*(kCathodeFrameDepth)-0.1" y="0" z="0"/>
   </physvol>
 EOF
 
 
     # Wire Plane physical volumes 
     # Center = ( -0.5*(kTPCWidth-kWireFrameWidth) , 0 , 0 )
     print GDML <<EOF;
   <physvol>
    <volumeref ref="volTPCPlaneVert"/>
    <position name="posTPCPlaneVert" unit="cm" x="-0.5*(kTPCWidth-kWireFrameWidth)-(2*$TPCWirePlaneSpacing)" y="0" z="0"/>
   </physvol>
   <physvol>
    <volumeref ref="volTPCPlane"/>
    <position name="posTPCPlane" unit="cm" x="-0.5*(kTPCWidth-kWireFrameWidth)-$TPCWirePlaneSpacing" y="0" z="0"/>
   </physvol>
   <physvol>
    <volumeref ref="volTPCPlane"/>
    <position name="posTPCPlane2" unit="cm" x="-0.5*(kTPCWidth-kWireFrameWidth)" y="0" z="0"/>
    <rotationref ref="rPlus180AboutY"/>
   </physvol>
 EOF
 
 
 $space=0;
 $i=1;
 while ( $space < ( $field_cage_width / 2 ) ) {
         $xPos=$space+2;
         print GDML <<EOF;
   <physvol>
    <volumeref ref="volFieldCageTubeTop"/>
    <position name="posFieldCageTubeTopA$i" unit="cm" x="$xPos" y="0.5*(kFieldCageLoopHeight-2*kFieldCageTubeRadius)" z="0"/>
   </physvol>
   <physvol>
    <volumeref ref="volFieldCageTubeBot"/>
    <position name="posFieldCageTubeBotA$i" unit="cm" x="$xPos" y="-0.5*(kFieldCageLoopHeight-2*kFieldCageTubeRadius)" z="0"/>
   </physvol>
   <physvol>
    <volumeref ref="volFieldCageTubeFront"/>
    <position name="posFieldCageTubeFrontA$i" unit="cm" x="$xPos" y="0" z="0.5*(kFieldCageLoopLength-2*kFieldCageTubeRadius)"/>
    <rotation name="rFieldCageVertPlusA$i" unit="deg" x="90" y="0" z="0"/>
   </physvol>
   <physvol>
    <volumeref ref="volFieldCageTubeBack"/>
    <position name="posFieldCageTubeBackA$i" unit="cm" x="$xPos" y="0" z="-0.5*(kFieldCageLoopLength-2*kFieldCageTubeRadius)"/>
    <rotation name="rFieldCageVertMinusA$i" unit="deg" x="-90" y="0" z="0"/>
   </physvol>
 EOF
         print GDML <<EOF;
   <physvol>
    <volumeref ref="volFieldCageTubeTop"/>
    <position name="posFieldCageTubeTopB$i" unit="cm" x="-$xPos" y="0.5*(kFieldCageLoopHeight-2*kFieldCageTubeRadius)" z="0"/>
   </physvol>
   <physvol>
    <volumeref ref="volFieldCageTubeBot"/>
    <position name="posFieldCageTubeBotB$i" unit="cm" x="-$xPos" y="-0.5*(kFieldCageLoopHeight-2*kFieldCageTubeRadius)" z="0"/>
   </physvol>
   <physvol>
    <volumeref ref="volFieldCageTubeFront"/>
    <position name="posFieldCageTubeFrontB$i" unit="cm" x="-$xPos" y="0" z="0.5*(kFieldCageLoopLength-2*kFieldCageTubeRadius)"/>
    <rotation name="rFieldCageVertFrontB$i" unit="deg" x="90" y="0" z="0"/>
   </physvol>
   <physvol>
    <volumeref ref="volFieldCageTubeBack"/>
    <position name="posFieldCageTubeBackB$i" unit="cm" x="-$xPos" y="0" z="-0.5*(kFieldCageLoopLength-2*kFieldCageTubeRadius)"/>
    <rotation name="rFieldCageVertBackB$i" unit="deg" x="-90" y="0" z="0"/>
   </physvol>
 EOF
         $space+=4*$field_cage_loop_interval;
         $i++;
 }
 
 
 
     # Closes TPC volume definition space
     print GDML <<EOF;
  </volume> 
 </structure>
 </gdml>
 EOF
 
    close(GDML);
 }

sub gen_tpcplane ( )

Definition at line 351 of file generate_gdml_simple.pl.

 {
 
 #### temporary edits
 ####    - my $NumberWires = $TPCLength / $TPCWirePitch - 1;
 ####    - my $NumberWiresPerEdge = int( $TPCLength / $TPCYWirePitch );
 ####    - my $NumberSideWires = int( $TanUVAngle * $TPCWidth / $TPCYWirePitch );
 ####    - <tube name="TPCWireCommon" rmax="0.5*kTPCWireThickness" z="kTPCWidth/kCosUVAngle" deltaphi="2*kPi" aunit="rad" lunit="cm"/>
 ####    - <box name="TPCPlane"  x="kTPCWirePlaneThickness"  y="kTPCWidth"  z="kTPCLength"  lunit="cm"/>
 
     my $NumberWires = ( $wires_plength ) / $TPCWirePitch - 1;
 
     $GDML = "microboone/micro-plane" . $suffix . ".gdml";
     push (@gdmlFiles, $GDML); # Add file to list of GDML fragments
     $GDML = ">" . $GDML;
     open(GDML) or die("Could not open file $GDML for writing");
 
     # Calculate the number of wire ends on a given y-edge of the plane.
     my $TPCYWirePitch = $TPCWirePitch / $CosUVAngle;
     my $NumberWiresPerEdge = int( ( $wires_plength ) / $TPCYWirePitch );
 
     # How many side wires will be "cut off" by the lower or higher
     # z-edge?
     my $NumberSideWires = int( $TanUVAngle * ( $wires_pwidth ) / $TPCYWirePitch );
 
     # The number of full-length "center" wires.
     my $NumberCenterWires = $NumberWiresPerEdge - $NumberSideWires;
 
     # define the solids
     print GDML <<EOF;
 <?xml version='1.0'?>
 <gdml>
 <solids>
 EOF
 
     # wires on either end of the tpc
     for($i = 0; $i < $NumberSideWires; ++$i)
     {
         print GDML <<EOF;
 <tube name="TPCWire$i"
   rmax="0.5*kTPCWireThickness"
   z="kTPCWireXPitch*($i+1)/kSinUVAngle"
   deltaphi="2*kPi"
   aunit="rad"
   lunit="cm"/>
 EOF
     }
 
     # The solids for the middle wire and the TPC wire plane, and start off the structures.
     print GDML <<EOF;
 <tube name="TPCWireCommon"
   rmax="0.5*kTPCWireThickness"
   z="($wires_pwidth)/kCosUVAngle"
   deltaphi="2*kPi"
   aunit="rad"
   lunit="cm"/>
 <box name="TPCPlane"
   x="kTPCWirePlaneThickness"
   y="($wires_pwidth)"
   z="($wires_plength)"
   lunit="cm"/>
 </solids>
 <structure>
 EOF
  
     # the wires at either end of the plane
     for ($i = 0; $i < $NumberSideWires; ++$i)
     {
         print GDML <<EOF;
     <volume name="volTPCWire$i">
         <materialref ref="Titanium"/>
         <solidref ref="TPCWire$i"/>
     </volume>
 EOF
     }
 
   
     # The wires in the middle of the plane, and the plane itself.
     print GDML <<EOF;
     <volume name="volTPCWireCommon">
       <materialref ref="Titanium"/>
       <solidref ref="TPCWireCommon"/>
     </volume>
     <volume name="volTPCPlane">
       <materialref ref="LAr"/>
       <solidref ref="TPCPlane"/>
 EOF
 
   # the wires at the -z end
   for ($i = 0; $i < $NumberSideWires; ++$i)
   {
       print GDML <<EOF;
     <physvol>
      <volumeref ref="volTPCWire$i"/> 
      <position name="posTPCWire$i" unit="cm" y="(-0.5*kTPCWirePlaneWidth)+0.5*($i+1)*kTPCWireXPitch/kTanUVAngle" z="-0.5*kTPCWirePlaneLength+0.5*kTPCWireXPitch*($i+1)" x="0"/>
      <rotationref ref="rPlusUVAngleAboutX"/>
     </physvol>
 EOF
   }
 
   # The wires in the middle.
   for ($i = 0; $i < $NumberCenterWires - 1 ; ++$i)
   {
       my $j = $NumberSideWires+$i;
       print GDML <<EOF;
     <physvol>
      <volumeref ref="volTPCWireCommon"/>
      <position name="posTPCWire$j" unit="cm" y="0" z="(-0.5*kTPCWirePlaneLength)+kTPCWireXPitch*(0.5*$NumberSideWires + $i+1)" x="0"/>
      <rotationref ref="rPlusUVAngleAboutX"/>
     </physvol>
 EOF
   }
 
   # the wires at the +z end
   for ($i = 0; $i < $NumberSideWires; ++$i)
   {
       my $j = $NumberSideWires-$i-1;
       my $k = $NumberCenterWires+$NumberSideWires+$i;
 
       print GDML <<EOF;
     <physvol>
      <volumeref ref="volTPCWire$j"/>
      <position name="posTPCWire$k" unit="cm" y="0.5*kTPCWirePlaneWidth-0.5*($j+1)*kTPCWireXPitch/kTanUVAngle" z="0.5*kTPCWirePlaneLength-0.5*kTPCWireXPitch*($j+1)" x="0"/>
      <rotationref ref="rPlusUVAngleAboutX"/>
     </physvol>
 EOF
   }
 
       print GDML <<EOF;
   </volume>
 </structure>
 </gdml>
 EOF
 
   close(GDML);
 }

sub gen_tpcplanevert ( )

Definition at line 280 of file generate_gdml_simple.pl.

 {
 
 
 ##### temporary edits:
 #####   - TPCPlaneVert   y="TPCWidth" z="kTPCLength"
 #####   - TPCWireVert    y="TPCWidth"
 #####   - (above)my $NumberWires = int($TPCLength / $TPCWirePitch ) - 1
 
 
 
     my $NumberWires = int( ( $wires_plength ) / $TPCWirePitch ) - 1;
 
     $GDML = "microboone/micro-vertplane" . $suffix . ".gdml";
     push (@gdmlFiles, $GDML); # Add file to list of GDML fragments
     $GDML = ">" . $GDML;
     open(GDML) or die("Could not open file $GDML for writing");
 
     # Define the solids and structures: the wires and the TPC wire plane.
 
     print GDML <<EOF;
 <?xml version='1.0'?>
 <gdml>
 <solids>
 <tube name="TPCWireVert"
   rmax="0.5*kTPCWireThickness"
   z="$wires_pwidth"             
   deltaphi="2*kPi"
   aunit="rad"
   lunit="cm"/>
 <box name="TPCPlaneVert"
   x="kTPCWirePlaneThickness" 
   y="$wires_pwidth" 
   z="$wires_plength"
   lunit="cm"/>
 </solids>
 <structure>
   <volume name="volTPCWireVert">
     <materialref ref="Titanium"/>
     <solidref ref="TPCWireVert"/>
   </volume>
   <volume name="volTPCPlaneVert">
     <materialref ref="LAr"/>       
     <solidref ref="TPCPlaneVert"/>
 EOF
 
     # the wires 
     for ( $i = 0; $i < ( $NumberWires / $wire_int ); ++$i )
     {
         $j=($wire_int*$i);
         print GDML <<EOF;
             <physvol>
              <volumeref ref="volTPCWireVert"/>
              <position name="posTPCWireVert$i" unit="cm" z="(-0.5*kTPCWirePlaneLength)+kTPCWirePitch*($j+1)" x="0" y="0"/>
              <rotationref ref="rPlus90AboutX"/>
             </physvol>
 EOF
     }
 
     print GDML <<EOF;
   </volume>
 </structure>
 </gdml>
 EOF
 
     close(GDML);
 }

sub gen_world ( )

Definition at line 1067 of file generate_gdml_simple.pl.

 {
     # Set up the output file.
     $GDML = "micro-world" . $suffix . ".gdml";
     push (@gdmlFiles, $GDML); # Add file to list of GDML fragments
     $GDML = ">" . $GDML;
     open(GDML) or die("Could not open file $GDML for writing");
 
     print GDML <<EOF;
 <?xml version='1.0'?>
 <gdml>
 <solids>
   <box name="World" 
     lunit="cm" 
     x="kWorldW" 
     y="kWorldH" 
     z="kWorldL"/>
   <tube name="Ground"
     rmin="620*2.54"
     rmax="((50*12)+620)*2.54"
     z="41*12*2.54"
     deltaphi="360" 
     lunit="cm"
     aunit="deg"/>
   <tube name="ConcreteEnclosure"
     rmin="584*2.54"
     rmax="620*2.54"
     z="38*12*2.54"
     deltaphi="360" 
     lunit="cm"
     aunit="deg"/>
   <tube name="ConcreteEnclosureBottom"
     rmin="0"
     rmax="620*2.54"
     z="36*2.54"
     deltaphi="360" 
     lunit="cm"
     aunit="deg"/>
   <tube name="Overburden"
     rmin="0"
     rmax="584*2.54"
     z="10*12*2.54"
     deltaphi="360" 
     lunit="cm"
     aunit="deg"/>
 </solids>
 
 <structure>
   <volume name="volGround" >
     <materialref ref="Dirt" />
     <solidref ref="Ground" />
   </volume>
   <volume name="volOverburden" >
     <materialref ref="Dirt" />
     <solidref ref="Overburden" />
   </volume>
   <volume name="volConcreteEnclosure" >
     <materialref ref="Concrete" />
     <solidref ref="ConcreteEnclosure" />
   </volume>
   <volume name="volConcreteEnclosureBottom" >
     <materialref ref="Concrete" />
     <solidref ref="ConcreteEnclosureBottom" />
   </volume>
   <volume name="volWorld" >
     <materialref ref="Air"/> 
     <solidref ref="World"/>
     <physvol>
       <volumeref ref="volConcreteEnclosure"/>
       <position name="posConcreteEnclosure" unit="cm" x="0.5*kTPCDepth" y="36*2.54/2" z="0.5*kTPCLength"/>
       <rotationref ref="rPlus90AboutX"/>
     </physvol>
     <physvol>
       <volumeref ref="volConcreteEnclosureBottom"/>
       <position name="posConcreteEnclosureBottom" unit="cm" x="0.5*kTPCDepth" y="-38*12*2.54/2" z="0.5*kTPCLength"/>
       <rotationref ref="rPlus90AboutX"/>
     </physvol>
     <physvol>
       <volumeref ref="volGround"/>
       <position name="posGround" unit="cm" x="0.5*kTPCDepth" y="0" z="0.5*kTPCLength"/>
       <rotationref ref="rPlus90AboutX"/>
     </physvol>
     <physvol>
       <volumeref ref="volOverburden"/>
       <position name="posOverburden" unit="cm" x="0.5*kTPCDepth" y="(41-10)*12*2.54/2" z="0.5*kTPCLength"/>
       <rotationref ref="rPlus90AboutX"/>
     </physvol>
     <physvol>
       <volumeref ref="volDetEnclosure"/>
       <position name="posDetEnclosure" unit="cm" x="0.5*kTPCDepth" y="0" z="0.5*kTPCLength"/>
     </physvol>
   </volume> 
 </structure>
 </gdml>
 EOF
 
    close(GDML);
 }

GetOptions	(	"input\|i:s"	,
		\	$input,
		"help\|h"	,
		\	$help,
		"suffix\|s:s"	,
		\	$suffix,
		"output\|o:s"	,
		\	$output
	)

if ( defined $help )

Definition at line 20 of file generate_gdml_simple.pl.

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

if ( !defined $suffix )

Definition at line 27 of file generate_gdml_simple.pl.

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

if ( $pmt_switch eq"on" )

Definition at line 103 of file generate_gdml_simple.pl.

103 { gen_pmt(); } # physical volumes defined in gen_cryostat()

gen_pmt

sub gen_pmt

Definition: generate_gdml_simple.pl:780

gen_cryostat

gen_cryostat()

defined

then echo echo For and will not be changed by echo further linking echo echo B echo The symbol is in the uninitialized data multiple common symbols may appear with the echo same name If the symbol is defined the common echo symbols are treated as undefined references For more echo details on common see the discussion of warn common echo in *Note Linker see the discussion of warn common echo in *Note Linker such as a global int variable echo as opposed to a large global array echo echo I echo The symbol is an indirect reference to another symbol This echo is a GNU extension to the a out object file format which is echo rarely used echo echo N echo The symbol is a debugging symbol echo echo R echo The symbol is in a read only data section echo echo S echo The symbol is in an uninitialized data section for small echo objects echo echo T echo The symbol is in the the normal defined echo symbol is used with no error When a weak undefined symbol echo is linked and the symbol is not defined

Definition: find_global_symbol.sh:107

in

if &&[-z"$BASH_VERSION"] then echo Attempting to switch to bash bash shellSwitch exit fi &&["$1"= 'shellSwitch'] shift declare a IncludeDirectives for Dir in

Definition: quiet_Math_Functor_test.sh:23

open ( PMT )

push	(	@	gdmlFiles,
			$PMT
	)

sub usage ( )

Definition at line 114 of file generate_gdml_simple.pl.

 {
     print "Usage: $0 [-h|--help] -i|--input <parameters-file> [-o|--output <fragments-file>] [-s|--suffix <string>]\n";
     print "       -i/--input can be omitted; <parameters-file> contains geometry and material parameters\n";
     print "       if -o is omitted, output goes to STDOUT; <fragments-file> is input to make_gdml.pl\n";
     print "       -s <string> appends the string to the file names; useful for multiple detector versions\n";
     print "       -h prints this message, then quits\n";
 }

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

$CosUVAngle = cos( deg2rad($UVAngle) )

Definition at line 65 of file generate_gdml_simple.pl.

$electronics_height =12

Definition at line 87 of file generate_gdml_simple.pl.

$field_cage_cross_length = sqrt(($field_cage_width)**2+($field_cage_height-50)**2)

Definition at line 83 of file generate_gdml_simple.pl.

$field_cage_height = 180

Definition at line 82 of file generate_gdml_simple.pl.

$field_cage_length = 1000

Definition at line 84 of file generate_gdml_simple.pl.

$field_cage_loop_interval = 1

Definition at line 85 of file generate_gdml_simple.pl.

$field_cage_width = 200

Definition at line 81 of file generate_gdml_simple.pl.

$NumberOfTPCPlanes =3

Definition at line 72 of file generate_gdml_simple.pl.

$parser = new XML::LibXML

Definition at line 41 of file generate_gdml_simple.pl.

$PMT = ">" . $PMT

Definition at line 784 of file generate_gdml_simple.pl.

$pmt_switch ="off"

Definition at line 88 of file generate_gdml_simple.pl.

$SinUVAngle = sin( deg2rad($UVAngle) )

Definition at line 64 of file generate_gdml_simple.pl.

$spacers_on_off = "off"

Definition at line 86 of file generate_gdml_simple.pl.

$TanUVAngle = tan( deg2rad($UVAngle) )

Definition at line 66 of file generate_gdml_simple.pl.

$tpc_neg_height =200

Definition at line 75 of file generate_gdml_simple.pl.

$tpc_neg_length =1000

Definition at line 73 of file generate_gdml_simple.pl.

$tpc_neg_width =200

Definition at line 74 of file generate_gdml_simple.pl.

$wire_frame_width =9.5

Definition at line 76 of file generate_gdml_simple.pl.

$wire_int =20

Definition at line 71 of file generate_gdml_simple.pl.

$wire_plane_height =240

Definition at line 77 of file generate_gdml_simple.pl.

$wire_plane_length =1042

Definition at line 78 of file generate_gdml_simple.pl.

$wires_on =1

Definition at line 70 of file generate_gdml_simple.pl.

$wires_plength =($wire_plane_length - 2*$wire_frame_width)

Definition at line 79 of file generate_gdml_simple.pl.

$wires_pwidth =($wire_plane_height - 2*$wire_frame_width)

Definition at line 80 of file generate_gdml_simple.pl.

$xmldata = $parser->parse_file($input)

Definition at line 45 of file generate_gdml_simple.pl.

else

Initial value:

{
    $suffix = "-" . $suffix

Definition at line 34 of file generate_gdml_simple.pl.

exit

Definition at line 110 of file generate_gdml_simple.pl.

sub gen_pmt

Initial value:

{
    $PMT = "microboone/micro-pmtdef" . $suffix . ".gdml"

Definition at line 780 of file generate_gdml_simple.pl.

Macros

Functions

Variables

Macro Definition Documentation

Function Documentation

Variable Documentation