[Rivet] New Analysis on Rivet 1.2.1

Holger Schulz holger.schulz at physik.hu-berlin.de
Tue May 18 14:31:59 BST 2010


On 18/05/10 14:05, Flavia Dias wrote:
> Hello all,
>
> I'm trying to make my own analysis on Rivet 1.2.1, installed in my 
> home directory.
>
> First, my problem is, I am able to compile it with no errors using the 
> command:
>
>  g++ -fPIC -o libRivetMuonAnalysis.so -shared 
> -I/home/fladias/RIVET/local/include 
> -I/afs/cern.ch/sw/lcg/external/HepMC/2.05.00/slc4_amd64_gcc34/include 
> <http://cern.ch/sw/lcg/external/HepMC/2.05.00/slc4_amd64_gcc34/include> -I/afs/cern.ch/sw/lcg/external/Boost/1.34.1/slc4_amd64_gcc34/include/boost-1_34_1 
> <http://cern.ch/sw/lcg/external/Boost/1.34.1/slc4_amd64_gcc34/include/boost-1_34_1> 
> MuonAnalysis.cc
>
>
> But then, if I run rivet --list-analysis, my analysis won't be listed.
You need to copy the RivetMuonAnalysis.so to some place where Rivet can 
find it.
You can use 'rivet-config --libdir' to see where Rivet looks for 
libraries by
default.

>
> The code I use is attached.
>
>
>
> The other thing is: I'll have some data from D0 to create a AIDA file 
> and use the automatic histogram booking to use in this analysis. The 
> data will be public released around next week. How should I proceed to 
> use the data as a reference file for my analysis?

The file has to have the same name as the analysis, so in your case 
maybe D0_NEW_DATA.aida

You can edit a flat file for convenience and convert it to aida later 
with flat2aida.

The flatfile should look like this:

# BEGIN HISTOGRAM /REF/D0_NEW_DATA/d01-x01-y01
AidaPath=/REF/D0_NEW_DATA/d01-x01-y01
## xlow         xhigh           yval            yerrminus       yerrplus
1.000000e+00    1.250000e+00    2.644330e-01    1.448710e-02    1.448710e-02
1.250000e+00    1.500000e+00    3.543480e-01    1.907530e-02    1.907530e-02
1.500000e+00    1.750000e+00    4.518360e-01    2.404250e-02    2.404250e-02
# END HISTOGRAM

the next histogram could for instance be /REF/D0_NEW_DATA/d02-x01-y01
and so on

The histograms should be booked if you use
h1= bookHistogram1D(1, 1, 1);
h2= bookHistogram1D(2, 1, 1)

Hope that helps,
Holger


>
> Thank you very much,
>
> Best Regards,
>
> Flavia
>
>
>
>
>
> _______________________________________________
> Rivet mailing list
> Rivet at projects.hepforge.org
> http://www.hepforge.org/lists/listinfo/rivet

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://www.hepforge.org/lists-archive/rivet/attachments/20100518/676d48ee/attachment.htm 


More information about the Rivet mailing list