Warning, /analysis/ParticleID/RICHAnalysis_macros/macros_analyze/readme.txt is written in an unsupported language. File is not indexed.
0001 HOW TO USE THE CODE IN ROOT
0002
0003 .L bnl_rich_analyzer_pid.cpp++
0004 eic_bnl_rich *a = new eic_bnl_rich()
0005
0006 How to obtain the relevant results
0007 ___________________________________________________________
0008
0009 a->acq("../eval_RICH_e_10k.root",10000,1,0)
0010
0011
0012 where the arguments are:
0013
0014 1: .root file name
0015 2: number of events to be processed <= events in the file
0016 3: 0 (pizel size effects off) , 1 (pizel size effects on 3mm)
0017 4: particle type (0 e, 1 pi , 2 K , 3 p)
0018
0019 The output on the sceen is like this
0020
0021 0 0 0 0 0
0022 0.0343412 0.00136802 3 1 0
0023 0.0334033 0.00188749 30 2 0
0024 0.0333956 0.00157126 28 3 0
0025 0.0334581 0.00146258 28 4 0
0026 0.0334142 0.00141108 29 5 0
0027 0.0334099 0.00137278 29 6 0
0028 0.0334005 0.00133371 29 7 0
0029 0.0334324 0.00134786 30 8 0
0030 0.0334299 0.00132511 28 9 0
0031 0.0334531 0.00131175 29 10 0
0032 0.0334187 0.0013297 29 11 0
0033 0.0334165 0.00132093 29 12 0
0034 0.0334248 0.00129988 28 13 0
0035 0.0334396 0.00128306 30 14 0
0036 0.0334506 0.00130121 29 15 0
0037 0.0334684 0.00128705 29 16 0
0038 0.0334535 0.00131191 29 17 0
0039 0.0334197 0.00131452 29 18 0
0040 0.0334226 0.00130962 28 19 0
0041 0.033431 0.00128808 29 20 0
0042 0.0334568 0.00128365 29 21 0
0043 0.0334206 0.00128394 29 22 0
0044 0.0334155 0.00127789 28 23 0
0045 0.0334248 0.0012976 28 24 0
0046 0.0334277 0.0012828 27 25 0
0047 0.0334087 0.0012881 30 26 0
0048 0.0334435 0.0012975 29 27 0
0049 0.0334481 0.00130363 29 28 0
0050 0.0334212 0.00128097 29 29 0
0051 0.0334251 0.00127178 27 30 0
0052 0.0334329 0.00128094 30 31 0
0053 0.0334164 0.00128659 29 32 0
0054 0.0334172 0.00126684 28 33 0
0055 0.0334375 0.00128586 29 34 0
0056 0.0333946 0.00128502 27 35 0
0057 0.0334234 0.00128111 30 36 0
0058 0.0333789 0.00132507 29 37 0
0059 0.0334289 0.00126741 28 38 0
0060 0.033438 0.00129313 28 39 0
0061 0.0334268 0.00129156 28 40 0
0062 0.0334032 0.00131998 28 41 0
0063 0.0334517 0.00129523 29 42 0
0064 0.0334286 0.00127945 29 43 0
0065 0.033442 0.00127814 29 44 0
0066 0.0334369 0.00128676 29 45 0
0067 0.0334225 0.00131866 28 46 0
0068 0.0334726 0.00130807 28 47 0
0069 0.0334476 0.00131119 29 48 0
0070 0.0334379 0.00128868 28 49 0
0071 0.0334272 0.00127614 28 50 0
0072 0.0334564 0.00127255 28 51 0
0073 0.0334011 0.00130214 28 52 0
0074 0.0334446 0.001267 28 53 0
0075 0.0334112 0.00128135 29 54 0
0076 0.0334346 0.001309 27 55 0
0077 0.0334288 0.00129947 29 56 0
0078 0.0334127 0.0012866 29 57 0
0079 0.0334443 0.00127512 30 58 0
0080 0.0334045 0.00128447 29 59 0
0081
0082 Three line in the code can be uncommented to print in the out file
0083
0084
0085 How to reproduce the plots
0086 __________________________________________________________________
0087
0088 Run the simulation for all the particles:
0089
0090 a->acq("../eval_RICH_e_10k.root",10000,1,0)
0091 a->acq("../eval_RICH_pi_10k.root",10000,1,1)
0092 a->acq("../eval_RICH_k_10k.root",10000,1,2)
0093 a->acq("../eval_RICH_p_10k.root",10000,1,3)
0094
0095 and save the outputs in a single file, i.e. pid.txt
0096
0097 Then use the raed funtion:
0098
0099 a->read() \\pid.txt called inside
0100
0101 3 plots in output
0102
0103
0104
0105
0106
