Project

General

Profile

dthet_calc.txt

Kyle Hazelwood, 01/15/2015 11:46 AM

 
1

    
2
--------------
3
 if (spars.if_3bump==0) sts = alg_3bump1();
4
 else if (spars.if_3bump==1) sts = alg_svd();
5
 else sts = alg_lssol();
6
--------------
7
alg_lssol calls:
8
 sts = sorbit_ls_lssol(nced,nbpm,betac, betab,psic,psib,
9
           tune[pl],cor->stepcut,dx1,lb,ub,
10
           dthet,dx2,&status);
11
 input dx1 in mm
12
 output dthet in mRad(milli-radians)
13
--------------
14
Read g-table in cards before apply corrections
15
In corr3_g1
16
    sts = c453_td(GET,cor->plane,MDAT_GTABLE,spars.gtable_num,
17
                  spars.slot+1,1,cor->g1.x.shorts,&ndis);
18
--------------
19
add dthet in uRad:
20
 for (i = 0; i < RING_SIZE; i++) {
21
    if (cor->useced) cor->g2[i] = (short)(cor->g1.x.shorts[i] + 
22
                  cor->dthet[i] * 1000. * p_scale);
23
    /* '1000' to convert to uR,p scaled to 10 Gev/c where s.f. is defined */ 
24
  }
25
--------------