[Rivet] Thoughts on a Rivet 2.2.2 release?

Christian Gutschow christian.guetschow at tu-dresden.de
Fri May 22 10:05:51 BST 2015


Hi all,

just to add my two cents...

Yes, I agree. But I also think that YODA does need a physics-agnostic
> merging utility... or at least some merging helper functions.
> Unfortunately there's no way to tell if scale() is being called for
> normalization or some other purpose, so maybe that is an impossible goal.


I suppose you won't ever be able to come up with a script that can deal
with every possible case, but if you still want to provide a script to the
user, then at least it should cover the vast majority of cases -- which I
find is people plotting normalised or unnormalised differential cross
sections. Both of these can be tested for easily without the heuristic that
is currently implemented. One would simply have to check whether the
analysis objects have a ScaledBy attribute (most likely a differential
cross section then) and whether or not their integrals are unity.

I've got a privately modified version of the yodamerge script that does
just that -- I've attached this just in case this is useful to anyone. The
flag "--add" can also be passed to this script which will force a simple
stacking of the analysis objects without messing with the normalisation
(which is ever so handy when it comes to efficiencies and all that).

Cheers,
Chris


On 11 May 2015 at 09:21, Chris Pollard <cpollard at cern.ch> wrote:

> Hi Frank, all,
>
> Chris, I just saw your reply to Leif in the other mail:
>>
>> > If I'm not mistaken, this is more or less what we've begun
>> implementing; I
>> > think we have a ~working implementation for 1D histograms. We have not
>> > expanded this to other types yet.
>>
>> Can you iterate a bit on which work has started on this re-entrant
>> histogramming? I pinged Andy about it a few months ago when I realized
>> that it doesn't even have to be that interruptive because the basic
>> book keeping of analysis objects already exists. I would be interested
>> to hear which direction you are taking now. If you want, we could also
>> have a chat about it on Skype/Vidyo/...
>>
>
> Sorry: I didn't respond very clearly to Leif's email. The mechanism for
> storing temporary histograms that deal with the NLO counter events has been
> started. So far there are no (significant) changes to finalize() to deal
> with multiple runs.
>
> In any case, a skype chat might be nice to re-kickstart this effort a bit.
> I'm not certain I understand your statement about non-interruptive book
> keeping.
>
> Chris
>
> _______________________________________________
> Rivet mailing list
> Rivet at projects.hepforge.org
> https://www.hepforge.org/lists/listinfo/rivet
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://www.hepforge.org/lists-archive/rivet/attachments/20150522/592c700b/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: yodacombine
Type: application/octet-stream
Size: 9357 bytes
Desc: not available
URL: <https://www.hepforge.org/lists-archive/rivet/attachments/20150522/592c700b/attachment.obj>


More information about the Rivet mailing list