Hi Sheila, I had a look at your plots and have a few (maybe a lot) comments. Some were already mentioned yesterday in the meeting. Plotting: --------- - when you plot number of leptons (or e, m, ..) you should fix the maximum of your plot. Either to 1 or test for the maximum of both histograms. - I like histograms better (with error bars). Maybe you can try something like h1->Draw("histo") h1->Draw("e1same") h2->Draw("histosame") h2->Draw("e1same") Depends a bit if you can still see something, but I really appreciate the histo option for the spectra - the binning for the spectra could be finer. Maybe 50 MeV? - you should say more clearly what you have plotted in the captions to your plots. Number of particle is Number of leptons per BReco event. You say it in the text but maybe you can also add it to the plot caption. - Your normalisation makes only sense, if you look at a quantity, which you count for all events. Clearly, the spectra are only for events with at least one lepton (or e, m). Here you should normalise to the number of entries. - for the spectra, do you take only the most energetic particle in an event or all? The most energetic particle should be better simulated than all, so take only that. - how do you get the gif plots? Do you know that you can use canvas.Print("plot.gif")? Number of Breco vs run number: ------------------------------ - I was thinking of something like Toyoko showed in one of the Wednesday meetings on page 4: http://www.slac.stanford.edu/BFROOT/www/Physics/WednesdayMeetings/W072005/orimoto.pdf There, the mes yield/fb-1 is plotted, but for us the number of BReco events is still fine. The binning is 200 runs, maybe 100 is also fine. For data you clearly need the luminosity, which you can get from the lumiscript with adjusted tcl files. I am not sure how this is done for MC. I guess there is always the same number of events in a run and the runnumbers are evenly distributed. But then you would expect a more or less straight line for MC. Maybe this needs some further investigation, but I am not sure if we learn something from this. R2 with R3 data comparision: ---------------------------- - the muon momentum for R2 is missing. - the number of leptons for R3 is missing. - the number of leptons for R2 looks stange, you don't have entries without a lepton (as for electrons and muons seperately). - apart from the normalisation the electron and muon spectrum look good and could be compatible. I am not sure for the lepton spectrum. - the number of leptons (e, m, l) look not as expected. I would expect a lower number of muons in R3 compared to R2 because of the decreasing efficiency. I don't expect such a difference for electrons. This can be cross checked with the PID control histograms from the PID group. - the number of Kaons looks good. R2 data/MC comparison: ---------------------- - the muon spectrum for MC is missing. - the number of lepton for data is missing. - for data you only should use OnPeak. You can check if the OffPeak data is important to the analysis. Interesting numbers are number of Breco events per fb-1, number of leptons per fb-1 and number of events with BReco, lepton and lepton momentum > 1 GeV. This number should be very low. Already the number of BRecos should be quite low in the OffPeak data. This is one reason why we use the BReco sample. - The same argument should hold for ccbar and uds events. If you still have a discrepancy between data and MC you should also look at these samples. - the number of ... plots don't look too promissing. For example the number of electrons differ by absolute ~5% in the second bin. Maybe this is coming from the OffPeak. Let's wait for the updated plots. Numbers: -------- A few numbers would be helpful: - number of events with at least one lepton/electron/muon per fb-1 (N/intLum), which is the average of the ... vs run number plots. Same for BReco. Cheers, Wolfgang