[Rivet] D0 2001 W/Z analysis

Lars Sonnenschein sonne at mail.cern.ch
Mon 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