[Rivet] JADE_OPAL_2000_S4300807 etc.

[Frank Siegert]

Hi Christoph,
(cc'ing to the Rivet list, to test their German ;-))

> gerade beschäftige ich mich wieder mit deiner Analyse - Andy Buckley hatte
> mich nämlich nach dem Integrationspunkt für die integrierten Binraten
> gefragt. Ich habe das mit Stefan Kluth diskutiert, und Andy hat es
> entsprechend geändert.

Ok, good.

> Woher kommen eigentlich in
>         00032       FastJets jadeJets = FastJets(fs, FastJets::JADE, 0.7);
>         00033       FastJets durhamJets = FastJets(fs, FastJets::DURHAM,
> 0.7);
> Die Zahlen 0.7 ? In e+e- verwendet man ja keinen Radius. Ich hab das mal auf
> pi/2 gesetzt, da kommt aber in einem schnellen Test mit 100 events nichts
> anderes raus (vielleicht sollte man mehr prozessieren).

Indeed, in the ee cases the R parameter will be ignored.

> Ich möchte auch ähnliche Rivet-Analysen schreiben:
>
> 1) Weinzierls "inclusive" kt Algorithmus. Das klappt einfach, kt mit einem
> Schnitt auf minimale Jetenergie.

Ok.

> 2) ee-SISCone, in meinem Messprogramm:
>
>   fastjet::SISConeSphericalPlugin * plugin =new
> fastjet::SISConeSphericalPlugin(Rparam, f, npass, Emin);
>   fastjet::JetDefinition jet_def(plugin);
>
> SISConeSpherical ist im Gegensatz zu SISCone nicht im "Wrapper enum for
> selected Fastjet jet algorithms", also habe ich in meiner Rivet-Analyse
> probiert
>
>   fastjet::SISConeSphericalPlugin * plugin = new
> fastjet::SISConeSphericalPlugin(R, f, npass, Emin);
>   addProjection(fs, plugin);
>
> das gibt
>
>   "error: no matching function for call to
> 'Rivet::MyAnalysis::addProjection(const Rivet::FinalState&,
> fastjet::SISConeSphericalPlugin*&)'
>
> Die Typen verstehe ich nicht gut, was sollte ich ändern?

You have to use the constructor of the FastJets projection which takes
a plugin as argument, cf. the documentation of the FastJets class:
http://rivet.hepforge.org/code/dev/a00209.html#aea0b268f9f0032bccb9011ff9f1610a6

> 3) ee-antikt, in meinem Messprogramm:
>
>   fastjet::JetDefinition ee_antikt(fastjet::ee_genkt_algorithm, R, -1,
> fastjet::E_scheme);
>
> Das ist wieder nicht im enum und es gibt in Rivet kein
>
>   FastJets(const FinalState& fsp, fastjet::JetAlgorithm type,
> fastjet::RecombinationScheme recom, double rparameter, int exponent);
>
> mit passenden Parametern.
> Sollte ich ein plugin schreiben wie in 2) ?

Indeed, there is currently no option in Rivet to initialise the
FastJets projection with a jet definition that needs an extra
parameter (making use of
http://fastjet.fr/repo/doxygen-3.0.6/classfastjet_1_1JetDefinition.html#a66e84c3d69756360862c1e22e12d7890).
I'm not sure whether it's also possible to use ee_genkt_algorithm as a
plugin? That would probably be the easiest solution for you.
Alternatively we'll have to extend the interface of the FastJets
projection to allow for such algorithms as well.

Cheers,
Frank