Final calibration of Tower, SMD, pre, post
for 2005 run
using MIPS with 00xx00 signature in SMD

1- just SMD


  1. used ~360K+400K CuCu200 minB events from runs: 6049091 6049092 6049126 6049129 6049130 6049131 + 6050016 6050017 6050018 6050019 6050020
    Goal: produce relative SMD gains
  2. Iteration 0 , INPUT=raw spectra, some SMD strips masked
  3. QA of SMD slopes , Weihong
  4. OUTPUT=gain= task- > doSlopesOnly(760./1.4); // for CuCu200-minB
    Goal: produce absolute SMD gains
  5. Iteration 1 , INPUT=smd slopes, masks; no pre/post/tw gains 
      // MIP cut .........
  float thrMipSmdE=0.5/1000.;
  int emptyStripCount=7;
  float offCenter=0.7; // fiducial area of tower =offCenter^2
  float twMipRelEneLow=0.5, twMipRelEneHigh=2.; 
  int thrMipPresAdc=12; // thres over pedestal for pre/post
    Iteration 1c , SMD thrs=15 ADC over ped, OUTPUT:
  6. Summary plot of MIP distribution, ~350K events
    Only UxV condition

    Also pre/post/tw cuts included
  7. MIP response for pair of strips, for 12x2 SMD plains:
    * Landau fit for each group of strips: PDF
    * Mean of Landau fit for each plane: PDF, red=average, blue=goal=1.3 MeV
    (sum over 12 consecutive pairs of strips, 24 groups per SMD plane)
    Iteration 1d SMD thrs=0.3 MeV, OUTPUT not posted, MIP response is at 1.2 MeV instead of 1.4 MeV seen with ADCthr>15. The peaks in energy sums for large strip ID are even more pronounce.
    Iteration 1e SMD thrs=0.5 MeV, OUTPUT:
  8. Summary plot of MIP distribution, ~670K events
    Oulu UxV condition

    Also pre/post/tw cuts included
  9. MIP response for pair of strips, for 12x2 SMD plains:
    * Landau fit for each groupe of strips: PDF, green=goal=1.3 MeV, red=Landau fit
    * Mean of Landau fit for each plane: PDF, red=average, blue=goal=1.3 MeV
    (sum over 12 consecutive pairs of strips, 24 groups per SMD plane)
    Goal: average MIP correction for SMD over many plains
    Method: MIP corrections were averaged over 8 SMD plains in each layer, outliers dropped, low degree polynomial fitted, old smd gains from slopes corrected by the ratio of smooth fit function (blue) to 1.3 MeV.
  10. Absolute SMD gains (iter2-inp) red lines are averages for each SMD plain, 1.3 MeV is the goal, ignore chi2/dof on plots
    * SMD-layer-1 (closest to IP)

    * SMD-layer-2

    * SMD-layer-3

    Used script reCalSmd.C
  11. Final absolute calibration for SMD,
    X-axis =strip ID, one plot per SMD plain.
    Y-axis shows gains in ch/GeV of deposited energy for normal incident MIP traversing a strip through the apex. Grid is at 10Kch/GeV, 20K, etc. Energy/GeV = (rawAdc-ped)/gain;
    * SMD-layer-1 (closest to IP), sectors: 01V , 02U , 04V , 05U , 07V , 08U , 10V , 11U

    * SMD-layer-2 , sectors: 02V , 03U , 05V , 06U , 08V , 09U , 11V , 12U

    * SMD-layer-3 , sectors: 03V , 04U , 06V , 07U , 09V , 10U , 12V , 01U
    Goal: Verify absolute SMD gains
    Method: reco the same events once more and calculate next iteration of corrections:
  12. Iteration 2 (output) ,
    * SMD-layer-1 (closest to IP)

    * SMD-layer-2

    * SMD-layer-3

    Note, the SMD correction from iteration 2 were NOT used.

    SMD CALIBRATION IS DONE