KaliVeda Toolkit C o r e d e v e l o p me n t t e a - - PowerPoint PPT Presentation

kaliveda toolkit
SMART_READER_LITE
LIVE PREVIEW

KaliVeda Toolkit C o r e d e v e l o p me n t t e a - - PowerPoint PPT Presentation

KaliVeda Toolkit C o r e d e v e l o p me n t t e a m: J o h n F r a n k l a n d , E r i c B o n n e t ( G A N I L ) D i e g o G r u y e r ( G A N I L / I N F N F i r e n z


slide-1
SLIDE 1

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

KaliVeda Toolkit

C

  • r

e d e v e l

  • p

me n t t e a m: J

  • h

n F r a n k l a n d , E r i c B

  • n

n e t ( G A N I L ) D i e g

  • G

r u y e r ( G A N I L / I N F N F i r e n z e ) G u i l a i n A d e ma r d ( I P N O r s a y )

I n c l u d i n g c

  • n

t r i b u t i

  • n

s f r

  • m:

D a n i e l C u s s

  • l

( L P C C a e n ) R

  • l

a n d D a y r a s , E n r i c

  • d

e F i l i p p

  • (

C E A / S P h N S a c l a y / I N F N C a t a n i a ) L a u r e n t T a s s a n

  • G
  • t

, M a r i e

  • F

r a n c e R i v e t ( I P N O r s a y ) M a u r i c e M

  • r

j e a n , M a u r y c y R e j mu n d ( G A N I L ) , . . .

slide-2
SLIDE 2

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

KaliVeda Toolkit

I N D R A i n D 5 I N T R O D U C T I O N

  • S

i mu l a t i

  • n

& a n a l y s i s

  • f

h e a v y

  • i
  • n

c

  • l

l i s i

  • n

d a t a i n 5

  • 1

M e V / n u c l e

  • n

e n e r g y r a n g e

– charged partjcle detectjon & identjfjcatjon (from protons to uranium...) – no treatment for neutrons or γ (only identjfjcatjon veto)

  • C

+ + f r a me w

  • r

k b a s e d

  • n

R O O T

– ~450 classes in v1.9.2

  • O

p e n

  • s
  • u

r c e c

  • l

l a b

  • r

a t i v e p r

  • j

e c t

– distributed version-control (bazaar VCS) – web-based code-hostjng, bug reports etc. (launchpad.net/kaliveda)

slide-3
SLIDE 3

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

KaliVeda Toolkit

  • F

i r s t d e v e l

  • p

e d f

  • r

a n a l y s i s & ma n a g e me n t

  • f

I N D R A d a t a ( 2 2

  • 2

9 )

– regroup existjng tools for simulatjon, data reductjon, data analysis – unifjed analysis environment for database of ~100 difgerent

reactjons (1993-2011...)

– including INDRA-VAMOS coupling (2007)

  • G

e n e r a l i s e d f

  • r

n e w a n d f u t u r e c h a r g e d

  • p

a r t i c l e a r r a y s ( 2 1

  • .

. . )

– FAZIA, arbitrary detector geometries, energy-loss calculatjons, ...

I N T R O D U C T I O N

slide-4
SLIDE 4

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Main Ingredients

  • N

u c l e i & e v e n t s ( c

  • l

l e c t i

  • n

s

  • f

n u c l e i )

– nuclear propertjes, reactjon kinematjcs, systematjcs

  • D

e t e c t

  • r

s i mu l a t i

  • n

– geometry, energy-loss calculatjons

  • P

a r t i c l e i d e n t i f i c a t i

  • n

– Z & A identjfjcatjon grids, semi-automatjc identjfjcatjon

I N T R O D U C T I O N 1 . N U C L E I & E V E N T S 2 . D E T E C T O R S I M U L A T I O N 3 . P A R T I C L E I D E N T I F I C A T I O N

slide-5
SLIDE 5

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Nuclei

1 . N U C L E I & E V E N T S All nuclei in tables (~3000) Nuclei with t1/2>1000s.

N u c l e a r p r

  • p

e r t i e s f r

  • m

N N D C

  • N

U D A T 2 d a t a b a s e : a b u n d a n c e , l i f e t i me / w i d t h , ma s s e x c e s s

http://www.nndc.bnl.gov/nudat2

A l s

  • a

v a i l a b l e : g . s . c h a r g e r a d i i

A t . D a t a N u c l . D a t a T a b . 9 9 ( 2 1 3 ) 6 9 _ 9 5 KVNuclearChart

slide-6
SLIDE 6

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Events

1 . N U C L E I & E V E N T S

C

  • l

l e c t i

  • n

s

  • f

n u c l e i w i t h f u l l r e l a t i v i s t i c k i n e ma t i c s

  • h

a n d l e L

  • r

e n t z i a n b

  • s

t / r

  • t

a t i

  • n

KVEventViewer

C

  • n

t a i n e r f

  • r

w r i t i n g d a t a i n R O O T T r e e s

  • u

s e r a n a l y s i s b a s e d

  • n

T S e l e c t

  • r

f r a me w

  • r

k ( = > c a n u s e P R O O F )

E n t r y

  • p
  • i

n t f

  • r

c

  • n

v e r s i

  • n
  • f

e x i s t i n g d a t a ( e . g . s i mu l a t i

  • n

) i n

  • r

d e r t

  • u

s e t

  • l

k i t ( e . g . f i l t e r w i t h d e t e c t

  • r

r e s p

  • n

s e )

slide-7
SLIDE 7

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Reaction systematics

1 . N U C L E I & E V E N T S

S

  • me

u s e f u l s y s t e ma t i c s

  • n

k i n e ma t i c s

  • f

h e a v y

  • i
  • n

r e a c t i

  • n

s :

root [0] KV2Body twob(new KVNucleus("181Ta",39.),new KVNucleus("66Zn")) root [1] twob.CalculateKinematics() root [2] twob.Print()

***** REACTION 181Ta + 66Zn ---> 181Ta + 66Zn ******

E.LAB = 7059 MEV QGG = 0 MEV E.EXC = 0 MEV ==> Q-REACTION = 0 MEV AVAILABLE ENERGY IN C.M. : ECM = 1877.39 MEV (7.60079 MEV/A) PROJECTILE VELOCITY IN LAB 8.41439 CM/NS ( 0.280674 * C ) VELOCITY OF C.M. 6.23431 CM/NS ENERGY - VELOCITY OF NUCLEUS 1 IN CM : 504.887 MEV 2.31522 CM/NS (K=2.69275) ENERGY - VELOCITY OF NUCLEUS 2 IN CM : 1372.51 MEV 6.23431 CM/NS (K=1) ENERGY - VELOCITY OF NUCLEUS 3 IN CM : 485.656 MEV 2.27089 CM/NS (K=2.74531) ENERGY - VELOCITY OF NUCLEUS 4 IN CM : 1319.18 MEV 6.11238 CM/NS (K=1.01995) MAXIMUM SCATTERING ANGLE IN LABORATORY THETA #3# 20.9357 DEG. THETA #4# 78.4028 DEG. GRAZING ANGLE IN LABORATORY : PROJECTILE 2.008 DEG. GRAZING ANGLE IN LABORATORY : TARGET 86.2389 DEG.

Based on cibi.for (?), an old code from Saclay (?)... (via M. Morjean - GANIL)

slide-8
SLIDE 8

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Reaction systematics

1 . N U C L E I & E V E N T S

S

  • me

u s e f u l s y s t e ma t i c s

  • n

k i n e ma t i c s

  • f

h e a v y

  • i
  • n

r e a c t i

  • n

s :

root [3] twob.GetELabVsThetaLabFunc(3)->Draw() root [4] twob.GetELabVsThetaLabFunc(4)->Draw("same")

Based on cibi.for (?), an old code from Saclay (?)... (via M. Morjean - GANIL) Elastjc scatuering kinematjcs

slide-9
SLIDE 9

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Reaction systematics

1 . N U C L E I & E V E N T S

S

  • me

u s e f u l s y s t e ma t i c s

  • n

k i n e ma t i c s

  • f

h e a v y

  • i
  • n

r e a c t i

  • n

s :

root [0] KVNucleus proj("34Si") root [1] KVNucleus targ("132Sn") root [2] KVWilckeReactionParameters wrp(proj,targ) root [3] wrp.Print()

  • PARAMETERS INDEPENDENT OF BOMBARDING ENERGY
  • ATOMIC NUMBERS: ZP= 14 ZT= 50 ZC= 64 (166Gd)

NEUTRON NUMBERS: NP= 20 NT= 82 NC=102 ... [etc. etc.] LIQUID DROP PARAMETERS: GAMMA= 0.863 MeV/fm**2 PROX-FACTOR= 22.98 MeV FUSION RELATED PARAMETERS: R-BARRIER=10.97 fm V(RB)= 85.2 MeV L-CRITICAL= 106 HBAR

  • EL/u ELAB ECM ECM/VC k ETA LMAX SGMAR SGFUS QP-CM QP-LP EP-QP
  • 1.0 34 27 0.33 5.9 110.2 0 0 0 180.0 -0.0 12

2.0 68 54 0.65 8.4 77.9 0 0 0 180.0 -0.0 24 3.0 102 81 0.98 10.2 63.6 0 0 0 180.0 -0.0 36 4.0 136 108 1.31 11.8 55.1 69 1100 802 76.8 63.5 102 4.5 153 122 1.47 12.5 52.0 86 1498 1132 62.2 50.7 126 ... [etc. etc.]

W.W. Wilcke et al.,

  • At. Data and Nucl. Data

Tables 25 (1980) 389-619

slide-10
SLIDE 10

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Detector geometries

2 . D E T E C T O R S I M U L A T I O N

R O O T g e

  • me

t r y p a c k a g e u s e d f

  • r

mo d e l l i n g a n d p r

  • p

a g a t i

  • n
  • f

e v e n t s / n u c l e i

KVINDRAFAZIA

FAZIA demonstrator INDRA (rings 6-17)

slide-11
SLIDE 11

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Detector geometries

2 . D E T E C T O R S I M U L A T I O N

A n y R O O T g e

  • me

t r y c a n b e i mp

  • r

t e d i n t

  • t

h e f r a me w

  • r

k

Silicon wafer 5x5cm2

root [2] KVGeoImport gimp(gGeoManager, KVMaterial::GetRangeTable(), new KVMultiDetArray) root [3] gimp.ImportGeometry() Info in <KVGeoImport::ImportGeometry>: Tested 650161 directions - Theta=[0.000000,180.000000:0.100000] Phi=[0.000000,360.000000:1.000000] Info in <KVGeoImport::ImportGeometry>: Imported 1 detectors into array Info in <KVGeoImport::ImportGeometry>: Created 0 identification telescopes in array

default range table ROOT geometry KaliVeda class handling detector geometries

slide-12
SLIDE 12

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Detector geometries

2 . D E T E C T O R S I M U L A T I O N

A n y R O O T g e

  • me

t r y c a n b e i mp

  • r

t e d i n t

  • t

h e f r a me w

  • r

k

Silicon wafer 5x5cm2

root [5] KVDetector* det = gMultiDetArray->GetDetector("SILICON_1") for(int i=0;i<100000;i++){ TVector3 point = det->GetRandomPointOnEntranceWindow(); ... }

slide-13
SLIDE 13

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Energy loss/range tables

2 . D E T E C T O R S I M U L A T I O N

A n a l y t i c i n t e r p

  • l

a t i

  • n
  • f

N

  • r

t h c l i f f e / S c h i l l i n g & H u b e r t / B i mb

  • t

/ G a u v i n t a b l e s

Nuclear Data Tables A7 (1970) 233 Atomic Data and Nuclear Data Tables 46 (1990) 1

slide-14
SLIDE 14

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Energy loss/range tables

2 . D E T E C T O R S I M U L A T I O N

A n a l y t i c i n t e r p

  • l

a t i

  • n
  • f

N

  • r

t h c l i f f e / S c h i l l i n g & H u b e r t / B i mb

  • t

/ G a u v i n t a b l e s

Nuclear Data Tables A7 (1970) 233 Atomic Data and Nuclear Data Tables 46 (1990) 1

slide-15
SLIDE 15

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Detection simulation

2 . D E T E C T O R S I M U L A T I O N

G e

  • me

t r y + r a n g e t a b l e s = d e t e c t

  • r

r e s p

  • n

s e f

  • r

s i mu l a t e d e v e n t s

gMultiDetArray->DetectEvent(KVEvent*) "Large" angle particles hitting frame inside ionisation chamber ∆E<∆Etheo Collimated IoCh response ∆E=∆Etheo

θ<5°

IoCh: 2cm C4H10 2.5µm windows aluminium frame ∆E in gas ∆E in gas

slide-16
SLIDE 16

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Z & A Identification

3 . P A R T I C L E I D E N T I F I C A T I O N

T

  • l

s f

  • r

c h a r g e d p a r t i c l e i d e n t i f i c a t i

  • n

i n ∆E

  • E
  • r

P S A ma p s

Linearisation algorithms:

  • L. Tassan-Got

(IPNO) KVIDZAGrid KVIDGCsI

slide-17
SLIDE 17

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Z & A Identification

3 . P A R T I C L E I D E N T I F I C A T I O N

∆E

  • E

g r i d s c a n b e g e n e r a t e d s e mi

  • a

u t

  • ma

t i c a l l y . . .

  • D. Gruyer, PhD Thesis (2014)
slide-18
SLIDE 18

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Z & A Identification

3 . P A R T I C L E I D E N T I F I C A T I O N

...t h e n e x t r a p

  • l

a t e d t

  • c
  • v

e r a l l d a t a

  • D. Gruyer, PhD Thesis (2014)

Automatjc grid Extrapolatjon Identjfjcatjon

Extrapolatjon via fjt using Tassan-Got functjonals

"A new functional for charge and mass identification in ∆E−E telescopes",

  • Nucl. Instr. and Meth. in Phys. Res.

B194 (2002) 503

slide-19
SLIDE 19

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Z & A Identification

3 . P A R T I C L E I D E N T I F I C A T I O N

F

  • r

h i g h

  • r

e s

  • l

u t i

  • n

d e t e c t

  • r

s ( e . g . F A Z I A ) , a s e c

  • n

d l i n e a r i s a t i

  • n

a l l

  • w

s a u t

  • ma

t i c i s

  • t
  • p

i c i d e n t i f i c a t i

  • n

First linearisatjon Second linearisatjon

  • D. Gruyer, PhD Thesis (2014)
slide-20
SLIDE 20

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Other tools (GUI)

M I S C E L L A N Y

G r a p h i c a l i n t e r f a c e s f

  • r

i d e n t i f i c a t i

  • n

g r i d s

slide-21
SLIDE 21

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Other tools (GUI)

M I S C E L L A N Y

E x p e r i me n t d a t a b a s e ma n a g e r

slide-22
SLIDE 22

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Other tools (GUI)

M I S C E L L A N Y

I mp r

  • v

e d ( ? ) T T r e e b r

  • w

s e r / a n a l y s i s G U I

slide-23
SLIDE 23

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Other tools (GUI)

M I S C E L L A N Y

S i mu l a t i

  • n

f i l t e r i n g / a n a l y s i s G U I

slide-24
SLIDE 24

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

Other tools (GUI)

M I S C E L L A N Y

U s e r a n a l y s i s s u b mi s s i

  • n

G U I i n t e r f a c e d w i t h C C I N 2 P 3 b a t c h s y s t e m

slide-25
SLIDE 25

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

KaliVeda Toolkit

I N D R A i n D 5 S U M M A R Y

  • F
  • r

mo r e i n f

  • r

ma t i

  • n

:

http://indra.in2p3.fr/KaliVedaDoc

slide-26
SLIDE 26

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

KaliVeda Toolkit

I N D R A i n D 5 S U M M A R Y

  • F
  • r

mo r e i n f

  • r

ma t i

  • n

:

https://launchpad.net/kaliveda

*or links on main website

slide-27
SLIDE 27

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

KaliVeda Toolkit

I N D R A i n D 5 S U M M A R Y

  • A

v a i l a b l e f

  • r

d

  • w

n l

  • a

d a s :

– release tarballs – bazaar branches – Ubuntu packages (experimental)

  • S

u p p

  • r

t e d p l a t f

  • r

ms :

– Linux 32/64 bit (tested Ubuntu/Fedora/Scientjfjc Linux) – original support for MacOSX & Windows abandoned, currently

no developers working on these platgorms

slide-28
SLIDE 28

ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14

KaliVeda Toolkit

  • P

r

  • s

:

– regroup many useful tools for analysis & simulatjon of HIC in ~GANIL energy

range

– contjnual development/improvement through interactjon with users (&

developers are also users!)

  • C
  • n

s :

– only charged partjcles are treated (INDRA/FAZIA data focus) – default range tables validated for GANIL energies

  • b

u t c

  • u

l d b e r e p l a c e d ( A B C i n t e r f a c e )

  • R

O O T V i r t u a l M

  • n

t e C a r l

  • a

l s

  • p
  • s

s i b l e ( = > G 3 , G 4 , . . . )

– for the moment, problems with ROOT v6 (in progress)

S U M M A R Y