TWiki> CALICE Web>AnalysisTips (revision 24)EditAttachPDF
This page is to allow users (as well as core software developers) to exchange hints/tips which they think others analysing data might find useful. Developers are welcome to add to these pages, but are encouraged not to change the content from analysers without checking with the contributer first.

Run lists

Important analysis-related tasks

This is the evolving list of critical jobs which need to be carried out, each related to publishing a particular paper. Please look at the list to see where you can help!

Runs for Nov. 2009 reprocessing tests

These are a small number of ~20-30 "typical" runs which allow technical checks that all steps necessary for reprocessing can be performed by the people who will be responsible, and that the analysers can carry out rapid quality checks of the reprocessed data. The list of such runs (required per subsystem and per year) is maintained here.

Getting Started - Basic Information

Sharing results/ideas

  • Analysis Meetings (approximately) bi-weekly, organised by telephone/web, details of the presentations can be found on Indico; or older meetings on CDS

This is the forum for discussing analysis results and software developments.

Example Analysis Job

Useful Marlin processors for analysis

These may be useful as is, or may be useful starting points for your own work.
  • Naive ECAL digitisation procesor (DRW) ; SimpleEcalDigi.h ; snippet of steering file
    • Average noise (default 0.12 MIP) applied to all cells containing signal.
    • No noise contributions in empty cells.
    • No simulation of dead cells.
    • Can disable layers at the back of the calorimeter (for DESY'06) or lower slabs (DESY/CERN'06).
  • Electron selection procesor (DRW) ; ElectronSelection.h ; snippet of steering file
    • Uses ECAL information only.
    • Loose cut on visible energy.
    • Optional cuts on Cerenkov trigger (data only) and double shower removal.
    • Only essential input parameter is beam energy (not reliably included in event data yet).
    • Selected event flagged in bit 1 of EventParameter "EventClass"
  • Hadron selection procesor (DRW) ; HadronSelection.h ; snippet of steering file
    • Simple cuts on ECAL, HCAL, TCMT energy to identify muons (flagged in bit 3 of EventParameter "EventClass")
    • Cut on total ECAL+HCAL energy against double beam particles
    • Cut on energy in first two layers to reduce upstream showering
    • Optionally demand or veto on Cerenkov trigger (data only)
    • Selected event flagged in bit 2 of EventParameter "EventClass"
  • Simple trigger procesor (DRW) ; TriggerSim.h ; snippet of steering file
    • Sets trigger bits for MC events depending on hits in beam scintillators.
    • Only set up and tested for 2006 (CERN) and 2007 (CERN).
  • Comparison of track and Calorimeter positions (DRW) ; TkProfile.h ; snippet of steering file
    • Fills profile histograms of shower positions in ECAL and HCAL and compares with extrapolated track.
    • Prints estimate of calorimeter displacements and of drift velocity corrections.
    • Also creates a small root file which may be checked.
  • Notes on how to find out the composition of MC hits (DRW)
  • Mip Finder procesor (Philippe Doublet) ; MipFinder.h ; snippet of steering file
    • Looks for MIPs in a given layer range (developped for the Si-W ECAL).
    • Cuts are distance and angle between hits and/or straight clusters.
    • Ouput is a cluster of Calorimeter Hits.
  • Interaction Finder procesor (Philippe Doublet) ; InteractionFinder.h ; snippet of steering file
    • Finds the first interaction layer (developped for the Si-W ECAL).
    • Can distinguish events with peaked layers and Fireball events.
    • Parameters are written on the Calorimeter Hit collection.

Useful executables (based on Marlin) for analysis

* ASCII to LCIO file: small program used to convert an ASCII file (e.g. output from GEANT3) to an LCIO file; uses some assumptions on the content of the ASCII file. Here is the program: and a small example input file: example.evt

Useful runs

A list of 2007 and 2008 runs that seem to useful for HCAL and / or combined analysis is listed in the flcWiki. It contains pion and proton data taken at CERN and Fermilab with the Si-W ECAL, AHCAL and TCMT, as well as HCAL only runs.

ConfigurationChanges for 2007 data

