createSTTwithPassive.C File Reference

#include <iostream>
#include <iomanip>
#include <fstream>
#include <math.h>
#include <string>
#include <sstream>
#include <vector>

Go to the source code of this file.

Macros
#define	outerDiam   84.4000
#define	innerDiam   30.0000
#define	tubeInnerDiam   1.0000
#define	tubeOuterDiam   1.0060
#define	wireDiam   0.0020
#define	tubeLength   119.6000
#define	endcapInnerDiam   0.9000
#define	endcapOuterDiam   0.9950
#define	endcapLength   0.7200
#define	topcapInnerDiam   0.1000
#define	topcapLength   0.0800
#define	crimppinLength   1.2000
#define	beltInnerDiam   0.3100
#define	beltLength   0.0800
#define	beltDistance   0.2000
#define	gassupplyInnerDiam   0.0500
#define	gassupplyOuterDiam   0.0800
#define	gassupplyLength   3.0000
#define	gassupplyTransY   0.4000
#define	ringInnerDiam   0.9500
#define	ringLength   0.0250
#define	springInnerDiam   0.5000
#define	springOuterDiam   0.5500
#define	springLength   0.1000
#define	noseOuterDiam   0.3000
#define	noseLength   0.5000
#define	safety   0.2000
#define	sttCenterX   0.
#define	sttCenterY   0.
#define	sttCenterZ   20.
#define	innerCoverThickness   1.0000
#define	outerCoverThickness   1.2000
#define	panelthickness   0.1000
#define	pipeDiam   4.0800
#define	noSupport   0
#define	skewangle   3.
#define	pi   3.141592653589793238512808959406186204433
#define	cave   "cave"
#define	sttassembly   "stt01assembly"
#define	innerCylinder   "stt01innerCylinder"
#define	outerCylinder   "stt01outerCylinder"
#define	panel1   "stt01box#1"
#define	panel2   "stt01box#2"
#define	panel3   "stt01box#3"
#define	panel4   "stt01box#4"
#define	air   "air"
#define	AlBe   "carbon"
#define	mylar   "mylar"
#define	GasMixture   "argon"
#define	W   "tungsten"
#define	ABS   "carbon"
#define	PVC   "carbon"
#define	Cu   "copper"
#define	CuBe   "copper"
#define	GFK   "carbon"

Functions
void	writename (char const *name, bool support=false, bool leftside=false)
void	writemother (char const *name, bool original=true, bool support=false, bool firsttube=false)
void	writemedium (char *name)
void	writetube (double inner, double outer, double length)
void	writehalftube (double inner, double outer, double length)
void	writepanel (char const *name, bool firstone, double xthick, double ythick, double length, int side)
void	writetrans (double x, double y, double z)
void	writerot (double x00, double x01, double x02, double x10, double x11, double x12, double x20, double x21, double x22)
bool	putStraw (double posX, double posY, double posZ)
void	plotrotright (double x, double y, double z)
void	plotrotleft (double x, double y, double z)
bool	putStrawRotatedShortLeft (double posX, double posY, double posZ, double xvector, double yvector, double zvector, double length)
bool	putStrawRotatedShortRight (double posX, double posY, double posZ, double xvector, double yvector, double zvector, double length)
bool	putStrawRotatedLeft (double posX, double posY, double posZ, double xvector, double yvector, double zvector, double length)
bool	putStrawRotatedRight (double posX, double posY, double posZ, double xvector, double yvector, double zvector, double length)
void	placeSingleLayerStraightExact (double ringteller)
void	placeSingleLayerSkewedRight (double ringposition)
void	placeSingleLayerSkewedLeft (double ringposition)
void	makeDoubleLayerStraightExact (double &startRadius)
void	makeSingleLayerStraightExact (double &startRadius)
void	makeDoubleLayerSkewedRight (double startRadius)
void	makeDoubleLayerSkewedLeft (double startRadius)
int	main ()

Variables
static int	counter1
static int	counter2
static int	counter3
static int	counter4
static int	counter5
static int	counter6
bool	passive
double	tubeSeperation = tubeOuterDiam
static int	tubeteller
static int	axialtubeteller
static int	skewedtubeteller
static int	shortskewedtubeteller
static double	maximumradius = 0
static double	minimumradius = 100000

Macro Definition Documentation

#define ABS   "carbon"

Definition at line 98 of file createSTTwithPassive.C.

Referenced by geom(), and putStraw().

#define air   "air"

Definition at line 93 of file createSTTwithPassive.C.

Referenced by main(), putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

#define AlBe   "carbon"

Definition at line 94 of file createSTTwithPassive.C.

Referenced by main(), and writepanel().

#define beltDistance   0.2000

Definition at line 52 of file createSTTwithPassive.C.

Referenced by putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

#define beltInnerDiam   0.3100

Definition at line 50 of file createSTTwithPassive.C.

Referenced by putStraw().

#define beltLength   0.0800

Definition at line 51 of file createSTTwithPassive.C.

Referenced by putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

#define cave   "cave"

Definition at line 83 of file createSTTwithPassive.C.

Referenced by main().

#define crimppinLength   1.2000

Definition at line 49 of file createSTTwithPassive.C.

Referenced by putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

#define Cu   "copper"

Definition at line 100 of file createSTTwithPassive.C.

Referenced by hypGeGeoCOSYsetup2014(), hypGeGeoCOSYsetup2014Actives(), hypGeGeoCOSYsetup2014Passives(), hypGeGeoCOSYTarget(), putStraw(), PndSttSingleStraw::STUrban(), and PndFtsSingleStraw::STUrban().

#define CuBe   "copper"

Definition at line 101 of file createSTTwithPassive.C.

Referenced by putStraw().

#define endcapInnerDiam   0.9000

Definition at line 44 of file createSTTwithPassive.C.

Referenced by putStraw().

#define endcapLength   0.7200

Definition at line 46 of file createSTTwithPassive.C.

Referenced by putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

#define endcapOuterDiam   0.9950

Definition at line 45 of file createSTTwithPassive.C.

Referenced by putStraw().

#define GasMixture   "argon"

Definition at line 96 of file createSTTwithPassive.C.

Referenced by putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

#define gassupplyInnerDiam   0.0500

Definition at line 53 of file createSTTwithPassive.C.

Referenced by putStraw().

#define gassupplyLength   3.0000

Definition at line 55 of file createSTTwithPassive.C.

Referenced by putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

#define gassupplyOuterDiam   0.0800

Definition at line 54 of file createSTTwithPassive.C.

Referenced by putStraw().

#define gassupplyTransY   0.4000

Definition at line 56 of file createSTTwithPassive.C.

Referenced by putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

#define GFK   "carbon"

Definition at line 102 of file createSTTwithPassive.C.

Referenced by putStraw().

#define innerCoverThickness   1.0000

Definition at line 72 of file createSTTwithPassive.C.

Referenced by main(), makeSingleLayerStraightExact(), and putStraw().

#define innerCylinder   "stt01innerCylinder"

Definition at line 85 of file createSTTwithPassive.C.

Referenced by main().

#define innerDiam   30.0000

Definition at line 32 of file createSTTwithPassive.C.

Referenced by main(), and putStraw().

#define mylar   "mylar"

Definition at line 95 of file createSTTwithPassive.C.

Referenced by putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

#define noseLength   0.5000

Definition at line 63 of file createSTTwithPassive.C.

Referenced by putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

#define noseOuterDiam   0.3000

Definition at line 62 of file createSTTwithPassive.C.

Referenced by putStraw().

#define noSupport   0

Definition at line 77 of file createSTTwithPassive.C.

Referenced by main().

#define outerCoverThickness   1.2000

Definition at line 73 of file createSTTwithPassive.C.

Referenced by main(), and putStraw().

#define outerCylinder   "stt01outerCylinder"

Definition at line 86 of file createSTTwithPassive.C.

Referenced by main().

#define outerDiam   84.4000

Definition at line 31 of file createSTTwithPassive.C.

Referenced by main(), and putStraw().

#define panel1   "stt01box#1"

Definition at line 87 of file createSTTwithPassive.C.

Referenced by main().

#define panel2   "stt01box#2"

Definition at line 88 of file createSTTwithPassive.C.

Referenced by main().

#define panel3   "stt01box#3"

Definition at line 89 of file createSTTwithPassive.C.

Referenced by main().

#define panel4   "stt01box#4"

Definition at line 90 of file createSTTwithPassive.C.

Referenced by main().

#define panelthickness   0.1000

Definition at line 75 of file createSTTwithPassive.C.

Referenced by main(), placeSingleLayerSkewedLeft(), placeSingleLayerSkewedRight(), and putStraw().

#define pi   3.141592653589793238512808959406186204433

Definition at line 81 of file createSTTwithPassive.C.

#define pipeDiam   4.0800

Definition at line 76 of file createSTTwithPassive.C.

Referenced by main(), placeSingleLayerSkewedLeft(), placeSingleLayerSkewedRight(), putStraw(), and writehalftube().

#define PVC   "carbon"

Definition at line 99 of file createSTTwithPassive.C.

Referenced by hypGeGeoCOSYInBeamStuff(), and putStraw().

#define ringInnerDiam   0.9500

Definition at line 57 of file createSTTwithPassive.C.

Referenced by putStraw().

#define ringLength   0.0250

Definition at line 58 of file createSTTwithPassive.C.

Referenced by putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

#define safety   0.2000

Definition at line 65 of file createSTTwithPassive.C.

Referenced by placeSingleLayerSkewedLeft(), and placeSingleLayerSkewedRight().

#define skewangle   3.

Definition at line 79 of file createSTTwithPassive.C.

Referenced by placeSingleLayerSkewedLeft(), placeSingleLayerSkewedRight(), plotrotleft(), plotrotright(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

#define springInnerDiam   0.5000

Definition at line 59 of file createSTTwithPassive.C.

Referenced by putStraw(), and putStrawRotatedRight().

#define springLength   0.1000

Definition at line 61 of file createSTTwithPassive.C.

Referenced by putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

#define springOuterDiam   0.5500

Definition at line 60 of file createSTTwithPassive.C.

Referenced by putStraw(), and putStrawRotatedRight().

#define sttassembly   "stt01assembly"

Definition at line 84 of file createSTTwithPassive.C.

Referenced by main(), putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), putStrawRotatedShortRight(), and writepanel().

#define sttCenterX   0.

Definition at line 67 of file createSTTwithPassive.C.

#define sttCenterY   0.

Definition at line 68 of file createSTTwithPassive.C.

#define sttCenterZ   20.

Definition at line 70 of file createSTTwithPassive.C.

Referenced by main().

#define topcapInnerDiam   0.1000

Definition at line 47 of file createSTTwithPassive.C.

Referenced by putStraw().

#define topcapLength   0.0800

Definition at line 48 of file createSTTwithPassive.C.

Referenced by putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

#define tubeInnerDiam   1.0000

Definition at line 33 of file createSTTwithPassive.C.

Referenced by putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

#define tubeLength   119.6000

Definition at line 38 of file createSTTwithPassive.C.

Referenced by main(), placeSingleLayerSkewedLeft(), placeSingleLayerSkewedRight(), and putStraw().

#define tubeOuterDiam   1.0060

Definition at line 34 of file createSTTwithPassive.C.

Referenced by putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

#define W   "tungsten"

Definition at line 97 of file createSTTwithPassive.C.

Referenced by putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

#define wireDiam   0.0020

Definition at line 37 of file createSTTwithPassive.C.

Referenced by putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

Function Documentation

int main

(

void

)

Definition at line 2032 of file createSTTwithPassive.C.

References air, AlBe, axialtubeteller, cave, counter1, counter2, counter3, counter4, counter5, counter6, innerCoverThickness, innerCylinder, innerDiam, makeDoubleLayerSkewedLeft(), makeDoubleLayerSkewedRight(), makeDoubleLayerStraightExact(), maximumradius, minimumradius, noSupport, outerCoverThickness, outerCylinder, outerDiam, panel1, panel2, panel3, panel4, panelthickness, passive, pipeDiam, shortskewedtubeteller, skewedtubeteller, sqrt(), sttassembly, sttCenterZ, tubeLength, tubeSeperation, tubeteller, writehalftube(), writemedium(), writemother(), writename(), writepanel(), writerot(), and writetrans().

