Skip to content
Snippets Groups Projects
Commit 7165dd7a authored by Dey sdey Sourav's avatar Dey sdey Sourav
Browse files

added topoana instructions

parent 68cf2c9a
Branches main
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
PartialReco_Dpluspi.card 0 → 100644
+ 121
− 0
View file @ 7165dd7a
% Names of input root files
{
dummy.root
}
% TTree name
{
dpluspi
}
% Component analysis --- decay trees
{
Y 100
}
% Process charge conjugate objects together (Two options: Y and N. Default: N)
{
Y
}
% PDG code of ISR photons (Default: 222222222)
{
222222222
}
% PDG code of FSR photons (Default: -22)
{
-22
}
% Ignore ISR photons (Three options: Ys, Yg and N. Default: N)
{
Yg
}
% Ignore FSR photons (Three options: Ys, Yg and N. Default: N)
{
Yg
}
% Signal identification --- decay branches
{
D+ --> K- pi+ pi+ & DplusisSig
D0 --> K- pi+ & D0isSig
}
% Signal identification --- cascade decay branches
{
0 & B- --> D0 pi- & -1 & B2D0pim
1 & D0 --> K- pi+ & 0
0 & B- --> D0 K- & -1 & B2D0Km
1 & D0 --> K- pi+ & 0
0 & B- --> D*0 pi- & -1 & B2Dst0pim_pi0
1 & D*0 --> D0 pi0 & 0
2 & D0 --> K- pi+ & 1
0 & B- --> D*0 pi- & -1 & B2Dst0pim_gamma
1 & D*0 --> D0 gamma & 0
2 & D0 --> K- pi+ & 1
0 & B- --> D*0 K- & -1 & B2Dst0Km_pi0
1 & D*0 --> D0 pi0 & 0
2 & D0 --> K- pi+ & 1
0 & B- --> D*0 K- & -1 & B2Dst0Km_gamma
1 & D*0 --> D0 gamma & 0
2 & D0 --> K- pi+ & 1
0 & B- --> D0 rho- & -1 & B2D0rhom
1 & D0 --> K- pi+ & 0
0 & anti-B0 --> D+ pi- & -1 & B2Dppim
1 & D+ --> K- pi+ pi+ & 0
0 & anti-B0 --> D+ K- & -1 & B2DpKm
1 & D+ --> K- pi+ pi+ & 0
0 & anti-B0 --> D+ rho- & -1 & B2Dprhom
1 & D+ --> K- pi+ pi+ & 0
0 & anti-B0 --> D*+ pi- & -1 & B2Dstppim_pi0
1 & D*+ --> D+ pi0 & 0
2 & D+ --> K- pi+ pi+ & 1
0 & anti-B0 --> D*+ pi- & -1 & B2Dstppim_gamma
1 & D*+ --> D+ gamma & 0
2 & D+ --> K- pi+ pi+ & 1
0 & anti-B0 --> D*+ K- & -1 & B2DstpKm_pi0
1 & D*+ --> D+ pi0 & 0
2 & D+ --> K- pi+ pi+ & 1
0 & anti-B0 --> D*+ K- & -1 & B2DstpKm_gamma
1 & D*+ --> D+ gamma & 0
2 & D+ --> K- pi+ pi+ & 1
0 & anti-B0 --> D+ K*- & -1 & B2DpKstm
1 & D+ --> K- pi+ pi+ & 0
0 & anti-B0 --> D*+ pi- & -1 & B2Dstppim_piplus
1 & D*+ --> D0 pi+ & 0
2 & D0 --> K- pi+ & 1
0 & anti-B0 --> D*+ K- & -1 & B2DstpKm_piplus
1 & D*+ --> D0 pi+ & 0
2 & D0 --> K- pi+ & 1
}
% Main name of output files (Default: Main name of the card file)
{
}
README.md
+ 24
− 1
View file @ 7165dd7a
Measurement of the Ratio of Branching Fractions of B0->D+pi- and B0->D*+pi-
\ No newline at end of file
Measurement of the Ratio of Branching Fractions of B0->D+pi- and B0->D*+pi-
Analysis Directories:
Ntuplizers: /home/belle2/sdey/BRratio/
Ntuples Location: /group/belle2/users2022/sdey/BRratio
After ntuple production jobs are done, we apply topoana module on charged and mixed samples.
Topoana installation instructions:
git clone git@gitlab.desy.de:belle2/software/topoana.git
cd topoana
./Configure
make
./Setup Belle_II
then add path to bash_profile
currently, the topoana is installed in /group/belle2/users2022/sdey/BRratio/mytopoana
example script: topoana.sh
At the end, the outputs are: /group/belle2/users2022/sdey/BRratio/MC15ri/*.root
\ No newline at end of file
continuumsupression/README.txt~ deleted 100644 → 0
+ 0
− 3
View file @ 68cf2c9a
We borrowed these continuum supression scripts from soft pi0 efficiency analysis
After completion of all the jobs. We take equal number of continuum events and non-continuum events for training purpose. The truncation of these root files according to their cross section is taken care of by
continuumsupression/checkjobs.sh~ deleted 100644 → 0
+ 0
− 179
View file @ 68cf2c9a
echo "sdey_MC15ri_b_uubar1" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_uubar1 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_uubar2" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_uubar2 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_uubar3" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_uubar3 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_uubar4" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_uubar4 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_uubar5" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_uubar5 | tail -2 >> jobstatus.txt
#echo "sdey_MC15ri_b_uubar6" >> jobstatus.txt
#gb2_job_status -p sdey_MC15ri_b_uubar6 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_ddbar1" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_ddbar1 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_ddbar2" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_ddbar2 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_ddbar3" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_ddbar3 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_ddbar4" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_ddbar4 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_ddbar5" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_ddbar5 | tail -2 >> jobstatus.txt
#echo "sdey_MC15ri_b_ddbar6" >> jobstatus.txt
#gb2_job_status -p sdey_MC15ri_b_ddbar6 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_ssbar1" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_ssbar1 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_ssbar2" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_ssbar2 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_ssbar3" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_ssbar3 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_ssbar4" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_ssbar4 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_ssbar5" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_ssbar5 | tail -2 >> jobstatus.txt
#echo "sdey_MC15ri_b_ssbar6" >> jobstatus.txt
#gb2_job_status -p sdey_MC15ri_b_ssbar6 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_ccbar1" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_ccbar1 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_ccbar2" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_ccbar2 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_ccbar3" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_ccbar3 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_ccbar4" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_ccbar4 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_ccbar5" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_ccbar5 | tail -2 >> jobstatus.txt
#echo "sdey_MC15ri_b_ccbar6" >> jobstatus.txt
#gb2_job_status -p sdey_MC15ri_b_ccbar6 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_charged1" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_charged1 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_charged2" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_charged2 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_charged3" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_charged3 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_charged4" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_charged4 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_charged5" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_charged5 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_charged6" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_charged6 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_charged7" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_charged7 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_charged8" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_charged8 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_charged9" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_charged9 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_charged10" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_charged10 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_charged11" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_charged11 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_charged12" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_charged12 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_charged13" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_charged13 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_charged14" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_charged14 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_charged15" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_charged15 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_mixed1" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_mixed1 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_mixed2" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_mixed2 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_mixed3" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_mixed3 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_mixed4" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_mixed4 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_mixed5" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_mixed5 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_mixed6" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_mixed6 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_mixed7" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_mixed7 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_mixed8" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_mixed8 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_mixed9" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_mixed9 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_mixed10" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_mixed10 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_mixed11" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_mixed11 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_mixed12" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_mixed12 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_mixed13" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_mixed13 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_mixed14" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_mixed14 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_mixed15" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_mixed15 | tail -2 >> jobstatus.txt
echo "sdey_MC15ri_b_mixed16" >> jobstatus.txt
gb2_job_status -p sdey_MC15ri_b_mixed16 | tail -2 >> jobstatus.txt
continuumsupression/mytestntuplemaker_mc.py~ deleted 100644 → 0
+ 0
− 247
View file @ 68cf2c9a
#!/usr/bin/env python3
# -*- coding: utf-8 -*-
import basf2 as b2
import modularAnalysis as ma
from variables import variables
from stdPhotons import stdPhotons
from modularAnalysis import fillParticleList
from modularAnalysis import matchMCTruth
from modularAnalysis import fillParticleListFromMC
from variables.MCGenTopo import mc_gen_topo
import stdCharged as stdc
import stdV0s as stdv0
import vertex as vx
import variables.utils as vu
import variables.collections as vc
import sys
import os
from vertex import kFit
# create path
my_path = b2.create_path()
ma.inputMdstList('default', '', path=my_path)
ma.fillParticleList(decayString='pi+:goodtracks',
cut='pt> 0.1',
path=my_path)
ma.fillParticleList(decayString='gamma:goodclusters',
cut='E > 0.1',
path=my_path)
# Second, create event shape
ma.buildEventShape(inputListNames=['pi+:goodtracks', 'gamma:goodclusters'],
allMoments=True,
foxWolfram=True,
harmonicMoments=True,
cleoCones=True,
thrust=True,
collisionAxis=True,
jets=True,
sphericity=True,
checkForDuplicates=False,
path=my_path)
# Apply a selection at the event level, to avoid
# processing useless events
ma.applyEventCuts(cut='foxWolframR2 < 0.4', path=my_path)
variables.addAlias("CMclusterE", "useCMSFrame(clusterE)")
variables.addAlias('CMpx', 'useCMSFrame(px)')
variables.addAlias('CMpy', 'useCMSFrame(py)')
variables.addAlias('CMpz', 'useCMSFrame(pz)')
variables.addAlias('CMpt', 'useCMSFrame(pt)')
variables.addAlias('CMp', 'useCMSFrame(p)')
variables.addAlias('CME', 'useCMSFrame(E)')
stdc.stdPi('all', path=my_path)
stdc.stdK('all', path=my_path)
# Calling standard particle lists
ma.cutAndCopyList('pi+:slow', 'pi+:all', 'abs(d0) < 0.5 and abs(z0) < 2', path=my_path)
ma.cutAndCopyList('pi-:slow', 'pi-:all', 'abs(d0) < 0.5 and abs(z0) < 2', path=my_path)
ma.cutAndCopyList('pi+:sig', 'pi+:all', 'abs(d0) < 0.5 and abs(z0) < 2 and thetaInCDCAcceptance and nCDCHits>0', path=my_path)
ma.cutAndCopyList('pi-:sig', 'pi-:all', 'abs(d0) < 0.5 and abs(z0) < 2 and thetaInCDCAcceptance and nCDCHits>0', path=my_path)
ma.cutAndCopyList('K+:sig', 'K+:all', 'abs(d0) < 0.5 and abs(z0) < 2 and thetaInCDCAcceptance and nCDCHits>0', path=my_path)
ma.cutAndCopyList('K-:sig', 'K-:all', 'abs(d0) < 0.5 and abs(z0) < 2 and thetaInCDCAcceptance and nCDCHits>0', path=my_path)
ma.matchMCTruth('pi+:sig', path=my_path)
ma.matchMCTruth('pi-:sig', path=my_path)
ma.matchMCTruth('K+:sig', path=my_path)
ma.matchMCTruth('K-:sig', path=my_path)
ma.matchMCTruth('pi-:slow', path=my_path)
ma.matchMCTruth('pi+:slow', path=my_path)
ma.reconstructDecay('D0:kpi -> K-:sig pi+:sig', '1.8 < M < 1.95', path=my_path)
ma.matchMCTruth('D0:kpi', path=my_path)
ma.reconstructDecay('D+:kpipi -> K-:sig pi+:sig pi+:sig', '1.8 < M < 1.95', path=my_path)
ma.matchMCTruth('D+:kpipi', path=my_path)
#ma.reconstructDecay('D*+:kpi -> D0:kpi pi+:slow', '0.135<massDifference(0)<0.16', path=my_path)
ma.reconstructDecay('D*+:kpi -> D0:kpi pi+:slow', '', path=my_path)
ma.matchMCTruth('D*+:kpi', path=my_path)
ma.reconstructDecay('anti-B0:Dstpi -> D*+:kpi pi-:sig', 'Mbc > 5.2 and abs(deltaE) < 0.5', path=my_path)
ma.matchMCTruth('anti-B0:Dstpi', path=my_path)
ma.reconstructDecay('anti-B0:Dpluspi -> D+:kpipi pi-:sig', 'Mbc > 5.2 and abs(deltaE) < 0.5', path=my_path)
ma.matchMCTruth('anti-B0:Dpluspi', path=my_path)
#ma.reconstructDecay('B-:Dpi -> D0:kpi pi-:sig', 'Mbc > 5.2 and abs(deltaE) < 0.5', path=my_path)
#ma.matchMCTruth('B-:Dpi', path=my_path)
#building continuum supression
cleanMask = ('cleanMask', 'p >= 0.05 and useCMSFrame(p)<=3.5','')
#ma.buildRestOfEvent('anti-B0:Dstpi',path=my_path)
ma.buildRestOfEvent('anti-B0:Dpluspi',path=my_path)
#ma.buildRestOfEvent('B-:Dpi',path=my_path)
ma.appendROEMasks(list_name='anti-B0:Dpluspi', mask_tuples=[cleanMask], path=my_path)
ma.buildContinuumSuppression(list_name='anti-B0:Dpluspi', roe_mask='cleanMask', path=my_path)
cont_suppr_base = [
'R2',
'thrustBm',
'thrustOm',
'cosTBTO',
'cosTBz'
]
#: Transverse energy and missing momentum squared (calculated within the KSFW moments)
et_mm2 = [
'KSFWVariables(et)',
'KSFWVariables(mm2)'
]
#: KSFW Moments calculated using primary particles
ksfw_moments = [
'KSFWVariables(hso10)',
'KSFWVariables(hso12)'
]
#: All available continuum suppression variables
cont_suppr_all = cont_suppr_base + ksfw_moments
trainVars = cont_suppr_all
# Save target variable necessary for training.
targetVar = ['isNotContinuumEvent']
vextra = ['Mbc', 'deltaE']
vids = ['kaonID','pionID']
vkin = ["px","py","pz","pt","InvM","genMotherPDG(0)","genMotherPDG(1)","mcPDG"]
cms_kinematics = ['CMpx', 'CMpy', 'CMpz', 'CMpt', 'CMp', 'CME']
v4mom = ["px","py","pz","pt","E"]
vmcgen = ['genMotherID(0)','genMotherPDG(0)','genNStepsToDaughter(0)', 'genMotherID(1)','genMotherPDG(1)','genNStepsToDaughter(1)', 'genMotherID(2)','genMotherPDG(2)','genNStepsToDaughter(2)','genMotherID(3)','genMotherPDG(3)','genNStepsToDaughter(3)','genMotherID(4)','genMotherPDG(4)','genNStepsToDaughter(4)','genMotherID(5)','genMotherPDG(5)','genNStepsToDaughter(5)','genParticleID']
V_dpluspi = vu.create_aliases_for_selected(list_of_variables=vkin,
decay_string='anti-B0:Dpluspi -> [^D+:kpipi -> K-:sig pi+:sig pi+:sig] ^pi-:sig') + \
vu.create_aliases_for_selected(list_of_variables=vids,
decay_string='anti-B0:Dpluspi -> [D+:kpipi -> ^K-:sig ^pi+:sig ^pi+:sig] ^pi-:sig') + \
vu.create_aliases_for_selected(list_of_variables=vc.mc_truth,
decay_string='^anti-B0:Dpluspi -> [^D+:kpipi -> ^K-:sig ^pi+:sig ^pi+:sig] ^pi-:sig')
ma.variablesToNtuple(decayString='anti-B0:Dpluspi', variables=V_dpluspi + trainVars + targetVar + vextra + mc_gen_topo(160), filename='output_MC15.root', treename='dpluspi', path=my_path)
#----------------------- soft pi0 dependent things
mylist = ['eff60_May2020Fit','eff50_May2020Fit','eff40_May2020Fit','eff30_May2020Fit','eff20_May2020Fit','eff10_May2020Fit']
ma.fillParticleList('gamma:eff60_May2020','[[clusterNHits>1.5] and [0.2967< clusterTheta<2.6180]] and [[clusterReg==1 and E>0.0225] or [clusterReg==2 and E>0.020] or [clusterReg==3 and E>0.020]]', path=my_path)
ma.reconstructDecay('pi0:eff60_May2020 -> gamma:eff60_May2020 gamma:eff60_May2020', '0.03<InvM', path=my_path)
ma.cutAndCopyList('pi0:eff60_May2020Fit', 'pi0:eff60_May2020', '', True, my_path)
kFit('pi0:eff60_May2020Fit', 0.0, 'mass', path=my_path)
ma.matchMCTruth('pi0:eff60_May2020Fit', my_path)
ma.fillParticleList('gamma:eff50_May2020','[[clusterNHits>1.5] and [0.2967< clusterTheta<2.6180]] and [[clusterReg==1 and E>0.025] or [clusterReg==2 and E>0.025] or [clusterReg==3 and E>0.040]]', path=my_path)
ma.reconstructDecay('pi0:eff50_May2020 -> gamma:eff50_May2020 gamma:eff50_May2020', '0.105<InvM<0.150', path=my_path)
ma.cutAndCopyList('pi0:eff50_May2020Fit', 'pi0:eff50_May2020', '', True, my_path)
kFit('pi0:eff50_May2020Fit', 0.0, 'mass', path=my_path)
ma.matchMCTruth('pi0:eff50_May2020Fit', my_path)
ma.fillParticleList('gamma:eff40_May2020','[[clusterNHits>1.5] and [0.2967< clusterTheta<2.6180]] and [[clusterReg==1 and E>0.080] or [clusterReg==2 and E>0.030] or [clusterReg==3 and E>0.060]]', path=my_path)
ma.reconstructDecay('pi0:eff40_May2020 -> gamma:eff40_May2020 gamma:eff40_May2020', '0.120<InvM<0.145', path=my_path)
ma.cutAndCopyList('pi0:eff40_May2020Fit', 'pi0:eff40_May2020', '', True, my_path)
kFit('pi0:eff40_May2020Fit', 0.0, 'mass', path=my_path)
ma.matchMCTruth('pi0:eff40_May2020Fit', my_path)
ma.fillParticleList('gamma:eff30_May2020','[[clusterNHits>1.5] and [0.2967< clusterTheta<2.6180]] and [[clusterReg==1 and E>0.080] or [clusterReg==2 and E>0.030] or [clusterReg==3 and E>0.060 ]] ', path=my_path)
ma.reconstructDecay('pi0:eff30_May2020 -> gamma:eff30_May2020 gamma:eff30_May2020', '0.120<InvM<0.145 and -1.5<daughterDiffOfPhi(0,1)<1.5 and daughterAngle(0,1)<1.4', path=my_path)
ma.cutAndCopyList('pi0:eff30_May2020Fit', 'pi0:eff30_May2020', '', True, my_path)
kFit('pi0:eff30_May2020Fit', 0.0, 'mass', path=my_path)
ma.matchMCTruth('pi0:eff30_May2020Fit', my_path)
ma.fillParticleList('gamma:eff20_May2020','[[clusterNHits>1.5] and [0.2967< clusterTheta<2.6180]] and [[clusterReg==1 and E>0.120] or [clusterReg==2 and E>0.030] or [clusterReg==3 and E>0.080 ]] and clusterE1E9\
>0.4', path=my_path)
ma.reconstructDecay('pi0:eff20_May2020 -> gamma:eff20_May2020 gamma:eff20_May2020', '0.121<InvM<0.142 and -1.0<daughterDiffOfPhi(0,1)<1.0 and daughterAngle(0,1)<0.9', path=my_path)
ma.cutAndCopyList('pi0:eff20_May2020Fit', 'pi0:eff20_May2020', '', True, my_path)
kFit('pi0:eff20_May2020Fit', 0.0, 'mass', path=my_path)
ma.matchMCTruth('pi0:eff20_May2020Fit', my_path)
ma.fillParticleList('gamma:eff10_May2020','[[clusterNHits>1.5] and [0.2967< clusterTheta<2.6180]] and [[clusterReg==1 and E>0.200] or [clusterReg==2 and E>0.100] or [clusterReg==3 and E>0.180]] and clusterE1E9>\
0.5', path=my_path)
ma.reconstructDecay('pi0:eff10_May2020 -> gamma:eff10_May2020 gamma:eff10_May2020', '0.127<InvM<0.139 and -0.9<daughterDiffOfPhi(0,1)<0.9 and daughterAngle(0,1)<0.8', path=my_path)
ma.cutAndCopyList('pi0:eff10_May2020Fit', 'pi0:eff10_May2020', '', True, my_path)
kFit('pi0:eff10_May2020Fit', 0.0, 'mass', path=my_path)
ma.matchMCTruth('pi0:eff10_May2020Fit', my_path)
for element in mylist:
ma.reconstructDecay('D*+:dpi0'+element+' -> D+:kpipi pi0:'+element, '', path=my_path)
ma.matchMCTruth('D*+:dpi0'+element, path=my_path)
ma.reconstructDecay('D*0:kpi'+element+' -> D0:kpi pi0:'+element, '', path=my_path)
ma.matchMCTruth('D*0:kpi'+element, path=my_path)
ma.reconstructDecay('B-:Dst0pi'+element+' -> D*0:kpi'+element+' pi-:sig', 'Mbc > 5.2 and abs(deltaE) < 0.5', path=my_path)
ma.matchMCTruth('B-:Dst0pi'+element, path=my_path)
ma.reconstructDecay('anti-B0:Dstpi_ch2'+element+' -> D*+:dpi0'+element+' pi-:sig', 'Mbc > 5.2 and abs(deltaE) < 0.5', path=my_path)
ma.matchMCTruth('anti-B0:Dstpi_ch2'+element, path=my_path)
ma.buildRestOfEvent('anti-B0:Dstpi_ch2'+element,path=my_path)
V_dstpi_ch2 = vu.create_aliases_for_selected(list_of_variables=["isSignal","mcErrors"],
decay_string='^anti-B0:Dstpi_ch2'+element+' -> [D*+:dpi0'+element+' -> [D+:kpipi -> K-:sig pi+:sig pi+:sig] pi0:'+element+'] pi-:sig') + \
vu.create_aliases_for_selected(list_of_variables=["InvM","genMotherPDG","isSignal","mcErrors","M"],
decay_string='anti-B0:Dstpi_ch2'+element+' -> [^D*+:dpi0'+element+' -> [D+:kpipi -> K-:sig pi+:sig pi+:sig] pi0:'+element+'] pi-:sig') + \
vu.create_aliases_for_selected(list_of_variables=["InvM","isSignal","mcErrors","M","genMotherPDG(0)","genMotherPDG(1)","px","py","pz","pt","genMotherID(1)"],
decay_string='anti-B0:Dstpi_ch2'+element+' -> [D*+:dpi0'+element+' -> [^D+:kpipi -> K-:sig pi+:sig pi+:sig] pi0:'+element+'] pi-:sig') + \
vu.create_aliases_for_selected(list_of_variables=vids,
decay_string='anti-B0:Dstpi_ch2'+element+' -> [D*+:dpi0'+element+' -> [D+:kpipi -> ^K-:sig ^pi+:sig ^pi+:sig] pi0:'+element+'] pi-:sig') + \
vu.create_aliases_for_selected(list_of_variables=cms_kinematics +["mcPDG","InvM","px","py","pz","pt","p","E","isSignal","genMotherPDG(0)","genMotherPDG(1)","mcE"],
decay_string='anti-B0:Dstpi_ch2'+element+' -> [D*+:dpi0'+element+' -> [D+:kpipi -> K-:sig pi+:sig pi+:sig] ^pi0:'+element+'] pi-:sig') + \
vu.create_aliases_for_selected(list_of_variables=vids+["px","py","pz","pt","genMotherPDG(0)","genMotherPDG(1)","mcPDG","genMotherID(0)"],
decay_string='anti-B0:Dstpi_ch2'+element+' -> [D*+:dpi0'+element+' -> [D+:kpipi -> K-:sig pi+:sig pi+:sig] pi0:'+element+'] ^pi-:sig')
ma.variablesToNtuple(decayString='anti-B0:Dstpi_ch2'+element, variables=V_dstpi_ch2 + vextra, filename='output_MC15.root', treename='dstpi_ch2'+element, path=my_path)
# Process the events
b2.process(my_path)
# print out the summary
print(b2.statistics)
topoana.sh 0 → 100644
+ 64
− 0
View file @ 7165dd7a
hadd MC15ri_mixed_part2.root MC15ri_mixed_part2/sub0*/*.root
hadd MC15ri_mixed_part3.root MC15ri_mixed_part3/sub0*/*.root
hadd MC15ri_mixed_part4.root MC15ri_mixed_part4/sub0*/*.root
hadd MC15ri_mixed_part5.root MC15ri_mixed_part5/sub0*/*.root
hadd MC15ri_mixed_part6.root MC15ri_mixed_part6/sub0*/*.root
hadd MC15ri_mixed_part7.root MC15ri_mixed_part7/sub0*/*.root
hadd MC15ri_mixed_part8.root MC15ri_mixed_part8/sub0*/*.root
hadd MC15ri_mixed_part9.root MC15ri_mixed_part9/sub0*/*.root
hadd MC15ri_mixed_part10.root MC15ri_mixed_part10/sub0*/*.root
hadd MC15ri_mixed_part11.root MC15ri_mixed_part11/sub0*/*.root
hadd MC15ri_mixed_part12.root MC15ri_mixed_part12/sub0*/*.root
hadd MC15ri_mixed_part13.root MC15ri_mixed_part13/sub0*/*.root
hadd MC15ri_mixed_part14.root MC15ri_mixed_part14/sub0*/*.root
hadd MC15ri_mixed_part15.root MC15ri_mixed_part15/sub0*/*.root
hadd MC15ri_mixed_part16.root MC15ri_mixed_part16/sub0*/*.root
topoana.exe PartialReco_Dpluspi.card -i MC15ri_mixed_part2.root -o mixed2Partial
topoana.exe PartialReco_Dpluspi.card -i MC15ri_mixed_part3.root -o mixed3Partial
topoana.exe PartialReco_Dpluspi.card -i MC15ri_mixed_part4.root -o mixed4Partial
topoana.exe PartialReco_Dpluspi.card -i MC15ri_mixed_part5.root -o mixed5Partial
topoana.exe PartialReco_Dpluspi.card -i MC15ri_mixed_part6.root -o mixed6Partial
topoana.exe PartialReco_Dpluspi.card -i MC15ri_mixed_part7.root -o mixed7Partial
topoana.exe PartialReco_Dpluspi.card -i MC15ri_mixed_part8.root -o mixed8Partial
topoana.exe PartialReco_Dpluspi.card -i MC15ri_mixed_part9.root -o mixed9Partial
topoana.exe PartialReco_Dpluspi.card -i MC15ri_mixed_part10.root -o mixed10Partial
topoana.exe PartialReco_Dpluspi.card -i MC15ri_mixed_part11.root -o mixed11Partial
topoana.exe PartialReco_Dpluspi.card -i MC15ri_mixed_part12.root -o mixed12Partial
topoana.exe PartialReco_Dpluspi.card -i MC15ri_mixed_part13.root -o mixed13Partial
topoana.exe PartialReco_Dpluspi.card -i MC15ri_mixed_part14.root -o mixed14Partial
topoana.exe PartialReco_Dpluspi.card -i MC15ri_mixed_part15.root -o mixed15Partial
topoana.exe PartialReco_Dpluspi.card -i MC15ri_mixed_part16.root -o mixed16Partial
rm -f MC15ri_mixed_part2.root
rm -f MC15ri_mixed_part3.root
rm -f MC15ri_mixed_part4.root
rm -f MC15ri_mixed_part5.root
rm -f MC15ri_mixed_part6.root
rm -f MC15ri_mixed_part7.root
rm -f MC15ri_mixed_part8.root
rm -f MC15ri_mixed_part9.root
rm -f MC15ri_mixed_part10.root
rm -f MC15ri_mixed_part11.root
rm -f MC15ri_mixed_part12.root
rm -f MC15ri_mixed_part13.root
rm -f MC15ri_mixed_part14.root
rm -f MC15ri_mixed_part15.root
rm -f MC15ri_mixed_part16.root
mv mixed2Partial.root MC15ri_mixed_part2.root
mv mixed3Partial.root MC15ri_mixed_part3.root
mv mixed4Partial.root MC15ri_mixed_part4.root
mv mixed5Partial.root MC15ri_mixed_part5.root
mv mixed6Partial.root MC15ri_mixed_part6.root
mv mixed7Partial.root MC15ri_mixed_part7.root
mv mixed8Partial.root MC15ri_mixed_part8.root
mv mixed9Partial.root MC15ri_mixed_part9.root
mv mixed10Partial.root MC15ri_mixed_part10.root
mv mixed11Partial.root MC15ri_mixed_part11.root
mv mixed12Partial.root MC15ri_mixed_part12.root
mv mixed13Partial.root MC15ri_mixed_part13.root
mv mixed14Partial.root MC15ri_mixed_part14.root
mv mixed15Partial.root MC15ri_mixed_part15.root
mv mixed16Partial.root MC15ri_mixed_part16.root
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment

imprint | privacy | Conditions of Use