LISTSERV mailing list manager LISTSERV 16.5

Help for HPS-SOFTWARE Archives


HPS-SOFTWARE Archives

HPS-SOFTWARE Archives


HPS-SOFTWARE@LISTSERV.SLAC.STANFORD.EDU


View:

Message:

[

First

|

Previous

|

Next

|

Last

]

By Topic:

[

First

|

Previous

|

Next

|

Last

]

By Author:

[

First

|

Previous

|

Next

|

Last

]

Font:

Proportional Font

LISTSERV Archives

LISTSERV Archives

HPS-SOFTWARE Home

HPS-SOFTWARE Home

HPS-SOFTWARE  May 2015

HPS-SOFTWARE May 2015

Subject:

Re: Trigger acceptance in MC

From:

Nathan Baltzell <[log in to unmask]>

Reply-To:

Software for the Heavy Photon Search Experiment <[log in to unmask]>

Date:

Tue, 26 May 2015 19:30:54 -0400

Content-Type:

text/plain

Parts/Attachments:

Parts/Attachments

text/plain (252 lines)

Very nice picture explanation!  

Does a driver exist in hps-java to take care of this and merge it back into a
format like real data, independently of the clustering algorithm?  Shouldn't it?



On May 26, 2015, at 7:02 PM, Kyle McCarty <[log in to unmask]> wrote:

