[Rivet] dividing Profile1D histograms

[Andy Buckley]

Hi Sercan,

The code means what it says -- there is no current implementation of
profile division. It will be supported in a forthcoming release with
completely new histogramming... or we'll happily accept patches for the
current version!

To get the same effect now I think you'll need to explicitly make a 2D
DataPointSet with its x and y values and errors set to vectors made by
looping over the numerator and denominator histo bin IDs. This is
exactly what would be implemented in that AIDA placeholder.

Andy


On 08/10/12 21:14, Sercan Sen wrote:
> 
> Hello, 
> 
> No problem with dividing two Histogram1D but is it possible to divide two Profile1D by using IHistogramFactory ? 
> 
> when I try to do this, it throws the error below.
> 
> ==========
> error: no matching function for call to ‘AIDA::IHistogramFactory::divide(const std::string, AIDA::IProfile1D&, AIDA::IProfile1D&)’
> /afs/cern.ch/sw/lcg/external/MCGenerators_hepmc2.06.08/rivet/1.8.1/x86_64-slc5-gcc43-opt/include/LWH/AIHistogramFactory.h:47: note: candidates are: virtual AIDA::IDataPointSet* AIDA::IHistogramFactory::divide(const std::string&, const AIDA::IHistogram1D&, const AIDA::IHistogram1D&)
> ==========
> 
> /afs/cern.ch/sw/lcg/external/MCGenerators_hepmc2.06.08/rivet/1.8.1/x86_64-slc5-gcc43-opt/include/LWH/AIHistogramFactory.h
> 
> this has some commented lines such as,
> 
>     // virtual IDataPointSet * divide(const std::string &, const IProfile1D &,
>     //                                const IProfile1D &) = 0;
> 
> Thanks, 
> Sercan
> _______________________________________________
> Rivet mailing list
> Rivet at projects.hepforge.org
> http://www.hepforge.org/lists/listinfo/rivet
> 


-- 
Dr Andy Buckley, SUPA Advanced Research Fellow
Particle Physics Expt Group, University of Edinburgh

The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.