[Pythia] Extraire l'abscisse d'un max

[Gwyddon]

Salut à tous,

Si parmi vous il y a des spécialistes de Pythia, ce serait cool

Voilà j'aurais besoin d'extraire d'un histogramme l'abscisse du maximum de l'histogramme. Or je ne vois pas comment faire. J'avais pensé à créer des tableaux, et à chercher le max sur ce tableau, sauf qu'a priori je ne connais pas le découpage à l'avance de mon tableau (et ça risque de trop discrétiser pour obtenir des choses intéressantes).

Si quelqu'un le sait, je lui en serais reconnaissant
-----

[Gwyddon]

Salut,

Je fais un petit up (oui je sais pas bien )

Au pire si pas de réponses je le déplace en info ..

[shamrock]

Salut Gwyddon

Je pense pas avoir la réponse à ton problème, j'utilise Pythia en amateur
Tu fais quoi exactement, tu remplis un pybook ?

[Gwyddon]

Lool, enfin une réponse

En effet je remplis un PYBOOK. Le truc c'est que je ne sais pas comment Pythia stocke ses pybook. Ce que je pourrais faire à la rigueur c'est extraire d'un fichier dans lequel je mets le contenu du pybook (via PYDUMP) en permanence le maximum (encore que je ne sais pas comment faire là tout de suite), mais ça prendrait un temps fou...

En plus les simu que je fais tourner commencent à être un peu grosse : étude du spectre de l'énergie manquante dans la désintégration de charginos produits à l'ILC en fonction du splitting de masses neutralinos-charginos, scénario AMSB

EDIT : par curiosité, tu t'en sers pour quoi ?

EDIT 2 : je devrais dormir à cette heure ci...

[A voir en vidéo sur Futura]

[shamrock]

La seule idée qui me vient comme ça, c'est effectivement de faire un PYDUMP, puis, par exemple remplir un histo sous ROOT, là tu dois avoir plus de fonctions pour manipuler ton histo.
En général ce que je fais, j'utilise mon prog fortran qui appelle PYTHIA pour générer les datas qui m'intéressent, je mets tout ça dans un fichier et j'exploite dans un prog en C++ qui utilise ROOT. Mais bon, c'est peut être pas la meilleure méthode.

Je m'intéresse pour ma thèse au gravitino dans le modèle GMSB, l'aspect cosmo, matière noire, nucléosynthèse primordiale. Et je profite d'encadrer un stage pour regarde un canal de désintégration au LHC qui produit des gravitinos :
stau -> tau + neutralino et neutralino -> photon + gravitino

[Gwyddon]

Pas mal, je vais regarder ça en tout cas

Bon courage pour ta thèse

Ce sera dans un an et demi pour moi ce truc

[Gwyddon]

Re,

Il y a un moyen de remplir un histo 2D à taille variable ? Le manuel n'est pas très clair sur ça...