Come usare OSCAR
Tommaso Boccali
SNS Pisa
4/9/2003
TISB
Outline

Cosa spiego:





Dove trovare OSCAR
Come lanciarlo su ntuple CMKIN / particle gun
Come scegliere la fisica
Come scegliere la geometria
Cosa non spiego:


Validazione
Confronto di fisica fra Cmsim e OSCAR
4/9/2003
TISB
2
Stato Attuale

OSCAR è in fase avanzata di validazione


Traduzione: non abbiamo ancora controllato tutto
ma quello che è stato controllato è ok.
OSCAR è integrato con ORCA_7_3_x

Questo vuol dire che non funziona ancora con
POOL, ma con plain ROOT
4/9/2003
TISB
3
Scema operativo
The “old” way
(Geant 3)
The “new” way
(Geant 4)
Phys Sim
Pythia
HEPEVT
ntuple
Pythia
HEPEVT
ntuple
Det Sim
CmSim
.fz events
.rz geometry
OSCAR
.fz events
.rz geometry
Conversion
to OO
ORCA
OBJY
ooHits
ORCA
OBJY
ooHits
Digitization
ORCA
OBJY
Digis
ORCA
OBJY
Digis
ORCA
Ntuple
Tags
ORCA
Ntuple
Tags
Analysis
4/9/2003
TISB
4
Dove trovare OSCAR

Inutile dirlo, no?


scram list OSCAR
Se volete giocarci, al momento usate la 2_3_2
scram project OSCAR OSCAR_2_3_2
 cd OSCAR_2_3_2/src
 project OSCAR
 cvs co –r OSCAR_2_3_2
OscarApplication/G4SimApplication/test


OscarApplication/G4SimApplication/test contiene
tutti gli script, e l’eseguibile!
4/9/2003
TISB
5
Cosa serve per girare?
1.
2.
Una geometria (CMS, per esempio?)
Un evento di fisica
1.
2.
3.
CMKIN ntuple
Particle gun per studi veloci
Tante altro cose volendo fare studi specifici
1.
2.
3.
4.
Set di tagli per i materiali
La possibilità di spegnere i calorimetri
Restringere le zone in eta
Ecc ecc
4/9/2003
TISB
6
Geometria

Facciamola semplice, volete girare tutto CMS
con l’ultima geometria disponibile

Non fate nulla!
4/9/2003
TISB
7
Cosa modificare
Alla fine tecnicamente il file da modificare è
solo 1 (o al max 2): .orcarc
 Linee rilevanti


CMSRandom:Seeds = 40 3


OutputDataSet = /System/TestG4/testOSCAR


Solita sequenza iniziale di numeri casuali
Ricordatevi! Scriviamo direttamente un ROOT DB!
GoPersistent = 1

4/9/2003
Immagino che useremo sempre la persistenza!
Altrimenti, OSCAR può essere usato direttamente per
generare istogrammi ecc ecc
TISB
8
Senza persistenza…
Per esempio, per
calcolare il Material
Budget, a che serve
girare ORCA?
4/9/2003
TISB
9
.orcarc

ExtraPackages =
CMSProductionPhysics:MantisNtplReader

1.
2.

4/9/2003
Questa è la prima selezione importante; sono
disponibili 2 opzioni
MantisNtplReader: serve per leggere CMKIN ntuples
MantisGunReader: serve per accendere il particle gun
e sparare singole particelle
A questo punto le possibilità si biforcano
TISB
10
MantisNtplReader

In generale questo è il production mode, le
card rilevanti sono

EventNtplReader:NtplFileName =
/afs/cern.ch/cms/geant4rep/genntpl/minbias.ntpl


EventNtplReader:NtplID = 101


Numero della ntupla nel file
HepEventCmkinNtupleReader:FirstEvent = 0


Nome della ntupla
Da quale evento cominciare
NumberOfEventsToBeProcessed = 10

4/9/2003
Indovinate…
TISB
11
MantisNtplReader

Consiglio spassionato:

Se volete fare questo, lasciate stare tutte le altre
opzioni! O(100)
4/9/2003
TISB
12
MantisGunReader

Uso più specialistico, 2 possibili casi
E Fissata
ParticleGun:Gun =
FlatRandomEGun
ParticleGun:PID = 13
ParticleGun:MinEta = 0.
ParticleGun:MaxEta = 2.5
ParticleGun:MinPhi =
ParticleGun:MaxPhi =
ParticleGun:MinE = 10
ParticleGun:MaxE = 10
Pt Fissato
ParticleGun:Gun =
FlatRandomPtGun
ParticleGun:PID = 13
ParticleGun:MinEta = 0.
ParticleGun:MaxEta = 2.5
ParticleGun:MinPhi =
ParticleGun:MaxPhi =
ParticleGun:MinPt = 19.8
ParticleGun:MaxPt = 20.2
4/9/2003
TISB
13
Inoltre, (volendo)

Cambiare tipo di smearing del primario
ParticleGun:VertexGenerator =
FlatEventVertexGenerator
VertexGenerator:MinX = 0.
VertexGenerator:MinY = 0.
VertexGenerator:MinZ = 0.
VertexGenerator:MaxX = 0.
VertexGenerator:MaxY = 0.
VertexGenerator:MaxZ = 0.
4/9/2003
TISB
ParticleGun:VertexGenerator =
GaussianEventVertexGenerator
VertexGenerator:MeanX = 0.
VertexGenerator:MeanY = 0.
VertexGenerator:MeanZ = 0.
VertexGenerator:SigmaX = 0.015
VertexGenerator:SigmaY = 0.015
VertexGenerator:SigmaZ = 53.0
14
Mi fermo qui

Sono possibili altri modi






Interattivo (comandi interattivi!)
Visualizzazione di G4
Macro mode (come interattivo ma da file)
Prerelease con POOL
…
Cehdete e vi sarà spiegato …
4/9/2003
TISB
15
Run!
`eval scram runtime –csh`
 source writeTrigger.csh




Semplicemente, setta DETINPUT all’ultima
geometria disponibile, dovrebbe essere la
cms133.rz adesso
Serve solamente per leggere il campo magnetico,
che ancora non è disponibile in C++
oscar
4/9/2003
TISB
16
Per rileggere, semplicemente
Usare ORCA_7_3_0
 Fare tutto nel modo standard, semplicemente
settare


InputCollections =
/System/TestG4/testOSCAR//testOSCAR


(nel caso in esame ora)
Ora, quali sono le performance che possiamo
aspettarci? (dal talk di Tony)
4/9/2003
TISB
17
Performance?
4/9/2003
TISB
18
Physics channels - timing
Sample
CMS132 or
CMS133
OSCAR 232pre3 or
OSCAR 232pre4
udsg_pt170
361
597
ZZ->4e
260
535
(no info)
780
bt_b_pt170up
173*
432
bt_c_pt30-40
(no info)
277
eg_jets_2g_pt170u
p
*) Normalisation not certain due to running conditions
4/9/2003
TISB
19
Summary





OSCAR 232pre3(+) suitable for Production running
Production has started at outside institutes (Wisconsin, Pisa,
hopefully others too by now)
OSCAR is stable enough, physics events show ~95% success
for jobs (1/5000 events fail)
Performance is acceptable, factor 2 w.r.t. CMSIM confirmed in
several channels
Production will continue with 232pre4 to gain experience and
shake out bugs.
4/9/2003
TISB
20