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  January 2016

HPS-SOFTWARE January 2016

Subject:

Re: Pass4 MC

From:

Bradley T Yale <[log in to unmask]>

Reply-To:

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

Date:

Sun, 10 Jan 2016 01:19:32 +0000

Content-Type:

text/plain

Parts/Attachments:

Parts/Attachments

text/plain (203 lines)

Hi Jeremy,

Thanks a lot for finding this.
We could probably just replace the mass value in the copy of 'particle.tbl' with an awk command in the script as needed, the way it is done in the MadGraph runcards. I'll try to work it into the standard procedure.

Brad

________________________________________
From: McCormick, Jeremy I. <[log in to unmask]>
Sent: Friday, January 8, 2016 8:16 PM
To: Bradley T Yale; Omar Moreno
Cc: Graham, Mathew Thomas; Uemura, Sho; Nathan Baltzell; hps-software
Subject: RE: Pass4 MC

Hi, Bradley.

I figured out what is going on with this, so let me explain for future reference.

Geant4 has trouble with particles that it does not know about e.g. the PDG 622 particle in your stdhep events representing the A prime is not something in G4's particle database.  Since this particle is marked as INTERMEDIATE, SLIC will create a G4PrimaryParticle for it, but the primary will be immediately decayed into its assigned products via the G4UnknownDecay process.  This is because it does not know how to track this particle in the simulation.

I'm able to hack around this by using the particle.tbl input to SLIC which is designed for introducing new particle definitions into the simulation.

I copied and then modified slic/HEAD/data/particle.tbl and added a row for the A prime:

622 A_prime                 0      0.05        0.0       99999999999

This has an infinite lifetime assigned with the mass and charge from the stdhep file, but the actual lifetime will be set through the preassigned decay and overrides this dummy value.

So then when I run slic with the new particle definition...

slic -P particle.tbl [...]

I can see the particle then travels properly before it decays...

*********************************************************************************************************
* G4Track Information:   Particle = not defined,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0        0        0        0       676        0        0         0 tracking_volume initStep
    1   0.0658    -0.25     7.82       676        0     7.82      7.82 tracking_volume Decay

And the daughters both have the proper displacement...

*********************************************************************************************************
* G4Track Information:   Particle = e-,   Track ID = 3,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0   0.0658    -0.25     7.82       371        0        0         0 tracking_volume initStep
    1     98.3    -27.4 1.32e+03       371 2.42e-23 1.31e+03  1.31e+03 world_volume Transportation
    2      187    -51.9  2.5e+03       371 2.18e-23 1.19e+03   2.5e+03  OutOfWorld Transportation

*********************************************************************************************************
* G4Track Information:   Particle = e+,   Track ID = 4,   Parent ID = 1
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
    0   0.0658    -0.25     7.82       355        0        0         0 tracking_volume initStep
    1    -80.1    -57.5 1.32e+03       355 2.41e-23 1.31e+03  1.31e+03 world_volume Transportation
    2     -152     -109  2.5e+03       355 2.17e-23 1.19e+03   2.5e+03  OutOfWorld Transportation

This is something we dealt with in ILC studies, and there is an extension in SLIC that is activated which will read additional particle definitions, but you have to provide this supplemental table on the command line for it to work.

I'm not sure how to solve this properly in all generality, as each StdHep record can potentially define a different mass and charge for some unknown PDG ID.  But Geant4 (from what I know) needs to have a fixed G4ParticleDefinition before event processing starts.

Hope that helps somewhat.

--Jeremy


-----Original Message-----
From: Bradley T Yale [mailto:[log in to unmask]]
Sent: Friday, January 08, 2016 11:16 AM
To: McCormick, Jeremy I.; Omar Moreno
Cc: Graham, Mathew Thomas; Uemura, Sho; Nathan Baltzell; hps-software
Subject: Re: Pass4 MC

The A' samples are finished, but the daughters which are supposed to have displaced vertices still do not appear to:

~/hps_soft/lcio/trunk/bin/dumpevent /cache/mss/hallb/hps/production/pass4/slic/ap/1pt05/50/apsignalv2_displaced_HPS-EngRun2015-Nominal-v4-4-fieldmap_4.slcio 5000





The ct was set to be 1 mm, so applying the boost to the A' in the above example should have given (E/m)*ct = (0.931/0.05)*1 mm = 18.62 mm for the daughters.





As usual, they appear to still be correctly displaced in both vertex and creation time at the stdhep level.


See the attached plots, where only the daughter electrons are shown.




Representative stdhep file:


/work/hallb/hps/mc_production/lhe/ap/1pt05/50/apsignalv2_displaced_4.stdhep




Could there still be a problem with the stdhep file that makes SLIC unable to handle this properly?



________________________________

From: [log in to unmask] <[log in to unmask]> on behalf of Bradley T Yale <[log in to unmask]>
Sent: Thursday, January 7, 2016 2:27 AM
To: Omar Moreno
Cc: Mathew Thomas Graham; Uemura, Sho; Nathan Baltzell; [log in to unmask]
Subject: Re: Pass4 MC

A' samples are in progress.
These are being made using the version of 'lhe_tridents' that displaces the creation time of the daughters as well as the vertex.


Also, I'm making two different versions in this batch: one which decays at the target, and one that is supposed to have ct=1 mm (before the boost).


________________________________

From: [log in to unmask] <[log in to unmask]> on behalf of Omar Moreno <[log in to unmask]>
Sent: Wednesday, January 6, 2016 5:05 PM
To: Bradley T Yale
Cc: Omar Moreno; Mathew Thomas Graham; Uemura, Sho; Nathan Baltzell; [log in to unmask]
Subject: Re: Pass4 MC

Are there pure A' samples as well?


On Wed, Jan 6, 2016 at 1:26 PM, Bradley T Yale <[log in to unmask]> wrote:


        Cc'ing the software list...



        Samples of everything are here:

        /mss/hallb/hps/production/pass4/





        which used the same jar as pass4 data (hps-distribution-3.5-20151218.205540-15-bin.jar)




        Everything was made with dst_maker v0.10, and there are also some tritrig made with v0.9 for comparison, labeled 'OLDDST'.


        Let me know if you need anything else. I'll be traveling tomorrow, so will try to submit things tonight.


        Brad




________________________________

        From: [log in to unmask] <[log in to unmask]> on behalf of Omar Moreno <[log in to unmask]>
        Sent: Wednesday, January 6, 2016 4:05 PM
        To: Bradley T Yale; Mathew Thomas Graham; Uemura, Sho; Nathan Baltzell
        Subject: Pass4 MC

        Hi Bradley,


        Are you planning to run pass 4 MC? Both the recon output and DST structure have changed so we will need pass4 MC files soon.


        Thanks,


        --Omar Moreno




________________________________

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

March 2024
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