Newer
Older
Presentations / tau23mu_lhcb / 20_08_2013 / group_meeting.tex
@mchrzasz mchrzasz on 10 Oct 2013 12 KB update before changing laptops
% This program can be redistributed and/or modified under the terms
% of the GNU Public License, version 3.
%
% Seth Brown, Ph.D.
% sethbrown@drbunsen.org
%
% Compiled with XeLaTeX
% Dependencies:
%   Fontin Sans font (http://www.exljbris.com/fontinsans.fsrtml)
%
\documentclass{beamer}
\usepackage{pgf}
\usepackage{tikz}
\usepackage{times}
\usepackage[T1]{fontenc}
%\usepackage[mathscr]{eucal}
%\usepackage{mathptmx}
%\usepackage{mathrsfs}
\usepackage{hyperref}
\usepackage{color}
\usepackage{graphicx}
\usepackage{wasysym}
% \usepackage{pgfpagfes}
% \setbeameroption{show notes on second screen}
%\pgfpagesuselayout{4 on 1}[a4paper,border shrink=5mm,landscape]
\usepackage{xcolor}
\usepackage{xcolor,multirow}
%\usepackage[table]{xcolor}
%\usepackage[dvipsnames]{xcolor}
%\usepackage{amsfonts}
%\usepackage{amsmath}
% \usepackage[amssymb]{SIunits}
%\usepackage{natbib}
%\usepackage{amssymb}
\usepackage{hepparticles}
\usepackage{hepnicenames}
\usepackage{hepunits}
\usepackage{tikz}
\usepackage[english]{babel}
%%\usepackage{lmodern}
%\usepackage{feynmp}
% suppress navigation bar
\beamertemplatenavigationsymbolsempty
\usepackage[mathscr]{eucal}
\usepackage{mathrsfs}
\mode<presentation>
{
  \usetheme{bunsen}
  \setbeamercovered{transparent}
  \setbeamertemplate{items}[circle]
}
\newcommand{\Simley}[1]{%
\begin{tikzpicture}[scale=0.15]
    \newcommand*{\SmileyRadius}{1.0}%
    \draw [fill=brown!10] (0,0) circle (\SmileyRadius)% outside circle
        %node [yshift=-0.22*\SmileyRadius cm] {\tiny #1}% uncomment this to see the smile factor
        ;  

    \pgfmathsetmacro{\eyeX}{0.5*\SmileyRadius*cos(30)}
    \pgfmathsetmacro{\eyeY}{0.5*\SmileyRadius*sin(30)}
    \draw [fill=cyan,draw=none] (\eyeX,\eyeY) circle (0.15cm);
    \draw [fill=cyan,draw=none] (-\eyeX,\eyeY) circle (0.15cm);

    \pgfmathsetmacro{\xScale}{2*\eyeX/180}
    \pgfmathsetmacro{\yScale}{1.0*\eyeY}
    \draw[color=red, domain=-\eyeX:\eyeX]   
        plot ({\x},{
            -0.1+#1*0.15 % shift the smiley as smile decreases
            -#1*1.75*\yScale*(sin((\x+\eyeX)/\xScale))-\eyeY});
\end{tikzpicture}%
}%
% set fonts
\usepackage{amsfonts}
\usepackage{amsmath}
\usepackage{verbatim}

\usepackage{fancyvrb}
\DefineVerbatimEnvironment{code}{Verbatim}{fontsize=\small}
\DefineVerbatimEnvironment{example}{Verbatim}{fontsize=\small}
\usepackage{listings}
  \usepackage{courier}
 \lstset{
         basicstyle=\footnotesize\ttfamily, % Standardschrift
         %numbers=left,               % Ort der Zeilennummern
         numberstyle=\tiny,          % Stil der Zeilennummern
         %stepnumber=2,               % Abstand zwischen den Zeilennummern
         numbersep=5pt,              % Abstand der Nummern zum Text
         tabsize=2,                  % Groesse von Tabs
         extendedchars=true,         %
         breaklines=true,            % Zeilen werden Umgebrochen
         keywordstyle=\color{red},
    		frame=b,         
 %        keywordstyle=[1]\textbf,    % Stil der Keywords
 %        keywordstyle=[2]\textbf,    %
 %        keywordstyle=[3]\textbf,    %
 %        keywordstyle=[4]\textbf,   \sqrt{\sqrt{}} %
         stringstyle=\color{white}\ttfamily, % Farbe der String
         showspaces=false,           % Leerzeichen anzeigen ?
         showtabs=false,             % Tabs anzeigen ?
         xleftmargin=17pt,
         framexleftmargin=17pt,
         framexrightmargin=5pt,
         framexbottommargin=4pt,
         %backgroundcolor=\color{lightgray},
         showstringspaces=false      % Leerzeichen in Strings anzeigen ?        
 }

    %\DeclareCaptionFont{blue}{\color{blue}} 

  %\captionsetup[lstlisting]{singlelinecheck=false, labelfont={blue}, textfont={blue}}
  \usepackage{caption}
\DeclareCaptionFont{white}{\color{white}}
\DeclareCaptionFormat{listing}{\colorbox[cmyk]{0.43, 0.35, 0.35,0.01}{\parbox{\textwidth}{\hspace{15pt}#1#2#3}}}
\captionsetup[lstlisting]{format=listing,labelfont=white,textfont=white, singlelinecheck=false, margin=0pt, font={bf,footnotesize}}
\usetikzlibrary{arrows}
\usetikzlibrary{shapes}

%\usepackage{gfsartemisia-euler}
%\usepackage[T1]{fontenc}
\setbeamerfont{frametitle}{size=\LARGE,series=\bfseries}
\tikzstyle{decision} = [diamond, draw, fill=gray!20, 
    text width=4.5em, text badly centered, node distance=3cm, inner sep=0pt]
\tikzstyle{block} = [rectangle, draw, fill=blue!10, 
    text width=5em, text centered, rounded corners, minimum height=2em]
\tikzstyle{line} = [draw, -latex']
\tikzstyle{cloud} = [draw, ellipse,fill=red!10, node distance=3cm,
    minimum height=2em]

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


\renewcommand{\PKs}{{\HepParticle{K}{S}{}\xspace}}
% color definitions
\usepackage{color}
\definecolor{uipoppy}{RGB}{225, 64, 5}
\definecolor{uipaleblue}{RGB}{96,123,139}
\definecolor{uiblack}{RGB}{0, 0, 0}

% caption styling
%\DeclareCaptionFont{uiblack}{\color{uiblack}}
%\DeclareCaptionFont{uipoppy}{\color{uipoppy}}
%\captionsetup{labelfont={uipoppy},textfont=uiblack}

% see the macros.tex file for definitions
\include{macros }

% title slide definition
\title{Optimisation of isolation and binning for $\tau \to 3 \mu$.}
%\subtitle{a bias report}
\author{ Marcin Chrz\k{a}szcz$^{1,2}$ ,  Nicola Serra$^{1}$ }
\institute[UTH, IFJ]
{
 %\begin{tiny}
$ ^1$ University of Zurich , $ ^2$ Institute of Nuclear Physics, Krakow,  
 %\end{tiny}smallsmall
}
  

\date{ \begin{small} $21^{th}$ August 2013 \end{small}}

%--------------------------------------------------------------------
%                           Introduction
%--------------------------------------------------------------------

\begin{document}




\setbeamertemplate{background}
{\includegraphics[width=\paperwidth,height=\paperheight]{frontpage_bg_mine}}
\setbeamertemplate{footline}[default]

\begin{frame}
\vspace{1.1cm}
\begin{columns}
\column{2.75in}
  \titlepage

  \begin{center}
  \includegraphics[height=1.0cm ]{pic/uzh.jpg}
  % \hspace{0.5cm}
  % \includegraphics[height=1.5cm]{pic/babar.jpg}
   \hspace{1cm}
   \includegraphics[height=1.0cm]{pic/ifj.png}
   \hspace{1cm}
   %\includegraphics[height=1.0cm]{pic/SNS.jpg}
   \end{center}      
  \vspace{10cm}
\column{2.0in}
\end{columns}
\end{frame}

%--------------------------------------------------------------------
%                          OUTLINE
%--------------------------------------------------------------------




\section[Outline]{}
\begin{frame}
\tableofcontents
\end{frame}




%-------------------------------------------------------------------
%                          Introduction
%-------------------------------------------------------------------
%
% Set the background for the rest of the slides.
% Insert infoline
\setbeamertemplate{background}
 {\includegraphics[width=\paperwidth,height=\paperheight]{slide_bg}}
\setbeamertemplate{footline}[bunsentheme]

\title{Update on analysis}


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\setbeamertemplate{background}
 {\includegraphics[width=\paperwidth,height=\paperheight]{slide_bg}}
\setbeamertemplate{footline}[bunsentheme]
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%2>%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Binning}
\subsection{Binning optimisation}
\begin{frame}\frametitle{Binning optimisation}
For 2011 we did 2 times one dimensional binning optimisation. This method has disadvantages: 
\begin{itemize}
\item The best bin of one classifier is split to very small pieces by the other optimisation.
\item You end up having best bins in the middle. See Pauls presentation.
\end{itemize}

	\textref {M.Chrz\k{a}szcz, N.Serra 2013}

\end{frame}


\begin{frame}\frametitle{How to fix?}
\only<1>
{
\begin{itemize}
\item Perform a simultaneous optimisation in 2D.
\item Brutal force method is not good enough, cuz number of combination explodes.
\item Use MC methods for optimisation.
\item Takes 3 hours to optimise.
\item Use Punzi FOM instead of Cls method.
\item Use 2011 data to optimise.
\item Apply Fine-tuning. You need a given number of events to perform a stable fit etc.
\end{itemize}
}
\only<2>
{
  \includegraphics[scale=0.23]{inflaton/punzi1.png}

}
	\textref {M.Chrz\k{a}szcz, N.Serra 2013}

\end{frame}








\section{Isolation}
\subsection{Isolation optimisation}

\begin{frame}\frametitle{Iso optimisation}
{~}


\only<1>
{	
	Till now every analysis that used track isolation parameter used the ones develeloped and optimised for $B_s \to \mu \mu$.
	This is based on an abstract definitions of isolating and non-isolating tracks:
	\begin{itemize}
	\item Non-isolating track to a given track($\mu$ from $B_s \to \mu \mu$ for example) will be a track that has the same primary mother as muon. 
	\item Isolating is the negation of non-isolating. 
\end{itemize}	
{~} \\
	\begin{Large}
	Many thanks to Giampi for discussion and advises !
	\end{Large}
	
}

\only<2>
{	
	This definition has potentially dangerous implications.
	  	\begin{itemize}
	  	\item Imagine a long chain of decays. Every of this decay is non-isolating.
	  	\item Why very long living particles ($\Lambda$, $K_s$) have to be considered non-isolating?
	  	\item When we do our analysis we are operating on basis of signal and bck hypothesis.
	  	\item There isnt a 1:1 correspondence between isolating and bck etc.
		\item $B_s \to \mu \mu$ does not really suffer from this. 	
	
	  	
	  		  		\end{itemize}	
	  	


}

	\textref {M.Chrz\k{a}szcz, N.Serra 2013}

\end{frame}



\begin{frame}\frametitle{How to train?}
{~}



\only<1>
{	
	\begin{enumerate}
	
	
	\item The main point of isolation variable is to fight again combinatorial bck.(example two decays trees are close and one picks something from the other).
	\item We build our bck sample taking from MC truth the candidates that are combinatorial bck. 
		
	\end{enumerate}
		
}
\only<2>
{	
	
	Now I will loose you all :P
	\begin{enumerate}
	
	
	\item We need to swap our signal and bck sample.
	\item Why? Our signal sample contains: signal candidate(3 tracks)+ tracks surrounding this candidate. Our selection should be optimised in a way that we should end up with our single signal candidate without any tracks nearby.
	\item That is why our signal sample is our background sample.
				
		
	\end{enumerate}
		
}
\only<3>
{	
	\begin{enumerate}

	\item We define the training variables as Giampi did:+tckchi2+IP.
	\item We put everything inside tmva.
	\item Then we scan the BDT response space and write how many tracks survive the cut.
	\item Optimisation of the cut has to be done inside the BDT that we will use.
	
	
		
	\end{enumerate}
		
}
\only<4>
{	
	\begin{enumerate}

	\item In practice what we do is to scan BDT form 0. to 0.5 and count the tracks for each of the BDT value.
	\item Then our new ntuple will have like 100 isolation parameters.
	\item How to choose the best one?
	\item Well isolation parameter on its own is useless. It has to be combined with other variables in TMVA. Than you can choose the best cut on the BDT.
	
	
		
	\end{enumerate}
		
}

	\textref {M.Chrz\k{a}szcz, N.Serra 2013}

\end{frame}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5555
\begin{frame}\frametitle{$\tau \to 3 \mu$ specifics}
{~}


\only<1>
{	
In case of $\tau \to 3 \mu$ we want different isolation parameters for different kinds of decays:
	\begin{itemize}
	\item $D \to \tau$
	\item $Ds \to \tau$		
	\item $B \to D \to \tau$
	\item $B \to Ds \to \tau$			
	\item $B \to \tau$	
	\end{itemize}
	
	
}
\only<2>
{	
\begin{enumerate}

\item Does it make any sense to make my life so complicated?
\item YES! 
\item Example: $B \to \tau$ is in $99 \%$ $B \to D \tau X$.
\item This means we if you have D and tau close to each other track from D can go to $\tau$ etc. 
\item In their approach this truck would be considered non-isolating which is nonsense because it forms a bck candidate!
\item From first looks the problem can be reduced to 3 chains:  $B \to \tau$, 
 $B \to Dx \to \tau$, $D \to \tau$.
\end{enumerate}
	
	
}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5



\only<3>
{	

	\begin{columns}
\column{2.5in}

  \includegraphics[scale=0.21]{tmva_tau23mu/mva_BDT.png}


\column{2.5in}

 \includegraphics[scale=0.21]{tmva_tau23mu/rejBvsS.png}
\end{columns}


	\begin{columns}
\column{2.5in}

  \includegraphics[scale=0.21]{tmva_tau23mu/mva_BDT_21513000.png}


\column{2.5in}

 \includegraphics[scale=0.21]{tmva_tau23mu/rejBvsS_21513000.png}
\end{columns}





	
}







\end{frame}



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55





\end{document}