Batreco
Marcin Byszewski
/afs/cern.ch/user/b/byszew/public/batreco
10/19/10
M.Byszewski - batreco
1
Two modes: Command line and Graphics
Raw data and alignment of detectors (on clusters positions)
Advice needed:
Root does not pass over (to Argc,Argv) command
line options that it itself recognizes.
The –h command line option never reaches my
program
$ ./batreco -H
Usage: batreco -c configfile
Options:
-b (--no-gui) run without gui
-c filename specify configuration file
-f filename specify root file name to analyze
-n N
limit number of events to analyze
-v N
verbosity 0-3, default 1
$./batreco -b -c mumegas.offline.myconfig_new_frame_diagram
batch mode. Gui disabled
Warning in <CBatConfig::KeywordValue>: Setting value of keyword 'unknown' to -1
Warning in <CBatConfig::KeywordValue>: Setting value of keyword 'unknown' to -1
Warning in <CBatConfig::KeywordValue>: Setting value of keyword 'unknown' to -1
Warning in <CBatConfig::ReadConfigFile>: Module 6: mmega.batmod6xy and single axis inversion are specified. Resulting in x: inverted
,y:inverted
Warning in <TClass::TClass>: no dictionary for class RunInfo is available
10/19/10
M.Byszewski - batreco
2
Configuration file
Based on mumegas.offline.myconfig_new_frame_diagram
############################################################
# PROJECT : MAMMA (Muon ATLAS Micromegas Analysis)
# PROGRAM : Analysis configuration file
# PURPOSE : To steer Micromegas/BAT telescope analysis
# AUTHOR : Venkatesh Kaushik ([email protected])
# MODIFIED : 2009 August 1
# MODIFIED : 2010 August 27 by Theodoros Alexopoulos <[email protected]>
# MODIFIED: 2010 October 4 by Marcin Byszewski ([email protected]): added Batreco section and
mmega.MuMegaXY key
############################################################
…
mmega.RunNumbers: 4786
…
############################ BATRECO - BAT RECONSTRUCTION #######################################
## cluster identification: charge tresholds
mmega.batreco.bat.treshold.qclusterStripMin: 20
mmega.batreco.bat.treshold.qclusterStripMax: -1
mmega.batreco.bat.treshold.qclusterTotalMin: 400
mmega.batreco.bat.treshold.qclusterTotalMax: -1
mmega.batreco.bat.treshold.qsingleStripMin: 100
mmega.batreco.bat.treshold.qsingleStripMax: -1
mmega.batreco.mm.treshold.tauMin: -1
mmega.batreco.mm.treshold.tauMax: -1
############################ END BATRECO PROPERTIES #############################################
10/19/10
M.Byszewski - batreco
3
Options tab in GUI mode
change values of filters
mmega.batreco.bat.treshold.qclusterStripMin: 20
mmega.batreco.bat.treshold.qclusterStripMax: -1
mmega.batreco.bat.treshold.qclusterTotalMin: 400
mmega.batreco.bat.treshold.qclusterTotalMax: -1
mmega.batreco.bat.treshold.qsingleStripMin: 100
mmega.batreco.bat.treshold.qsingleStripMax: -1
mmega.batreco.mm.treshold.tauMin: -1
mmega.batreco.mm.treshold.tauMax: -1
and redo analysis before saving data to the file
At present:
batreco saves temporary file after finishing reconstruction
original root file is replaced with the new version on exit
10/19/10
M.Byszewski - batreco
4
Configuration file
Added keys
###################### ALIGNMENT
#################################################
# Coordinate system: righhand coordinate
# Z downstream beam, Z=0 at last BAT detector
# see drawing below
## REFERNCE FRAME PARAMETERS (LEFT EDGE OF not inverted MOD6 is 0,0,0)
## Z position of BAT Mod1
mmega.BATMod1Z: 0, -41, -42, -814
## Was BAT Mod1 in TB
mmega.BATMod1In: true,true,true,true
## XY plane was upright/inverted
#mmega.BATMod1XY: upright,upright,upright,upright
## X,Y inversion defined separately (true/false) - optional to and overriding mmega.BATMod1XY
mmega.BATMod1XInv: false, false, false, true
mmega.BATMod1YInv: false, false, false, false
…
mmega.MuMega1Z:-1760,-1760, -1667, -1787
#mmega.MuMega1XY: upright,upright,upright,upright
mmega.MuMega1XInv: false, false, false, false
mmega.MuMega1YInv: false, false, false, false
mmega.MuMega1VH: unknown, unknown, unknown, vertical
## Run Range validity [includes boundaries]
mmega.RunRangeMin: 0, 2535, 3551, 4135
mmega.RunRangeMax: 2534,3550, 4134, 10000
Warning in <CBatConfig::ReadConfigFile>: Module 6: mmega.batmod6xy and single axis inversion are specified. Resulting in x: inverted , y:inverted
10/19/10
M.Byszewski - batreco
5
Output: SimpleFastReco root file
Alignment to RunInfo tree:
bat_align1x, 3,6
bat_align1y, 3,6
mm_align1x,
mm_align1y
other data to ntp tree:
Bat Clusters: (x, y data)
bat_nlcu1x,
-number of clusters (vector)
bat_clq1x
-cluster charge (vector)
bat_clsz1x,
-cluster size (vector)
bat_clpos1x
-cluster position (vector)
Tracks: (x, y data)
bat_ntracksx
bat_trackax
bat_trackbx
bat_trackmm1x
- number of tracks
- track fit: 1. parameter (offset) (vector)
- track fit 2. parameter (slope) (vector)
- track extrapolated position to micromega plane (vector)
At present:
batreco saves temporary file after finishing reconstruction
original root file is replaced with the new version on exit
10/19/10
M.Byszewski - batreco
6
GUI mode: Raw Data
Raw data and alignment of detectors (on clusters positions)
Bat 1
Bat 3
Bat 6
Black – raw data
Red – aligned
mm
Bat 6 taken as reference
10/19/10
M.Byszewski - batreco
7
Event Display (x axis)
Raw data (black line) and found clusters (filled bars). Variable x, and y scales
Bat 1
Bat 3
Bat 6
mm
10/19/10
M.Byszewski - batreco
8
Cluster Statistics (x axis)
Size, number, charge of clusters and eta distribution for position correction
Bat 1
Bat 3
Bat 6
10/19/10
M.Byszewski - batreco
9
Reconstruction 1 (x axis)
Comparison with mm data (clusters)
Raw Mm cluster q weighted position
Fitted track Extrapolated to
Mm plane
Position differences =
Mm cluster q-weighted position
–
Fitted track Extrapolated to
Mm plane
10/19/10
M.Byszewski - batreco
10
Reconstruction 2
Beam profile
Beam profile
(extrapolated tracks to mm)
10/19/10
M.Byszewski - batreco
Track angle vs position on mm plane
11
Thank you
Fitting procedure, track identification
per event , lsq fit to all combination of points
select the one with the smallest error
remove fitted points from further fitting
continue looking for tracks if more than three points
277
ijk=(0 0 0) x1=29.443 x2=29.3564 x3=31.0352 corr=0.251963 error=0.51075 0.00178584
ijk=(0 0 1) x1=29.443 x2=29.3564 x3=30.3301 corr=0.222541 error=0.296875 0.00103803
ijk=(0 0 2) x1=29.443 x2=29.3564 x3=28.9253 corr=0.443144 error=0.129278 0.00045202
ijk=(0 0 3) x1=29.443 x2=29.3564 x3=28.7484 corr=0.401535 error=0.182941 0.000639655
ijk=(0 0 4) x1=29.443 x2=29.3564 x3=17.0234 corr=0.298253 error=3.73962 0.0130756
ijk=(0 1 0) x1=29.443 x2=17.0427 x3=31.0352 corr=0.131718 error=4.45946 0.0155925
ijk=(0 1 1) x1=29.443 x2=17.0427 x3=30.3301 corr=0.163135 error=4.24559 0.0148447
ijk=(0 1 2) x1=29.443 x2=17.0427 x3=28.9253 corr=0.240393 error=3.81943 0.0133547
ijk=(0 1 3) x1=29.443 x2=17.0427 x3=28.7484 corr=0.251553 error=3.76577 0.013167
ijk=(0 1 4) x1=29.443 x2=17.0427 x3=17.0234 corr=0.997818 error=0.209094 0.000731099
ijk=(1 0 0) x1=16.985 x2=29.3564 x3=31.0352 corr=0.996218 error=0.294805 0.00103079
ijk=(1 0 1) x1=16.985 x2=29.3564 x3=30.3301 corr=0.999697 error=0.0809304 0.000282974
ijk=(1 0 2) x1=16.985 x2=29.3564 x3=28.9253 corr=0.993806 error=0.345223 0.00120707
ijk=(1 0 3) x1=16.985 x2=29.3564 x3=28.7484 corr=0.991618 error=0.398886 0.00139471
ijk=(1 0 4) x1=16.985 x2=29.3564 x3=17.0234 corr=0.211781 error=3.95556 0.0138306
ijk=(1 1 0) x1=16.985 x2=17.0427 x3=31.0352 corr=0.295974 error=4.24352 0.0148375
ijk=(1 1 1) x1=16.985 x2=17.0427 x3=30.3301 corr=0.296146 error=4.02964 0.0140897
ijk=(1 1 2) x1=16.985 x2=17.0427 x3=28.9253 corr=0.29655 error=3.60349 0.0125996
ijk=(1 1 3) x1=16.985 x2=17.0427 x3=28.7484 corr=0.296608 error=3.54983 0.012412
ijk=(1 1 4) x1=16.985 x2=17.0427 x3=17.0234 corr=0.860665 error=0.00685069 2.39534e-05
ijk=(1 1 4) x1=16.985 x2=17.0427 x3=17.0234 a=17.034 b=5.95327e-05<-277
ijk=(0 0 0) x1=29.443 x2=29.3564 x3=31.0352 corr=0.251963 error=0.51075 0.00178584
ijk=(0 0 1) x1=29.443 x2=29.3564 x3=30.3301 corr=0.222541 error=0.296875 0.00103803
ijk=(0 0 2) x1=29.443 x2=29.3564 x3=28.9253 corr=0.443144 error=0.129278 0.00045202
ijk=(0 0 3) x1=29.443 x2=29.3564 x3=28.7484 corr=0.401535 error=0.182941 0.000639655
ijk=(0 0 2) x1=29.443 x2=29.3564 x3=28.9253 a=29.1262 b=-0.000403236<-Data saved to temporary root file.
10/19/10
M.Byszewski - batreco
12