101def common_mixing_matrix(mtype):
102 presusylog.info(f"Will set mixing to {mtype}")
103
104
105 if mtype == 'higgsino':
106
107 param_blocks['VMIX']={}
108 param_blocks['VMIX']['1 1']='0.00E+00'
109 param_blocks['VMIX']['1 2']='1.00E+00'
110 param_blocks['VMIX']['2 1']='1.00E+00'
111 param_blocks['VMIX']['2 2']='0.00E+00'
112
113 param_blocks['UMIX']={}
114 param_blocks['UMIX']['1 1']='0.00E+00'
115 param_blocks['UMIX']['1 2']='1.00E+00'
116 param_blocks['UMIX']['2 1']='1.00E+00'
117 param_blocks['UMIX']['2 2']='0.00E+00'
118
119 param_blocks['NMIX']={}
120 param_blocks['NMIX']['1 1']=' 0.00E+00'
121 param_blocks['NMIX']['1 2']=' 0.00E+00'
122 param_blocks['NMIX']['1 3']=' 7.07E-01'
123 param_blocks['NMIX']['1 4']='-7.07E-01'
124 param_blocks['NMIX']['2 1']=' 0.00E+00'
125 param_blocks['NMIX']['2 2']=' 0.00E+00'
126 param_blocks['NMIX']['2 3']='-7.07E-01'
127 param_blocks['NMIX']['2 4']='-7.07E-01'
128 param_blocks['NMIX']['3 1']=' 1.00E+00'
129 param_blocks['NMIX']['3 2']=' 0.00E+00'
130 param_blocks['NMIX']['3 3']=' 0.00E+00'
131 param_blocks['NMIX']['3 4']=' 0.00E+00'
132 param_blocks['NMIX']['4 1']=' 0.00E+00'
133 param_blocks['NMIX']['4 2']='-1.00E+00'
134 param_blocks['NMIX']['4 3']=' 0.00E+00'
135 param_blocks['NMIX']['4 4']=' 0.00E+00'
136 if masses['1000022'] * masses['1000023'] > 0:
137 presusylog.warning("Expected N1 and N2 masses to have opposite sign for a higgsino signal. Possibly set after the mixing.")
138 elif mtype == 'winobino':
139
140 param_blocks['VMIX']={}
141 param_blocks['VMIX']['1 1']='9.72557835E-01'
142 param_blocks['VMIX']['1 2']='-2.32661249E-01'
143 param_blocks['VMIX']['2 1']='2.32661249E-01'
144 param_blocks['VMIX']['2 2']='9.72557835E-01'
145
146 param_blocks['UMIX']={}
147 param_blocks['UMIX']['1 1']='9.16834859E-01'
148 param_blocks['UMIX']['1 2']='-3.99266629E-01'
149 param_blocks['UMIX']['2 1']='3.99266629E-01'
150 param_blocks['UMIX']['2 2']='9.16834859E-01'
151
152 param_blocks['NMIX']={}
153 param_blocks['NMIX']['1 1']='9.86364430E-01'
154 param_blocks['NMIX']['1 2']='-5.31103553E-02'
155 param_blocks['NMIX']['1 3']='1.46433995E-01'
156 param_blocks['NMIX']['1 4']='-5.31186117E-02'
157 param_blocks['NMIX']['2 1']='9.93505358E-02'
158 param_blocks['NMIX']['2 2']='9.44949299E-01'
159 param_blocks['NMIX']['2 3']='-2.69846720E-01'
160 param_blocks['NMIX']['2 4']='1.56150698E-01'
161 param_blocks['NMIX']['3 1']='-6.03388002E-02'
162 param_blocks['NMIX']['3 2']='8.77004854E-02'
163 param_blocks['NMIX']['3 3']='6.95877493E-01'
164 param_blocks['NMIX']['3 4']='7.10226984E-01'
165 param_blocks['NMIX']['4 1']='-1.16507132E-01'
166 param_blocks['NMIX']['4 2']='3.10739017E-01'
167 param_blocks['NMIX']['4 3']='6.49225960E-01'
168 param_blocks['NMIX']['4 4']='-6.84377823E-01'
169 elif mtype == 'stau_maxmix':
170 param_blocks['selmix'] = {}
171
172 param_blocks['selmix'][ '3 3' ] = '0.70710678118'
173 param_blocks['selmix'][ '3 6' ] = '0.70710678118'
174 param_blocks['selmix'][ '6 3' ] = '-0.70710678118'
175 param_blocks['selmix'][ '6 6' ] = '0.70710678118'
176 elif mtype == 'stau_nomix':
177 param_blocks['selmix'] = {}
178
179 param_blocks['selmix'][ '3 3' ] = '1.0'
180 param_blocks['selmix'][ '3 6' ] = '0.0'
181 param_blocks['selmix'][ '6 3' ] = '0.0'
182 param_blocks['selmix'][ '6 6' ] = '1.0'
183 elif mtype == 'stop_maxmix':
184 param_blocks['USQMIX']={}
185
186 param_blocks['USQMIX']['1 1'] = '1.0'
187 param_blocks['USQMIX']['2 2'] = '1.0'
188 param_blocks['USQMIX']['3 3'] = ' 0.70710678118'
189 param_blocks['USQMIX']['3 6'] = ' 0.70710678118'
190 param_blocks['USQMIX']['4 4'] = '1.0'
191 param_blocks['USQMIX']['5 5'] = '1.0'
192 param_blocks['USQMIX']['6 3'] = ' 0.70710678118'
193 param_blocks['USQMIX']['6 6'] = ' -0.70710678118'
194 elif mtype == 'stop_pureL':
195 param_blocks['USQMIX']={}
196
197 param_blocks['USQMIX']['1 1'] = '1.0'
198 param_blocks['USQMIX']['2 2'] = '1.0'
199 param_blocks['USQMIX']['3 3'] = '1.0'
200 param_blocks['USQMIX']['3 6'] = '0.0'
201 param_blocks['USQMIX']['4 4'] = '1.0'
202 param_blocks['USQMIX']['5 5'] = '1.0'
203 param_blocks['USQMIX']['6 3'] = '0.0'
204 param_blocks['USQMIX']['6 6'] = '1.0'
205 elif mtype == 'stop_pureR':
206 param_blocks['USQMIX']={}
207
208 param_blocks['USQMIX']['1 1'] = '1.0'
209 param_blocks['USQMIX']['2 2'] = '1.0'
210 param_blocks['USQMIX']['3 3'] = '0.0'
211 param_blocks['USQMIX']['3 6'] = '1.0'
212 param_blocks['USQMIX']['4 4'] = '1.0'
213 param_blocks['USQMIX']['5 5'] = '1.0'
214 param_blocks['USQMIX']['6 3'] = '1.0'
215 param_blocks['USQMIX']['6 6'] = '0.0'