--- a/system_settings/optic_input_example_lidar.py Mon Nov 14 23:42:20 2016 +0100
+++ b/system_settings/optic_input_example_lidar.py Tue Nov 15 03:37:02 2016 +0100
@@ -12,56 +12,64 @@
# => internal calib with LinPol before the receiver
print(" Lidar system :", EID, ", ", LID)
-
-# --- IL Laser IL and +-Uncertainty
+# +++ IL Laser and +-Uncertainty
bL = 1. #degree of linear polarization; default 1
-RotL, dRotL, nRotL = 90, 1., 0 #alpha; rotation of laser polarization in degrees; default 0
-# --- ME Emitter and +-Uncertainty
-DiE, dDiE, nDiE = 0., 0.1, 0 # Diattenuation
-TiE = 1. # Unpolarized transmittance
-RetE, dRetE, nRetE = 0., 180.0, 0 # Retardance in degrees
+RotL, dRotL, nRotL = 0., 1., 1 #alpha; rotation of laser polarization in degrees; default 0
+
+# +++ ME Emitter optics and +-Uncertainty; default = no emitter optics
+DiE, dDiE, nDiE = 0.0, 0.1, 0 # Diattenuation
+TiE = 1.0 # Unpolarized transmittance
+RetE, dRetE, nRetE = 0., 180., 0 # Retardance in degrees
RotE, dRotE, nRotE = 0., 1.0, 0 # beta: Rotation of optical element in degrees
-# --- MO Receiver Optics including telescope
-DiO, dDiO, nDiO = 0.0, 0.0022, 0
+# +++ MO Receiver optics including telescope
+DiO, dDiO, nDiO = 0.0, 0.01, 1
TiO = 1.0
-RetO, dRetO, nRetO = 0., 180.0, 0
-RotO, dRotO, nRotO = 0., 0.5, 0 #gamma
-
-# --- PBS MT transmitting path defined with TS, TP, PolFilter extinction ratio ERaT, and +-Uncertainty
-# --- Polarizing beam splitter
-TP, dTP, nTP = 0.512175, 0.0024, 1
-TS, dTS, nTS = 1-TP, 0.02, 0
+RetO, dRetO, nRetO = 0., 180., 2
+RotO, dRotO, nRotO = 0., 0.5, 0 #gamma: Rotation of optical element in degrees
+
+# +++++ PBS MT Transmitting path defined with TS, TP, PolFilter extinction ratio ERaT, and +-Uncertainty
+# --- Polarizing beam splitter transmitting path
+TP, dTP, nTP = 0.95, 0.01, 1
+TS, dTS, nTS = 0.02, 0.01, 1
+RetT, dRetT, nRetT = 0.0, 180., 0 # Retardance in degrees
+# --- Pol.Filter behind transmitted path of PBS
+ERaT, dERaT, nERaT = 0.0001, 0.0001, 1 # Extinction ratio
+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°
+# --
TiT = 0.5 * (TP + TS)
DiT = (TP-TS)/(TP+TS)
-RetT, dRetT, nRetT = 0., 180., 0 # Retardance in degrees
-# --- Pol.Filter
-ERaT, dERaT, nERaT = 0.0001, 0.0001, 1 # Extinction ratio
-RotaT, dRotaT, nRotaT = 90., 2., 0 # Rotation of the pol.-filter in degrees
DaT = (1-ERaT)/(1+ERaT)
-TaT = 0.5*(1+ERaT)
+TaT = 0.5*(1+ERaT)
-# --- PBS MR reflecting path defined with RS, RP, PolFilter extinction ratio ERaR and +-Uncertainty
-# --- Polarizing beam splitter
-RP, dRP, nRP = 1-TP, 0.02, 0
-RS, dRS, nRS = 1-TS, 0.00, 0
-RetR, dRetR, nRetR = 0., 180., 0
+# +++++ PBS MR Reflecting path defined with RS, RP, PolFilter extinction ratio ERaR and +-Uncertainty
+# ---- 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.
+RS_RP_depend_on_TS_TP = True
+# --- Polarizing beam splitter reflecting path
+if(RS_RP_depend_on_TS_TP):
+ RP, dRP, nRP = 1-TP, 0.00, 0 # do not change this
+ RS, dRS, nRS = 1-TS, 0.00, 0 # do not change this
+else:
+ RP, dRP, nRP = 0.05, 0.01, 1 # change this if RS_RP_depend_on_TS_TP = False
+ RS, dRS, nRS = 0.98, 0.01, 1 # change this if RS_RP_depend_on_TS_TP = False
+RetR, dRetR, nRetR = 0.0, 180., 0
+# --- Pol.Filter behind reflected path of PBS
+ERaR, dERaR, nERaR = 0.0001, 0.0001, 1 # Extinction ratio
+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°
+# --
TiR = 0.5 * (RP + RS)
DiR = (RP-RS)/(RP+RS)
-# --- Pol.Filter
-ERaR, dERaR, nERaR = 1, 0.003, 0
-RotaR, dRotaR, nRotaR = 0., 2., 0
DaR = (1-ERaR)/(1+ERaR)
-TaR = 0.5*(1+ERaR)
+TaR = 0.5*(1+ERaR)
-# --- Parallel signal detected in the transmitted channel => Y = 1, or in the reflected channel => Y = -1
-Y = -1.
+# +++ Parallel signal detected in the transmitted channel => Y = +1, or in the reflected channel => Y = -1
+Y = +1.
-# --- Calibrator Location
+# +++ Calibrator Location
LocC = 3 #location of calibrator: 1 = behind laser; 2 = behind emitter; 3 = before receiver; 4 = before PBS
# --- Calibrator Type used; defined by matrix values below
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°
-# --- MC Calibrator
+# --- MC Calibrator parameters
if TypeC == 1: #mechanical rotator
DiC, dDiC, nDiC = 0., 0., 0
TiC = 1.
@@ -100,7 +108,7 @@
sys.exit()
# --- LDRCal assumed atmospheric linear depolarization ratio during the calibration measurements (first guess)
-LDRCal,dLDRCal,nLDRCal= 0.006, 0.02, 1
+LDRCal,dLDRCal,nLDRCal= 0.008, 0.003, 0
# ====================================================
# NOTE: there is no need to change anything below.