{
  // reset counters for the numbers of straws
  counter1 = 0;
  counter2 = 0;
  counter3 = 0;
  counter4 = 0;
  counter5 = 0;
  counter6 = 0;
  tubeteller = 0;
  axialtubeteller = 0;
  skewedtubeteller = 0;

  cout << setiosflags(ios::fixed) << setprecision(6);

  //------------------------------------------------------ volumes

  //                                    //
  //         STT Support:               //
  //                                    //

  if(!noSupport)
    {
      // container for all tubes!
      // for X > 0     
      cout << sttassembly << endl;
      cout << cave << endl;
      cout << "ASSEMBLY" << endl;
      cout << air  << endl;
      counter4++;
      writetrans(0., 0., sttCenterZ);
      writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
      // cylinder inner x > 0
      writename(innerCylinder, 1, 1);
      writemother(sttassembly, 1, 1);
      writemedium(AlBe);
      writehalftube((innerDiam / 2.), (innerDiam / 2.) + innerCoverThickness, (tubeLength + 0.4) / 2.);
      writetrans(0., 0., 0.);
      writerot(0., 1., 0., -1., 0., 0., 0., 0., 1.);
      //      writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);

      // cylinder outer x > 0
      writename(outerCylinder, 1, 1);
      writemother(sttassembly, 1, 1);
      writemedium(AlBe);
      writehalftube((outerDiam / 2.) - outerCoverThickness, (outerDiam / 2.), (tubeLength + 0.4) / 2.);
      writetrans(0., 0., 0.);
      writerot(0., 1., 0., -1., 0., 0., 0., 0., 1.);
      //      writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);

      // cylinder inner x < 0
      writename(innerCylinder, 1, 0);
      writemother(sttassembly, 1, 1);
      writemedium(AlBe);
      writehalftube((innerDiam / 2.), (innerDiam / 2.) + innerCoverThickness, (tubeLength + 0.4) / 2.);
      writetrans(0., 0., 0.);
      writerot(0., -1., 0., 1., 0., 0., 0., 0., 1.);
      //      writerot(1., 0., 0., 0., -1., 0., 0., 0., 1.);

      // cylinder outer x < 0
      writename(outerCylinder, 1, 0);
      writemother(sttassembly, 1, 1);
      writemedium(AlBe);
      writehalftube((outerDiam / 2.) - outerCoverThickness, (outerDiam / 2.), (tubeLength + 0.4) / 2.);
      writetrans(0., 0., 0.);
      writerot(0., -1., 0., 1., 0., 0., 0., 0., 1.);
      //      writerot(1., 0., 0., 0., -1., 0., 0., 0., 1.);

      // around the pipe CHECK why (pipeDiam/2.) "-" panelthickness/2?? -----------------
      // panel up x > 0
      writepanel(panel1, true, panelthickness, (outerDiam/2. - innerDiam/2.),  (tubeLength + 0.4), 1);
      writetrans((pipeDiam/2.) + panelthickness/2.,  (outerDiam/2. + innerDiam/2.)/2., 0.);
      writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      // panel up x < 0
      writepanel(panel2,  false, panelthickness, (outerDiam/2. - innerDiam/2.),  (tubeLength + 0.4), 2);
      writetrans(-((pipeDiam/2.) + panelthickness/2.), (outerDiam/2. + innerDiam/2.)/2., 0.);
      writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      // panel down x > 0
      writepanel(panel3, false, panelthickness, (outerDiam/2. - innerDiam/2.),  (tubeLength + 0.4), 1); 
      writetrans((pipeDiam/2.) + panelthickness/2., -(outerDiam/2. + innerDiam/2.)/2., 0.);
      writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      // panel down x < 0
      writepanel(panel4, false, panelthickness, (outerDiam/2. - innerDiam/2.),  (tubeLength + 0.4), 2); 
      writetrans(-((pipeDiam/2.) + panelthickness/2.), -(outerDiam/2. + innerDiam/2.)/2., 0.);
      writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);

    }

  //                                    //
  //  Layers Geometry and               //
  //  Placement of the logical tubes :  //
  //                                    //

  double
    startRadius = innerCoverThickness + innerDiam / 2.; // CHECK added CoverThickness // CHECK

  bool
    wintzdesign = true;

  passive = true;//false;

  if (wintzdesign)
    {

      // NEW ONE!
      // 4 double layers of straight straws
      makeDoubleLayerStraightExact(startRadius);
      makeDoubleLayerStraightExact(startRadius);
      makeDoubleLayerStraightExact(startRadius);
      makeDoubleLayerStraightExact(startRadius);

      cerr << "start radius before skewed " << startRadius << endl;
      // startRadius = (ringteller + 1) * sqrt(3.) * radius;
  
      makeDoubleLayerSkewedLeft((26 * sqrt(3.) + 2) * (tubeSeperation / 2.));
      makeDoubleLayerSkewedRight((27 * sqrt(3.) + 4) * (tubeSeperation / 2.) + 0.244);
      makeDoubleLayerSkewedLeft((28 * sqrt(3.) + 6) * (tubeSeperation / 2.) + 2*0.244);
      makeDoubleLayerSkewedRight((29 * sqrt(3.) + 8) * (tubeSeperation / 2.) + 3*0.244);

      
      cerr << "start radius after " << startRadius << endl;
      startRadius += 10 * tubeSeperation * sqrt(3.) /2.;

      cerr << "start radius " << startRadius << endl;
    
      // then two double layers of straight straws
      makeDoubleLayerStraightExact(startRadius);
      makeDoubleLayerStraightExact(startRadius);

    }
  else
    {
      makeDoubleLayerStraightExact(startRadius);
      startRadius += tubeSeperation;
      makeDoubleLayerSkewedLeft(startRadius);
      startRadius += tubeSeperation + sqrt(3.) * (tubeSeperation / 2.);
      makeDoubleLayerSkewedRight(startRadius);
      startRadius += (tubeSeperation / 2.) + sqrt(3.) * (tubeSeperation / 2.);
      makeDoubleLayerStraightExact(startRadius);
      startRadius += tubeSeperation;
      makeDoubleLayerSkewedLeft(startRadius);
      startRadius += tubeSeperation + sqrt(3.) * (tubeSeperation / 2.);
      makeDoubleLayerSkewedRight(startRadius);
      startRadius += (tubeSeperation / 2.) + sqrt(3.) * (tubeSeperation / 2.);
      makeDoubleLayerStraightExact(startRadius);
      startRadius += tubeSeperation;
      makeDoubleLayerSkewedLeft(startRadius);
      startRadius += tubeSeperation + sqrt(3.) * (tubeSeperation / 2.);
      makeDoubleLayerSkewedRight(startRadius);
      startRadius += (tubeSeperation / 2.) + sqrt(3.) * (tubeSeperation / 2.);
      makeDoubleLayerStraightExact(startRadius);
    }

  // fill the rest of the volume with straight straws
  // straws outside the outerDiam will not be placed
  bool
    fillup = true;
   
  if (fillup)
    {
      makeDoubleLayerStraightExact(startRadius);
      makeDoubleLayerStraightExact(startRadius);
      makeDoubleLayerStraightExact(startRadius);
      makeDoubleLayerStraightExact(startRadius);
      makeDoubleLayerStraightExact(startRadius);
      makeDoubleLayerStraightExact(startRadius);
      makeDoubleLayerStraightExact(startRadius);
      makeDoubleLayerStraightExact(startRadius);
      makeDoubleLayerStraightExact(startRadius);
      makeDoubleLayerStraightExact(startRadius);
    }

  // output of all counters
  cerr << "tubes: " << tubeteller << endl; 
  cerr << "axial tubes: " << axialtubeteller << endl; 
  cerr << "skewed tubes full length: " << skewedtubeteller << endl; 
  cerr << "skewed tubes short: " << shortskewedtubeteller << endl; 
  cerr << "minimum radius: " << minimumradius << endl;
  cerr << "maximum radius: " << maximumradius << endl;

  return 0;
}

void makeDoubleLayerSkewedLeft

(

double

startRadius

)

Definition at line 2022 of file createSTTwithPassive.C.

References placeSingleLayerSkewedLeft(), sqrt(), and tubeSeperation.

{
  // place first skewed left layer at the specified radius
  placeSingleLayerSkewedLeft(startRadius);

  // the straw separation is the separation between the centers of two straws next to each other within one layer. 
  // Close packing then dictates that separation between the layers is sqrt(3) * (separation / 2)
  placeSingleLayerSkewedLeft(startRadius + sqrt(3.) * (tubeSeperation / 2.));
}

void makeDoubleLayerSkewedRight

(

double

startRadius

)

Definition at line 2015 of file createSTTwithPassive.C.

References placeSingleLayerSkewedRight(), sqrt(), and tubeSeperation.

{
  // see comments at makeDoubleLayerSkewedLeft()
  placeSingleLayerSkewedRight(startRadius);
  placeSingleLayerSkewedRight(startRadius + sqrt(3.) * (tubeSeperation / 2.));
}

void makeDoubleLayerStraightExact

(

double &

startRadius

)

Definition at line 1959 of file createSTTwithPassive.C.

References placeSingleLayerStraightExact(), sqrt(), and tubeSeperation.

{
  // close packed straight layers are only possible at specific positions
  // these are indicated by the integer ringteller. ringteller = 0 would be 
  // a single straw at the origin. ringteller = 1 would be the first layer of
  // 6 straws closepacked around that single straw. And so on ...

  // look for the first ring number which is completely outside the startradius 
  // given:
  double
    radius = tubeSeperation / 2.;

  int
    ringteller = 0;
    
  while (ringteller * sqrt(3.) * radius - radius <= startRadius)
    {
      ringteller++;
    }
  cerr << "startRadius " << startRadius << endl;
  // move the startRadius to the next double layer
  startRadius = (ringteller + 1) * sqrt(3.) * radius;
    
  // place the first layer of this double layer
  cerr << "positioning ring at: " << ringteller * sqrt(3.) * radius << endl;
  placeSingleLayerStraightExact(ringteller);
  // place the second layer
  cerr << "positioning ring at: " << (ringteller + 1) * sqrt(3.) * radius << endl;
  placeSingleLayerStraightExact(ringteller + 1);

  cerr << "ringteller " << ringteller << endl;
}

void makeSingleLayerStraightExact

(

double &

startRadius

)

Definition at line 1992 of file createSTTwithPassive.C.

References innerCoverThickness, placeSingleLayerStraightExact(), sqrt(), and tubeSeperation.

{
  // see comments at makeDoubleLayerStraightExact()

  double 
    radius = tubeSeperation / 2.;
 
  int
    ringteller = 0;
  
  while (ringteller * sqrt(3.) * radius - radius + innerCoverThickness <= startRadius) // CHECK
    {
      ringteller++;
    }

  startRadius = (ringteller + 1) * sqrt(3.) * radius;
  
  cerr << "positioning ring at: " << ringteller * sqrt(3.) * radius << endl;
  placeSingleLayerStraightExact(ringteller);
  cerr << "ringteller " << ringteller << " radius " << radius << endl;
 
}

void placeSingleLayerSkewedLeft

(

double

ringposition

)

Definition at line 1547 of file createSTTwithPassive.C.

References cos(), fabs(), i, limit, panelthickness, pi, pipeDiam, putStrawRotatedLeft(), putStrawRotatedShortLeft(), safety, sin(), skewangle, sqrt(), tubeLength, tubeSeperation, and zpos.

{  
  // calculation of some general properties of the layer layout
  // see layerfront.ps and layertop.ps if description is not clear

  // skewed straws become ellipses in the x-y projection
  // the short ellipse radius is still the straw radius, the
  // long ellipse radius (within the straw plane) is given by:

  double 
    radius = tubeSeperation / 2.;

  double
    newradius = radius / cos((skewangle / 180.) * pi);
    
  
  double
    xpos = 0.,  
    ypos = ringposition / cos((30. / 180.) * pi),
    zpos = 0.,
    tmpxpos,
    tmpypos,
    vectorx,
    vectory,
    vectorz;
    
 
  // the available space for a layer at the radius ringposition
  double
    availableSpace = 2. * ringposition * tan((30. / 180.) * pi);
     
  // the number of skewed full-length straws that fit in the available space
  int
    possibleStraws = int((availableSpace - tubeLength * sin((skewangle / 180.) * pi)) / (newradius * 2));

  // the number of short straws that fit on either side of the pack of full-length straws 
  int
    extraStraws = int(availableSpace / (newradius * 2)) - possibleStraws;

  // the amount of space that is left in the layer when considering the sum of all
  // full length and short straws together
  double
    extraspace = (availableSpace - (possibleStraws + extraStraws) * (newradius * 2)) / 2.;
    
  //  cerr << possibleStraws << "  " <<  extraStraws << endl;
    
  double
    skewanglerad = (skewangle / 180.) * pi,
    sixtyrad = (60. / 180.) * pi,
    // translation between the front end of the straw and the center
    translationToLeft = - 0.5 * tubeLength * tan(skewanglerad);
       
  // calculate the lenghts of the different short straws in this layer
  vector<double>
    lengthsShort;
    
  double
    // start at the right edge of the layer
    xpos2 = xpos - availableSpace * sin((60. / 180.) * pi),
    ypos2 = ypos - availableSpace * cos((60. / 180.) * pi),
    // shift to the center of the rightmost straw
    tmpxpos2 = xpos2,
    tmpypos2 = ypos2  - extraspace - newradius,
    // the limit on the left side of the plane beyond which no center of the short straws may be
    limit = ypos2 - availableSpace + newradius + safety + extraspace;

  //   cerr << "xpos2/ypos2 " << xpos2 << " " << ypos2 << endl;
  //   cerr << "tmpxpos2/tmpypos2 " << tmpxpos2 << " " << tmpypos2 << endl;
  

  // shift to the center of of the leftmost full-length straw 
  for(int i = 0; i < possibleStraws; i++) 
    {
      tmpypos2 -= 2 * newradius;
    }
  int extracounter = 0; // how many short tubes can be really placed 
  // calculate the maximum allowable length for this short straw and shift to the next 
  for (int i = 0; i < extraStraws; i++)
    {
     if((tmpypos2 - limit) < 0) break; 
     extracounter++; 
     double
        lengthShort = fabs((tmpypos2 - newradius - limit) / sin(skewangle * (pi / 180.)));
        
      if (lengthShort > tubeLength)
        lengthShort = tubeLength;
        
      //      cerr << "length: " << lengthShort << endl;
      lengthsShort.push_back(lengthShort);
        
      tmpypos2 -= 2 * newradius;
    } 
  extraStraws = extracounter; // update the number 
  // =================================================    
  // planes intersecting pipe ========================
  // pipe
  double pipespace = (panelthickness + pipeDiam/2.) * 1./(sqrt(3)/2.);  

  // the available space for a layer at the radius ringposition
  double
    availableSpacePipe = 2. * ringposition * tan((30. / 180.) * pi) - pipespace;
  
  // the number of skewed full-length straws that fit in the available space
  int
    possibleStrawsPipe = int((availableSpacePipe - tubeLength * sin((skewangle / 180.) * pi)) / (newradius * 2));
  
  // the number of short straws that fit on either side of the pack of full-length straws 
  int
    extraStrawsPipe = int(availableSpacePipe / (newradius * 2)) - possibleStrawsPipe;
  
  // the amount of space that is left in the layer when considering the sum of all
  // full length and short straws together
  double
    extraspacePipe = (availableSpacePipe - (possibleStrawsPipe + extraStrawsPipe) * (newradius * 2)) / 2.;
    
  //  cerr << "layer with pipe " << possibleStraws << "  " <<  extraStraws << endl;
    
  // calculate the lenghts of the different short straws in this layer
  vector<double>
    lengthsShortPipe;
    
  double
    // start at the right edge of the layer
    xpos2Pipe = xpos - availableSpacePipe * sin((60. / 180.) * pi),
    ypos2Pipe = ypos - availableSpacePipe * cos((60. / 180.) * pi),
    // shift to the center of the rightmost straw
    tmpxpos2Pipe = xpos2Pipe, 
    tmpypos2Pipe = ypos2Pipe - extraspacePipe - newradius,
    // the limit on the left side of the plane beyond which no center of the short straws may be
    limitPipe = ypos2Pipe - availableSpacePipe + newradius + safety + extraspacePipe;

  // shift to the center of of the leftmost full-length straw 
  for(int i = 0; i < possibleStrawsPipe; i++) 
    {
      tmpypos2Pipe -= 2 * newradius;
    }
  
  // calculate the maximum allowable length for this short straw and shift to the next 
  extracounter = 0; // how many short tubes can be really placed 
  for (int i = 0; i < extraStrawsPipe; i++)
    {
      if((tmpypos2Pipe - limitPipe) < 0) break; 
      extracounter++;
      double
        lengthShort = fabs((tmpypos2Pipe - newradius - limitPipe) / sin(skewangle * (pi / 180.)));
        
      if (lengthShort > tubeLength)
        lengthShort = tubeLength;
        
      lengthsShortPipe.push_back(lengthShort);
        
      tmpypos2Pipe -= 2 * newradius;
    } 
  extraStrawsPipe = extracounter; // update the number 
  // =================================================
  // plane 1
  // place full straws:
  vectorx = -cos((60. / 180.) * pi);
  vectory = sin((60. / 180.) * pi);
  vectorz = 0;
  tmpxpos = xpos - (extraspacePipe + newradius + pipespace) * sin((60. / 180.) * pi);
  tmpypos = ypos - (extraspacePipe + newradius + pipespace) * cos((60. / 180.) * pi);
    
  for(int i = 0; i < possibleStrawsPipe; i++) 
    {
      // place a straw
      putStrawRotatedLeft(tmpxpos + translationToLeft * sin(sixtyrad), 
                          tmpypos + translationToLeft * cos(sixtyrad), 
                          zpos, vectorx, vectory, vectorz, tubeLength);
      
      tmpxpos -= 2 * newradius * sin((60. / 180.) * pi);
      tmpypos -= 2 * newradius * cos((60. / 180.) * pi);
    }
    
  for (int i = 0; i < extraStrawsPipe; i++)
    {
      double
        lengthShort = lengthsShortPipe[i];
        
      double
        translationToFrontPar = ((tubeLength - lengthShort) / 2.) * cos(skewanglerad),
        translationToFrontPer = ((tubeLength - lengthShort) / 2.) * sin(skewanglerad),
        // translation to the other side of the straw layer
        switchSides = 2 * newradius * (possibleStrawsPipe + 1 + (2 * i)) * cos(skewanglerad); 
        
        
      putStrawRotatedShortLeft(tmpxpos + (translationToLeft + translationToFrontPer) * sin(sixtyrad), 
                               tmpypos + (translationToLeft + translationToFrontPer) * cos(sixtyrad), 
                               zpos + translationToFrontPar, 
                               vectorx, vectory, vectorz, lengthShort);
        
      putStrawRotatedShortLeft(tmpxpos + (translationToLeft - translationToFrontPer + switchSides) * sin(sixtyrad), 
                               tmpypos + (translationToLeft - translationToFrontPer + switchSides) * cos(sixtyrad),
                               zpos - translationToFrontPar, 
                               vectorx, vectory, vectorz, lengthShort);
    
        
      tmpxpos -= 2 * newradius * sin((60. / 180.) * pi);
      tmpypos -= 2 * newradius * cos((60. / 180.) * pi);
    }

  // plane 2
  vectorx = -1.;
  vectory = 0.;
  vectorz = 0.;
  xpos -= availableSpace * sin((60. / 180.) * pi);
  ypos -= availableSpace * cos((60. / 180.) * pi);
  tmpxpos = xpos;
  tmpypos = ypos - extraspace - newradius;

  for(int i = 0; i < possibleStraws; i++) 
    {
      putStrawRotatedLeft(tmpxpos, tmpypos + translationToLeft, zpos, vectorx, vectory, vectorz, tubeLength);
      tmpypos -= 2 * newradius;
    }
    
  for (int i = 0; i < extraStraws; i++)
    {
      double
        lengthShort = lengthsShort[i];
        
      double
        translationToFrontPar = ((tubeLength - lengthShort) / 2.) * cos(skewanglerad),
        translationToFrontPer = ((tubeLength - lengthShort) / 2.) * sin(skewanglerad),
        switchSides = 2 * newradius * (possibleStraws + 1 + (2 * i)) * cos(skewanglerad); 
        
      putStrawRotatedShortLeft(tmpxpos,
                               tmpypos + (translationToLeft + translationToFrontPer), 
                               zpos + translationToFrontPar, 
                               vectorx, vectory, vectorz, lengthShort);
        
      putStrawRotatedShortLeft(tmpxpos,
                               tmpypos + (translationToLeft - translationToFrontPer + switchSides), 
                               zpos - translationToFrontPar, 
                               vectorx, vectory, vectorz, lengthShort);
        
      tmpypos -= 2 * newradius;
    }
    
  // plane 3
  vectorx =  -1 * cos((60. / 180.) * pi);
  vectory =  -1 * sin((60. / 180.) * pi);
  vectorz = 0;
  ypos -= availableSpace;
  tmpxpos = xpos + (extraspacePipe + newradius) * sin((60. / 180.) * pi);
  tmpypos = ypos - (extraspacePipe + newradius) * cos((60. / 180.) * pi);

  for(int i = 0; i < possibleStrawsPipe; i++) 
    {
        
      putStrawRotatedLeft(tmpxpos - translationToLeft * sin(sixtyrad), 
                          tmpypos + translationToLeft * cos(sixtyrad), 
                          zpos, vectorx, vectory, vectorz, tubeLength);
      
      tmpxpos += 2 * newradius * sin((60. / 180.) * pi);
      tmpypos -= 2 * newradius * cos((60. / 180.) * pi);
    }

  for (int i = 0; i < extraStrawsPipe; i++)
    {
      double
        lengthShort = lengthsShortPipe[i];
        
      double
        translationToFrontPar = ((tubeLength - lengthShort) / 2.) * cos(skewanglerad),
        translationToFrontPer = ((tubeLength - lengthShort) / 2.) * sin(skewanglerad),
        switchSides = 2 * newradius * (possibleStrawsPipe + 1 + (2 * i)) * cos(skewanglerad); 
        
      putStrawRotatedShortLeft(tmpxpos - (translationToLeft + translationToFrontPer) * sin(sixtyrad), 
                               tmpypos + (translationToLeft + translationToFrontPer) * cos(sixtyrad), 
                               zpos + translationToFrontPar, 
                               vectorx, vectory, vectorz, lengthShort);

      putStrawRotatedShortLeft(tmpxpos - (translationToLeft - translationToFrontPer + switchSides) * sin(sixtyrad), 
                               tmpypos + (translationToLeft - translationToFrontPer + switchSides) * cos(sixtyrad), 
                               zpos - translationToFrontPar, 
                               vectorx, vectory, vectorz, lengthShort);
        
      tmpxpos += 2 * newradius * sin((60. / 180.) * pi);
      tmpypos -= 2 * newradius * cos((60. / 180.) * pi);
    }
    
  // plane 4
  vectorx = 1 * cos((60. / 180.) * pi);
  vectory = -1 * sin((60. / 180.) * pi);
  vectorz = 0;
  xpos += availableSpace * sin((60. / 180.) * pi);
  ypos -= availableSpace * cos((60. / 180.) * pi);
  tmpxpos = xpos + (extraspacePipe + newradius + pipespace) * sin((60. / 180.) * pi);
  tmpypos = ypos + (extraspacePipe + newradius + pipespace) * cos((60. / 180.) * pi);
    
  for(int i = 0; i < possibleStrawsPipe; i++) 
    {
        
      putStrawRotatedLeft(tmpxpos - translationToLeft * sin(sixtyrad), 
                          tmpypos - translationToLeft * cos(sixtyrad), 
                          zpos, vectorx, vectory, vectorz, tubeLength);
        
      tmpxpos += 2 * newradius * sin((60. / 180.) * pi);
      tmpypos += 2 * newradius * cos((60. / 180.) * pi);
    }
    
  for (int i = 0; i < extraStrawsPipe; i++)
    {
      double
        lengthShort = lengthsShortPipe[i];
        
      double
        translationToFrontPar = ((tubeLength - lengthShort) / 2.) * cos(skewanglerad),
        translationToFrontPer = ((tubeLength - lengthShort) / 2.) * sin(skewanglerad),
        switchSides = 2 * newradius * (possibleStrawsPipe + 1 + (2 * i)) * cos(skewanglerad); 

        
      putStrawRotatedShortLeft(tmpxpos - (translationToLeft + translationToFrontPer) * sin(sixtyrad), 
                               tmpypos - (translationToLeft + translationToFrontPer) * cos(sixtyrad), 
                               zpos + translationToFrontPar, 
                               vectorx, vectory, vectorz, lengthShort);
      
      putStrawRotatedShortLeft(tmpxpos - (translationToLeft - translationToFrontPer + switchSides) * sin(sixtyrad), 
                               tmpypos - (translationToLeft - translationToFrontPer + switchSides) * cos(sixtyrad), 
                               zpos - translationToFrontPar, 
                               vectorx, vectory, vectorz, lengthShort);

      tmpxpos += 2 * newradius * sin((60. / 180.) * pi);
      tmpypos += 2 * newradius * cos((60. / 180.) * pi);
    }
    
  // plane 5
  vectorx = 1.;
  vectory = 0.;
  vectorz = 0.;
  xpos += availableSpace * sin((60. / 180.) * pi);
  ypos += availableSpace * cos((60. / 180.) * pi);
  tmpxpos = xpos;
  tmpypos = ypos + extraspace + newradius;

  for(int i = 0; i < possibleStraws; i++) 
    {

      putStrawRotatedLeft(tmpxpos, tmpypos - translationToLeft, zpos, vectorx, vectory, vectorz, tubeLength);

      tmpypos += 2 * newradius;
    }
    
  for (int i = 0; i < extraStraws; i++)
    {
      double
        lengthShort = lengthsShort[i];
        
      double
        translationToFrontPar = ((tubeLength - lengthShort) / 2.) * cos(skewanglerad),
        translationToFrontPer = ((tubeLength - lengthShort) / 2.) * sin(skewanglerad),
        switchSides = 2 * newradius * (possibleStraws + 1 + (2 * i)) * cos(skewanglerad); 
        
      putStrawRotatedShortLeft(tmpxpos, 
                               tmpypos - (translationToLeft + translationToFrontPer), 
                               zpos + translationToFrontPar, 
                               vectorx, vectory, vectorz, lengthShort);
        
      putStrawRotatedShortLeft(tmpxpos,
                               tmpypos - (translationToLeft - translationToFrontPer + switchSides), 
                               zpos - translationToFrontPar, 
                               vectorx, vectory, vectorz, lengthShort);

      tmpypos += 2 * newradius;
    }
    
  // plane 6
  vectorx = cos((60. / 180.) * pi);
  vectory = sin((60. / 180.) * pi);
  vectorz = 0;
  ypos += availableSpace;
  tmpxpos = xpos - (extraspacePipe + newradius) * sin((60. / 180.) * pi);
  tmpypos = ypos + (extraspacePipe + newradius) * cos((60. / 180.) * pi);

  for(int i = 0; i < possibleStrawsPipe; i++) 
    {

      putStrawRotatedLeft(tmpxpos + translationToLeft * sin(sixtyrad), 
                          tmpypos - translationToLeft * cos(sixtyrad), 
                          zpos, vectorx, vectory, vectorz, tubeLength);

      tmpxpos -= 2 * newradius * sin((60. / 180.) * pi);
      tmpypos += 2 * newradius * cos((60. / 180.) * pi);
    }
    
  for (int i = 0; i < extraStrawsPipe; i++)
    {
      double
        lengthShort = lengthsShortPipe[i];
        
      double
        translationToFrontPar = ((tubeLength - lengthShort) / 2.) * cos(skewanglerad),
        translationToFrontPer = ((tubeLength - lengthShort) / 2.) * sin(skewanglerad),
        switchSides = 2 * newradius * (possibleStrawsPipe + 1 + (2 * i)) * cos(skewanglerad); 
        
      putStrawRotatedShortLeft(tmpxpos + (translationToLeft + translationToFrontPer) * sin(sixtyrad), 
                               tmpypos - (translationToLeft + translationToFrontPer) * cos(sixtyrad), 
                               zpos + translationToFrontPar, 
                               vectorx, vectory, vectorz, lengthShort);
        
      putStrawRotatedShortLeft(tmpxpos + (translationToLeft - translationToFrontPer + switchSides) * sin(sixtyrad), 
                               tmpypos - (translationToLeft - translationToFrontPer + switchSides) * cos(sixtyrad), 
                               zpos - translationToFrontPar, 
                               vectorx, vectory, vectorz, lengthShort);
        
      tmpxpos -= 2 * newradius * sin((60. / 180.) * pi);
      tmpypos += 2 * newradius * cos((60. / 180.) * pi);
    }
}

