Changeset 582
- Timestamp:
- Oct 21, 2011, 9:09:14 AM (13 years ago)
- Location:
- BAORadio/AmasNancay/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
BAORadio/AmasNancay/trunk/etude_mergeAna.pic
r577 r582 1 2 3 set source $1 4 set ncycles $2 5 6 7 #set fcalib $3 8 #set f1420low $4 9 #set f1420high $5 10 11 set fcalib 1410 12 set f1420low 1420.2 13 set f1420high 1420.6 14 15 16 17 nbins1420 = floor((${f1420high}-${f1420low})/250*8192) 18 19 echo $nbins1420 20 21 22 23 clearscript diffOnOff 24 clearscript evolAt1420Freq 25 clearscript evolAt1420SideFreq 26 clearscript IntAt1420 27 clearscript sigmaRaw1420Side 28 ############################################################################ 29 1 30 defscript evolAtCalibFreq 2 31 … … 56 85 57 86 58 #Follow up of the intensity arround the (1420.4 +/- 2)MHz frequency in different calibration conditions87 #Follow up of the intensity arround the [1420,1421]MHz frequency in different calibration conditions 59 88 # NO calibration 60 89 # Calibration coeff. mean per Run … … 68 97 settitle "Raw ON-OFF ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt 69 98 n/pl onoffevol.onoffRaw11420%cycle ! ! "red cpts notit nsta" 70 setaxelabels "cycle" "I[ 1420.4MHz +/- 2MHz]Raw (a.u)" $axedefatt99 setaxelabels "cycle" "I[${f1420low},${f1420high}]MHz Raw (a.u)" $axedefatt 71 100 72 101 … … 75 104 settitle "ON-OFF Calib Run ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 76 105 n/pl onoffevol.onoffRun11420%cycle ! ! "red cpts notit nsta" 77 setaxelabels "cycle" "I[ 1420.4MHz +/- 2MHz]Run (a.u)" $axedefatt106 setaxelabels "cycle" "I[${f1420low},${f1420high}]MHz Run (a.u)" $axedefatt 78 107 79 108 … … 82 111 settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 83 112 n/pl onoffevol.onoffCycle11420%cycle ! ! "red cpts notit nsta" 84 setaxelabels "cycle" "I[ 1420.4MHz +/- 2MHz]Cycle (a.u)" $axedefatt85 86 87 newh1d hSigmaRaw0 -0.0 2 0.02 5088 newh1d hSigmaRaw1 -0.0 2 0.02 5089 newh1d hSigmaRun0 -0.0 2 0.02 5090 newh1d hSigmaRun1 -0.0 2 0.02 5091 newh1d hSigmaCycl0 -0.0 2 0.02 5092 newh1d hSigmaCycl1 -0.0 2 0.02 50113 setaxelabels "cycle" "I[${f1420low},${f1420high}]MHz Cycle (a.u)" $axedefatt 114 115 116 newh1d hSigmaRaw0 -0.01 0.01 25 117 newh1d hSigmaRaw1 -0.01 0.01 25 118 newh1d hSigmaRun0 -0.01 0.01 25 119 newh1d hSigmaRun1 -0.01 0.01 25 120 newh1d hSigmaCycl0 -0.01 0.01 25 121 newh1d hSigmaCycl1 -0.01 0.01 25 93 122 94 123 graphicatt "" … … 98 127 n/proj hSigmaRaw1 onoffevol.onoffRaw11420 ! ! ! "same red cpts notit" 99 128 settitle "Raw ON-OFF ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt 100 setaxelabels "(a.u)" "I[ 1420.4MHz +/- 2MHz]Raw" $axedefatt129 setaxelabels "(a.u)" "I[${f1420low},${f1420high}]MHz Raw" $axedefatt 101 130 102 131 newwin 1 1 … … 104 133 n/proj hSigmaRun1 onoffevol.onoffRun11420 ! ! ! "same red cpts notit" 105 134 settitle "ON-OFF Calib Run ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 106 setaxelabels "(a.u)" "I[ 1420.4MHz +/- 2MHz]Run" $axedefatt135 setaxelabels "(a.u)" "I[${f1420low},${f1420high}]MHz Run" $axedefatt 107 136 108 137 … … 111 140 n/proj hSigmaCycl1 onoffevol.onoffCycle11420 ! ! ! "same red cpts notit" 112 141 settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 113 setaxelabels "(a.u)" "I[1420.4MHz +/- 2MHz] Cycle" $axedefatt 142 setaxelabels "(a.u)" "I[${f1420low},${f1420high}]MHz Cycle" $axedefatt 143 144 145 endscript 146 ################################################## 147 defscript sigmaRaw1420Side 148 149 150 newh1d hSigRaw0All -0.01 0.01 25 151 newh1d hSigRaw1All -0.01 0.01 25 152 153 n/proj hSigRaw0All onoffevol.onoffRaw01420side 154 n/proj hSigRaw1All onoffevol.onoffRaw11420side 155 156 echo "Resultats sur ${ncycles} cycles (Ch0, Ch1): mean (${hSigRaw0All.mean}, ${hSigRaw1All.mean}) sigma (${hSigRaw0All.sigma}, ${hSigRaw1All.sigma})" 157 158 newwin 1 1 159 disp hSigRaw0All "blue cpts notit" 160 disp hSigRaw1All "red same cpts notit" 161 settitle "Mean and Sigma for 1 cycle" 162 163 #Ts les 10cycles 164 n10cycles = floor(${ncycles}/10) 165 x10min = ${hSigRaw0All.xmin} 166 x10max = ${hSigRaw0All.xmax} 167 n10bin = ${hSigRaw0All.nbin} 168 169 170 for i 0:${n10cycles} 171 ic = ${i}*10+1 172 jc = ${ic}+9 173 newh1d hSigRaw010c${i} ${x10min} ${x10max} ${n10bin} 174 newh1d hSigRaw110c${i} ${x10min} ${x10max} ${n10bin} 175 176 n/proj hSigRaw010c${i} onoffevol.onoffRaw01420side (cycle>=${ic})&&(cycle<=${jc}) 177 cp hSigRaw010c${i} h0 178 mSig010c[i] = ${h0.mean} 179 180 n/proj hSigRaw110c${i} onoffevol.onoffRaw11420side (cycle>=${ic})&&(cycle<=${jc}) 181 cp hSigRaw110c${i} h1 182 mSig110c[i] = ${h1.mean} 183 end 184 185 186 line2vec hSigRaw010c $mSig010c 187 line2vec hSigRaw110c $mSig110c 188 189 190 191 192 newwin 1 1 193 disp hSigRaw010c "blue cpts notit" 194 disp hSigRaw110c "red same cpts notit" 195 settitle "Mean and Sigma for 10 cycles" 196 197 198 endscript 199 ########################################### 200 defscript evolAt1420SideFreq 201 202 203 #Follow up of the intensity arround the [${f1420low},${f1420high}]MHz frequency in different calibration conditions 204 # NO calibration 205 # Calibration coeff. mean per Run 206 # Calibration coeff. per cycles 207 208 209 graphicatt "xylimits=0,500,-0.01,0.01" 210 211 newwin 1 2 212 n/pl onoffevol.onoffRaw01420side%cycle ! ! "blue cpts notit nsta" 213 settitle "Raw ON-OFF ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt 214 n/pl onoffevol.onoffRaw11420side%cycle ! ! "red cpts notit nsta" 215 setaxelabels "cycle" "I[1418,1419]U[1422,1423]MHz Raw (a.u)" $axedefatt 216 217 218 newwin 1 2 219 n/pl onoffevol.onoffRun01420side%cycle ! ! "blue cpts notit nsta" 220 settitle "ON-OFF Calib Run ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 221 n/pl onoffevol.onoffRun11420side%cycle ! ! "red cpts notit nsta" 222 setaxelabels "cycle" "I[1418,1419]U[1422,1423]MHz Run (a.u)" $axedefatt 223 224 225 newwin 1 2 226 n/pl onoffevol.onoffCycle01420side%cycle ! ! "blue cpts notit nsta" 227 settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 228 n/pl onoffevol.onoffCycle11420side%cycle ! ! "red cpts notit nsta" 229 setaxelabels "cycle" "I[1418,1419]U[1422,1423]MHz Cycle (a.u)" $axedefatt 230 231 232 newh1d hSigmaSideRaw0 -0.01 0.01 25 233 newh1d hSigmaSideRaw1 -0.01 0.01 25 234 newh1d hSigmaSideRun0 -0.01 0.01 25 235 newh1d hSigmaSideRun1 -0.01 0.01 25 236 newh1d hSigmaSideCycl0 -0.01 0.01 25 237 newh1d hSigmaSideCycl1 -0.01 0.01 25 238 239 graphicatt "" 240 241 newwin 1 1 242 n/proj hSigmaSideRaw0 onoffevol.onoffRaw01420side ! ! ! "blue cpts notit" 243 n/proj hSigmaSideRaw1 onoffevol.onoffRaw11420side ! ! ! "same red cpts notit" 244 settitle "Raw ON-OFF ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt 245 setaxelabels "(a.u)" "I[1418,1419]U[1422,1423]MHz Raw" $axedefatt 246 247 newwin 1 1 248 n/proj hSigmaSideRun0 onoffevol.onoffRun01420side ! ! ! "blue cpts notit" 249 n/proj hSigmaSideRun1 onoffevol.onoffRun11420side ! ! ! "same red cpts notit" 250 settitle "ON-OFF Calib Run ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 251 setaxelabels "(a.u)" "I[1418,1419]U[1422,1423]MHz Run" $axedefatt 252 253 254 newwin 1 1 255 n/proj hSigmaSideCycl0 onoffevol.onoffCycle01420side ! ! ! "blue cpts notit" 256 n/proj hSigmaSideCycl1 onoffevol.onoffCycle11420side ! ! ! "same red cpts notit" 257 settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 258 setaxelabels "(a.u)" "I[1418,1419]U[1422,1423]MHz Cycle" $axedefatt 259 260 261 endscript 262 263 ########################################### 264 defscript IntAt1420 265 266 267 #Follow up of the intensity arround the [${f1420low},${f1420high}]MHz 268 #frequency in different calibration conditions with baseline removed 269 # NO calibration 270 # Calibration coeff. mean per Run 271 # Calibration coeff. per cycles 272 273 274 graphicatt "xylimits=0,500,-0.1,0.1" 275 276 newwin 1 2 277 n/pl onoffevol.${nbins1420}*(onoffRaw01420-onoffRaw01420side)%cycle ! ! "blue cpts notit nsta" 278 settitle "Raw ON-OFF ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt 279 n/pl onoffevol.${nbins1420}*(onoffRaw11420-onoffRaw11420side)%cycle ! ! "red cpts notit nsta" 280 setaxelabels "cycle" "I 1420MHz Raw (a.u)" $axedefatt 281 282 283 newwin 1 2 284 n/pl onoffevol.${nbins1420}*(onoffRun01420-onoffRun01420side)%cycle ! ! "blue cpts notit nsta" 285 settitle "ON-OFF Calib Run ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 286 n/pl onoffevol.${nbins1420}*(onoffRun11420-onoffRun11420side)%cycle ! ! "red cpts notit nsta" 287 setaxelabels "cycle" "I 1420MHz Run (a.u)" $axedefatt 288 289 290 newwin 1 2 291 n/pl onoffevol.${nbins1420}*(onoffCycle01420-onoffCycle01420side)%cycle ! ! "blue cpts notit nsta" 292 settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 293 n/pl onoffevol.${nbins1420}*(onoffCycle11420-onoffCycle11420side)%cycle ! ! "red cpts notit nsta" 294 setaxelabels "cycle" "I 1420MHz Cycle (a.u)" $axedefatt 295 296 297 newh1d h1420IntRaw0 -0.1 0.1 25 298 newh1d h1420IntRaw1 -0.1 0.1 25 299 newh1d h1420IntRun0 -0.1 0.1 25 300 newh1d h1420IntRun1 -0.1 0.1 25 301 newh1d h1420IntCycl0 -0.1 0.1 25 302 newh1d h1420IntCycl1 -0.1 0.1 25 303 304 graphicatt "" 305 306 newwin 1 1 307 n/proj h1420IntRaw0 onoffevol.${nbins1420}*(onoffRaw01420-onoffRaw01420side) ! ! ! "blue cpts notit" 308 n/proj h1420IntRaw1 onoffevol.${nbins1420}*(onoffRaw11420-onoffRaw11420side) ! ! ! "same red cpts notit" 309 settitle "Raw ON-OFF ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt 310 setaxelabels "(a.u)" "I 1420MHz Raw" $axedefatt 311 312 newwin 1 1 313 n/proj h1420IntRun0 onoffevol.${nbins1420}*(onoffRun01420-onoffRun01420side) ! ! ! "blue cpts notit" 314 n/proj h1420IntRun1 onoffevol.${nbins1420}*(onoffRun11420-onoffRun11420side) ! ! ! "same red cpts notit" 315 settitle "ON-OFF Calib Run ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 316 setaxelabels "(a.u)" "I 1420MHz Run" $axedefatt 317 318 319 320 newwin 1 1 321 n/proj h1420IntCycl0 onoffevol.${nbins1420}*(onoffCycle01420-onoffCycle01420side) ! ! ! "blue cpts notit" 322 n/proj h1420IntCycl1 onoffevol.${nbins1420}*(onoffCycle11420-onoffCycle11420side) ! ! ! "same red cpts notit" 323 settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt 324 setaxelabels "(a.u)" "I 1420MHz Cycle" $axedefatt 325 114 326 115 327 … … 164 376 #use minuscule for source label 165 377 166 set source $1 167 set ncycles $2 168 set fcalib $3 378 379 169 380 170 381 openppf onoffsurvey_${source}-${ncycles}Cycles.ppf … … 181 392 182 393 183 184 394 #diffOnOff 185 186 187 evolAt1420Freq 395 #evolAt1420Freq 396 #evolAt1420SideFreq 397 #IntAt1420 -
BAORadio/AmasNancay/trunk/mergeAnaFiles.cc
r577 r582 251 251 //Compute the mean of Diff ON-OFF BAO-calibrated spectra and also the mean/sigma of rebinned spectra 252 252 // 253 void meanCalibBAODiffOnOffCycles() throw(string) {253 void 1meanCalibBAODiffOnOffCycles() throw(string) { 254 254 255 255 list<string> listOfFiles; … … 286 286 sa_size_t chCalibHigh = freqToChan(para.rcalibFreq_ + (para.rcalibBandFreq_*0.5)); 287 287 //Lower and Higher freq. just arround 1420.4MHz Freq. bin to perform mean follow up 288 sa_size_t ch1420Low = freqToChan(1420 );289 sa_size_t ch1420High = freqToChan(142 1);288 sa_size_t ch1420Low = freqToChan(1420.4-0.2); 289 sa_size_t ch1420High = freqToChan(1420.4+0.2); 290 290 291 291 //Lower and Higher freq. on the sides of 1420.4Mhz Freq. bin to perform mean follow up … … 296 296 297 297 298 if (para.debuglev_> 90){298 if (para.debuglev_>0){ 299 299 cout << "freq. band for follow up [" << chCalibLow << ", "<< chCalibHigh << "]" << endl; 300 300 cout << "freq. band for follow up [" << ch1420Low << ", "<< ch1420High << "]" << endl; 301 cout << "freq. band for follow up [" << ch1420aLow << ", "<< ch1420aHigh << "]" << endl; 302 cout << "freq. band for follow up [" << ch1420bLow << ", "<< ch1420bHigh << "]" << endl; 301 303 } 302 304
Note: See TracChangeset
for help on using the changeset viewer.