Reconstruction & Analysis instruction for the data including FE65p2

Instruction : installation of normal EUTelescope

################################################################
# you can set any TBRECBASEDIR as you want.

export LANG=C
export LC_ALL=C
export TBRECBASEDIR=/Path/To/Your/TBREC/DIR/

################################################################

export INSTALLDIR=$TBRECBASEDIR/ilcsoft20170515
mkdir -p $INSTALLDIR
cd $INSTALLDIR/
export ILCSOFT=$PWD/ILCSOFT
mkdir -p $ILCSOFT
cd $ILCSOFT
svn co https://svnsrv.desy.de/public/ilctools/ilcinstall/branches/eutel-git ilcinstall_eutel-git
 

# modify following file to newer Eudaq and Eutelescope_version :

emacs -nw $ILCSOFT/ilcinstall_eutel-git/releases/v01-17-05/release-versions.py

# Eutelescope_version = "tags/v1.0"
# Eudaq_version = "tags/v1.5.1"

# fix root repository path

emacs -nw $ILCSOFT/ilcinstall_eutel-git/ilcsoft/root.py 

# self.download.url = 'https://root.cern.ch/download/root_v%s.source.tar.gz' % self.version

source /afs/cern.ch/sw/lcg/contrib/gcc/4.8/x86_64-slc6/setup.sh
cd $ILCSOFT/ilcinstall_eutel-git
./ilcsoft-install -i examples/eutelescope/release-standalone.cfg

# take 1hr
# when it successfully finished make script for the next login

cat > $TBRECBASEDIR/setupITKTB.sh <<EOF
export TB_ANALYSIS=$TBRECBASEDIR/Reconstruction
export TB_RECDIR=$TB_ANALYSIS/TBReconstruction/TBReconstruction-trunk
export ILCSOFT=$TBRECBASEDIR/ilcsoft20170515/ILCSOFT
source /afs/cern.ch/sw/lcg/contrib/gcc/4.8/x86_64-slc6/setup.sh
source \$ILCSOFT/v01-17-05/Eutelescope/v1.0/build_env.sh
EOF

##
## better to exit once and lonin again
## (in case you loged out once : do export TBRECBASEDIR=/Path/To/Your/TBREC/DIR/ again)

source $TBRECBASEDIR/setupITKTB.sh
export WORKDIR=$TBRECBASEDIR/work/run_cern_2016_nov_itk
mkdir -p $WORKDIR

See http://atlaspc5.kek.jp/do/view/Main/HowToInstallEUTelescope#x5_B9F_Reconstruction for test run.

YARR Converter and necessary geo files for FE65p2 Reconstruction

# copy files for source files

cp /afs/cern.ch/work/k/kojin/public/ITKTB/ilcsoft20170515/ILCSOFT/v01-17-05/Eutelescope/v1.0/external/eudaq/v1.5.1/main/lib/plugins/YARRConverterPlugin.cc $ILCSOFT/v01-17-05/Eutelescope/v1.0/external/eudaq/v1.5.1/main/lib/plugins/
cp /afs/cern.ch/work/k/kojin/public/ITKTB/ilcsoft20170515/ILCSOFT/v01-17-05/Eutelescope/v1.0/pixgeo/include/FE65p2.h $ILCSOFT/v01-17-05/Eutelescope/v1.0/pixgeo/include/
cp /afs/cern.ch/work/k/kojin/public/ITKTB/ilcsoft20170515/ILCSOFT/v01-17-05/Eutelescope/v1.0/pixgeo/src/FE65p2.cc $ILCSOFT/v01-17-05/Eutelescope/v1.0/pixgeo/src
cp /afs/cern.ch/work/k/kojin/public/ITKTB/ilcsoft20170515/ILCSOFT/v01-17-05/Eutelescope/v1.0/pixgeo/include/FEI4FourChipKEKTypeB.h $ILCSOFT/v01-17-05/Eutelescope/v1.0/pixgeo/include/
cp /afs/cern.ch/work/k/kojin/public/ITKTB/ilcsoft20170515/ILCSOFT/v01-17-05/Eutelescope/v1.0/pixgeo/src/FEI4FourChipKEKTypeB.cc   $ILCSOFT/v01-17-05/Eutelescope/v1.0/pixgeo/src

# copy config files

cd $WORKDIR
tar -zxvf /afs/cern.ch/work/k/kojin/public/ITKTB/work/run_cern_2016_nov_itk/config.tar.gz

# change some hard coding histogram range :

emacs -nw $ILCSOFT/v01-17-05/Eutelescope/v1.0/src/EUTelMille.cc

### l.3072-3074
# const int NBin = 600;
# const double Min = -300.;
# const double Max = 300.;

emacs -nw $ILCSOFT/v01-17-05/Eutelescope/v1.0/src/EUTelPreAlignment.cc

### l.129 133
# AIDA::IHistogram1D * histo1Da = AIDAProcessor::histogramFactory(this)->createHistogram1D( (basePath + tempHistoName).c_str(), 100 , -20., 20.);
# AIDA::IHistogram1D * histo1Db = AIDAProcessor::histogramFactory(this)->createHistogram1D( (basePath + tempHistoName).c_str(), 100 , -20., 20.);

## bit odd but this is necessary way to compile...

source compile.sh

## change your config_fe65p2.cfg files to have right BasePath

BasePath = $TBRECBASEDIR/work/run_%(TestBeamDate)s
## note $TBRECBASEDIR env have to be extracted since it doesn't work in the config file

## run converter as the first try

./run.sh -c config_fe65.cfg -t 11111 -r 8174 

# -t option correspond to the step executing
# converter/clustering/hitmaker/align/fitter
#(e.g)
# 10000 : converter only
# 11000 : converter + clustering
#

TBmon2 Analysis

You can find Ntuple.root in the working directory.

This correspond to the tbtrack file (Ntuple of TBmon2 input)

In case you skipped Reconstruction please use following file instead.

/home/kojin/work/PPSTB/PPSdata/v17/cern_2016_nov_itk/tbtrack/tbtrack008174.root

-- Koji Nakamura - 2017-05-22