diff --git a/productions/local_test_data/Ds_eepi_data_local.py b/productions/local_test_data/Ds_eepi_data_local.py index ae7920e..4524000 100644 --- a/productions/local_test_data/Ds_eepi_data_local.py +++ b/productions/local_test_data/Ds_eepi_data_local.py @@ -38,8 +38,8 @@ 'L0HadronDecision', 'Hlt1TrackMVADecision', 'Hlt1TrackMVALooseDecision', - 'Hlt2RareCharmD2piMuMuOSDecision', - 'Hlt2RareCharmD2piEEOSDecision', + 'Hlt2RareCharmD2PiMuMuOSDecision', + 'Hlt2RareCharmD2PiEEOSDecision', ] #StdPions = DataOnDemand(Location = "Phys/StdAllLoosePions/Particles") @@ -79,7 +79,7 @@ RequiredSelections=[StrippingSel]) pion_Filter = FilterDesktop('pionFilter', Code="PPINFO(PROBNNpi)>0.2") -pionFilter_sel = Selection("Filtered_pion_sel", +pion_Filter_sel = Selection("Filtered_pion_sel", Algorithm=pion_Filter, RequiredSelections=[pion_from_Ds_sel]) @@ -97,7 +97,7 @@ phi_sel = Selection("phi_sel", Algorithm=phi, - RequiredSelections=[electron_from_Ds_sel]) + RequiredSelections=[electron_Filter_sel]) ############ #Remake the Ds @@ -108,7 +108,7 @@ Ds.MotherCut = "(M<6500) & (VFASPF(VCHI2/VDOF) < 6.0) & (BPVDIRA>0.999)" Ds_sel = Selection("Ds_sel", Algorithm=Ds, - RequiredSelections=[phi_sel, pionFilter_sel]) #the order matters!! It makes sense to put the rarer decay first in order to + RequiredSelections=[phi_sel, pion_Filter_sel]) #the order matters!! It makes sense to put the rarer decay first in order to #save time -> put the most distinctive decay first Ds_phi_pi_sel_seq = SelectionSequence('Ds_phi_pi_selection_seq', @@ -142,6 +142,14 @@ }) + +Ds_OfflineTree.Ds.addTupleTool('TupleToolDecayTreeFitter/ConsD') + +#Ds_OfflineTree.Ds.ConsD.constrainToOriginVertex = True +Ds_OfflineTree.Ds.ConsD.Verbose = True +Ds_OfflineTree.Ds.ConsD.daughtersToConstrain = ['phi(1020)'] +Ds_OfflineTree.Ds.ConsD.UpdateDaughters = True + decay_trees.append(Ds_OfflineTree) # Lb_SL_OfflineTree_WS = DecayTreeTuple("Lb_SL_OfflineTree_WS") @@ -249,7 +257,7 @@ #DaVinci().UserAlgorithms.extend(decay_trees) DaVinci().TupleFile = "Ds_eepi.root" -#DaVinci().Input = ['PFN:/disk/data12/lhcb/davide/Rphipi/00070442_00006584_1.charm.mdst'] +DaVinci().Input = ['PFN:/disk/gangadir/davide/workspace/davide/local_data/00069601_00000009_1.charm.mdst'] DaVinci().DataType = "2016" diff --git a/productions/local_test_data/Ds_mumupi_data_local.py b/productions/local_test_data/Ds_mumupi_data_local.py index 1b9f222..7aef226 100644 --- a/productions/local_test_data/Ds_mumupi_data_local.py +++ b/productions/local_test_data/Ds_mumupi_data_local.py @@ -38,8 +38,8 @@ 'L0HadronDecision', 'Hlt1TrackMVADecision', 'Hlt1TrackMVALooseDecision', - 'Hlt2RareCharmD2piMuMuOSDecision', - 'Hlt2RareCharmD2piEEOSDecision', + 'Hlt2RareCharmD2PiMuMuOSDecision', + 'Hlt2RareCharmD2PiEEOSDecision', ] #StdPions = DataOnDemand(Location = "Phys/StdAllLoosePions/Particles") @@ -82,7 +82,7 @@ # pion_Filter = FilterDesktop('pionFilter', Code="ALL") pion_Filter = FilterDesktop('pionFilter', Code="PPINFO(PROBNNpi)>0.2") -pionFilter_sel = Selection("Filtered_pion_sel", +pion_Filter_sel = Selection("Filtered_pion_sel", Algorithm=pion_Filter, RequiredSelections=[pion_from_Ds_sel]) @@ -100,7 +100,7 @@ phi_sel = Selection("phi_sel", Algorithm=phi, - RequiredSelections=[muon_from_Ds_sel]) + RequiredSelections=[muon_Filter_sel]) ############ #Remake the Ds @@ -111,7 +111,7 @@ Ds.MotherCut = "(M<6500) & (VFASPF(VCHI2/VDOF) < 6.0) & (BPVDIRA>0.999)" Ds_sel = Selection("Ds_sel", Algorithm=Ds, - RequiredSelections=[phi_sel, pionFilter_sel]) #the order matters!! It makes sense to put the rarer decay first in order to + RequiredSelections=[phi_sel, pion_Filter_sel]) #the order matters!! It makes sense to put the rarer decay first in order to #save time -> put the most distinctive decay first Ds_phi_pi_sel_seq = SelectionSequence('Ds_phi_pi_selection_seq', @@ -143,6 +143,14 @@ }) + +Ds_OfflineTree.Ds.addTupleTool('TupleToolDecayTreeFitter/ConsD') + +#Ds_OfflineTree.Ds.ConsD.constrainToOriginVertex = True +Ds_OfflineTree.Ds.ConsD.Verbose = True +Ds_OfflineTree.Ds.ConsD.daughtersToConstrain = ['phi(1020)'] +Ds_OfflineTree.Ds.ConsD.UpdateDaughters = True + decay_trees.append(Ds_OfflineTree) # Lb_SL_OfflineTree_WS = DecayTreeTuple("Lb_SL_OfflineTree_WS") @@ -250,7 +258,7 @@ #DaVinci().UserAlgorithms.extend(decay_trees) DaVinci().TupleFile = "Ds_mumupi.root" -#DaVinci().Input = ['PFN:/disk/data12/lhcb/davide/Rphipi/00070442_00006584_1.charm.mdst'] +DaVinci().Input = ['PFN:/disk/gangadir/davide/workspace/davide/local_data/00069601_00000009_1.charm.mdst'] #remove this for the grid and specify in ganga j.inputdata DaVinci().DataType = "2016" diff --git a/productions/local_test_mc/Ds_eepi_mc_local.py b/productions/local_test_mc/Ds_eepi_mc_local.py index 083d51a..461fd6f 100644 --- a/productions/local_test_mc/Ds_eepi_mc_local.py +++ b/productions/local_test_mc/Ds_eepi_mc_local.py @@ -38,8 +38,8 @@ 'L0HadronDecision', 'Hlt1TrackMVADecision', 'Hlt1TrackMVALooseDecision', - 'Hlt2RareCharmD2piMuMuOSDecision', - 'Hlt2RareCharmD2piEEOSDecision', + 'Hlt2RareCharmD2PiMuMuOSDecision', + 'Hlt2RareCharmD2PiEEOSDecision', ] #StdPions = DataOnDemand(Location = "Phys/StdAllLoosePions/Particles") @@ -78,7 +78,7 @@ RequiredSelections=[StrippingSel]) pion_Filter = FilterDesktop('pionFilter', Code="PPINFO(PROBNNpi)>0.2") -pionFilter_sel = Selection("Filtered_pion_sel", +pion_Filter_sel = Selection("Filtered_pion_sel", Algorithm=pion_Filter, RequiredSelections=[pion_from_Ds_sel]) @@ -96,7 +96,7 @@ phi_sel = Selection("phi_sel", Algorithm=phi, - RequiredSelections=[electron_from_Ds_sel]) + RequiredSelections=[electron_Filter_sel]) ############ #Remake the Ds @@ -107,7 +107,7 @@ Ds.MotherCut = "(M<6500) & (VFASPF(VCHI2/VDOF) < 6.0) & (BPVDIRA>0.999)" Ds_sel = Selection("Ds_sel", Algorithm=Ds, - RequiredSelections=[phi_sel, pionFilter_sel]) #the order matters!! It makes sense to put the rarer decay first in order to + RequiredSelections=[phi_sel, pion_Filter_sel]) #the order matters!! It makes sense to put the rarer decay first in order to #save time -> put the most distinctive decay first Ds_phi_pi_sel_seq = SelectionSequence('Ds_phi_pi_selection_seq',