Analysed simon's runs 149663 - 150958? using the jbot0g? sample that was
reprocessed with 4.8.2. Files from nglas08:/cdf/scratch/waschke/temp482a/
Ran on 4.8.2 with JPsiModules checked out from the head
Used alignment table 100030 1 GOOD and no extra retracking was required
Require SVX beam position algorithm for all events
Apply the following cuts:
if( Beamtype != 3 ) return kFALSE;
if( fitprb1 < 0.001 ) return kFALSE;
if( abs(ptmu1)<1.5 ) return kFALSE;
if( abs(ptmu2)<1.5 ) return kFALSE;
if( !((mu1type==1)||(mu1type==3)) ) return kFALSE;
if( !((mu2type==1)||(mu2type==3)) ) return kFALSE;
if( cmu1x2xp>9. ) return kFALSE;
if( cmu2x2xp>9. ) return kFALSE;
if( ! ( ncthtax1>=36&&ncthtax2>=36
&&ncthtst1>=24&&ncthtst2>=24
&&nsihit1>=3&&nsihit2>=3)
) return kFALSE;
if( silayer1[1]<1 ) return kFALSE;
if( silayer2[1]<1 ) return kFALSE;
The root files were:
- jchain.C is used to chain
Simon's output files from JPsi and runs jpsi482sel
- jpsiHist.C and
jpsiHist.h define the
histograms for the running over the chain. These are used by jpsi482sel
- jpsi482sel.C and
jpsi482sel.h:
Main program that reads the jpsi ntuple, applies cuts, and creates
an output ntuple with ct, cterr, ptjpsi and mass
- masspt.C: Reads the output of
jpsi482sel and looks at mass vs pt.
To do the analysis of the JpsiNtuple (simon's output for example) start
root and :
- .L jpsiHistC+
- .L jpsi482sel.C+
- .x jchain.C
You get a file called j482.root.
To analyze that file, you can run masspt.C. Type root then:
- TFile * f = new TFile("j482.root")
- .L masspt.C+
- .x masspt.C(Fitter);
Runs
Plot 1
All Runs for this data sample(right) and Runs with Beam Info (Left)
The following files were excluded as they had no candidates:
87-116, 119-122, 129-131, 139, 160-162, 172,
178-179, 214, 272-315, 318-340.
The left and right plots are not the same because there may be runs
within a file that don't have a valid beamline.
Quantities upon which Cuts are Placed
The following plots are those quantities upon which cuts are placed.
All plots are made before any cuts are made.
Plot 2a
Uncorrected D0 (no cut used)
Plot 2b
Corrected D0 (no cut used)
Plot 2c
Raw Track Z0 ( no cut used)
Plot 2d
COT Track Axial and Stereo Hits (ax>36, st>24, si>3 for both tracks)
Plot 2e
Silicon Hits (>3 for both tracks)
Plot 2f
Lxy Errors (no cut made)
Plot 2g
Muon Chi2 (< 9)
Plot 2h
Muon Pt and Type (Pt> 1.5 for both and both are Type 1 or 3)
Plot 2j
Di-Muon Pt (no cut made: will study this during fits)
Plot 2k
Fit Probability and Silicon Layer 0 (< 0.001 and layer0>1)
Rough Results
Plot 3a
Mass and D0/SigD0. M1 is before any selection. JpsiMass is after the
cuts listed at the top of the page.
Plot 3b
Mass and Lifetime. The mass (h3) is the same as JpsiMass but with 70
instead of 100 bins.
Fit Studies
Study of the window cuts to use was done. Plots include pt dependence
of the mass. Binning in Pt of the dimuon pair was
done to match
Simon's plots. The results are given in plot and
tabular form.
Plot 4a
(surface)
Plot 4b (box)
Plot of mass vs pt in surface and box form.
low pt.
Fits of the mass using guassian plus a flat background are given for
dimu pt bins of 2 GeV:
Plot 5a
>1.5,1.5-2,2-3,3-4 GeV
Plot 5b
4-5,5-7,7-10,10-15 GeV
Plot 5c
Plots of the Mass, sigma and background vs pt. The point at 20GeV is
actually anything above 1.5GeV.
The definition of the signal window was +- 50Mev in the Run 1
analysis. The Sidebands were at 3.0 and 3.2 GeV.
The value for three standard deviations and 7 standard deviations
are given in the table below:
Pt |
Signal |
SideBands |
>1.5 | 3.053 to 3.141 | 2.994 and 3.200 |
1.5-2 | 3.056 to 3.138 | 3.000 and 3.193 |
2-3 | 3.056 to 3.138 | 3.002 and 3.192 |
3-4 | 3.054 to 3.140 | 2.997 and 3.196 |
4-5 | 3.052 to 3.142 | 2.993 and 3.201 |
5-7 | 3.050 to 3.144 | 2.987 and 3.207 |
7-10 | 3.045 to 3.149 | 2.976 and 3.218 |
10-15 | 3.037 to 3.157 | 2.957 and 3.237 |