anchor.diff:        x  = x + sin_mb(tiltangle) * c_length
         xp = xp - tiltangle
       end if
  +
  +    ! CRY Only x_in0 (i.e. b) have to be assigned after the change of reference frame
  +    x_in0 = x
   
       ! For selected collimator, first turn reset particle distribution to simple pencil beam
       nprim = 3
  
positive.code:   real(kind=fPrec) x_Dump,xpDump,y_Dump,ypDump,s_Dump
    real(kind=fPrec) cRot,sRot,cRRot,sRRot
    real(kind=fPrec) s_impact,xinn,xpinn,yinn,ypinn
    real(kind=fPrec) x_in0,xp_in0,s_cry
  
    ! Initilaisation
  
  
negative.code:     zlm   = -one * length
  
      if (cdb_isCrystal(icoll)) then
        call cry_doCrystal(ie,iturn,j,mat,x,xp,z,zp,s,p,x_in0,xp_in0,zlm,nhit,nabs,lhit_pos,lhit_turn,&
          part_abs_pos_local,part_abs_turn_local,impact,indiv,c_length)
  !      p = sqrt(p**2+pmap**2)
      else
  