Project

General

Profile

Support #1256

Projecting beams backwards (-ve Pz)

Added by Carlisle, Timothy almost 10 years ago. Updated almost 10 years ago.

Status:
Closed
Priority:
Normal
Assignee:
Rogers, Chris
Category:
Simulation
Target version:
Future MAUS release
Start date:
17 April 2013
Due date:
% Done:

100%

Estimated time:
Workflow:
New Issue

Description

I'm trying to track a bunch of particles back through my geometry, which I do by flipping px, py, pz and the PID (particle -> anti-particle). This all seemed fine in G4MICE, but for some reason Maus doesn't seem to let me...

TRANSFORM/MERGE/OUTPUT: Processed 1 events so far, 1 events in buffer.
TRANSFORM/MERGE/OUTPUT: Processed 2 events so far, 1 events in buffer.
Traceback (most recent call last):
File "/home/carlisle/maus0.5/src/common_py/ErrorHandler.py", line 162, in HandleCppException
raise CppError(error_message)
ErrorHandler.CppError: Particle is outside world volume at MAUSPrimaryGeneratorAction::GeneratePrimaries
Traceback (most recent call last):
File "/home/carlisle/maus0.5/src/common_py/ErrorHandler.py", line 162, in HandleCppException
raise CppError(error_message)
ErrorHandler.CppError: In branch errors
Failed to convert json value to string at StringProcessor::JsonToCpp
TRANSFORM/MERGE/OUTPUT: Processed 3 events so far, 1 events in buffer.

I can input beams (in icool for003 format) and send them down the axis, but back up doesn't work...am I missing something basic/subtle? Cards & for003 file attached.

Thanks Chris

Files

for003_TRef1.dat (927 KB) for003_TRef1.dat Carlisle, Timothy, 17 April 2013 19:50
Cards_31.py (844 Bytes) Cards_31.py Carlisle, Timothy, 17 April 2013 19:50
StepIV_none__8m.dat (684 Bytes) StepIV_none__8m.dat Carlisle, Timothy, 17 April 2013 19:51
#1

Updated by Carlisle, Timothy almost 10 years ago

#3

Updated by Rogers, Chris almost 10 years ago

  • Category set to Simulation
  • Target version set to Future MAUS release

I ran with --verbose_level 0 and got a bit more debugging output (below). Note the first error

ErrorHandler.CppError: Particle is outside world volume at MAUSPrimaryGeneratorAction::GeneratePrimaries

I tried setting the world volume to 6000*6000*31000 m and it looks like it works. You may have some bug in your for003 file (should be in [m]/[GeV] as per ICOOL specifications). Note I was missing some files so I couldn't run with the full geometry that you provided.

Stack trace:/home/cr67/MAUS/maus_release/build/libMausCpp.so(_ZN4MAUS26MAUSPrimaryGeneratorAction17GeneratePrimariesEP7G4Event+0x247) [0x7fd6f108c387]
/home/cr67/MAUS/maus_release/third_party/build/geant4.9.2.p04/lib/Linux-g++/libG4run.so(_ZN12G4RunManager13GenerateEventEi+0x10a) [0x7fd6ed19617a]
/home/cr67/MAUS/maus_release/third_party/build/geant4.9.2.p04/lib/Linux-g++/libG4run.so(_ZN12G4RunManager11DoEventLoopEiPKci+0xb6) [0x7fd6ed195aa6]
/home/cr67/MAUS/maus_release/third_party/build/geant4.9.2.p04/lib/Linux-g++/libG4run.so(_ZN12G4RunManager6BeamOnEiPKci+0x7c) [0x7fd6ed19566c]
/home/cr67/MAUS/maus_release/build/libMausCpp.so(_ZN4MAUS17MAUSGeant4Manager16RunManyParticlesEN4Json5ValueE+0x242) [0x7fd6f1091a72]
/home/cr67/MAUS/maus_release/build/_MapCppSimulation.so(_ZN4MAUS16MapCppSimulation7processESs+0xe0) [0x7fd6dfffc1a0]
/home/cr67/MAUS/maus_release/build/_MapCppSimulation.so(+0x77f2) [0x7fd6dffff7f2]
/home/cr67/MAUS/maus_release/third_party/install/lib/libpython2.7.so.1.0(PyEval_EvalFrameEx+0x52be) [0x7fd6fa497f0e]
/home/cr67/MAUS/maus_release/third_party/install/lib/libpython2.7.so.1.0(PyEval_EvalCodeEx+0x855) [0x7fd6fa499ce5]
/home/cr67/MAUS/maus_release/third_party/install/lib/libpython2.7.so.1.0(PyEval_EvalFrameEx+0x5245) [0x7fd6fa497e95]
/home/cr67/MAUS/maus_release/third_party/install/lib/libpython2.7.so.1.0(PyEval_EvalFrameEx+0x629a) [0x7fd6fa498eea]

Traceback (most recent call last):
  File "/home/cr67/MAUS/maus_release/src/common_py/ErrorHandler.py", line 162, in HandleCppException
    raise CppError(error_message)