void placeSingleLayerSkewedRight

(

double

ringposition

)

Definition at line 1159 of file createSTTwithPassive.C.

References cos(), fabs(), i, limit, panelthickness, pi, pipeDiam, putStrawRotatedRight(), putStrawRotatedShortRight(), safety, sin(), skewangle, sqrt(), tubeLength, tubeSeperation, and zpos.

{  
  // see void placeSingleLayerSkewedLeft(double ringposition)
    
  double radius = tubeSeperation / 2.;
    
  double
    newskewangle = skewangle * -1.;
    
  double
    newradius = radius / cos((newskewangle / 180.) * pi);
    
  double
    xpos = 0., 
    ypos = ringposition / cos((30. / 180.) * pi),  
    zpos = 0.,
    tmpxpos,
    tmpypos,
    vectorx,
    vectory,
    vectorz;
 
  double
    availableSpace = 2. * ringposition * tan((30. / 180.) * pi);
  
  int
    possibleStraws = int((availableSpace - fabs(tubeLength * sin((newskewangle / 180.) * pi))) / (newradius * 2)),
    extraStraws = -1 * possibleStraws + int(availableSpace / (newradius * 2));

  double
    extraspace = (availableSpace - (possibleStraws + extraStraws) * (newradius * 2)) / 2.;

  cerr << "extra space: " << extraspace * 2 << endl;

  double
    skewanglerad = (newskewangle / 180.) * pi,
    sixtyrad = (60. / 180.) * pi,
    translationToLeft = - 0.5 * tubeLength * tan(skewanglerad);


  vector<double>
    lengthsShort;

  double
    xpos2 = xpos - availableSpace * sin((60. / 180.) * pi),
    ypos2 = ypos - availableSpace * cos((60. / 180.) * pi),
    tmpxpos2 = xpos2,
    tmpypos2 = ypos2  - extraspace - newradius,
    limit = ypos2 - availableSpace + newradius + safety + extraspace; 

  for(int i = 0; i < possibleStraws; i++) 
    {
      tmpypos2 -= 2 * newradius;
    }
  int extracounter = 0; // how many short tubes can be really placed 
  for (int i = 0; i < extraStraws; i++)
    {
      if((tmpypos2 - limit) < 0) break; 
      extracounter++; 
      double
        lengthShort = fabs((tmpypos2 - newradius - limit) / sin(newskewangle * (pi / 180.)));

      if (lengthShort > tubeLength)
        lengthShort = tubeLength;

      //      cerr << "length: " << lengthShort << endl;
      lengthsShort.push_back(fabs(lengthShort));

      tmpypos2 -= 2 * newradius;
    }
  extraStraws = extracounter; // update the number 

  double
    additionalShift = fabs(tubeLength * sin((newskewangle / 180.) * pi));
  
  // =================================================    
  // planes intersecting pipe ========================
  // pipe
  double pipespace = (panelthickness + pipeDiam/2.) * 1./(sqrt(3)/2.);  

  // the available space for a layer at the radius ringposition
  double
    availableSpacePipe = 2. * ringposition * tan((30. / 180.) * pi) - pipespace;
  
  // the number of skewed full-length straws that fit in the available space
  int
    possibleStrawsPipe = int((availableSpacePipe - fabs(tubeLength * sin((skewangle / 180.) * pi))) / (newradius * 2));
  
  // the number of short straws that fit on either side of the pack of full-length straws 
  int
    extraStrawsPipe = -1 * possibleStrawsPipe + int(availableSpacePipe / (newradius * 2));
  // the amount of space that is left in the layer when considering the sum of all
  // full length and short straws together
  double
    extraspacePipe = (availableSpacePipe - (possibleStrawsPipe + extraStrawsPipe) * (newradius * 2)) / 2.;
    
  //  cerr << "layer with pipe " << possibleStraws << "  " <<  extraStraws << endl;
    
  // calculate the lenghts of the different short straws in this layer
  vector<double>
    lengthsShortPipe;
    
  double
    // start at the right edge of the layer
    xpos2Pipe = xpos - availableSpacePipe * sin((60. / 180.) * pi),
    ypos2Pipe = ypos - availableSpacePipe * cos((60. / 180.) * pi),
    // shift to the center of the rightmost straw
    tmpxpos2Pipe = xpos2Pipe,
    tmpypos2Pipe = ypos2Pipe - extraspacePipe - newradius,
    // the limit on the left side of the plane beyond which no center of the short straws may be
    limitPipe = ypos2Pipe - availableSpacePipe + newradius + safety + extraspacePipe;


  // shift to the center of of the leftmost full-length straw 
  for(int i = 0; i < possibleStrawsPipe; i++) 
    {
      tmpypos2Pipe -= 2 * newradius;
    }
  extracounter = 0; // how many short tubes can be really placed 
  // calculate the maximum allowable length for this short straw and shift to the next 
  for (int i = 0; i < extraStrawsPipe; i++)
    {
      if((tmpypos2Pipe - limitPipe) < 0) break; 
      extracounter++;
      double
        lengthShort = fabs((tmpypos2Pipe - newradius - limitPipe) / sin(skewangle * (pi / 180.)));
        
      if (lengthShort > tubeLength)
        lengthShort = tubeLength;
      lengthsShortPipe.push_back(lengthShort);
      tmpypos2Pipe -= 2 * newradius;
    }
  extraStrawsPipe = extracounter; // update the number 
  // =================================================

  // plane 1
  vectorx = -cos((60. / 180.) * pi);
  vectory = sin((60. / 180.) * pi);
  vectorz = 0;
  tmpxpos = xpos - (extraspacePipe + newradius + pipespace) * sin((60. / 180.) * pi);
  tmpypos = ypos - (extraspacePipe + newradius + pipespace) * cos((60. / 180.) * pi);
  
  for(int i = 0; i < possibleStrawsPipe; i++) 
    {
      
      putStrawRotatedRight(tmpxpos + (translationToLeft - additionalShift) * sin(sixtyrad),
                           tmpypos + (translationToLeft - additionalShift) * cos(sixtyrad), 
                           zpos, vectorx, vectory, vectorz, tubeLength);
      
      tmpxpos -= 2 * newradius * sin((60. / 180.) * pi);
      tmpypos -= 2 * newradius * cos((60. / 180.) * pi);
    }
  
  for (int i = 0; i < extraStrawsPipe; i++)
    {
      double
        lengthShort = lengthsShortPipe[i];

      double
        translationToFrontPar = ((tubeLength - lengthShort) / 2.) * cos(skewanglerad),
        translationToFrontPer = ((tubeLength - lengthShort) / 2.) * sin(skewanglerad),
        switchSides = 2 * newradius * (possibleStrawsPipe + 1 + (2 * i)) * cos(skewanglerad); 

      putStrawRotatedShortRight(tmpxpos + (translationToLeft - additionalShift - translationToFrontPer) * sin(sixtyrad), 
                                tmpypos + (translationToLeft - additionalShift - translationToFrontPer) * cos(sixtyrad), 
                                zpos - translationToFrontPar, 
                                vectorx, vectory, vectorz, lengthShort);

      putStrawRotatedShortRight(tmpxpos + (translationToLeft - additionalShift + translationToFrontPer + switchSides) * sin(sixtyrad), 
                                tmpypos + (translationToLeft - additionalShift + translationToFrontPer + switchSides) * cos(sixtyrad), 
                                zpos + translationToFrontPar, 
                                vectorx, vectory, vectorz, lengthShort);
                           

      tmpxpos -= 2 * newradius * sin((60. / 180.) * pi);
      tmpypos -= 2 * newradius * cos((60. / 180.) * pi);
    }
  

  // plane 2
  vectorx = -1.;
  vectory = 0.;
  vectorz = 0.;
  xpos -= availableSpace * sin((60. / 180.) * pi);
  ypos -= availableSpace * cos((60. / 180.) * pi);
  tmpxpos = xpos;
  tmpypos = ypos - extraspace - newradius;

  for(int i = 0; i < possibleStraws; i++) 
    {
      putStrawRotatedRight(tmpxpos, tmpypos - additionalShift + translationToLeft, zpos, vectorx, vectory, vectorz, tubeLength);
      tmpypos -= 2 * newradius;
    }
  
  for (int i = 0; i < extraStraws; i++)
    {
      double
        lengthShort = lengthsShort[i];
     
      double
        translationToFrontPar = ((tubeLength - lengthShort) / 2.) * cos(skewanglerad),
        translationToFrontPer = ((tubeLength - lengthShort) / 2.) * sin(skewanglerad),
        switchSides = 2 * newradius * (possibleStraws + 1 + (2 * i)) * cos(skewanglerad); 
     
      putStrawRotatedShortRight(tmpxpos, 
                                tmpypos - additionalShift + (translationToLeft - translationToFrontPer), 
                                zpos - translationToFrontPar, 
                                vectorx, vectory, vectorz, lengthShort);
      
      putStrawRotatedShortRight(tmpxpos, 
                                tmpypos - additionalShift + (translationToLeft + translationToFrontPer + switchSides), 
                                zpos + translationToFrontPar, 
                                vectorx, vectory, vectorz, lengthShort);
      
      tmpypos -= 2 * newradius;
    }
  
  // plane 3
  vectorx = -1 * cos((60. / 180.) * pi);
  vectory =  -1 * sin((60. / 180.) * pi);
  vectorz = 0;
  ypos -= availableSpace;
  tmpxpos = xpos + (extraspacePipe + newradius) * sin((60. / 180.) * pi);
  tmpypos = ypos - (extraspacePipe + newradius) * cos((60. / 180.) * pi);

  for(int i = 0; i < possibleStrawsPipe; i++) 
    {
      
      putStrawRotatedRight(tmpxpos - (translationToLeft - additionalShift) * sin(sixtyrad),
                           tmpypos + (translationToLeft - additionalShift) * cos(sixtyrad), 
                           zpos, vectorx, vectory, vectorz, tubeLength);
        
      tmpxpos += 2 * newradius * sin((60. / 180.) * pi);
      tmpypos -= 2 * newradius * cos((60. / 180.) * pi);
    }
  
  for (int i = 0; i < extraStrawsPipe; i++)
    {
      double
        lengthShort = lengthsShortPipe[i];

      double
        translationToFrontPar = ((tubeLength - lengthShort) / 2.) * cos(skewanglerad),
        translationToFrontPer = ((tubeLength - lengthShort) / 2.) * sin(skewanglerad),
        switchSides = 2 * newradius * (possibleStrawsPipe + 1 + (2 * i)) * cos(skewanglerad); 
            
      putStrawRotatedShortRight(tmpxpos - (translationToLeft - translationToFrontPer - additionalShift) * sin(sixtyrad), 
                                tmpypos + (translationToLeft - translationToFrontPer - additionalShift) * cos(sixtyrad), 
                                zpos - translationToFrontPar, 
                                vectorx, vectory, vectorz, lengthShort);
 
      putStrawRotatedShortRight(tmpxpos - (translationToLeft + translationToFrontPer - additionalShift + switchSides) * sin(sixtyrad), 
                                tmpypos + (translationToLeft + translationToFrontPer - additionalShift + switchSides) * cos(sixtyrad), 
                                zpos + translationToFrontPar, 
                                vectorx, vectory, vectorz, lengthShort);
      
      tmpxpos += 2 * newradius * sin((60. / 180.) * pi);
      tmpypos -= 2 * newradius * cos((60. / 180.) * pi);
    }
  
  // plane 4
  vectorx = 1 * cos((60. / 180.) * pi);
  vectory = -1 * sin((60. / 180.) * pi);
  vectorz = 0;
  xpos += availableSpace * sin((60. / 180.) * pi);
  ypos -= availableSpace * cos((60. / 180.) * pi);
  tmpxpos = xpos + (extraspacePipe + newradius + pipespace) * sin((60. / 180.) * pi);
  tmpypos = ypos + (extraspacePipe + newradius + pipespace) * cos((60. / 180.) * pi);
    
  for(int i = 0; i < possibleStrawsPipe; i++) 
    {
      
      putStrawRotatedRight(tmpxpos - (translationToLeft - additionalShift) * sin(sixtyrad), 
                           tmpypos - (translationToLeft - additionalShift) * cos(sixtyrad), 
                           zpos, vectorx, vectory, vectorz, tubeLength);
      
      tmpxpos += 2 * newradius * sin((60. / 180.) * pi);
      tmpypos += 2 * newradius * cos((60. / 180.) * pi);
    }
  
  for (int i = 0; i < extraStrawsPipe; i++)
    {
      double
        lengthShort = lengthsShortPipe[i];

      double
        translationToFrontPar = ((tubeLength - lengthShort) / 2.) * cos(skewanglerad),
        translationToFrontPer = ((tubeLength - lengthShort) / 2.) * sin(skewanglerad),
        switchSides = 2 * newradius * (possibleStrawsPipe + 1 + (2 * i)) * cos(skewanglerad); 
      
      putStrawRotatedShortRight(tmpxpos - (translationToLeft - additionalShift - translationToFrontPer) * sin(sixtyrad), 
                                tmpypos - (translationToLeft - additionalShift - translationToFrontPer) * cos(sixtyrad), 
                                zpos - translationToFrontPar, 
                                vectorx, vectory, vectorz, lengthShort);
 
      putStrawRotatedShortRight(tmpxpos - (translationToLeft - additionalShift + translationToFrontPer + switchSides) * sin(sixtyrad), 
                                tmpypos - (translationToLeft - additionalShift + translationToFrontPer + switchSides) * cos(sixtyrad), 
                                zpos + translationToFrontPar, 
                                vectorx, vectory, vectorz, lengthShort);
      
      tmpxpos += 2 * newradius * sin((60. / 180.) * pi);
      tmpypos += 2 * newradius * cos((60. / 180.) * pi);
    }
  
  // plane 5
  vectorx = 1.;
  vectory = 0.;
  vectorz = 0.;
  xpos += availableSpace * sin((60. / 180.) * pi);
  ypos += availableSpace * cos((60. / 180.) * pi);
  tmpxpos = xpos;
  tmpypos = ypos + extraspace + newradius;

  for(int i = 0; i < possibleStraws; i++) 
    {
      putStrawRotatedRight(tmpxpos, tmpypos - (translationToLeft - additionalShift), zpos, vectorx, vectory, vectorz, tubeLength);
      tmpypos += 2 * newradius;
    }
  
  for (int i = 0; i < extraStraws; i++)
    {
      double
        lengthShort = lengthsShort[i];

      double
        translationToFrontPar = ((tubeLength - lengthShort) / 2.) * cos(skewanglerad),
        translationToFrontPer = ((tubeLength - lengthShort) / 2.) * sin(skewanglerad),
        switchSides = 2 * newradius * (possibleStraws + 1 + (2 * i)) * cos(skewanglerad); 
      
      putStrawRotatedShortRight(tmpxpos, tmpypos - (translationToLeft - additionalShift - translationToFrontPer), 
                                zpos - translationToFrontPar, 
                                vectorx, vectory, vectorz, lengthShort);
 
      putStrawRotatedShortRight(tmpxpos, tmpypos - (translationToLeft - additionalShift + translationToFrontPer + switchSides), 
                                zpos + translationToFrontPar, 
                                vectorx, vectory, vectorz, lengthShort);

      tmpypos += 2 * newradius;
    }
  
  // plane 6
  vectorx = cos((60. / 180.) * pi);
  vectory = sin((60. / 180.) * pi);
  vectorz = 0;
  ypos += availableSpace;
  tmpxpos = xpos - (extraspacePipe + newradius) * sin((60. / 180.) * pi);
  tmpypos = ypos + (extraspacePipe + newradius) * cos((60. / 180.) * pi);

  for(int i = 0; i < possibleStrawsPipe; i++) 
    {
      
      putStrawRotatedRight(tmpxpos + (translationToLeft - additionalShift) * sin(sixtyrad), 
                           tmpypos - (translationToLeft - additionalShift) * cos(sixtyrad), 
                           zpos, vectorx, vectory, vectorz, tubeLength);
      
      tmpxpos -= 2 * newradius * sin((60. / 180.) * pi);
      tmpypos += 2 * newradius * cos((60. / 180.) * pi);
    }
  
  for (int i = 0; i < extraStrawsPipe; i++)
    {
      double
        lengthShort = lengthsShortPipe[i];

      double
        translationToFrontPar = ((tubeLength - lengthShort) / 2.) * cos(skewanglerad),
        translationToFrontPer = ((tubeLength - lengthShort) / 2.) * sin(skewanglerad),
        switchSides = 2 * newradius * (possibleStrawsPipe + 1 + (2 * i)) * cos(skewanglerad); 

      putStrawRotatedShortRight(tmpxpos + (translationToLeft - additionalShift - translationToFrontPer) * sin(sixtyrad), 
                                tmpypos - (translationToLeft - additionalShift - translationToFrontPer) * cos(sixtyrad), 
                                zpos - translationToFrontPar, 
                                vectorx, vectory, vectorz, lengthShort);
      
      putStrawRotatedShortRight(tmpxpos + (translationToLeft - additionalShift + translationToFrontPer + switchSides) * sin(sixtyrad), 
                                tmpypos - (translationToLeft - additionalShift + translationToFrontPer + switchSides) * cos(sixtyrad), 
                                zpos + translationToFrontPar, 
                                vectorx, vectory, vectorz, lengthShort);

      tmpxpos -= 2 * newradius * sin((60. / 180.) * pi);
      tmpypos += 2 * newradius * cos((60. / 180.) * pi);
    }
  
}

