[Halld-offline] Geant3/Geant4 comparison MC

Justin Stevens jrsteven at jlab.org
Tue Jul 17 18:38:11 EDT 2018


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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.jlab.org/pipermail/halld-offline/attachments/20180717/e9d95574/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: PastedGraphic-4.png
Type: image/png
Size: 121367 bytes
Desc: not available
URL: <https://mailman.jlab.org/pipermail/halld-offline/attachments/20180717/e9d95574/attachment-0003.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: PastedGraphic-5.png
Type: image/png
Size: 115599 bytes
Desc: not available
URL: <https://mailman.jlab.org/pipermail/halld-offline/attachments/20180717/e9d95574/attachment-0004.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: PastedGraphic-6.png
Type: image/png
Size: 112782 bytes
Desc: not available
URL: <https://mailman.jlab.org/pipermail/halld-offline/attachments/20180717/e9d95574/attachment-0005.png>


More information about the Halld-offline mailing list