[Halld-offline] Geant3/Geant4 comparison MC

Wed Jul 18 00:12:12 EDT 2018

On Tue, Jul 17, 2018 at 11:10 PM Sean Dobbs <sdobbs at fsu.edu> wrote:

> Hi all,
> Please see some updates here:
> https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_JeffersonLab_HDGeant4_issues_54&d=DwIFaQ&c=lz9TcOasaINaaC3U7FbMev2lsutwpI4--09aP8Lu18s&r=oyKV5joTkJsuRYv6hh48IMTw3i-IrYD-ZUAHHU0DdAY&m=NNKFjuV4zWK9n8qpp8yLRSMBw8Qqj6gEO4OQe-9Avjo&s=-p4UHylXcMnyC_179sXsOlW0UuUe1Yd3QqnRyPOec5E&e= 
>
> ---Sean
>
> On Tue, Jul 17, 2018 at 7:22 PM Justin Stevens <jrsteven at jlab.org> wrote:
>
>> Hi Thomas,
>>
>> Thanks for putting these together.  First a technical point: there seems
>> to be an issue with the p=1 GeV and theta=70 degrees single track samples,
>> since the histograms are almost all empty.
>>
>> Now onto the other single track MC samples: The forward tracking momentum
>> resolution is pretty poor from G4 (see pages 2 and 3 for both WireBased and
>> TimeBased tracks in the single pi- at theta=4 and p=4 GeV
>> https://halldweb.jlab.org/gluex_sim/G3G4_Comparisons/ver01/pim_4_4_Hist_Reconstruction.pdf).
>> This is likely related to the # of reconstructed hits on tracks (page 17 of
>> the same PDF) where very few of the tracks have the maximum number of hits
>> expected for this region in G4.
>>
>> Interestingly the number of FDC wire and cathode hits, as well as pseudo
>> points, for this sample are basically the same (see below where G3 is left
>> and G4 is right).  Meaning the FDC hit-level objects are there, but there’s
>> a problem with tracking through the FDC.  Any ideas from the experts on why
>> this might be the case?
>>
>> A similar conclusion can be drawn from the single-particle samples at 15
>> degrees.
>>
>> -Justin
>>
>> FDCWireHits
>>
>>
>>
>> FDCCathodeHits
>>
>>
>> FDCPseudoHits
>>
>>
>> On Jul 17, 2018, at 2:04 PM, Thomas Britton <tbritton at jlab.org> wrote:
>>
>> Initial comparisons can be found
>> https://halldweb.jlab.org/gluex_sim/G3G4_Comparisons/ver01/ and also at
>> /volatile/halld/home/tbritton/REQUESTED_MC/
>>
>> *_Reconstruction are all the graphs from the
>> /Independent/Hist_Reconstruction/ directory in hd_root.root and *_Neutrals
>> come from /Independent/Hist_Neutrals/
>>
>> The prefixes should line up with the directories below.
>>
>> I haven't looked at all of the plots yet so if something looks wildly off
>> please let me know and I can try to regenerate things by Thursday.
>>
>> Thomas
>>
>> ----- Original Message -----
>> From: "Thomas Britton" <tbritton at jlab.org>
>> To: "GlueX Offline Software Email List" <halld-offline at jlab.org>
>> Sent: Monday, July 16, 2018 10:02:49 AM
>> Subject: [Halld-offline] Geant3/Geant4 comparison MC
>>
>> initial runs of particle gun comparisons can be found:
>>
>> +----------------------------------------------------------------------------+
>> | OutputLocation
>>                                                             |
>>
>> +----------------------------------------------------------------------------+
>> |
>> /volatile/halld/home/tbritton/REQUESTED_MC/prot_g3_4_4_20180712023219pm/   |
>> |
>> /volatile/halld/home/tbritton/REQUESTED_MC/prot_g3_4_15_20180712023438pm/  |
>> |
>> /volatile/halld/home/tbritton/REQUESTED_MC/prot_g3_1_70_20180712023505pm/  |
>> |
>> /volatile/halld/home/tbritton/REQUESTED_MC/pim_g3_4_15_20180712023845pm/   |
>> |
>> /volatile/halld/home/tbritton/REQUESTED_MC/pim_g3_1_70_20180712023858pm/   |
>> |
>> /volatile/halld/home/tbritton/REQUESTED_MC/gamma_g3_4_4_20180712023935pm/  |
>> |
>> /volatile/halld/home/tbritton/REQUESTED_MC/gamma_g3_4_15_20180712024002pm/ |
>> |
>> /volatile/halld/home/tbritton/REQUESTED_MC/gamma_g3_1_70_20180712024104pm/ |
>> |
>> /volatile/halld/home/tbritton/REQUESTED_MC/gamma_g4_1_70_20180712024133pm/ |
>> |
>> /volatile/halld/home/tbritton/REQUESTED_MC/gamma_g4_4_15_20180712024211pm/ |
>> |
>> /volatile/halld/home/tbritton/REQUESTED_MC/gamma_g4_4_4_20180712024219pm/  |
>> |
>> /volatile/halld/home/tbritton/REQUESTED_MC/pim_g4_1_70_20180712024321pm/   |
>> |
>> /volatile/halld/home/tbritton/REQUESTED_MC/pim_g4_4_15_20180712024344pm/   |
>> | /volatile/halld/home/tbritton/REQUESTED_MC/pim_g4_4_4_20180712024352pm/
>>    |
>> |
>> /volatile/halld/home/tbritton/REQUESTED_MC/prot_g4_4_4_20180712024455pm/   |
>> |
>> /volatile/halld/home/tbritton/REQUESTED_MC/prot_g4_4_15_20180712024515pm/  |
>> |
>> /volatile/halld/home/tbritton/REQUESTED_MC/prot_g4_1_70_20180712024531pm/  |
>>
>> +----------------------------------------------------------------------------+
>>
>> where prot_g3_4_15_20180712023438pm is formed from [particle
>> user]_g[3/4]_[P GeV]_[Theta degrees]_[Time Stamp]
>> thus the example above is for 4 GeV protons shot through geant3 at 15
>> degrees.
>>
>> The signal MC can be found:
>>
>>
>> +---------------------------------------------------------------------------+
>> | OutputLocation
>>                                                            |
>>
>> +---------------------------------------------------------------------------+
>> | /volatile/halld/home/tbritton/REQUESTED_MC/RHO_geant3_20180713024643pm/
>>   |
>> | /volatile/halld/home/tbritton/REQUESTED_MC/RHO_geant4_20180713024658pm/
>>   |
>> |
>> /volatile/halld/home/tbritton/REQUESTED_MC/OMEGA_geant4_20180713024842pm/ |
>> |
>> /volatile/halld/home/tbritton/REQUESTED_MC/OMEGA_geant3_20180713024852pm/ |
>> | /volatile/halld/home/tbritton/REQUESTED_MC/PHI_geant3_20180713024922pm/
>>   |
>> | /volatile/halld/home/tbritton/REQUESTED_MC/PHI_geant4_20180713024930pm/
>>   |
>> | /volatile/halld/home/tbritton/REQUESTED_MC/PI0_geant3_20180713025004pm/
>>   |
>>
>> +---------------------------------------------------------------------------+
>>
>> These use the same 4 vectors and are of the form
>> [Channel]_geant[3/4]_[Time Stamp].
>>                                                  RHO->pi+ pi-
>>                                                  OMEGA->pi+ pi- pi0
>>                                                  PHI->K+ K-
>>                                                  PI0->gg
>>
>>
>> +-------------------------------------------------------------------------+
>> | OutputLocation
>>                                                          |
>>
>> +-------------------------------------------------------------------------+
>> | /volatile/halld/home/tbritton/REQUESTED_MC/PI0_geant3_20180713025004pm/
>> |
>>
>> +-------------------------------------------------------------------------+
>>
>> contains two straggler jobs which are currently running....
>>
>> NOTE:  In All cases the smearing is actually smearing and the software is
>> "the latest and greatest": version_2.37_jlab.xml
>>
>> I will begin to process them using Sean's comparison scripts but they
>> should exist for others to start to look at them before the Geant4 meeting
>> on Friday.
>>
>> Thomas
>>
>>
>>
>>
>>
>> _______________________________________________
>> Halld-offline mailing list
>> Halld-offline at jlab.org
>> https://mailman.jlab.org/mailman/listinfo/halld-offline
>> _______________________________________________
>> Halld-offline mailing list
>> Halld-offline at jlab.org
>> https://mailman.jlab.org/mailman/listinfo/halld-offline
>>
>>
>> _______________________________________________
>> Halld-offline mailing list
>> Halld-offline at jlab.org
>> https://mailman.jlab.org/mailman/listinfo/halld-offline
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.jlab.org/pipermail/halld-offline/attachments/20180717/431e4c4a/attachment-0002.html>