> Hello Holly,
> 
> You aren't using running pedestals in MC- it's not even using this driver
> 
> Do you mean that it's not used by anything in the driver (and could therefore just be removed entirely) or that it isn't present in the driver? It is definitely in my version of that steering.
> 
> Holly, what does your driver actually look at for hits? If it looks at raw hits, you are fine with Monte Carlo. If you look at FADC hits, though, you would need to be able to look across multiple events to form your clusters. Consider the graphic below:
> 
> <HitDist.png>
> 
> The raw hits (the unprocessed SLiC output) is converted into FADC hits, but the FADC hits that correspond to the cluster are actually spread across six events (FADC hits are only written out every other event to simulate a clock-cycle) different events. Thusly, to correctly form this cluster, you would need to retain all six of those events. Otherwise, you are probably forming three different clusters from the same cluster (or possible just losing energy, since the latter hits may below threshold).
> 
> - Kyle
> 
> On Tue, May 26, 2015 at 6:48 PM, Holly Vance <[log in to unmask]> wrote:
> Hi Pelle,
> 
> A few things to check:
> What does the timing of the hits in the cluster look like? There is a time cut set in the ReconClusterer, but it is probably going to be slightly different for MC (not sure).
> 
> In general, these problems tend to arise when making the hits into readout hits. (You aren't using running pedestals in MC- it's not even using this driver). I suspect the issue is in EcalRawConverter. 
> 
> On Tue, May 26, 2015 at 6:32 PM, Hansson Adrian, Per Ola <[log in to unmask]> wrote:
> Thanks Holly,
> 
> 
> 
> On May 26, 2015, at 3:11 PM, Holly Vance <[log in to unmask]> wrote:
> 
>> Hi Pelle,
>> 
>> I have a simple steering file here that can take in your non-spaced events:
>> org.hps.steering.users.holly.EcalReadoutNoPileUp
>> 
> 
> Cool. That is good to have. 
> 
>> This will treat each event independently and can tell you if you are having problems in your readout of the data (it will show you where the peaks are and if you did spacing correctly, etc). This uses the full Recon Clusterer which works fine with Monte Carlo (as least I have never had any problems with it). 
>> 
> 
> Right now I’m trying to use our official simulation (see other email, didn’t see your reply before sending it) which includes spaced out events. The energy of the recon clusters are a little weird. Positions doesn’t seem crazy though. Any ideas what can go wrong?
> 
> /pelle
> 
> 
>> -Holly
>> 
>> 
>> On Tue, May 26, 2015 at 5:50 PM, Kyle McCarty <[log in to unmask]> wrote:
>> Hello Pelle,
>> 
>> Your cluster energy distributions definitely look different from what I saw (see the attached PDF). What clustering algorithm are you using? You need to use the GTPClusterer algorithm for it work right with Monte Carlo; you might be using GTPOnlineClusterer, which only works for the EvIO readout. I'm not sure why you would get more When I run reconstruction for Monte Carlo, I used the following drivers (note that I didn't include tracking):
>> 	• org.hps.conditions.ConditionsDriver
>> 	• org.hps.readout.ecal.FADCEcalReadoutDriver
>> 	• org.hps.recon.ecal.EcalRawConverterDriver
>> 	• org.hps.recon.ecal.cluster.GTPClusterDriver
>> 	• org.hps.readout.ecal.FADCPrimaryTriggerDriver
>> You need the GTPClusterDriver to properly handle Monte Carlo hit readout, GTP clustering to make sure that you are simulating the hardware, and then the FADCPrimaryTriggerDriver simulates the hardware pair trigger. Are you using these same drivers?
>> 
>> - Kyle
>> 
>> On Tue, May 26, 2015 at 5:42 PM, Hansson Adrian, Per Ola <[log in to unmask]> wrote:
>> Hi Again,
>> 
>> following up on this topic.
>> 
>> Running recon:
>> 
>> java -jar target/hps-distribution-3.3.1-SNAPSHOT-bin.jar  /org/hps/steering/recon/EngineeringRun2015FullRecon.lcsim  -DoutputFile=outfile -i bunched-readout.slcio
>> 
>> I see ECal energies that don’t look right. Any idea?
>> 
>> I see about 180 tracks in these 1342 triggered events. They look ok but the number seems low?
>> 
>> Am I forgetting something again?
>> 
>> /Pelle
>> 
>> 
>> <Screen Shot 2015-05-26 at 2.22.35 PM.png>
>> <Screen Shot 2015-05-26 at 2.31.44 PM.png>
>> 
>> 
>> On May 26, 2015, at 1:43 PM, Hansson Adrian, Per Ola <[log in to unmask]> wrote:
>> 
>>> 
>>> Ok,
>>> 
>>> using 150 bunches on 10k events I get ~13% acceptance.
>>> 
>>> Thanks,
>>> pelle
>>> 
>>> 
>>> Trigger Processing Results
>>> Single-Cluster Cuts
>>> Total Clusters Processed     :: 20295
>>> Passed Seed Energy Cut       :: 14119
>>> Passed Hit Count Cut         :: 11938
>>> Passed Total Energy Cut      :: 10991
>>> 
>>> Cluster Pair Cuts
>>> Total Pairs Processed        :: 3168
>>> Passed Energy Sum Cut        :: 2942
>>> Passed Energy Difference Cut :: 2941
>>> Passed Energy Slope Cut      :: 2755
>>> Passed Coplanarity Cut       :: 1342
>>> 
>>> Trigger Count :: 1342
>>> 
>>> Trigger Module Cut Values:
>>> Seed Energy Low        :: 0.050
>>> Seed Energy High       :: 6.600
>>> Cluster Energy Low     :: 0.060
>>> Cluster Energy High    :: 0.630
>>> Cluster Hit Count      :: 2
>>> Pair Energy Sum Low    :: 0.200
>>> Pair Energy Sum High   :: 0.860
>>> Pair Energy Difference :: 0.540
>>> Pair Energy Slope      :: 0.6
>>> Pair Coplanarity       :: 30.0
>>> FADCPrimaryTriggerDriver: Trigger count: 1342
>>> 
>>> 
>>> 
>>> On May 25, 2015, at 3:23 PM, Kyle McCarty <[log in to unmask]> wrote:
>>> 
>>>> Hello Pelle,
>>>> 
>>>> If you are reconstructing everything from a SLiC file, you need to space out the events because each one contains an A' event and will produce weird pile-up otherwise. I usually insert around 150 empty events between each real event to ensure that the hits and clusters (which are displaced by around 60ish events from the source) do not overlap at all. You can do this with the command:
>>>> 
>>>> java -cp $HPS_JAVA org.hps.users.meeg.FilterMCBunches $INPUT $OUTPUT -e150 -a
>>>> 
>>>> You should definitely get better acceptance that that as well; I got somewhere between 15% - 20% acceptance for 40 MeV (I don't have the exact value on hand).
>>>> 
>>>> Note also that there are (or at least were last I checked) differences in how the Monte Carlo simulation is stored versus the EvIO. In the EvIO data, all the hits and clusters for an event are stored in the same event, but in Monte Carlo readout, they are spaced across several events with each event representing 2 ns of time. As such, some drivers that work for EvIO readout do not for the Monte Carlo (this mainly affects clustering). Sho can correct me if I am misrepresenting something here.
>>>> 
>>>> - Kyle
>>>> 
>>>> On Mon, May 25, 2015 at 6:08 PM, Hansson Adrian, Per Ola <[log in to unmask]> wrote:
>>>> Just realized, how is the bunch spacing simulated in the readout step? I suppose I need to add some amount of “time” between the events or do a “NoPileUp” simulation...
>>>> 
>>>> /Pelle
>>>> 
>>>> On May 25, 2015, at 3:04 PM, Hansson Adrian, Per Ola <[log in to unmask]> wrote:
>>>> 
>>>>> Hi,
>>>>> 
>>>>> Have some issues getting the MC simulation to run. 
>>>>> 
>>>>> I’m running readout and recon over some 1.1GeV 40MeV A' MC files locally with the trunk. Slic file is here:
>>>>> 
>>>>> /nfs/slac/g/hps3/users/phansson/ap1.1gev40mevall_1_200ux40u_beamspot_gammactau_0cm_30mrad_SLIC-v04-00-00_Geant4-v10-00-02_QGSP_BERT_HPS-EngRun2015-Nominal-v1.slcio  
>>>>> 
>>>>> I’m using Nominal-v1 detector and I simulate readout with:
>>>>> 
>>>>> java -jar target/hps-distribution-3.3.1-SNAPSHOT-bin.jar  ../kepler2/hps-java/steering-files/src/main/resources/org/hps/steering/readout/EngineeringRun2015TrigPairs1.lcsim -Ddetector=HPS-EngRun2015-Nominal-v1 -Drun=2000
>>>>> 
>>>>> 
>>>>> I see 245 accepted events out of 10k events which is much lower than I would naively expect?  
>>>>> 
>>>>> 
>>>>> I then run recon with:
>>>>> 
>>>>> java -jar target/hps-distribution-3.3.1-SNAPSHOT-bin.jar  ../kepler2/hps-java/steering-files/src/main/resources/org/hps/steering/recon/EngineeringRun2015FullRecon.lcsim  -DoutputFile=outfiles/ap1.1gev40mevall_1_200ux40u_beamspot_gammactau_0cm_30mrad_SLIC-v04-00-00_Geant4-v10-00-02_QGSP_BERT_HPS-EngRun2015-Nominal-v1-recon -i outfiles/ap1.1gev40mevall_1_200ux40u_beamspot_gammactau_0cm_30mrad_SLIC-v04-00-00_Geant4-v10-00-02_QGSP_BERT_HPS-EngRun2015-Nominal-v1-readout.slcio
>>>>> 
>>>>> And I basically get to the 2nd event and then it just hangs there for 5-10mins.
>>>>> 
>>>>> 
>>>>> 
>>>>> Are these the right steps with the latest and greatest sw? Any ideas?
>>>>> 
>>>>> /pelle
>>>>> 
>>>>> 
>>>>> 
>>>>> 
>>>>> 
>>>>> 
>>>>> 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
>>> 
>> 
>> 
>> 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
>> 
>> NOTE: This message was trained as non-spam. If this is wrong, please correct the training as soon as possible. 
>> Spam
>> Not spam
>> Forget previous vote
> 
> NOTE: This message was trained as non-spam. If this is wrong, please correct the training as soon as possible. 
> Spam
> Not spam
> Forget previous vote
> 
> 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

Top of Message | Previous Page | Permalink

Advanced Options


Options

Log In

Log In

Get Password

Get Password


Search Archives

Search Archives


Subscribe or Unsubscribe

Subscribe or Unsubscribe


Archives

February 2024
January 2024
December 2023
November 2023
October 2023
September 2023
August 2023
July 2023
June 2023
May 2023
April 2023
March 2023
February 2023
January 2023
December 2022
November 2022
October 2022
September 2022
August 2022
June 2022
April 2022
March 2022
February 2022
January 2022
December 2021
November 2021
October 2021
September 2021
August 2021
July 2021
June 2021
May 2021
April 2021
March 2021
February 2021
January 2021
December 2020
November 2020
October 2020
September 2020
August 2020
July 2020
June 2020
May 2020
April 2020
March 2020
February 2020
January 2020
December 2019
November 2019
October 2019
September 2019
August 2019
July 2019
June 2019
May 2019
April 2019
March 2019
February 2019
January 2019
December 2018
November 2018
October 2018
September 2018
August 2018
July 2018
June 2018
May 2018
April 2018
March 2018
February 2018
January 2018
December 2017
November 2017
October 2017
September 2017
August 2017
July 2017
June 2017
May 2017
April 2017
March 2017
February 2017
January 2017
December 2016
November 2016
October 2016
September 2016
August 2016
July 2016
June 2016
May 2016
April 2016
March 2016
February 2016
January 2016
December 2015
November 2015
October 2015
September 2015
August 2015
July 2015
June 2015
May 2015
April 2015
March 2015
February 2015
January 2015
December 2014
November 2014
October 2014
September 2014
August 2014
July 2014
June 2014
May 2014
April 2014
March 2014
February 2014
January 2014
December 2013
November 2013
October 2013
September 2013
August 2013
July 2013
June 2013
May 2013
April 2013
March 2013
February 2013
January 2013
December 2012
November 2012
October 2012
September 2012
August 2012
July 2012
June 2012
May 2012
April 2012
March 2012
February 2012
January 2012
December 2011
November 2011
October 2011
September 2011
August 2011
July 2011
June 2011
May 2011
April 2011
March 2011

ATOM RSS1 RSS2



LISTSERV.SLAC.STANFORD.EDU

Secured by F-Secure Anti-Virus CataList Email List Search Powered by the LISTSERV Email List Manager

Privacy Notice, Security Notice and Terms of Use