Newer
Older
Presentations / Kstarmumu_Run2 / Filter1 / mchrzasz.tex
\documentclass[11 pt,xcolor={dvipsnames,svgnames,x11names,table}]{beamer}

\usepackage[english]{babel}
\usepackage{polski}


\usetheme[
	bullet=circle,		% Other option: square
	bigpagenumber,		% circled page number on lower right
	topline=true,			% colored bar at the top of the frame
	shadow=false,			% Shading for beamer blocks
	watermark=BG_lower,	% png file for the watermark
	]{Flip}

%\logo{\kern+1.em\includegraphics[height=1cm]{SHiP-3_LightCharcoal}}


\usepackage[lf]{berenis}
\usepackage[LY1]{fontenc}
\usepackage[utf8]{inputenc}

\usepackage{emerald}
\usefonttheme{professionalfonts}
\usepackage[no-math]{fontspec}
\defaultfontfeatures{Mapping=tex-text}	% This seems to be important for mapping glyphs properly

\setmainfont{Gillius ADF}			% Beamer ignores "main font" in favor of sans font
\setsansfont{Gillius ADF}			% This is the font that beamer will use by default
% \setmainfont{Gill Sans Light}		% Prettier, but harder to read

\setbeamerfont{title}{family=\fontspec{Gillius ADF}}

\input t1augie.fd

%\newcommand{\handwriting}{\fontspec{augie}} % From Emerald City, free font
%\newcommand{\handwriting}{\usefont{T1}{fau}{m}{n}} % From Emerald City, free font
% \newcommand{\handwriting}{}	% If you prefer no special handwriting font or don't have augie

%% Gill Sans doesn't look very nice when boldfaced
%% This is a hack to use Helvetica instead
%% Usage: \textbf{\forbold some stuff}
%\newcommand{\forbold}{\fontspec{Arial}}

\usepackage{graphicx}
\usepackage[export]{adjustbox}

\usepackage{amsmath}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{bm}
\usepackage{colortbl}
\usepackage{mathrsfs} 			% For Weinberg-esque letters
\usepackage{cancel}				% For "SUSY-breaking" symbol
\usepackage{slashed}            % for slashed characters in math mode
\usepackage{bbm}                % for \mathbbm{1} (unit matrix)
\usepackage{amsthm}				% For theorem environment
\usepackage{multirow}			% For multi row cells in table
\usepackage{arydshln} 			% For dashed lines in arrays and tables
\usepackage{siunitx}
\usepackage{xhfill}
\usepackage{grffile}
\usepackage{textpos}
\usepackage{subfigure}
\usepackage{tikz}

%\usepackage{hepparticles}
\usepackage[italic]{hepparticles}

\usepackage{hepnicenames}

% Drawing a line
\tikzstyle{lw} = [line width=20pt]
\newcommand{\topline}{%
  \tikz[remember picture,overlay] {%
    \draw[crimsonred] ([yshift=-23.5pt]current page.north west)
             -- ([yshift=-23.5pt,xshift=\paperwidth]current page.north west);}}



% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %
\usepackage{tikzfeynman}		% For Feynman diagrams
\usetikzlibrary{arrows,shapes}
\usetikzlibrary{trees}
\usetikzlibrary{matrix,arrows} 				% For commutative diagram
% http://www.felixl.de/commu.pdf
\usetikzlibrary{positioning}				% For "above of=" commands
\usetikzlibrary{calc,through}				% For coordinates
\usetikzlibrary{decorations.pathreplacing}  % For curly braces
% http://www.math.ucla.edu/~getreuer/tikz.html
\usepackage{pgffor}							% For repeating patterns

\usetikzlibrary{decorations.pathmorphing}	% For Feynman Diagrams
\usetikzlibrary{decorations.markings}
\tikzset{
	% >=stealth', %%  Uncomment for more conventional arrows
	vector/.style={decorate, decoration={snake}, draw},
	provector/.style={decorate, decoration={snake,amplitude=2.5pt}, draw},
	antivector/.style={decorate, decoration={snake,amplitude=-2.5pt}, draw},
	fermion/.style={draw=gray, postaction={decorate},
		decoration={markings,mark=at position .55 with {\arrow[draw=gray]{>}}}},
	fermionbar/.style={draw=gray, postaction={decorate},
		decoration={markings,mark=at position .55 with {\arrow[draw=gray]{<}}}},
	fermionnoarrow/.style={draw=gray},
	gluon/.style={decorate, draw=black,
		decoration={coil,amplitude=4pt, segment length=5pt}},
	scalar/.style={dashed,draw=black, postaction={decorate},
		decoration={markings,mark=at position .55 with {\arrow[draw=black]{>}}}},
	scalarbar/.style={dashed,draw=black, postaction={decorate},
		decoration={markings,mark=at position .55 with {\arrow[draw=black]{<}}}},
	scalarnoarrow/.style={dashed,draw=black},
	electron/.style={draw=black, postaction={decorate},
		decoration={markings,mark=at position .55 with {\arrow[draw=black]{>}}}},
	bigvector/.style={decorate, decoration={snake,amplitude=4pt}, draw},
}

