[Prev][Next][Index][Thread]

istruzioni 6pioni



 Scopo e piano di massima del lavoro 6pioni focus
----------------------------------------------------
 1) riprodurre i plot usati nella pubblicazione 6pioni E687
 che si possono prendere da http://arXiv.org/abs/hep-ex/0106029
 (se clicchi source ti da' tutto il tar file comprende
 quindi anche i singoli postscript). 
 2) utilizzare la stessa catena di kumac per girare i
 dati  6pioni focus
 3) ottimizzare i tagli sui dati focus e pubblicare



 1) Riprodurre i plot 6pioni E687
  la kumac che produce i plot usati per la pubblicazione
  (piu' altri) si trova in
  /afs/..../e831/user/adriano/lebrun3/seie68712_2001.kumac
  le ntuple di input sono
   set6pi1990_500000.nut
             _500k-10000k.ntu
             _1000k-last.ntu
   set6pi1991_500000.ntu
             _500k-10000k.ntu
             _1000k-last.ntu   

 
 2) Dati E831
 per la produzione degli istogrammi a partire dall'ntupla
 vale quanto sopra.
 Per il fit la directory e'
 /afs/lnf.infn.it/project/e831/user/adriano/6pioni831/fit
  adriano ha scritto delle istruzioni nel file
  aaareadme.txt che sono anche utilizzabili per il
  fit E687, e che riporto in fondo.












 ============================aaareadme.txt
Informazioni sul direttorio: 8-Agosto-2001

Questo direttorio viene usato per fittare i dati di E831 contenuti 
negli istogrammi del file e831_1-19_21-38_01.his che si trova nel
direttorio superiore. La kumac e831.kumac produce il file e831.ps
con gli istogrammi piu' interessanti.
E' possibile fittare  gli spettri di massa 6 pioni con 
con ID=6003(5MeV/bin), 7003(10MeV/bin),8003(15MeV/bin). 
A questo scopo si traduce l'istogramma in file di dati (x,y,dy) con
il programma garbwrite_e831.kumac, ottenendo il file 6003.dat. 
Si puo' usare questo programma per creare 7003.dat e 8003.dat.
Per fittare questi dati si usa il programma fit1y.f di Garbincius.
Attualmente e' settato per il file 6003.dat (300 punti).Agire su read-data
per cambiare il file.Volendo analizzare 7003.dat e 8003.dat bisogna 
inoltre dimensionare opportunamente i connon di input in fit1y.f
Il comando link_cern fit1y produce fit1y.exe che poi si esegue.
L'intervallo di fit viene dato dai parametri nbeg, nend che si settano nella
routine fnct. In numero dei punti npt viene settato nella routine fit.
In questa routine viene fatto il settaggio dei parametri.
In fit viene chiamato migrad e poi minos.
I parametri del fit sono nel file di uscita fort.8. Scrivendo materialmente
questi dati nella function garb_e831.f e runnando la kumac
fig6_garb_e831.kumac si produce l'istogramma per e831 equivalente,
a parte la sottrazione del fondo, e l'eliminazione della risoluzione,
della figura 6.