-\section{Run by run QA, how to and code}\r
-\subsection{Online - Francesco, Michael}\r
-\r
-\subsubsection{Creation and checking of online QA histograms (AMORE)}\r
-The QA histograms for the EMCal are created and filled by the class:\\ \r
-\texttt{\$ALICE\_ROOT/EMCAL/AliEMCALQADataMakerRec}\\\r
-It is run both during the offline reconstruction and by the online data quality monitoring framework. Its main methods are:\\\r
-\r
-\texttt{InitRaws()}:\r
-All the QA histograms are created here. Their titles should be self-explaining. Each of them has 2 important flags:\\\r
-\texttt{expert} - if true, the histogram is shipped only to the AMORE expert agent (run in the EMCal station), otherwise it is also shipped to the AMORE shifter agent (at the moment we have 4 histograms monitored by the DQM shifter).\\\r
-\texttt{image} - if true, the plots is saved to the loogbok (there should be 9 plots in the logbook).\\\r
-Each of those histograms is replicated 4 times by the amore framework and filled according to the event species (Calib, Cosmic, LowMultiplicity, HighMultiplicity).\\\r
-\r
-\texttt{MakeRaws(AliRawReader* rawReader)}:\r
-Here we loop over the input raw data stream and we fill the histograms.\\\r
-\r
-\texttt{MakeRawsSTU(AliRawReader* rawReader)}:\r
-STU raw data decoding, provided EMCal trigger experts.\\\r
-\r
-\texttt{GetCalibRefFromOCDB()}:\r
-Get the reference histograms from the OCDB. The ratio between histograms from current run and the reference run is calculated in the \texttt{MakeRaws} method.\\\r
-\r
-\r
-The QA histograms are then analyzed by the class:\\\r
-\texttt{\$ALICE\_ROOT/EMCAL/AliEMCALQAChecker}\\\r
-It checks the data contained in the 4 non-expert histograms. They are at the moment: \r
-\begin{enumerate}\r
-\item Ratio distribution of towers' amplitude in the current run w.r.t. the reference run.\r
-\item Number of hits of L1 Gamma patch.\r
-\item Number of hits of L1 Jet patch.\r
-\item Number of links between TRU and STU. \r
-\end{enumerate}\r
-The main method is:\r
-\texttt{CheckRaws()}: For the first histogram, we check how many\r
-channels have the ratio in the region $0.8-1.2$. If there are more than\r
-the threshold (default = 90\%), everything is ok, otherwise a red box with a ``call\r
-expert'' message is displayed. For L1 Jet/Gamma patch, we check if the\r
-rate of a single patch is higher than the average rate of all other\r
-patches times a certain threshold value (default = 0.5):\r
-\begin{equation}\r
-Rate_{patch}/Rate_{Total} > Threshold /(1 + Threshold)\r
-\end{equation}\r
-\r
- If so red box with a ``hot spot - call expert'' message is\r
- displayed. Finally, if there is a number (default = 1) of missing\r
- STU-TRU links, another red box with a warning message is shown.\r
-If one wants to change the thresholds, this can be done by updating\r
-them in the database (type \texttt{amoreConfigFileBrowser d emc} at\r
-any machine at P2 and edit the file \texttt{QAThresholds.configfile} accordingly).\r
-\r
-\subsubsection{How to's for EMCal AMORE experts}\r
-\r
-\paragraph{How to run AMORE at point 2 (P2):}\r
-\r
-Some useful scripts for running AMORE at P2 are located in the \texttt{\textasciitilde/amore} directory in the aldaqacr51 machine. \texttt{StartAmore.sh} is used to run the agent and the GUI (in the expert mode). It simply launches the \texttt{configMonitor.sh} and the \texttt{checkLogs.sh} (used to check the logs of the agents). \texttt{configMonitor.sh} setup the agent and the GUI using some kdialog commands and runs them.\r
-Basically the agent is run by the following line in the script:\\\r
-\texttt{sshdqm \$agentName nohup \$scriptsPATH/startAmoreAgent.sh \$onlineMode \&}\\\r
-The agent is actually run on the EMCal DQM node via the sshdqm command, so this line basically is used to connect to our dqm node and run the agent there.\\\r
-The GUI is launched in our machine by the command:\\\r
-\texttt{amore -c EMC -m EMCUIQA(TRU) -g configFile.txt}\\\r
-The \texttt{configFile.txt} file contains just one line, i.e. \texttt{runType i}, where \texttt {i} runs from 1 to 4 (1=LowMultiplicity, 2=HighMultiplicity, 3=Cosmic, 4=Calib) and it is used to select the event specie you want to display in the GUI (this file is automatically created each time you run \texttt{configMonitor.sh}).\\\r
-Some other useful scripts (which are launched by \texttt{configMonitor.sh} using the sshdqm) are located in our DQM node. One can login to it by using:\\\r
-\texttt{ssdqm EMCQA}\\\r
-They are located in the \texttt{scripts} directory.\\\r
-\texttt{startAmoreAgent.sh} basically runs the agent, it simply issues the command:\\\r
-\texttt{amoreAgent -a EMCQA -s @\$1: -u}\\\r
-The only parameter here is \$1, which is the gdc to monitor, defined in \texttt{configMonitor.sh}. \r
-\texttt{killAgent.sh} kills any running agent. Be aware that there cannot be 2 running agents of the same kind at the same time (i.e. just one EMCQA and one EMCQAshifter at a time). It can happen that some time the agent does not start when you try to run it. Most of the time this is because the dqm shifter is already running the EMCQA agent instead of just the EMCQAshifter one (they can run it from the dqm station even if they shouldn't). To check if the dqm shifter is running it, simply do \texttt{ps aux | grep EMCQA} in the dmq node, and see if there is an EMCQA process belonging to the user daq. If so, kindly ask the dqm shifter to please kill it.\\\r
-The agent is always running using the same parameters, which are\r
-stored in a database. In order to change them, you should run the\r
-\texttt{amoreConfigFileBrowser} command in the aldaqacr51 machine. A\r
-window will appear, there you can browse a lot of configuration\r
-files. Our files are \texttt{EMCAL\_config.txt} (it contains just the\r
-librares to be loaded by the agent, and the event species to monitor),\r
-\texttt{QAdescriptions\-EMC.configfile} (it contains the descriptions\r
-of the histograms shipped to the dqm shifter), and\r
-\texttt{QAThresholds.configfile} (it contains the thresholds for the\r
-QA checker - see above). You can edit them by pressing the edit button.\r
-\r
-\paragraph{How to run AMORE in the test machine:}\r
-\r
-The EMCal AMORE test machine can be reached via \texttt{ssh\r
- emcal@pcaldbl601}. The standard setup there is identical to the\r
-setup at P2. In the home directory there are some symbolic links which\r
-need to be changed in order to change the setup and test new\r
-features. \texttt{alirootLink} usually links to\r
-\texttt{/opt/aliroot-<some-ver>}, which is the version currently at P2\r
-(daq team updates the software in the \texttt{/opt} directory when\r
-needed). The same holds for \texttt{rootLink}, \texttt{amoreLink} and\r
-\texttt{amoreSiteLink} (you should not need to change\r
-\texttt{rootLink} and \texttt{amoreLink}). For testing some changes to\r
-the EMCAL QA classes, one has to compile an own version of aliroot\r
-with the new version of those classes, and then make a symbolic link\r
-to the path of this aliroot version. There is an aliroot trunk version\r
-which was used to be updated from time to time for tests in \texttt{fblanco/alisoft}. In the directory \texttt{myamoreStuff} there are 2 version of the AMORE modules: \texttt{current\_deploy} and \texttt{trunk} (of course you should do svn up from time to time). Let's suppose you want to make some modification to our AMORE GUI (the expert one) and test them. Remove the \texttt{amoreSiteLink} (it usually points to \texttt{/amoreSite}). Then do:\\\r
-\texttt{ln -s myamoreStuff/amoreSite amoreSiteLink}\\\r
-(or any other directory you would like to use). Then:\\\r
-\texttt{cd trunk[current\_deploy]/amoreEMC}\\\r
-At this point you can modify either the \texttt{src/ui/EMCUIQA} or \texttt{src/ui/EMCUIQATRU} class. This class contains just some manipulations of canvas/histograms to be shown in the expert panel and should be easy to understand. If there are some doubts, just ask. The first one also contains the hack we use in order to use our own reference file at P2 in used to calculate the ratio to reference (next section will explain how to create a new reference file). \texttt{make install} will install the modified libraries into \texttt{amoreSiteLink}. You can use some of the scripts in the \texttt{fblanco} directory to run the agent and the GUI.\r
-In order to run the expert agent with the expert GUI you should do:\\\r
-\texttt{amoreAgent -a EMCQA -s <some-raw-data>\footnote{If you want to create a raw data file, you should first download a chunk of raw data from alien. Then do:\\ \texttt{deroot file.root file.raw}.\\ Keep in mind that the test machine is shared with other detectors, so avoid to store a lot of raw data there and clean up the space from time to time.} -g emcal\_amore.cfg}\\\r
-After doing ssh to the test machine in another terminal, you do \\\r
-\texttt{amore -d EMC -m EMCUIQA(TRU) -g configGUI.txt}\\\r
-If you want to test the shifter agent, you simply do:\\\r
-\texttt{amoreAgent -a EMCQAshifter -s <some-raw-data> -g emcal\_amore.cfg}\\\r
-and type \texttt{amoreGui} in another terminal window.\\\r
-You can commit changes to the trunk (and not to the current\_deploy) with the usual \texttt{svn commit}. Once you feel that your changes are ready to be deployed at P2, send an email with a request to date-support.\r
-\r
-\r
-\r
-\paragraph{How to create a new reference file:}\r
-\r
-In order to create a new reference file, you have to use the \texttt{doReco.sh} script. The only thing you should do there is the path to the raw data in alien and the number of chunks to analyze (check it in the alien path).\r
-Remember to do \texttt{alien-token-init} and \texttt{source /tmp/gclient\_env\_\$UID} before running the script.\r
-After the script is executed you will have some directories (chunck10, chunk11...). Each of them contains an \texttt{EMCAL.QA.<RunNumber>.root} file. You can do :\\\r
-\texttt{hadd EMCAL.QA.0.root chunck10/EMCAL.QA.<RunNumber>.root \newline chunck11/EMCAL.QA.<RunNumber>.root}\\\r
-to merge them (the output files should always be called \texttt{EMCAL.QA.0.root}).\r
-At this point the output file can be already copied to the aldaqacr51 machine in the emcal station, in which we can change the QARef file whenever we want.\r
-In order to do that, you should copy it to your lxplus area, then from the \texttt{\textasciitilde/amore/QARef} directory in the aldaqacr51 you can do:\\\r
-\texttt{scp your-afs-account@aldaqgw01:/afs/cern.ch/<path-to-file>/<file> .}\\\r
-Do:\\\r
-\texttt{ln -s new-file QA.Ref.root}\\\r
-and add a note to the notes file in the directory.\r
-Then you have to create an OCDB file. In order to do that, you must do\\ \r
-\texttt{aliroot Save2OCDB.C}\\\r
-Standard parameters of the macro are ``EMCAL'' (detector name), 0 (run number) and ``12'' (year). You may need to change only the last one. The macro will create a directory named \texttt{QARef/EMCAL/QA/Calib/}, and the file \texttt{Run0\_999999999\_v0\_s0.root} in it. This file has to be committed to the \texttt{QARef/EMCAL/QA/Calib/} directory of aliroot. You can check it with the \texttt{checkCDB.C} macro (it just displays a couple of histograms).\r
-%\subsection{}\r
-\r
-\subsubsection{Some more informations}\r
-Further details about AMORE can be found here:\\\r
-\href{https://ph-dep-aid.web.cern.ch/ph-dep-aid/amore/}{https://ph-dep-aid.web.cern.ch/ph-dep-aid/amore/}\\\r
-The Twiki page with the general EMCal informations for the DQM shifter is:\\\r
-\href{https://twiki.cern.ch/twiki/bin/viewauth/ALICE/EVEEMC}{https://twiki.cern.ch/twiki/bin/viewauth/ALICE/EVEEMC}\r
-The Twiki page (DQM blackboard) with temporary informations for the DQM shifter is:\\\r
-\href{https://twiki.cern.ch/twiki/bin/viewauth/ALICE/DQMBlackboard}{https://twiki.cern.ch/twiki/bin/viewauth/ALICE/DQMBlackboard}\r
-\r
-\subsection{Offline - Marie}\r
-Analysis code, what we control, how\r
-\r
-\subsection{Event display}\r
-\r
+\section{Run by run QA, how to and code}
+\subsection{Online - Francesco, Michael}
+
+\subsubsection{Creation and checking of online QA histograms (AMORE)}
+The QA histograms for the EMCal are created and filled by the class:\\
+\texttt{\$ALICE\_ROOT/EMCAL/AliEMCALQADataMakerRec}\\
+It is run both during the offline reconstruction and by the online data quality monitoring framework. Its main methods are:\\
+
+\texttt{InitRaws()}:
+All the QA histograms are created here. Their titles should be self-explaining. Each of them has 2 important flags:\\
+\texttt{expert} - if true, the histogram is shipped only to the AMORE expert agent (run in the EMCal station), otherwise it is also shipped to the AMORE shifter agent (at the moment we have 4 histograms monitored by the DQM shifter).\\
+\texttt{image} - if true, the plots is saved to the loogbok (there should be 9 plots in the logbook).\\
+Each of those histograms is replicated 4 times by the amore framework and filled according to the event species (Calib, Cosmic, LowMultiplicity, HighMultiplicity).\\
+
+\texttt{MakeRaws(AliRawReader* rawReader)}:
+Here we loop over the input raw data stream and we fill the histograms.\\
+
+\texttt{MakeRawsSTU(AliRawReader* rawReader)}:
+STU raw data decoding, provided EMCal trigger experts.\\
+
+\texttt{GetCalibRefFromOCDB()}:
+Get the reference histograms from the OCDB. The ratio between histograms from current run and the reference run is calculated in the \texttt{MakeRaws} method.\\
+
+
+The QA histograms are then analyzed by the class:\\
+\texttt{\$ALICE\_ROOT/EMCAL/AliEMCALQAChecker}\\
+It checks the data contained in the 4 non-expert histograms. They are at the moment:
+\begin{enumerate}
+\item Ratio distribution of towers' amplitude in the current run w.r.t. the reference run.
+\item Number of hits of L1 Gamma patch.
+\item Number of hits of L1 Jet patch.
+\item Number of links between TRU and STU.
+\end{enumerate}
+The main method is:
+\texttt{CheckRaws()}: For the first histogram, we check how many
+channels have the ratio in the region $0.8-1.2$. If there are more than
+the threshold (default = 90\%), everything is ok, otherwise a red box with a ``call
+expert'' message is displayed. For L1 Jet/Gamma patch, we check if the
+rate of a single patch is higher than the average rate of all other
+patches times a certain threshold value (default = 0.5):
+\begin{equation}
+Rate_{patch}/Rate_{Total} > Threshold /(1 + Threshold)
+\end{equation}
+
+ If so red box with a ``hot spot - call expert'' message is
+ displayed. Finally, if there is a number (default = 1) of missing
+ STU-TRU links, another red box with a warning message is shown.
+If one wants to change the thresholds, this can be done by updating
+them in the database (type \texttt{amoreConfigFileBrowser d emc} at
+any machine at P2 and edit the file \texttt{QAThresholds.configfile} accordingly).
+
+\subsubsection{How to's for EMCal AMORE experts}
+
+\paragraph{How to run AMORE at point 2 (P2):}
+
+Some useful scripts for running AMORE at P2 are located in the \texttt{\textasciitilde/amore} directory in the aldaqacr51 machine. \texttt{StartAmore.sh} is used to run the agent and the GUI (in the expert mode). It simply launches the \texttt{configMonitor.sh} and the \texttt{checkLogs.sh} (used to check the logs of the agents). \texttt{configMonitor.sh} setup the agent and the GUI using some kdialog commands and runs them.
+Basically the agent is run by the following line in the script:\\
+\texttt{sshdqm \$agentName nohup \$scriptsPATH/startAmoreAgent.sh \$onlineMode \&}\\
+The agent is actually run on the EMCal DQM node via the sshdqm command, so this line basically is used to connect to our dqm node and run the agent there.\\
+The GUI is launched in our machine by the command:\\
+\texttt{amore -c EMC -m EMCUIQA(TRU) -g configFile.txt}\\
+The \texttt{configFile.txt} file contains just one line, i.e. \texttt{runType i}, where \texttt {i} runs from 1 to 4 (1=LowMultiplicity, 2=HighMultiplicity, 3=Cosmic, 4=Calib) and it is used to select the event specie you want to display in the GUI (this file is automatically created each time you run \texttt{configMonitor.sh}).\\
+Some other useful scripts (which are launched by \texttt{configMonitor.sh} using the sshdqm) are located in our DQM node. One can login to it by using:\\
+\texttt{ssdqm EMCQA}\\
+They are located in the \texttt{scripts} directory.\\
+\texttt{startAmoreAgent.sh} basically runs the agent, it simply issues the command:\\
+\texttt{amoreAgent -a EMCQA -s @\$1: -u}\\
+The only parameter here is \$1, which is the gdc to monitor, defined in \texttt{configMonitor.sh}.
+\texttt{killAgent.sh} kills any running agent. Be aware that there cannot be 2 running agents of the same kind at the same time (i.e. just one EMCQA and one EMCQAshifter at a time). It can happen that some time the agent does not start when you try to run it. Most of the time this is because the dqm shifter is already running the EMCQA agent instead of just the EMCQAshifter one (they can run it from the dqm station even if they shouldn't). To check if the dqm shifter is running it, simply do \texttt{ps aux | grep EMCQA} in the dmq node, and see if there is an EMCQA process belonging to the user daq. If so, kindly ask the dqm shifter to please kill it.\\
+The agent is always running using the same parameters, which are
+stored in a database. In order to change them, you should run the
+\texttt{amoreConfigFileBrowser} command in the aldaqacr51 machine. A
+window will appear, there you can browse a lot of configuration
+files. Our files are \texttt{EMCAL\_config.txt} (it contains just the
+librares to be loaded by the agent, and the event species to monitor),
+\texttt{QAdescriptions\-EMC.configfile} (it contains the descriptions
+of the histograms shipped to the dqm shifter), and
+\texttt{QAThresholds.configfile} (it contains the thresholds for the
+QA checker - see above). You can edit them by pressing the edit button.
+
+\paragraph{How to run AMORE in the test machine:}
+
+The EMCal AMORE test machine can be reached via \texttt{ssh
+ emcal@pcaldbl601}. The standard setup there is identical to the
+setup at P2. In the home directory there are some symbolic links which
+need to be changed in order to change the setup and test new
+features. \texttt{alirootLink} usually links to
+\texttt{/opt/aliroot-<some-ver>}, which is the version currently at P2
+(daq team updates the software in the \texttt{/opt} directory when
+needed). The same holds for \texttt{rootLink}, \texttt{amoreLink} and
+\texttt{amoreSiteLink} (you should not need to change
+\texttt{rootLink} and \texttt{amoreLink}). For testing some changes to
+the EMCAL QA classes, one has to compile an own version of aliroot
+with the new version of those classes, and then make a symbolic link
+to the path of this aliroot version. There is an aliroot trunk version
+which was used to be updated from time to time for tests in \texttt{fblanco/alisoft}. In the directory \texttt{myamoreStuff} there are 2 version of the AMORE modules: \texttt{current\_deploy} and \texttt{trunk} (of course you should do svn up from time to time). Let's suppose you want to make some modification to our AMORE GUI (the expert one) and test them. Remove the \texttt{amoreSiteLink} (it usually points to \texttt{/amoreSite}). Then do:\\
+\texttt{ln -s myamoreStuff/amoreSite amoreSiteLink}\\
+(or any other directory you would like to use). Then:\\
+\texttt{cd trunk[current\_deploy]/amoreEMC}\\
+At this point you can modify either the \texttt{src/ui/EMCUIQA} or \texttt{src/ui/EMCUIQATRU} class. This class contains just some manipulations of canvas/histograms to be shown in the expert panel and should be easy to understand. If there are some doubts, just ask. The first one also contains the hack we use in order to use our own reference file at P2 in used to calculate the ratio to reference (next section will explain how to create a new reference file). \texttt{make install} will install the modified libraries into \texttt{amoreSiteLink}. You can use some of the scripts in the \texttt{fblanco} directory to run the agent and the GUI.
+In order to run the expert agent with the expert GUI you should do:\\
+\texttt{amoreAgent -a EMCQA -s <some-raw-data>\footnote{If you want to create a raw data file, you should first download a chunk of raw data from alien. Then do:\\ \texttt{deroot file.root file.raw}.\\ Keep in mind that the test machine is shared with other detectors, so avoid to store a lot of raw data there and clean up the space from time to time.} -g emcal\_amore.cfg}\\
+After doing ssh to the test machine in another terminal, you do \\
+\texttt{amore -d EMC -m EMCUIQA(TRU) -g configGUI.txt}\\
+If you want to test the shifter agent, you simply do:\\
+\texttt{amoreAgent -a EMCQAshifter -s <some-raw-data> -g emcal\_amore.cfg}\\
+and type \texttt{amoreGui} in another terminal window.\\
+You can commit changes to the trunk (and not to the current\_deploy) with the usual \texttt{svn commit}. Once you feel that your changes are ready to be deployed at P2, send an email with a request to date-support.
+
+
+
+\paragraph{How to create a new reference file:}
+
+In order to create a new reference file, you have to use the \texttt{doReco.sh} script. The only thing you should do there is the path to the raw data in alien and the number of chunks to analyze (check it in the alien path).
+Remember to do \texttt{alien-token-init} and \texttt{source /tmp/gclient\_env\_\$UID} before running the script.
+After the script is executed you will have some directories (chunck10, chunk11...). Each of them contains an \texttt{EMCAL.QA.<RunNumber>.root} file. You can do :\\
+\texttt{hadd EMCAL.QA.0.root chunck10/EMCAL.QA.<RunNumber>.root \newline chunck11/EMCAL.QA.<RunNumber>.root}\\
+to merge them (the output files should always be called \texttt{EMCAL.QA.0.root}).
+At this point the output file can be already copied to the aldaqacr51 machine in the emcal station, in which we can change the QARef file whenever we want.
+In order to do that, you should copy it to your lxplus area, then from the \texttt{\textasciitilde/amore/QARef} directory in the aldaqacr51 you can do:\\
+\texttt{scp your-afs-account@aldaqgw01:/afs/cern.ch/<path-to-file>/<file> .}\\
+Do:\\
+\texttt{ln -s new-file QA.Ref.root}\\
+and add a note to the notes file in the directory.
+Then you have to create an OCDB file. In order to do that, you must do\\
+\texttt{aliroot Save2OCDB.C}\\
+Standard parameters of the macro are ``EMCAL'' (detector name), 0 (run number) and ``12'' (year). You may need to change only the last one. The macro will create a directory named \texttt{QARef/EMCAL/QA/Calib/}, and the file \texttt{Run0\_999999999\_v0\_s0.root} in it. This file has to be committed to the \texttt{QARef/EMCAL/QA/Calib/} directory of aliroot. You can check it with the \texttt{checkCDB.C} macro (it just displays a couple of histograms).
+%\subsection{}
+
+\subsubsection{Some more informations}
+Further details about AMORE can be found here:\\
+\href{https://ph-dep-aid.web.cern.ch/ph-dep-aid/amore/}{https://ph-dep-aid.web.cern.ch/ph-dep-aid/amore/}\\
+The Twiki page with the general EMCal informations for the DQM shifter is:\\
+\href{https://twiki.cern.ch/twiki/bin/viewauth/ALICE/EVEEMC}{https://twiki.cern.ch/twiki/bin/viewauth/ALICE/EVEEMC}
+The Twiki page (DQM blackboard) with temporary informations for the DQM shifter is:\\
+\href{https://twiki.cern.ch/twiki/bin/viewauth/ALICE/DQMBlackboard}{https://twiki.cern.ch/twiki/bin/viewauth/ALICE/DQMBlackboard}
+
+\subsection{Offline - Marie}
+\label{sec:QAOffline}
+
+The aim of the offline QA is to check EMCAL data quality but also to get information of the global run quality with respect to EMCAL. The offline QA is automaticclly processed for each official data reconstruction pass, it can be calibration passes, validation passes, production passes or specific early dedicated production such as the ``muoncalo'' production for calorimeters and muon detectors. The user can then analyse the outputs created during official reconstruction passes which are stored on alien. The output are in that case stored in /alice/data/year/period/pass/QAresults.root files.
+But each analyser can also add the QA task in its ianalysis in order to perform the QA of his dataset.
+
+
+The following sections describes how to run the QA task then how to analyse the ouputs run by run and how to look at EMCAL stability via some trending observables for whole periods/year.
+
+\subsubsection{Creation of offline QA histograms }
+
+The offline QA histograms for the EMCal are created and filled by the class:\\
+\texttt{\$ALICE\_ROOT/PWGGA/CaloTrackCorrelations/AliAnaCalorimeterQA}\\
+This class takes the information on AliESD/AliAOD objects and fills histograms related to EMCAL cluster or cells but also on correlation between EMCAL and the other detectors in the run.
+Its main methods are:
+\begin{itemize}
+
+\item \texttt{ MakeAnalysisFillHistograms()}: main loop where the histogram filling once the clusters and cells are loaded.
+\item \texttt{ ClusterLoopHistograms(const TObjArray $\*$ clusters, AliVCaloCells $\*$ cells)}: fill of cluster related histograms.
+\item \texttt{ CellHistograms(AliVCaloCells $\*$ cells)}: fill cells related histograms.
+\item \texttt{ Correlate()}: this methods allows to correlate clusters observables to other detector observables: tracks, VO signal, PHOS clusters.
+\item \texttt{ InvariantMassHistograms(const Int$\_$t iclus, const TLorentzVector mom, const Int$\_$t nModule,const TObjArray$\*$ caloClusters, AliVCaloCells $\*$ cells)}: fill Invariant mass histograms in cluster loop.
+\end{itemize}
+
+The way to initialize the macro is the following: \\
+\texttt{AddTaskCalorimeterQA(Bool$\_t$ kSimulation = kFALSE,const char *suffix = default, TString outputFile='' '' ,Int$\_t$ year = 2012, Bool$\_t$ kPrintSettings = kFALSE)}
+
+where
+\begin{itemize}
+
+\item kSimulation: TRUE if you want to run all the MC histograms.
+\item suffix: suffix added to $CaloQA\_$ to the name of the TDirectory and TList where the output histograms will be stored. This is important in case you want to add the task for different \textbf{AliVEvent} offline trigger type set by SelectCollisionCandidates(kTriggerMask).
+\item outputFile: this will be let to the Analysis manager to define it;
+\item year: by default to 2012 will set all the SuperModules histograms (12) even if data from previous years are analysed. (if sets to 2010 or 2011 the associated nb of SM (10 for 2011 to 2013 ; 4 for 2010)
+\item kPrintSettings: to have lots of printings.
+
+\end{itemize}
+
+An example of configuration macro for this class can be found in \\
+\texttt{\$ALICE\_ROOT/PWGGA/CaloTrackCorrelations/AliAnaCalorimeterQA/macros/\\QA/AddTaskCalorimeterQA.C}\\
+
+Lots of settings for this class can be configured via the macro
+such as the filling of MonteCarlo info, the range of the histograms, the different histogram class filling.
+
+The ingredient to be passed to the \texttt{AddTaskCalorimeterQA.C} macro are the ones of\\
+\texttt{\$ALICE\_ROOT/PWG/CaloTrackCorrBase/AliAnalysisTaskCaloTrackCorrelation.cxx}.
+
+An example of how to call it can be found in
+\texttt{\$ALICE\_ROOT/PWGGA/CaloTrackCorrelations/\\AliAnaCalorimeterQA/macros/ana.C}\\
+
+\subsubsection{Run by Run Quality Assurance}
+
+ The Run by Run Quality assurance is made by analyzing different histograms.
+A macro \\ \texttt{\$ALICE\_ROOT/PWGGA/CaloTrackCorrelations/AliAnaCalorimeterQA/macros/QA/QAplots.C} is used systematically on all the runs.
+This macro runs on the .root file generated by the\\ \texttt{\$ALICE\_ROOT/PWGGA/CaloTrackCorrelations/AliAnaCalorimeterQA/macros/ \\ QA/AddTaskCalorimeterQA.C} and looks at several histograms generated by the \texttt{AliAnaCalorimeterQA} task.
+To run this macro it is supposed that the output (.root files) of the QA task are in a directory named {\bf period/pass} where period is the period (e.g. LHC12h) and pass the reconstruction pass you are looking at. The .root files are supposed to be renamed {\bf runnumber.root} in case of production passes, validation passes or dedicated calo production passes and {\bf runnumber$\_$outer.root} or {\bf runnumber$\_$barrel.root} in case of cpass (calibration pass). In the calibration pass case, the Minimum Bias triggers are reconstructed within the $\_$barrel.root extansion and the specific triggers (muon, EMCAL, ...) are reconstructed with the $\_$outer.root extension. The QA task is called twice at the end of reconstruction to separate both kind of triggers.
+
+To run the macro you must have in the period/pass repository a text file called ``runlist.txt containing the list of run you whant to plot, the runlist.txt beeing of the type:\\
+1 ~~~~~run1\\
+2 ~~~~~run2\\
+.....\\
+The output plots are saved automatically in repository period/pass/runnumber so that you should create those repository prior to running the macro.\\
+
+
+
+The arguments to pass to the macro:
+\texttt{ QAplots(TString fCalorimeter = "EMCAL", TString period = "LHC11h", TString pass = "pass", TString fTrigger= "default", TString system = "PbPb")} are the following:
+\begin{itemize}
+\item TString fCalorimeter : name of the calorimeter: "EMCAL".
+\item TString period: path(name) of the repository for the period.
+\item TString pass: name of the repository where you saved the .root files (in the example period/pass).
+\item TString fTtrigger type: "default" for the minbias; "trigEMC" or "EMC7" for EMCAL triggers. This corresponds to the TString suffix of the AddTaskCalorimeterQA() when initializing the task.
+\item TString system: "PbPb" or "pp": a flag for setting/adjusting scales for some drawings to the multiplicities/energies reached in the system.
+\end{itemize}
+
+
+For typical automatic QA of official reconstruction passes the histograms we check automatically with runnning of the QAplots.C macro are the following: but the user can add as much histograms verifications as he wants.
+\begin{itemize}
+\item Occupancy map: $\eta$ $\phi$ distributions of cells with signsl (in row/col units).
+\item Summed energy map: $\eta$ $\phi$ distribution of mean energy/event (in row/col units).
+\item Time vs E of custers: Time vs Energy of clusters.
+\item EMCAL Cluster versus track correlation (multiplicities/energies).
+\item EMCAL Clusters vs VO signal correlation.
+\item EMCAL vs PHOS clusters correlations.
+\end{itemize}
+
+On top of those default settings and checks, the user can add as much histograms checks as he wants.
+
+\subsubsection{Period Quality Assurance and stability of EMCAL}
+
+The stability of EMCAL is checked also systematically at the end of each reconstruction passes (calibration pass, validation pass, production pass).
+This step is done via a set of macros called ``trending'' macros which can be found in \texttt{\$ALICE\_ROOT/PWGGA/CaloTrackCorrelations/\\ AliAnaCalorimeterQA/macros/QA/}. The goal of these two macros (trendingClusters.C and trendingPi0.C) is to look at averages observables concerning clusters and $\pi^0$ respectively, as a function of runnumber.
+Those macros uses the output of .root file generated by the\\ \texttt{\$ALICE\_ROOT/PWGGA/CaloTrackCorrelations/AliAnaCalorimeterQA/macros/ \\ QA/AddTaskCalorimeterQA.C}.
+It requires the same repository, naming conventions as in the previous section: period/pass/runnumber.root files and a runlist.txt file containing the list of runs to check. It can be run depending on the ``fTrigger" option EMCAL triggered events and on the MinBias events.
+
+{\bf Cluster Averages: (trendingCluster.C):}
+
+This macro computes some averages over histograms and store them in a TGraph as a fuction of runnumber.
+The folling histograms are used:
+\begin{itemize}
+\item EMCAL\_hE: cluster energy distribution in EMCAL.
+\item EMCAL\_hNClusters\_Mod\_xx: 2D histogram with number of cluster per supermodule.
+\item EMCAL\_hNCellsPerCluster\_Mod\_xx: 2D histogram with number of cells per supermodule.
+\item EMCAL\_hE\_Mod\_xx: 2D histogram with Cluster Energy per supermodule.
+\end{itemize}
+
+
+and the following averages are computed and displayed as a funcion of runnumber:
+
+\begin{itemize}
+\item Mean number of clusters per event is computed over EMCAL and also per supermodule. An example is displayed on figure \ref{fig:trendingCluster} for LHC12i period.
+\item Average energy of clusters (with given threshold energies) per event is computed over EMCAL and also per supermodule.
+\item Number of cells of clusters (with given threshold energies)per event is computed over EMCAL and also per supermodule.
+\end{itemize}
+
+
+\begin{figure}[ht]
+\begin{center}
+\includegraphics[width=0.9\textwidth]{figures/trendingClusterLHC12iMB.pdf}
+\end{center}
+\caption{\label{fig:trendingCluster}
+Example of trending plot for period LHC12i: mean number of cluster per event as a function of runnumber for MinBias events.}
+\end{figure}
+
+
+
+{\bf $\pi^0$ Averages: (trendingPi0.C):}
+
+This macro computes some fits over invariant mass histograms and store the fit parameters and $\pi^0$ computed values in a TGraph as a funcion of runnumber. The following histograms are used:
+\begin{itemize}
+\item EMCAL\_hIM: Invariant mass histogram filled over the whole cluster pairs.
+\item EMCAL$\_$hIM$\_$Mod$\_$xx: Invariant mass histogram filled over the cluster pairs in the same supermodule.
+\end{itemize}
+
+The fit is done by a gaussian plus a 2$^{nd}$ order polinom for background. This is a very ``raw'' fitting wich may not be appropriate especially for the PbPb cases but allows to check the stabilities of averages (mass, width, number of $\pi^O$).
+
+
+
+\begin{itemize}
+\item Mean number of $\pi^0$ per event is computed over EMCAL and also per supermodule.
+\item Mass of the $\pi^0$ per event is computed over EMCAL and also per supermodule.
+\item Witdh of the $\pi^0$ per event is computed over EMCAL and also per supermodule.
+\end{itemize}
+
+
+
+The use of these set of QA macros allows to identify problematic runs (subperiod with missing part of detectors) or more general features. Also it can be used to check global observables for different dataset used for analysis purpose.
+
+
+
+\subsection{Event display}
+
\subsection{Logbook tips}
\ No newline at end of file
\begin{figure}[ht]
\begin{center}
-\includegraphics[width=0.8\textwidth]{figures/EMCalGeometryStructure.eps}
+\includegraphics[width=0.8\textwidth]{figures/EMCalGeometryStructure.pdf}
\end{center}
\caption{\label{fig:GeometryDirectory}
Here is shown a typical hierarchical geometry structure. This
wide defaults are used.
\begin{tabular*}{\textwidth}[ht]{p{.5\textwidth}p{0.5\textwidth}}
-\includegraphics[width=0.5\textwidth]{figures/EMCalMCStep.eps}
+\includegraphics[width=0.5\textwidth]{figures/EMCalMCStep.pdf}
\label{fig:emcalStepManager_ParticleStep}
&
\vspace{-10.5cm}
--- /dev/null
+%!PS-Adobe-2.0 EPSF-2.0
+%%BoundingBox: 0 0 567 269
+%%Title: /scratch/alicehp2/germain/oldQA/LHC12i/vpass1/LHC12ivpass1MBClAvEnt.eps: ClusterAverages
+%%Creator: ROOT Version 5.34/02
+%%CreationDate: Thu Apr 11 15:59:10 2013
+%%EndComments
+%%BeginProlog
+80 dict begin
+/s {stroke} def /l {lineto} def /m {moveto} def /t {translate} def
+/r {rotate} def /rl {roll} def /R {repeat} def
+/d {rlineto} def /rm {rmoveto} def /gr {grestore} def /f {eofill} def
+/c {setrgbcolor} def /black {0 setgray} def /sd {setdash} def
+/cl {closepath} def /sf {scalefont setfont} def /lw {setlinewidth} def
+/box {m dup 0 exch d exch 0 d 0 exch neg d cl} def
+/NC{systemdict begin initclip end}def/C{NC box clip newpath}def
+/bl {box s} def /bf {box f} def /Y { 0 exch d} def /X { 0 d} def
+/K {{pop pop 0 moveto} exch kshow} bind def
+/ita {/ang 15 def gsave [1 0 ang dup sin exch cos div 1 0 0] concat} def
+/mp {newpath /y exch def /x exch def} def
+/side {[w .77 mul w .23 mul] .385 w mul sd w 0 l currentpoint t -144 r} def
+/mr {mp x y w2 0 360 arc} def /m24 {mr s} def /m20 {mr f} def
+/mb {mp x y w2 add m w2 neg 0 d 0 w neg d w 0 d 0 w d cl} def
+/mt {mp x y w2 add m w2 neg w neg d w 0 d cl} def
+/m21 {mb f} def /m25 {mb s} def /m22 {mt f} def /m26{mt s} def
+/m23 {mp x y w2 sub m w2 w d w neg 0 d cl f} def
+/m27 {mp x y w2 add m w3 neg w2 neg d w3 w2 neg d w3 w2 d cl s} def
+/m28 {mp x w2 sub y w2 sub w3 add m w3 0 d 0 w3 neg d w3 0 d 0 w3 d w3 0 d 0 w3 d w3 neg 0 d 0 w3 d w3 neg 0 d 0 w3 neg d w3 neg 0 d cl s } def
+/m29 {mp gsave x w2 sub y w2 add w3 sub m currentpoint t 4 {side} repeat cl fill gr} def
+/m30 {mp gsave x w2 sub y w2 add w3 sub m currentpoint t 4 {side} repeat cl s gr} def
+/m31 {mp x y w2 sub m 0 w d x w2 sub y m w 0 d x w2 sub y w2 add m w w neg d x w2 sub y w2 sub m w w d s} def
+/m32 {mp x y w2 sub m w2 w d w neg 0 d cl s} def
+/m33 {mp x y w2 add m w3 neg w2 neg d w3 w2 neg d w3 w2 d cl f} def
+/m34 {mp x w2 sub y w2 sub w3 add m w3 0 d 0 w3 neg d w3 0 d 0 w3 d w3 0 d 0 w3 d w3 neg 0 d 0 w3 d w3 neg 0 d 0 w3 neg d w3 neg 0 d cl f } def
+/m2 {mp x y w2 sub m 0 w d x w2 sub y m w 0 d s} def
+/m5 {mp x w2 sub y w2 sub m w w d x w2 sub y w2 add m w w neg d s} def
+/reEncode {exch findfont dup length dict begin {1 index /FID eq {pop pop} {def} ifelse } forall /Encoding exch def currentdict end dup /FontName get exch definefont pop } def [/Times-Bold /Times-Italic /Times-BoldItalic /Helvetica /Helvetica-Oblique
+ /Helvetica-Bold /Helvetica-BoldOblique /Courier /Courier-Oblique /Courier-Bold /Courier-BoldOblique /Times-Roman /AvantGarde-Book /AvantGarde-BookOblique /AvantGarde-Demi /AvantGarde-DemiOblique /Bookman-Demi /Bookman-DemiItalic /Bookman-Light
+ /Bookman-LightItalic /Helvetica-Narrow /Helvetica-Narrow-Bold /Helvetica-Narrow-BoldOblique /Helvetica-Narrow-Oblique /NewCenturySchlbk-Roman /NewCenturySchlbk-Bold /NewCenturySchlbk-BoldItalic /NewCenturySchlbk-Italic /Palatino-Bold
+ /Palatino-BoldItalic /Palatino-Italic /Palatino-Roman ] {ISOLatin1Encoding reEncode } forall
+%%EndProlog
+%%BeginSetup
+%%EndSetup
+newpath gsave .25 .25 scale gsave 0 0 t black[ ] 0 sd 3 lw 1 1 1 c 2268 1076 0 0 bf black 1 1 1 c 2041 861 181 108 bf black 2041 861 181 108 bl 1 1 1 c 2041 861 181 108 bf black 2041 861 181 108 bl 181 108 m 2041 X s[ 4 8] 0 sd 181 969 m -861 Y
+ s 277 969 m -861 Y s 374 969 m -861 Y s 470 969 m -861 Y s 566 969 m -861 Y s 662 969 m -861 Y s 758 969 m -861 Y s 854 969 m -861 Y s 950 969 m -861 Y s 1046 969 m -861 Y s 1142 969 m -861 Y s 1238 969 m -861 Y s 1334 969 m -861 Y s 1430 969 m -861
+ Y s 1526 969 m -861 Y s 1622 969 m -861 Y s 1718 969 m -861 Y s 1814 969 m -861 Y s 1910 969 m -861 Y s 2006 969 m -861 Y s 2102 969 m -861 Y s 2198 969 m -861 Y s 2198 969 m -861 Y s[ ] 0 sd 181 108 m 861 Y s[ 4 8] 0 sd 2222 165 m -2041 X s 2222
+ 280 m -2041 X s 2222 395 m -2041 X s 2222 510 m -2041 X s 2222 624 m -2041 X s 2222 739 m -2041 X s 2222 854 m -2041 X s 2222 969 m -2041 X s 2222 165 m -2041 X s 0 0 0.6 c[ ] 0 sd 1 1 1 c black 0 0 0.6 c 181 108 m 2041 X s black 181 108 m 2041 X s
+ gsave 2268 1076 0 0 C 212.26 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 192772) show NC gr
+ gsave 2268 1076 0 0 C 262.076 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 192775) show NC gr
+ gsave 2268 1076 0 0 C 309.726 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 192778) show NC gr
+ gsave 2268 1076 0 0 C 359.542 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 192779) show NC gr
+ gsave 2268 1076 0 0 C 407.193 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 192820) show NC gr
+ gsave 2268 1076 0 0 C 457.009 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 192822) show NC gr
+ gsave 2268 1076 0 0 C 504.659 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 192824) show NC gr
+ gsave 2268 1076 0 0 C 554.475 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193004) show NC gr
+ gsave 2268 1076 0 0 C 602.125 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193005) show NC gr
+ gsave 2268 1076 0 0 C 649.776 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193007) show NC gr
+ gsave 2268 1076 0 0 C 699.592 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193008) show NC gr
+ gsave 2268 1076 0 0 C 747.242 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193010) show NC gr
+ gsave 2268 1076 0 0 C 797.058 15.1614 t 90 r /Helvetica findfont 23.8251 sf 0 0 m
+ 69 58 45 32 19 6
+( 193011) K NC gr
+ gsave 2268 1076 0 0 C 844.708 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193014) show NC gr
+ gsave 2268 1076 0 0 C 894.524 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193047) show NC gr
+ gsave 2268 1076 0 0 C 942.174 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193049) show NC gr
+ gsave 2268 1076 0 0 C 989.825 12.9955 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193051) show NC gr
+ gsave 2268 1076 0 0 C 1039.64 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193092) show NC gr
+ gsave 2268 1076 0 0 C 1087.29 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193093) show NC gr
+ gsave 2268 1076 0 0 C 1137.11 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193094) show NC gr
+ gsave 2268 1076 0 0 C 1184.76 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193097) show NC gr
+ gsave 2268 1076 0 0 C 1234.57 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193148) show NC gr
+ gsave 2268 1076 0 0 C 1282.22 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193150) show NC gr
+ gsave 2268 1076 0 0 C 1329.87 12.9955 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193151) show NC gr
+ gsave 2268 1076 0 0 C 1379.69 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193152) show NC gr
+ gsave 2268 1076 0 0 C 1427.34 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193155) show NC gr
+ gsave 2268 1076 0 0 C 1477.16 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193156) show NC gr
+ gsave 2268 1076 0 0 C 1524.81 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193184) show NC gr
+ gsave 2268 1076 0 0 C 1574.62 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193187) show NC gr
+ gsave 2268 1076 0 0 C 1622.27 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193188) show NC gr
+ gsave 2268 1076 0 0 C 1672.09 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193189) show NC gr
+ gsave 2268 1076 0 0 C 1719.74 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193192) show NC gr
+ gsave 2268 1076 0 0 C 1767.39 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193194) show NC gr
+ gsave 2268 1076 0 0 C 1817.21 12.9955 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193301) show NC gr
+ gsave 2268 1076 0 0 C 1864.86 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193750) show NC gr
+ gsave 2268 1076 0 0 C 1914.67 12.9955 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193751) show NC gr
+ gsave 2268 1076 0 0 C 1962.32 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193752) show NC gr
+ gsave 2268 1076 0 0 C 2012.14 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193757) show NC gr
+ gsave 2268 1076 0 0 C 2059.79 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193758) show NC gr
+ gsave 2268 1076 0 0 C 2107.44 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193759) show NC gr
+ gsave 2268 1076 0 0 C 2157.25 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193760) show NC gr
+ gsave 2268 1076 0 0 C 2204.9 8.66367 t 90 r /Helvetica findfont 23.8251 sf 0 0 m ( 193766) show NC gr 181 137 m -29 Y s 277 137 m -29 Y s 374 137 m -29 Y s 470 137 m -29 Y s 566 137 m -29 Y s 662 137 m -29 Y s 758 137 m -29 Y s 854 137 m -29 Y s
+ 950 137 m -29 Y s 1046 137 m -29 Y s 1142 137 m -29 Y s 1238 137 m -29 Y s 1334 137 m -29 Y s 1430 137 m -29 Y s 1526 137 m -29 Y s 1622 137 m -29 Y s 1718 137 m -29 Y s 1814 137 m -29 Y s 1910 137 m -29 Y s 2006 137 m -29 Y s 2102 137 m -29 Y s
+ 2198 137 m -29 Y s 2198 137 m -29 Y s 181 108 m 861 Y s
+ gsave 2268 1076 0 0 C 62.8116 855.538 t 90 r /Helvetica findfont 34.6547 sf 0 0 m (>/event) show NC gr
+ gsave 2268 1076 0 0 C 71.4753 773.233 t 90 r /Helvetica findfont 23.8251 sf 0 0 m (clusters) show NC gr
+ gsave 2268 1076 0 0 C 62.8116 729.914 t 90 r /Helvetica findfont 34.6547 sf 0 0 m (<N) show NC gr 236 165 m -55 X s 209 194 m -28 X s 209 222 m -28 X s 209 251 m -28 X s 236 280 m -55 X s 209 309 m -28 X s 209 337 m -28 X s 209 366 m -28 X s 236
+ 395 m -55 X s 209 423 m -28 X s 209 452 m -28 X s 209 481 m -28 X s 236 510 m -55 X s 209 538 m -28 X s 209 567 m -28 X s 209 596 m -28 X s 236 624 m -55 X s 209 653 m -28 X s 209 682 m -28 X s 209 710 m -28 X s 236 739 m -55 X s 209 768 m -28 X s
+ 209 797 m -28 X s 209 825 m -28 X s 236 854 m -55 X s 209 883 m -28 X s 209 911 m -28 X s 209 940 m -28 X s 236 969 m -55 X s 236 165 m -55 X s 209 136 m -28 X s 209 108 m -28 X s
+ gsave 2268 1076 0 0 C 121.291 153.78 t 0 r /Helvetica findfont 34.6547 sf 0 0 m (0.6) show NC gr
+ gsave 2268 1076 0 0 C 121.291 268.574 t 0 r /Helvetica findfont 34.6547 sf 0 0 m (0.8) show NC gr
+ gsave 2268 1076 0 0 C 155.946 381.202 t 0 r /Helvetica findfont 34.6547 sf 0 0 m (1) show NC gr
+ gsave 2268 1076 0 0 C 121.291 495.995 t 0 r /Helvetica findfont 34.6547 sf 0 0 m (1.2) show NC gr
+ gsave 2268 1076 0 0 C 121.291 610.789 t 0 r /Helvetica findfont 34.6547 sf 0 0 m (1.4) show NC gr
+ gsave 2268 1076 0 0 C 121.291 727.749 t 0 r /Helvetica findfont 34.6547 sf 0 0 m (1.6) show NC gr
+ gsave 2268 1076 0 0 C 121.291 842.542 t 0 r /Helvetica findfont 34.6547 sf 0 0 m (1.8) show NC gr
+ gsave 2268 1076 0 0 C 149.448 955.17 t 0 r /Helvetica findfont 34.6547 sf 0 0 m (2) show NC gr 221 457 m -16 X s 205 453 m 8 Y s 238 457 m 15 X s 253 453 m 8 Y s 269 460 m -16 X s 253 456 m 8 Y s 286 460 m 15 X s 301 456 m 8 Y s 317 461 m -16 X s
+ 301 457 m 8 Y s 334 461 m 15 X s 349 457 m 8 Y s 365 456 m -16 X s 349 452 m 9 Y s 382 456 m 16 X s 398 452 m 9 Y s 413 460 m -15 X s 398 456 m 8 Y s 430 460 m 16 X s 446 456 m 8 Y s 461 460 m -15 X s 446 456 m 8 Y s 478 460 m 16 X s 494 456 m 8 Y s
+ 509 457 m -15 X s 494 452 m 9 Y s 526 457 m 16 X s 542 452 m 9 Y s 557 462 m -15 X s 542 458 m 9 Y s 574 462 m 16 X s 590 458 m 9 Y s 605 458 m -15 X s 590 454 m 9 Y s 622 458 m 16 X s 638 454 m 9 Y s 653 461 m -15 X s 638 457 m 8 Y s 670 461 m 16 X
+ s 686 457 m 8 Y s 701 459 m -15 X s 686 455 m 9 Y s 718 459 m 16 X s 734 455 m 9 Y s 749 461 m -15 X s 734 457 m 9 Y s 766 461 m 16 X s 782 457 m 9 Y s 797 458 m -15 X s 782 454 m 9 Y s 814 458 m 16 X s 830 454 m 9 Y s 845 460 m -15 X s 830 456 m 8
+ Y s 862 460 m 16 X s 878 456 m 8 Y s 893 462 m -15 X s 878 457 m 9 Y s 910 462 m 16 X s 926 457 m 9 Y s 941 461 m -15 X s 926 457 m 8 Y s 958 461 m 16 X s 974 457 m 8 Y s 989 459 m -15 X s 974 455 m 9 Y s 1006 459 m 16 X s 1022 455 m 9 Y s 1037 467
+ m -15 X s 1022 463 m 9 Y s 1054 467 m 16 X s 1070 463 m 9 Y s 1085 465 m -15 X s 1070 460 m 9 Y s 1103 465 m 15 X s 1118 460 m 9 Y s 1133 462 m -15 X s 1118 457 m 9 Y s 1151 462 m 15 X s 1166 457 m 9 Y s 1181 466 m -15 X s 1166 461 m 9 Y s 1199 466
+ m 15 X s 1214 461 m 9 Y s 1229 460 m -15 X s 1214 456 m 9 Y s 1247 460 m 15 X s 1262 456 m 9 Y s 1277 465 m -15 X s 1262 460 m 9 Y s 1295 465 m 15 X s 1310 460 m 9 Y s 1325 477 m -15 X s 1310 472 m 9 Y s 1343 477 m 15 X s 1358 472 m 9 Y s 1373 461 m
+ -15 X s 1358 457 m 9 Y s 1391 461 m 15 X s 1406 457 m 9 Y s 1421 464 m -15 X s 1406 459 m 9 Y s 1439 464 m 15 X s 1454 459 m 9 Y s 1469 461 m -15 X s 1454 457 m 8 Y s 1487 461 m 15 X s 1502 457 m 8 Y s 1517 464 m -15 X s 1502 460 m 8 Y s 1535 464 m
+ 15 X s 1550 460 m 8 Y s 1565 468 m -15 X s 1550 463 m 9 Y s 1583 468 m 15 X s 1598 463 m 9 Y s 1613 464 m -15 X s 1598 459 m 9 Y s 1631 464 m 15 X s 1646 459 m 9 Y s 1661 464 m -15 X s 1646 460 m 8 Y s 1679 464 m 15 X s 1694 460 m 8 Y s 1757 461 m
+ -15 X s 1742 457 m 8 Y s 1775 461 m 15 X s 1790 457 m 8 Y s 1854 462 m -16 X s 1838 458 m 8 Y s 1871 462 m 15 X s 1886 458 m 8 Y s 1902 468 m -16 X s 1886 463 m 9 Y s 1919 468 m 15 X s 1934 463 m 9 Y s 1950 475 m -16 X s 1934 471 m 9 Y s 1967 475 m
+ 15 X s 1982 471 m 9 Y s 1998 491 m -16 X s 1982 487 m 9 Y s 2015 491 m 15 X s 2030 487 m 9 Y s 2006 500 m 19 Y s 2002 519 m 9 X s 2006 483 m -19 Y s 2002 464 m 9 X s 2094 475 m -16 X s 2078 471 m 9 Y s 2111 475 m 15 X s 2126 471 m 9 Y s 2142 480 m
+ -16 X s 2126 475 m 9 Y s 2159 480 m 15 X s 2174 475 m 9 Y s 2190 477 m -16 X s 2174 472 m 9 Y s 2207 477 m 15 X s 2222 472 m 9 Y s /w 17 def /w2 {w 2 div} def /w3 {w 3 div} def 229 457 277 460 325 461 374 456 422 460 470 460 518 457 566 462 614 458
+ 662 461 710 459 758 461 806 458 854 460 902 462 950 461 998 459 1046 467 1094 465 1142 462 1190 466 1238 460 1286 465 1334 477 1382 461 1430 464 1478 461 1526 464 1574 468 1622 464 1670 464 31 { m20} R 1766 461 m20 1862 462 1910 468 1958 475 2006
+ 491 4 { m20} R 2102 475 2150 480 2198 477 3 { m20} R 221 470 m -16 X s 205 466 m 8 Y s 238 470 m 15 X s 253 466 m 8 Y s 269 469 m -16 X s 253 464 m 9 Y s 286 469 m 15 X s 301 464 m 9 Y s 317 456 m -16 X s 301 452 m 9 Y s 334 456 m 15 X s 349 452 m 9
+ Y s 365 460 m -16 X s 349 456 m 8 Y s 382 460 m 16 X s 398 456 m 8 Y s 413 464 m -15 X s 398 459 m 9 Y s 430 464 m 16 X s 446 459 m 9 Y s 461 460 m -15 X s 446 456 m 9 Y s 478 460 m 16 X s 494 456 m 9 Y s 470 469 m 2 Y s 465 471 m 9 X s 470 452 m -2
+ Y s 465 450 m 9 X s 509 456 m -15 X s 494 452 m 8 Y s 526 456 m 16 X s 542 452 m 8 Y s 557 464 m -15 X s 542 460 m 8 Y s 574 464 m 16 X s 590 460 m 8 Y s 605 459 m -15 X s 590 455 m 8 Y s 622 459 m 16 X s 638 455 m 8 Y s 653 459 m -15 X s 638 455 m
+ 9 Y s 670 459 m 16 X s 686 455 m 9 Y s 701 466 m -15 X s 686 461 m 9 Y s 718 466 m 16 X s 734 461 m 9 Y s 749 467 m -15 X s 734 463 m 9 Y s 766 467 m 16 X s 782 463 m 9 Y s 797 464 m -15 X s 782 460 m 9 Y s 814 464 m 16 X s 830 460 m 9 Y s 845 466 m
+ -15 X s 830 461 m 9 Y s 862 466 m 16 X s 878 461 m 9 Y s 893 463 m -15 X s 878 459 m 9 Y s 910 463 m 16 X s 926 459 m 9 Y s 941 462 m -15 X s 926 458 m 8 Y s 958 462 m 16 X s 974 458 m 8 Y s 989 466 m -15 X s 974 462 m 8 Y s 1006 466 m 16 X s 1022
+ 462 m 8 Y s 1037 469 m -15 X s 1022 465 m 8 Y s 1054 469 m 16 X s 1070 465 m 8 Y s 1085 467 m -15 X s 1070 463 m 9 Y s 1103 467 m 15 X s 1118 463 m 9 Y s 1133 464 m -15 X s 1118 460 m 9 Y s 1151 464 m 15 X s 1166 460 m 9 Y s 1181 466 m -15 X s 1166
+ 462 m 8 Y s 1199 466 m 15 X s 1214 462 m 8 Y s 1229 462 m -15 X s 1214 458 m 8 Y s 1247 462 m 15 X s 1262 458 m 8 Y s 1277 468 m -15 X s 1262 464 m 9 Y s 1295 468 m 15 X s 1310 464 m 9 Y s 1325 477 m -15 X s 1310 473 m 9 Y s 1343 477 m 15 X s 1358
+ 473 m 9 Y s 1334 486 m 2 Y s 1330 488 m 8 X s 1334 469 m -2 Y s 1330 467 m 8 X s 1373 463 m -15 X s 1358 458 m 9 Y s 1391 463 m 15 X s 1406 458 m 9 Y s 1421 476 m -15 X s 1406 472 m 8 Y s 1439 476 m 15 X s 1454 472 m 8 Y s 1469 465 m -15 X s 1454
+ 461 m 8 Y s 1487 465 m 15 X s 1502 461 m 8 Y s 1517 478 m -15 X s 1502 474 m 8 Y s 1535 478 m 15 X s 1550 474 m 8 Y s 1565 465 m -15 X s 1550 461 m 8 Y s 1583 465 m 15 X s 1598 461 m 8 Y s 1613 469 m -15 X s 1598 464 m 9 Y s 1631 469 m 15 X s 1646
+ 464 m 9 Y s 1661 448 m -15 X s 1646 443 m 9 Y s 1679 448 m 15 X s 1694 443 m 9 Y s 1670 456 m 1 Y s 1666 457 m 8 X s 1670 439 m -1 Y s 1666 438 m 8 X s 1757 468 m -15 X s 1742 464 m 9 Y s 1775 468 m 15 X s 1790 464 m 9 Y s 1854 470 m -16 X s 1838
+ 465 m 9 Y s 1871 470 m 15 X s 1886 465 m 9 Y s 1902 430 m -16 X s 1886 426 m 9 Y s 1919 430 m 15 X s 1934 426 m 9 Y s 1910 439 m 8 Y s 1906 447 m 9 X s 1910 422 m -8 Y s 1906 414 m 9 X s 1950 496 m -16 X s 1934 491 m 9 Y s 1967 496 m 15 X s 1982 491
+ m 9 Y s 1998 392 m -16 X s 1982 388 m 8 Y s 2015 392 m 15 X s 2030 388 m 8 Y s 2094 487 m -16 X s 2078 483 m 8 Y s 2111 487 m 15 X s 2126 483 m 8 Y s 2142 488 m -16 X s 2126 484 m 8 Y s 2159 488 m 15 X s 2174 484 m 8 Y s 2190 486 m -16 X s 2174 481
+ m 9 Y s 2207 486 m 15 X s 2222 481 m 9 Y s 1 0 0 c 229 470 277 469 325 456 374 460 422 464 470 460 518 456 566 464 614 459 662 459 710 466 758 467 806 464 854 466 902 463 950 462 998 466 1046 469 1094 467 1142 464 1190 466 1238 462 1286 468 1334 477
+ 1382 463 1430 476 1478 465 1526 478 1574 465 1622 469 1670 448 31 { m21} R 1766 468 m21 1862 470 1910 430 1958 496 2006 392 4 { m21} R 2102 487 2150 488 2198 486 3 { m21} R black 224 459 m -19 X s 205 455 m 9 Y s 235 459 m 18 X s 253 455 m 9 Y s 272
+ 468 m -19 X s 253 464 m 8 Y s 283 468 m 18 X s 301 464 m 8 Y s 320 458 m -19 X s 301 454 m 8 Y s 331 458 m 18 X s 349 454 m 8 Y s 368 458 m -19 X s 349 454 m 9 Y s 379 458 m 19 X s 398 454 m 9 Y s 416 467 m -18 X s 398 463 m 9 Y s 427 467 m 19 X s
+ 446 463 m 9 Y s 464 466 m -18 X s 446 462 m 8 Y s 475 466 m 19 X s 494 462 m 8 Y s 470 475 m s 465 475 m 9 X s 470 457 m s 465 457 m 9 X s 512 457 m -18 X s 494 453 m 8 Y s 523 457 m 19 X s 542 453 m 8 Y s 560 466 m -18 X s 542 462 m 9 Y s 571 466 m
+ 19 X s 590 462 m 9 Y s 608 461 m -18 X s 590 457 m 8 Y s 619 461 m 19 X s 638 457 m 8 Y s 656 459 m -18 X s 638 455 m 9 Y s 667 459 m 19 X s 686 455 m 9 Y s 704 461 m -18 X s 686 457 m 9 Y s 715 461 m 19 X s 734 457 m 9 Y s 752 468 m -18 X s 734 464
+ m 8 Y s 763 468 m 19 X s 782 464 m 8 Y s 801 455 m -19 X s 782 450 m 9 Y s 811 455 m 19 X s 830 450 m 9 Y s 849 462 m -19 X s 830 458 m 8 Y s 859 462 m 19 X s 878 458 m 8 Y s 897 467 m -19 X s 878 462 m 9 Y s 907 467 m 19 X s 926 462 m 9 Y s 945 460
+ m -19 X s 926 455 m 9 Y s 955 460 m 19 X s 974 455 m 9 Y s 993 464 m -19 X s 974 460 m 8 Y s 1003 464 m 19 X s 1022 460 m 8 Y s 1041 465 m -19 X s 1022 460 m 9 Y s 1051 465 m 19 X s 1070 460 m 9 Y s 1089 469 m -19 X s 1070 465 m 8 Y s 1099 469 m 19
+ X s 1118 465 m 8 Y s 1137 459 m -19 X s 1118 454 m 9 Y s 1147 459 m 19 X s 1166 454 m 9 Y s 1185 473 m -19 X s 1166 468 m 9 Y s 1195 473 m 19 X s 1214 468 m 9 Y s 1233 470 m -19 X s 1214 465 m 9 Y s 1243 470 m 19 X s 1262 465 m 9 Y s 1281 463 m -19
+ X s 1262 459 m 8 Y s 1291 463 m 19 X s 1310 459 m 8 Y s 1329 484 m -19 X s 1310 479 m 9 Y s 1339 484 m 19 X s 1358 479 m 9 Y s 1334 492 m 2 Y s 1330 494 m 8 X s 1334 475 m -2 Y s 1330 473 m 8 X s 1377 464 m -19 X s 1358 459 m 9 Y s 1387 464 m 19 X s
+ 1406 459 m 9 Y s 1425 471 m -19 X s 1406 467 m 8 Y s 1435 471 m 19 X s 1454 467 m 8 Y s 1473 460 m -19 X s 1454 455 m 9 Y s 1483 460 m 19 X s 1502 455 m 9 Y s 1521 468 m -19 X s 1502 464 m 9 Y s 1531 468 m 19 X s 1550 464 m 9 Y s 1569 463 m -19 X s
+ 1550 459 m 9 Y s 1579 463 m 19 X s 1598 459 m 9 Y s 1617 456 m -19 X s 1598 452 m 9 Y s 1627 456 m 19 X s 1646 452 m 9 Y s 1665 469 m -19 X s 1646 464 m 9 Y s 1675 469 m 19 X s 1694 464 m 9 Y s 1670 477 m 4 Y s 1666 481 m 8 X s 1670 460 m -4 Y s
+ 1666 456 m 8 X s 1761 464 m -19 X s 1742 460 m 9 Y s 1771 464 m 19 X s 1790 460 m 9 Y s 1857 456 m -19 X s 1838 451 m 9 Y s 1867 456 m 19 X s 1886 451 m 9 Y s 1905 412 m -19 X s 1886 407 m 9 Y s 1915 412 m 19 X s 1934 407 m 9 Y s 1910 420 m 6 Y s
+ 1906 426 m 9 X s 1910 403 m -5 Y s 1906 398 m 9 X s 1953 466 m -19 X s 1934 462 m 8 Y s 1963 466 m 19 X s 1982 462 m 8 Y s 2001 430 m -19 X s 1982 425 m 9 Y s 2011 430 m 19 X s 2030 425 m 9 Y s 2006 438 m 28 Y s 2002 466 m 9 X s 2006 421 m -28 Y s
+ 2002 393 m 9 X s 2097 464 m -19 X s 2078 460 m 8 Y s 2108 464 m 18 X s 2126 460 m 8 Y s 2145 473 m -19 X s 2126 469 m 9 Y s 2156 473 m 18 X s 2174 469 m 9 Y s 2193 468 m -19 X s 2174 464 m 8 Y s 2204 468 m 18 X s 2222 464 m 8 Y s 0 1 0 c 229 459 277
+ 468 325 458 374 458 422 467 470 466 518 457 566 466 614 461 662 459 710 461 758 468 806 455 854 462 902 467 950 460 998 464 1046 465 1094 469 1142 459 1190 473 1238 470 1286 463 1334 484 1382 464 1430 471 1478 460 1526 468 1574 463 1622 456 1670 469
+ 31 { m22} R 1766 464 m22 1862 456 1910 412 1958 466 2006 430 4 { m22} R 2102 464 2150 473 2198 468 3 { m22} R black 224 451 m -19 X s 205 447 m 9 Y s 235 451 m 18 X s 253 447 m 9 Y s 272 460 m -19 X s 253 456 m 8 Y s 283 460 m 18 X s 301 456 m 8 Y s
+ 320 466 m -19 X s 301 461 m 9 Y s 331 466 m 18 X s 349 461 m 9 Y s 368 457 m -19 X s 349 453 m 9 Y s 379 457 m 19 X s 398 453 m 9 Y s 416 463 m -18 X s 398 458 m 9 Y s 427 463 m 19 X s 446 458 m 9 Y s 464 465 m -18 X s 446 461 m 8 Y s 475 465 m 19 X
+ s 494 461 m 8 Y s 470 474 m s 465 474 m 9 X s 470 456 m s 465 456 m 9 X s 512 457 m -18 X s 494 453 m 8 Y s 523 457 m 19 X s 542 453 m 8 Y s 560 457 m -18 X s 542 453 m 9 Y s 571 457 m 19 X s 590 453 m 9 Y s 608 464 m -18 X s 590 460 m 8 Y s 619 464
+ m 19 X s 638 460 m 8 Y s 656 466 m -18 X s 638 462 m 8 Y s 667 466 m 19 X s 686 462 m 8 Y s 704 463 m -18 X s 686 458 m 9 Y s 715 463 m 19 X s 734 458 m 9 Y s 752 465 m -18 X s 734 461 m 8 Y s 763 465 m 19 X s 782 461 m 8 Y s 801 455 m -19 X s 782
+ 450 m 9 Y s 811 455 m 19 X s 830 450 m 9 Y s 849 463 m -19 X s 830 458 m 9 Y s 859 463 m 19 X s 878 458 m 9 Y s 897 469 m -19 X s 878 465 m 9 Y s 907 469 m 19 X s 926 465 m 9 Y s 945 455 m -19 X s 926 451 m 8 Y s 955 455 m 19 X s 974 451 m 8 Y s 993
+ 459 m -19 X s 974 455 m 9 Y s 1003 459 m 19 X s 1022 455 m 9 Y s 1041 471 m -19 X s 1022 467 m 9 Y s 1051 471 m 19 X s 1070 467 m 9 Y s 1089 470 m -19 X s 1070 465 m 9 Y s 1099 470 m 19 X s 1118 465 m 9 Y s 1137 457 m -19 X s 1118 453 m 9 Y s 1147
+ 457 m 19 X s 1166 453 m 9 Y s 1185 469 m -19 X s 1166 465 m 8 Y s 1195 469 m 19 X s 1214 465 m 8 Y s 1233 468 m -19 X s 1214 464 m 8 Y s 1243 468 m 19 X s 1262 464 m 8 Y s 1281 466 m -19 X s 1262 462 m 8 Y s 1291 466 m 19 X s 1310 462 m 8 Y s 1329
+ 479 m -19 X s 1310 475 m 8 Y s 1339 479 m 19 X s 1358 475 m 8 Y s 1334 488 m 1 Y s 1330 489 m 8 X s 1334 470 m -1 Y s 1330 469 m 8 X s 1377 469 m -19 X s 1358 465 m 8 Y s 1387 469 m 19 X s 1406 465 m 8 Y s 1425 461 m -19 X s 1406 457 m 8 Y s 1435
+ 461 m 19 X s 1454 457 m 8 Y s 1473 462 m -19 X s 1454 458 m 9 Y s 1483 462 m 19 X s 1502 458 m 9 Y s 1521 457 m -19 X s 1502 452 m 9 Y s 1531 457 m 19 X s 1550 452 m 9 Y s 1569 468 m -19 X s 1550 463 m 9 Y s 1579 468 m 19 X s 1598 463 m 9 Y s 1617
+ 471 m -19 X s 1598 466 m 9 Y s 1627 471 m 19 X s 1646 466 m 9 Y s 1665 465 m -19 X s 1646 461 m 8 Y s 1675 465 m 19 X s 1694 461 m 8 Y s 1670 474 m 3 Y s 1666 477 m 8 X s 1670 456 m -2 Y s 1666 454 m 8 X s 1761 458 m -19 X s 1742 454 m 9 Y s 1771
+ 458 m 19 X s 1790 454 m 9 Y s 1857 453 m -19 X s 1838 449 m 9 Y s 1867 453 m 19 X s 1886 449 m 9 Y s 1905 446 m -19 X s 1886 442 m 9 Y s 1915 446 m 19 X s 1934 442 m 9 Y s 1910 455 m 15 Y s 1906 470 m 9 X s 1910 438 m -15 Y s 1906 423 m 9 X s 1953
+ 461 m -19 X s 1934 456 m 9 Y s 1963 461 m 19 X s 1982 456 m 9 Y s 2001 518 m -19 X s 1982 514 m 8 Y s 2011 518 m 19 X s 2030 514 m 8 Y s 2006 527 m 70 Y s 2002 597 m 9 X s 2006 509 m -70 Y s 2002 439 m 9 X s 2097 468 m -19 X s 2078 464 m 9 Y s 2108
+ 468 m 18 X s 2126 464 m 9 Y s 2145 474 m -19 X s 2126 469 m 9 Y s 2156 474 m 18 X s 2174 469 m 9 Y s 2193 474 m -19 X s 2174 470 m 8 Y s 2204 474 m 18 X s 2222 470 m 8 Y s 0 0 1 c 229 451 277 460 325 466 374 457 422 463 470 465 518 457 566 457 614
+ 464 662 466 710 463 758 465 806 455 854 463 902 469 950 455 998 459 1046 471 1094 470 1142 457 1190 469 1238 468 1286 466 1334 479 1382 469 1430 461 1478 462 1526 457 1574 468 1622 471 1670 465 31 { m23} R 1766 458 m23 1862 453 1910 446 1958 461
+ 2006 518 4 { m23} R 2102 468 2150 474 2198 474 3 { m23} R black 221 459 m -16 X s 205 454 m 9 Y s 238 459 m 15 X s 253 454 m 9 Y s 269 453 m -16 X s 253 449 m 8 Y s 286 453 m 15 X s 301 449 m 8 Y s 317 468 m -16 X s 301 464 m 8 Y s 334 468 m 15 X s
+ 349 464 m 8 Y s 365 457 m -16 X s 349 453 m 9 Y s 382 457 m 16 X s 398 453 m 9 Y s 413 464 m -15 X s 398 460 m 9 Y s 430 464 m 16 X s 446 460 m 9 Y s 461 459 m -15 X s 446 455 m 9 Y s 478 459 m 16 X s 494 455 m 9 Y s 470 468 m 1 Y s 465 469 m 9 X s
+ 470 451 m -1 Y s 465 450 m 9 X s 509 462 m -15 X s 494 458 m 8 Y s 526 462 m 16 X s 542 458 m 8 Y s 557 454 m -15 X s 542 450 m 8 Y s 574 454 m 16 X s 590 450 m 8 Y s 605 459 m -15 X s 590 455 m 8 Y s 622 459 m 16 X s 638 455 m 8 Y s 653 464 m -15 X
+ s 638 459 m 9 Y s 670 464 m 16 X s 686 459 m 9 Y s 701 456 m -15 X s 686 451 m 9 Y s 718 456 m 16 X s 734 451 m 9 Y s 749 458 m -15 X s 734 454 m 9 Y s 766 458 m 16 X s 782 454 m 9 Y s 797 456 m -15 X s 782 452 m 8 Y s 814 456 m 16 X s 830 452 m 8 Y
+ s 845 461 m -15 X s 830 457 m 8 Y s 862 461 m 16 X s 878 457 m 8 Y s 893 464 m -15 X s 878 459 m 9 Y s 910 464 m 16 X s 926 459 m 9 Y s 941 459 m -15 X s 926 455 m 9 Y s 958 459 m 16 X s 974 455 m 9 Y s 989 459 m -15 X s 974 454 m 9 Y s 1006 459 m
+ 16 X s 1022 454 m 9 Y s 1037 468 m -15 X s 1022 464 m 8 Y s 1054 468 m 16 X s 1070 464 m 8 Y s 1085 464 m -15 X s 1070 459 m 9 Y s 1103 464 m 15 X s 1118 459 m 9 Y s 1133 459 m -15 X s 1118 455 m 8 Y s 1151 459 m 15 X s 1166 455 m 8 Y s 1181 468 m
+ -15 X s 1166 464 m 9 Y s 1199 468 m 15 X s 1214 464 m 9 Y s 1229 457 m -15 X s 1214 453 m 8 Y s 1247 457 m 15 X s 1262 453 m 8 Y s 1277 467 m -15 X s 1262 462 m 9 Y s 1295 467 m 15 X s 1310 462 m 9 Y s 1325 474 m -15 X s 1310 469 m 9 Y s 1343 474 m
+ 15 X s 1358 469 m 9 Y s 1334 482 m 1 Y s 1330 483 m 8 X s 1334 465 m -1 Y s 1330 464 m 8 X s 1373 463 m -15 X s 1358 459 m 8 Y s 1391 463 m 15 X s 1406 459 m 8 Y s 1421 456 m -15 X s 1406 452 m 9 Y s 1439 456 m 15 X s 1454 452 m 9 Y s 1469 463 m -15
+ X s 1454 459 m 9 Y s 1487 463 m 15 X s 1502 459 m 9 Y s 1517 458 m -15 X s 1502 454 m 9 Y s 1535 458 m 15 X s 1550 454 m 9 Y s 1565 465 m -15 X s 1550 460 m 9 Y s 1583 465 m 15 X s 1598 460 m 9 Y s 1613 476 m -15 X s 1598 472 m 9 Y s 1631 476 m 15 X
+ s 1646 472 m 9 Y s 1661 446 m -15 X s 1646 442 m 9 Y s 1679 446 m 15 X s 1694 442 m 9 Y s 1670 455 m 1 Y s 1666 456 m 8 X s 1670 438 m -1 Y s 1666 437 m 8 X s 1757 457 m -15 X s 1742 453 m 9 Y s 1775 457 m 15 X s 1790 453 m 9 Y s 1854 464 m -16 X s
+ 1838 460 m 8 Y s 1871 464 m 15 X s 1886 460 m 8 Y s 1902 489 m -16 X s 1886 485 m 9 Y s 1919 489 m 15 X s 1934 485 m 9 Y s 1910 498 m 17 Y s 1906 515 m 9 X s 1910 481 m -17 Y s 1906 464 m 9 X s 1950 482 m -16 X s 1934 477 m 9 Y s 1967 482 m 15 X s
+ 1982 477 m 9 Y s 1998 487 m -16 X s 1982 482 m 9 Y s 2015 487 m 15 X s 2030 482 m 9 Y s 2006 495 m 53 Y s 2002 548 m 9 X s 2006 478 m -53 Y s 2002 425 m 9 X s 2094 479 m -16 X s 2078 475 m 8 Y s 2111 479 m 15 X s 2126 475 m 8 Y s 2142 486 m -16 X s
+ 2126 481 m 9 Y s 2159 486 m 15 X s 2174 481 m 9 Y s 2190 483 m -16 X s 2174 479 m 9 Y s 2207 483 m 15 X s 2222 479 m 9 Y s 1 1 0 c 229 459 277 453 325 468 374 457 422 464 470 459 518 462 566 454 614 459 662 464 710 456 758 458 806 456 854 461 902
+ 464 950 459 998 459 1046 468 1094 464 1142 459 1190 468 1238 457 1286 467 1334 474 1382 463 1430 456 1478 463 1526 458 1574 465 1622 476 1670 446 31 { m24} R 1766 457 m24 1862 464 1910 489 1958 482 2006 487 4 { m24} R 2102 479 2150 486 2198 483 3 {
+ m24} R black 221 460 m -16 X s 205 456 m 9 Y s 238 460 m 15 X s 253 456 m 9 Y s 269 459 m -16 X s 253 454 m 9 Y s 286 459 m 15 X s 301 454 m 9 Y s 317 462 m -16 X s 301 457 m 9 Y s 334 462 m 15 X s 349 457 m 9 Y s 365 457 m -16 X s 349 453 m 8 Y s
+ 382 457 m 16 X s 398 453 m 8 Y s 413 459 m -15 X s 398 454 m 9 Y s 430 459 m 16 X s 446 454 m 9 Y s 461 459 m -15 X s 446 455 m 9 Y s 478 459 m 16 X s 494 455 m 9 Y s 470 468 m s 465 468 m 9 X s 470 451 m -1 Y s 465 450 m 9 X s 509 461 m -15 X s 494
+ 457 m 8 Y s 526 461 m 16 X s 542 457 m 8 Y s 557 467 m -15 X s 542 462 m 9 Y s 574 467 m 16 X s 590 462 m 9 Y s 605 453 m -15 X s 590 449 m 9 Y s 622 453 m 16 X s 638 449 m 9 Y s 653 459 m -15 X s 638 455 m 9 Y s 670 459 m 16 X s 686 455 m 9 Y s 701
+ 451 m -15 X s 686 446 m 9 Y s 718 451 m 16 X s 734 446 m 9 Y s 749 459 m -15 X s 734 455 m 8 Y s 766 459 m 16 X s 782 455 m 8 Y s 797 457 m -15 X s 782 453 m 8 Y s 814 457 m 16 X s 830 453 m 8 Y s 845 458 m -15 X s 830 453 m 9 Y s 862 458 m 16 X s
+ 878 453 m 9 Y s 893 456 m -15 X s 878 451 m 9 Y s 910 456 m 16 X s 926 451 m 9 Y s 941 464 m -15 X s 926 459 m 9 Y s 958 464 m 16 X s 974 459 m 9 Y s 989 460 m -15 X s 974 455 m 9 Y s 1006 460 m 16 X s 1022 455 m 9 Y s 1037 468 m -15 X s 1022 463 m
+ 9 Y s 1054 468 m 16 X s 1070 463 m 9 Y s 1085 463 m -15 X s 1070 459 m 8 Y s 1103 463 m 15 X s 1118 459 m 8 Y s 1133 466 m -15 X s 1118 462 m 9 Y s 1151 466 m 15 X s 1166 462 m 9 Y s 1181 467 m -15 X s 1166 462 m 9 Y s 1199 467 m 15 X s 1214 462 m 9
+ Y s 1229 449 m -15 X s 1214 444 m 9 Y s 1247 449 m 15 X s 1262 444 m 9 Y s 1277 471 m -15 X s 1262 467 m 8 Y s 1295 471 m 15 X s 1310 467 m 8 Y s 1286 480 m s 1282 480 m 8 X s 1286 462 m s 1282 462 m 8 X s 1325 474 m -15 X s 1310 470 m 8 Y s 1343
+ 474 m 15 X s 1358 470 m 8 Y s 1334 483 m 1 Y s 1330 484 m 8 X s 1334 465 m -1 Y s 1330 464 m 8 X s 1373 466 m -15 X s 1358 462 m 8 Y s 1391 466 m 15 X s 1406 462 m 8 Y s 1421 463 m -15 X s 1406 458 m 9 Y s 1439 463 m 15 X s 1454 458 m 9 Y s 1469 463
+ m -15 X s 1454 459 m 9 Y s 1487 463 m 15 X s 1502 459 m 9 Y s 1517 460 m -15 X s 1502 456 m 9 Y s 1535 460 m 15 X s 1550 456 m 9 Y s 1565 464 m -15 X s 1550 460 m 8 Y s 1583 464 m 15 X s 1598 460 m 8 Y s 1613 454 m -15 X s 1598 450 m 9 Y s 1631 454
+ m 15 X s 1646 450 m 9 Y s 1661 459 m -15 X s 1646 454 m 9 Y s 1679 459 m 15 X s 1694 454 m 9 Y s 1670 467 m 3 Y s 1666 470 m 8 X s 1670 450 m -2 Y s 1666 448 m 8 X s 1757 458 m -15 X s 1742 453 m 9 Y s 1775 458 m 15 X s 1790 453 m 9 Y s 1854 463 m
+ -16 X s 1838 458 m 9 Y s 1871 463 m 15 X s 1886 458 m 9 Y s 1902 500 m -16 X s 1886 496 m 9 Y s 1919 500 m 15 X s 1934 496 m 9 Y s 1910 509 m 19 Y s 1906 528 m 9 X s 1910 492 m -19 Y s 1906 473 m 9 X s 1950 486 m -16 X s 1934 481 m 9 Y s 1967 486 m
+ 15 X s 1982 481 m 9 Y s 1998 427 m -16 X s 1982 423 m 9 Y s 2015 427 m 15 X s 2030 423 m 9 Y s 2006 436 m 26 Y s 2002 462 m 9 X s 2006 419 m -26 Y s 2002 393 m 9 X s 2094 474 m -16 X s 2078 470 m 8 Y s 2111 474 m 15 X s 2126 470 m 8 Y s 2142 479 m
+ -16 X s 2126 475 m 8 Y s 2159 479 m 15 X s 2174 475 m 8 Y s 2190 474 m -16 X s 2174 470 m 8 Y s 2207 474 m 15 X s 2222 470 m 8 Y s 1 0 1 c 229 460 277 459 325 462 374 457 422 459 470 459 518 461 566 467 614 453 662 459 710 451 758 459 806 457 854
+ 458 902 456 950 464 998 460 1046 468 1094 463 1142 466 1190 467 1238 449 1286 471 1334 474 1382 466 1430 463 1478 463 1526 460 1574 464 1622 454 1670 459 31 { m25} R 1766 458 m25 1862 463 1910 500 1958 486 2006 427 4 { m25} R 2102 474 2150 479 2198
+ 474 3 { m25} R black 224 458 m -19 X s 205 454 m 9 Y s 235 458 m 18 X s 253 454 m 9 Y s 272 462 m -19 X s 253 457 m 9 Y s 283 462 m 18 X s 301 457 m 9 Y s 320 461 m -19 X s 301 456 m 9 Y s 331 461 m 18 X s 349 456 m 9 Y s 368 456 m -19 X s 349 452 m
+ 9 Y s 379 456 m 19 X s 398 452 m 9 Y s 416 454 m -18 X s 398 450 m 9 Y s 427 454 m 19 X s 446 450 m 9 Y s 464 460 m -18 X s 446 456 m 9 Y s 475 460 m 19 X s 494 456 m 9 Y s 512 457 m -18 X s 494 452 m 9 Y s 523 457 m 19 X s 542 452 m 9 Y s 560 464 m
+ -18 X s 542 460 m 8 Y s 571 464 m 19 X s 590 460 m 8 Y s 608 455 m -18 X s 590 450 m 9 Y s 619 455 m 19 X s 638 450 m 9 Y s 656 459 m -18 X s 638 454 m 9 Y s 667 459 m 19 X s 686 454 m 9 Y s 704 456 m -18 X s 686 451 m 9 Y s 715 456 m 19 X s 734 451
+ m 9 Y s 752 463 m -18 X s 734 458 m 9 Y s 763 463 m 19 X s 782 458 m 9 Y s 801 458 m -19 X s 782 454 m 9 Y s 811 458 m 19 X s 830 454 m 9 Y s 849 459 m -19 X s 830 455 m 9 Y s 859 459 m 19 X s 878 455 m 9 Y s 897 461 m -19 X s 878 456 m 9 Y s 907
+ 461 m 19 X s 926 456 m 9 Y s 945 465 m -19 X s 926 461 m 8 Y s 955 465 m 19 X s 974 461 m 8 Y s 993 457 m -19 X s 974 453 m 9 Y s 1003 457 m 19 X s 1022 453 m 9 Y s 1041 468 m -19 X s 1022 464 m 9 Y s 1051 468 m 19 X s 1070 464 m 9 Y s 1089 464 m
+ -19 X s 1070 460 m 8 Y s 1099 464 m 19 X s 1118 460 m 8 Y s 1137 467 m -19 X s 1118 463 m 8 Y s 1147 467 m 19 X s 1166 463 m 8 Y s 1185 463 m -19 X s 1166 459 m 8 Y s 1195 463 m 19 X s 1214 459 m 8 Y s 1233 456 m -19 X s 1214 452 m 9 Y s 1243 456 m
+ 19 X s 1262 452 m 9 Y s 1281 465 m -19 X s 1262 461 m 8 Y s 1291 465 m 19 X s 1310 461 m 8 Y s 1286 474 m s 1282 474 m 8 X s 1286 456 m s 1282 456 m 8 X s 1329 476 m -19 X s 1310 472 m 9 Y s 1339 476 m 19 X s 1358 472 m 9 Y s 1334 485 m 3 Y s 1330
+ 488 m 8 X s 1334 468 m -3 Y s 1330 465 m 8 X s 1377 462 m -19 X s 1358 458 m 9 Y s 1387 462 m 19 X s 1406 458 m 9 Y s 1425 469 m -19 X s 1406 465 m 9 Y s 1435 469 m 19 X s 1454 465 m 9 Y s 1473 460 m -19 X s 1454 456 m 9 Y s 1483 460 m 19 X s 1502
+ 456 m 9 Y s 1521 459 m -19 X s 1502 454 m 9 Y s 1531 459 m 19 X s 1550 454 m 9 Y s 1569 467 m -19 X s 1550 463 m 9 Y s 1579 467 m 19 X s 1598 463 m 9 Y s 1617 457 m -19 X s 1598 452 m 9 Y s 1627 457 m 19 X s 1646 452 m 9 Y s 1665 462 m -19 X s 1646
+ 458 m 8 Y s 1675 462 m 19 X s 1694 458 m 8 Y s 1670 471 m 2 Y s 1666 473 m 8 X s 1670 453 m -2 Y s 1666 451 m 8 X s 1761 458 m -19 X s 1742 454 m 9 Y s 1771 458 m 19 X s 1790 454 m 9 Y s 1857 463 m -19 X s 1838 459 m 8 Y s 1867 463 m 19 X s 1886 459
+ m 8 Y s 1905 502 m -19 X s 1886 498 m 8 Y s 1915 502 m 19 X s 1934 498 m 8 Y s 1910 511 m 19 Y s 1906 530 m 9 X s 1910 493 m -20 Y s 1906 473 m 9 X s 1953 470 m -19 X s 1934 466 m 8 Y s 1963 470 m 19 X s 1982 466 m 8 Y s 2001 562 m -19 X s 1982 558
+ m 9 Y s 2011 562 m 19 X s 2030 558 m 9 Y s 2006 571 m 106 Y s 2002 677 m 9 X s 2006 554 m -106 Y s 2002 448 m 9 X s 2097 473 m -19 X s 2078 469 m 9 Y s 2108 473 m 18 X s 2126 469 m 9 Y s 2145 477 m -19 X s 2126 473 m 8 Y s 2156 477 m 18 X s 2174 473
+ m 8 Y s 2193 475 m -19 X s 2174 470 m 9 Y s 2204 475 m 18 X s 2222 470 m 9 Y s 0 1 1 c 229 458 277 462 325 461 374 456 422 454 470 460 518 457 566 464 614 455 662 459 710 456 758 463 806 458 854 459 902 461 950 465 998 457 1046 468 1094 464 1142 467
+ 1190 463 1238 456 1286 465 1334 476 1382 462 1430 469 1478 460 1526 459 1574 467 1622 457 1670 462 31 { m26} R 1766 458 m26 1862 463 1910 502 1958 470 2006 562 4 { m26} R 2102 473 2150 477 2198 475 3 { m26} R black 225 459 m -20 X s 205 455 m 9 Y s
+ 234 459 m 19 X s 253 455 m 9 Y s 273 456 m -20 X s 253 452 m 8 Y s 282 456 m 19 X s 301 452 m 8 Y s 321 461 m -20 X s 301 457 m 8 Y s 330 461 m 19 X s 349 457 m 8 Y s 369 456 m -20 X s 349 452 m 9 Y s 378 456 m 20 X s 398 452 m 9 Y s 417 462 m -19 X
+ s 398 458 m 9 Y s 426 462 m 20 X s 446 458 m 9 Y s 465 464 m -19 X s 446 459 m 9 Y s 474 464 m 20 X s 494 459 m 9 Y s 470 472 m 2 Y s 465 474 m 9 X s 470 455 m -2 Y s 465 453 m 9 X s 513 454 m -19 X s 494 450 m 9 Y s 522 454 m 20 X s 542 450 m 9 Y s
+ 561 460 m -19 X s 542 456 m 8 Y s 570 460 m 20 X s 590 456 m 8 Y s 609 459 m -19 X s 590 455 m 9 Y s 618 459 m 20 X s 638 455 m 9 Y s 657 453 m -19 X s 638 449 m 8 Y s 666 453 m 20 X s 686 449 m 8 Y s 705 461 m -19 X s 686 457 m 8 Y s 714 461 m 20 X
+ s 734 457 m 8 Y s 753 461 m -19 X s 734 457 m 9 Y s 762 461 m 20 X s 782 457 m 9 Y s 801 457 m -19 X s 782 453 m 9 Y s 810 457 m 20 X s 830 453 m 9 Y s 849 465 m -19 X s 830 461 m 9 Y s 858 465 m 20 X s 878 461 m 9 Y s 897 462 m -19 X s 878 458 m 9
+ Y s 906 462 m 20 X s 926 458 m 9 Y s 945 460 m -19 X s 926 456 m 9 Y s 954 460 m 20 X s 974 456 m 9 Y s 993 456 m -19 X s 974 452 m 9 Y s 1002 456 m 20 X s 1022 452 m 9 Y s 1041 465 m -19 X s 1022 461 m 8 Y s 1050 465 m 20 X s 1070 461 m 8 Y s 1090
+ 464 m -20 X s 1070 460 m 9 Y s 1098 464 m 20 X s 1118 460 m 9 Y s 1138 465 m -20 X s 1118 460 m 9 Y s 1146 465 m 20 X s 1166 460 m 9 Y s 1186 460 m -20 X s 1166 455 m 9 Y s 1194 460 m 20 X s 1214 455 m 9 Y s 1234 460 m -20 X s 1214 455 m 9 Y s 1242
+ 460 m 20 X s 1262 455 m 9 Y s 1282 464 m -20 X s 1262 460 m 8 Y s 1290 464 m 20 X s 1310 460 m 8 Y s 1330 490 m -20 X s 1310 486 m 8 Y s 1338 490 m 20 X s 1358 486 m 8 Y s 1334 499 m 4 Y s 1330 503 m 8 X s 1334 481 m -4 Y s 1330 477 m 8 X s 1378 452
+ m -20 X s 1358 448 m 9 Y s 1386 452 m 20 X s 1406 448 m 9 Y s 1426 470 m -20 X s 1406 466 m 9 Y s 1434 470 m 20 X s 1454 466 m 9 Y s 1474 458 m -20 X s 1454 453 m 9 Y s 1482 458 m 20 X s 1502 453 m 9 Y s 1522 461 m -20 X s 1502 456 m 9 Y s 1530 461
+ m 20 X s 1550 456 m 9 Y s 1570 474 m -20 X s 1550 470 m 9 Y s 1578 474 m 20 X s 1598 470 m 9 Y s 1618 468 m -20 X s 1598 464 m 8 Y s 1626 468 m 20 X s 1646 464 m 8 Y s 1666 459 m -20 X s 1646 455 m 8 Y s 1674 459 m 20 X s 1694 455 m 8 Y s 1670 468 m
+ 2 Y s 1666 470 m 8 X s 1670 450 m -2 Y s 1666 448 m 8 X s 1762 462 m -20 X s 1742 458 m 9 Y s 1770 462 m 20 X s 1790 458 m 9 Y s 1858 464 m -20 X s 1838 460 m 8 Y s 1867 464 m 19 X s 1886 460 m 8 Y s 1906 499 m -20 X s 1886 495 m 8 Y s 1915 499 m 19
+ X s 1934 495 m 8 Y s 1910 508 m 22 Y s 1906 530 m 9 X s 1910 490 m -21 Y s 1906 469 m 9 X s 1954 468 m -20 X s 1934 464 m 8 Y s 1963 468 m 19 X s 1982 464 m 8 Y s 2002 605 m -20 X s 1982 601 m 8 Y s 2011 605 m 19 X s 2030 601 m 8 Y s 2006 614 m 132
+ Y s 2002 746 m 9 X s 2006 596 m -132 Y s 2002 464 m 9 X s 2098 478 m -20 X s 2078 474 m 8 Y s 2107 478 m 19 X s 2126 474 m 8 Y s 2146 480 m -20 X s 2126 476 m 8 Y s 2155 480 m 19 X s 2174 476 m 8 Y s 2194 480 m -20 X s 2174 475 m 9 Y s 2203 480 m 19
+ X s 2222 475 m 9 Y s 0.35 0.83 0.33 c 229 459 277 456 325 461 374 456 422 462 470 464 518 454 566 460 614 459 662 453 710 461 758 461 806 457 854 465 902 462 950 460 998 456 1046 465 1094 464 1142 465 1190 460 1238 460 1286 464 1334 490 1382 452
+ 1430 470 1478 458 1526 461 1574 474 1622 468 1670 459 31 { m27} R 1766 462 m27 1862 464 1910 499 1958 468 2006 605 4 { m27} R 2102 478 2150 480 2198 480 3 { m27} R black 221 452 m -16 X s 205 448 m 8 Y s 238 452 m 15 X s 253 448 m 8 Y s 269 453 m
+ -16 X s 253 449 m 8 Y s 286 453 m 15 X s 301 449 m 8 Y s 317 456 m -16 X s 301 451 m 9 Y s 334 456 m 15 X s 349 451 m 9 Y s 365 455 m -16 X s 349 451 m 8 Y s 382 455 m 16 X s 398 451 m 8 Y s 413 459 m -15 X s 398 454 m 9 Y s 430 459 m 16 X s 446 454
+ m 9 Y s 461 465 m -15 X s 446 460 m 9 Y s 478 465 m 16 X s 494 460 m 9 Y s 470 473 m 2 Y s 465 475 m 9 X s 470 456 m -2 Y s 465 454 m 9 X s 509 451 m -15 X s 494 447 m 9 Y s 526 451 m 16 X s 542 447 m 9 Y s 557 467 m -15 X s 542 463 m 8 Y s 574 467
+ m 16 X s 590 463 m 8 Y s 605 458 m -15 X s 590 453 m 9 Y s 622 458 m 16 X s 638 453 m 9 Y s 653 467 m -15 X s 638 463 m 8 Y s 670 467 m 16 X s 686 463 m 8 Y s 701 464 m -15 X s 686 459 m 9 Y s 718 464 m 16 X s 734 459 m 9 Y s 749 455 m -15 X s 734
+ 451 m 9 Y s 766 455 m 16 X s 782 451 m 9 Y s 797 461 m -15 X s 782 457 m 9 Y s 814 461 m 16 X s 830 457 m 9 Y s 845 459 m -15 X s 830 455 m 9 Y s 862 459 m 16 X s 878 455 m 9 Y s 893 455 m -15 X s 878 450 m 9 Y s 910 455 m 16 X s 926 450 m 9 Y s 941
+ 454 m -15 X s 926 450 m 9 Y s 958 454 m 16 X s 974 450 m 9 Y s 989 457 m -15 X s 974 452 m 9 Y s 1006 457 m 16 X s 1022 452 m 9 Y s 1037 467 m -15 X s 1022 462 m 9 Y s 1054 467 m 16 X s 1070 462 m 9 Y s 1085 454 m -15 X s 1070 450 m 9 Y s 1103 454 m
+ 15 X s 1118 450 m 9 Y s 1133 468 m -15 X s 1118 464 m 9 Y s 1151 468 m 15 X s 1166 464 m 9 Y s 1181 457 m -15 X s 1166 452 m 9 Y s 1199 457 m 15 X s 1214 452 m 9 Y s 1229 460 m -15 X s 1214 456 m 9 Y s 1247 460 m 15 X s 1262 456 m 9 Y s 1277 467 m
+ -15 X s 1262 463 m 8 Y s 1295 467 m 15 X s 1310 463 m 8 Y s 1325 482 m -15 X s 1310 478 m 8 Y s 1343 482 m 15 X s 1358 478 m 8 Y s 1334 491 m 2 Y s 1330 493 m 8 X s 1334 473 m -3 Y s 1330 470 m 8 X s 1373 458 m -15 X s 1358 453 m 9 Y s 1391 458 m 15
+ X s 1406 453 m 9 Y s 1421 465 m -15 X s 1406 461 m 8 Y s 1439 465 m 15 X s 1454 461 m 8 Y s 1469 456 m -15 X s 1454 452 m 8 Y s 1487 456 m 15 X s 1502 452 m 8 Y s 1517 463 m -15 X s 1502 459 m 8 Y s 1535 463 m 15 X s 1550 459 m 8 Y s 1565 475 m -15
+ X s 1550 470 m 9 Y s 1583 475 m 15 X s 1598 470 m 9 Y s 1613 459 m -15 X s 1598 455 m 8 Y s 1631 459 m 15 X s 1646 455 m 8 Y s 1661 478 m -15 X s 1646 474 m 8 Y s 1679 478 m 15 X s 1694 474 m 8 Y s 1670 487 m 4 Y s 1666 491 m 8 X s 1670 469 m -4 Y s
+ 1666 465 m 8 X s 1757 461 m -15 X s 1742 457 m 8 Y s 1775 461 m 15 X s 1790 457 m 8 Y s 1854 463 m -16 X s 1838 459 m 9 Y s 1871 463 m 15 X s 1886 459 m 9 Y s 1902 484 m -16 X s 1886 479 m 9 Y s 1919 484 m 15 X s 1934 479 m 9 Y s 1910 492 m 20 Y s
+ 1906 512 m 9 X s 1910 475 m -19 Y s 1906 456 m 9 X s 1950 476 m -16 X s 1934 472 m 9 Y s 1967 476 m 15 X s 1982 472 m 9 Y s 1998 444 m -16 X s 1982 439 m 9 Y s 2015 444 m 15 X s 2030 439 m 9 Y s 2006 452 m 41 Y s 2002 493 m 9 X s 2006 435 m -41 Y s
+ 2002 394 m 9 X s 2094 480 m -16 X s 2078 476 m 8 Y s 2111 480 m 15 X s 2126 476 m 8 Y s 2142 481 m -16 X s 2126 477 m 9 Y s 2159 481 m 15 X s 2174 477 m 9 Y s 2190 476 m -16 X s 2174 472 m 8 Y s 2207 476 m 15 X s 2222 472 m 8 Y s 0.35 0.33 0.85 c
+ 229 452 277 453 325 456 374 455 422 459 470 465 518 451 566 467 614 458 662 467 710 464 758 455 806 461 854 459 902 455 950 454 998 457 1046 467 1094 454 1142 468 1190 457 1238 460 1286 467 1334 482 1382 458 1430 465 1478 456 1526 463 1574 475 1622
+ 459 1670 478 31 { m28} R 1766 461 m28 1862 463 1910 484 1958 476 2006 444 4 { m28} R 2102 480 2150 481 2198 476 3 { m28} R black 224 445 m -19 X s 205 440 m 9 Y s 235 445 m 18 X s 253 440 m 9 Y s 229 449 m s 225 449 m 9 X s 229 440 m s 225 440 m 9 X
+ s 272 457 m -19 X s 253 453 m 8 Y s 283 457 m 18 X s 301 453 m 8 Y s 277 461 m s 273 461 m 9 X s 277 453 m s 273 453 m 9 X s 320 460 m -19 X s 301 456 m 8 Y s 331 460 m 18 X s 349 456 m 8 Y s 325 464 m 1 Y s 321 465 m 9 X s 325 456 m -1 Y s 321 455
+ m 9 X s 368 452 m -19 X s 349 448 m 8 Y s 379 452 m 19 X s 398 448 m 8 Y s 416 451 m -18 X s 398 447 m 9 Y s 427 451 m 19 X s 446 447 m 9 Y s 464 450 m -18 X s 446 445 m 9 Y s 475 450 m 19 X s 494 445 m 9 Y s 470 454 m 4 Y s 465 458 m 9 X s 470 445
+ m -4 Y s 465 441 m 9 X s 512 452 m -18 X s 494 448 m 9 Y s 523 452 m 19 X s 542 448 m 9 Y s 560 464 m -18 X s 542 460 m 9 Y s 571 464 m 19 X s 590 460 m 9 Y s 566 469 m s 561 469 m 9 X s 566 460 m s 561 460 m 9 X s 608 457 m -18 X s 590 453 m 8 Y s
+ 619 457 m 19 X s 638 453 m 8 Y s 614 461 m s 609 461 m 9 X s 614 453 m s 609 453 m 9 X s 656 469 m -18 X s 638 465 m 9 Y s 667 469 m 19 X s 686 465 m 9 Y s 662 474 m 4 Y s 657 478 m 9 X s 662 465 m -4 Y s 657 461 m 9 X s 704 461 m -18 X s 686 457 m
+ 8 Y s 715 461 m 19 X s 734 457 m 8 Y s 752 455 m -18 X s 734 450 m 9 Y s 763 455 m 19 X s 782 450 m 9 Y s 801 461 m -19 X s 782 457 m 8 Y s 811 461 m 19 X s 830 457 m 8 Y s 849 450 m -19 X s 830 445 m 9 Y s 859 450 m 19 X s 878 445 m 9 Y s 897 457 m
+ -19 X s 878 453 m 8 Y s 907 457 m 19 X s 926 453 m 8 Y s 945 466 m -19 X s 926 461 m 9 Y s 955 466 m 19 X s 974 461 m 9 Y s 950 470 m 1 Y s 945 471 m 9 X s 950 461 m -1 Y s 945 460 m 9 X s 993 457 m -19 X s 974 452 m 9 Y s 1003 457 m 19 X s 1022 452
+ m 9 Y s 1041 467 m -19 X s 1022 462 m 9 Y s 1051 467 m 19 X s 1070 462 m 9 Y s 1046 471 m 1 Y s 1041 472 m 9 X s 1046 462 m s 1041 462 m 9 X s 1089 459 m -19 X s 1070 455 m 8 Y s 1099 459 m 19 X s 1118 455 m 8 Y s 1094 463 m 1 Y s 1090 464 m 8 X s
+ 1094 455 m s 1090 455 m 8 X s 1137 456 m -19 X s 1118 452 m 8 Y s 1147 456 m 19 X s 1166 452 m 8 Y s 1142 460 m 2 Y s 1138 462 m 8 X s 1142 452 m -2 Y s 1138 450 m 8 X s 1185 464 m -19 X s 1166 460 m 8 Y s 1195 464 m 19 X s 1214 460 m 8 Y s 1233 463
+ m -19 X s 1214 459 m 8 Y s 1243 463 m 19 X s 1262 459 m 8 Y s 1238 467 m 2 Y s 1234 469 m 8 X s 1238 459 m -2 Y s 1234 457 m 8 X s 1281 454 m -19 X s 1262 450 m 9 Y s 1291 454 m 19 X s 1310 450 m 9 Y s 1286 459 m 3 Y s 1282 462 m 8 X s 1286 450 m -3
+ Y s 1282 447 m 8 X s 1329 461 m -19 X s 1310 457 m 9 Y s 1339 461 m 19 X s 1358 457 m 9 Y s 1334 466 m 5 Y s 1330 471 m 8 X s 1334 457 m -5 Y s 1330 452 m 8 X s 1377 457 m -19 X s 1358 452 m 9 Y s 1387 457 m 19 X s 1406 452 m 9 Y s 1382 461 m 3 Y s
+ 1378 464 m 8 X s 1382 452 m -2 Y s 1378 450 m 8 X s 1425 451 m -19 X s 1406 446 m 9 Y s 1435 451 m 19 X s 1454 446 m 9 Y s 1430 455 m 2 Y s 1426 457 m 8 X s 1430 446 m -2 Y s 1426 444 m 8 X s 1473 459 m -19 X s 1454 454 m 9 Y s 1483 459 m 19 X s
+ 1502 454 m 9 Y s 1521 471 m -19 X s 1502 466 m 9 Y s 1531 471 m 19 X s 1550 466 m 9 Y s 1526 475 m 4 Y s 1522 479 m 8 X s 1526 466 m -4 Y s 1522 462 m 8 X s 1569 468 m -19 X s 1550 464 m 8 Y s 1579 468 m 19 X s 1598 464 m 8 Y s 1574 472 m 2 Y s 1570
+ 474 m 8 X s 1574 464 m -2 Y s 1570 462 m 8 X s 1617 461 m -19 X s 1598 457 m 8 Y s 1627 461 m 19 X s 1646 457 m 8 Y s 1622 465 m 3 Y s 1618 468 m 8 X s 1622 457 m -3 Y s 1618 454 m 8 X s 1665 489 m -19 X s 1646 484 m 9 Y s 1675 489 m 19 X s 1694 484
+ m 9 Y s 1670 493 m 9 Y s 1666 502 m 8 X s 1670 484 m -9 Y s 1666 475 m 8 X s 1761 460 m -19 X s 1742 456 m 8 Y s 1771 460 m 19 X s 1790 456 m 8 Y s 1857 460 m -19 X s 1838 455 m 9 Y s 1867 460 m 19 X s 1886 455 m 9 Y s 1905 462 m -19 X s 1886 458 m
+ 8 Y s 1915 462 m 19 X s 1934 458 m 8 Y s 1910 466 m 19 Y s 1906 485 m 9 X s 1910 458 m -19 Y s 1906 439 m 9 X s 1953 468 m -19 X s 1934 464 m 9 Y s 1963 468 m 19 X s 1982 464 m 9 Y s 1958 473 m 1 Y s 1954 474 m 9 X s 1958 464 m -2 Y s 1954 462 m 9 X
+ s 2001 514 m -19 X s 1982 509 m 9 Y s 2011 514 m 19 X s 2030 509 m 9 Y s 2006 518 m 58 Y s 2002 576 m 9 X s 2006 509 m -58 Y s 2002 451 m 9 X s 2097 474 m -19 X s 2078 470 m 9 Y s 2108 474 m 18 X s 2126 470 m 9 Y s 2102 479 m 2 Y s 2098 481 m 9 X s
+ 2102 470 m -3 Y s 2098 467 m 9 X s 2145 477 m -19 X s 2126 473 m 8 Y s 2156 477 m 18 X s 2174 473 m 8 Y s 2150 481 m 1 Y s 2146 482 m 9 X s 2150 473 m -1 Y s 2146 472 m 9 X s 2193 473 m -19 X s 2174 468 m 9 Y s 2204 473 m 18 X s 2222 468 m 9 Y s
+ 2198 477 m 1 Y s 2194 478 m 9 X s 2198 468 m -1 Y s 2194 467 m 9 X s 0 1 1 c 229 445 277 457 325 460 374 452 422 451 470 450 518 452 566 464 614 457 662 469 710 461 758 455 806 461 854 450 902 457 950 466 998 457 1046 467 1094 459 1142 456 1190 464
+ 1238 463 1286 454 1334 461 1382 457 1430 451 1478 459 1526 471 1574 468 1622 461 1670 489 31 { m29} R 1766 460 m29 1862 460 1910 462 1958 468 2006 514 4 { m29} R 2102 474 2150 477 2198 473 3 { m29} R black 224 446 m -19 X s 205 441 m 9 Y s 235 446 m
+ 18 X s 253 441 m 9 Y s 229 451 m 1 Y s 225 452 m 9 X s 229 440 m -1 Y s 225 439 m 9 X s 272 460 m -19 X s 253 456 m 9 Y s 283 460 m 18 X s 301 456 m 9 Y s 277 465 m 1 Y s 273 466 m 9 X s 277 455 m -1 Y s 273 454 m 9 X s 320 464 m -19 X s 301 459 m 9
+ Y s 331 464 m 18 X s 349 459 m 9 Y s 325 469 m 1 Y s 321 470 m 9 X s 325 459 m -2 Y s 321 457 m 9 X s 368 447 m -19 X s 349 443 m 8 Y s 379 447 m 19 X s 398 443 m 8 Y s 416 447 m -18 X s 398 442 m 9 Y s 427 447 m 19 X s 446 442 m 9 Y s 422 452 m s
+ 417 452 m 9 X s 422 441 m s 417 441 m 9 X s 464 441 m -18 X s 446 437 m 9 Y s 475 441 m 19 X s 494 437 m 9 Y s 470 447 m 6 Y s 465 453 m 9 X s 470 436 m -6 Y s 465 430 m 9 X s 512 455 m -18 X s 494 451 m 8 Y s 523 455 m 19 X s 542 451 m 8 Y s 518
+ 460 m 1 Y s 513 461 m 9 X s 518 450 m -1 Y s 513 449 m 9 X s 560 459 m -18 X s 542 455 m 8 Y s 571 459 m 19 X s 590 455 m 8 Y s 566 464 m 2 Y s 561 466 m 9 X s 566 454 m -2 Y s 561 452 m 9 X s 608 459 m -18 X s 590 455 m 9 Y s 619 459 m 19 X s 638
+ 455 m 9 Y s 614 465 m s 609 465 m 9 X s 614 454 m -1 Y s 609 453 m 9 X s 656 453 m -18 X s 638 448 m 9 Y s 667 453 m 19 X s 686 448 m 9 Y s 662 458 m 5 Y s 657 463 m 9 X s 662 447 m -5 Y s 657 442 m 9 X s 704 455 m -18 X s 686 451 m 8 Y s 715 455 m
+ 19 X s 734 451 m 8 Y s 752 456 m -18 X s 734 452 m 8 Y s 763 456 m 19 X s 782 452 m 8 Y s 758 461 m 1 Y s 753 462 m 9 X s 758 451 m -1 Y s 753 450 m 9 X s 801 458 m -19 X s 782 453 m 9 Y s 811 458 m 19 X s 830 453 m 9 Y s 806 463 m s 801 463 m 9 X s
+ 806 452 m s 801 452 m 9 X s 849 450 m -19 X s 830 446 m 8 Y s 859 450 m 19 X s 878 446 m 8 Y s 897 464 m -19 X s 878 459 m 9 Y s 907 464 m 19 X s 926 459 m 9 Y s 902 469 m 1 Y s 897 470 m 9 X s 902 458 m -1 Y s 897 457 m 9 X s 945 473 m -19 X s 926
+ 469 m 8 Y s 955 473 m 19 X s 974 469 m 8 Y s 950 478 m 4 Y s 945 482 m 9 X s 950 468 m -3 Y s 945 465 m 9 X s 993 454 m -19 X s 974 449 m 9 Y s 1003 454 m 19 X s 1022 449 m 9 Y s 998 459 m 1 Y s 993 460 m 9 X s 998 448 m s 993 448 m 9 X s 1041 465 m
+ -19 X s 1022 460 m 9 Y s 1051 465 m 19 X s 1070 460 m 9 Y s 1046 470 m 2 Y s 1041 472 m 9 X s 1046 460 m -2 Y s 1041 458 m 9 X s 1089 478 m -19 X s 1070 474 m 8 Y s 1099 478 m 19 X s 1118 474 m 8 Y s 1094 483 m 3 Y s 1090 486 m 8 X s 1094 473 m -3 Y
+ s 1090 470 m 8 X s 1137 447 m -19 X s 1118 443 m 8 Y s 1147 447 m 19 X s 1166 443 m 8 Y s 1142 452 m 3 Y s 1138 455 m 8 X s 1142 442 m -3 Y s 1138 439 m 8 X s 1185 469 m -19 X s 1166 465 m 8 Y s 1195 469 m 19 X s 1214 465 m 8 Y s 1190 474 m 1 Y s
+ 1186 475 m 8 X s 1190 464 m -1 Y s 1186 463 m 8 X s 1233 453 m -19 X s 1214 448 m 9 Y s 1243 453 m 19 X s 1262 448 m 9 Y s 1238 458 m 3 Y s 1234 461 m 8 X s 1238 448 m -4 Y s 1234 444 m 8 X s 1281 451 m -19 X s 1262 446 m 9 Y s 1291 451 m 19 X s
+ 1310 446 m 9 Y s 1286 456 m 5 Y s 1282 461 m 8 X s 1286 446 m -6 Y s 1282 440 m 8 X s 1329 449 m -19 X s 1310 445 m 9 Y s 1339 449 m 19 X s 1358 445 m 9 Y s 1334 455 m 7 Y s 1330 462 m 8 X s 1334 444 m -7 Y s 1330 437 m 8 X s 1377 447 m -19 X s 1358
+ 443 m 9 Y s 1387 447 m 19 X s 1406 443 m 9 Y s 1382 452 m 5 Y s 1378 457 m 8 X s 1382 442 m -4 Y s 1378 438 m 8 X s 1425 443 m -19 X s 1406 438 m 9 Y s 1435 443 m 19 X s 1454 438 m 9 Y s 1430 448 m 3 Y s 1426 451 m 8 X s 1430 437 m -3 Y s 1426 434 m
+ 8 X s 1473 464 m -19 X s 1454 459 m 9 Y s 1483 464 m 19 X s 1502 459 m 9 Y s 1478 469 m 1 Y s 1474 470 m 8 X s 1478 459 m -1 Y s 1474 458 m 8 X s 1521 468 m -19 X s 1502 464 m 9 Y s 1531 468 m 19 X s 1550 464 m 9 Y s 1526 474 m 7 Y s 1522 481 m 8 X
+ s 1526 463 m -7 Y s 1522 456 m 8 X s 1569 468 m -19 X s 1550 463 m 9 Y s 1579 468 m 19 X s 1598 463 m 9 Y s 1574 473 m 4 Y s 1570 477 m 8 X s 1574 463 m -4 Y s 1570 459 m 8 X s 1617 469 m -19 X s 1598 465 m 8 Y s 1627 469 m 19 X s 1646 465 m 8 Y s
+ 1622 474 m 5 Y s 1618 479 m 8 X s 1622 464 m -5 Y s 1618 459 m 8 X s 1665 479 m -19 X s 1646 475 m 9 Y s 1675 479 m 19 X s 1694 475 m 9 Y s 1670 484 m 13 Y s 1666 497 m 8 X s 1670 474 m -13 Y s 1666 461 m 8 X s 1761 462 m -19 X s 1742 458 m 8 Y s
+ 1771 462 m 19 X s 1790 458 m 8 Y s 1766 467 m 1 Y s 1762 468 m 8 X s 1766 457 m -1 Y s 1762 456 m 8 X s 1857 457 m -19 X s 1838 453 m 8 Y s 1867 457 m 19 X s 1886 453 m 8 Y s 1905 437 m -19 X s 1886 432 m 9 Y s 1915 437 m 19 X s 1934 432 m 9 Y s
+ 1910 442 m 20 Y s 1906 462 m 9 X s 1910 432 m -20 Y s 1906 412 m 9 X s 1953 469 m -19 X s 1934 465 m 8 Y s 1963 469 m 19 X s 1982 465 m 8 Y s 1958 474 m 4 Y s 1954 478 m 9 X s 1958 464 m -4 Y s 1954 460 m 9 X s 2001 554 m -19 X s 1982 550 m 9 Y s
+ 2011 554 m 19 X s 2030 550 m 9 Y s 2006 559 m 92 Y s 2002 651 m 9 X s 2006 549 m -92 Y s 2002 457 m 9 X s 2097 467 m -19 X s 2078 462 m 9 Y s 2108 467 m 18 X s 2126 462 m 9 Y s 2102 472 m 4 Y s 2098 476 m 9 X s 2102 461 m -4 Y s 2098 457 m 9 X s
+ 2145 476 m -19 X s 2126 472 m 8 Y s 2156 476 m 18 X s 2174 472 m 8 Y s 2150 481 m 2 Y s 2146 483 m 9 X s 2150 471 m -2 Y s 2146 469 m 9 X s 2193 474 m -19 X s 2174 470 m 9 Y s 2204 474 m 18 X s 2222 470 m 9 Y s 2198 480 m 3 Y s 2194 483 m 9 X s 2198
+ 469 m -3 Y s 2194 466 m 9 X s 0.754 0.715 0.676 c 229 446 277 460 325 464 374 447 422 447 470 441 518 455 566 459 614 459 662 453 710 455 758 456 806 458 854 450 902 464 950 473 998 454 1046 465 1094 478 1142 447 1190 469 1238 453 1286 451 1334 449
+ 1382 447 1430 443 1478 464 1526 468 1574 468 1622 469 1670 479 31 { m30} R 1766 462 m30 1862 457 1910 437 1958 469 2006 554 4 { m30} R 2102 467 2150 476 2198 474 3 { m30} R black 1 1 1 c 794 269 907 646 bf black
+ gsave 2268 1076 0 0 C 1104.62 896.69 t 0 r /Helvetica findfont 19.4933 sf 0 0 m (<# of clusters>) show NC gr
+ gsave 2268 1076 0 0 C 1104.62 877.197 t 0 r /Helvetica findfont 19.4933 sf 0 0 m (det = EMCAL) show NC gr
+ gsave 2268 1076 0 0 C 1104.62 859.87 t 0 r /Helvetica findfont 19.4933 sf 0 0 m (average) show NC gr 1006 863 m20
+ gsave 2268 1076 0 0 C 1104.62 836.044 t 0 r /Helvetica findfont 19.4933 sf 0 0 m (SM_ 0) show NC gr 1 0 0 c 1006 843 m21 black
+ gsave 2268 1076 0 0 C 1104.62 816.551 t 0 r /Helvetica findfont 19.4933 sf 0 0 m (SM_ 1) show NC gr 0 1 0 c 1006 822 m22 black
+ gsave 2268 1076 0 0 C 1104.62 794.892 t 0 r /Helvetica findfont 19.4933 sf 0 0 m (SM_ 2) show NC gr 0 0 1 c 1006 801 m23 black
+ gsave 2268 1076 0 0 C 1104.62 775.399 t 0 r /Helvetica findfont 19.4933 sf 0 0 m (SM_ 3) show NC gr 1 1 0 c 1006 780 m24 black
+ gsave 2268 1076 0 0 C 1104.62 753.74 t 0 r /Helvetica findfont 19.4933 sf 0 0 m (SM_ 4) show NC gr 1 0 1 c 1006 760 m25 black
+ gsave 2268 1076 0 0 C 1104.62 734.246 t 0 r /Helvetica findfont 19.4933 sf 0 0 m (SM_ 5) show NC gr 0 1 1 c 1006 739 m26 black
+ gsave 2268 1076 0 0 C 1104.62 712.587 t 0 r /Helvetica findfont 19.4933 sf 0 0 m (SM_ 6) show NC gr 0.35 0.83 0.33 c 1006 718 m27 black
+ gsave 2268 1076 0 0 C 1104.62 693.094 t 0 r /Helvetica findfont 19.4933 sf 0 0 m (SM_ 7) show NC gr 0.35 0.33 0.85 c 1006 698 m28 black
+ gsave 2268 1076 0 0 C 1104.62 671.435 t 0 r /Helvetica findfont 19.4933 sf 0 0 m (SM_ 8) show NC gr 0 1 1 c 1006 677 m29 black
+ gsave 2268 1076 0 0 C 1104.62 649.776 t 0 r /Helvetica findfont 19.4933 sf 0 0 m (SM_ 9) show NC gr 0.754 0.715 0.676 c 1006 656 m30
+ gr gr showpage
+end
+%%EOF
\subsubsection{Time calibration - Marie }
-The time of the amplitude measured by a given cell is a good candidate to reject noisy towers, identify pile up events, or even identify heavy hadrons at low energy. The average time is around 650 ns. The aim of the time calibration is to move this mean value to 0, with as small spread as possible (negative values are unavoidable for the moment).
+The time of the amplitude measured by a given cell is a good candidate to reject noisy towers, identify pile up events when coming from different Bunch Crossing, or even identify heavy hadrons at low energy. The average time is around 580 ns. The aim of the time calibration is to do a relative calibration between cells to align all cells to a mean value of 0 ns, with as small spread as possible (negative values are unavoidable for the moment). The time calibration coefficient for each cell is the result of the average time of the cell when belonging to a cluster with enough energy (>1GeV).
+The calibration coefficient have to be subtracted to the cell time.
+
+\paragraph*{Time Calibration Procedure\\}
+
+Since the some variations of mean time have been observed depending on the bunch cross numbers (BC) $\%$ 4 the computation of the time coefficients is done for each bunch cross numbers BC$\%$ 4 scheme.
+
+The time calibration coefficient computing is done in 2 iterations.
+\begin{itemize}
+\item 1$^{rst}$ iteration:\\
+Get Bunch Cross Number for the event.
+Loop on all cluster of the event. \\ Loop on all cells in the cluster.
+If cell amplitude is > 0.9 GeV and 500ns < cell time < 700ns then
+compute average per cell per BC$\%$4 and fill in 1D histogam with calibration coefficients: hAveragesBC$x$ where $x$ stands for the result of BC$\%$4.
+\item 2$^{nd}$iteration:\\
+Get Bunch Cross Number for the event.\\
+Loop on all cluster of the event. \\ Loop on all cells in the cluster.\\
+Get Calibration coefficient for BC$\%$4 from hAveragesBC$x$. If cell amplitude is > 0.9 GeV and -20ns < (cell time-cell Calibration Coefficient) < 20ns.\\
+Compute average time per cell per BC$\%$4 and fill in 1D histogam with those new calibration coefficients: hAllAveragesBC$x$.
+\end{itemize}
+
+
+\paragraph*{Acces time calibration coefficients\\}
+
+The time calibration coefficient are stored in OADB in TH1D histograms named hAllTimeAvBCx where x stands for the BC$\%$4 value.
+They have the usual structure:
+runrange/pass/
+
+To use them in your analysis you may do the following:\\
+\begin{DDbox}{\linewidth}
+\begin{lstlisting}
+ AliOADBContainer *contTRF=new AliOADBContainer("");
+contTRF->InitFromFile(Form("%s/EMCALTimeCalib.root",fOADBFilePathEMCAL.Data()),"AliEMCALTimeCalib");
+TObjArray trecal=(TObjArray)contTRF->GetObject(runnumber);
+if(trecal)
+{TObjArray trecalpass=(TObjArray)trecal->FindObject(pass);
+if(trecalpass)
+{printf("AliCalorimeterUtils::SetOADBParameters() - Time Recalibrate EMCAL \n");
+for (Int_t ibc = 0; ibc < 4; ++ibc)
+{
+TH1F *h = GetEMCALChannelTimeRecalibrationFactors(ibc);
+if (h)
+delete h;
+h = (TH1F*)trecalpass->FindObject(Form("hAllTimeAvBC%d",ibc));
+}}
+\end{lstlisting}
+\end{DDbox}
+
+As already mentioned the time calibration of the cells is not done at the reconstruction level but offline during analysis. The methods to recalibrate cells is implemented in the class
+\texttt{\$ALICE\_ROOT/\\EMCAL/AliEMCALRecoUtils}.
+The method \texttt{RecalibrateCellTime(absId,bc,time)}: is called by the method \texttt{SwitchOnTimeRecalibration} , modifies the provided time with the calibration parameters. The inputs are the bunch crossing number that can be recovered from the event with InputEvent()->GetBunchCrossNumber(), and the absolute ID of the cell.
+The way to pass the calibration parameters to AliEMCALRecoUtils is the following.\\
+\begin{DDbox}{\linewidth}
+\begin{lstlisting}
+AliCalorimeterUtils *cu = new AliCalorimeterUtils ;
+TGeoManager::Import("geometry.root") ; //need file "geometry.root" in local dir!!!!
+AliEMCALRecoUtils * reco = cu->GetEMCALRecoUtils();
+and then
+for(Int_t i =0; i< 4; i++) reco->SetEMCALChannelTimeRecalibrationFactors( i, (TH1F*) file->Get(Form("hAllTimeAvBC%d",i)))
+\end{lstlisting}
+\end{DDbox}
+, where TFile *file is in the file containing the time calibration coefficients.
+
+Some more details on time recalibration can be found in twikis:
+https://twiki.cern.ch/twiki/bin/\\viewauth/ALICE/EMCalCode:HowTo\#Energy\_and\_Time\_calibration and in \\
+https://twiki.cern.ch/twiki/bin/viewauth/ALICE/EMCALTimeCalibration
\subsection{Alignment - Marco}
\subsection{Bad channel finding - Alexis}
-The analysis is done on the output of QA histograms. The idea is to check distributions over the cells of:
+The analysis is done on the output of offline Quality assurance see section \ref{sec:QAOffline} histograms \texttt{TH2F EMCAL\_hAmpId} containing the distribution of amplitudes (energy) of cell versus cell Absolute ID number (AbsId). The idea is to check distributions over the cells of different observables extracted from this histograms. Then each cell is tested regarding to the distribution over all the cells for each obbservable. The different tests are the following:
+
+\begin{enumerate}
+\item average energy (average computed for Emin< E <Emax)
+\item average number of hit per event (average computed for Emin< E <Emax)
+\item Shape criteria : A fit of the cell energy (amplitude) distribution is performed with the function:$A*e^{-B*x}/x^2$ between Emin and Emax. Then the $\chi^{2}/ndf$, A and B which are parameters from the fit of each cell amplitude.
+\end{enumerate}
+
+Each criteria is tested at least once (they can be tested also for different energy i nterval). At the end of each test the marked cells are excluded (if above nsigma from mean value, usually nsigma is taken equal 4 or 5) before computing the next distribution.
+
+The typical nsigma used is 4 or 5.
+The min energy considered is 0.1 GeV/0.3 GeV. And the maximum energy depends on the data (minbias or triggered data).
+
+
+There are different levels of classification for cells wich will enter the deadMap:
\begin{itemize}
-\item average energy (criteria 1) and
-\item average number of hit per event (criteria 2) (average computed for E > Emin)
-\item Shape criteria : $\chi^{2}/ndf$ (criteria 3), A (criteria 4) and B (criteria 5) which are parameters from the fit of each cell amplitude (the fit function is $A*e^{-B*x}/x^2$ and the fit range is from Emin to Emax).
+\item kAlive = 0: cell is OK
+\item kDead = 1: cell is dead
+\item kHot = 2: cell is bad
+\item kWarning=3: cell may have problems (warm or miscalibrated)
\end{itemize}
-Each criteria is ran once, at each step, and the marked cells are excluded (above nsigma from mean value) to compute the next distribution. \footnote{For each criteria we have some parameters Emin (min energy) Emax, (max energy for the Energy distribution fit), and nsigma, nb of sigma we use for excluding the cell;}
+The distinction of the bad/warm status is done by visual check of the energy distribution of the cells detected by the different tests described above.
-The typical nsigma used is 4 or 5.
-The min energy considered is 0.1 GeV -0.3 GeV. And the max energy for the fit depends on the data. Bad/warm channels are not detected automatically. The distinction is made by a visual check, so it is at some point subjective. (??????)
+In the reconstruction pass the only the cells marked as kHot and kDead are not reconstructed.
+
+The macro to compute the BadCellAnalysis can be found on \\ \texttt{\$ALICE\_ROOT/PWGGA/CaloTrackCorrelations/macros/QA/BadChannelAnalysis.C}. This macro allows first to merge all the individual run histograms in order to get the energy distribution over a large statistics. Usually this is done at the end of a data taking period. Then it computes the different tests and produces the text file with the dead and bad cells detected. Finally the macro produce a pdf file containing the individual energy distribution for all the detected cells.
+
+All the results and OCDB created are updated on the following twiki:\\
+https://twiki.cern.ch/twiki/bin/viewauth/ALICE/EMCalQABadChannels.
-The cells are then marked bad or warm and passed through OCDB, in the reconstruction pass, the bad ones are excluded.
git://git.uio.no / u / mrichter / AliRoot.git / commitdiff