#
# Configuration of SelDSTWriter
# per-event an per-line selective writing of the raw event is active (selectiveRawEvent=True)
#
SelDSTWriterElements = {
'default': stripDSTElements(pack=enablePacking),
charmMicroDSTname: mdstElements,
leptonicMicroDSTname: mdstElements,
pidMicroDSTname: mdstElements,
bhadronMicroDSTname: mdstElements
}
SelDSTWriterConf = {
'default':
stripDSTStreamConf(pack=enablePacking, selectiveRawEvent=True),
charmMicroDSTname:
mdstStreamConf,
leptonicMicroDSTname:
mdstStreamConf,
bhadronMicroDSTname:
mdstStreamConf,
pidMicroDSTname:
stripCalibMicroDSTStreamConf(pack=enablePacking, selectiveRawEvent=True)
}
dstWriter = SelDSTWriter("MyDSTWriter",
StreamConf=SelDSTWriterConf,
MicroDSTElements=SelDSTWriterElements,
OutputFileSuffix='20151106',
SelectionSequences=sc.activeStreams())
chic2ConstNoPV = tup.Lambda_b0.addTupleTool(
'TupleToolDecayTreeFitter/ConstChic2FitNoPV')
chic2ConstNoPV.UpdateDaughters = True
chic2ConstNoPV.Substitutions = {
'Beauty -> ^chi_c1(1P) Baryon Meson': 'chi_c2(1P)'
}
chic2ConstNoPV.daughtersToConstrain += ['J/psi(1S)', 'chi_c2(1P)']
##################################################################
# If we want to write a DST do this
##################################################################
from DSTWriters.microdstelements import *
from DSTWriters.Configuration import (SelDSTWriter, stripDSTStreamConf,
stripDSTElements)
SelDSTWriterElements = {'default': stripDSTElements()}
SelDSTWriterConf = {'default': stripDSTStreamConf()}
if MODE == 'MC':
dstWriter = SelDSTWriter("MyDSTWriter",
StreamConf=SelDSTWriterConf,
MicroDSTElements=SelDSTWriterElements,
OutputFileSuffix='MC',
SelectionSequences=sc.activeStreams())
###################### DAVINCI SETTINGS ############################################
DaVinci().SkipEvents = 0 #1945
DaVinci().PrintFreq = 10000
DaVinci().EvtMax = EVTMAX
DaVinci().TupleFile = "DVTuples1.root"
DaVinci().HistogramFile = 'DVHistos.root'
DaVinci().RootInTES = rootInTES
DaVinci().InputType = "MDST"
particle.addTool(TISTOSTool, name = "TISTOSTool")
particle.ToolList += [ "TupleToolTISTOS/TISTOSTool" ]
##################################################################
# If we want to write a DST do this
##################################################################
from DSTWriters.microdstelements import *
from DSTWriters.Configuration import (SelDSTWriter,
stripDSTStreamConf,
stripDSTElements
)
SelDSTWriterElements = {
'default' : stripDSTElements()
}
SelDSTWriterConf = {
'default' : stripDSTStreamConf()
}
if MODE == 'MC':
dstWriter = SelDSTWriter( "MyDSTWriter",
StreamConf = SelDSTWriterConf,
MicroDSTElements = SelDSTWriterElements,
OutputFileSuffix ='MC',
SelectionSequences = sc.activeStreams()
)
###################### DAVINCI SETTINGS ############################################
DaVinci().SkipEvents = 0 #1945
DaVinci().PrintFreq = 1000
DaVinci().EvtMax = EVTMAX
DaVinci().TupleFile = "DVTuples1.root"
DaVinci().HistogramFile = 'DVHistos.root'
sc = StrippingConf(
Streams=[MyStream],
MaxCandidates=2000,
AcceptBadEvents=False,
BadEventSelection=filterBadEvents,
HDRLocation='NonexistingPlaceholder',
)
#
# Configuration of SelDSTWriter
#
SelDSTWriterElements = {
'default': stripDSTElements(),
}
SelDSTWriterConf = {
'default': stripDSTStreamConf(),
}
dstWriter = SelDSTWriter("MyDSTWriter",
StreamConf=SelDSTWriterConf,
OutputFileSuffix='TAGGING',
SelectionSequences=sc.activeStreams())
# DaVinci Config
DaVinci().appendToMainSequence([
sc.sequence(),
dstWriter.sequence(),
])
from DSTWriters.microdstelements import *
from DSTWriters.Configuration import ( SelDSTWriter,
stripDSTStreamConf,
stripDSTElements,
stripMicroDSTStreamConf,
stripMicroDSTElements,
stripCalibMicroDSTStreamConf )
#
# Configuration of SelDSTWriter
# per-event an per-line selective writing of the raw event is active (selectiveRawEvent=True)
#
SelDSTWriterElements = { 'default' : stripDSTElements(pack=enablePacking) }
SelDSTWriterConf = { 'default' : stripDSTStreamConf(pack=enablePacking,selectiveRawEvent=True) }
dstWriter = SelDSTWriter( "MyDSTWriter",
StreamConf = SelDSTWriterConf,
MicroDSTElements = SelDSTWriterElements,
OutputFileSuffix ='000000',
SelectionSequences = sc.activeStreams() )
#
# Add stripping TCK
#
#from Configurables import StrippingTCK
#stck = StrippingTCK(HDRLocation = '/Event/Strip/Phys/DecReports', TCK=0x36112100)
#
#Configure DaVinci
#
mdstElements = stripMicroDSTElements(pack=enablePacking)
#
# Configuration of SelDSTWriter
# per-event an per-line selective writing of the raw event is active (selectiveRawEvent=True)
#
SelDSTWriterElements = {
'default' : stripDSTElements(pack=enablePacking),
charmMicroDSTname : mdstElements,
leptonicMicroDSTname : mdstElements,
pidMicroDSTname : mdstElements,
bhadronMicroDSTname : mdstElements
}
SelDSTWriterConf = {
'default' : stripDSTStreamConf(pack=enablePacking, selectiveRawEvent=True),
charmMicroDSTname : mdstStreamConf,
leptonicMicroDSTname : mdstStreamConf,
bhadronMicroDSTname : mdstStreamConf,
pidMicroDSTname : stripCalibMicroDSTStreamConf(pack=enablePacking, selectiveRawEvent=True)
}
dstWriter = SelDSTWriter( "MyDSTWriter",
StreamConf = SelDSTWriterConf,
MicroDSTElements = SelDSTWriterElements,
OutputFileSuffix ='000000',
SelectionSequences = sc.activeStreams() )
#
# Add stripping TCK
#
