|
[Rivet] dividing histos in RivetHannes Jung hannes.jung at cern.chThu Mar 24 05:15:27 GMT 2011
Hi Andy thanks a lot. This example did help indeed. BTW, is there any reason, why I have to book still a AIDA histo, to get the LWH histos in the output, or am I doing something weird ? Anyway, at the end I got the histos.. Cheers Hannes On 23.03.2011, at 10:50, Andy Buckley wrote: > On 22/03/11 20:54, Hannes Jung wrote: >> dear riveties >> >> I was wondering, whether there is a simple way to divide two histograms in rivet/aida ? >> Just couldn't find it in the manual > > Hi Hannes, > > To do divisions or other histogram combinations you need to go via the > HistogramFactory, and mix and match histogram paths with histogram > pointers... it's just a bit cumbersome. To avoid getting duplicate > numerator and denominator histograms in the output, I tend to use raw > (non-persistent) LWH histograms for those. You'll find an example (well, > a few examples) in the src/Analyses/MC_GENERIC.cc source file. > > This round-about way of combining histograms is an AIDA monstrosity: a > good example of how not to design a histogramming interface! Leif may > know about a neater way to do this with the LWH histograms, whose > interface is less awkward than the AIDA ones (i.e. LWH::Histogram1D > rather than AIDA::IHistogram1D). It's all going to be nicer in Rivet > version 2.0. > > Hope that helps, > Andy > > -- > Dr Andy Buckley > SUPA Advanced Research Fellow > Particle Physics Experiment Group, University of Edinburgh > > The University of Edinburgh is a charitable body, registered in > Scotland, with registration number SC005336. > *********************************************************************** Hannes Jung Email: Hannes.Jung at cern.ch mobile :+49 40 8998 93741 http://www.desy.de/~jung Tel: +49 (0) 40 8998 3741 (DESY) Tel: +41 22 76 71691 (CERN) CERN - PH 40-3 B11 CH-1211 Genève 23 Switzerland ***********************************************************************
More information about the Rivet mailing list |