---------- Forwarded message ----------
From: Sebouh Paul <[log in to unmask]>
Date: Tue, Oct 14, 2014 at 8:54 PM
Subject: Re: no tracks found
To: "McCormick, Jeremy I." <[log in to unmask]>
I typed:
ifarm1101> java -Xmx1500m -jar
~/hps-distribution-3.0.3-20141007.165651-110-bin.jar ~/SimpleMCRecon.lcsim
-i /work/hallb/hps/sebouh/trident_slic/radmuon_12.lcio.slcio
-DoutputFile=/work/hallb/hps/sebouh/test -DrunNumber=100
1) I'm running montecarlo muon radiative and bethe-heitler trident events
for 4.4 GeV beam energy, so that I can figure out how far out the threshold
in z should be set for the muon vertexing
2) Stepan set me up to this.
3) I couldnt find existing events for this
4, 5) I have modified nothing: my command is
/group/hps/hps_soft/slic/v00-02/slic/v03-01-04/build/bin/slic -P
/group/hps/hps_soft/slic/v00-02/slic/v03-01-04/data/particle
.tbl -g
~/detector-data/detectors/HPS-Proposal2014-v7-6pt6/HPS-Proposal2014-v7-6pt6.lcdd
-i
trident/radmuon_${N}.stdhep -o trident_slic/radmuon_${N}.lcio -r 10000
5) I set EventMarkerDriver eventInterval to 1, and found that it was event
#1778 that gets the error:
....
Event 1788
getPathLengthToPlaneApprox ERROR t is negative (-0.02849407469347215)!
p 0.5016171795755748 rho 8.964778426857066E-4 a 4.4968868700000003E-4 A
3.226783361630766E-4 B 0.6932096975604128 C 394.38172820812906
track params: 0.33654476569411973, 0.3715470696473806,
-0.017612227732148794, 0.5016171795755748, 1289.7084037010472,
313.646930637243, -48.658595699839786,
xp [ 902.35, 77.493, -30.330]
eta [ 0.99953, 0.030515, 0.0000]
h [ 0.0000, 0.0000, 1.0000]
java.lang.RuntimeException: Problem in calculating the approximate path
length to the plane.
at org.hps.recon.tracking.WTrack.getPathLengthToPlaneApprox(WTrack.java:163)
at org.hps.recon.tracking.WTrack.getHelixAndPlaneIntercept(WTrack.java:288)
at
org.hps.recon.tracking.TrackUtils.getHelixPlaneIntercept(TrackUtils.java:141)
at
org.hps.recon.tracking.MultipleScattering.getHelixIntersection(MultipleScattering.java:288)
at
org.hps.recon.tracking.MultipleScattering.getHelixIntersection(MultipleScattering.java:157)
at
org.hps.recon.tracking.MultipleScattering.FindHPSScatterPoints(MultipleScattering.java:96)
at
org.hps.recon.tracking.MultipleScattering.FindHPSScatters(MultipleScattering.java:59)
at
org.hps.recon.tracking.MultipleScattering.FindScatters(MultipleScattering.java:44)
at
org.lcsim.recon.tracking.seedtracker.HelixFitter.FitCandidate(HelixFitter.java:110)
at
org.lcsim.recon.tracking.seedtracker.SeedTrackFinder.FindTracks(SeedTrackFinder.java:140)
at
org.lcsim.recon.tracking.seedtracker.SeedTracker.process(SeedTracker.java:111)
at org.lcsim.util.Driver.doProcess(Driver.java:273)
at org.lcsim.util.Driver.processChildren(Driver.java:284)
at org.lcsim.util.Driver.process(Driver.java:198)
at
org.hps.recon.tracking.TrackerReconDriver.process(TrackerReconDriver.java:172)
at org.lcsim.util.Driver.doProcess(Driver.java:273)
at org.lcsim.util.Driver.processChildren(Driver.java:284)
at org.lcsim.util.Driver.process(Driver.java:198)
at org.lcsim.util.DriverAdapter.recordSupplied(DriverAdapter.java:74)
at
org.freehep.record.loop.DefaultRecordLoop.consumeRecord(DefaultRecordLoop.java:832)
at
org.freehep.record.loop.DefaultRecordLoop.loop(DefaultRecordLoop.java:668)
at
org.freehep.record.loop.DefaultRecordLoop.execute(DefaultRecordLoop.java:566)
at org.lcsim.util.loop.LCSimLoop.loop(LCSimLoop.java:151)
at org.lcsim.job.JobControlManager.run(JobControlManager.java:418)
at org.lcsim.job.JobControlManager.main(JobControlManager.java:180)
On Tue, Oct 14, 2014 at 8:36 PM, McCormick, Jeremy I. <
[log in to unmask]> wrote:
> Hi, there.
>
> Can we backup a minute here so I can get some info from you....
>
> What are you trying to accomplish overall with this?
>
> Who requested you to do this task?
>
> Are there already existing simulated events at JLab that you can use for
> this?
>
> What detector model exactly are you using? Have you modified at all any
> of the detector description files?
>
> What version of SLIC are you using?
>
> I'm not sure about that error. I've actually never seen it before...
>
> Can you set "eventInterval" in EventMarkerDriver to "1" so that we can
> exactly which event is causing the problem?
>
> Thanks.
>
> --Jeremy
>
> -----Original Message-----
> From: Sebouh Paul [mailto:[log in to unmask]]
> Sent: Tuesday, October 14, 2014 5:32 PM
> To: McCormick, Jeremy I.
> Subject: Re: no tracks found
>
> actually, it has now reached the same type of error as before, but prints
> out other stuff before the error appears:
>
> GBLOutput: perPar covariance matrix
> [[ 0.045286, -0.00037628, -1.0901e-06, 0.0000, 0.0000]
> [ -0.00037628, 3.6293e-06, 1.1360e-08, 0.0000, 0.0000]
> [ -1.0901e-06, 1.1360e-08, 4.0390e-11, 0.0000, 0.0000]
> [ 0.0000, 0.0000, 0.0000, 0.0016188, -1.6808e-05]
> [ 0.0000, 0.0000, 0.0000, -1.6808e-05, 1.7911e-07]
> GBLOutput: truth perPar chi2
> 58110239054383830000000000000000000000000.000000
> 5 hits
> GBLOutput: layer 1
> trkpos at intercept [88.7457316196 2.2479630992 1.5571328860] trkposTruth
> at intercept [ -1.0000e+06, -1.0000e+06, -1.0000e+06] layer 1 uRes
> -0.0002020983
> GBLOutput: layer 2
> trkpos at intercept [96.1937023167 2.5592111482 1.6849016328] trkposTruth
> at intercept [ -1.0000e+06, -1.0000e+06, -1.0000e+06] layer 2 uRes
> 0.0013486757
> GBLOutput: layer 3
> trkpos at intercept [188.6087034940 8.2624687543 3.2719860485] trkposTruth
> at intercept [ -1.0000e+06, -1.0000e+06, -1.0000e+06] layer 3 uRes
> 0.0266609766
> GBLOutput: layer 4
> trkpos at intercept [196.0476245141 8.8701410158 3.3999129666] trkposTruth
> at intercept [ -1.0000e+06, -1.0000e+06, -1.0000e+06] layer 4 uRes
> -0.0294352594
> GBLOutput: layer 5
> trkpos at intercept [288.3499734145 18.2633493836 4.9902245570]
> trkposTruth at intercept [ -1.0000e+06, -1.0000e+06, -1.0000e+06] layer 5
> uRes 0.0195566949
> GBLOutput: layer 6
> trkpos at intercept [295.7797894214 19.1692584111 5.1185138278]
> trkposTruth at intercept [ -1.0000e+06, -1.0000e+06, -1.0000e+06] layer 6
> uRes -0.0179816191
> GBLOutput: layer 7
> trkpos at intercept [487.9857199620 50.4835337402 8.4571671653]
> trkposTruth at intercept [ -1.0000e+06, -1.0000e+06, -1.0000e+06] layer 7
> uRes -0.0938543443
> GBLOutput: layer 8
> trkpos at intercept [495.3950348020 51.9985052331 8.5867894143]
> trkposTruth at intercept [ -1.0000e+06, -1.0000e+06, -1.0000e+06] layer 8
> uRes 0.1020305513
> GBLOutput: layer 9
> trkpos at intercept [686.5898944397 99.2563507857 11.9633261476]
> trkposTruth at intercept [ -1.0000e+06, -1.0000e+06, -1.0000e+06] layer 9
> uRes -0.2937698932
> GBLOutput: layer 10
> trkpos at intercept [693.9798596028 101.4051253157 12.0952351167]
> trkposTruth at intercept [ -1.0000e+06, -1.0000e+06, -1.0000e+06] layer 10
> uRes 0.2852378959 getPathLengthToPlaneApprox ERROR t is negative
> (-0.02849407469347215)!
> p 0.5016171795755748 rho 8.964778426857066E-4 a 4.4968868700000003E-4 A
> 3.226783361630766E-4 B 0.6932096975604128 C 394.38172820812906 track
> params: 0.33654476569411973, 0.3715470696473806, -0.017612227732148794,
> 0.5016171795755748, 1289.7084037010472, 313.646930637243,
> -48.658595699839786,
> xp [ 902.35, 77.493, -30.330]
> eta [ 0.99953, 0.030515, 0.0000]
> h [ 0.0000, 0.0000, 1.0000]
> java.lang.RuntimeException: Problem in calculating the approximate path
> length to the plane.
> at
> org.hps.recon.tracking.WTrack.getPathLengthToPlaneApprox(WTrack.java:163)
> at org.hps.recon.tracking.WTrack.getHelixAndPlaneIntercept(WTrack.java:288)
> at
> org.hps.recon.tracking.TrackUtils.getHelixPlaneIntercept(TrackUtils.java:141)
> at
> org.hps.recon.tracking.MultipleScattering.getHelixIntersection(MultipleScattering.java:288)
> at
> org.hps.recon.tracking.MultipleScattering.getHelixIntersection(MultipleScattering.java:157)
> at
> org.hps.recon.tracking.MultipleScattering.FindHPSScatterPoints(MultipleScattering.java:96)
> at
> org.hps.recon.tracking.MultipleScattering.FindHPSScatters(MultipleScattering.java:59)
> at
> org.hps.recon.tracking.MultipleScattering.FindScatters(MultipleScattering.java:44)
> at
> org.lcsim.recon.tracking.seedtracker.HelixFitter.FitCandidate(HelixFitter.java:110)
> at
> org.lcsim.recon.tracking.seedtracker.SeedTrackFinder.FindTracks(SeedTrackFinder.java:140)
> at
> org.lcsim.recon.tracking.seedtracker.SeedTracker.process(SeedTracker.java:111)
> at org.lcsim.util.Driver.doProcess(Driver.java:273)
> at org.lcsim.util.Driver.processChildren(Driver.java:284)
> at org.lcsim.util.Driver.process(Driver.java:198)
> at
> org.hps.recon.tracking.TrackerReconDriver.process(TrackerReconDriver.java:172)
> at org.lcsim.util.Driver.doProcess(Driver.java:273)
> at org.lcsim.util.Driver.processChildren(Driver.java:284)
> at org.lcsim.util.Driver.process(Driver.java:198)
> at org.lcsim.util.DriverAdapter.recordSupplied(DriverAdapter.java:74)
> at
> org.freehep.record.loop.DefaultRecordLoop.consumeRecord(DefaultRecordLoop.java:832)
> at
> org.freehep.record.loop.DefaultRecordLoop.loop(DefaultRecordLoop.java:668)
> at
> org.freehep.record.loop.DefaultRecordLoop.execute(DefaultRecordLoop.java:566)
> at org.lcsim.util.loop.LCSimLoop.loop(LCSimLoop.java:151)
> at org.lcsim.job.JobControlManager.run(JobControlManager.java:418)
> at org.lcsim.job.JobControlManager.main(JobControlManager.java:180)
>
> On Tue, Oct 14, 2014 at 8:25 PM, Sebouh Paul <[log in to unmask]>
> wrote:
>
>
> thanks! i am running SimpleMCRecon.lcsim and it looks like it is
> working. I will update you if I need any more help.
>
> On Tue, Oct 14, 2014 at 8:05 PM, McCormick, Jeremy I. <
> [log in to unmask]> wrote:
>
>
> Hi,
>
> It could be that the recon steering configuration you are
> running expects the events to be run through the readout simulation first
> to simulate properly the timing and pile-up. This information is not
> present when you try to reconstruct directly the SLIC simulation output.
> (At least I think this is what that particular steering file expects as
> input.)
>
> What is the result if instead you use the following
> steering file to reconstruct directly your single particle events?
>
>
> steering-files/src/main/resources/org/hps/steering/recon/SimpleMCRecon.lcsim
>
> If this file isn't present in your trunk then 'svn up' to
> get it, as I just added it.
>
> I believe this should treat each event separately, and it
> should not depend on running the readout simulation first, so you should
> see reconstructed tracks.
>
> This is essentially what I use in the "MCReconTest" in
> integration-tests to test reconstruction of SLIC output directly without
> running the trigger simulation first.
>
> Your other option is running the readout simulation,
> either with the "NoPileup" steering file to treat each event separately, or
> inserting empty events so that your physics events will be properly spaced
> to create triggered events. (This is all covered both on the Confluence
> and in a bunch of previous emails to this list.)
>
> But let's see if the simpler option I suggested works
> first before you try to run the readout simulation...
>
> --Jeremy
>
> -----Original Message-----
> From: [log in to unmask] [mailto:
> [log in to unmask]] On Behalf Of Sebouh Paul
> Sent: Tuesday, October 14, 2014 2:22 PM
> To: hps-software
> Subject: no tracks found
>
> hey, i'm trying to do reconstruction, but I am not getting
> any tracks. below is the command i am using, and the output (it is done
> going through the entire file after only a few seconds), as well as the
> steering file I am using. The input I am using is muon radiative trident
> events.
>
>
> ifarm1101> java -Xmx1500m -jar
> ifarm1101>
> ~/hps-distribution-3.0.3-20141007.165651-110-bin.jar
> ifarm1101> ~/HPS2014OfflineRecon.lcsim -i
> ifarm1101>
> /work/hallb/hps/sebouh/trident_slic/radmuon_1.slcio
> ifarm1101> -DoutputFile=/work/hallb/hps/sebouh/test
> -DrunNumber=100
>
> No input files in XML file.
> Got ConditionsEvent with run: 0
> Reading calibrations calibSVT/base for run: 0 Use this
> calibration from run -1: calibSVT/default.base Reading calibrations
> calibSVT/tp for run: 0 Use this calibration from run -1: calibSVT/
> default.tp Loading the SVT bad channels for run 0 File
> daqmap/svt0.badchannels was not found! Continuing with only QA bad channels
> Loading SVT gains ...
> Loading SVT t0 shifts ...
> Loading fieldmap for run 0
> reading ECal DAQ map
> reading ECal bad channels
> reading pedestals for ECal
> reading pedestals for ECal
> name ""
> >> Event 0
> >> Event 1000
> >> Event 2000
> >> Event 3000
> >> Event 4000
> >> Event 5000
> GBLOutputDriver: Total Number of Events = 5636
> GBLOutputDriver: Total Number of Tracks = 0
> GBLOutputDriver: Total Number of Tracks Processed = 0 End
> of file reached
>
> and the contents of the driver file are:
>
>
> ifarm1101> more ~/HPS2014OfflineRecon.lcsim
> <!--
> Offline reconstruction for 2014 (electron run) data.
> @author Sho Uemura <[log in to unmask]>
> @version $Id: HPS2014OfflineRecon.lcsim,v 1.7 2013/10/30
> 16:23:32 phansson Exp $
> -->
> <lcsim xmlns:xs="http://www.w3.org/2001/XMLSchema-instance
> "
> xs:noNamespaceSchemaLocation="
> http://www.lcsim.org/schemas/lcsim/1.0/lcsim.xsd">
> <execute>
> <driver name="EventMarkerDriver"/>
> <driver name="CalibrationDriver"/>
> <driver name="RawTrackerHitSensorSetup"/>
> <driver name="RawTrackerHitFitterDriver" />
> <driver name="TrackerHitDriver"/>
> <driver name="HelicalTrackHitDriver"/>
> <driver name="TrackerReconDriver"/>
> <driver name="EcalRawConverter" />
> <driver name="EcalClusterer" />
> <driver name="ReconParticle" />
> <driver name="TrackDataDriver" />
> <driver name="GBLDriver"/>
> <driver name="LCIOWriter"/>
> <driver name="CleanupDriver"/>
> </execute>
> <drivers>
> <driver name="CalibrationDriver"
> type="org.hps.conditions.deprecated.CalibrationDriver">
> <!--
> <runNumber>${runNumber}</runNumber>-->
> </driver>
> <driver name="EventMarkerDriver"
> type="org.lcsim.job.EventMarkerDriver">
> <eventInterval>1000</eventInterval>
> </driver>
> <driver name="RawTrackerHitSensorSetup"
> type="org.lcsim.recon.tracking.digitization.sisim.config.RawTrackerHitSensorSetup"/>
> <driver name="RawTrackerHitFitterDriver"
> type="org.hps.recon.tracking.RawTrackerHitFitterDriver">
> <fitAlgorithm>Pileup</fitAlgorithm>
> <useTimestamps>true</useTimestamps>
> <correctT0Shift>false</correctT0Shift>
> <useTruthTime>false</useTruthTime>
> <subtractTOF>true</subtractTOF>
> <debug>false</debug>
> </driver>
> <driver name="TrackerHitDriver"
> type="org.hps.recon.tracking.DataTrackerHitDriver">
> <neighborDeltaT>8.0</neighborDeltaT>
> </driver>
> <driver name="HelicalTrackHitDriver"
> type="org.hps.recon.tracking.HelicalTrackHitDriver">
> <debug>false</debug>
> <clusterTimeCut>12.0</clusterTimeCut>
> <maxSeperation>20.0</maxSeperation>
> <tolerance>1.0</tolerance>
> </driver>
> <driver name="TrackerReconDriver"
> type="org.hps.recon.tracking.TrackerReconDriver">
>
> <strategyResource>HPS-Full.xml</strategyResource>
> <debug>false</debug>
> </driver>
> <driver name="GBLDriver"
> type="org.hps.recon.tracking.gbl.GBLOutputDriver">
> <debug>0</debug>
> <isMC>false</isMC>
> <gblFileName></gblFileName>
> </driver>
>
> <driver name="EcalRawConverter"
> type="org.hps.recon.ecal.EcalRawConverterDriver">
>
> <ecalCollectionName>EcalCalHits</ecalCollectionName>
> <use2014Gain>true</use2014Gain>
> <useTimestamps>true</useTimestamps>
> <useTruthTime>false</useTruthTime>
> </driver>
> <driver name="EcalClusterer"
> type="org.hps.recon.ecal.EcalClusterICBasic">
> <ecalName>Ecal</ecalName>
>
> <ecalCollectionName>EcalCalHits</ecalCollectionName>
> <timeCut>true</timeCut>
> </driver>
> <driver name="ReconParticle"
> type="org.hps.recon.particle.HpsReconParticleDriver">
> </driver>
> <driver name="TrackDataDriver"
> type="org.hps.recon.tracking.TrackDataDriver" />
> <driver name="LCIOWriter"
> type="org.lcsim.util.loop.LCIODriver">
>
> <outputFilePath>${outputFile}.slcio</outputFilePath>
> </driver>
> <driver name="CleanupDriver"
> type="org.lcsim.recon.tracking.digitization.sisim.config.ReadoutCleanupDriver"/>
> </drivers>
> </lcsim>
>
>
> ________________________________
>
>
> Use REPLY-ALL to reply to list
>
> To unsubscribe from the HPS-SOFTWARE list, click the
> following link:
>
> https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=HPS-SOFTWARE&A=1
>
>
>
> ########################################################################
>
> Use REPLY-ALL to reply to list
>
> To unsubscribe from the HPS-SOFTWARE list, click the
> following link:
>
> https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=HPS-SOFTWARE&A=1
>
>
>
>
>
########################################################################
Use REPLY-ALL to reply to list
To unsubscribe from the HPS-SOFTWARE list, click the following link:
https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=HPS-SOFTWARE&A=1
