How to decode MWPC data » History » Version 1
Jixie Zhang, 01/21/2021 04:07 AM
| 1 | 1 | Jixie Zhang | h1. How to decode MWPC data |
|---|---|---|---|
| 2 | |||
| 3 | The decoder of MWPC can finally create root tree. |
||
| 4 | |||
| 5 | The decoder is located at /home/a-compton/MWPC/analysis/track/anatka.C |
||
| 6 | |||
| 7 | MWPC calibration is stored at /home/a-compton/MWPC/analysis/track/anatk.inp |
||
| 8 | It includes the reference positions of each chamber, also the pedestal of each TDC module. |
||
| 9 | Most importantly, it tells the time cut value which will be used to determine if a hit is out of time. |
||
| 10 | |||
| 11 | Here is how is looks: |
||
| 12 | ---------------------------------------------------------------------------------------------------------------------------- |
||
| 13 | [a-compton@soliddaq track]> cat anatk.inp |
||
| 14 | /home/nfs/ftbf_user/experiments/EDIT2018/MWPC/data/exp_20180312155107.dat ifile |
||
| 15 | 1 |
||
| 16 | 5000000 |
||
| 17 | 198.8 475.0 758.8 1534.5 zz |
||
| 18 | 6.362 6.498 6.45 6.48 xi |
||
| 19 | 6.44 6.49 6.41 6.42 yi |
||
| 20 | 0.45 |
||
| 21 | 59 58 59 60 50 51 51 52 |
||
| 22 | 50 50 53 54 50 50 49 48 |
||
| 23 | 18 | itwin (time window) |
||
| 24 | |||
| 25 | ======================================================================== |
||
| 26 | line 1 file name of data to process |
||
| 27 | line 2 file flag. if 1 use the file selection, if 2 use the name from liine 1 |
||
| 28 | line 3 event limit |
||
| 29 | line 4 Z locations of chambers. 0 is upstream face of hall. units cm |
||
| 30 | line 5 x0 location of first wire in x |
||
| 31 | line 6 y0 |
||
| 32 | line 7 Allowable window for finding possible hits on a track |
||
| 33 | line 8 Center of time distributatioin TCDs 0 - 7 |
||
| 34 | line 8 8 - 15 |
||
| 35 | line 9 Require time in the tdcs be within +/- this value |
||
| 36 | ---------------------------------------------------------------------------------------------------------------------------- |
||
| 37 | |||
| 38 | |||
| 39 | How to run the decoder: |
||
| 40 | step 1): synchronize the raw data file /home/a-compton/MWPC/cpdata.csh |
||
| 41 | step 2) run decoder: /home/a-compton/MWPC/analysis/track/wctree.csh |
||
| 42 | |||
| 43 | |||
| 44 | A) It will ask you which file you want to decode, you make your choice here, usually |
||
| 45 | the last file is the latest run. |
||
| 46 | B) After you select which file to decode, it will read the first line to tell you what is the |
||
| 47 | possible run number. Some times this value is wrong. You can type in the correct run |
||
| 48 | number if you are sure, otherwise just copy the run number. |
||
| 49 | C) What is the output? |
||
| 50 | root tree will be stored at /home/a-compton/MWPC/wcnt_xxx.root, where xxx is the run number. |
||
| 51 | beam position figure will be store at /home/a-compton/MWPC/log/BeamPosition_xxx.png |
||
| 52 | |||
| 53 | ======================================================= |
||
| 54 | The root file includes a lot of monitoring histogram, also the tree named 'M': |
||
| 55 | Here is the tree branches: |
||
| 56 | |||
| 57 | spill -->spill number |
||
| 58 | evn -->event number |
||
| 59 | etime -->event time relative to spill start |
||
| 60 | nhit_raw -->number of raw hits |
||
| 61 | wire_raw -->wire id of raw hits |
||
| 62 | tdc_raw --> tdc of raw hits |
||
| 63 | nhit_tcut -->number of hits after time cut |
||
| 64 | wire_tcut -->wire id of hits after time cut |
||
| 65 | tdc_tcut --> tdc values after subtracted the pedestal and pass time cut |
||
| 66 | |||
| 67 | nx0 -->number of hits in chamber0 X plane |
||
| 68 | x0 -->X positions of chamber0 X plane |
||
| 69 | ny0 -->number of hits in chamber0 Y plane |
||
| 70 | y0 -->Y positions of chamber0 Y plane |
||
| 71 | |||
| 72 | nx1 -->number of hits in chamber1 X plane |
||
| 73 | x1 -->X positions of chamber1 X plane |
||
| 74 | ny1 -->number of hits in chamber1 Y plane |
||
| 75 | y1 -->Y positions of chamber1 Y plane |
||
| 76 | |||
| 77 | nx2 -->number of hits in chamber2 X plane |
||
| 78 | x2 -->X positions of chamber2 X plane |
||
| 79 | ny2 -->number of hits in chamber2 Y plane |
||
| 80 | y2 -->Y positions of chamber2 Y plane |
||
| 81 | |||
| 82 | nx3 -->number of hits in chamber3 X plane |
||
| 83 | x3 -->X positions of chamber3 X plane |
||
| 84 | ny3 -->number of hits in chamber3 Y plane |
||
| 85 | y3 -->Y positions of chamber3 Y plane |
||
| 86 | |||
| 87 | |||
| 88 | ======================================================= |
||
| 89 | The description of monitoring histogram is here: |
||
| 90 | OBJ: TH1F h100 All Wire Hits : 0 at: 0x1ed4740 |
||
| 91 | OBJ: TH1F h200 M TDC 0 : 0 at: 0x1f7e800 |
||
| 92 | KEY: TH1F h301;1 X Dev 2 1234 |
||
| 93 | KEY: TH1F h302;1 X Dev 3 1234 |
||
| 94 | KEY: TH1F h311;1 Y Dev 2 1234 |
||
| 95 | KEY: TH1F h312;1 Y Dev 3 1234 |
||
| 96 | KEY: TH1F h305;1 X Dev 2 123 |
||
| 97 | KEY: TH1F h306;1 X Dev 2 124 |
||
| 98 | KEY: TH1F h307;1 X Dev 3 134 |
||
| 99 | KEY: TH1F h308;1 X Dev 3 234 |
||
| 100 | KEY: TH1F h315;1 Y Dev 2 123 |
||
| 101 | KEY: TH1F h316;1 Y Dev 2 124 |
||
| 102 | KEY: TH1F h317;1 Y Dev 3 134 |
||
| 103 | KEY: TH1F h318;1 Y Dev 3 234 |
||
| 104 | KEY: TH1F h380;1 Raw Hit/Event |
||
| 105 | KEY: TH1F h381;1 Hits/Event after Tcut |
||
| 106 | KEY: TH1F h390;1 Times after Tcut |
||
| 107 | KEY: TH1F h391;1 Times before Tcut |
||
| 108 | KEY: TH1F h392;1 Offset Times before Tcut |
||
| 109 | KEY: TH1F h100;1 All Wire Hits |
||
| 110 | KEY: TH1F h120;1 Event Time |
||
| 111 | KEY: TH1F h121;1 Time bet. Ev |
||
| 112 | KEY: TH1F h130;1 Events per Spill |
||
| 113 | KEY: TH1F h150;1 h150 Wire Chamber 1 X |
||
| 114 | KEY: TH1F h151;1 h151 Wire Chamber 1 Y |
||
| 115 | KEY: TH1F h152;1 h152 Wire Chamber 2 X |
||
| 116 | KEY: TH1F h153;1 h153 Wire Chamber 2 Y |
||
| 117 | KEY: TH1F h154;1 h154 Wire Chamber 3 X |
||
| 118 | KEY: TH1F h155;1 h155 Wire Chamber 3 Y |
||
| 119 | KEY: TH1F h156;1 h156 Wire Chamber 4 X |
||
| 120 | KEY: TH1F h157;1 h157 Wire Chamber 4 Y |
||
| 121 | KEY: TH1F h160;1 Hits 1X |
||
| 122 | KEY: TH1F h161;1 Hits 1Y |
||
| 123 | KEY: TH1F h162;1 Hits 2X |
||
| 124 | KEY: TH1F h163;1 Hits 2Y |
||
| 125 | KEY: TH1F h164;1 Hits 3X |
||
| 126 | KEY: TH1F h165;1 Hits 3Y |
||
| 127 | KEY: TH1F h166;1 Hits 4X |
||
| 128 | KEY: TH1F h167;1 Hits 4Y |
||
| 129 | KEY: TH1F h110;1 All Time Hits |
||
| 130 | KEY: TH2F h181;1 Chamber 1 |
||
| 131 | KEY: TH2F h182;1 Chamber 2 |
||
| 132 | KEY: TH2F h183;1 Chamber 3 |
||
| 133 | KEY: TH2F h184;1 Chamber 4 |
||
| 134 | KEY: TH1F h200;1 M TDC 0 |
||
| 135 | KEY: TH1F h201;1 M TDC 1 |
||
| 136 | KEY: TH1F h202;1 M TDC 2 |
||
| 137 | KEY: TH1F h203;1 M TDC 3 |
||
| 138 | KEY: TH1F h204;1 M TDC 4 |
||
| 139 | KEY: TH1F h205;1 M TDC 5 |
||
| 140 | KEY: TH1F h206;1 M TDC 6 |
||
| 141 | KEY: TH1F h207;1 M TDC 7 |
||
| 142 | KEY: TH1F h208;1 M TDC 8 |
||
| 143 | KEY: TH1F h209;1 M TDC 9 |
||
| 144 | KEY: TH1F h210;1 M TDC 10 |
||
| 145 | KEY: TH1F h211;1 M TDC 11 |
||
| 146 | KEY: TH1F h212;1 M TDC 12 |
||
| 147 | KEY: TH1F h213;1 M TDC 13 |
||
| 148 | KEY: TH1F h214;1 M TDC 14 |
||
| 149 | KEY: TH1F h215;1 M TDC 15 |
||
| 150 | KEY: TH2F h320;1 X(6) Y(5) Tracks |
||
| 151 | KEY: TH1F h321;1 X Tracks |
||
| 152 | KEY: TH1F h322;1 Y Tracks |
||
| 153 | KEY: TH1F h341;1 X Track Slope |
||
| 154 | KEY: TH1F h342;1 X Track Intercept |
||
| 155 | KEY: TH1F h343;1 Y Track Slope |
||
| 156 | KEY: TH1F h344;1 Y Track Intercept |
||
| 157 | KEY: TNtuple trk;1 StLine |
||
| 158 | ======================================================= |
||
| 159 | |||
| 160 |