<parameter name="LCIOInputFiles"> simjob.slcio </parameter> ..
<parameter name="GearXMLFile"> gear_ldc.xml </parameter> ..
</global>

Marlin -c mysteer.xml

LCIO Available Collections:
ECAL007                                 SimCalorimeterHit             simjob.slcio
FinalMCParticles                        MCParticle                    simjob.slcio
MCParticle                              MCParticle                    simjob.slcio
SiliconRawHitExample                    TrackerPulse                  simjob.slcio
TPC4711                                 SimTrackerHit                 simjob.slcio
TPC4711UserFloatExtension               LCFloatVec                    simjob.slcio
TPC4711UserIntExtension                 LCIntVec                      simjob.slcio
TPCRawFADCMCTruth                       LCRelation                    simjob.slcio
TrackerDataExample                      TrackerData                   simjob.slcio
TrackerPulseExample                     TrackerPulse                  simjob.slcio
TrackerRawDataExample                   TrackerRawData                simjob.slcio
AdditionalExtension                     LCFloatVec                    simjob.slcio

Active Processors:
MyAIDAProcessor                         AIDAProcessor                  [ Active ] 
MyTestProcessor                         TestProcessor                  [ Active ] 
MyLCIOOutputProcessor                   LCIOOutputProcessor            [ Active ] 

Inactive Processors:
MyAbsCalibr                             AbsCalibr                      [ Inactive ] 
MyCEDViewer                             CEDViewer                      [ Inactive ] 
MyCalorimeterFinalStage                 CalorimeterFinalStage          [ Inactive ] 
MyClusterCheater                        ClusterCheater                 [ Inactive ] 
MyClusterCheater5_3                     ClusterCheater5_3              [ Inactive ] 
MyConditionsProcessor                   ConditionsProcessor            [ Inactive ] 
MyConversionTagger                      ConversionTagger               [ Inactive ] 
MyConvertToMIPScale                     ConvertToMIPScale              [ Inactive ] 
MyCurlKillerProcessor                   CurlKillerProcessor            [ Inactive ] 
MyCutOnCosThetaQQbar                    CutOnCosThetaQQbar             [ Inactive ] 
MyCutOnGEANT4Bug                        CutOnGEANT4Bug                 [ Inactive ] 
MyDSTAIDAPlotProcessor                  DSTAIDAPlotProcessor           [ Inactive ] 
MyDSTCollectionProcessor                DSTCollectionProcessor         [ Inactive ] 
MyDSTPlotProcessor                      DSTPlotProcessor               [ Inactive ]
MyDrawMCParticles                       DrawMCParticles                [ Inactive ]
MyEMShowerFinder                        EMShowerFinder                 [ Inactive ]
MyETDDigiProcessor                      ETDDigiProcessor               [ Inactive ]
....

No Errors were found :)

Esto es debido a que la variable de entorno MARLIN_DLL carga las librerias dinámicas que le digas, pero eso no quiere decir que ejecute nada de ellas. Para ejecutar algún módulo que esté contenido en alguna librería hay que activarlo en la sección "execute" del fichero steer de Marlin. Por tanto, que un procesador esté inactivo, no es ningún error, está inactivo porque esta cargado en memoria pero no se hace uso de él. Por ejemplo, supón que cargas la librería que contiene el procesador del visor:

 export MARLIN_DLL=/mnt/data2/ilcsoft/ilcinstall_prueba/ilcinstall/my_v01-04/CEDViewer/v00-06/lib/libCEDViewer.so

<execute>
  <processor name="MyAIDAProcessor"/>
  <processor name="MyTestProcessor"/>  
  <processor name="MyLCIOOutputProcessor"/>  
 </execute>

<execute>
  <processor name="MyAIDAProcessor"/>
  <processor name="MyTestProcessor"/>  
  <processor name="MyCEDViewer" />
  <processor name="MyLCIOOutputProcessor"/>  
 </execute>

Ahora, si lanzo Marlin, veré el procesador activo:

Active Processors:
MyAIDAProcessor                         AIDAProcessor                  [ Active ]
MyTestProcessor                         TestProcessor                  [ Active ]
MyCEDViewer                             CEDViewer                      [ Active ]
MyLCIOOutputProcessor                   LCIOOutputProcessor            [ Active ]

