AnalysisTips < CALICE

CALICE Web>AnalysisTips (2011-09-19, AngelaLucaciTimoce) (raw view)

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
   * [[2011RunList: W-HCAL][2011: W-HCAL]]
   * [[2010RunList: W-HCAL][2010: W-HCAL]]
   * [[2009RunList][2009]]
   * [[2008RunList][2008]]
   * [[2007RunList][2007]]
   * [[2006RunList][2006]]
   * [[HowToUseElogInfo][How to use the e-log information]]

---++ Centrally produced runs
Follow [[CentrallyProducedRuns][this link]].

---++ Important analysis-related tasks
  This is the [[SoftwareAnalysisOrganisation][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 [[DataReprocessingRequest][here]].


---++ Getting Started - Basic Information
   * There are many ILC software packages at the [[http://ilcsoft.desy.de/portal/software_packages][ILCSoft portal]].
 Look at [[http://ilcsoft.desy.de/portal/software_packages/ilcinstal][lcinstall]] as a useful tool for installing them.
   * CALICE-specific ones can be found from the [[SoftwareMain][CALICE Software pages]]
   * AHCAL-specific ones can be found from the [[http://www-flc.desy.de/flc/flcwiki/HCAL_Documentation][flcWiki]] pages

---++ Sharing results/ideas
   * *Analysis Meetings* (approximately) bi-weekly, organised by telephone/web, details of the presentations can be found on  <a href="http://ilcagenda.linearcollider.org/categoryDisplay.py?categId=158">Indico</a>; or older meetings on  <a href="http://nicadd.niu.edu/cdsagenda//displayLevel.php?fid=30"> CDS</a> 

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

   * The [[MailingLists#CaliceSW][ *CALICE-SW mailing list* ]] is used for discussion of analysis/software: [[MailingLists#CaliceSW][how to join and view the archive of previous postings]].

---++ Example Analysis Job
<br />

---++ 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)  [[%ATTACHURL%/SimpleEcalDigi.cc][SimpleEcalDigi.cc]] ;  [[%ATTACHURL%/SimpleEcalDigi.h][SimpleEcalDigi.h]] ;  [[%ATTACHURL%/SimpleEcalDigi.xml][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)  [[%ATTACHURL%/ElectronSelection.cc][ElectronSelection.cc]] ;  [[%ATTACHURL%/ElectronSelection.h][ElectronSelection.h]] ;  [[%ATTACHURL%/ElectronSelection.xml][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)  [[%ATTACHURL%/HadronSelection.cc][HadronSelection.cc]] ;  [[%ATTACHURL%/HadronSelection.h][HadronSelection.h]] ;  [[%ATTACHURL%/HadronSelection.xml][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)  [[%ATTACHURL%/TriggerSim.cc][TriggerSim.cc]] ;  [[%ATTACHURL%/TriggerSim.h][TriggerSim.h]] ;  [[%ATTACHURL%/TriggerSim.xml][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)  [[%ATTACHURL%/TkProfile.cc][TkProfile.cc]] ;  [[%ATTACHURL%/TkProfile.h][TkProfile.h]] ;  [[%ATTACHURL%/TkProfile.xml][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.  
   * [[HitComposition][Notes on how to find out the composition of MC hits]] (DRW)
   * Mip Finder procesor (Philippe Doublet)  [[%ATTACHURL%/MipFinder.cc][MipFinder.cc]] ;  [[%ATTACHURL%/MipFinder.hh][MipFinder.h]] ;  [[%ATTACHURL%/MipFinder.xml][snippet of steering file]]  
      * Looks for MIPs in a given layer range (developed 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)  [[%ATTACHURL%/InteractionFinder.cc][InteractionFinder.cc]] ;  [[%ATTACHURL%/InteractionFinder.hh][InteractionFinder.h]] ;  [[%ATTACHURL%/InteractionFinder.xml][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.
   * Overlaying several CALICE events  (Vasiliy Morgunov)  [[%ATTACHURL%/Mixture.cc][Mixture.cc]] ; [[%ATTACHURL%/Mixture.h][Mixture.h]]  [[%ATTACHURL%/mix.xml][snippet of steering file]]
      *  Uses input files prepared by !CALICEOverlayPreparation.cc
      *  Allows shifting of events in X-Y plane with coordinates given in the steering file   
      *  Allows the primary track to be cut off

---++ 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: [[%ATTACHURL%/CreateSimHits.cc][CreateSimHits.cc]] and a small example input file: [[%ATTACHURL%/example.evt][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 [[http://www-flc.desy.de/flc/flcwiki/HcalRunLists][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").
      * Details of how to do this are at <a target="_self" href="http://grid.desy.de/certs/">http://grid.desy.de/certs/</a><br />
   * To check the status of a particular grid site in detail, you could look at  [[http://gstat.gridops.org/gstat][http://gstat.gridops.org/gstat]].

---++ 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.<br />
   * 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 [[http://ganga.web.cern.ch/ganga][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 [[MonteCarloSampleWishList]["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 [[MonteCarloProduction][centrally produced MC samples]].

---++ Mokka Tips
   * [[MokkaBeamProfile][How to get correct beam profile with Mokka/GEANT4]] (Angela) 

---++ CALICE Internal Notes
   The [[CaliceInternalNotes][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.

Attachments

Topic attachments
I	Attachment	History	Action	Size	Date	Who	Comment
cc	CreateSimHits.cc	r1	manage	7.2 K	2010-01-11 - 15:46	AngelaLucaciTimoceExternal	program to create an LCIO file with SimCalorimeterHits from an ASCII file
cc	ElectronSelection.cc	r1	manage	8.6 K	2009-04-30 - 12:09	DavidWard
h	ElectronSelection.h	r1	manage	1.6 K	2009-04-30 - 12:09	DavidWard
xml	ElectronSelection.xml	r4 r3 r2 r1	manage	1.1 K	2008-09-03 - 17:12	DavidWard
cc	HadronSelection.cc	r1	manage	7.9 K	2008-09-03 - 15:27	DavidWard
h	HadronSelection.h	r1	manage	1.6 K	2008-09-03 - 15:29	DavidWard
xml	HadronSelection.xml	r2 r1	manage	1.5 K	2008-09-03 - 17:14	DavidWard
cc	InteractionFinder.cc	r1	manage	5.7 K	2010-02-22 - 09:43	PhilippeDoubletExternal
hh	InteractionFinder.hh	r1	manage	1.9 K	2010-02-22 - 09:43	PhilippeDoubletExternal
xml	InteractionFinder.xml	r1	manage	1.2 K	2010-02-22 - 09:43	PhilippeDoubletExternal
cc	MipFinder.cc	r1	manage	19.1 K	2010-02-22 - 09:43	PhilippeDoubletExternal
hh	MipFinder.hh	r1	manage	2.6 K	2010-02-22 - 09:44	PhilippeDoubletExternal
xml	MipFinder.xml	r1	manage	2.3 K	2010-02-22 - 09:44	PhilippeDoubletExternal
cc	Mixture.cc	r1	manage	19.7 K	2010-12-02 - 16:00	DavidWard
h	Mixture.h	r1	manage	4.9 K	2010-12-02 - 16:02	DavidWard
cc	SimpleEcalDigi.cc	r1	manage	6.4 K	2008-06-09 - 13:25	DavidWard
h	SimpleEcalDigi.h	r1	manage	3.1 K	2008-06-09 - 13:26	DavidWard
xml	SimpleEcalDigi.xml	r2 r1	manage	1.0 K	2008-09-03 - 17:09	DavidWard
cc	TkProfile.cc	r1	manage	14.0 K	2009-01-30 - 11:59	DavidWard
h	TkProfile.h	r1	manage	5.2 K	2009-01-30 - 12:00	DavidWard
xml	TkProfile.xml	r1	manage	1.1 K	2009-01-30 - 12:02	DavidWard
cc	TriggerSim.cc	r1	manage	6.6 K	2008-09-03 - 15:28	DavidWard
h	TriggerSim.h	r1	manage	1.1 K	2008-09-03 - 15:29	DavidWard
xml	TriggerSim.xml	r1	manage	0.1 K	2008-09-03 - 17:10	DavidWard
evt	example.evt	r1	manage	1.8 K	2010-01-11 - 15:47	AngelaLucaciTimoceExternal	Small example ASCII file with 2 events
xml	mix.xml	r1	manage	3.9 K	2010-12-02 - 16:01	DavidWard

Topic revision: r28 - 2011-09-19 - AngelaLucaciTimoce

CAlorimeter for the LInear Collider Experiment

- Cern Search
- TWiki Search
- Google Search
CALICE All webs

Copyright &© 2008-2022 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