% TIKZ - for block diagrams,
% from http://www.texample.net/tikz/examples/control-system-principles/
% \usetikzlibrary{shapes,arrows}
\tikzstyle{block} = [draw, rectangle,
minimum height=3em, minimum width=6em]




\usetikzlibrary{backgrounds}
\usetikzlibrary{mindmap,trees}	% For mind map
\newcommand{\degree}{\ensuremath{^\circ}}
\newcommand{\E}{\mathrm{E}}
\newcommand{\Var}{\mathrm{Var}}
\newcommand{\Cov}{\mathrm{Cov}}
\newcommand\Ts{\rule{0pt}{2.6ex}}       % Top strut
\newcommand\Bs{\rule[-1.2ex]{0pt}{0pt}} % Bottom strut

\graphicspath{{images/}}	% Put all images in this directory. Avoids clutter.

% SOME COMMANDS THAT I FIND HANDY
% \renewcommand{\tilde}{\widetilde} % dinky tildes look silly, dosn't work with fontspec
\newcommand{\comment}[1]{\textcolor{comment}{\footnotesize{#1}\normalsize}} % comment mild
\newcommand{\Comment}[1]{\textcolor{Comment}{\footnotesize{#1}\normalsize}} % comment bold
\newcommand{\COMMENT}[1]{\textcolor{COMMENT}{\footnotesize{#1}\normalsize}} % comment crazy bold
\newcommand{\Alert}[1]{\textcolor{Alert}{#1}} % louder alert
\newcommand{\ALERT}[1]{\textcolor{ALERT}{#1}} % loudest alert
%% "\alert" is already a beamer pre-defined
\newcommand*{\Scale}[2][4]{\scalebox{#1}{$#2$}}%

\def\Put(#1,#2)#3{\leavevmode\makebox(0,0){\put(#1,#2){#3}}}

\usepackage{gmp}
\usepackage[final]{feynmp-auto}

\usepackage[backend=bibtex,style=numeric-comp,firstinits=true]{biblatex}
\bibliography{bib}
\setbeamertemplate{bibliography item}[text]

\makeatletter\let\frametextheight\beamer@frametextheight\makeatother

% suppress frame numbering for backup slides
% you always need the appendix for this!
\newcommand{\backupbegin}{
	\newcounter{framenumberappendix}
	\setcounter{framenumberappendix}{\value{framenumber}}
}
\newcommand{\backupend}{
	\addtocounter{framenumberappendix}{-\value{framenumber}}
	\addtocounter{framenumber}{\value{framenumberappendix}}
}


\definecolor{links}{HTML}{2A1B81}
%\hypersetup{colorlinks,linkcolor=,urlcolor=links}

% For shapo's formulas:
\def\lsi{\raise0.3ex\hbox{$<$\kern-0.75em\raise-1.1ex\hbox{$\sim$}}}
\def\gsi{\raise0.3ex\hbox{$>$\kern-0.75em\raise-1.1ex\hbox{$\sim$}}}
\newcommand{\lsim}{\mathop{\lsi}}
\newcommand{\gsim}{\mathop{\gsi}}
\newcommand{\wt}{\widetilde}
%\newcommand{\ol}{\overline}
\newcommand{\Tr}{\rm{Tr}}
\newcommand{\tr}{\rm{tr}}
\newcommand{\eqn}[1]{&\hspace{-0.7em}#1\hspace{-0.7em}&}
\newcommand{\vev}[1]{\rm{$\langle #1 \rangle$}}
\newcommand{\abs}[1]{\rm{$\left| #1 \right|$}}
\newcommand{\eV}{\rm{eV}}
\newcommand{\keV}{\rm{keV}}
\newcommand{\GeV}{\rm{GeV}}
\newcommand{\im}{\rm{Im}}
\newcommand{\disp}{\displaystyle}
\def\be{\begin{equation}}
\def\ee{\end{equation}}
\def\ba{\begin{eqnarray}}
\def\ea{\end{eqnarray}}
\def\d{\partial}
\def\l{\left(}
\def\r{\right)}
\def\la{\langle}
\def\ra{\rangle}
\def\e{{\rm e}}
\def\Br{{\rm Br}}

\def\ARROW{{\color{JungleGreen}{$\Rrightarrow$}}\xspace}
\def\ARROWR{{\color{WildStrawberry}{$\Rrightarrow$}}\xspace}

\author{ {\fontspec{Trebuchet MS}Marcin Chrz\k{a}szcz} (Universit\"{a}t Z\"{u}rich)}
\institute{UZH}
\title[$\PBzero \to \PKstar \Pmuon \APmuon$ MC Filter]{$\PBzero \to \PKstar \Pmuon \APmuon$ MC Filter}
\date{25 September 2014}


\begin{document}
\tikzstyle{every picture}+=[remember picture]

{
\setbeamertemplate{sidebar right}{\llap{\includegraphics[width=\paperwidth,height=\paperheight]{bubble2}}}
\begin{frame}[c]%{\phantom{title page}}
\begin{center}
\begin{center}
	\begin{columns}
		\begin{column}{0.75\textwidth}
			\flushright\fontspec{Trebuchet MS}\bfseries \Huge {$\PBzero \to \PKstar \Pmuon \APmuon$\\ MC Filter}
		\end{column}
                \begin{column}{0.02\textwidth}
                  {~}
                  \end{column}
                \begin{column}{0.23\textwidth}
                 % \hspace*{-1.cm}
                  \vspace*{-3mm}
                  \includegraphics[width=0.6\textwidth]{lhcb-logo}
                  \end{column}

	\end{columns}
\end{center}
	\quad
	\vspace{3em}
\begin{columns}
\begin{column}{0.44\textwidth}
\flushright \vspace{-1.8em} {\fontspec{Trebuchet MS} \Large Marcin ChrzÄ…szcz\\\vspace{-0.1em}\small \href{mailto:mchrzasz@cern.ch}{mchrzasz@cern.ch}\\ \Large Thomas Blake \\\vspace{-0.1em}\small \href{mailto:thomas.blake@cern.ch}{thomas.blake@cern.ch} }

\end{column}
\begin{column}{0.53\textwidth}
\includegraphics[height=1.3cm]{uzh-transp}
\end{column}
\end{columns}

\vspace{1em}
%		\footnotesize\textcolor{gray}{With N. Serra, B. Storaci\\Thanks to the theory support from M. Shaposhnikov, D. Gorbunov}\normalsize\\
\vspace{0.5em}

	\textcolor{normal text.fg!50!Comment}{$\PBzero \to \PKstar \Pmuon \APmuon$ meeting, CERN\\October 10, 2016}
\end{center}
\end{frame}
}


\begin{frame}[c]{A glimpse in the Run1 analysis}

	\begin{minipage}{\textwidth}
\ARROW In the Run1 we have asked for a filtered MC to correct for detector acceptance.\\
\ARROW Asked for $5.5M$ events (after stripping in DST), which means we generated around $110M$ events. \\
\ARROW After our full selection we ended up with with only $1.4M$ events.\\
\\{~}\\{~}
\begin{alertblock}{Warning!}
The stripping line has a PID cut inside: ${\rm PID_{\mu}}>-3$.\\
This essentially means we model that efficiency from MC.
\end{alertblock}



	\end{minipage}
		\vspace*{2.cm}
\end{frame}



\begin{frame}[c]{Run2 options}
	\begin{minipage}{\textwidth}


\begin{enumerate}
\item Repeat what we did in Run1 and keep the PID cuts.
\item Filter on stripping removing the PID cut.
\item Filter on MC truth: 
\begin{itemize}
\item 4 charge tracks on $\rm{StdAllNoPiDPions/Kaons/Muons}$
\item And truth matched the decay channel: $\rm{mcMatch( '[ B0 => K*(892)0 mu+ mu- ]CC}')$
\end{itemize}
\end{enumerate}

\begin{exampleblock}{Why MCTruth?}
\ARROWR We are using a very old stripping line that for sure can be (and should be) optimized for the final analysis of Run2!\\
\ARROWR Producing an MCTRUTH match sample would allow the sample to be reused for future analysis even if the stripping line will change!
\end{exampleblock}

	\end{minipage}
		\vspace*{2.cm}
\end{frame}



\begin{frame}[c]{Retentions}
	\begin{minipage}{\textwidth}

\ARROW To study the solution I have used 2012 Physics MC.\\
\ARROW I have taken 17.250 simulated events.\\
\ARROW Here is the results:\\{~}\\
\begin{tabular}{|c|c|c|c|}
\hline 
Type & Filter retention & Events in the ntuple & Truth Matched \\ \hline
Strip & $3447~(20~\%)$ & $4975$ & $1648$ \\
Strip no PID$_{\mu}$ & $3504~(20.3~\%)$ & $5176$ & $1660$ \\
MCTruth & $5009~(29~\%)$ & $4456$ & $1660$ \\
\hline
\end{tabular}
{~}\\
\ARROW Now I have cross check this running the same algorithms on stripped and non stripped MC always getting the same numbers.\\
\ARROW For speed purpose I have put a cut on the $m_{\PKstar}<1300~\rm MeV$ (can be adjusted if needed).\\
\ARROW Other option to consider is to remove \texttt{ISMUON} form stripping to get all efficiencies from \texttt{PIDCalib}.



	\end{minipage}
		\vspace*{2.cm}
\end{frame}


\begin{frame}[c]{Plans}
	\begin{minipage}{\textwidth}

\ARROW With Tom we feel that it would be best to ask for $200M$ generated events. \\
\ARROW Also we noticed that we have $50M$ events of some old MC10 (Stripping 12) MC, which we propose to delete.\\
\ARROW For PPG: The $R(\PDstar)$ have already got green light for more then $1000M$ generated events, so we getting the $200M$ should not be a problem. \\
\ARROW To discuss: Do we want a flat $m(K\pi)$ sample or we can keep the $\PKstar$?


	\end{minipage}
		\vspace*{2.cm}
\end{frame}



\begin{frame}[c]{Plans 2}
	\begin{minipage}{\textwidth}
\ARROW Besides the normal $\PB \to \PKstar \Pmu \Pmu$ PHSP we should ask for other MC channels.\\
\ARROW I proposed to scale the old numbers by factor: $\frac{5}{3}$.
{~}\\
\begin{small}

\begin{tabular}{lrrr}\hline                                                                                                                                                                                                                   
Decay & DecFile event type & N. of events & N. of events Run2\\ \hline\hline                                                                                                                                                                               
$\PB \to \PKstar \PJpsi$ (physics) & {\tt 11144001} & 2M & 3.5M \\                                                                                                                                                                             
$\PB \to \PKstar \PJpsi$ (PHSP?) & {\tt xxxxxxxx} & 0 & 3.5M \\                                                                                                                                                                             
 $\PB \to \PKstar \mu \mu$ (physics)                & {\tt 11114001} & 1M & 1.5M\\                                                                                                                                                                             
$\PLambdab \to \PLambda(1530)\mu\mu$          & {\tt 15114000} & 1M & 1.5M\\                                                                                                                                                                             
$\PLambdab \to p K \mu\mu$              & {\tt 15114011} & 2M & 3.5M\\                                                                                                                                                                             
$\PBs \to \phi\mu\mu$              & {\tt 13114002} & 0.6M & 1M\\                                                                                                                                                                           
$\PB_u \to  \PK \mu\mu$               & {\tt 12113001} & 1M & 1.5M\\\hline                                                                                                                                                                       
\hline\end{tabular} 
\end{small}

\ARROW This would be unfiltered production and this MC will be needed for other analysis as well.\\
\ARROW Do we want to simulate a flat $q^2$ in the $\PB \to \PKstar \Pmu \Pmu$?\\
\ARROW Do we want to have a flat $\PK \Ppi$ mass distribution in the simulation?

	\end{minipage}
		\vspace*{2.cm}
\end{frame}
\begin{frame}
\begin{center}
\begin{Huge}
This is not related to MC requests.
\end{Huge}
\end{center}


\end{frame}


\begin{frame}\frametitle{MCmatching studies.}
\ARROW Let's look how the candidates that have been matched by: $\rm{mcMatch( '[ B0 => K*(892)0 mu+ mu- ]CC}')$ look like:

\begin{center}
\includegraphics[angle=-90,width=0.5\textwidth]{{images/BKGCAT_mctruthmatched}.pdf}
{~}\\


\end{center}
\ARROW \texttt{BKGCAT==10} is the pure signal. The mcMatch is not changing anything in that number of entries.\\
\ARROW \texttt{BKGCAT==30} is the K=$\PK \leftrightarrows \pi$ swaps. This goes away with some PID selection

\end{frame}

\begin{frame}\frametitle{MCmatching studies.}

\ARROW Now all \texttt{BKGCAT==10} have true mcMatch:
\begin{center}
\includegraphics[angle=-90,width=0.4\textwidth]{{images/BKGCAT10}.pdf}

\end{center}



\ARROW How does \texttt{BKKCAT==50,40} (missID +FSR, FSR)look like:

\includegraphics[angle=-90,width=0.45\textwidth]{{images/BKGCAT40MCMatch}.pdf}
\includegraphics[angle=-90,width=0.45\textwidth]{{images/BKGCAT50MCMatch}.pdf}\\
\ARROW We need to consider which \texttt{BKGCAT} we should use for the analysis.

\end{frame}

\backupbegin

\begin{frame}\frametitle{Backup}
\topline

\end{frame}

\backupend

\end{document}