//we have to know that the muon has a muon stub in CMU, CMP, CMX, BMU
        cut &= (uber.muon.det[j] != 0);
  
  // Cosmic cut
  
  // From Objects.C
        float pt = 1/fabs(my_correctedCurv(data, uber, trind));
        float p, px, py, pz;
        PVector(pt, uber.muon.eveta[j], uber.muon.phi0[j], p, px, py, pz);
  
  // Track Pt > 18 GeV
        if (!cut) return cut;
        cut &= pt> 18;
  
  // Energy deposition
        if (!cut) return cut;
        cut &= uber.muon.em[j] < ( 2+  0.5*fabs((p-100)*0.0115) +0.5*((p-100)*0.0115) );
        if (!cut) return cut;
        cut &= uber.muon.had[j] < ( 6+  0.5*fabs((p-100)*0.028) +0.5*((p-100)*0.028) );
  
  // Track cut
        if (!cut) return cut;
        cut &= (3<=  naxialsl( uber, trind));
        if (!cut) return cut;
        cut &= (3<=  nstereosl( uber, trind));
  
   //vertex cut
        if (!cut) return cut;
        cut &= fabs(uber.muon.ztrack[j]) <= 60;
  
  
  // NO fiducial cuts !