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

Re: PROBLEMA trasporto secondarie

Ciao Monica, 

il problema era stato riscontrato da Francesca e da Barbara forse un paio
di anni fa con il passaggio delle cernlib da 2001 a 2002. Il problema
nasce in quanto e' stato cambiato un parametro di geant, MXGKIN da 100 a
250. Questo crea un problema nel programma fidamc in runtime. In pratica
il nostro codice e' stato scritto per accettare 100 come parametro della
cinematica e non 250.  La conseguenza e' che non vengono copiate le
particelle secondarie dal trasporto/tracking alla cinematica. Il problema
non viene riscontrato ricostruendo i dati di Montecarlo o guardando
l'interattivo di Geant ma solo controllando la presenza nelle banche della
cinematica delle particelle secondarie.

Per risolvere il problema si possono ricompilare le cernlib dopo aver
modificato il file geant321/gckmax.inc settando PARAMETER (MXGKIN=100)
invece di 250. Questo ovviamente e' gia' stato fatto sulle cernlib 2002
delle macchine della produzione di Finuda in counting room per la
distribuzione RH 9.0: 

fnddiskfe:/data02/cern-rh90

non mi ricordo di averlo fatto sulle 2003. 

Su Alpha funziona in quanto usa le cernlib 2001 o precedenti, su linux
gentoo posso tranquillamente usare le 2001. Le distribuzioni RedHat o
derivate dalla 9.0 in poi non possono ricompilare ed utilizzare le 2001.

All'epoca avevo spedito un mail ad Aldo (non c'era la mailing list
offline) forwardando una botta e risposta che ho avuto con Mc Lauren che
mantiene il codice di geant3 e delle cernlib. Se non sbaglio avevo fatto 
circolare un e-mail con tutti gli upgrade da fare su RH9.0 per risolvere 
il problema. Ne avevo discusso anche con Luigi che e' un esperto di 
ricompilazione di cernlib (se non mi sbaglio mi aveva aiutato in croom).

On Fri, 28 Oct 2005, Germano Bonomi wrote:

> 
> Ciao Monica,
> noi abbiamo Fedora Core 3, cern 2003, gcc 3.4.4-2 e il MC mi funziona 
> bene.
> 

Ottimo, ora siamo piu' tranquilli, possiamo installare Fedora 3 o distro
derivate come SLC per le nostre simulazioni, ma hai controllato la
presenza delle secondarie nella banca della cinematica?

ciao!
Stefano


> Ciao, Germano
> 
> On Fri, 28 Oct 2005, Monica Bertani wrote:
> > 
> > Caro Stefano,
> > noi usiamo la RH9, cernlib2003.
> > 
> > Il problema e' che dopo aver generato eventi con il MC, la routine everec 
> > non riempie l'ntupla 110  dalla var n.4 in poi, con i valori delle
> > tracce generate. 
> > 
> > Ad es. la stampa delle variabili          
> >          JK     = LQ(JKINE -2)        ! pointer to K- KINE bank
> >          IVERT  = IFIX (Q(JK + 8))     ! pi- vertex number
> >          JV     = LQ(JVERTX - IVERT)   ! pointer to pi- vertex bank
> >          ITRACK = IFIX (Q(JV + 8))     ! pi- track number
> >          THN    = Q(JV+4)              ! pi- vertex time = HN formation time
> > su linux e' per ogni evento: 
> > jk,ivert,jv,itrack,thn= 4601784 0 0 0  0.0058546993
> > 
> > Ho fatto il test su alpha e in effetti, come dice Francesca, li' le stesse
> > variabili sono piu' sensate e l'ntupla e' riempita correttamente, almeno
> > la parte del pi-. 
> > 
> > Se tu o altri  avete avuto questo problema e lo avete risolto fatelo
> > sapere per favore, 
> > perche' fare generazioni di numeri elevati di eventi su alfa non e'
> > possibile! 
> > 
> > grazie,
> > ciao  Monica
> > 
> > ============================================================================
> > Dr. Monica Bertani                       
> > Laboratori Nazionali di Frascati         phone : +39 06 94032760
> > dell'INFN - Alte Energie                 fax   : +39 06 94032427
> > Via Enrico Fermi, 40                     e-mail: monica.bertani@lnf.infn.it
> > 00044 Frascati (Italy)
> > ============================================================================
> > 
> > 
> > On Wed, 26 Oct 2005, Stefano Piano wrote:
> > 
> > > Ciao,
> > > riuscite a mandarmi una descrizione del problema? Non vedete le tracce in 
> > > ricostruzione della simulazione o non vedete le tracce nelle banche di 
> > > transporto di geant ?
> > > Per piacere specificate anche quale versione di Linux, cernlib e di 
> > > compilatore state usando. 
> > > saluti,
> > > Stefano
> > > 
> > > On Wed, 26 Oct 2005, Francesca De Mori wrote:
> > > 
> > > > Ciao Monica,
> > > > tu mi dici che a Frascati non riuscite a vedere lo spettro generato col
> > > > MC.
> > > > Io mi ricordo di aver riscontrato un problema simile  nel passare alle
> > > > CERNLIB2003(ma c'era anche per le 2002 almeno du RH9), insieme a Barbara e
> > > > Stefano.
> > > > Avendo riscontrato dei problemi e avendo a disposizione delle alpha qui a
> > > > Torino ho sempre girato i MC su alpha.
> > > > Sono sicura che pero' Barbara e Stefano avevano risolto il problema e
> > > > quindi senza girarti i mail che sono circolati tra di noi ho pensato che
> > > > sia piu' efficace chiedere a loro come hanno risolto il problema.
> > > > Invio questo mail a tutto il gruppo perche' sia a tutti presente (se
> > > > ancora esiste,
> > > > ma tu mi dici di si') e cosi' uno non debba scervellarsi per capire cosa
> > > > non funziona. Probabilmente anche Germano non ha questo problema(giri su
> > > > linux vero?).
> > > > Direi anzi che sarebbe utile che una comunicazione di questo tipo non sia
> > > > punto a punto, ma magari pubblica sulla pagina web dell'offline.
> > > > Bug fixes e simili potrebbero essere depositati da qualche parte perche'
> > > > siano disponibili per tutti...gia' siamo pochi se poi ciascuno deve
> > > > dedicare del tempo a risolvere problemi gia' noti...e' un po' dura.
> > > > 
> > > > Fra
> > > > 
> > > >   Francesca De Mori
> > > >   Dipartimento di Fisica Sperimentale
> > > >   Universita' degli Studi di Torino
> > > >   via Pietro Giuria 1
> > > >   10125 Torino
> > > >   Italy
> > > >   tel 0039 0116707323
> > > >   fax 0039 0116707325
> > > > 
> > > > 
> > > 
> > > -- 
> > > Stefano Piano
> > > INFN sez. Trieste
> > > Via A. Valerio, 2
> > > I-34127   Trieste
> > > --
> > > tel. +39 040 5583386
> > > fax. +39 040 5583350
> > > mob. +39 349 3534904
> > > --
> > > stefano.piano@ts.infn.it
> > > 
> > 
> > 
> 
> 

-- 
Stefano Piano
INFN sez. Trieste
Via A. Valerio, 2
I-34127   Trieste
--
tel. +39 040 5583386
fax. +39 040 5583350
mob. +39 349 3534904
--
stefano.piano@ts.infn.it
Follow-Ups: