JeffersonLab / gluex_MCwrapper

1 stars 4 forks source link

geant=0 flag not working #135

Open psharp15 opened 1 week ago

psharp15 commented 1 week ago

I am running simulation on the container, and I want to save the generated files before they go into geant. I've been testing the flag "geant=0", expecting it to stop after the generation stage. However, it keeps going and ends up running into errors with mcsmear. When I run with "mcsmear=0" and "recon=0", I get the same error (See below). I see that gen_gcf_090303_000_geant4_smeared.hddm is made, but gen_gcf_090303_000_geant4.hddm is not made. (here: /work/halld2/home/psharp/simulation/SRC-CT_Simulation/geant0_test/90303_0/). Here is my run script: /work/halld2/home/psharp/simulation/SRC-CT_Simulation/geant0_test/run_AV18_C12.sh

Thank you for your help! -Phoebe

`BeamProperties: Parsing config file gen_gcf_090303_000_beam.conf

BeamProperties: Using flux from input ROOT histogram file

BeamProperties: Using fixed polarization = 0.4 Working on reservoir event 0 opening HDDM file: gen_gcf_090303_000.hddm Total cross section: 120.13 nb RUNNING MCSMEAR running MCsmear without folding in random background mcsmear -PTHREAD_TIMEOUT_FIRST_EVENT=3600 -PTHREAD_TIMEOUT=3000 -ogen_gcf_090303_000_geant4_smeared.hddm gen_gcf_090303_000_geant4.hddm JANA >> Enabling z compression of output HDDM file stream (default) JANA >> Enabling CRC data integrity check in output HDDM file stream JANA ERROR>> JANA ERROR>> xxxxxxxxxxxx Unable to open event source "gen_gcf_090303_000_geant4.hddm"! xxxxxxxxxxxx JANA ERROR>> JANA ERROR>> xxxxxxxxxxxx NO VALID EVENT SOURCES GIVEN !!! xxxxxxxxxxxx
JANA ERROR>> JANA >> JANA >>Telling all threads to quit ... JANA >>Launching threads . JANA >> JANA >>Telling all threads to quit ...

Closed HDDM file 0 event written to gen_gcf_090303_000_geant4_smeared.hddm JANA >>Merging event reader thread ... JANA >> 0 events processed (0 events read) Average rate: 0.0Hz

Something went wrong with mcsmear status code: 66 ending gluex_MC.py`

T-Britton commented 1 week ago

I think the error is less of an error and more of not properly handling the case as well as it should. Geant=0 merely tells the system to not run Geant. Smearing requires Geant to have been run, so it crashes. I think think we may want to consider adding a control branch that just ends after generation.

Geant=0 works as intended but obviously causes some confusion Thomas Britton

On Oct 14, 2024, at 2:50 PM, Phoebe Sharp @.***> wrote:



I am running simulation on the container, and I want to save the generated files before they go into geant. I've been testing the flag "geant=0", expecting it to stop after the generation stage. However, it keeps going and ends up running into errors with mcsmear. When I run with "mcsmear=0" and "recon=0", I get the same error (See below). I see that gen_gcf_090303_000_geant4_smeared.hddm is made, but gen_gcf_090303_000_geant4.hddm is not made. (here: /work/halld2/home/psharp/simulation/SRC-CT_Simulation/geant0_test/90303_0/). Here is my run script: /work/halld2/home/psharp/simulation/SRC-CT_Simulation/geant0_test/run_AV18_C12.sh

Thank you for your help! -Phoebe

`BeamProperties: Parsing config file gen_gcf_090303_000_beam.conf

BeamProperties: Using flux from input ROOT histogram file

BeamProperties: Using fixed polarization = 0.4 Working on reservoir event 0 opening HDDM file: gen_gcf_090303_000.hddm Total cross section: 120.13 nb RUNNING MCSMEAR running MCsmear without folding in random background mcsmear -PTHREAD_TIMEOUT_FIRST_EVENT=3600 -PTHREAD_TIMEOUT=3000 -ogen_gcf_090303_000_geant4_smeared.hddm gen_gcf_090303_000_geant4.hddm JANA >> Enabling z compression of output HDDM file stream (default) JANA >> Enabling CRC data integrity check in output HDDM file stream JANA ERROR>> JANA ERROR>> xxxxxxxxxxxx Unable to open event source "gen_gcf_090303_000_geant4.hddm"! xxxxxxxxxxxx JANA ERROR>> JANA ERROR>> xxxxxxxxxxxx NO VALID EVENT SOURCES GIVEN !!! xxxxxxxxxxxx JANA ERROR>> JANA >> JANA >>Telling all threads to quit ... JANA >>Launching threads . JANA >> JANA >>Telling all threads to quit ...

Closed HDDM file 0 event written to gen_gcf_090303_000_geant4_smeared.hddm JANA >>Merging event reader thread ... JANA >> 0 events processed (0 events read) Average rate: 0.0Hz

Something went wrong with mcsmear status code: 66 ending gluex_MC.py`

— Reply to this email directly, view it on GitHubhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_JeffersonLab_gluex-5FMCwrapper_issues_135&d=DwMCaQ&c=CJqEzB1piLOyyvZjb8YUQw&r=ccFffx721N71hPpKcJGvJIqY9RM4gBTuzp9ir7rze5Q&m=2SJE4WkQl5hS0ujwWbgq_AC9sgDw8-Vx_CG9ytkD1PJIXiyJZY7OViGuv1l1OTG2&s=G5FJPAb4D6hK2WG0xoewjJ_-QfNXM1sGK475tXloA_k&e=, or unsubscribehttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_notifications_unsubscribe-2Dauth_AFRO2BDTHJZFXL3M6VUNOPDZ3QG5RAVCNFSM6AAAAABP5SI53OVHI2DSMVQWIX3LMV43ASLTON2WKOZSGU4DMNZZHE3TAOA&d=DwMCaQ&c=CJqEzB1piLOyyvZjb8YUQw&r=ccFffx721N71hPpKcJGvJIqY9RM4gBTuzp9ir7rze5Q&m=2SJE4WkQl5hS0ujwWbgq_AC9sgDw8-Vx_CG9ytkD1PJIXiyJZY7OViGuv1l1OTG2&s=CJGOVhjPMrEsSRUJIBs0JnQfb9WEbgpr9N1DHLAJ7mo&e=. You are receiving this because you are subscribed to this thread.Message ID: @.***>

aaust commented 1 week ago

The usage text gives the impression that the process should end after this step:

geant=[0/1] where 0 means that the geant step and any subsequent step will not run (default is 1)

s6pepaul commented 5 days ago

I am working on this fix. It should do what the usage text says. Run generation and then stop.