//==========================================
// Garfield V8.1 Input file for:  
//
// Vertical Drift Chamber (VDC) Simulation
// for the Qweak experiment at JLab
//
// => generates the CELL definition
//    file needed for additional Garfield
//    simulations: 
//    FIELD, XT, ARRIVAL-TIME-DISTRIBUTION, 
//    TIMING, SIGNAL
//
//================================================================
// Simplified VDC layout for Garfield
//
// VDC consists of one wire plane and two HV planes.
// Garfield forces us to use this orientation of the VDC since 
// in Garfiled  you can only define a X-RANGE **only** for a XT-PLOT ... 
// (There is no YT-PLOT or Y-RANGE for XT-PLOT defined in Garfield)
//
//     |     | / 
//     |  *  |/
//     |     |            |
//     |  * /|          Y |
//     |   / |            |
//     |  *  |            ------> X
//     | /   |        used coordinate system
//     |/ *  |
//     |     |          ------> Gravity vector
//    /|  *  |
//   / |     |
//  /  |  *  |
//
//=================================================================
// 
// Author: Klaus Grimm
//         Department of Physics
//         The College of William and Mary
//         Williamsburg, VA 23185 
//
//  email: grimm@jlab.org
//
//==========================================
// Last Change: 05/02/2004
//==========================================


GLOBAL anode      = -4100         // high voltage on planes

GLOBAL wnum       = 41            // Number of wires
GLOBAL plane      = 1.3           // distance wire plane to HV plane (13mm)
GLOBAL pitch      = 0.423         // perp. distance between signal wires   (4.23mm)

GLOBAL wdiameter  = 0.0025         // 25um wire diameter
GLOBAL wpot       = 0.0           // voltage on signal wire = 0.0 V
GLOBAL wtension   = 70.0          // wire tention = 70g
GLOBAL wlength    = 50.0*1.141    // overall wire length (45deg) 

GLOBAL wxpos      =  0.0           // all signal wires with x=0
GLOBAL wystart    =  -({wnum}*0.5-1)*pitch     // position of first signal wire:  y=-8.46mm



// Delete some result files, otherwise results from different settings/simulations
// will be saved in the same file ...   


//=============================================
& CELL
//=============================================

//-----------------------------------------------------------------------
// Simplified VDC layout for Garfield
//
// Garfield forces us to use this orientation of the VDC since 
// in Garfiled  you can only define a X-RANGE **only** for a XT-PLOT ... 
//
//   |     |            |
//   |  *  |          Y |
//   |     |            |
//   |  *  |            ------> X
//   |     |        used coordinate system
//   |  *  |
//   |     |          ------> Gravity vector
//
//-------------------------------------------------------------------------
 opt cell-pr

// Location of the HV planes
 plane x={-plane} v={anode}
 plane x={+plane} v={anode}


// Location of the signal wires
 ROWS CARTESIAN 
   S {wnum} ... 
   {wdiameter} ... 
   {wxpos} ...
   {wystart}+({pitch}*I) ...
   {wpot} ...
   {wtension} ...
   {wlength} ...
   TUNGSTEN

// Gravity points perp. to HV plane = x directection
 GRAVITY  1 0 0

// Save this &CELL definition
 WRITE DATASET vdc_cell.dat


//=============================================
& OPTIMISE
//=============================================

 //Prints the current potential settings
 DISPLAY

 //forces acting on a wire and the wire displacement
 SELECT 21
 FORCES  PRINT-SAG PLOT-SAG KEEP-SAG DETAILED ...
 ELECTROSTATICS GRAVITY ...
 SCANNING-GRID 25 25 ...
 SCAN-AREA -0.15 -0.15 +0.15 +0.15




//=============================================
& QUIT
//=============================================