-- SajanEaso - 21 Jan 2014

Options for Target study

The tool fires particles into blocks of material and returns information on probability of hadronic interaction and multiplicity and composition of daughters.
To run the tool ensure that you have the latest versions of Sim/GaussTools and Gen/LbPGuns

getpack Sim/SimChecks
cd scripts/Target

Remember to setup Gauss "SetupProject Gauss (your version)"

Single Target

This options can be used to obtain rootfiles for one specific combination of initial particle ID and energy and target material and thickness.
In the next secttion an option to run many combinations at once is described.
To run the tool setup Gauss and then:

cd SingleTarget

gaudirun.py Gauss-LHCb-Target-2013.py Gauss-Job-SingleTarget.py

Inside Gauss-Job-SingleTarget.py you can instead choose number of events generate and name of output files.

Target properties:

Inside Gauss-LHCb-Target-2013.py you can change the target properties by changing the line

target = 'Target_(thickness)mm(material)'

Available targets are 1, 5, 10 mm and availble maretials are 'Al' 'Be' 'Si'. (e.g. target = 'Target_1mmAl')

Incoming particle properties:

You can copy the file $LBPGUNSROOT/options/TargetMaterialGun.py. This contains the definition of the particle gun.
In particular the following two lines modify ID and energy of the particle.

ParticleGun.MaterialEval.PdgCode = 2212 (proton)

ParticleGun.MaterialEval.ModP = 3 * GeV

Before running change the following line in "Gauss-Job-SingleTarget.py" to pick up your pgun option file.

importOptions ("$LBPGUNSROOT/options/TargetMaterialGun.py")

Hadronic interaction model:

You can choose the hadronic interaction model to test by changing the following line in replacint QGSP_BERT with one of the models available in you Gauss version. (es. FTFP_BERT, LHEP)

Gauss.PhysicsList = {"Em":'NoCuts', "Hadron":'QGSP_BERT', "GeneralPhys":True, "LHCbPhys":True}

Multiple Target

This options file allow to run the tool at once for many patricle types and energies, for many target materials and thicknesses and many hadronic interaction models.

cd MultiTarget

In the file Targer_CreateEvents.py you'll find the following lines which you can modify to run with your choice of energies, materials, etc.

models=['QGSP_BERT','FTFP_BERT','LHEP'] #any present in the version of Gauss you are using
energies=['1', '5', '10', '50', '100'] #any
materialsTodo=['Al'] # 'Al' 'Be' 'Si'
thicks=['1'] #in mm 1, 5, 10
particlesTodo=['p', 'pbar', 'Piminus', 'Piplus', 'Kminus', 'Kplus'] # Available: 'Piminus' 'Piplus' 'Kminus' 'Kplus' 'p' 'pbar'

Finally to run the tool: python Target_CreateEvents.py

This creates a "TargetOutput.root" file containing a tree with information about multiplicity and cross-section for all targets and pguns.

To get histograms you can look at the script "TargetPlot.py". This works in a similar way as Targer_CreateEvents.py and you need to set the same options and then run

python TargetPlots.py

In the MultiTarget folder there is a subfolder "options" containing some option files for specific options combinations. E.g. MakeEvents_1mm_Al_AllPguns.py will genetare all needed data for all pguns but onlu with the 1mm Al target. And will also produce plots in the end. Options beginning with "MakeEvents" make events and automatically make plots out of them. Instead the options beginning with "Plot" will perform plots on already existing data.

Edit | Attach | Watch | Print version | History: r8 < r7 < r6 < r5 < r4 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r7 - 2014-11-12 - LucaPescatore
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    LHCb All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2022 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback