Fair simulation analysis framework fairroot
This presentation is the property of its rightful owner.
Sponsored Links
1 / 33

FAIR Simulation & Analysis Framework FairRoot PowerPoint PPT Presentation


  • 83 Views
  • Uploaded on
  • Presentation posted in: General

FAIR Simulation & Analysis Framework FairRoot. M. Al-Turany, D. Bertini, F. Uhlig GSI-IT. FairRoot FairRoot new features Integrated Track follower (Geane) Geant4 configuration classes New reader (ROOT converted CAD Step Format ) Fast Simulation CMake/CTest

Download Presentation

FAIR Simulation & Analysis Framework FairRoot

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript


Fair simulation analysis framework fairroot

FAIR Simulation & Analysis FrameworkFairRoot

M. Al-Turany, D. Bertini, F. Uhlig

GSI-IT


Overview

FairRoot

FairRoot new features

Integrated Track follower (Geane)

Geant4 configuration classes

New reader (ROOT converted CAD Step Format )

Fast Simulation

CMake/CTest

FAIR experiments design studies

CBM

PANDA

Summary

Overview

IT-Palaver


Fairroot

FairRoot (Former CbmRoot) has started end of 2003

First released in March 2004

Oct 04 release was used to produce data for the CBM technical report

June 05 release ( Hades initialization scheme adapted )

Sept. 06 PANDA collaboration decided to use CbmRoot as simulation and analysis framework

Oct. 06 CbmRoot was renamed to FairRoot

FairRoot

IT-Palaver


Geane integration in fairroot

The integration into the VMC (TGeant3) is done

In FairRoot:

Geane can be used in the analysis or from macro

Propagation to

Length

Plane

Volume (Enter or Exit point)

CbmPoints and/or CbmTrackPar can be used as input for propagation

Geane Integration in FairRoot

IT-Palaver


What is geane 1

Package to calculate the average trajectories of particles through dense materials and to calculate the transport matrix as well as the propagated errors covariance matrix in a given track representation.

Geane is a tool to calculate extrapolated track parameters and propagated errors through dense materials.

With VMC it is straight forward to use it.

What is GEANE? (1)

IT-Palaver


What is geane 2

It is a track follower:

it predicts the trajectory of a charged particle in terms of mean values and errors both in forward and in backward direction. Three effects are taken into account:

energy loss (affects mean values and errors)

Coulomb multiple scattering (affects errors only)

magnetic field (affects mean values only)

Geometry and magnetic fields are handled by Geant3

(In VMC applications TGeoManager handles the geometry)

What is Geane? (2)

IT-Palaver


Geane muon absorber in cbm

Geane: Muon Absorber in CBM

IT-Palaver


Fair simulation analysis framework fairroot

Geane : Panda detector

IT-Palaver


Fair simulation analysis framework fairroot

Pulls for the whole Panda detector

s=1.41

s=1.13

s=1.19

s=1.05

s=0.96

IT-Palaver


Geane extrapolation hades mdc

Module 1 layer 1(D1S1) To Module 4 layer 6 (D4S6)

Propagate MC points (No input errors)

Geane::PropagateToVolume is used:

Helix track representation

Internal representation for Geane (can be transformed)

Pull distributions where not calculated :

Error are in the volume frame (have to be transformed to lab)

Geane Extrapolation (Hades MDC):

IT-Palaver


Geane hades example

Geane: Hades example

  • Hades Simulation

  • Geometry: 28 02 2003

  • Field map

  • Electrons:

    • 0.05-0.7 GeV

    • Polar angle range (20., 85.) degree

IT-Palaver


Hades simulation in fairroot

Hades Simulation in FairRoot

IT-Palaver


Geane red vs mc black

Geane (Red) vs. MC (Black)

IT-Palaver


Geane red vs mc black1

Geane (Red) vs. MC (Black)

IT-Palaver


New geometry reader

A standalone STEP to ROOT geometry converter has been implemented by Tobias Stockmanns (PANDA collaboration)

The output ROOT file has the geometry but not the full material properties needed by the simulation engines.

A reader has been introduced to the framework that read this geometry and replace the media definition by a proper one.

New Geometry Reader

IT-Palaver


Step to root geometry

STEP to ROOT Geometry

PANDA MVD detector

9345 volumes

IT-Palaver


Geant4 geant3 configuration

gconfig/g3Config.C

gconfig/g4Config.C

geant3->SetPAIR(1);

geant3->SetCOMP(1);

geant3->SetPHOT(1);

geant3->SetPFIS(0);

geant3->SetDRAY(1);

geant3->SetANNI(1);

geant3->SetBREM(1);

geant3->SetHADR(3);

geant3->SetMUNU(1);

geant3->SetDCAY(1);

geant3->SetLOSS(1);

geant3->SetMULS(1);

geant3->SetCKOV(1);

geant3->SetRAYL(1);

geant3->SetCUTS(…)

geant4->SetProcess("PAIR",1); /** pair production*/

geant4->SetProcess("COMP",1); /**Compton scattering*/

geant4->SetProcess("PHOT",1); /** photo electric effect */

geant4->SetProcess("PFIS",0); /**photofission*/

geant4->SetProcess("DRAY",1); /**delta-ray*/

geant4->SetProcess("ANNI",1); /**annihilation*/

geant4->SetProcess("BREM",1); /**bremsstrahlung*/

geant4->SetProcess("HADR",1); /**hadronic process*/

geant4->SetProcess("MUNU",1); /**muon nuclear interaction*/

geant4->SetProcess("DCAY",0); /**decay*/

geant4->SetProcess("LOSS",1); /**energy loss*/

geant4->SetProcess("MULS",1); /**multiple scattering*/

geant4->SetProcess("CKOV",1); /**Cerenkov photon generation*/

