407 const int runno,
double& scale)
415 if (expno == 7) {scale = 1.001932;
return;}
416 if (expno == 9) {scale = 1.001413;
return;}
417 if (expno == 11) {scale = 1.001366;
return;}
419 if (expno == 13) {scale = 1.001460;
return;}
420 if (expno >= 15) {B2ERROR(
"scale_momenta mode=1: not ready for exp >=15");}
425 if (runno >= 0 && runno <= 99) {scale = 1.00262126;
return;}
426 if (runno >= 100 && runno <= 199) {scale = 1.00253537;
return;}
427 if (runno >= 200 && runno <= 299) {scale = 1.00246555;
return;}
428 if (runno >= 300 && runno <= 399) {scale = 1.0026452 ;
return;}
429 if (runno >= 400 && runno <= 499) {scale = 1.00290386;
return;}
430 if (runno >= 500 && runno <= 599) {scale = 1.00287753;
return;}
431 if (runno >= 600 && runno <= 699) {scale = 1.00266845;
return;}
432 if (runno >= 700 && runno <= 799) {scale = 1.00184472;
return;}
433 if (runno >= 800 && runno <= 899) {scale = 1.00200021;
return;}
434 if (runno >= 900 && runno <= 999) {scale = 1.00178485;
return;}
435 if (runno >= 1000 && runno <= 1099) {scale = 1.00196071;
return;}
436 if (runno >= 1100 && runno <= 1199) {scale = 1.00224453;
return;}
437 if (runno >= 1200 && runno <= 1299) {scale = 1.00209138;
return;}
438 if (runno >= 1300 && runno <= 1399) {scale = 1.00145455;
return;}
439 if (runno >= 1400 && runno <= 1499) {scale = 1.00119243;
return;}
440 if (runno >= 1500 && runno <= 1599) {scale = 1.00115572;
return;}
441 if (runno >= 1600 && runno <= 1699) {scale = 1.00133916;
return;}
442 if (runno >= 1700 && runno <= 1799) {scale = 1.00186441;
return;}
443 if (runno >= 1800 && runno <= 1899) {scale = 1.00206786;
return;}
444 if (runno >= 1900 && runno <= 1999) {scale = 1.00184794;
return;}
445 if (runno >= 2000 && runno <= 2099) {scale = 1.00187219;
return;}
446 if (runno >= 2100 && runno <= 2199) {scale = 1.00214844;
return;}
447 if (runno >= 2200 && runno <= 2299) {scale = 1.00154502;
return;}
448 if (runno >= 2300 && runno <= 2399) {scale = 1.00179518;
return;}
449 if (runno >= 2400 && runno <= 2499) {scale = 1.00198983;
return;}
450 if (runno >= 2500 && runno <= 2599) {scale = 1.00219242;
return;}
451 if (runno >= 2600 && runno <= 2699) {scale = 1.00227363;
return;}
452 if (runno >= 2700 && runno <= 2799) {scale = 1.00202328;
return;}
453 if (runno >= 2800 && runno <= 2899) {scale = 1.00193147;
return;}
454 }
else if (expno == 9) {
455 if (runno >= 0 && runno <= 99) {scale = 1.00126754;
return;}
456 if (runno >= 100 && runno <= 199) {scale = 1.00119683;
return;}
457 if (runno >= 200 && runno <= 299) {scale = 1.00123771;
return;}
458 if (runno >= 300 && runno <= 399) {scale = 1.00142625;
return;}
459 if (runno >= 400 && runno <= 499) {scale = 1.00134481;
return;}
460 if (runno >= 500 && runno <= 599) {scale = 1.00131973;
return;}
461 if (runno >= 600 && runno <= 699) {scale = 1.00145963;
return;}
462 if (runno >= 700 && runno <= 799) {scale = 1.00161304;
return;}
463 if (runno >= 800 && runno <= 899) {scale = 1.00141795;
return;}
464 if (runno >= 900 && runno <= 999) {scale = 1.00137157;
return;}
465 if (runno >= 1000 && runno <= 1099) {scale = 1.00126209;
return;}
466 if (runno >= 1100 && runno <= 1199) {scale = 1.00147941;
return;}
467 if (runno >= 1200 && runno <= 1299) {scale = 1.00145555;
return;}
468 }
else if (expno == 11) {
469 if (runno >= 0 && runno <= 99) {scale = 1.00149576;
return;}
470 if (runno >= 100 && runno <= 199) {scale = 1.00155153;
return;}
471 if (runno >= 200 && runno <= 299) {scale = 1.0013597 ;
return;}
472 if (runno >= 300 && runno <= 399) {scale = 1.00147504;
return;}
473 if (runno >= 400 && runno <= 499) {scale = 1.00138242;
return;}
474 if (runno >= 500 && runno <= 599) {scale = 1.00144554;
return;}
475 if (runno >= 600 && runno <= 699) {scale = 1.00136961;
return;}
476 if (runno >= 700 && runno <= 799) {scale = 1.00141443;
return;}
477 if (runno >= 800 && runno <= 899) {scale = 1.00136617;
return;}
478 if (runno >= 900 && runno <= 999) {scale = 1.00126262;
return;}
479 if (runno >= 1000 && runno <= 1099) {scale = 1.00133337;
return;}
480 if (runno >= 1100 && runno <= 1199) {scale = 1.00137906;
return;}
481 if (runno >= 1200 && runno <= 1299) {scale = 1.00127944;
return;}
483 if (runno >= 1300 && runno <= 1399) {scale = 1.00127944;
return;}
484 }
else if (expno == 13) {
485 if (runno >= 0 && runno <= 99) {scale = 1.00110691;
return;}
486 if (runno >= 100 && runno <= 199) {scale = 1.00106123;
return;}
487 if (runno >= 200 && runno <= 299) {scale = 1.00109934;
return;}
488 if (runno >= 300 && runno <= 399) {scale = 1.00105759;
return;}
489 if (runno >= 400 && runno <= 499) {scale = 1.000986887;
return;}
490 if (runno >= 500 && runno <= 599) {scale = 1.000928764;
return;}
491 if (runno >= 600 && runno <= 699) {scale = 1.00103925;
return;}
492 if (runno >= 700 && runno <= 799) {scale = 1.0010591 ;
return;}
493 if (runno >= 800 && runno <= 899) {scale = 1.00127043;
return;}
494 if (runno >= 900 && runno <= 999) {scale = 1.00154033;
return;}
495 if (runno >= 1000 && runno <= 1099) {scale = 1.00180656;
return;}
496 if (runno >= 1100 && runno <= 1199) {scale = 1.00202059;
return;}
497 if (runno >= 1200 && runno <= 1299) {scale = 1.00184308;
return;}
498 if (runno >= 1300 && runno <= 1399) {scale = 1.00180153;
return;}
499 if (runno >= 1400 && runno <= 1499) {scale = 1.00189577;
return;}
500 if (runno >= 1500 && runno <= 1599) {scale = 1.00176026;
return;}
502 if (runno >= 1600 && runno <= 1699) {scale = 1.00176026;
return;}
504 }
else if (expno == 15) {
505 if (runno >= 0 && runno <= 99) {scale = 1.00178719;
return;}
506 if (runno >= 100 && runno <= 199) {scale = 1.00205712;
return;}
507 if (runno >= 200 && runno <= 299) {scale = 1.00197622;
return;}
508 if (runno >= 300 && runno <= 399) {scale = 1.00191558;
return;}
509 if (runno >= 400 && runno <= 499) {scale = 1.00207795;
return;}
510 if (runno >= 500 && runno <= 599) {scale = 1.00191871;
return;}
511 if (runno >= 600 && runno <= 699) {scale = 1.00187917;
return;}
512 if (runno >= 700 && runno <= 799) {scale = 1.0019573 ;
return;}
513 if (runno >= 800 && runno <= 899) {scale = 1.00196562;
return;}
514 if (runno >= 900 && runno <= 999) {scale = 1.00205215;
return;}
515 if (runno >= 1000 && runno <= 1099) {scale = 1.00200559;
return;}
516 if (runno >= 1100 && runno <= 1199) {scale = 1.00198568;
return;}
517 if (runno >= 1200 && runno <= 1299) {scale = 1.00212539;
return;}
518 if (runno >= 1300 && runno <= 1399) {scale = 1.00199158;
return;}
519 if (runno >= 1400 && runno <= 1499) {scale = 1.00208257;
return;}
521 }
else if (expno >= 17) { B2ERROR(
"scale_momenta mode=2: not ready for exp >=17");}
529 const int runno,
double& scale)
538 B2ERROR(
"scale_momenta mode=1 not ready for exp " << expno);
544 if (runno >= 0 && runno <= 99) {scale = 1.00284313;
return;}
545 if (runno >= 100 && runno <= 199) {scale = 1.00272563;
return;}
546 if (runno >= 200 && runno <= 299) {scale = 1.00276953;
return;}
547 if (runno >= 300 && runno <= 399) {scale = 1.00286722;
return;}
548 if (runno >= 400 && runno <= 499) {scale = 1.00318492;
return;}
549 if (runno >= 500 && runno <= 599) {scale = 1.00304208;
return;}
550 if (runno >= 600 && runno <= 699) {scale = 1.00292413;
return;}
551 if (runno >= 700 && runno <= 799) {scale = 1.00201286;
return;}
552 if (runno >= 800 && runno <= 899) {scale = 1.00225237;
return;}
553 if (runno >= 900 && runno <= 999) {scale = 1.00210375;
return;}
554 if (runno >= 1000 && runno <= 1099) {scale = 1.00213222;
return;}
555 if (runno >= 1100 && runno <= 1199) {scale = 1.00253683;
return;}
556 if (runno >= 1200 && runno <= 1299) {scale = 1.00234023;
return;}
557 if (runno >= 1300 && runno <= 1399) {scale = 1.00168421;
return;}
558 if (runno >= 1400 && runno <= 1499) {scale = 1.00136169;
return;}
559 if (runno >= 1500 && runno <= 1599) {scale = 1.00138431;
return;}
560 if (runno >= 1600 && runno <= 1699) {scale = 1.00154521;
return;}
561 if (runno >= 1700 && runno <= 1799) {scale = 1.00204405;
return;}
562 if (runno >= 1800 && runno <= 1899) {scale = 1.00223498;
return;}
563 if (runno >= 1900 && runno <= 1999) {scale = 1.00198975;
return;}
564 if (runno >= 2000 && runno <= 2099) {scale = 1.00201868;
return;}
565 if (runno >= 2100 && runno <= 2199) {scale = 1.00233875;
return;}
566 if (runno >= 2200 && runno <= 2299) {scale = 1.00175269;
return;}
567 if (runno >= 2300 && runno <= 2399) {scale = 1.00192858;
return;}
568 if (runno >= 2400 && runno <= 2499) {scale = 1.00217693;
return;}
569 if (runno >= 2500 && runno <= 2599) {scale = 1.00209337;
return;}
570 if (runno >= 2600 && runno <= 2699) {scale = 1.00244277;
return;}
571 if (runno >= 2700 && runno <= 2799) {scale = 1.00221485;
return;}
572 if (runno >= 2800 && runno <= 2899) {scale = 1.00214678;
return;}
573 }
else if (expno == 9) {
574 if (runno >= 0 && runno <= 99) {scale = 1.00127702;
return;}
575 if (runno >= 100 && runno <= 199) {scale = 1.00119957;
return;}
576 if (runno >= 200 && runno <= 299) {scale = 1.00132677;
return;}
577 if (runno >= 300 && runno <= 399) {scale = 1.0014785 ;
return;}
578 if (runno >= 400 && runno <= 499) {scale = 1.00144872;
return;}
579 if (runno >= 500 && runno <= 599) {scale = 1.00145376;
return;}
580 if (runno >= 600 && runno <= 699) {scale = 1.00151659;
return;}
581 if (runno >= 700 && runno <= 799) {scale = 1.00167384;
return;}
582 if (runno >= 800 && runno <= 899) {scale = 1.00153754;
return;}
583 if (runno >= 900 && runno <= 999) {scale = 1.0014984 ;
return;}
584 if (runno >= 1000 && runno <= 1099) {scale = 1.0013764 ;
return;}
585 if (runno >= 1100 && runno <= 1199) {scale = 1.00145545;
return;}
586 if (runno >= 1200 && runno <= 1299) {scale = 1.0017164 ;
return;}
587 }
else if (expno == 11) {
588 if (runno >= 0 && runno <= 99) {scale = 1.00159624;
return;}
589 if (runno >= 100 && runno <= 199) {scale = 1.00153365;
return;}
590 if (runno >= 200 && runno <= 299) {scale = 1.00143809;
return;}
591 if (runno >= 300 && runno <= 399) {scale = 1.00148495;
return;}
592 if (runno >= 400 && runno <= 499) {scale = 1.00143193;
return;}
593 if (runno >= 500 && runno <= 599) {scale = 1.00141905;
return;}
594 if (runno >= 600 && runno <= 699) {scale = 1.0013992 ;
return;}
595 if (runno >= 700 && runno <= 799) {scale = 1.00148836;
return;}
596 if (runno >= 800 && runno <= 899) {scale = 1.00139019;
return;}
597 if (runno >= 900 && runno <= 999) {scale = 1.00151389;
return;}
598 if (runno >= 1000 && runno <= 1099) {scale = 1.00142988;
return;}
599 if (runno >= 1100 && runno <= 1199) {scale = 1.00143519;
return;}
600 if (runno >= 1200 && runno <= 1299) {scale = 1.00132641;
return;}
601 }
else if (expno == 13) {
602 if (runno >= 0 && runno <= 99) {scale = 1.00126317;
return;}
603 if (runno >= 100 && runno <= 199) {scale = 1.00118837;
return;}
604 if (runno >= 200 && runno <= 299) {scale = 1.00122276;
return;}
605 if (runno >= 300 && runno <= 399) {scale = 1.00113233;
return;}
606 if (runno >= 400 && runno <= 499) {scale = 1.00111677;
return;}
607 if (runno >= 500 && runno <= 599) {scale = 1.00100963;
return;}
608 if (runno >= 600 && runno <= 699) {scale = 1.0011474 ;
return;}
609 if (runno >= 700 && runno <= 799) {scale = 1.00116245;
return;}
610 if (runno >= 800 && runno <= 899) {scale = 1.00139826;
return;}
611 if (runno >= 900 && runno <= 999) {scale = 1.0017014 ;
return;}
612 if (runno >= 1000 && runno <= 1099) {scale = 1.00190706;
return;}
613 if (runno >= 1100 && runno <= 1199) {scale = 1.00214688;
return;}
614 if (runno >= 1200 && runno <= 1299) {scale = 1.00207336;
return;}
615 if (runno >= 1300 && runno <= 1399) {scale = 1.00192885;
return;}
616 if (runno >= 1400 && runno <= 1499) {scale = 1.00196289;
return;}
617 if (runno >= 1500 && runno <= 1599) {scale = 1.00185017;
return;}
618 if (runno >= 1600 && runno <= 1699) {scale = 1.00191256;
return;}
619 }
else if (expno == 15) {
620 if (runno >= 0 && runno <= 99) {scale = 1.001953;
return;}
621 if (runno >= 100 && runno <= 199) {scale = 1.00210508;
return;}
622 if (runno >= 200 && runno <= 299) {scale = 1.002139;
return;}
623 if (runno >= 300 && runno <= 399) {scale = 1.00207803;
return;}
624 if (runno >= 400 && runno <= 499) {scale = 1.0022331;
return;}
625 if (runno >= 500 && runno <= 599) {scale = 1.00206421;
return;}
626 if (runno >= 600 && runno <= 699) {scale = 1.00205153;
return;}
627 if (runno >= 700 && runno <= 799) {scale = 1.00214664;
return;}
628 if (runno >= 800 && runno <= 899) {scale = 1.0021456;
return;}
629 if (runno >= 900 && runno <= 999) {scale = 1.00221055;
return;}
630 if (runno >= 1000 && runno <= 1099) {scale = 1.00216472;
return;}
631 if (runno >= 1100 && runno <= 1199) {scale = 1.00214208;
return;}
632 if (runno >= 1200 && runno <= 1299) {scale = 1.00228874;
return;}
633 if (runno >= 1300 && runno <= 1399) {scale = 1.00215047;
return;}
634 if (runno >= 1400 && runno <= 1499) {scale = 1.00219808;
return;}
635 }
else if (expno == 17) {
636 if (runno >= 0 && runno <= 99) {scale = 1.00216742;
return;}
637 if (runno >= 100 && runno <= 199) {scale = 1.00229433;
return;}
638 if (runno >= 200 && runno <= 299) {scale = 1.0022749 ;
return;}
639 if (runno >= 300 && runno <= 399) {scale = 1.00231283;
return;}
640 if (runno >= 400 && runno <= 499) {scale = 1.00223654;
return;}
641 if (runno >= 500 && runno <= 599) {scale = 1.00236524;
return;}
642 if (runno >= 600 && runno <= 699) {scale = 1.00229312;
return;}
643 if (runno >= 700 && runno <= 799) {scale = 1.002262 ;
return;}
644 if (runno >= 800 && runno <= 899) {scale = 1.00220104;
return;}
645 if (runno >= 900 && runno <= 999) {scale = 1.00215475;
return;}
646 }
else if (expno == 19) {
647 if (runno >= 0 && runno <= 99) {scale = 1.00227682;
return;}
648 if (runno >= 100 && runno <= 199) {scale = 1.00235795;
return;}
649 if (runno >= 200 && runno <= 299) {scale = 1.00251509;
return;}
650 if (runno >= 300 && runno <= 399) {scale = 1.0022106 ;
return;}
651 if (runno >= 400 && runno <= 499) {scale = 1.00230963;
return;}
652 if (runno >= 500 && runno <= 599) {scale = 1.002229 ;
return;}
653 if (runno >= 600 && runno <= 699) {scale = 1.00215402;
return;}
654 if (runno >= 700 && runno <= 799) {scale = 1.00217892;
return;}
655 if (runno >= 800 && runno <= 899) {scale = 1.0022305 ;
return;}
656 if (runno >= 900 && runno <= 999) {scale = 1.0021458 ;
return;}
657 if (runno >= 1000 && runno <= 1099) {scale = 1.00215792;
return;}
658 if (runno >= 1100 && runno <= 1199) {scale = 1.0022667 ;
return;}
659 if (runno >= 1200 && runno <= 1299) {scale = 1.00219581;
return;}
660 if (runno >= 1300 && runno <= 1399) {scale = 1.00218378;
return;}
661 if (runno >= 1400 && runno <= 1499) {scale = 1.00232172;
return;}
662 if (runno >= 1500 && runno <= 1599) {scale = 1.00217172;
return;}
663 if (runno >= 1600 && runno <= 1700) {scale = 1.00208639;
return;}
664 if (runno >= 1700 && runno <= 1799) {scale = 1.00188253;
return;}
666 B2ERROR(
"scale_momenta not ready for this exp,run "
667 << expno <<
"," << runno);
669 }
else if (expno == 21) {
670 if (runno >= 0 && runno <= 99) {scale = 1.00228182;
return;}
671 if (runno >= 100 && runno <= 199) {scale = 1.00233446;
return;}
672 if (runno >= 200 && runno <= 299) {scale = 1.00238846;
return;}
673 if (runno >= 300 && runno <= 399) {scale = 1.00238061 ;
return;}
675 B2ERROR(
"scale_momenta not ready for this exp,run "
676 << expno <<
"," << runno);
678 }
else if (expno == 23) {
679 if (runno >= 0 && runno <= 99) {scale = 1.00216089;
return;}
680 if (runno >= 100 && runno <= 199) {scale = 1.0019783;
return;}
681 if (runno >= 200 && runno <= 299) {scale = 1.00172162;
return;}
682 if (runno >= 300 && runno <= 399) {scale = 1.0017938 ;
return;}
683 if (runno >= 400 && runno <= 499) {scale = 1.00177832 ;
return;}
684 if (runno >= 500 && runno <= 599) {scale = 1.0017609 ;
return;}
685 if (runno >= 600 && runno <= 699) {scale = 1.0017756 ;
return;}
687 B2ERROR(
"scale_momenta not ready for this exp,run "
688 << expno <<
"," << runno);
690 }
else if (expno == 25) {
691 if (runno >= 0 && runno <= 99) {scale = 1.00190068;
return;}
692 if (runno >= 100 && runno <= 199) {scale = 1.00199038;
return;}
693 if (runno >= 200 && runno <= 299) {scale = 1.00196171;
return;}
694 if (runno >= 300 && runno <= 399) {scale = 1.00200167;
return;}
695 if (runno >= 400 && runno <= 499) {scale = 1.0019667 ;
return;}
696 if (runno >= 500 && runno <= 599) {scale = 1.00212109;
return;}
697 if (runno >= 600 && runno <= 699) {scale = 1.00202115;
return;}
698 if (runno >= 700 && runno <= 799) {scale = 1.00196456;
return;}
699 if (runno >= 800 && runno <= 899) {scale = 1.00200392;
return;}
700 if (runno >= 900 && runno <= 999) {scale = 1.00199495;
return;}
701 if (runno >= 1000 && runno <= 1099) {scale = 1.00202212;
return;}
702 if (runno >= 1100 && runno <= 1199) {scale = 1.00205312;
return;}
703 if (runno >= 1200 && runno <= 1299) {scale = 1.00200891;
return;}
704 if (runno >= 1300 && runno <= 1399) {scale = 1.00196866;
return;}
705 if (runno >= 1400 && runno <= 1499) {scale = 1.0020257 ;
return;}
706 if (runno >= 1500 && runno <= 1599) {scale = 1.00196995;
return;}
707 if (runno >= 1600 && runno <= 1699) {scale = 1.00189232;
return;}
708 if (runno >= 1700 && runno <= 1799) {scale = 1.00179028;
return;}
709 if (runno >= 1800 && runno <= 1899) {scale = 1.00179459;
return;}
710 if (runno >= 1900 && runno <= 1999) {scale = 1.001762 ;
return;}
711 if (runno >= 2000 && runno <= 2099) {scale = 1.00179514;
return;}
712 if (runno >= 2100 && runno <= 2199) {scale = 1.00175211;
return;}
714 B2ERROR(
"scale_momenta not ready for this exp,run "
715 << expno <<
"," << runno);
717 }
else if (expno == 27) {
718 if (runno >= 0 && runno <= 99) {scale = 1.00192647;
return;}
719 if (runno >= 100 && runno <= 199) {scale = 1.0019717 ;
return;}
720 if (runno >= 200 && runno <= 299) {scale = 1.00196428;
return;}
721 if (runno >= 300 && runno <= 399) {scale = 1.00190998;
return;}
722 if (runno >= 400 && runno <= 499) {scale = 1.00204645;
return;}
723 if (runno >= 500 && runno <= 599) {scale = 1.0020687 ;
return;}
724 if (runno >= 600 && runno <= 699) {scale = 1.00198209;
return;}
726 if (runno >= 700 && runno <= 800) {scale = 1.00187505;
return;}
727 if (runno >= 801 && runno <= 900) {scale = 1.00186805;
return;}
728 if (runno >= 901 && runno <= 1000) {scale = 1.00185576;
return;}
729 if (runno >= 1001 && runno <= 1079) {scale = 1.00177176;
return;}
730 if (runno >= 1080 && runno <= 1100) {scale = 1.00177184;
return;}
731 if (runno >= 1101 && runno <= 1200) {scale = 1.00174057;
return;}
732 if (runno >= 1201 && runno <= 1251) {scale = 1.00189649;
return;}
733 if (runno >= 1252 && runno <= 1300) {scale = 1.00189999;
return;}
734 if (runno >= 1301 && runno <= 1400) {scale = 1.0018818 ;
return;}
735 if (runno >= 1401 && runno <= 1500) {scale = 1.00200148;
return;}
736 if (runno >= 1501 && runno <= 1600) {scale = 1.00220137;
return;}
737 if (runno >= 1601 && runno <= 1632) {scale = 1.00217034;
return;}
739 B2ERROR(
"scale_momenta not ready for this exp,run "
740 << expno <<
"," << runno);
742 }
else if (expno == 31) {
743 if (runno >= 0 && runno <= 99) {scale = 1.000373374;
return;}
744 if (runno >= 100 && runno <= 199) {scale = 1.000684843;
return;}
745 if (runno >= 200 && runno <= 299) {scale = 1.000741216;
return;}
746 if (runno >= 300 && runno <= 399) {scale = 1.00092523 ;
return;}
747 if (runno >= 400 && runno <= 499) {scale = 1.00104576 ;
return;}
748 if (runno >= 500 && runno <= 599) {scale = 1.00103982 ;
return;}
749 if (runno >= 600 && runno <= 699) {scale = 1.0010601 ;
return;}
750 if (runno >= 700 && runno <= 799) {scale = 1.000982515;
return;}
751 if (runno >= 800 && runno <= 899) {scale = 1.00101577 ;
return;}
752 if (runno >= 900 && runno <= 999) {scale = 1.000972193;
return;}
753 if (runno >= 1000 && runno <= 1099) {scale = 1.00102536 ;
return;}
754 if (runno >= 1100 && runno <= 1199) {scale = 1.00127571 ;
return;}
755 if (runno >= 1200 && runno <= 1299) {scale = 1.00117938 ;
return;}
756 if (runno >= 1300 && runno <= 1399) {scale = 1.00121084 ;
return;}
757 if (runno >= 1400 && runno <= 1499) {scale = 1.00125911 ;
return;}
758 if (runno >= 1500 && runno <= 1599) {scale = 1.0012907 ;
return;}
759 if (runno >= 1600 && runno <= 1699) {scale = 1.00125362 ;
return;}
760 if (runno >= 1700 && runno <= 1799) {scale = 1.00131686 ;
return;}
762 B2ERROR(
"scale_momenta not ready for this exp,run "
763 << expno <<
"," << runno);
765 }
else if (expno == 33) {
766 if (runno >= 0 && runno <= 99) {scale = 1.00136407;
return;}
767 if (runno >= 100 && runno <= 199) {scale = 1.0013701 ;
return;}
768 if (runno >= 200 && runno <= 299) {scale = 1.00156244;
return;}
769 if (runno >= 300 && runno <= 399) {scale = 1.00138501;
return;}
770 if (runno >= 400 && runno <= 499) {scale = 1.00148468;
return;}
771 if (runno >= 500 && runno <= 599) {scale = 1.0013797 ;
return;}
772 if (runno >= 600 && runno <= 699) {scale = 1.00152298;
return;}
773 if (runno >= 700 && runno <= 799) {scale = 1.001524 ;
return;}
774 if (runno >= 800 && runno <= 899) {scale = 1.0014635 ;
return;}
776 B2ERROR(
"scale_momenta not ready for this exp,run "
777 << expno <<
"," << runno);
779 }
else if (expno == 35) {
780 if (runno >= 0 && runno <= 99) {scale = 1.00135452;
return;}
781 if (runno >= 100 && runno <= 199) {scale = 1.00135798;
return;}
782 if (runno >= 200 && runno <= 299) {scale = 1.00139979;
return;}
783 if (runno >= 300 && runno <= 399) {scale = 1.00136061;
return;}
784 if (runno >= 400 && runno <= 499) {scale = 1.00135776;
return;}
785 if (runno >= 500 && runno <= 599) {scale = 1.00118666;
return;}
786 if (runno >= 600 && runno <= 699) {scale = 1.0011346 ;
return;}
788 B2ERROR(
"scale_momenta not ready for this exp,run "
789 << expno <<
"," << runno);
791 }
else if (expno == 37) {
792 if (runno >= 0 && runno <= 99) {scale = 1.00076522 ;
return;}
793 if (runno >= 100 && runno <= 199) {scale = 1.000875154;
return;}
794 if (runno >= 200 && runno <= 299) {scale = 1.000946763;
return;}
795 if (runno >= 300 && runno <= 399) {scale = 1.000868444;
return;}
796 if (runno >= 400 && runno <= 499) {scale = 1.000781409;
return;}
797 if (runno >= 500 && runno <= 599) {scale = 1.000781171;
return;}
798 if (runno >= 600 && runno <= 699) {scale = 1.000857962;
return;}
799 if (runno >= 700 && runno <= 799) {scale = 1.000756082;
return;}
800 if (runno >= 800 && runno <= 899) {scale = 1.000727355;
return;}
801 if (runno >= 900 && runno <= 999) {scale = 1.000740236;
return;}
802 if (runno >= 1000 && runno <= 1099) {scale = 1.000499606;
return;}
803 if (runno >= 1100 && runno <= 1199) {scale = 1.000481664;
return;}
804 if (runno >= 1200 && runno <= 1299) {scale = 1.000706924;
return;}
805 if (runno >= 1300 && runno <= 1399) {scale = 1.000673738;
return;}
806 if (runno >= 1400 && runno <= 1499) {scale = 1.000662648;
return;}
807 if (runno >= 1500 && runno <= 1599) {scale = 1.000671198;
return;}
808 if (runno >= 1600 && runno <= 1699) {scale = 1.000604364;
return;}
809 if (runno >= 1700 && runno <= 1799) {scale = 1.000717372;
return;}
810 if (runno >= 1800 && runno <= 1899) {scale = 1.000512399;
return;}
811 if (runno >= 1900 && runno <= 1999) {scale = 1.000436958;
return;}
813 B2ERROR(
"scale_momenta not ready for this exp,run "
814 << expno <<
"," << runno);
816 }
else if (expno == 39) {
817 if (runno >= 0 && runno <= 99) {scale = 1.000504342;
return;}
818 if (runno >= 100 && runno <= 199) {scale = 1.000704544;
return;}
819 if (runno >= 200 && runno <= 299) {scale = 1.00094335 ;
return;}
820 if (runno >= 300 && runno <= 399) {scale = 1.000928819;
return;}
821 if (runno >= 400 && runno <= 499) {scale = 1.000884638;
return;}
822 if (runno >= 500 && runno <= 599) {scale = 1.00083459 ;
return;}
823 if (runno >= 600 && runno <= 699) {scale = 1.000767604;
return;}
824 if (runno >= 700 && runno <= 799) {scale = 1.000882219;
return;}
825 if (runno >= 800 && runno <= 899) {scale = 1.000781437;
return;}
826 if (runno >= 900 && runno <= 999) {scale = 1.000853168;
return;}
827 if (runno >= 1000 && runno <= 1099) {scale = 1.000926527;
return;}
828 if (runno >= 1100 && runno <= 1199) {scale = 1.000942882;
return;}
829 if (runno >= 1200 && runno <= 1299) {scale = 1.000932802;
return;}
830 if (runno >= 1300 && runno <= 1399) {scale = 1.000898892;
return;}
832 B2ERROR(
"scale_momenta not ready for this exp,run "
833 << expno <<
"," << runno);
835 }
else if (expno == 41) {
836 if (runno >= 0 && runno <= 99) {scale = 1.00178427;
return;}
837 if (runno >= 100 && runno <= 199) {scale = 1.00188559;
return;}
838 if (runno >= 200 && runno <= 299) {scale = 1.0019292 ;
return;}
839 if (runno >= 300 && runno <= 399) {scale = 1.00196352;
return;}
840 if (runno >= 400 && runno <= 499) {scale = 1.0019078 ;
return;}
841 if (runno >= 500 && runno <= 599) {scale = 1.00185598;
return;}
842 if (runno >= 600 && runno <= 699) {scale = 1.00191314;
return;}
843 if (runno >= 700 && runno <= 799) {scale = 1.00179647;
return;}
844 if (runno >= 800 && runno <= 899) {scale = 1.00189776;
return;}
845 if (runno >= 900 && runno <= 999) {scale = 1.00184798;
return;}
846 if (runno >= 1000 && runno <= 1099) {scale = 1.00177963;
return;}
847 if (runno >= 1100 && runno <= 1199) {scale = 1.00176408;
return;}
848 if (runno >= 1200 && runno <= 1299) {scale = 1.00171371;
return;}
850 B2ERROR(
"scale_momenta not ready for this exp,run "
851 << expno <<
"," << runno);
853 }
else if (expno == 43) {
854 if (runno >= 0 && runno <= 99) {scale = 1.00142307;
return;}
855 if (runno >= 100 && runno <= 199) {scale = 1.000979455;
return;}
856 if (runno >= 200 && runno <= 299) {scale = 1.000974458;
return;}
857 if (runno >= 300 && runno <= 399) {scale = 1.00103301;
return;}
858 if (runno >= 400 && runno <= 499) {scale = 1.001111994;
return;}
859 if (runno >= 500 && runno <= 599) {scale = 1.00100635;
return;}
860 if (runno >= 600 && runno <= 699) {scale = 1.00105078;
return;}
861 if (runno >= 700 && runno <= 799) {scale = 1.00103593;
return;}
862 if (runno >= 800 && runno <= 899) {scale = 1.00105158;
return;}
863 if (runno >= 900 && runno <= 999) {scale = 1.000955608;
return;}
864 if (runno >= 1000 && runno <= 1099) {scale = 1.00099199;
return;}
865 if (runno >= 1100 && runno <= 1199) {scale = 1.0011439;
return;}
867 B2ERROR(
"scale_momenta not ready for this exp,run "
868 << expno <<
"," << runno);
870 }
else if (expno == 45) {
871 if (runno >= 0 && runno <= 99) {scale = 1.00126261;
return;}
872 if (runno >= 100 && runno <= 199) {scale = 1.00138601;
return;}
873 if (runno >= 200 && runno <= 299) {scale = 1.00135372;
return;}
874 if (runno >= 300 && runno <= 399) {scale = 1.00141286;
return;}
875 if (runno >= 400 && runno <= 499) {scale = 1.00147822;
return;}
877 B2ERROR(
"scale_momenta not ready for this exp,run "
878 << expno <<
"," << runno);
880 }
else if (expno == 47) {
881 if (runno >= 0 && runno <= 99) {scale = 1.00156977;
return;}
882 if (runno >= 100 && runno <= 199) {scale = 1.00155614;
return;}
883 if (runno >= 200 && runno <= 299) {scale = 1.0016555;
return;}
884 if (runno >= 300 && runno <= 399) {scale = 1.00167046;
return;}
885 if (runno >= 400 && runno <= 499) {scale = 1.00168705;
return;}
886 if (runno >= 500 && runno <= 599) {scale = 1.00169555;
return;}
887 if (runno >= 600 && runno <= 699) {scale = 1.00175653;
return;}
888 if (runno >= 700 && runno <= 799) {scale = 1.00174358;
return;}
889 if (runno >= 800 && runno <= 899) {scale = 1.00174004;
return;}
891 B2ERROR(
"scale_momenta not ready for this exp,run "
892 << expno <<
"," << runno);
894 }
else if (expno == 49) {
895 if (runno >= 0 && runno <= 99) {scale = 1.00158837;
return;}
896 if (runno >= 100 && runno <= 199) {scale = 1.00163884;
return;}
897 if (runno >= 200 && runno <= 299) {scale = 1.00160595;
return;}
898 if (runno >= 300 && runno <= 399) {scale = 1.00149916;
return;}
899 if (runno >= 400 && runno <= 499) {scale = 1.0014956 ;
return;}
900 if (runno >= 500 && runno <= 599) {scale = 1.00156212;
return;}
901 if (runno >= 600 && runno <= 699) {scale = 1.00121868;
return;}
902 if (runno >= 700 && runno <= 799) {scale = 1.00134613;
return;}
903 if (runno >= 800 && runno <= 899) {scale = 1.00138985;
return;}
904 if (runno >= 900 && runno <= 999) {scale = 1.00129356;
return;}
905 if (runno >= 1000 && runno <= 1099) {scale = 1.00119732;
return;}
906 if (runno >= 1100 && runno <= 1199) {scale = 1.00121481;
return;}
907 if (runno >= 1200 && runno <= 1299) {scale = 1.00121108;
return;}
909 B2ERROR(
"scale_momenta not ready for this exp,run "
910 << expno <<
"," << runno);
912 }
else if (expno == 51) {
913 if (runno >= 0 && runno <= 99) {scale = 1.00160252;
return;}
914 if (runno >= 100 && runno <= 199) {scale = 1.00156099;
return;}
915 if (runno >= 200 && runno <= 299) {scale = 1.00154760;
return;}
916 if (runno >= 300 && runno <= 399) {scale = 1.00146316;
return;}
917 if (runno >= 400 && runno <= 499) {scale = 1.00145525;
return;}
918 if (runno >= 500 && runno <= 599) {scale = 1.00134429;
return;}
919 if (runno >= 600 && runno <= 699) {scale = 1.00135581;
return;}
920 if (runno >= 700 && runno <= 799) {scale = 1.00134382;
return;}
921 if (runno >= 800 && runno <= 899) {scale = 1.00126462;
return;}
922 if (runno >= 900 && runno <= 999) {scale = 1.00130752;
return;}
923 if (runno >= 1000 && runno <= 1099) {scale = 1.00130452;
return;}
924 if (runno >= 1100 && runno <= 1199) {scale = 1.00131440;
return;}
925 if (runno >= 1200 && runno <= 1299) {scale = 1.00130864;
return;}
926 if (runno >= 1300 && runno <= 1399) {scale = 1.00105290;
return;}
927 if (runno >= 1400 && runno <= 1499) {scale = 1.00126645;
return;}
928 if (runno >= 1500 && runno <= 1599) {scale = 1.00126383;
return;}
929 if (runno >= 1600 && runno <= 1699) {scale = 1.00141111;
return;}
930 if (runno >= 1700 && runno <= 1799) {scale = 1.00126220;
return;}
931 if (runno >= 1800 && runno <= 1899) {scale = 1.00105098;
return;}
933 B2ERROR(
"scale_momenta not ready for this exp,run "
934 << expno <<
"," << runno);
936 }
else if (expno == 53) {
937 if (runno >= 0 && runno <= 99) {scale = 1.0011516;
return;}
938 if (runno >= 100 && runno <= 199) {scale = 1.00115527;
return;}
939 if (runno >= 200 && runno <= 299) {scale = 1.00114844;
return;}
941 B2ERROR(
"scale_momenta not ready for this exp,run "
942 << expno <<
"," << runno);
944 }
else if (expno == 55) {
945 if (runno >= 0 && runno <= 99) {scale = 1.00114284;
return;}
946 if (runno >= 100 && runno <= 199) {scale = 1.00111458;
return;}
947 if (runno >= 200 && runno <= 299) {scale = 1.00109686;
return;}
948 if (runno >= 300 && runno <= 399) {scale = 1.00119475;
return;}
949 if (runno >= 400 && runno <= 499) {scale = 1.00117818;
return;}
950 if (runno >= 500 && runno <= 599) {scale = 1.00115789;
return;}
951 if (runno >= 600 && runno <= 699) {scale = 1.00122261;
return;}
952 if (runno >= 700 && runno <= 799) {scale = 1.00118454;
return;}
953 if (runno >= 800 && runno <= 899) {scale = 1.00118042;
return;}
954 if (runno >= 900 && runno <= 999) {scale = 1.00124759;
return;}
955 if (runno >= 1000 && runno <= 1099) {scale = 1.00128055;
return;}
956 if (runno >= 1100 && runno <= 1199) {scale = 1.00119131;
return;}
957 if (runno >= 1200 && runno <= 1299) {scale = 1.00122238;
return;}
958 if (runno >= 1300 && runno <= 1399) {scale = 1.00129538;
return;}
959 if (runno >= 1400 && runno <= 1499) {scale = 1.00130387;
return;}
960 if (runno >= 1500 && runno <= 1599) {scale = 1.00130858;
return;}
961 if (runno >= 1600 && runno <= 1699) {scale = 1.00111854;
return;}
962 if (runno >= 1700 && runno <= 1799) {scale = 1.00136261;
return;}
964 B2ERROR(
"scale_momenta not ready for this exp,run "
965 << expno <<
"," << runno);
967 }
else if (expno == 61) {
968 if (runno >= 0 && runno <= 99) {scale = 1.0009992;
return;}
969 if (runno >= 100 && runno <= 199) {scale = 1.00113704;
return;}
970 if (runno >= 200 && runno <= 299) {scale = 1.00129904;
return;}
971 if (runno >= 300 && runno <= 399) {scale = 1.00141879;
return;}
972 if (runno >= 400 && runno <= 499) {scale = 1.00146707;
return;}
973 if (runno >= 500 && runno <= 599) {scale = 1.00150101;
return;}
974 if (runno >= 600 && runno <= 699) {scale = 1.00147322;
return;}
975 if (runno >= 700 && runno <= 799) {scale = 1.00153929;
return;}
976 if (runno >= 800 && runno <= 899) {scale = 1.00159997;
return;}
977 if (runno >= 900 && runno <= 999) {scale = 1.00164032;
return;}
978 if (runno >= 1000 && runno <= 1099) {scale = 1.00165878;
return;}
979 if (runno >= 1100 && runno <= 1199) {scale = 1.00163475;
return;}
980 if (runno >= 1200 && runno <= 1207) {scale = 1.00166193;
return;}
981 if (runno >= 1208 && runno <= 1299) {scale = 1.00235824;
return;}
982 if (runno >= 1300 && runno <= 1399) {scale = 1.00242282;
return;}
984 B2ERROR(
"scale_momenta not ready for this exp,run "
985 << expno <<
"," << runno);
987 }
else if (expno == 65) {
989 B2ERROR(
"scale_momenta not ready for this exp,run "
990 << expno <<
"," << runno);
992 if (runno >= 1000 && runno <= 1336) {scale = 1.00145234;
return;}
994 B2ERROR(
"scale_momenta not ready for this exp,run "
995 << expno <<
"," << runno);
998 B2ERROR(
"scale_momenta mode=2 not ready for exp " << expno);
1005 const int runno,
double& scale)
1014 B2ERROR(
"scale_momenta mode=1 not ready for exp " << expno);
1021 }
else if (expno == 9) {
1023 }
else if (expno == 11) {
1025 }
else if (expno == 13) {
1027 }
else if (expno == 15) {
1029 }
else if (expno == 17) {
1031 }
else if (expno == 19) {
1032 B2ERROR(
"scale_momenta not ready for this exp,run "
1033 << expno <<
"," << runno);
1035 }
else if (expno == 21) {
1036 B2ERROR(
"scale_momenta not ready for this exp,run "
1037 << expno <<
"," << runno);
1039 }
else if (expno == 23) {
1040 B2ERROR(
"scale_momenta not ready for this exp,run "
1041 << expno <<
"," << runno);
1043 }
else if (expno == 25) {
1044 B2ERROR(
"scale_momenta not ready for this exp,run "
1045 << expno <<
"," << runno);
1047 }
else if (expno == 27) {
1048 B2ERROR(
"scale_momenta not ready for this exp,run "
1049 << expno <<
"," << runno);
1051 }
else if (expno == 31) {
1052 if (runno >= 0 && runno <= 137) {
1053 B2ERROR(
"scale_momenta not ready for this exp,run "
1054 << expno <<
"," << runno);
1056 if (runno >= 138 && runno <= 199) {scale = 1.000931841;
return;}
1057 if (runno >= 200 && runno <= 299) {scale = 1.000916397;
return;}
1058 if (runno >= 300 && runno <= 399) {scale = 1.00108023 ;
return;}
1059 if (runno >= 400 && runno <= 499) {scale = 1.00118662 ;
return;}
1060 if (runno >= 500 && runno <= 599) {scale = 1.00117739 ;
return;}
1061 if (runno >= 600 && runno <= 699) {scale = 1.00119542 ;
return;}
1062 if (runno >= 700 && runno <= 799) {scale = 1.00110396 ;
return;}
1063 if (runno >= 800 && runno <= 899) {scale = 1.00109603 ;
return;}
1064 if (runno >= 900 && runno <= 999) {scale = 1.00112795 ;
return;}
1065 if (runno >= 1000 && runno <= 1099) {scale = 1.00118365 ;
return;}
1066 if (runno >= 1100 && runno <= 1199) {scale = 1.00142214 ;
return;}
1067 if (runno >= 1200 && runno <= 1299) {scale = 1.00133150 ;
return;}
1068 if (runno >= 1300 && runno <= 1399) {scale = 1.00132831 ;
return;}
1069 if (runno >= 1400 && runno <= 1499) {scale = 1.00136554 ;
return;}
1070 if (runno >= 1500 && runno <= 1599) {scale = 1.00141187 ;
return;}
1071 if (runno >= 1600 && runno <= 1699) {scale = 1.00136628 ;
return;}
1072 if (runno >= 1700 && runno <= 1799) {scale = 1.00139273 ;
return;}
1073 if (runno >= 1800) {
1074 B2ERROR(
"scale_momenta not ready for this exp,run "
1075 << expno <<
"," << runno);
1077 }
else if (expno == 33) {
1078 if (runno >= 0 && runno <= 99) {scale = 1.00149319;
return;}
1079 if (runno >= 100 && runno <= 199) {scale = 1.00150915;
return;}
1080 if (runno >= 200 && runno <= 299) {scale = 1.00173040;
return;}
1081 if (runno >= 300 && runno <= 399) {scale = 1.00150449;
return;}
1082 if (runno >= 400 && runno <= 499) {scale = 1.00161519;
return;}
1083 if (runno >= 500 && runno <= 599) {scale = 1.00151670;
return;}
1084 if (runno >= 600 && runno <= 699) {scale = 1.00164347;
return;}
1085 if (runno >= 700 && runno <= 799) {scale = 1.00164165;
return;}
1086 if (runno >= 800 && runno <= 899) {scale = 1.00161369;
return;}
1088 B2ERROR(
"scale_momenta not ready for this exp,run "
1089 << expno <<
"," << runno);
1091 }
else if (expno == 35) {
1092 if (runno >= 0 && runno <= 99) {scale = 1.00147034;
return;}
1093 if (runno >= 100 && runno <= 199) {scale = 1.00148523;
return;}
1094 if (runno >= 200 && runno <= 299) {scale = 1.00153372;
return;}
1095 if (runno >= 300 && runno <= 399) {scale = 1.00148256;
return;}
1096 if (runno >= 400 && runno <= 499) {scale = 1.00144902;
return;}
1097 if (runno >= 500 && runno <= 599) {scale = 1.00131501;
return;}
1098 if (runno >= 600 && runno <= 699) {scale = 1.00126371;
return;}
1100 B2ERROR(
"scale_momenta not ready for this exp,run "
1101 << expno <<
"," << runno);
1103 }
else if (expno == 37) {
1104 if (runno >= 0 && runno <= 99) {scale = 1.000916277;
return;}
1105 if (runno >= 100 && runno <= 199) {scale = 1.001035310;
return;}
1106 if (runno >= 200 && runno <= 299) {scale = 1.001123403;
return;}
1107 if (runno >= 300 && runno <= 399) {scale = 1.001017718;
return;}
1108 if (runno >= 400 && runno <= 499) {scale = 1.000932890;
return;}
1109 if (runno >= 500 && runno <= 599) {scale = 1.000928479;
return;}
1110 if (runno >= 600 && runno <= 699) {scale = 1.000997938;
return;}
1111 if (runno >= 700 && runno <= 799) {scale = 1.000899663;
return;}
1112 if (runno >= 800 && runno <= 899) {scale = 1.000860910;
return;}
1113 if (runno >= 900 && runno <= 999) {scale = 1.000882920;
return;}
1114 if (runno >= 1000 && runno <= 1099) {scale = 1.000616966;
return;}
1115 if (runno >= 1100 && runno <= 1199) {scale = 1.000613018;
return;}
1116 if (runno >= 1200 && runno <= 1299) {scale = 1.000832338;
return;}
1117 if (runno >= 1300 && runno <= 1399) {scale = 1.000803640;
return;}
1118 if (runno >= 1400 && runno <= 1499) {scale = 1.000770454;
return;}
1119 if (runno >= 1500 && runno <= 1599) {scale = 1.000786608;
return;}
1120 if (runno >= 1600 && runno <= 1699) {scale = 1.000718089;
return;}
1121 if (runno >= 1700 && runno <= 1799) {scale = 1.000826042;
return;}
1122 if (runno >= 1800 && runno <= 1899) {scale = 1.000638150;
return;}
1123 if (runno >= 1900 && runno <= 1999) {scale = 1.000529173;
return;}
1124 if (runno >= 2000) {
1125 B2ERROR(
"scale_momenta not ready for this exp,run "
1126 << expno <<
"," << runno);
1128 }
else if (expno == 39) {
1129 if (runno >= 0 && runno <= 99) {scale = 1.000610857;
return;}
1130 if (runno >= 100 && runno <= 199) {scale = 1.000838583;
return;}
1131 if (runno >= 200 && runno <= 299) {scale = 1.00105918 ;
return;}
1132 if (runno >= 300 && runno <= 399) {scale = 1.00105841 ;
return;}
1133 if (runno >= 400 && runno <= 499) {scale = 1.001025523;
return;}
1134 if (runno >= 500 && runno <= 599) {scale = 1.000967373;
return;}
1135 if (runno >= 600 && runno <= 699) {scale = 1.000898585;
return;}
1136 if (runno >= 700 && runno <= 799) {scale = 1.001003199;
return;}
1137 if (runno >= 800 && runno <= 899) {scale = 1.000897072;
return;}
1138 if (runno >= 900 && runno <= 999) {scale = 1.000972551;
return;}
1139 if (runno >= 1000 && runno <= 1099) {scale = 1.001044677;
return;}
1140 if (runno >= 1100 && runno <= 1199) {scale = 1.00106451 ;
return;}
1141 if (runno >= 1200 && runno <= 1299) {scale = 1.00108570 ;
return;}
1142 if (runno >= 1300 && runno <= 1399) {scale = 1.00102381 ;
return;}
1143 if (runno >= 1400) {
1144 B2ERROR(
"scale_momenta not ready for this exp,run "
1145 << expno <<
"," << runno);
1147 }
else if (expno == 41) {
1148 if (runno >= 0 && runno <= 99) {scale = 1.00189378;
return;}
1149 if (runno >= 100 && runno <= 199) {scale = 1.00197304;
return;}
1150 if (runno >= 200 && runno <= 299) {scale = 1.00204049;
return;}
1151 if (runno >= 300 && runno <= 399) {scale = 1.00205065;
return;}
1152 if (runno >= 400 && runno <= 499) {scale = 1.00199205;
return;}
1153 if (runno >= 500 && runno <= 599) {scale = 1.00195618;
return;}
1154 if (runno >= 600 && runno <= 699) {scale = 1.00200889;
return;}
1155 if (runno >= 700 && runno <= 799) {scale = 1.00190365;
return;}
1156 if (runno >= 800 && runno <= 899) {scale = 1.00204192;
return;}
1157 if (runno >= 900 && runno <= 999) {scale = 1.00196542;
return;}
1158 if (runno >= 1000 && runno <= 1099) {scale = 1.00189706;
return;}
1159 if (runno >= 1100 && runno <= 1199) {scale = 1.00187422;
return;}
1160 if (runno >= 1200 && runno <= 1299) {scale = 1.00183714;
return;}
1161 if (runno >= 1300) {
1162 B2ERROR(
"scale_momenta not ready for this exp,run "
1163 << expno <<
"," << runno);
1165 }
else if (expno == 43) {
1166 if (runno >= 0 && runno <= 99) {scale = 1.00151737;
return;}
1167 if (runno >= 100 && runno <= 199) {scale = 1.00110489;
return;}
1168 if (runno >= 200 && runno <= 299) {scale = 1.00108144;
return;}
1169 if (runno >= 300 && runno <= 399) {scale = 1.00114918;
return;}
1170 if (runno >= 400 && runno <= 499) {scale = 1.00122723;
return;}
1171 if (runno >= 500 && runno <= 599) {scale = 1.00111069;
return;}
1172 if (runno >= 600 && runno <= 699) {scale = 1.00115667;
return;}
1173 if (runno >= 700 && runno <= 799) {scale = 1.00113759;
return;}
1174 if (runno >= 800 && runno <= 899) {scale = 1.00115609;
return;}
1175 if (runno >= 900 && runno <= 999) {scale = 1.00105426;
return;}
1179 if (runno >= 1000 && runno <= 1099) {scale = 1.00111210;
return;}
1180 if (runno >= 1100 && runno <= 1199) {scale = 1.00123104;
return;}
1181 if (runno >= 1200) {
1182 B2ERROR(
"scale_momenta not ready for this exp,run "
1183 << expno <<
"," << runno);
1185 }
else if (expno == 45) {
1186 if (runno >= 0 && runno <= 99) {scale = 1.00136477;
return;}
1187 if (runno >= 100 && runno <= 199) {scale = 1.00151600;
return;}
1188 if (runno >= 200 && runno <= 299) {scale = 1.00146757;
return;}
1189 if (runno >= 300 && runno <= 399) {scale = 1.00153299;
return;}
1190 if (runno >= 400 && runno <= 499) {scale = 1.00159018;
return;}
1192 B2ERROR(
"scale_momenta not ready for this exp,run "
1193 << expno <<
"," << runno);
1195 }
else if (expno == 47) {
1196 if (runno >= 0 && runno <= 99) {scale = 1.00166672;
return;}
1197 if (runno >= 100 && runno <= 199) {scale = 1.00165120;
return;}
1198 if (runno >= 200 && runno <= 299) {scale = 1.00175597;
return;}
1199 if (runno >= 300 && runno <= 399) {scale = 1.00177319;
return;}
1200 if (runno >= 400 && runno <= 499) {scale = 1.00179552;
return;}
1201 if (runno >= 500 && runno <= 599) {scale = 1.00179413;
return;}
1202 if (runno >= 600 && runno <= 699) {scale = 1.00186237;
return;}
1203 if (runno >= 700 && runno <= 799) {scale = 1.00183016;
return;}
1204 if (runno >= 800 && runno <= 899) {scale = 1.00184324;
return;}
1206 B2ERROR(
"scale_momenta not ready for this exp,run "
1207 << expno <<
"," << runno);
1209 }
else if (expno == 49) {
1210 if (runno >= 0 && runno <= 99) {scale = 1.00171645;
return;}
1211 if (runno >= 100 && runno <= 199) {scale = 1.00177728;
return;}
1212 if (runno >= 200 && runno <= 299) {scale = 1.00173301;
return;}
1213 if (runno >= 300 && runno <= 399) {scale = 1.00162075;
return;}
1214 if (runno >= 400 && runno <= 499) {scale = 1.00163153;
return;}
1215 if (runno >= 500 && runno <= 599) {scale = 1.00168559;
return;}
1216 if (runno >= 600 && runno <= 699) {scale = 1.00139227;
return;}
1217 if (runno >= 700 && runno <= 799) {scale = 1.00148583;
return;}
1218 if (runno >= 800 && runno <= 899) {scale = 1.00150403;
return;}
1219 if (runno >= 900 && runno <= 999) {scale = 1.00142759;
return;}
1220 if (runno >= 1000 && runno <= 1099) {scale = 1.00134573;
return;}
1221 if (runno >= 1100 && runno <= 1199) {scale = 1.00138313;
return;}
1222 if (runno >= 1200 && runno <= 1299) {scale = 1.00151369;
return;}
1223 if (runno >= 1300) {
1224 B2ERROR(
"scale_momenta not ready for this exp,run "
1225 << expno <<
"," << runno);
1227 }
else if (expno == 51) {
1228 if (runno >= 0 && runno <= 99) {scale = 1.00165035;
return;}
1229 if (runno >= 100 && runno <= 199) {scale = 1.00161504;
return;}
1230 if (runno >= 200 && runno <= 299) {scale = 1.00160162;
return;}
1231 if (runno >= 300 && runno <= 399) {scale = 1.00152725;
return;}
1232 if (runno >= 400 && runno <= 499) {scale = 1.00149943;
return;}
1233 if (runno >= 500 && runno <= 599) {scale = 1.00141294;
return;}
1234 if (runno >= 600 && runno <= 699) {scale = 1.00140154;
return;}
1235 if (runno >= 700 && runno <= 799) {scale = 1.00140759;
return;}
1236 if (runno >= 800 && runno <= 899) {scale = 1.00133671;
return;}
1237 if (runno >= 900 && runno <= 999) {scale = 1.00136792;
return;}
1238 if (runno >= 1000 && runno <= 1099) {scale = 1.00135251;
return;}
1239 if (runno >= 1100 && runno <= 1199) {scale = 1.00138229;
return;}
1240 if (runno >= 1200 && runno <= 1299) {scale = 1.00134938;
return;}
1241 if (runno >= 1300 && runno <= 1399) {scale = 1.00106240;
return;}
1242 if (runno >= 1400 && runno <= 1499) {scale = 1.00132666;
return;}
1243 if (runno >= 1500 && runno <= 1599) {scale = 1.00132654;
return;}
1244 if (runno >= 1600 && runno <= 1699) {scale = 1.00146619;
return;}
1245 if (runno >= 1700 && runno <= 1799) {scale = 1.00131902;
return;}
1246 if (runno >= 1800 && runno <= 1899) {scale = 1.00114243;
return;}
1247 if (runno >= 1900) {
1248 B2ERROR(
"scale_momenta not ready for this exp,run "
1249 << expno <<
"," << runno);
1251 }
else if (expno == 53) {
1252 if (runno >= 0 && runno <= 99) {scale = 1.00125475;
return;}
1253 if (runno >= 100 && runno <= 199) {scale = 1.00124954;
return;}
1254 if (runno >= 200 && runno <= 299) {scale = 1.00122914;
return;}
1256 B2ERROR(
"scale_momenta not ready for this exp,run "
1257 << expno <<
"," << runno);
1259 }
else if (expno == 55) {
1260 if (runno >= 0 && runno <= 99) {scale = 1.00119352;
return;}
1261 if (runno >= 100 && runno <= 199) {scale = 1.00117130;
return;}
1262 if (runno >= 200 && runno <= 299) {scale = 1.00115825;
return;}
1263 if (runno >= 300 && runno <= 399) {scale = 1.00125005;
return;}
1264 if (runno >= 400 && runno <= 499) {scale = 1.00124720;
return;}
1265 if (runno >= 500 && runno <= 599) {scale = 1.00122234;
return;}
1266 if (runno >= 600 && runno <= 699) {scale = 1.00128709;
return;}
1267 if (runno >= 700 && runno <= 799) {scale = 1.00123081;
return;}
1268 if (runno >= 800 && runno <= 899) {scale = 1.00124198;
return;}
1269 if (runno >= 900 && runno <= 999) {scale = 1.00131118;
return;}
1270 if (runno >= 1000 && runno <= 1099) {scale = 1.00132496;
return;}
1271 if (runno >= 1100 && runno <= 1199) {scale = 1.00126186;
return;}
1272 if (runno >= 1200 && runno <= 1299) {scale = 1.00127849;
return;}
1273 if (runno >= 1300 && runno <= 1399) {scale = 1.00135312;
return;}
1274 if (runno >= 1400 && runno <= 1499) {scale = 1.00136637;
return;}
1275 if (runno >= 1500 && runno <= 1599) {scale = 1.00136270;
return;}
1276 if (runno >= 1600 && runno <= 1699) {scale = 1.00118422;
return;}
1277 if (runno >= 1700 && runno <= 1799) {scale = 1.00142667;
return;}
1278 if (runno >= 1800) {
1279 B2ERROR(
"scale_momenta not ready for this exp,run "
1280 << expno <<
"," << runno);
1282 }
else if (expno == 61) {
1283 if (runno >= 0 && runno <= 99) {scale = 1.00103013;
return;}
1284 if (runno >= 100 && runno <= 199) {scale = 1.00116185;
return;}
1285 if (runno >= 200 && runno <= 299) {scale = 1.00133560;
return;}
1286 if (runno >= 300 && runno <= 399) {scale = 1.00145027;
return;}
1287 if (runno >= 400 && runno <= 499) {scale = 1.00147949;
return;}
1288 if (runno >= 500 && runno <= 599) {scale = 1.00151022;
return;}
1289 if (runno >= 600 && runno <= 699) {scale = 1.00150439;
return;}
1290 if (runno >= 700 && runno <= 799) {scale = 1.00155006;
return;}
1291 if (runno >= 800 && runno <= 899) {scale = 1.00162396;
return;}
1292 if (runno >= 900 && runno <= 999) {scale = 1.00168542;
return;}
1293 if (runno >= 1000 && runno <= 1099) {scale = 1.00168249;
return;}
1294 if (runno >= 1100 && runno <= 1207) {scale = 1.00166891;
return;}
1295 if (runno >= 1208 && runno <= 1299) {scale = 1.00249956;
return;}
1296 if (runno >= 1300 && runno <= 1399) {scale = 1.00255134;
return;}
1297 if (runno >= 1400) {
1298 B2ERROR(
"scale_momenta not ready for this exp,run "
1299 << expno <<
"," << runno);
1301 }
else if (expno == 63) {
1302 if (runno >= 0 && runno <= 99) {scale = 1.00129667;
return;}
1303 if (runno >= 100 && runno <= 199) {scale = 1.00123725;
return;}
1304 if (runno >= 200 && runno <= 299) {scale = 1.00126795;
return;}
1305 if (runno >= 300 && runno <= 399) {scale = 1.00122458;
return;}
1306 if (runno >= 400 && runno <= 499) {scale = 1.00116489;
return;}
1307 if (runno >= 500 && runno <= 599) {scale = 1.00116968;
return;}
1308 if (runno >= 600 && runno <= 699) {scale = 1.000918379;
return;}
1309 if (runno >= 700 && runno <= 799) {scale = 1.0010429;
return;}
1311 B2ERROR(
"scale_momenta not ready for this exp,run "
1312 << expno <<
"," << runno);
1314 }
else if (expno == 65) {
1315 if (runno >= 0 && runno <= 99) {scale = 1.00116975;
return;}
1316 if (runno >= 100 && runno <= 199) {scale = 1.00111926;
return;}
1317 if (runno >= 200 && runno <= 299) {scale = 1.00110162;
return;}
1318 if (runno >= 300 && runno <= 399) {scale = 1.00109524;
return;}
1319 if (runno >= 400 && runno <= 499) {scale = 1.00106913;
return;}
1320 if (runno >= 500 && runno <= 599) {scale = 1.00110941;
return;}
1321 if (runno >= 600 && runno <= 699) {scale = 1.000897865;
return;}
1322 if (runno >= 700 && runno <= 999) {scale = 1.00104385;
return;}
1323 if (runno >= 1000 && runno <= 1299) {scale = 1.000876489;
return;}
1324 if (runno >= 1000) {
1325 B2ERROR(
"scale_momenta not ready for this exp,run "
1326 << expno <<
"," << runno);
1328 }
else if (expno == 67) {
1329 if (runno >= 0 && runno <= 199) {scale = 1.000826364;
return;}
1330 if (runno >= 200 && runno <= 299) {scale = 1.000836576;
return;}
1331 if (runno >= 300 && runno <= 399) {scale = 1.000904815;
return;}
1332 if (runno >= 400 && runno <= 499) {scale = 1.000966045;
return;}
1333 if (runno >= 500 && runno <= 599) {scale = 1.000988147;
return;}
1334 if (runno >= 600 && runno <= 699) {scale = 1.000988147;
return;}
1335 if (runno >= 700 && runno <= 742) {scale = 1.000837414;
return;}
1336 if (runno >= 1000 && runno <= 1099) {scale = 1.000984865;
return;}
1337 if (runno >= 1100 && runno <= 1123) {scale = 1.00105248 ;
return;}
1338 if (runno >= 1124) {
1339 B2ERROR(
"scale_momenta not ready for this exp,run "
1340 << expno <<
"," << runno);
1342 }
else if (expno == 69) {
1343 if (runno >= 0 && runno <= 99) {scale = 1.000791450;
return;}
1344 if (runno >= 100 && runno <= 199) {scale = 1.000891748;
return;}
1345 if (runno >= 200 && runno <= 299) {scale = 1.000866165;
return;}
1346 if (runno >= 300 && runno <= 399) {scale = 1.000838834;
return;}
1347 if (runno >= 400 && runno <= 499) {scale = 1.000811878;
return;}
1348 if (runno >= 500 && runno <= 599) {scale = 1.000779810;
return;}
1349 if (runno >= 600 && runno <= 699) {scale = 1.000799086;
return;}
1350 if (runno >= 700 && runno <= 799) {scale = 1.000833797;
return;}
1351 if (runno >= 800 && runno <= 899) {scale = 1.000875203;
return;}
1352 if (runno >= 900 && runno <= 999) {scale = 1.000891998;
return;}
1353 if (runno >= 1000 && runno <= 1099) {scale = 1.000921074;
return;}
1354 if (runno >= 1100 && runno <= 1199) {scale = 1.000900829;
return;}
1355 if (runno >= 1200 && runno <= 1299) {scale = 1.000958405;
return;}
1356 if (runno >= 1300 && runno <= 1399) {scale = 1.000836841;
return;}
1357 if (runno >= 1400) {
1358 B2ERROR(
"scale_momenta not ready for this exp,run "
1359 << expno <<
"," << runno);
1361 }
else if (expno == 71) {
1362 if (runno >= 0 && runno <= 99) {scale = 1.000962999;
return;}
1363 if (runno >= 100 && runno <= 199) {scale = 1.001478932;
return;}
1364 if (runno >= 200 && runno <= 300) {scale = 1.001486524;
return;}
1365 if (runno >= 301 && runno <= 384) {scale = 1.001430843;
return;}
1366 if (runno >= 385 && runno <= 499) {scale = 1.001505696;
return;}
1367 if (runno >= 500 && runno <= 599) {scale = 1.001523980;
return;}
1368 if (runno >= 600 && runno <= 699) {scale = 1.001480830;
return;}
1369 if (runno >= 1000 && runno <= 1013) {scale = 1.001480830;
return;}
1370 if (runno >= 2000 && runno <= 2099) {scale = 1.001617882;
return;}
1371 if (runno >= 2100 && runno <= 2199) {scale = 1.001644395;
return;}
1372 if (runno >= 2200 && runno <= 2299) {scale = 1.001722184;
return;}
1373 if (runno >= 700 && runno <= 999) {
1374 B2ERROR(
"scale_momenta not ready for this exp,run "
1375 << expno <<
"," << runno);
1377 if (runno >= 1014 && runno <= 1999) {
1378 B2ERROR(
"scale_momenta not ready for this exp,run "
1379 << expno <<
"," << runno);
1381 if (runno >= 2299) {
1382 B2ERROR(
"scale_momenta not ready for this exp,run "
1383 << expno <<
"," << runno);
1385 }
else if (expno == 73) {
1386 if (runno >= 0 && runno <= 99) {scale = 1.000721587;
return;}
1387 if (runno >= 100 && runno <= 199) {scale = 1.000707089;
return;}
1388 if (runno >= 200 && runno <= 299) {scale = 1.000722517;
return;}
1389 if (runno >= 300 && runno <= 399) {scale = 1.000722517;
return;}
1390 if (runno >= 400 && runno <= 499) {scale = 1.000750776;
return;}
1391 if (runno >= 500 && runno <= 599) {scale = 1.000729771;
return;}
1392 if (runno >= 600 && runno <= 699) {scale = 1.000751190;
return;}
1393 if (runno >= 700 && runno <= 799) {scale = 1.000702455;
return;}
1394 if (runno >= 800 && runno <= 899) {scale = 1.000771074;
return;}
1395 if (runno >= 900 && runno <= 999) {scale = 1.000868463;
return;}
1396 if (runno >= 1000) {
1397 B2ERROR(
"scale_momenta not ready for this exp,run "
1398 << expno <<
"," << runno);
1402 B2ERROR(
"scale_momenta mode=2 not ready for exp " << expno);
2780 int no_exp{0}, no_run{0}, no_evt{0}, no_frm{0}, expmc{0};
2781 get_event_id(&no_exp, &no_run, &no_evt, &no_frm, &expmc);
2786 static int first_exp73 =
true;
2787 if (no_exp == 73 && first_exp73) {
2788 first_exp73 =
false;
2790 "scale_error: warning: scale parameters for exp#71 are tentatively used for exp#73.");
2803 Belle::Mdst_trk_fit_Manager& fitmgr = Belle::Mdst_trk_fit_Manager::get_manager();
2804 for (std::vector<Belle::Mdst_trk_fit>::iterator it = fitmgr.begin(); it != fitmgr.end(); ++it) {
2805 Belle::Mdst_trk_fit& fit = *it;
2806 if (fit.helix(2) == 0.)
continue;
2808 const double pt = 1. / fabs(fit.helix(2));
2809 const double tanl = fit.helix(4);
2813 fit.error(0, scale[0]*scale[0]*fit.error(0));
2814 fit.error(1, scale[1]*scale[0]*fit.error(1));
2815 fit.error(2, scale[1]*scale[1]*fit.error(2));
2816 fit.error(3, scale[2]*scale[0]*fit.error(3));
2817 fit.error(4, scale[2]*scale[1]*fit.error(4));
2818 fit.error(5, scale[2]*scale[2]*fit.error(5));
2819 fit.error(6, scale[3]*scale[0]*fit.error(6));
2820 fit.error(7, scale[3]*scale[1]*fit.error(7));
2821 fit.error(8, scale[3]*scale[2]*fit.error(8));
2822 fit.error(9, scale[3]*scale[3]*fit.error(9));
2823 fit.error(10, scale[4]*scale[0]*fit.error(10));
2824 fit.error(11, scale[4]*scale[1]*fit.error(11));
2825 fit.error(12, scale[4]*scale[2]*fit.error(12));
2826 fit.error(13, scale[4]*scale[3]*fit.error(13));
2827 fit.error(14, scale[4]*scale[4]*fit.error(14));
2832 Belle::Mdst_vee_daughters_Manager& daumgr = Belle::Mdst_vee_daughters_Manager::get_manager();
2833 for (std::vector<Belle::Mdst_vee_daughters>::iterator it = daumgr.begin(); it != daumgr.end(); ++it) {
2834 Belle::Mdst_vee_daughters& dau = *it;
2835 if (dau.helix_p(2) == 0. || dau.helix_m(2) == 0.)
continue;
2838 const double pt_p = 1. / fabs(dau.helix_p(2));
2839 const double tanl_p = dau.helix_p(4);
2843 dau.error_p(0, scale[0]*scale[0]*dau.error_p(0));
2844 dau.error_p(1, scale[1]*scale[0]*dau.error_p(1));
2845 dau.error_p(2, scale[1]*scale[1]*dau.error_p(2));
2846 dau.error_p(3, scale[2]*scale[0]*dau.error_p(3));
2847 dau.error_p(4, scale[2]*scale[1]*dau.error_p(4));
2848 dau.error_p(5, scale[2]*scale[2]*dau.error_p(5));
2849 dau.error_p(6, scale[3]*scale[0]*dau.error_p(6));
2850 dau.error_p(7, scale[3]*scale[1]*dau.error_p(7));
2851 dau.error_p(8, scale[3]*scale[2]*dau.error_p(8));
2852 dau.error_p(9, scale[3]*scale[3]*dau.error_p(9));
2853 dau.error_p(10, scale[4]*scale[0]*dau.error_p(10));
2854 dau.error_p(11, scale[4]*scale[1]*dau.error_p(11));
2855 dau.error_p(12, scale[4]*scale[2]*dau.error_p(12));
2856 dau.error_p(13, scale[4]*scale[3]*dau.error_p(13));
2857 dau.error_p(14, scale[4]*scale[4]*dau.error_p(14));
2860 const double pt_m = 1. / fabs(dau.helix_m(2));
2861 const double tanl_m = dau.helix_m(4);
2865 dau.error_m(0, scale[0]*scale[0]*dau.error_m(0));
2866 dau.error_m(1, scale[1]*scale[0]*dau.error_m(1));
2867 dau.error_m(2, scale[1]*scale[1]*dau.error_m(2));
2868 dau.error_m(3, scale[2]*scale[0]*dau.error_m(3));
2869 dau.error_m(4, scale[2]*scale[1]*dau.error_m(4));
2870 dau.error_m(5, scale[2]*scale[2]*dau.error_m(5));
2871 dau.error_m(6, scale[3]*scale[0]*dau.error_m(6));
2872 dau.error_m(7, scale[3]*scale[1]*dau.error_m(7));
2873 dau.error_m(8, scale[3]*scale[2]*dau.error_m(8));
2874 dau.error_m(9, scale[3]*scale[3]*dau.error_m(9));
2875 dau.error_m(10, scale[4]*scale[0]*dau.error_m(10));
2876 dau.error_m(11, scale[4]*scale[1]*dau.error_m(11));
2877 dau.error_m(12, scale[4]*scale[2]*dau.error_m(12));
2878 dau.error_m(13, scale[4]*scale[3]*dau.error_m(13));
2879 dau.error_m(14, scale[4]*scale[4]*dau.error_m(14));