void placeSingleLayerStraightExact

(

double

ringteller

)

Definition at line 1108 of file createSTTwithPassive.C.

References i, putStraw(), sqrt(), tubeSeperation, and zpos.

{  
  double
    radius = tubeSeperation / 2.;

  // place a straight double layer  
  double
    xpos = 0.,  
    ypos = ((ringteller) * 2 * radius), 
    zpos = 0.;

  // 6 loops, one for each side of the hexagon
  //      1  / \  4
  //      2  | |  5
  //      3  \ /  6
  for(int i = 0; i < ringteller; i++) 
    {
      xpos -= sqrt(3.) * radius;
      ypos -= radius;
      putStraw(xpos, ypos, zpos);
    }
  for(int i = 0; i < ringteller; i++) 
    {
      ypos -= 2 * radius;
      putStraw(xpos, ypos, zpos);
    }
  for(int i = 0; i < ringteller; i++) 
    {
      xpos += sqrt(3.) *  radius;
      ypos -= radius;
      putStraw(xpos, ypos, zpos);
    }
  for(int i = 0; i < ringteller; i++) 
    {
      xpos += sqrt(3.) * radius;
      ypos += radius;
      putStraw(xpos, ypos, zpos);
    }
  for(int i = 0; i < ringteller; i++) 
    {
      ypos += 2 * radius;
      putStraw(xpos, ypos, zpos);
    }
  for(int i = 0; i < ringteller; i++) 
    {
      xpos -= sqrt(3.) * radius;
      ypos += radius;
      putStraw(xpos, ypos, zpos);
    }
}