Access to Data

  • You need to obtain a grid certificate - for details, contact your local grid Certificate Authority)
  • When you have a grid certificate, you then need to request membership of the CALICE VO ("virtual organisation").
  • To check the status of a particular grid site in detail, you could look at

Running Grid Jobs

  • Your Marlin job needs to be bundled together for submission to the grid. Assume you will make a binary as on your local PC, upload this to the grid, and then submit your job with the XML steering file each time you want to run over a large number of files.
  • Your shell script to run Marlin will need to download the files to be analysed to the grid machine ("Computing Element") where it happens to start running.
  • Grid jobs can be submitted using a variety of tools, one convenient one is Ganga.

Monte Carlo Samples

If you need specific Monte Carlo samples generated for a particular study and think these would be useful for others, please add them to the "wish list" and email the responsible people to discuss your needs or propose it in the next analysis/software meeting.

Further information is available on the pages for centrally produced MC samples.

Mokka Tips

CALICE Internal Notes

The CALICE Internal Notes describe technical details of analyses which have been performed. They are explicitly for use only by members of the Collaboration and their results are not approved for presentation as results of the Collaboration.
Topic attachments
I Attachment History Action Size Date Who Comment
Unknown file formatcc r1 manage 7.2 K 2010-01-11 - 15:46 AngelaLucaciTimoceExternal program to create an LCIO file with SimCalorimeterHits from an ASCII file
Unknown file formatcc r1 manage 8.6 K 2009-04-30 - 12:09 DavidWard  
Header fileh ElectronSelection.h r1 manage 1.6 K 2009-04-30 - 12:09 DavidWard  
XMLxml ElectronSelection.xml r4 r3 r2 r1 manage 1.1 K 2008-09-03 - 17:12 DavidWard  
Unknown file formatcc r1 manage 7.9 K 2008-09-03 - 15:27 DavidWard  
Header fileh HadronSelection.h r1 manage 1.6 K 2008-09-03 - 15:29 DavidWard  
XMLxml HadronSelection.xml r2 r1 manage 1.5 K 2008-09-03 - 17:14 DavidWard  
Unknown file formatcc r1 manage 5.7 K 2010-02-22 - 09:43 PhilippeDoubletExternal  
Unknown file formathh InteractionFinder.hh r1 manage 1.9 K 2010-02-22 - 09:43 PhilippeDoubletExternal  
XMLxml InteractionFinder.xml r1 manage 1.2 K 2010-02-22 - 09:43 PhilippeDoubletExternal  
Unknown file formatcc r1 manage 19.1 K 2010-02-22 - 09:43 PhilippeDoubletExternal  
Unknown file formathh MipFinder.hh r1 manage 2.6 K 2010-02-22 - 09:44 PhilippeDoubletExternal  
XMLxml MipFinder.xml r1 manage 2.3 K 2010-02-22 - 09:44 PhilippeDoubletExternal  
Unknown file formatcc r1 manage 6.4 K 2008-06-09 - 13:25 DavidWard  
Header fileh SimpleEcalDigi.h r1 manage 3.1 K 2008-06-09 - 13:26 DavidWard  
XMLxml SimpleEcalDigi.xml r2 r1 manage 1.0 K 2008-09-03 - 17:09 DavidWard  
Unknown file formatcc r1 manage 14.0 K 2009-01-30 - 11:59 DavidWard  
Header fileh TkProfile.h r1 manage 5.2 K 2009-01-30 - 12:00 DavidWard  
XMLxml TkProfile.xml r1 manage 1.1 K 2009-01-30 - 12:02 DavidWard  
Unknown file formatcc r1 manage 6.6 K 2008-09-03 - 15:28 DavidWard  
Header fileh TriggerSim.h r1 manage 1.1 K 2008-09-03 - 15:29 DavidWard  
XMLxml TriggerSim.xml r1 manage 0.1 K 2008-09-03 - 17:10 DavidWard  
Unknown file formatevt example.evt r1 manage 1.8 K 2010-01-11 - 15:47 AngelaLucaciTimoceExternal Small example ASCII file with 2 events
Edit | Attach | Watch | Print version | History: r28 | r26 < r25 < r24 < r23 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r24 - 2010-06-08 - AngelaLucaciTimoceExternal
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    CALICE All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2023 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback