[d2n-analysis-talk] Problem with Reconstructed Energy in One-Pass H2 Runs
Diana Parno
dparno at cmu.edu
Thu Feb 10 17:25:06 EST 2011
I've found the source of the problem that Matt observed with the
reconstructed energy in the one-pass H2 runs. It looks like there was
a mistake in the combination of StartType.pl and DBs that I had
checked out from the
Here's the location of the copy of StartType.pl used for the farm
replays:
/work/halla/e06014/disk1/diana/replay_32/20101118/d2n_analysis/d2n/DB/
StartType.pl
And here's the current d2n copy, which is what Matt used in his recent
test replay:
/usr/local/d2n_analysis/d2n/DB/StartType.pl
I'm looking at the subroutine adjBBps_sh. In the Jan-3 update to the
changelog, Matt wrote
> - Disabled the shower start type (for now)
>
> - Copied the db_BB.ts.ps.dat and db_BB.ts.sh.dat to DB DIR
Indeed, pretty much all of that subroutine is commented out. But in
the farm version, that subroutine contains a couple of commands:
> if ($runno>=1231 && $runno<=1294){
> print " Adjusting the BB ps/sh Energy coefficients. \n";
> dolink("$dbpath/db_BB.ts.ps.dat.1", "db_BB.ts.ps.dat");
> dolink("$dbpath/db_BB.ts.sh.dat.1", "db_BB.ts.sh.dat");
> return;
> }
Here's what happened. I checked out a CVS copy of the databases,
StartType.pl, etc, but that predated Matt's calculation of the shower/
preshower coefficients. I got the updated DB files from Matt and stuck
them in the main DB directory. I didn't realize that, for specific run
numbers, these lines in StartType.pl would use old DB files (from the
CVS copy) rather than the new ones. I'm pretty sure that the old
copies lack the coefficients entirely. I've attached the output of
some diff commands at the end of the email for general interest; I'm
not familiar with these DB files, but those look like energy
calibration coefficients to me.
So, this should be pretty easily fixable; all one has to do is disable
this section of StartType.pl. Then the correct DB file will be pulled,
and hooray. This is something that we can either do now (and re-replay
those few runs), or in the next Full Replay Round, whenever that
happens.
The good news is that this (the one-pass H2 dataset) is literally the
only portion of our farm-replayed runs that should be affected. I was
initially worried about the "if ($runno>=2087 && $runno<=2101)"
condition, but we don't have any production runs in that range. There
are a few five-pass runs that may be in the production range, but we
haven't done a farm replay of the five-pass runs yet, so nothing has
to be redone there.
Best,
Diana
[dseymour at d2n:~]$ diff /work/halla/e06014/disk1/diana/
replay_32/20101118/d2n_analysis/d2n/DB/20090210/db_BB.ts.ps.dat.1 /usr/
local/d2n_analysis/d2n/DB/db_BB.ts.ps.dat
16c16
< 0.0 0.0 0.77 - Meters
---
> -0.015 0.0 0.97 - Meters
31,36c31,36
< 336.0 329.0 313.0 321.0 418.0 419.0 429.0 434.0 367.0
< 340.0 328.0 330.0 484.0 416.0 350.0 409.0 371.0 303.0
< 244.0 300.0 399.0 443.0 415.0 471.0 395.0 376.0 299.0
< 337.0 383.0 445.0 323.0 451.0 374.0 621.0 418.0 540.0
< 432.0 482.0 307.0 533.0 559.0 482.0 414.0 498.0 405.0
< 440.0 298.0 669.0 653.0 567.0 471.0 550.0 319.0 385.0
---
> 337.0 329.0 313.0 321.0 419.0 419.0 429.0 440.0 366.0
> 340.0 328.0 330.0 484.0 417.0 351.0 409.0 371.0 303.0
> 244.0 300.0 399.0 443.0 415.0 471.0 395.0 376.0 299.0
> 337.0 383.0 455.0 323.0 451.0 377.0 621.0 418.0 540.0
> 432.0 482.0 307.0 533.0 559.0 482.0 414.0 449.0 405.0
> 441.0 298.0 669.0 653.0 567.0 471.0 550.0 319.0 385.0
38,42c38,46
< 1 1 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1 1 1
---
> 0.635367 0.697452 0.674238 0.668869 0.633355 0.632939
> 0.662406 0.710782 0.807697 0.879553 0.834624 0.832439
> 0.771791 0.772114 0.761963 0.735849 0.670249 0.657658
> 0.694703 0.665494 0.456888 0.599011 0.653604 0.725153
> 0.694023 0.586299 0.399231 0.650773 0.633457 0.591386
> 0.620168 0.550137 0.653063 0.688886 0.689436 0.863534
> 0.922742 0.865032 0.838056 0.83525 0.718808 0.724633
> 0.714448 0.640682 0.605775 0.643843 0.67314 0.635458
> 0.613033 0.646124 0.706132 0.724086 0.589599 0.311727
[dseymour at d2n:~]$ diff /work/halla/e06014/disk1/diana/
replay_32/20101118/d2n_analysis/d2n/DB/20090210/db_BB.ts.sh.dat.1 /usr/
local/d2n_analysis/d2n/DB/db_BB.ts.sh.dat
77,79c77,79
< 300.0 382.0 336.0 448.0 477.5 393.5 322.5
< 446.5 340.5 422.5 406.0 532.5 358.2 436.0
< 328.0 526.0 359.5 413.0 471.0 516.0 364.0
---
> 300.0 382.0 336.0 448.0 478.0 394.0 323.0
> 447.0 340.0 423.0 406.0 533.0 358.0 436.0
> 328.0 526.0 359.0 413.0 471.0 516.0 364.0
83,84c83,84
< 325.0 400.0 570.0 366.0 536.0 526.0 466.0
< 355.0 327.0 296.0 395.0 328.0 298.0 293.0
---
> 325.0 400.0 570.0 367.0 536.0 526.0 466.0
> 355.0 327.0 297.0 395.0 328.0 298.0 293.0
88,90c88,90
< 456.0 308.0 316.0 271.0 280.0 283.0 282.0
< 269.0 257.0 373.0 302.0 289.0 278.0 353.0
< 319.0 318.0 339.0 401.0 290.0 304.0 318.0
---
> 456.0 309.0 316.0 271.0 280.0 283.0 282.0
> 269.0 257.0 373.0 302.0 289.0 279.0 353.0
> 319.0 318.0 339.0 401.0 291.0 305.0 318.0
92c92
< 354.0 339.0 266.0 367.0 256.0 389.0 240.0
---
> 354.0 339.0 266.0 367.0 256.0 389.0 241.0
95c95
< 493.0 289.0 565.0 322.0 454.0 397.0 349.0
---
> 493.0 289.0 564.0 322.0 454.0 397.0 349.0
98c98
< 419.0 464.0 434.0 495.0 397.0 598.0 589.0
---
> 419.0 464.0 434.0 495.0 497.0 598.0 589.0
102,103c102,103
< 407.0 478.0 431.0 694.0 633.0 587.0 435.0
< 527.0 549.0 604.0 664.0 480.0 451.0 419.0
---
> 407.0 478.0 431.0 694.0 632.0 587.0 435.0
> 526.0 549.0 604.0 664.0 480.0 451.0 419.0
105,131c105,133
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
< 1 1 1 1 1 1 1
---
> 0.464380 0.480252 0.426431 0.419414 0.472768 0.430032 0.391296
> 0.395439 0.347502 0.339528 0.358516 0.366492 0.330258 0.409295
> 0.328884 0.322200 0.369123 0.322764 0.335218 0.339731 0.292846
> 0.386128 0.377557 0.420664 0.431206 0.403969 0.394914 0.393565
> 0.606290 0.463429 0.319667 0.465732 0.494974 0.493528 0.418442
> 0.510360 0.482645 0.480140 0.483297 0.503069 0.463225 0.439524
> 0.430395 0.446047 0.409424 0.409756 0.437363 0.362555 0.466882
> 0.423235 0.496694 0.483720 0.377215 0.385568 0.393544 0.417672
> 0.359835 0.564724 0.453687 0.482026 0.612437 0.513056 0.581336
> 0.582363 0.568828 0.570801 0.546519 0.492148 0.505277 0.489650
> 0.445866 0.432672 0.480349 0.471082 0.407719 0.503589 0.521233
> 0.547584 0.549862 0.471185 0.394080 0.543942 0.518961 0.621431
> 0.156339 0.496587 0.517081 0.587237 0.567797 0.603410 0.603324
> 0.550720 0.564284 0.553806 0.552481 0.528222 0.527040 0.507912
> 0.474999 1.409260 0.509266 0.490966 0.498538 0.487384 0.527075
> 0.519088 0.479362 0.416530 0.483688 0.402789 0.357542 0.469409
> 0.520351 0.569883 0.503047 0.578725 0.556899 0.551918 0.555821
> 0.588727 0.542551 0.539876 0.544335 0.474186 0.529514 0.451148
> 0.522177 0.466971 0.459936 1.255981 0.476657 0.509346 0.501497
> 0.517072 0.418300 0.410303 0.465548 0.432661 0.533832 0.372198
> 0.523918 0.533032 0.482876 0.493399 0.497118 0.496355 0.490289
> 0.484170 0.503562 0.433836 0.442722 0.439260 0.416772 0.422372
> 0.421690 0.433204 0.425137 0.427559 0.486981 0.478462 0.416606
> 0.392216 0.508214 0.517432 0.355693 0.475066 0.401692 0.513483
> 0.347694 0.454982 0.397289 0.415619 0.394961 0.376394 0.376200
> 0.402443 0.342085 0.362095 0.352506 0.340282 0.312670 0.370079
> 0.365096 0.353546 0.384143 0.413152 0.425236 0.369764 0.366190
>
>
More information about the d2n-analysis-talk
mailing list