Commit in hps-java/src/main on MAIN | |||
resources/org/lcsim/hps/steering/MuonAnalysis.lcsim | +21 | added 1.1 | |
java/org/lcsim/hps/analysis/muon/SimpleMuonAnalysis.java | +51 | added 1.1 | |
+72 |
skeleton analysis for new muon detector implementation
diff -N MuonAnalysis.lcsim --- /dev/null 1 Jan 1970 00:00:00 -0000 +++ MuonAnalysis.lcsim 28 Jan 2013 19:20:17 -0000 1.1 @@ -0,0 +1,21 @@
+<?xml version="1.0" encoding="UTF-8"?> +<lcsim xmlns:xs="http://www.w3.org/2001/XMLSchema-instance" xs:noNamespaceSchemaLocation="http://www.lcsim.org/schemas/lcsim/1.0/lcsim.xsd"> + <inputFiles> + <file>${inputFile}</file> + </inputFiles> + <execute> + <driver name="EventMarkerDriver" /> + <driver name="MuonAnalysis" /> + <driver name="AidaSaveDriver" /> + </execute> + <drivers> + <driver name="EventMarkerDriver" type="org.lcsim.job.EventMarkerDriver"> + <eventInterval>1</eventInterval> + </driver> + <driver name="MuonAnalysis" type="org.lcsim.hps.analysis.muon.SimpleMuonAnalysis"> + </driver> + <driver name="AidaSaveDriver" type="org.lcsim.job.AidaSaveDriver"> + <outputFileName>muon_plots</outputFileName> + </driver> + </drivers> +</lcsim>
diff -N SimpleMuonAnalysis.java --- /dev/null 1 Jan 1970 00:00:00 -0000 +++ SimpleMuonAnalysis.java 28 Jan 2013 19:20:18 -0000 1.1 @@ -0,0 +1,51 @@
+package org.lcsim.hps.analysis.muon; + +import hep.aida.ICloud1D; +import hep.aida.ICloud2D; +import hep.aida.IHistogram1D; + +import java.util.HashMap; +import java.util.List; + +import org.lcsim.event.CalorimeterHit; +import org.lcsim.event.EventHeader; +import org.lcsim.geometry.Detector; +import org.lcsim.units.clhep.SystemOfUnits; +import org.lcsim.util.Driver; +import org.lcsim.util.aida.AIDA; + +/** + * @author Jeremy McCormick <[log in to unmask]> + * @version $Id: SimpleMuonAnalysis.java,v 1.1 2013/01/28 19:20:18 jeremy Exp $ + */ +public class SimpleMuonAnalysis extends Driver +{ + AIDA aida = AIDA.defaultInstance(); + ICloud1D hite = aida.cloud1D("Hit Energy [MeV]"); + ICloud1D nhits = aida.cloud1D("Number Of Hits Per Event"); + ICloud2D xy = aida.cloud2D("Hit Position XY"); + IHistogram1D layer = aida.histogram1D("Hit Layer Number", 4, 0, 4); + + public void startOfData() + { + } + + public void endOfData() + { + } + + public void detectorChanged(Detector det) + { + } + + public void process(EventHeader event) + { + List<CalorimeterHit> hits = event.get(CalorimeterHit.class, "MUON_HITS"); + nhits.fill(hits.size()); + for (CalorimeterHit hit : hits) { + xy.fill(hit.getPosition()[0], hit.getPosition()[1]); + hite.fill(hit.getRawEnergy() * SystemOfUnits.MeV); + layer.fill(hit.getLayerNumber()); + } + } +}
\ No newline at end of file
Use REPLY-ALL to reply to list
To unsubscribe from the LCD-CVS list, click the following link:
https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=LCD-CVS&A=1