code: def _get_fault_type_dummy_variables(self, sites, rup, imt):
          
  
          Fn, Fr = 0, 0
          if rup.rake >= -135 and rup.rake <= -45:
              # normal
              Fn = 1
          elif rup.rake >= 45 and rup.rake <= 135:
              # reverse
              Fr = 1
          return Fn, Fr
docs: Same classification of SadighEtAl1997. Akkar and Bommer 2010 is based
          on Akkar and Bommer 2007b; read Strong-Motion Dataset and Record
          Processing on p. 514 (Akkar and Bommer 2007b).