\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 \\ Nicola Serra }

\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 9, 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}[c]%{Plans 2}
	\begin{minipage}{\textwidth}

\begin{center}
\begin{Huge}
MC model
\end{Huge}
\end{center}

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


\begin{frame}[c]{Acceptance correction}

	\begin{minipage}{\textwidth}

\begin{columns}

\column{0.6\textwidth}	
	
\ARROW The decay of $\PBzero \to \PKstar \Pmuon \APmuon$ is described by 3 helicity angles and the invariant mass squared of two leptons $(q^2)$.\\
\ARROW In order to model the detector acceptance we have used a large MC sample of PHSP simulated events.\\
\ARROW There is a caveat: the $q^2$ distribution.\\
\ARROW We had to reweight it to make it flat.
\column{0.4\textwidth}	


\includegraphics[width=0.99\textwidth]{images/q2PHSP.png}\\
\includegraphics[width=0.99\textwidth]{images/q2PHSPw.png}
\end{columns}

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


\begin{frame}[c]{Can we optimize it?}

	\begin{minipage}{\textwidth}

\ARROW It would be nice if we could generate not only the flat angle distributions but also a flat $q^2$.\\
\ARROW There exists already a model for it: \texttt{FLATQ2}.\\
\ARROW It basically reweighs the  distribution by $1/p^{\rm had}_T$.\\
\ARROW The problem is that it was design to generate the flat distribution of decays $\PB \to X \ell \nu$ :

\begin{center}
\includegraphics[width=0.7\textwidth]{images/flatq2.png}
\end{center}

\ARROW Will not work in current version for $\PB \to \PKstar \mu \mu$.

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


\begin{frame}[c]{Modifying the \texttt{FLATQ2} 1 }

	\begin{minipage}{\textwidth}

\ARROW I wrote a mirror model that requires that the two leptons are \texttt{DIRAC}, and called it \texttt{FLATQ2EWP}. \\
\ARROW And improves the situation a lot:
\begin{center}
\includegraphics[angle=-90,width=0.8\textwidth]{{images/ForTom_Q2_FLATQ2}.pdf}
\end{center}

\ARROW So much flatter but the and and the begging still not flat.\\
\ARROW End of the spectrum is due to $\PKstar$ width? $\twoheadrightarrow$ Lets test it with $\PB \to \PK \mu \mu$.
 
	\end{minipage}
		\vspace*{2.cm}
\end{frame}

\begin{frame}[c]{Modifying the \texttt{FLATQ2} 1 }

	\begin{minipage}{\textwidth}
{~}\\
\ARROW \texttt{FLATQ2EWP} use to simulate the $\PB  \to \PK \mu \mu$:
\begin{center}
\includegraphics[angle=-90,width=0.75\textwidth]{{images/ForTom_Q2_FLATQ2_Kmumu}.pdf}
\end{center}

\ARROW Oki so end of the spectrum is understood and not much can be done there.\\
\ARROW Now the low $q^2$: Can this be just Phase space suppression: $\sqrt{\lambda}=\sqrt{1-4m_{\ell}^2/q^2}$
	\end{minipage}
		\vspace*{2.cm}
\end{frame}

\begin{frame}[c]{Modifying the \texttt{FLATQ2} 2 }

	\begin{minipage}{\textwidth}
{~}\\
\ARROW \texttt{FLATQ2EWP} with phase space suppression factor. 
\begin{center}
\includegraphics[angle=-90,width=0.75\textwidth]{{images/ForTom_Q2_FLATQ2_Kmumu_lambda}.pdf}
\end{center}
\ARROW Now it's perfect.
	\end{minipage}
		\vspace*{2.cm}
\end{frame}




\begin{frame}[c]{Update since last week }
	\begin{minipage}{\textwidth}
{~}\\
\ARROW Discussion was made via: \href{https://its.cern.ch/jira/browse/LHCBGAUSS-911}{JIRA}\\
\ARROW It was suggested my Michal to incorporate the new model into the current one to save the code.\\
\ARROW Thanks to John for merging the two codes:
\begin{center}
\only<1>{
\includegraphics[width=0.7\textwidth]{images/evt1.png}
}
\only<2>{
\includegraphics[width=0.7\textwidth]{images/evt2.png}
}
\end{center}






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




\begin{frame}[c]{\texttt{FLATQ2} Conclusion }

	\begin{minipage}{\textwidth}
{~}\\
\ARROW The new model was tested by me and John.\\
\ARROW Changes won't have any influence on the existing \texttt{DEFILES} as the flag is by default switched off.\\
\ARROW The commit was merge to master by Gloria today.\\
\ARROW We thank all people involved action\\
{~}\\
\ARROW The whole things took $<$week and is already available for production!\\
\ARROW There is also other model \texttt{XLL}, see Biplab slides more suitable for $\PB \to \PK \Ppi \mu \mu$. 

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





\backupbegin

\begin{frame}\frametitle{Backup}


\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}




\backupend

\end{document}
