diff --git a/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_11_NEWJET.py b/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_11_NEWJET.py index 9e2251e..9636bd4 100644 --- a/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_11_NEWJET.py +++ b/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_11_NEWJET.py @@ -3,10 +3,10 @@ import math from Gaudi.Configuration import * #from GaudiPython import AppMgr -from Configurables import DecayTreeTuple,BTaggingTool,SubstitutePID,TrackScaleState -from Configurables import CombineParticles, FilterDesktop,TupleToolDecayTreeFitter,TupleToolDecay -from Configurables import TupleToolTrigger, BackgroundCategory, TupleToolTISTOS,TupleToolRecoStats, TupleToolTagging -from Configurables import DaVinci,HltSelReportsDecoder,HltVertexReportsDecoder,HltDecReportsDecoder,LoKi__Hybrid__TupleTool,TupleToolJets +from Configurables import DecayTreeTuple, BTaggingTool, SubstitutePID, TrackScaleState +from Configurables import CombineParticles, FilterDesktop, TupleToolDecayTreeFitter, TupleToolDecay +from Configurables import TupleToolTrigger, BackgroundCategory, TupleToolTISTOS, TupleToolRecoStats, TupleToolTagging +from Configurables import DaVinci, HltSelReportsDecoder, HltVertexReportsDecoder, HltDecReportsDecoder, LoKi__Hybrid__TupleTool, TupleToolJets from DecayTreeTuple.Configuration import * from StrippingConf.Configuration import StrippingConf, StrippingStream from StrippingSettings.Utils import strippingConfiguration @@ -23,17 +23,17 @@ 'L0DiMuonDecision', - # 'Hlt1SingleMuonNoIPDecision', - # 'Hlt1SingleMuonHighPTDecision', - # 'Hlt1DiMuonHighMassDecision', - # 'Hlt1DiMuonLowMassDecision', + # 'Hlt1SingleMuonNoIPDecision', + # 'Hlt1SingleMuonHighPTDecision', + # 'Hlt1DiMuonHighMassDecision', + # 'Hlt1DiMuonLowMassDecision', 'Hlt1GlobalDecision', 'Hlt1TrackAllL0Decision', 'Hlt1SingleMuonNoIPDecision', 'Hlt1SingleMuonHighPTDecision', 'Hlt1TrackMuonDecision', 'Hlt1MuTrackDecision', - # 'Hlt1GlobalDecision', + # 'Hlt1GlobalDecision', 'Hlt2TopoOSTF2BodyDecision', 'Hlt2TopoOSTF3BodyDecision', @@ -50,23 +50,23 @@ 'Hlt2TopoE2BodyBBDTDecision', 'Hlt2TopoE3BodyBBDTDecision', 'Hlt2TopoE4BodyBBDTDecision', - # 'Hlt2MuonFromHLT1Decision', + # 'Hlt2MuonFromHLT1Decision', 'Hlt2SingleMuonDecision', - # 'Hlt2SingleMuonHighPTDecision', - # 'Hlt2SingleMuonLowPTDecision', - # 'Hlt2DisplVerticesLowMassSingleDecision', - # 'Hlt2DisplVerticesHighMassSingleDecision', - # 'Hlt2DisplVerticesDoubleDecision', - # 'Hlt2DisplVerticesSinglePostScaledDecision', - # 'Hlt2DisplVerticesHighFDSingleDecision', - # 'Hlt2DisplVerticesSingleDownDecision', - # 'Hlt2DiMuonDecision', - # 'Hlt2DiMuonLowMassDecision', - # 'Hlt2DiMuonBDecision', + # 'Hlt2SingleMuonHighPTDecision', + # 'Hlt2SingleMuonLowPTDecision', + # 'Hlt2DisplVerticesLowMassSingleDecision', + # 'Hlt2DisplVerticesHighMassSingleDecision', + # 'Hlt2DisplVerticesDoubleDecision', + # 'Hlt2DisplVerticesSinglePostScaledDecision', + # 'Hlt2DisplVerticesHighFDSingleDecision', + # 'Hlt2DisplVerticesSingleDownDecision', + # 'Hlt2DiMuonDecision', + # 'Hlt2DiMuonLowMassDecision', + # 'Hlt2DiMuonBDecision', 'Hlt2DiMuonDetachedDecision', 'Hlt2MuTrackDecision', 'Hlt2GlobalDecision' - ] +] ########################################################################### @@ -76,7 +76,7 @@ eventNodeKiller.Nodes = [ '/Event/AllStreams', '/Event/Strip' ] -stripping='stripping20' +stripping = 'stripping20' #get the configuration dictionary from the database config = strippingConfiguration(stripping) #get the line builders from the archive @@ -88,17 +88,17 @@ for stream in streams: for line in stream.lines: if line.name() in MyLines: - print 'Strippingb2LcMuXB2DMuNuXLine' - MyStream.appendLines( [ line ] ) + print 'Strippingb2LcMuXB2DMuNuXLine' + MyStream.appendLines( [ line ] ) from Configurables import ProcStatusCheck filterBadEvents = ProcStatusCheck() sc = StrippingConf( Streams = [ MyStream ], - MaxCandidates = 2000, - AcceptBadEvents = False, - BadEventSelection = filterBadEvents ) + MaxCandidates = 2000, + AcceptBadEvents = False, + BadEventSelection = filterBadEvents ) ########################################################################### @@ -107,19 +107,18 @@ Lc_star = CombineParticles("Lc_star") -Lc_star.Inputs = ["Phys/StdAllLoosePions/Particles" ,"Phys/Lc2PKPiforB2DMuNuX/Particles"] -Lc_star.DecayDescriptors = ["[Lambda_c(2595)+ -> Lambda_c+ pi+ pi-]cc"] +Lc_star.Inputs = ["Phys/StdAllLoosePions/Particles", "Phys/Lc2PKPiforB2DMuNuX/Particles"] +Lc_star.DecayDescriptors = ["[Lambda_c(2595)+ -> Lambda_c+ pi+ pi-]cc"] Lc_star.MotherCut = "(M<3000) & (M>2300) & (VFASPF(VCHI2/VDOF) < 7.0) " -Lc_star.DaughtersCuts ={ "pi+" : "PT>0", "pi-": "PT>0", "Lambda_c+": "ADMASS('Lambda_c+') < 50.0 *MeV" } +Lc_star.DaughtersCuts = { "pi+": "PT>0", "pi-": "PT>0", "Lambda_c+": "ADMASS('Lambda_c+') < 50.0 *MeV" } ############################################################################# Lb_SL = CombineParticles("Lb_SL") -Lb_SL.Inputs = [ "Phys/MuforB2DMuNuX/Particles", "Phys/Lc_star/Particles"] -Lb_SL.DecayDescriptors = [ "[Lambda_b0 -> Lambda_c(2595)+ mu-]cc" ] +Lb_SL.Inputs = [ "Phys/MuforB2DMuNuX/Particles", "Phys/Lc_star/Particles"] +Lb_SL.DecayDescriptors = [ "[Lambda_b0 -> Lambda_c(2595)+ mu-]cc" ] Lb_SL.MotherCut = "(M<6500) & (VFASPF(VCHI2/VDOF) < 7.0) & (BPVDIRA>0.999)" - ''' ################################################################################ Lb_SLSel = CombineParticles("Lb_SLSel") @@ -134,8 +133,6 @@ ''' - - ######################################################################## #decay = "[ B*_s20 -> ^(B- -> ^(J/psi(1S) -> ^mu+ ^mu-) ^K-) ^[pi+]CC]CC" decay0 = "[Lambda_b0 -> ^(Lambda_c(2595)+ -> ^(Lambda_c+ -> ^K- ^p+ ^pi+ ) ^pi+ ^pi-) ^mu- ]CC" @@ -147,17 +144,17 @@ Lb_SL_OfflineTree.Inputs = ["Phys/Lb_SL/Particles"] Lb_SL_OfflineTree.Decay = decay0 Lb_SL_OfflineTree.Branches = { -# "H_10": "[^H_10 -> (Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) mu- ) CELLjet]CC", - "Lambda_b0": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) mu- ]CC", -# "CELLjet": "[H_10 -> (Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) mu- ) ^CELLjet]CC", - "Lambda_c_STAR": "[Lambda_b0 -> ^(Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) mu- ]CC", - "mu": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) ^mu- ]CC", - "Lambda_c": "[Lambda_b0 -> (Lambda_c(2595)+ -> ^(Lambda_c+ -> K- p+ pi+ ) pi+ pi-) mu- ]CC", - "pi_plus": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) ^pi+ pi-) mu- ]CC", - "pi_minus": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ ^pi-) mu- ]CC", - "p": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- ^p+ pi+ ) pi+ pi-) mu- ]CC", - "pi_Lc": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ ^pi+ ) pi+ pi-) mu- ]CC", - "K_Lc": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> ^K- p+ pi+ ) pi+ pi-) mu- ]CC" + # "H_10": "[^H_10 -> (Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) mu- ) CELLjet]CC", + "Lambda_b0": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) mu- ]CC", + # "CELLjet": "[H_10 -> (Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) mu- ) ^CELLjet]CC", + "Lambda_c_STAR": "[Lambda_b0 -> ^(Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) mu- ]CC", + "mu": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) ^mu- ]CC", + "Lambda_c": "[Lambda_b0 -> (Lambda_c(2595)+ -> ^(Lambda_c+ -> K- p+ pi+ ) pi+ pi-) mu- ]CC", + "pi_plus": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) ^pi+ pi-) mu- ]CC", + "pi_minus": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ ^pi-) mu- ]CC", + "p": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- ^p+ pi+ ) pi+ pi-) mu- ]CC", + "pi_Lc": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ ^pi+ ) pi+ pi-) mu- ]CC", + "K_Lc": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> ^K- p+ pi+ ) pi+ pi-) mu- ]CC" } Lb_SL_OfflineTree.ToolList = [ "TupleToolKinematic", @@ -180,9 +177,7 @@ 'TupleToolPropertime', 'TupleToolRecoStats', 'TupleToolTrackPosition' - ] - - +] Lb_SL_OfflineTree.addTool( TupleToolTrigger, name='TupleToolTrigger' ) @@ -198,13 +193,13 @@ ######################################################################## from Configurables import TupleToolNeutrinoReco -Lb_SL_OfflineTree.addTool( TupleToolNeutrinoReco,name = "TupleToolNeutrinoReco" ) +Lb_SL_OfflineTree.addTool( TupleToolNeutrinoReco, name = "TupleToolNeutrinoReco" ) Lb_SL_OfflineTree.TupleToolNeutrinoReco.Verbose = True -Lb_SL_OfflineTree.TupleToolNeutrinoReco.MotherMass=5619.5 #MeV -Lb_SL_OfflineTree.TupleToolNeutrinoReco.ExtraName='mLb' +Lb_SL_OfflineTree.TupleToolNeutrinoReco.MotherMass = 5619.5 # MeV +Lb_SL_OfflineTree.TupleToolNeutrinoReco.ExtraName = 'mLb' from Configurables import TupleToolKinematic -Lb_SL_OfflineTree.addTool( TupleToolKinematic,name = "TupleToolKinematic" ) +Lb_SL_OfflineTree.addTool( TupleToolKinematic, name = "TupleToolKinematic" ) Lb_SL_OfflineTree.TupleToolKinematic.Verbose = True from Configurables import TupleToolGeometry @@ -221,10 +216,7 @@ Lb_SL_OfflineTree.TrackIsolation.MinConeAngle = 0.5 Lb_SL_OfflineTree.TrackIsolation.MaxConeAngle = 1.5 Lb_SL_OfflineTree.TrackIsolation.StepSize = 0.1 -Lb_SL_OfflineTree.TrackIsolation.Verbose= True - - - +Lb_SL_OfflineTree.TrackIsolation.Verbose = True ################################################## @@ -239,7 +231,7 @@ DaVinci().SkipEvents = 0 # Events to skip DaVinci().PrintFreq = 10000 DaVinci().DataType = "2011" # Default is "DC06" -DaVinci().Simulation = False +DaVinci().Simulation = False DaVinci().Lumi = True DaVinci().appendToMainSequence( [ eventNodeKiller ] ) DaVinci().appendToMainSequence( [ sc.sequence() ] ) @@ -255,5 +247,3 @@ #appMgr = AppMgr() #appMgr.initialize() #appMgr.run( DaVinci().EvtMax ) - - diff --git a/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_12_NEWJET.py b/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_12_NEWJET.py index 11163ae..8445eb7 100644 --- a/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_12_NEWJET.py +++ b/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_12_NEWJET.py @@ -220,11 +220,11 @@ ################################################## -#DaVinci().Input = ([ -#'root://eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/2012/ALLSTREAMS.DST/00025296/0000/00025296_00000002_1.allstreams.dst' -#'root://eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/2012/ALLSTREAMS.DST/00025296/0000/00025296_00000003_1.allstreams.dst', -#'root://eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/2012/ALLSTREAMS.DST/00025296/0000/00025296_00000004_1.allstreams.dst' -#]) +DaVinci().Input = ([ + "DATAFILE='PFN:00021317_00001395_1.semileptonic.dst' TYP='POOL_ROOTTREE' OPT='READ'" + #'root://eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/2012/ALLSTREAMS.DST/00025296/0000/00025296_00000003_1.allstreams.dst', + #'root://eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/2012/ALLSTREAMS.DST/00025296/0000/00025296_00000004_1.allstreams.dst' +]) from Configurables import DaVinci #DaVinci().EventPreFilters = fltrs.filters('Filters') DaVinci().EvtMax = -1 # Number of events