void plotrotleft	(	double	x,
		double	y,
		double	z
	)

Definition at line 426 of file createSTTwithPassive.C.

References cos(), counter6, pi, sin(), and skewangle.

{
  cout << 1 + (1 - cos(skewangle * (pi / 180.))) * (x * x - 1) << "  " 
       << -z * sin(skewangle * (pi / 180.)) + (1 - cos(skewangle * (pi / 180.))) * x * y << "  " 
       << y * sin(skewangle * (pi / 180.)) + (1 - cos(skewangle * (pi / 180.))) * x * z << "  " 

       << z * sin(skewangle * (pi / 180.)) + (1 - cos(skewangle * (pi / 180.))) * x * y << "  " 
       << 1 + (1 - cos(skewangle * (pi / 180.))) * (y * y - 1) << "  "
       << -x * sin(skewangle * (pi / 180.)) + (1 - cos(skewangle * (pi / 180.))) * y * z << "  " 

       << -y * sin(skewangle * (pi / 180.)) + (1 - cos(skewangle * (pi / 180.))) * x * z << "  " 
       << x * sin(skewangle * (pi / 180.)) + (1 - cos(skewangle * (pi / 180.))) * y * z << "  " 
       << 1 + (1 - cos(skewangle * (pi / 180.))) * (z * z - 1) << endl;
  cout << "//----------------------------------------------------------" << endl;
  counter6++;
}

void plotrotright	(	double	x,
		double	y,
		double	z
	)

Definition at line 405 of file createSTTwithPassive.C.

References cos(), counter6, pi, sin(), and skewangle.

{
  double
    newskewangle = skewangle * -1.;

  cout << 1 + (1 - cos(newskewangle * (pi / 180.))) * (x * x - 1) << "  " 
       << -z * sin(newskewangle * (pi / 180.)) + (1 - cos(newskewangle * (pi / 180.))) * x * y << "  " 
       << y * sin(newskewangle * (pi / 180.)) + (1 - cos(newskewangle * (pi / 180.))) * x * z << "  " 

       << z * sin(newskewangle * (pi / 180.)) + (1 - cos(newskewangle * (pi / 180.))) * x * y << "  " 
       << 1 + (1 - cos(newskewangle * (pi / 180.))) * (y * y - 1) << "  "
       << -x * sin(newskewangle * (pi / 180.)) + (1 - cos(newskewangle * (pi / 180.))) * y * z << "  " 

       << -y * sin(newskewangle * (pi / 180.)) + (1 - cos(newskewangle * (pi / 180.))) * x * z << "  " 
       << x * sin(newskewangle * (pi / 180.)) + (1 - cos(newskewangle * (pi / 180.))) * y * z << "  " 
       << 1 + (1 - cos(newskewangle * (pi / 180.))) * (z * z - 1) << endl;
  cout << "//----------------------------------------------------------" << endl;
  counter6++;
}

bool putStraw	(	double	posX,
		double	posY,
		double	posZ
	)

Definition at line 199 of file createSTTwithPassive.C.

References ABS, air, axialtubeteller, beltDistance, beltInnerDiam, beltLength, crimppinLength, Cu, CuBe, endcapInnerDiam, endcapLength, endcapOuterDiam, GasMixture, gassupplyInnerDiam, gassupplyLength, gassupplyOuterDiam, gassupplyTransY, GFK, innerCoverThickness, innerDiam, maximumradius, minimumradius, mylar, noseLength, noseOuterDiam, outerCoverThickness, outerDiam, panelthickness, passive, pipeDiam, PVC, ringInnerDiam, ringLength, springInnerDiam, springLength, springOuterDiam, sqrt(), sttassembly, topcapInnerDiam, topcapLength, tubeInnerDiam, tubeLength, tubeOuterDiam, tubeteller, W, wireDiam, writemedium(), writemother(), writename(), writerot(), writetrans(), and writetube().

{

  // CHECK PIPE
  if(posX > -(pipeDiam/2. + panelthickness + tubeOuterDiam/2.) && posX < (pipeDiam/2. + panelthickness  + tubeOuterDiam/2.)) return false; 

  // check if the straws fit wihin the inner and outer diameter specified
  if (
      (sqrt(posX * posX + posY * posY) < ((outerDiam / 2.) - outerCoverThickness - (tubeOuterDiam / 2.))) &&
      (sqrt(posX * posX + posY * posY) > ((innerDiam / 2.) + innerCoverThickness + (tubeOuterDiam / 2.)))
      )
    {
      // convert int to string
      stringstream
        conv;
        
      string
        tubetellerStr;
        
      conv << tubeteller + 1;
      conv >> tubetellerStr;
        
      string
        nameItube = "stt01tube#" + tubetellerStr,
        nameIwall = "stt01wall#" + tubetellerStr,
        nameIgas = "stt01gas#" + tubetellerStr,
        nameIwire = "stt01wire#" + tubetellerStr,
        nameIendcap = "stt01endcap",
        nameItopcap = "stt01topcap",
        nameIcrimppin = "stt01crimppin",
        nameIbelt = "stt01belt",
        nameIgassupply = "stt01gassupply",
        nameIring = "stt01ring",
        nameIspring = "stt01spring",
        nameInose = "stt01nose";
        
      // keep track of minimum and maximum extend of the straw package 
      if (sqrt(posX * posX + posY * posY) - (tubeOuterDiam / 2.) < minimumradius)
        {
          minimumradius = sqrt(posX * posX + posY * posY) - (tubeOuterDiam / 2.);
        }
      if (sqrt(posX * posX + posY * posY) + (tubeOuterDiam / 2.) > maximumradius)
        {
          maximumradius = sqrt(posX * posX + posY * posY) + (tubeOuterDiam / 2.);
        }

      // place the tube
      if (tubeteller == 0)
        {
          // original tube
          writename(nameItube.c_str());
          writemother(sttassembly, 1, 0, 1);
          writemedium(air);
          writetrans(posX, posY, posZ);
          writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
         
          // Mylar wall
          writename(nameIwall.c_str());
          writemother(nameItube.c_str()); 
          writemedium(mylar);
          writetube(0, tubeOuterDiam / 2., tubeLength / 2.);
          writetrans(0., 0., 0.);
          writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
            
          // Gas filling
          writename(nameIgas.c_str());
          writemother(nameIwall.c_str());
          writemedium(GasMixture);
          writetube(0., tubeInnerDiam / 2., tubeLength / 2.);
          writetrans(0., 0., 0.);
          writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
            
          // Anode wire
          writename(nameIwire.c_str());
          writemother(nameIgas.c_str());
          writemedium(W);
          writetube(0., wireDiam / 2., tubeLength / 2.);
          writetrans(0., 0., 0.);
          writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);

          if (passive) {
            // Endcaps
            writename(nameIendcap.c_str(), 1, 1);
            writemother(nameIgas.c_str());
            writemedium(ABS);
            writetube(endcapInnerDiam / 2., endcapOuterDiam / 2., endcapLength / 2.);
            writetrans(0., 0., (tubeLength / 2. - topcapLength - endcapLength / 2.)); 
            writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
            
            writename(nameIendcap.c_str(), 1, 0);
            writemother(nameIgas.c_str(), false);
            writetrans(0., 0., -(tubeLength / 2. - topcapLength - endcapLength / 2.));
            writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
            
            // Top cap
            writename(nameItopcap.c_str(), 1, 1);
            writemother(nameIgas.c_str());
            writemedium(ABS);
            writetube(topcapInnerDiam / 2., endcapOuterDiam / 2., topcapLength / 2.);
            writetrans(0., 0., (tubeLength / 2. - topcapLength / 2.));
            writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
            
            writename(nameItopcap.c_str(), 1, 0);
            writemother(nameIgas.c_str(), false);
            writetrans(0., 0., -(tubeLength / 2. - topcapLength / 2.));
            writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
            
            // Crimp pin
            writename(nameIcrimppin.c_str(), 1, 1);
            writemother(nameItube.c_str());
            writemedium(Cu);
            writetube(wireDiam / 2., topcapInnerDiam / 2., crimppinLength / 2.);
            writetrans(0., 0., (tubeLength / 2. + crimppinLength / 2.));
            writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
            
            writename(nameIcrimppin.c_str(), 1, 0);
            writemother(nameItube.c_str(), false);
            writetrans(0., 0., -(tubeLength / 2. + crimppinLength / 2.));
            writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
            
            // Belt
            writename(nameIbelt.c_str(), 1, 1);
            writemother(nameItube.c_str());
            writemedium(GFK);
            writetube(beltInnerDiam / 2., tubeInnerDiam / 2., beltLength / 2.);
            writetrans(0., 0., (tubeLength / 2. + beltDistance + beltLength / 2.));
            writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
            
            writename(nameIbelt.c_str(), 1, 0);
            writemother(nameItube.c_str(), false);
            writetrans(0., 0., -(tubeLength / 2. + beltDistance + beltLength / 2.));
            writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
            
            // Gas supply
            writename(nameIgassupply.c_str(), 1, 1);
            writemother(nameItube.c_str());
            writemedium(PVC);
            writetube(gassupplyInnerDiam / 2., gassupplyOuterDiam / 2., gassupplyLength / 2.);
            writetrans(0., gassupplyTransY, (tubeLength / 2. + beltDistance + beltLength + gassupplyLength / 2.));
            writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
            
            writename(nameIgassupply.c_str(), 1, 0);
            writemother(nameItube.c_str(), false);
            writetrans(0., gassupplyTransY, -(tubeLength / 2. + beltDistance + beltLength + gassupplyLength / 2.));
            writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
            
            // Ring
            writename(nameIring.c_str(), 1, 1);
            writemother(nameItube.c_str());
            writemedium(CuBe);
            writetube(ringInnerDiam / 2., tubeInnerDiam / 2., ringLength / 2.);
            writetrans(0., 0., (tubeLength / 2. + ringLength / 2.));
            writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
            
            writename(nameIring.c_str(), 1, 0);
            writemother(nameItube.c_str(), false);
            writetrans(0., 0., -(tubeLength / 2. + ringLength / 2.));
            writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
            
            // Spring
            writename(nameIspring.c_str(), 1, 1);
            writemother(nameItube.c_str());
            writemedium(CuBe);
            writetube(springInnerDiam / 2., springOuterDiam / 2., springLength / 2.);
            writetrans(0., 0., (tubeLength / 2. + springLength / 2.));
            writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
            
            writename(nameIspring.c_str(), 1, 0);
            writemother(nameItube.c_str(), false);
            writetrans(0., 0., -(tubeLength / 2. + springLength / 2.));
            writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);

            // Nose
            writename(nameInose.c_str(), 1, 1);
            writemother(nameItube.c_str());
            writemedium(ABS);
            writetube( topcapInnerDiam / 2., noseOuterDiam / 2., noseLength / 2.);
            writetrans(0., 0., (tubeLength/2. + noseLength/2.));
            writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);

            writename(nameInose.c_str(), 1, 0);
            writemother(nameItube.c_str(), false);
            writetrans(0., 0., -(tubeLength/2. + noseLength/2.));
            writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
          }
        }
      else 
        {
          // copy 
          writename(nameItube.c_str());
          writemother(sttassembly, false);
          writetrans(posX, posY, posZ);
          writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        }
      tubeteller++;
      axialtubeteller++;
    }
  else
    {
      return false;
    }
    
  return true;
}