ErrorHandler.CppError: Particle is outside world volume at MAUSPrimaryGeneratorAction::GeneratePrimaries
mc numevts = 1000
Stack trace:/home/cr67/MAUS/maus_release/build/libMausCpp.so(_ZN4MAUS15StringProcessor9JsonToCppERKN4Json5ValueE+0x90) [0x7fd6f1135190]
/home/cr67/MAUS/maus_release/build/libMausCpp.so(_ZN4MAUS23ObjectMapValueProcessorISsE9JsonToCppERKN4Json5ValueE+0xd2) [0x7fd6f115e3c2]
/home/cr67/MAUS/maus_release/build/libMausCpp.so(_ZN4MAUS17ObjectProcessorNS9ValueItemINS_5SpillESt3mapISsSsSt4lessISsESaISt4pairIKSsSsEEEE12_SetCppChildERKN4Json5ValueERS2_+0x97) [0x7fd6f115c1a7]
/home/cr67/MAUS/maus_release/build/libMausCpp.so(_ZN4MAUS15ObjectProcessorINS_5SpillEE9JsonToCppERKN4Json5ValueE+0x60) [0x7fd6f115a9a0]
/home/cr67/MAUS/maus_release/build/libMausCpp.so(_ZNK4MAUS21JsonCppSpillConverter8_convertEPKN4Json5ValueE+0x76) [0x7fd6f120de56]
/home/cr67/MAUS/maus_release/build/_OutputCppRoot.so(_ZN4MAUS13OutputCppRoot11write_eventINS_21JsonCppSpillConverterENS_5SpillEEEbPNS_9MAUSEventIT0_EERKN4Json5ValueESs+0x1a4) [0x7fd6df2676d4]
/home/cr67/MAUS/maus_release/build/_OutputCppRoot.so(_ZN4MAUS13OutputCppRoot5_saveESs+0x23a) [0x7fd6df262f8a]
/home/cr67/MAUS/maus_release/build/_OutputCppRoot.so(_ZN4MAUS10OutputBaseISsE4saveESs+0x34) [0x7fd6df2638e4]
/home/cr67/MAUS/maus_release/build/_OutputCppRoot.so(+0x15dd7) [0x7fd6df26cdd7]
/home/cr67/MAUS/maus_release/third_party/install/lib/libpython2.7.so.1.0(PyEval_EvalFrameEx+0x52be) [0x7fd6fa497f0e]
/home/cr67/MAUS/maus_release/third_party/install/lib/libpython2.7.so.1.0(PyEval_EvalCodeEx+0x855) [0x7fd6fa499ce5]
/home/cr67/MAUS/maus_release/third_party/install/lib/libpython2.7.so.1.0(PyEval_EvalFrameEx+0x5245) [0x7fd6fa497e95]
/home/cr67/MAUS/maus_release/third_party/install/lib/libpython2.7.so.1.0(PyEval_EvalFrameEx+0x629a) [0x7fd6fa498eea]

Traceback (most recent call last):
  File "/home/cr67/MAUS/maus_release/src/common_py/ErrorHandler.py", line 162, in HandleCppException
    raise CppError(error_message)
ErrorHandler.CppError: In branch errors
Failed to convert json value to string at StringProcessor::JsonToCpp
#4

Updated by Carlisle, Timothy almost 10 years ago

Ah ok, never had volume be an issue before!

One (minor) thing, is there a way to apply a cut on the absolute value of a variable, i.e. for pz? Something like...cue incorrect pseudocode(!)...:

bunch.cut({ fabs('pz'):150.}, operator.le, global_cut=False)

If there isn't a simple Maus-way, I can work round it straightforwardly enough (if statement etc), just thought I'd ask.

Thanks

#5

Updated by Rogers, Chris almost 10 years ago

I think the pseudocode you posted was even correct. This will set statistical weight to 0. If you want to remove the particle altogether there is Bunch.conditional_remove

#6

Updated by Rogers, Chris almost 10 years ago

The documentation I linked to was crap. I rewrote it:

    """ 
    For when a cut is not good enough, remove hits from the bunch altogether if
    comparator(value, get([variable])) is True for any variable

    - variable_value_dict = dict of Bunch.get_variables() strings to the cut 
                            value; use hit['variable'] for the cut
    - comparator = comparator function; if comparator(Hit['variable']) is true,
                   remove particle from the bunch
    - value_is_nsigma_bool = boolean; if True, value is the number of standard
                             deviations

    e.g. bunch.cut({'energy':300, 'x':20.}, operator.ge) removes particles if
    they have energy >= 300 MeV OR x >= 20 mm. Here operator.ge is a functional 
    representation of the >= operator. ge(hit.get('energy'), 300) is the same as 
    (hit.get('energy') >= 300).

    A whole load of useful comparators can be found in the operator module. e.g.
    ge is >=; le is <=; etc. See also python built-in function "filter".
    """
#7

Updated by Rogers, Chris almost 10 years ago

  • Status changed from Open to Closed
  • % Done changed from 0 to 100

Closing the issue if that's okay.

#8

Updated by Carlisle, Timothy almost 10 years ago

yep, thanks.

Also available in: Atom PDF