|
[Rivet] D0 2001 W/Z analysisLars Sonnenschein sonne at mail.cern.chMon Oct 13 12:57:04 BST 2008
Hello all I have put high statistics plots with HerwigJimmy & ATLAS + W/Z Tune as well as with Pythia 6413 with Tune A at http://www-d0.fnal.gov/~sonne/D0_2001_S4674421/ D0_2001_WptZpt_HerwigJimmy6510_ATLAS_and_WZtune_datnormed_RootRatio_partons.eps and D0_2001_WptZpt_Pythia6413_TuneA_datnormed_RootRatio_partons.eps The simujation is normed to data in the first three plots and the ratio I have still computed with ROOT. Somehow I fail to manage to get the ratio written out correctly. If I put the ratio into a datapointset it is written out but all data points and errors are zero. If I write out a profile histogram, which I fill after the division there are problems due to the errors (their should be at least two entries per bin to get a mean and a variance.) Modifying the bins directly is forbidden (l-value). Adding a special function to the LWH histo class does not work either. It seems not be used but the AIDA class directly, therfore already the normalization and integral functions of the LWH class cannot be used. Either I simply remove the ratio histo from the analysis and everybody can do the divison on its own after running Rivet, or I would need to touch AIDA or I would need to touch the division. Lars On Fri, 10 Oct 2008, Lars Sonnenschein wrote: > Hello Andy et al > at > /afs/cern.ch/user/l/lsonnens/public/Pythia6test.pdf > I put a plot, where I make use of the scaling procedure. > (extracting the integrals and scaling afterwards the ratio) > Apart of the fact, that I scale to a wrong overal value via: > > _h_dsigdpt_wz_rat->scale(zpt_scaled_integral / wpt_integral * _brzee / > _brwenu); > > the shape stays the same and you can clearly see the three goups of > entries, first for 2GeV bins, then in the middle for 5GeV and rightmost > for 10GeV bins. It seems to indicate that somehow the binwidth has to be > multiplied to get the distribution about flat. > > 2. All bins are downshifted by one bin (the first group of bins should > consist of 10 measurements but there are only 9) > > So far > > Lars > -- ________________________________ Lars Sonnenschein ________________________________ Home Institution: PH/TH 53/1-052 CERN CH-1211 Geneve 23 Switzerland Tel.:+41(22)767-2801 -------------------------------- ________________________________ FNAL: D0, PK151 Mailstop #352 Fermilab, P.O.Box 500 Batavia, IL 60510-500 USA Tel.: +1(630)840-8740 ________________________________
More information about the Rivet mailing list |