2.a.ii. Dumpevent

La herramienta dumpevent muestra el contenido de un fichero en formato lcio, nos servimos de ella para comprobar que los ficheros que nos hemos descargado sólo contienen hasta la simulación del paso de las partículas a través del detector ILD.

En mi instalación, esta herramienta se encuentra en el subdirectorio: /ilcinstall/my_v01-04/lcio/v01-10-01/bin:

[iglesias@fpmac116 bin]$ dumpevent
 usage: dumpevent filename runNum evtNum
    or: dumpevent filename n
  where the first dumps the event with the specified run and event number
  and the second simply dumps the n-th event in the file

[iglesias@fpmac116]$ ./dumpevent ../Samples/SinglePar_muon_Theta1-40_1-50GeV_0001.slcio 2

SIOReader::getEventMap() recreating event map for direct access ...
 SIOReader::getEventMap() : done

============================================================================
        Event  : 1 - run:  0 - timestamp 0 - weight 1

============================================================================
 date:      01.01.1970  00:00:00.000000000
 detector : ILD_00
 event parameters:

collection name : ETDCollection
 parameters:

--------------- print out of SimTrackerHit collection ---------------

  flag:  0xb00
     LCIO::THBIT_BARREL : 0
     LCIO::THBIT_MOMENTUM : 0

 cellID[indices] | position (x,y,z) | dEdx | time  | PDG of MCParticle

    -----------------------------------------
       WARNING: CellIDDecoder - no CellIDEncoding parameter in collection !         -> using default : "byte0:8,byte1:8,byte2:8,byte3:8"
    ------------------------------------------
0: 1[byte0:1,byte1:0,byte2:0,byte3:0] | (370.141, 790.693, 2426) | 9.19724e-05 | 8.60581 | 13
1: 2[byte0:2,byte1:0,byte2:0,byte3:0] | (370.338, 791.39, 2428) | 8.59193e-05 |8.61291 | 13
2: 3[byte0:3,byte1:0,byte2:0,byte3:0] | (370.533, 792.088, 2430) | 7.71785e-05 | 8.62 | 13
3: 1[byte0:1,byte1:0,byte2:0,byte3:0] | (360.89, 759.687, 2426) | 0.00025424 | 8.78791 | 11

--------------------------------------------------------------------------------

• La colección MCParticle (que siempre existirá una única por suceso) contiene las partículas del Montecarlo truth que fueron generadas por el generador o que se añadieron (decay in flight) por el programa de simulación, Mokka.
• La clase genérica (y por tanto virtual) SimCalorimeterHit almacena colecciones de impactos en los calorímetros a partir del programa de simulación. Cada subdetector creará una colección de este tipo, por ejemplo, el ECAL de las tapas nos creará una colección de impactos en sus calorímetros llamada EcalEndcapCollection.
• La clase genérica SimTrackerHit almacena impactos de los subdetectores de trayectorias, por ejemplo el subdetector TPC, creará la colección de trazas simuladas TPCCollection.

Para este caso concreto del evento 2 de la muestra SinglePar_muon_Theta1-40_1-50GeV_0001.slcio , me aparece las siguientes colecciones para las clases SimTrackerHit, SimCalorimeterHit y MCParticle:

SimCalorimeterHit:

• EcalEndcapCollection
• EcalEndcapPreShowerCollection
• HcalEndCapsCollection
• MuonEndCapCollection

MCParticle:

• MCParticle

SimTrackerHit:

• ETDCollection
• TPCCollection
• VXDCollection

Pero si miramos en el resto de sucesos, encontraremos el resto de las colecciones para los diferentes detectores, por ejemplo SETCollection, SITCollection ... LumiCalCollection ... y otras. La informacion respecto al tipo de colecciones que corresponden a un suceso determinado puede verse en bash con:

dumpevent NombreFichero.slcio n|grep 'collection name'

dumpevent ../Samples/SinglePar_muon_Theta1-40_1-50GeV_0001.slcio 2|grep 'collection name'

 collection name : ETDCollection
 collection name : EcalEndcapCollection
 collection name : EcalEndcapPreShowerCollection
 collection name : HcalEndCapsCollection
 collection name : MCParticle
 collection name : MuonEndCapCollection
 collection name : SITCollection
 collection name : TPCCollection
 collection name : VXDCollection

