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 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 o n t r i b u t i o n s f r o m: D a n i e l C u s s o l ( L P C C a e n ) R o l a n d D a y r a s , E n r i c o d e F i l i p p o ( 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 o 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 o r j e a n , M a u r y c y R e j mu n d ( G A N I L ) , . . . ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14
KaliVeda Toolkit I N T R O D U C T I O N I N D R A i n D 5 ● S i mu l a t i o n & a n a l y s i s o f h e a v y - i o n c o l l i s i o n d a t a i n 5 - 1 0 0 M e V / n u c l e o 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 o r k b a s e d o n R O O T – ~450 classes in v1.9.2 ● O p e n - s o u r c e c o l l a b o r a t i v e p r o j e c t – distributed version-control (bazaar VCS) – web-based code-hostjng, bug reports etc. (launchpad.net/kaliveda) ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14
KaliVeda Toolkit I N T R O D U C T I O N ● F i r s t d e v e l o p e d f o r a n a l y s i s & ma n a g e me n t o f I N D R A d a t a ( 2 0 0 2 - 2 0 0 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 o 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 0 1 0 - . . . ) – FAZIA, arbitrary detector geometries, energy-loss calculatjons, ... ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14
Main Ingredients I N T R O D U C T I O N ● N u c l e i & e v e n t s ( c o l l e c t i o n s o f n u c l e i ) 1 . N U C L E I & E V E N T S – nuclear propertjes, reactjon kinematjcs, systematjcs ● D e t e c t o r s i mu l a t i o n 2 . D E T E C T O R S I M U L A T I O 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 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 – Z & A identjfjcatjon grids, semi-automatjc identjfjcatjon ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14
1 . N U C L E I Nuclei & E V E N T S N u c l e a r p r o p e r t i e s f r o 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 All nuclei in A l s o a v a i l a b l e : tables (~3000) 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 0 1 3 ) 6 9 _ 9 5 Nuclei with t 1/2 >1000s. KVNuclearChart http://www.nndc.bnl.gov/nudat2 ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14
1 . N U C L E I Events & E V E N T S C o l l e c t i o n s o 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 o r e n t z i a n b o o s t / r o t a t i o n C o n t a i n e r f o 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 o n T S e l e c t o r f r a me w o r k ( = > c a n u s e P R O O F ) E n t r y - p o i n t f o r c o n v e r s i o n o f KVEventViewer e x i s t i n g d a t a ( e . g . s i mu l a t i o n ) i n o r d e r t o u s e t o o l k i t ( e . g . f i l t e r w i t h d e t e c t o r r e s p o n s e ) ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14
1 . N U C L E I Reaction systematics & E V E N T S S o me u s e f u l s y s t e ma t i c s o n k i n e ma t i c s o f h e a v y - i o n r e a c t i o 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) ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14
1 . N U C L E I Reaction systematics & E V E N T S S o me u s e f u l s y s t e ma t i c s o n k i n e ma t i c s o f h e a v y - i o n r e a c t i o n s : root [3] twob.GetELabVsThetaLabFunc(3)->Draw() root [4] twob.GetELabVsThetaLabFunc(4)->Draw("same") Elastjc scatuering kinematjcs Based on cibi.for (?), an old code from Saclay (?)... (via M. Morjean - GANIL) ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14
1 . N U C L E I Reaction systematics & E V E N T S S o me u s e f u l s y s t e ma t i c s o n k i n e ma t i c s o f h e a v y - i o n r e a c t i o n s : root [0] KVNucleus proj("34Si") root [1] KVNucleus targ("132Sn") W.W. Wilcke et al., root [2] KVWilckeReactionParameters wrp(proj,targ) At. Data and Nucl. Data root [3] wrp.Print() Tables 25 (1980) 389-619 ------------------------------------------- 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.] ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14
2 . D E T E C T O R Detector geometries S I M U L A T I O N R O O T g e o me t r y p a c k a g e u s e d f o r mo d e l l i n g a n d p r o p a g a t i o n o f e v e n t s / n u c l e i KVINDRAFAZIA FAZIA demonstrator INDRA (rings 6-17) ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14
2 . D E T E C T O R Detector geometries S I M U L A T I O N A n y R O O T g e o me t r y c a n b e i mp o r t e d i n t o t h e f r a me w o r k Silicon wafer 5x5cm 2 KaliVeda class handling detector geometries root [2] KVGeoImport gimp(gGeoManager, KVMaterial::GetRangeTable(), new KVMultiDetArray) ROOT geometry default range table 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 ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14
2 . D E T E C T O R Detector geometries S I M U L A T I O N A n y R O O T g e o me t r y c a n b e i mp o r t e d i n t o t h e f r a me w o r k Silicon wafer 5x5cm 2 root [5] KVDetector* det = gMultiDetArray->GetDetector("SILICON_1") for(int i=0;i<100000;i++){ TVector3 point = det->GetRandomPointOnEntranceWindow(); ... } ICC/Spiral2 Common Analysis Tools WS, IPNO, 26/11/14
Recommend
More recommend