geant4->SetProcess("RAYL",1); /**Rayleigh scattering*/

...

geant4->SetCut("CUTGAM",cut1); /** gammas (GeV)*/ …

Geant4/Geant3 Configuration

IT-Palaver


Fair simulation analysis framework fairroot

Full Simulation-Analysis Chain

Determine particle properties

at target vertex

Event Generator

Transport particles through

the detector material

Simulation

Transport

SIM

Digitizer

Determine detector response

RAW

Determine physical space point

parameters from detector hits

Storage

Levels

Hit Finder

Determine momentum vector

and PID for all tracks

Analysis

Reconstruction

Physics Analysis

Calculate physics observables

IT-Palaver


Fair simulation analysis framework fairroot

Fast Simulation-Analysis Chain

Determine particle properties

at target vertex

Event Generator

No Transport

Put the events on the Stack

SIM

Parameterized detector

Response

Storage

Levels

Physics Analysis

Calculate physics observables

IT-Palaver


Fast simulation features

All generators available in FairRoot can be used

Events from different generators can be mixed

Detector response is implemented in Tasks

Give a unified output format for all different generators (CbmStack)

Fast analysis can be done with TTree::Draw()

The same code can be used to make fast simulation using the full simulation transport files (read only primary particles from stack)

The same physics analysis code can be used for fast and full simulation

Fast and full simulation results can be easily compared

Fast Simulation features

IT-Palaver


Fast simulation macro

// Load basic libraries

gROOT->LoadMacro("$VMCWORKDIR/gconfig/basiclibs.C");

basiclibs();

// Load this example libraries

gSystem->Load("libGeoBase");

gSystem->Load("libParBase");

gSystem->Load("libBase");

gSystem->Load("libMCStack");

gSystem->Load("libGen");

gSystem->Load("libPassive");

gSystem->Load("libPGen");

CbmRunSim *fRun = new CbmRunSim();

fRun->SetOutputFile("sim_fast.root");

Fast Simulation: Macro

IT-Palaver


Fast simulation macro1

// Create and Set Event Generator //

CbmPrimaryGenerator* primGen = new CbmPrimaryGenerator();

fRun->SetGenerator(primGen);

PndDpmGenerator *PndDpm = new

PndDpmGenerator("../../input/dpmevt_noelastic_36755.root");

primGen->AddGenerator(PndDpm);

/**switch off the transport of particles*/

primGen->DoTracking(kFALSE);

fRun->Init();

fRun->Run(1000);

Fast Simulation: Macro

IT-Palaver


Cmake ctest

Supports complex, large build environments.

CMake has been proved in large projects. (KDE 4)

Has powerful commands:

include the ability to locate include files, libraries, executables;

include external CMake files that encapsulate standard functionality; interfaces to testing systems;

supports recursive directory traversal with variable inheritance;

can run external programs;

supports conditional builds;

....

CMake & CTest

IT-Palaver


Cmake ctest1

Already in SVN (under testing)

Can run in parallel to Automake/Autoconf

Need to test:

Nightly Builds

E-mail to the user who has committed his code

Reports

See for an example of the web interface:

http://www.na-mic.org:8081/Insight/Dashboard/

http://lxg1417.gsi.de:8081

CMake & CTest

IT-Palaver


Cmake ctest dashboard

CMake/CTest : Dashboard

IT-Palaver


Dashboard cbm

Dashboard : CBM

IT-Palaver


Fair simulation analysis framework fairroot

Multi purpose detector at FAIR

Physics program

pp, pA collisions

1.515 GeV/c (p momentum)

  • Charmonium (cc) spectroscopy

  • Open charm spectroscopy

  • Search for gluonic excitations

    (hybrids - glueballs)

  • Charmed hadrons in nuclei

  • Single and double Hypernuclei

  • Other options (EFF, GPD, …)

The Panda experiment

IT-Palaver


Fair simulation analysis framework fairroot

PANDA Detector implementation: proposed geometry

IT-Palaver


Fair simulation analysis framework fairroot

Detector implementation: state of art

view from geometry manager

COILS

(dipole)

Muon Detector

DIRC

(Cherenkov)

TPC/STT

Micro Vertex

EMC

(Fwd EndCup)

COILS

(solenoid)

EMC

(barrel/Bkw EndCup)

IT-Palaver


Fair simulation analysis framework fairroot

crystal

crystal

cluster

two clusters

Reconstruction example: c in EMC barrel

c at rest

in lab frame

c

only barrel

not yet full

coverage

c

°

Clusterization

Jan Zhong

Dima Melnichuk

IT-Palaver


Fair simulation analysis framework fairroot

Task example - EMC reconstruction

Full reconstruction scheme

Migrated from Babar-like framework for PandaRoot

Dima Melnychuk (Warsaw)

IT-Palaver


Ongoing work

New Geane Interface is implemented and under testing

Check of Geant4 VMC interface (Physics list)

New G4 physics list for PANDA is needed

Energy cuts can be used directly, but the processes has to be adapted to the new list

Check of TFluka VMC interface

Use more intensively the TGeoManager services for the reconstruction:

First tests in PANDA MVD and CBM TOF has been made, and they look very promising

Ongoing work

IT-Palaver


Availability

The framework is available via SVN

https://subversion.gsi.de/fairroot/

Tested on

Red Hat 9.0 (gcc 3.2.2)

Suse 9.0 (gcc 3.3.1)

Suse 10.1 (gcc 4.1.0)

Debian (gcc 3.2.3)

Fedora Core 2(gcc 3.3.3)

Fedora Core 4 (gcc 4.0.0 )

Fedora Core 5 (gcc 4.1.0)

Gentoo (gcc 4.1.0)

SL 3 & 4

64 bit architectures (Debian, Suse)

Availability

IT-Palaver


  • Login