====== Stripping ======
===== Idea =====
Stripping is the idea of only fully reconstructing interesting decays, which can easily be accessed afterwards. It can cause a lot of pain, though...
===== Documentation =====
A lot of documentation can be found under: https://twiki.cern.ch/twiki/bin/view/LHCb/LHCbStripping
===== How to get the name of the lines =====
To know how to access your line in your options file, you have to know the name of it. To find this out, do the following:
* Find out in which stream it is (Dimuon, Dielectron, EW,...)
* Execute the following script:
SetupProject DaVinci vXrYpZ
python
from StrippingSettings.Stripping12 import StreamCalibration
locations = {}
for line in StreamCalibration.stream.lines :
locations[ line.name()] = str(StreamCalibration.stream.name()) + '/' + str(line.outputLocation())
for line, loc in locations.iteritems() :
print line, loc
* This will give you something like: ''StrippingDiMuonIncLine Dimuon/Phys/DiMuonIncLine'', where the first one is your line and the second one the location.
From Stripping 13 on, you have to do something like:
from StrippingSelections.Utils import buildStreams
streams = buildStreams('Stripping13')
for stream in streams:
if(stream.name() == "Calibration"):
for line in stream.lines :
print line.name(), "\t", line.outputLocation()
===== How to access the event =====
To access the event, you have to set the correct path, f.ex:
MyParticles = DataOnDemand(Location = '/Event/Dimuon/Phys/DiMuonIncLine')
or you can set the ''RootInTES'' ...
MyParticles = DataOnDemand(Location = 'Phys/DiMuonIncLine')
MyGaudiSequence.RootInTES = "/Event/Dimuon"
Note that you should not set ''RootInTES'' in the stripping, as it can cause problems (but can be useful for debugging sometimes).
===== How to test stripping lines on the GRID =====
First you have to make sure, the SDSTs as well as the RAW files are accessible and are staged in (as in the test-samples).
In ganga then you can do the following (in addition to the usual things):
j.inputdata = DaVinci().readInputData('/someDirectory/something_SDST.py')
j.inputdata.depth = 2
where it is crucial that you set ''j.inputdata.depth = 2'' to make sure the RAW files can be accessed.