Linear Collider Forum



Home » Analysis and Reconstruction » Analysis Tools » Simple analysis example of stdhep files ?
Re: Simple analysis example of stdhep files ? [message #1268 is a reply to message #1265] Wed, 07 November 2007 03:11 Go to previous messageGo to previous message
gaede
Messages: 233
Registered: January 2004
Location: DESY, Hamburg
Hi Graham,

in case you want to use a C++ tool, Marlin can read stdhep files and convert them to LCIO. So you can analyze the 4vectors as MCParticles. See $MARLIN/examples/mymarlin/src/MyProcessor.cc.
You can easily build this example by copying the full example directory and then use cmake:
 cp -rp $MARLIN/examples/mymarlin .
 cd mymarlin
 mkdir build
 cd build 
 cmake -C ../BuildSetup.cmake ..
 make install
 cd ..
 export MARLIN_DLL=$PWD/lib/libmymarlin.so
 #run
 $MARLIN/bin/Marlin readstdhep.xml


You have to modify the BuildSetup.cmake to sth. like:
# ILC_HOME
SET( ILC_HOME "/afs/desy.de/group/it/ilcsoft/v01-01"
    CACHE PATH "Path to ILC Software" FORCE )

# Path to Marlin
SET( Marlin_HOME "${ILC_HOME}/Marlin/v00-09-08"
    CACHE PATH "Path to Marlin" FORCE )

# Path to LCIO
SET( LCIO_HOME "${ILC_HOME}/lcio/v01-08-03"
    CACHE PATH "Path to LCIO" FORCE )

## Path to RAIDA
SET( RAIDA_HOME "${ILC_HOME}/RAIDA/v01-03"
    CACHE PATH "Path to RAIDA" FORCE )

# CMake Modules Path
SET( CMAKE_MODULE_PATH "${ILC_HOME}/CMakeModules/v01-01"
    CACHE PATH "Path to CMake Modules" FORCE )


SET( BUILD_WITH "RAIDA"
    CACHE STRING "Build mymarlin with these optional packages" FORCE )




I have attached an example steering file that shows how to read sthep files.

If this example doesn't work your stdhep file might be corrupted (or incompatible with the LCIO code). Could you then point us to the file that you are using.

Frank.

 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic:Neural Nets and lcsim.org
Next Topic:JAS3 support of HBOOK and ROOT files
Goto Forum:
  


Current Time: Tue Oct 16 13:21:23 Pacific Daylight Time 2018
.:: Contact :: Home ::.

Powered by: FUDforum 3.0.1.
Copyright ©2001-2010 FUDforum Bulletin Board Software