bool putStrawRotatedLeft	(	double	posX,
		double	posY,
		double	posZ,
		double	xvector,
		double	yvector,
		double	zvector,
		double	length
	)

Definition at line 766 of file createSTTwithPassive.C.

References air, beltDistance, beltLength, crimppinLength, endcapLength, GasMixture, gassupplyLength, gassupplyTransY, mylar, noseLength, passive, plotrotleft(), ringLength, sin(), skewangle, skewedtubeteller, springLength, sttassembly, topcapLength, tubeInnerDiam, tubeOuterDiam, tubeteller, W, wireDiam, writemedium(), writemother(), writename(), writerot(), writetrans(), and writetube().

{
  double extr1 = posX + (length/2.) * sin(skewangle);
  double extr2 = posX - (length/2.) * sin(skewangle);
  

  // convert int to string
  stringstream
    conv;
      
  string
    tubetellerStr;
    
  conv << tubeteller + 1;
  conv >> tubetellerStr;
    
  string
    nameItube = "stt01tube#" + tubetellerStr,
    nameIwall = "stt01wall#" + tubetellerStr,
    nameIgas = "stt01gas#" + tubetellerStr,
    nameIwire = "stt01wire#" + tubetellerStr,
    nameIendcap = "stt01endcap",
    nameItopcap = "stt01topcap",
    nameIcrimppin = "stt01crimppin",
    nameIbelt = "stt01belt",
    nameIgassupply = "stt01gassupply",
    nameIring = "stt01ring",
    nameIspring = "stt01spring",
    nameInose = "stt01nose";
  

  if (tubeteller == 0)
    {
      // original straw
      writename(nameItube.c_str());
      writemother(sttassembly, 1, 0, 1);
      writemedium(air);
      writetrans(posX, posY, posZ);
      plotrotleft(xvector, yvector, zvector);
        
      // Mylar wall
      writename(nameIwall.c_str());
      writemother(nameItube.c_str()); 
      writemedium(mylar);
      writetube(0, tubeOuterDiam / 2., length / 2.);
      writetrans(0., 0., 0.);
      writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);

      // Gas filling
      writename(nameIgas.c_str());
      writemother(nameItube.c_str());
      writemedium(GasMixture);
      writetube(0., tubeInnerDiam / 2., length / 2.);
      writetrans(0., 0., 0.);
      writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
      // Anode wire
      writename(nameIwire.c_str());
      writemother(nameIgas.c_str());
      writemedium(W);
      writetube(0., wireDiam / 2., length / 2.);
      writetrans(0., 0., 0.);
      writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);

      if (passive) {
        // Endcaps
        writename(nameIendcap.c_str(), 1, 1);
        writemother(nameIgas.c_str(), false);
        writetrans(0., 0., length / 2. - topcapLength - endcapLength / 2.);
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameIendcap.c_str(), 1, 0);
        writemother(nameIgas.c_str(), false);
        writetrans(0., 0., -(length / 2. - topcapLength - endcapLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        // Top cap
        writename(nameItopcap.c_str(), 1, 1);
        writemother(nameIgas.c_str(), false);
        writetrans(0., 0., (length / 2. - topcapLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameItopcap.c_str(), 1, 0);
        writemother(nameIgas.c_str(), false);
        writetrans(0., 0., -(length / 2. - topcapLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        // Crimp pin
        writename(nameIcrimppin.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length / 2. + crimppinLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameIcrimppin.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0.,-(length / 2. + crimppinLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        // Belt
        writename(nameIbelt.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length / 2. + beltDistance + beltLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameIbelt.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., -(length / 2. + beltDistance + beltLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        // Gas supply
        writename(nameIgassupply.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., gassupplyTransY, (length / 2. + beltDistance + beltLength + gassupplyLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameIgassupply.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., gassupplyTransY, -(length / 2. + beltDistance + beltLength + gassupplyLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        // Ring
        writename(nameIring.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length / 2. + ringLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameIring.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., -(length / 2. + ringLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        // Spring
        writename(nameIspring.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length / 2. + springLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameIspring.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., -(length / 2. + springLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        // Nose
        writename(nameInose.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length/2. + noseLength/2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameInose.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., -(length/2. + noseLength/2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      }
    }
  else
    {
      // copy volume

      writename(nameItube.c_str());
      writemother(sttassembly, false);
      writetrans(posX, posY, posZ);
      plotrotleft(xvector, yvector, zvector);
    }

  tubeteller++;
  skewedtubeteller++;
    
  return true;
}

bool putStrawRotatedRight	(	double	posX,
		double	posY,
		double	posZ,
		double	xvector,
		double	yvector,
		double	zvector,
		double	length
	)

Definition at line 936 of file createSTTwithPassive.C.

References air, beltDistance, beltLength, crimppinLength, endcapLength, GasMixture, gassupplyLength, gassupplyTransY, mylar, noseLength, passive, plotrotright(), ringLength, sin(), skewangle, skewedtubeteller, springInnerDiam, springLength, springOuterDiam, sttassembly, topcapLength, tubeInnerDiam, tubeOuterDiam, tubeteller, W, wireDiam, writemedium(), writemother(), writename(), writerot(), writetrans(), and writetube().

{


  double extr1 = posX + (length/2.) * sin(skewangle);
  double extr2 = posX - (length/2.) * sin(skewangle);
  

  // convert int to string
  stringstream
    conv;
    
  string
    tubetellerStr;
    
  conv << tubeteller + 1;
  conv >> tubetellerStr;
    
  string
    nameItube = "stt01tube#" + tubetellerStr,
    nameIwall = "stt01wall#" + tubetellerStr,
    nameIgas = "stt01gas#" + tubetellerStr,
    nameIwire = "stt01wire#" + tubetellerStr,
    nameIendcap = "stt01endcap",
    nameItopcap = "stt01topcap",
    nameIcrimppin = "stt01crimppin",
    nameIbelt = "stt01belt",
    nameIgassupply = "stt01gassupply",
    nameIring = "stt01ring",
    nameIspring = "stt01spring",
    nameInose = "stt01nose";
    
  if (tubeteller == 0)
    {
      // original volume
      writename(nameItube.c_str());
      writemother(sttassembly, 1, 0, 1);
      writemedium(air);
      writetrans(posX, posY, posZ);
      plotrotright(xvector, yvector, zvector);  
      
      // Mylar wall
      writename(nameIwall.c_str());
      writemother(nameItube.c_str()); 
      writemedium(mylar);
      writetube(0, tubeOuterDiam / 2., length / 2.);
      writetrans(0., 0., 0.);
      writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);

      // Gas filling
      writename(nameIgas.c_str());
      writemother(nameItube.c_str());
      writemedium(GasMixture);
      writetube(0., tubeInnerDiam / 2., length / 2.);
      writetrans(0., 0., 0.);
      writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
      // Anode wire
      writename(nameIwire.c_str());
      writemother(nameIgas.c_str());
      writemedium(W);
      writetube(0., wireDiam / 2., length / 2.);
      writetrans(0., 0., 0.);
      writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);

      if (passive) {
        // Endcaps
        writename(nameIendcap.c_str(), 1, 1);
        writemother(nameIgas.c_str(), false);
        writetrans(0., 0., length / 2. - topcapLength - endcapLength / 2.);
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        writename(nameIendcap.c_str(), 1, 0);
        writemother(nameIgas.c_str(), false);
        writetrans(0., 0., -(length / 2. - topcapLength - endcapLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        // Top cap
        writename(nameItopcap.c_str(), 1, 1);
        writemother(nameIgas.c_str(), false);
        writetrans(0., 0., (length / 2. - topcapLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        writename(nameItopcap.c_str(), 1, 0);
        writemother(nameIgas.c_str(), false);
        writetrans(0., 0., -(length / 2. - topcapLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        // Crimp pin
        writename(nameIcrimppin.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length / 2. + crimppinLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        writename(nameIcrimppin.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0.,-(length / 2. + crimppinLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        // Belt
        writename(nameIbelt.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length / 2. + beltDistance + beltLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        writename(nameIbelt.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., -(length / 2. + beltDistance + beltLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        // Gas supply
        writename(nameIgassupply.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., gassupplyTransY, (length / 2. + beltDistance + beltLength + gassupplyLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        writename(nameIgassupply.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., gassupplyTransY, -(length / 2. + beltDistance + beltLength + gassupplyLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        // Ring
        writename(nameIring.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length / 2. + ringLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        writename(nameIring.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., -(length / 2. + ringLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        // Spring
        writename(nameIspring.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length / 2. + springLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        writename(nameIspring.c_str(), 1, 0);
        writetube(springInnerDiam / 2., springOuterDiam / 2., springLength / 2.);
        writetrans(0., 0., -(length / 2. + springLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        // Nose
        writename(nameInose.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length/2. + noseLength/2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        writename(nameInose.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., -(length/2. + noseLength/2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      }
    }
  else
    {
      // copy volume 
      writename(nameItube.c_str());
      writemother(sttassembly, false);
      writetrans(posX, posY, posZ);
      plotrotright(xvector, yvector, zvector);
    }

  tubeteller++;
  skewedtubeteller++;
   
  return true;
}

bool putStrawRotatedShortLeft	(	double	posX,
		double	posY,
		double	posZ,
		double	xvector,
		double	yvector,
		double	zvector,
		double	length
	)

Definition at line 443 of file createSTTwithPassive.C.

References air, beltDistance, beltLength, crimppinLength, endcapLength, GasMixture, gassupplyLength, gassupplyTransY, mylar, noseLength, passive, plotrotleft(), ringLength, shortskewedtubeteller, sin(), skewangle, springLength, sttassembly, topcapLength, tubeInnerDiam, tubeOuterDiam, tubeteller, W, wireDiam, writemedium(), writemother(), writename(), writerot(), writetrans(), and writetube().

{
 
  double extr1 = posX + (length/2.) * sin(skewangle);
  double extr2 = posX - (length/2.) * sin(skewangle);
  
  if (length > 0.)
    {
      // convert int to string
      stringstream
        conv;
        
      string
        tubetellerStr;
        
      conv << tubeteller + 1;
      conv >> tubetellerStr;
        
      string
        nameItube = "stt01tube" + tubetellerStr,
        nameIwall = "stt01wall" + tubetellerStr,
        nameIgas = "stt01gas" + tubetellerStr,
        nameIwire = "stt01wire" + tubetellerStr,
        nameIendcap = "stt01endcap",
        nameItopcap = "stt01topcap",
        nameIcrimppin = "stt01crimppin",
        nameIbelt = "stt01belt",
        nameIgassupply = "stt01gassupply",
        nameIring = "stt01ring",
        nameIspring = "stt01spring",
        nameInose = "stt01nose";
      
      // N.B. we don't use copies here since the short straws are all different in length
      writename(nameItube.c_str());
      writemother(sttassembly, 1, 0, 1);
      writemedium(air);
      writetrans(posX, posY, posZ);
      plotrotleft(xvector, yvector, zvector);
        
      // Mylar wall
      writename(nameIwall.c_str());
      writemother(nameItube.c_str()); 
      writemedium(mylar);
      writetube(0, tubeOuterDiam / 2., length / 2.);
      writetrans(0., 0., 0.);
      writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);

      // Gas filling
      writename(nameIgas.c_str());
      writemother(nameItube.c_str());
      writemedium(GasMixture);
      writetube(0., tubeInnerDiam / 2., length / 2.);
      writetrans(0., 0., 0.);
      writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
      // Anode wire
      writename(nameIwire.c_str());
      writemother(nameIgas.c_str());
      writemedium(W);
      writetube(0., wireDiam / 2., length / 2.);
      writetrans(0., 0., 0.);
      writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);

      if (passive) {
        // Endcaps
        writename(nameIendcap.c_str(), 1, 1);
        writemother(nameIgas.c_str(), false);
        writetrans(0., 0., length / 2. - topcapLength - endcapLength / 2.);
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        writename(nameIendcap.c_str(), 1, 0);
        writemother(nameIgas.c_str(), false);
        writetrans(0., 0., -(length / 2. - topcapLength - endcapLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        // Top cap
        writename(nameItopcap.c_str(), 1, 1);
        writemother(nameIgas.c_str(), false);
        writetrans(0., 0., (length / 2. - topcapLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        writename(nameItopcap.c_str(), 1, 0);
        writemother(nameIgas.c_str(), false);
        writetrans(0., 0., -(length / 2. - topcapLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        // Crimp pin
        writename(nameIcrimppin.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length / 2. + crimppinLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        writename(nameIcrimppin.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., -(length / 2. + crimppinLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      
        // Belt
        writename(nameIbelt.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length / 2. + beltDistance + beltLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameIbelt.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., -(length / 2. + beltDistance + beltLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        // Gas supply
        writename(nameIgassupply.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., gassupplyTransY, (length / 2. + beltDistance + beltLength + gassupplyLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameIgassupply.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., gassupplyTransY, -(length / 2. + beltDistance + beltLength + gassupplyLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        // Ring
        writename(nameIring.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length / 2. + ringLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameIring.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., -(length / 2. + ringLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        // Spring
        writename(nameIspring.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length / 2. + springLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameIspring.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., -(length / 2. + springLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        // Nose
        writename(nameInose.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length/2. + noseLength/2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameInose.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., -(length/2. + noseLength/2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      }

      tubeteller++;
      shortskewedtubeteller++;

        
      return true;
    }
}

bool putStrawRotatedShortRight	(	double	posX,
		double	posY,
		double	posZ,
		double	xvector,
		double	yvector,
		double	zvector,
		double	length
	)

Definition at line 604 of file createSTTwithPassive.C.

References air, beltDistance, beltLength, crimppinLength, endcapLength, GasMixture, gassupplyLength, gassupplyTransY, mylar, noseLength, passive, plotrotright(), ringLength, shortskewedtubeteller, sin(), skewangle, springLength, sttassembly, topcapLength, tubeInnerDiam, tubeOuterDiam, tubeteller, W, wireDiam, writemedium(), writemother(), writename(), writerot(), writetrans(), and writetube().

{
  double extr1 = posX + (length/2.) * sin(skewangle);
  double extr2 = posX - (length/2.) * sin(skewangle);
  
  if (length > 0.)
    {



      // convert int to string
      stringstream
        conv;
        
      string
        tubetellerStr;
        
      conv << tubeteller + 1;
      conv >> tubetellerStr;
        
      string
        nameItube = "stt01tube" + tubetellerStr,
        nameIwall = "stt01wall" + tubetellerStr,
        nameIgas = "stt01gas" + tubetellerStr,
        nameIwire = "stt01wire" + tubetellerStr,
        nameIendcap = "stt01endcap",
        nameItopcap = "stt01topcap",
        nameIcrimppin = "stt01crimppin",
        nameIbelt = "stt01belt",
        nameIgassupply = "stt01gassupply",
        nameIring = "stt01ring",
        nameIspring = "stt01spring",
        nameInose = "stt01nose";
        
      // N.B. don't use copies here since all short straws have different length
      writename(nameItube.c_str());
      writemother(sttassembly, 1, 0, 1);
      writemedium(air);
      writetrans(posX, posY, posZ);
      plotrotright(xvector, yvector, zvector);  
      
      // Mylar wall
      writename(nameIwall.c_str());
      writemother(nameItube.c_str()); 
      writemedium(mylar);
      writetube(0, tubeOuterDiam / 2., length / 2.);
      writetrans(0., 0., 0.);
      writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);

      // Gas filling
      writename(nameIgas.c_str());
      writemother(nameItube.c_str());
      writemedium(GasMixture);
      writetube(0., tubeInnerDiam / 2., length / 2.);
      writetrans(0., 0., 0.);
      writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
      // Anode wire
      writename(nameIwire.c_str());
      writemother(nameIgas.c_str());
      writemedium(W);
      writetube(0., wireDiam / 2., length / 2.);
      writetrans(0., 0., 0.);
      writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);

      if (passive) {
        // Endcaps
        writename(nameIendcap.c_str(), 1, 1);
        writemother(nameIgas.c_str(), false);
        writetrans(0., 0., length / 2. - topcapLength - endcapLength / 2.);
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameIendcap.c_str(), 1, 0);
        writemother(nameIgas.c_str(), false);
        writetrans(0., 0., -(length / 2. - topcapLength - endcapLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        // Top cap
        writename(nameItopcap.c_str(), 1, 1);
        writemother(nameIgas.c_str(), false);
        writetrans(0., 0., (length / 2. - topcapLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameItopcap.c_str(), 1, 0);
        writemother(nameIgas.c_str(), false);
        writetrans(0., 0., -(length / 2. - topcapLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        // Crimp pin
        writename(nameIcrimppin.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length / 2. + crimppinLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameIcrimppin.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0.,-(length / 2. + crimppinLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        // Belt
        writename(nameIbelt.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length / 2. + beltDistance + beltLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameIbelt.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., -(length / 2. + beltDistance + beltLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        // Gas supply
        writename(nameIgassupply.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., gassupplyTransY, (length / 2. + beltDistance + beltLength + gassupplyLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameIgassupply.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., gassupplyTransY, -(length / 2. + beltDistance + beltLength + gassupplyLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        // Ring
        writename(nameIring.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length / 2. + ringLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameIring.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., -(length / 2. + ringLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        // Spring
        writename(nameIspring.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length / 2. + springLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameIspring.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., -(length / 2. + springLength / 2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        // Nose
        writename(nameInose.c_str(), 1, 1);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., (length/2. + noseLength/2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
        
        writename(nameInose.c_str(), 1, 0);
        writemother(nameItube.c_str(), false);
        writetrans(0., 0., -(length/2. + noseLength/2.));
        writerot(1., 0., 0., 0., 1., 0., 0., 0., 1.);
      }
        
      tubeteller++;
      shortskewedtubeteller++;
    
      return true;
    }
}

void writehalftube	(	double	inner,
		double	outer,
		double	length
	)

Definition at line 146 of file createSTTwithPassive.C.

References counter4, pi, and pipeDiam.

{
  cout << 0. << "    " << 0. << "    " << -1. * length * 10. << endl;  
  cout << inner * 10. << "   " << outer * 10. << endl;
  cout << 0. << "    " << 0. << "    " << length * 10. << endl;  
  cout << atan((pipeDiam/2.)/outer) * (180. / pi) << "    " << 180. - atan((pipeDiam/2.)/outer) * (180. / pi) << endl;
  counter4++;
}

void writemedium

(

char *

name

)

Definition at line 132 of file createSTTwithPassive.C.

References counter3.

{
  cout << name << endl;
  counter3++;
}

void writemother	(	char const *	name,
		bool	original = true,
		bool	support = false,
		bool	firsttube = false
	)

Definition at line 120 of file createSTTwithPassive.C.

References counter2.

{
  cout << name << endl;
  if ((original && !support && firsttube) || (original && support && firsttube))
    cout << "ASSEMBLY" << endl;
  if (original && !support && !firsttube)
    cout << "TUBE" << endl;
  if(original && support && !firsttube)
    cout << "TUBS" << endl;
  counter2++;
}

void writename	(	char const *	name,
		bool	support = false,
		bool	leftside = false
	)

Definition at line 112 of file createSTTwithPassive.C.

References counter1.

{
  if(support && leftside)  cout << name << "#1" << endl;
  else if(support && !leftside) cout << name << "#2" << endl;
  else cout << name << endl;
  counter1++;
}

void writepanel	(	char const *	name,
		bool	firstone,
		double	xthick,
		double	ythick,
		double	length,
		int	side
	)

Definition at line 155 of file createSTTwithPassive.C.

References AlBe, and sttassembly.

{
  xthick *= 10.;
  ythick *= 10.;
  length *= 10.;

  if(firstone) {
    cout << name << endl;
    cout << sttassembly << endl;
    cout << "BOX" << endl;
    cout << AlBe << endl;
    cout <<  xthick/2. << " " <<  ythick/2. << " " <<  length / 2. << endl;
    cout << -xthick/2. << " " <<  ythick/2. << " " <<  length / 2. << endl;
    cout <<  xthick/2. << " " <<  ythick/2. << " " << -length / 2. << endl;
    cout << -xthick/2. << " " <<  ythick/2. << " " << -length / 2. << endl;
    cout <<  xthick/2. << " " << -ythick/2. << " " <<  length / 2. << endl;
    cout << -xthick/2. << " " << -ythick/2. << " " <<  length / 2. << endl;
    cout <<  xthick/2. << " " << -ythick/2. << " " << -length / 2. << endl;
    cout << -xthick/2. << " " << -ythick/2. << " " << -length / 2. << endl;
  }
  else {
    cout << name << endl;
    cout << sttassembly << endl;
  }

}

void writerot	(	double	x00,
		double	x01,
		double	x02,
		double	x10,
		double	x11,
		double	x12,
		double	x20,
		double	x21,
		double	x22
	)

Definition at line 189 of file createSTTwithPassive.C.

References counter6.

{
  cout << x00 << "   "  << x01 << "   " << x02 << "   "  
       << x10 << "   "  << x11 << "   " << x12 << "   "  
       << x20 << "   "  << x21 << "   " << x22 << endl;
  cout << "//----------------------------------------------------------" << endl;
  counter6++;
}

void writetrans	(	double	x,
		double	y,
		double	z
	)

Definition at line 183 of file createSTTwithPassive.C.

References counter5.

{
  cout << x * 10. << "   " << y * 10. << "   " << z * 10. << endl;
  counter5++;
}

void writetube	(	double	inner,
		double	outer,
		double	length
	)

Definition at line 138 of file createSTTwithPassive.C.

References counter4.

{
  cout << 0. << "    " << 0. << "    " << -1. * length * 10. << endl;  
  cout << inner * 10. << "   " << outer * 10. << endl;
  cout << 0. << "    " << 0. << "    " << length * 10. << endl;  
  counter4++;
}

Variable Documentation

int axialtubeteller

static

Definition at line 105 of file createSTTwithPassive.C.

Referenced by main(), and putStraw().

int counter1

static

Definition at line 21 of file createSTTwithPassive.C.

Referenced by main(), and writename().

int counter2

static

Definition at line 22 of file createSTTwithPassive.C.

Referenced by main(), and writemother().

int counter3

static

Definition at line 23 of file createSTTwithPassive.C.

Referenced by main(), and writemedium().

int counter4

static

Definition at line 24 of file createSTTwithPassive.C.

Referenced by main(), writehalftube(), and writetube().

int counter5

static

Definition at line 25 of file createSTTwithPassive.C.

Referenced by main(), and writetrans().

int counter6

static

Definition at line 26 of file createSTTwithPassive.C.

Referenced by main(), plotrotleft(), plotrotright(), and writerot().

double maximumradius = 0

static

Definition at line 108 of file createSTTwithPassive.C.

Referenced by main(), and putStraw().

double minimumradius = 100000

static

Definition at line 109 of file createSTTwithPassive.C.

Referenced by main(), and putStraw().

bool passive

Definition at line 28 of file createSTTwithPassive.C.

Referenced by main(), putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

int shortskewedtubeteller

static

Definition at line 107 of file createSTTwithPassive.C.

Referenced by main(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().

int skewedtubeteller

static

Definition at line 106 of file createSTTwithPassive.C.

Referenced by main(), putStrawRotatedLeft(), and putStrawRotatedRight().

double tubeSeperation = tubeOuterDiam

Definition at line 36 of file createSTTwithPassive.C.

int tubeteller

static

Definition at line 104 of file createSTTwithPassive.C.

Referenced by main(), putStraw(), putStrawRotatedLeft(), putStrawRotatedRight(), putStrawRotatedShortLeft(), and putStrawRotatedShortRight().