for i in `seq 1 1000`; do dumpevent NombreFichero.slcio $i|grep 'collection name'; done

for i in `seq 1 1000`; do dumpevent ../Samples/SinglePar_muon_Theta1-40_1-50GeV_0001.slcio $i|grep 'collection name'; done

for i in `seq 1 1000`; do echo "$i :";dumpevent ../Samples/SinglePar_muon_Theta1-40_1-50GeV_0001.slcio $i |grep 'collection name' ; done

37 :
 collection name : ETDCollection
 collection name : EcalEndcapCollection
 collection name : EcalEndcapPreShowerCollection
 collection name : FTDCollection
 collection name : HcalEndCapsCollection
 collection name : MCParticle
 collection name : MuonEndCapCollection
 collection name : TPCCollection
 collection name : VXDCollection
38 :
 collection name : LHcalCollection
 collection name : LumiCalCollection
 collection name : MCParticle
 collection name : MuonEndCapCollection
39 :
 collection name : ETDCollection
 collection name : EcalEndcapCollection
 collection name : EcalEndcapPreShowerCollection
 collection name : HcalEndCapsCollection
 collection name : MCParticle
 collection name : MuonEndCapCollection
 collection name : TPCCollection
 collection name : VXDCollection
40 :
 collection name : ETDCollection
 collection name : EcalEndcapCollection
 collection name : EcalEndcapPreShowerCollection
 collection name : HcalEndCapsCollection
 collection name : MCParticle
 collection name : MuonEndCapCollection
 collection name : TPCCollection
 collection name : TPCSpacePointCollection

2.a.iii. Visor de Eventos

CED es el visor de sucesos dentro del entorno Marlin. Muestra una ventanita con lo que tengas descrito del detector y trazas, hits, etc... resultado de la reconstruccion. CED usa OpenGL y la herramienta GLUT (OpenGL Utility Toolkit), por lo que para correr CED GLUT previamente lo debes tener intalada en tu sistema.

Si ejecutamos Marlin -x, nos devolverá un ejemplo de fichero steer con todos los procesadores que tenemos cargados (tanto estática como dinámicamente); cada vez que carguemos un nuevo procesador, via la variable de entorno MARLIN_DLL, nos aparecerá en el fichero steer de ejemplo. Así, si lo ejecutamos después de incluir la librería del visor, veremos que disponemos de unos cuantos procesadores relacionados con el visor de sucesos: GenericViewer, CEDViewer, VertexViewer, DrawMCParticles, ... . Centrándonos en CEDViewer, activamos el procesador en la etiqueta execute del fichero steer (ver apartado sobre Marlin), y elegimos que colecciones vamos a visualizar (para extraer las colecciones de cada suceso ver apartado sobre Dumpevent). Para especificar que colecciones quiero visualiza voy alfichero steer, y bajo la definición del procesador CEDViewer:

<processor name="MyCEDViewer" type="CEDViewer...>
....
</processor>

<parameter name="DrawCollection" type="StringVec">TPCTracks 0 2</parameter>

Para utilizar el visor de sucesos, antes de ejecutar Marlin, lanzas

glced &

Marlin mysteer.xml

Ahora el suceso es mostrado en la ventana CED despues de finalizar la reconstruccion de trazas, clusters y objetos particle flow. Para manipular el evento en pantalla, se puede usar el raton (siempre que este colocado dentro de la ventana CED) o ciertos comnados del teclado:

• Usando el raton:
  • boton izquierdo de raton: rotar
  • boton central del raton: trasladas
  • boton derecho del raton: aumentar
• Usando los comandos del teclado:
  • ESC : front view. Projection of all points on the -plane.
  • s: side view. Projection of all points on the -plane.
  • c: centre view. A new origin for a rotation and zoom is set to the current position of the mouse cursor.
  • r: redisplay in some initial projection.
  • 0, 1, ... toggle layers. The assignment of the layers is done in the steering file

After an event is displayed, the system waits until the user hits a key to continue. Please note, that since the event display is realised as a Marlin processor, the reconstruction is also stopped as long as the computer waits for the user entry to continue.

-- Main.iglesias - 11 Dec 2008
to top