Muon tools

For the fiducial efficiencies, I have not been very careful about the scale factors. Normally we need to apply a few of them:
- the trigger scale factor
- the basic muon reconstruction scale factor (which I think is the one you're asking about)
- the additional scale factors that we derive

But all of these are quite small effects, so I have been neglecting them. But Louise applies them for the final background predictions and they should also be applied for the final fiducial efficiencies. We will send you the histogram with our own derived scale factors.

For the scale factor associated with the reconstruction, we just use sf.scaleFactor(p). For the luminosity vector, we use pretty much what is on the twiki:

  std::vector lumi_MuonEff(11);
  lumi_MuonEff[0] = 11.991; // luminosity for period B
  lumi_MuonEff[1] = 166.467; // luminosity for period D
  lumi_MuonEff[2] = 50.418; // luminosity for period E
  lumi_MuonEff[3] = 136.759; // luminosity for period F
  lumi_MuonEff[4] = 517.988; // luminosity for period G
  lumi_MuonEff[5] = 264.752; // luminosity for period H
  lumi_MuonEff[6] = 333.797; // luminosity for period I
  lumi_MuonEff[7] = 233.152; // luminosity for period J
  lumi_MuonEff[8] = 551.354; // luminosity for period K
  lumi_MuonEff[9] = 1415.898; // luminosity for period L
  lumi_MuonEff[10] = 1031.540; // luminosity for period M

these are the luminosities associated with the AllGood GRL, and should be very close to the luminosities for our GRL.

Hi Nick - here are the scale factors we calculated for the additional selection cuts we are imposing in this analysis:

///////////////////////////////////////
  // this function returns the scale factor for isolation + impact parameter
  // cuts to be used for MC.
  ///////////////////////////////////////

  float mySF=0;

  if      (pt < 22.0) mySF = 0.974444;
  else if (pt < 25.0) mySF = 0.974090;
  else if (pt < 30.0) mySF = 0.979768;
  else if (pt < 40.0) mySF = 0.986316;
  else if (pt < 50.0) mySF = 0.992907;
  else if (pt < 60.0) mySF = 0.994530;
  else if (pt < 80.0) mySF = 0.994577;
  else if (pt < 110.0) mySF = 0.994276;
  else mySF = 0.99;

  return mySF;
}