Mercurial > public > atmospheric_lidar_ghk / file comparison
comparison: system_settings/optic_input_example_lidar.py
system_settings/optic_input_example_lidar.py
- changeset 13
- f08818615e3a
- parent 12
- 8badc005e347
- child 16
- 313ac320b970
equal
deleted
inserted
replaced
| 12:8badc005e347 | 13:f08818615e3a |
|---|---|
| 10 LID = "example lidar" # Additional lidar ID (short descriptive text) | 10 LID = "example lidar" # Additional lidar ID (short descriptive text) |
| 11 # firet fit intern (FITLN1) => DO = 0, DT fixed -0.9998, eta and DR fitted, | 11 # firet fit intern (FITLN1) => DO = 0, DT fixed -0.9998, eta and DR fitted, |
| 12 # => internal calib with LinPol before the receiver | 12 # => internal calib with LinPol before the receiver |
| 13 print(" Lidar system :", EID, ", ", LID) | 13 print(" Lidar system :", EID, ", ", LID) |
| 14 | 14 |
| 15 # +++ IL Laser and +-Uncertainty | |
| 16 bL = 1. #degree of linear polarization; default 1 | |
| 17 RotL, dRotL, nRotL = 0., 1., 1 #alpha; rotation of laser polarization in degrees; default 0 | |
| 15 | 18 |
| 16 # --- IL Laser IL and +-Uncertainty | 19 # +++ ME Emitter optics and +-Uncertainty; default = no emitter optics |
| 17 bL = 1. #degree of linear polarization; default 1 | 20 DiE, dDiE, nDiE = 0.0, 0.1, 0 # Diattenuation |
| 18 RotL, dRotL, nRotL = 90, 1., 0 #alpha; rotation of laser polarization in degrees; default 0 | 21 TiE = 1.0 # Unpolarized transmittance |
| 19 # --- ME Emitter and +-Uncertainty | 22 RetE, dRetE, nRetE = 0., 180., 0 # Retardance in degrees |
| 20 DiE, dDiE, nDiE = 0., 0.1, 0 # Diattenuation | |
| 21 TiE = 1. # Unpolarized transmittance | |
| 22 RetE, dRetE, nRetE = 0., 180.0, 0 # Retardance in degrees | |
| 23 RotE, dRotE, nRotE = 0., 1.0, 0 # beta: Rotation of optical element in degrees | 23 RotE, dRotE, nRotE = 0., 1.0, 0 # beta: Rotation of optical element in degrees |
| 24 | 24 |
| 25 # --- MO Receiver Optics including telescope | 25 # +++ MO Receiver optics including telescope |
| 26 DiO, dDiO, nDiO = 0.0, 0.0022, 0 | 26 DiO, dDiO, nDiO = 0.0, 0.01, 1 |
| 27 TiO = 1.0 | 27 TiO = 1.0 |
| 28 RetO, dRetO, nRetO = 0., 180.0, 0 | 28 RetO, dRetO, nRetO = 0., 180., 2 |
| 29 RotO, dRotO, nRotO = 0., 0.5, 0 #gamma | 29 RotO, dRotO, nRotO = 0., 0.5, 0 #gamma: Rotation of optical element in degrees |
| 30 | 30 |
| 31 # --- PBS MT transmitting path defined with TS, TP, PolFilter extinction ratio ERaT, and +-Uncertainty | 31 # +++++ PBS MT Transmitting path defined with TS, TP, PolFilter extinction ratio ERaT, and +-Uncertainty |
| 32 # --- Polarizing beam splitter | 32 # --- Polarizing beam splitter transmitting path |
| 33 TP, dTP, nTP = 0.512175, 0.0024, 1 | 33 TP, dTP, nTP = 0.95, 0.01, 1 |
| 34 TS, dTS, nTS = 1-TP, 0.02, 0 | 34 TS, dTS, nTS = 0.02, 0.01, 1 |
| 35 RetT, dRetT, nRetT = 0.0, 180., 0 # Retardance in degrees | |
| 36 # --- Pol.Filter behind transmitted path of PBS | |
| 37 ERaT, dERaT, nERaT = 0.0001, 0.0001, 1 # Extinction ratio | |
| 38 RotaT, dRotaT, nRotaT = 0., 1., 1 # Rotation of the Pol.-filter in degrees; usually 0° because TP >> TS, but for PollyXTs it can also be 90° | |
| 39 # -- | |
| 35 TiT = 0.5 * (TP + TS) | 40 TiT = 0.5 * (TP + TS) |
| 36 DiT = (TP-TS)/(TP+TS) | 41 DiT = (TP-TS)/(TP+TS) |
| 37 RetT, dRetT, nRetT = 0., 180., 0 # Retardance in degrees | |
| 38 # --- Pol.Filter | |
| 39 ERaT, dERaT, nERaT = 0.0001, 0.0001, 1 # Extinction ratio | |
| 40 RotaT, dRotaT, nRotaT = 90., 2., 0 # Rotation of the pol.-filter in degrees | |
| 41 DaT = (1-ERaT)/(1+ERaT) | 42 DaT = (1-ERaT)/(1+ERaT) |
| 42 TaT = 0.5*(1+ERaT) | 43 TaT = 0.5*(1+ERaT) |
| 43 | 44 |
| 44 # --- PBS MR reflecting path defined with RS, RP, PolFilter extinction ratio ERaR and +-Uncertainty | 45 # +++++ PBS MR Reflecting path defined with RS, RP, PolFilter extinction ratio ERaR and +-Uncertainty |
| 45 # --- Polarizing beam splitter | 46 # ---- for PBS without absorption the change of RS and RP must depend on the change of TP and TS. Hence the values and uncertainties are not independent. |
| 46 RP, dRP, nRP = 1-TP, 0.02, 0 | 47 RS_RP_depend_on_TS_TP = True |
| 47 RS, dRS, nRS = 1-TS, 0.00, 0 | 48 # --- Polarizing beam splitter reflecting path |
| 48 RetR, dRetR, nRetR = 0., 180., 0 | 49 if(RS_RP_depend_on_TS_TP): |
| 50 RP, dRP, nRP = 1-TP, 0.00, 0 # do not change this | |
| 51 RS, dRS, nRS = 1-TS, 0.00, 0 # do not change this | |
| 52 else: | |
| 53 RP, dRP, nRP = 0.05, 0.01, 1 # change this if RS_RP_depend_on_TS_TP = False | |
| 54 RS, dRS, nRS = 0.98, 0.01, 1 # change this if RS_RP_depend_on_TS_TP = False | |
| 55 RetR, dRetR, nRetR = 0.0, 180., 0 | |
| 56 # --- Pol.Filter behind reflected path of PBS | |
| 57 ERaR, dERaR, nERaR = 0.0001, 0.0001, 1 # Extinction ratio | |
| 58 RotaR, dRotaR, nRotaR = 90., 1., 1 # Rotation of the Pol.-filter in degrees; usually 90° because RS >> RP, but for PollyXTs it can also be 0° | |
| 59 # -- | |
| 49 TiR = 0.5 * (RP + RS) | 60 TiR = 0.5 * (RP + RS) |
| 50 DiR = (RP-RS)/(RP+RS) | 61 DiR = (RP-RS)/(RP+RS) |
| 51 # --- Pol.Filter | |
| 52 ERaR, dERaR, nERaR = 1, 0.003, 0 | |
| 53 RotaR, dRotaR, nRotaR = 0., 2., 0 | |
| 54 DaR = (1-ERaR)/(1+ERaR) | 62 DaR = (1-ERaR)/(1+ERaR) |
| 55 TaR = 0.5*(1+ERaR) | 63 TaR = 0.5*(1+ERaR) |
| 56 | 64 |
| 57 # --- Parallel signal detected in the transmitted channel => Y = 1, or in the reflected channel => Y = -1 | 65 # +++ Parallel signal detected in the transmitted channel => Y = +1, or in the reflected channel => Y = -1 |
| 58 Y = -1. | 66 Y = +1. |
| 59 | 67 |
| 60 # --- Calibrator Location | 68 # +++ Calibrator Location |
| 61 LocC = 3 #location of calibrator: 1 = behind laser; 2 = behind emitter; 3 = before receiver; 4 = before PBS | 69 LocC = 3 #location of calibrator: 1 = behind laser; 2 = behind emitter; 3 = before receiver; 4 = before PBS |
| 62 # --- Calibrator Type used; defined by matrix values below | 70 # --- Calibrator Type used; defined by matrix values below |
| 63 TypeC = 3 #Type of calibrator: 1 = mechanical rotator; 2 = hwp rotator (fixed retardation); 3 = linear polarizer; 4 = qwp; 5 = circular polarizer; 6 = real HWP calibration +-22.5° | 71 TypeC = 3 #Type of calibrator: 1 = mechanical rotator; 2 = hwp rotator (fixed retardation); 3 = linear polarizer; 4 = qwp; 5 = circular polarizer; 6 = real HWP calibration +-22.5° |
| 64 # --- MC Calibrator | 72 # --- MC Calibrator parameters |
| 65 if TypeC == 1: #mechanical rotator | 73 if TypeC == 1: #mechanical rotator |
| 66 DiC, dDiC, nDiC = 0., 0., 0 | 74 DiC, dDiC, nDiC = 0., 0., 0 |
| 67 TiC = 1. | 75 TiC = 1. |
| 68 RetC, dRetC, nRetC = 0., 0., 0 | 76 RetC, dRetC, nRetC = 0., 0., 0 |
| 69 RotC, dRotC, nRotC = 0., 0.1, 1 #constant calibrator offset epsilon | 77 RotC, dRotC, nRotC = 0., 0.1, 1 #constant calibrator offset epsilon |
| 98 else: | 106 else: |
| 99 print ('calibrator not implemented yet') | 107 print ('calibrator not implemented yet') |
| 100 sys.exit() | 108 sys.exit() |
| 101 | 109 |
| 102 # --- LDRCal assumed atmospheric linear depolarization ratio during the calibration measurements (first guess) | 110 # --- LDRCal assumed atmospheric linear depolarization ratio during the calibration measurements (first guess) |
| 103 LDRCal,dLDRCal,nLDRCal= 0.006, 0.02, 1 | 111 LDRCal,dLDRCal,nLDRCal= 0.008, 0.003, 0 |
| 104 | 112 |
| 105 # ==================================================== | 113 # ==================================================== |
| 106 # NOTE: there is no need to change anything below. | 114 # NOTE: there is no need to change anything below. |
| 107 | 115 |
| 108 # --- LDRtrue for simulation of measurement => LDRsim | 116 # --- LDRtrue for simulation of measurement => LDRsim |
