tstGamma.f90
./Manager/emit/tests/tstGamma.f90
202 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 5.8041311241133610E-004 , new=",new,")"
225 insist(sum(this%binned(:, 11 ))== 0.0000000000000000 ," 40070 has zero counts if not has_gamma")
227 insist(sum(this%binned(:, 12 ))== 0.0000000000000000 ," 40080 has zero counts if not has_gamma")
229 insist(sum(this%binned(:, 13 ))== 0.0000000000000000 ," 40090 has zero counts if not has_gamma")
236 insist(sum(this%binned(:, 15 ))== 0.0000000000000000 ," 40110 has zero counts if not has_gamma")
238 insist(sum(this%binned(:, 16 ))== 0.0000000000000000 ," 50100 has zero counts if not has_gamma")
240 insist(sum(this%binned(:, 17 ))== 0.0000000000000000 ," 50110 has zero counts if not has_gamma")
242 insist(sum(this%binned(:, 18 ))== 0.0000000000000000 ," 50120 has zero counts if not has_gamma")
244 insist(sum(this%binned(:, 19 ))== 0.0000000000000000 ," 60120 has zero counts if not has_gamma")
246 insist(sum(this%binned(:, 20 ))== 0.0000000000000000 ," 60130 has zero counts if not has_gamma")
250 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 3.0992144909162689E-002 , new=",new,")"
263 insist(sum(this%binned(:, 24 ))== 0.0000000000000000 ," 70140 has zero counts if not has_gamma")
265 insist(sum(this%binned(:, 25 ))== 0.0000000000000000 ," 70150 has zero counts if not has_gamma")
272 insist(sum(this%binned(:, 27 ))== 0.0000000000000000 ," 80160 has zero counts if not has_gamma")
274 insist(sum(this%binned(:, 28 ))== 0.0000000000000000 ," 80170 has zero counts if not has_gamma")
276 insist(sum(this%binned(:, 29 ))== 0.0000000000000000 ," 80180 has zero counts if not has_gamma")
278 insist(sum(this%binned(:, 30 ))== 0.0000000000000000 ," 80190 has zero counts if not has_gamma")
280 insist(sum(this%binned(:, 31 ))== 0.0000000000000000 ," 90190 has zero counts if not has_gamma")
282 insist(sum(this%binned(:, 32 ))== 0.0000000000000000 ," 90200 has zero counts if not has_gamma")
284 insist(sum(this%binned(:, 33 ))== 0.0000000000000000 ," 100200 has zero counts if not has_gamma")
286 insist(sum(this%binned(:, 34 ))== 0.0000000000000000 ," 100210 has zero counts if not has_gamma")
288 insist(sum(this%binned(:, 35 ))== 0.0000000000000000 ," 100220 has zero counts if not has_gamma")
300 insist(sum(this%binned(:, 38 ))== 0.0000000000000000 ," 110230 has zero counts if not has_gamma")
307 insist(sum(this%binned(:, 40 ))== 0.0000000000000000 ," 110241 has zero counts if not has_gamma")
314 insist(sum(this%binned(:, 42 ))== 0.0000000000000000 ," 120240 has zero counts if not has_gamma")
316 insist(sum(this%binned(:, 43 ))== 0.0000000000000000 ," 120250 has zero counts if not has_gamma")
318 insist(sum(this%binned(:, 44 ))== 0.0000000000000000 ," 120260 has zero counts if not has_gamma")
330 insist(sum(this%binned(:, 47 ))== 0.0000000000000000 ," 130260 has zero counts if not has_gamma")
332 insist(sum(this%binned(:, 48 ))== 0.0000000000000000 ," 130270 has zero counts if not has_gamma")
349 insist(sum(this%binned(:, 52 ))== 0.0000000000000000 ," 140280 has zero counts if not has_gamma")
351 insist(sum(this%binned(:, 53 ))== 0.0000000000000000 ," 140290 has zero counts if not has_gamma")
353 insist(sum(this%binned(:, 54 ))== 0.0000000000000000 ," 140300 has zero counts if not has_gamma")
362 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 4.5932783156572995E-002 , new=",new,")"
365 insist(sum(this%binned(:, 57 ))== 0.0000000000000000 ," 150310 has zero counts if not has_gamma")
374 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 5.7196099764841087E-002 , new=",new,")"
377 insist(sum(this%binned(:, 60 ))== 0.0000000000000000 ," 150340 has zero counts if not has_gamma")
379 insist(sum(this%binned(:, 61 ))== 0.0000000000000000 ," 160320 has zero counts if not has_gamma")
381 insist(sum(this%binned(:, 62 ))== 0.0000000000000000 ," 160330 has zero counts if not has_gamma")
383 insist(sum(this%binned(:, 63 ))== 0.0000000000000000 ," 160340 has zero counts if not has_gamma")
387 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 3.1893720850639802E-002 , new=",new,")"
390 insist(sum(this%binned(:, 65 ))== 0.0000000000000000 ," 160360 has zero counts if not has_gamma")
392 insist(sum(this%binned(:, 66 ))== 0.0000000000000000 ," 160370 has zero counts if not has_gamma")
394 insist(sum(this%binned(:, 67 ))== 0.0000000000000000 ," 170350 has zero counts if not has_gamma")
401 insist(sum(this%binned(:, 69 ))== 0.0000000000000000 ," 170370 has zero counts if not has_gamma")
408 insist(sum(this%binned(:, 71 ))== 0.0000000000000000 ," 170381 has zero counts if not has_gamma")
410 insist(sum(this%binned(:, 72 ))== 0.0000000000000000 ," 180360 has zero counts if not has_gamma")
412 insist(sum(this%binned(:, 73 ))== 0.0000000000000000 ," 180370 has zero counts if not has_gamma")
414 insist(sum(this%binned(:, 74 ))== 0.0000000000000000 ," 180380 has zero counts if not has_gamma")
416 insist(sum(this%binned(:, 75 ))== 0.0000000000000000 ," 180390 has zero counts if not has_gamma")
418 insist(sum(this%binned(:, 76 ))== 0.0000000000000000 ," 180400 has zero counts if not has_gamma")
425 insist(sum(this%binned(:, 78 ))== 0.0000000000000000 ," 180420 has zero counts if not has_gamma")
427 insist(sum(this%binned(:, 79 ))== 0.0000000000000000 ," 190390 has zero counts if not has_gamma")
434 insist(sum(this%binned(:, 81 ))== 0.0000000000000000 ," 190410 has zero counts if not has_gamma")
446 insist(sum(this%binned(:, 84 ))== 0.0000000000000000 ," 190440 has zero counts if not has_gamma")
448 insist(sum(this%binned(:, 85 ))== 0.0000000000000000 ," 200400 has zero counts if not has_gamma")
450 insist(sum(this%binned(:, 86 ))== 0.0000000000000000 ," 200410 has zero counts if not has_gamma")
452 insist(sum(this%binned(:, 87 ))== 0.0000000000000000 ," 200420 has zero counts if not has_gamma")
454 insist(sum(this%binned(:, 88 ))== 0.0000000000000000 ," 200430 has zero counts if not has_gamma")
456 insist(sum(this%binned(:, 89 ))== 0.0000000000000000 ," 200440 has zero counts if not has_gamma")
460 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 5.8461036379483737E-002 , new=",new,")"
463 insist(sum(this%binned(:, 91 ))== 0.0000000000000000 ," 200460 has zero counts if not has_gamma")
470 insist(sum(this%binned(:, 93 ))== 0.0000000000000000 ," 200480 has zero counts if not has_gamma")
472 insist(sum(this%binned(:, 94 ))== 0.0000000000000000 ," 200490 has zero counts if not has_gamma")
474 insist(sum(this%binned(:, 95 ))== 0.0000000000000000 ," 210440 has zero counts if not has_gamma")
476 insist(sum(this%binned(:, 96 ))== 0.0000000000000000 ," 210441 has zero counts if not has_gamma")
478 insist(sum(this%binned(:, 97 ))== 0.0000000000000000 ," 210450 has zero counts if not has_gamma")
480 insist(sum(this%binned(:, 98 ))== 0.0000000000000000 ," 210451 has zero counts if not has_gamma")
484 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 9.2907327068596016E-002 , new=",new,")"
487 insist(sum(this%binned(:, 100 ))== 0.0000000000000000 ," 210461 has zero counts if not has_gamma")
499 insist(sum(this%binned(:, 103 ))== 0.0000000000000000 ," 210490 has zero counts if not has_gamma")
506 insist(sum(this%binned(:, 105 ))== 0.0000000000000000 ," 220440 has zero counts if not has_gamma")
508 insist(sum(this%binned(:, 106 ))== 0.0000000000000000 ," 220450 has zero counts if not has_gamma")
510 insist(sum(this%binned(:, 107 ))== 0.0000000000000000 ," 220460 has zero counts if not has_gamma")
512 insist(sum(this%binned(:, 108 ))== 0.0000000000000000 ," 220470 has zero counts if not has_gamma")
514 insist(sum(this%binned(:, 109 ))== 0.0000000000000000 ," 220480 has zero counts if not has_gamma")
516 insist(sum(this%binned(:, 110 ))== 0.0000000000000000 ," 220490 has zero counts if not has_gamma")
518 insist(sum(this%binned(:, 111 ))== 0.0000000000000000 ," 220500 has zero counts if not has_gamma")
525 insist(sum(this%binned(:, 113 ))== 0.0000000000000000 ," 230480 has zero counts if not has_gamma")
527 insist(sum(this%binned(:, 114 ))== 0.0000000000000000 ," 230490 has zero counts if not has_gamma")
529 insist(sum(this%binned(:, 115 ))== 0.0000000000000000 ," 230500 has zero counts if not has_gamma")
531 insist(sum(this%binned(:, 116 ))== 0.0000000000000000 ," 230510 has zero counts if not has_gamma")
543 insist(sum(this%binned(:, 119 ))== 0.0000000000000000 ," 230540 has zero counts if not has_gamma")
545 insist(sum(this%binned(:, 120 ))== 0.0000000000000000 ," 240480 has zero counts if not has_gamma")
547 insist(sum(this%binned(:, 121 ))== 0.0000000000000000 ," 240490 has zero counts if not has_gamma")
549 insist(sum(this%binned(:, 122 ))== 0.0000000000000000 ," 240500 has zero counts if not has_gamma")
551 insist(sum(this%binned(:, 123 ))== 0.0000000000000000 ," 240510 has zero counts if not has_gamma")
553 insist(sum(this%binned(:, 124 ))== 0.0000000000000000 ," 240520 has zero counts if not has_gamma")
555 insist(sum(this%binned(:, 125 ))== 0.0000000000000000 ," 240530 has zero counts if not has_gamma")
557 insist(sum(this%binned(:, 126 ))== 0.0000000000000000 ," 240540 has zero counts if not has_gamma")
564 insist(sum(this%binned(:, 128 ))== 0.0000000000000000 ," 250520 has zero counts if not has_gamma")
566 insist(sum(this%binned(:, 129 ))== 0.0000000000000000 ," 250530 has zero counts if not has_gamma")
568 insist(sum(this%binned(:, 130 ))== 0.0000000000000000 ," 250540 has zero counts if not has_gamma")
570 insist(sum(this%binned(:, 131 ))== 0.0000000000000000 ," 250550 has zero counts if not has_gamma")
587 insist(sum(this%binned(:, 135 ))== 0.0000000000000000 ," 260540 has zero counts if not has_gamma")
589 insist(sum(this%binned(:, 136 ))== 0.0000000000000000 ," 260550 has zero counts if not has_gamma")
591 insist(sum(this%binned(:, 137 ))== 0.0000000000000000 ," 260560 has zero counts if not has_gamma")
593 insist(sum(this%binned(:, 138 ))== 0.0000000000000000 ," 260570 has zero counts if not has_gamma")
595 insist(sum(this%binned(:, 139 ))== 0.0000000000000000 ," 260580 has zero counts if not has_gamma")
599 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 9.9839981208801470E-002 , new=",new,")"
602 insist(sum(this%binned(:, 141 ))== 0.0000000000000000 ," 260600 has zero counts if not has_gamma")
604 insist(sum(this%binned(:, 142 ))== 0.0000000000000000 ," 270550 has zero counts if not has_gamma")
606 insist(sum(this%binned(:, 143 ))== 0.0000000000000000 ," 270560 has zero counts if not has_gamma")
608 insist(sum(this%binned(:, 144 ))== 0.0000000000000000 ," 270570 has zero counts if not has_gamma")
610 insist(sum(this%binned(:, 145 ))== 0.0000000000000000 ," 270581 has zero counts if not has_gamma")
617 insist(sum(this%binned(:, 147 ))== 0.0000000000000000 ," 270590 has zero counts if not has_gamma")
621 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 7.7794668172793124E-002 , new=",new,")"
626 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 4.0713854684785566E-006 , new=",new,")"
639 insist(sum(this%binned(:, 152 ))== 0.0000000000000000 ," 280560 has zero counts if not has_gamma")
641 insist(sum(this%binned(:, 153 ))== 0.0000000000000000 ," 280570 has zero counts if not has_gamma")
643 insist(sum(this%binned(:, 154 ))== 0.0000000000000000 ," 280580 has zero counts if not has_gamma")
645 insist(sum(this%binned(:, 155 ))== 0.0000000000000000 ," 280590 has zero counts if not has_gamma")
647 insist(sum(this%binned(:, 156 ))== 0.0000000000000000 ," 280600 has zero counts if not has_gamma")
649 insist(sum(this%binned(:, 157 ))== 0.0000000000000000 ," 280610 has zero counts if not has_gamma")
651 insist(sum(this%binned(:, 158 ))== 0.0000000000000000 ," 280620 has zero counts if not has_gamma")
655 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 6.8060587681108660E-003 , new=",new,")"
658 insist(sum(this%binned(:, 160 ))== 0.0000000000000000 ," 280640 has zero counts if not has_gamma")
667 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 4.9393710887042759E-002 , new=",new,")"
675 insist(sum(this%binned(:, 164 ))== 0.0000000000000000 ," 290630 has zero counts if not has_gamma")
682 insist(sum(this%binned(:, 166 ))== 0.0000000000000000 ," 290650 has zero counts if not has_gamma")
699 insist(sum(this%binned(:, 170 ))== 0.0000000000000000 ," 300640 has zero counts if not has_gamma")
703 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.9869193073204290E-003 , new=",new,")"
706 insist(sum(this%binned(:, 172 ))== 0.0000000000000000 ," 300660 has zero counts if not has_gamma")
708 insist(sum(this%binned(:, 173 ))== 0.0000000000000000 ," 300670 has zero counts if not has_gamma")
710 insist(sum(this%binned(:, 174 ))== 0.0000000000000000 ," 300680 has zero counts if not has_gamma")
719 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 3.1907135371456243E-008 , new=",new,")"
722 insist(sum(this%binned(:, 177 ))== 0.0000000000000000 ," 300700 has zero counts if not has_gamma")
736 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 6.7099930788272200E-002 , new=",new,")"
739 insist(sum(this%binned(:, 181 ))== 0.0000000000000000 ," 310670 has zero counts if not has_gamma")
741 insist(sum(this%binned(:, 182 ))== 0.0000000000000000 ," 310680 has zero counts if not has_gamma")
743 insist(sum(this%binned(:, 183 ))== 0.0000000000000000 ," 310690 has zero counts if not has_gamma")
750 insist(sum(this%binned(:, 185 ))== 0.0000000000000000 ," 310710 has zero counts if not has_gamma")
757 insist(sum(this%binned(:, 187 ))== 0.0000000000000000 ," 310721 has zero counts if not has_gamma")
759 insist(sum(this%binned(:, 188 ))== 0.0000000000000000 ," 320680 has zero counts if not has_gamma")
761 insist(sum(this%binned(:, 189 ))== 0.0000000000000000 ," 320690 has zero counts if not has_gamma")
763 insist(sum(this%binned(:, 190 ))== 0.0000000000000000 ," 320700 has zero counts if not has_gamma")
765 insist(sum(this%binned(:, 191 ))== 0.0000000000000000 ," 320710 has zero counts if not has_gamma")
767 insist(sum(this%binned(:, 192 ))== 0.0000000000000000 ," 320711 has zero counts if not has_gamma")
769 insist(sum(this%binned(:, 193 ))== 0.0000000000000000 ," 320720 has zero counts if not has_gamma")
771 insist(sum(this%binned(:, 194 ))== 0.0000000000000000 ," 320730 has zero counts if not has_gamma")
773 insist(sum(this%binned(:, 195 ))== 0.0000000000000000 ," 320731 has zero counts if not has_gamma")
775 insist(sum(this%binned(:, 196 ))== 0.0000000000000000 ," 320740 has zero counts if not has_gamma")
782 insist(sum(this%binned(:, 198 ))== 0.0000000000000000 ," 320751 has zero counts if not has_gamma")
784 insist(sum(this%binned(:, 199 ))== 0.0000000000000000 ," 320760 has zero counts if not has_gamma")
796 insist(sum(this%binned(:, 202 ))== 0.0000000000000000 ," 330710 has zero counts if not has_gamma")
798 insist(sum(this%binned(:, 203 ))== 0.0000000000000000 ," 330720 has zero counts if not has_gamma")
800 insist(sum(this%binned(:, 204 ))== 0.0000000000000000 ," 330730 has zero counts if not has_gamma")
802 insist(sum(this%binned(:, 205 ))== 0.0000000000000000 ," 330740 has zero counts if not has_gamma")
804 insist(sum(this%binned(:, 206 ))== 0.0000000000000000 ," 330750 has zero counts if not has_gamma")
806 insist(sum(this%binned(:, 207 ))== 0.0000000000000000 ," 330751 has zero counts if not has_gamma")
818 insist(sum(this%binned(:, 210 ))== 0.0000000000000000 ," 340720 has zero counts if not has_gamma")
820 insist(sum(this%binned(:, 211 ))== 0.0000000000000000 ," 340730 has zero counts if not has_gamma")
822 insist(sum(this%binned(:, 212 ))== 0.0000000000000000 ," 340740 has zero counts if not has_gamma")
824 insist(sum(this%binned(:, 213 ))== 0.0000000000000000 ," 340750 has zero counts if not has_gamma")
826 insist(sum(this%binned(:, 214 ))== 0.0000000000000000 ," 340760 has zero counts if not has_gamma")
828 insist(sum(this%binned(:, 215 ))== 0.0000000000000000 ," 340770 has zero counts if not has_gamma")
830 insist(sum(this%binned(:, 216 ))== 0.0000000000000000 ," 340771 has zero counts if not has_gamma")
832 insist(sum(this%binned(:, 217 ))== 0.0000000000000000 ," 340780 has zero counts if not has_gamma")
836 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 3.6759974434876781E-002 , new=",new,")"
839 insist(sum(this%binned(:, 219 ))== 0.0000000000000000 ," 340791 has zero counts if not has_gamma")
841 insist(sum(this%binned(:, 220 ))== 0.0000000000000000 ," 340800 has zero counts if not has_gamma")
843 insist(sum(this%binned(:, 221 ))== 0.0000000000000000 ," 340810 has zero counts if not has_gamma")
845 insist(sum(this%binned(:, 222 ))== 0.0000000000000000 ," 340811 has zero counts if not has_gamma")
847 insist(sum(this%binned(:, 223 ))== 0.0000000000000000 ," 340820 has zero counts if not has_gamma")
859 insist(sum(this%binned(:, 226 ))== 0.0000000000000000 ," 350760 has zero counts if not has_gamma")
861 insist(sum(this%binned(:, 227 ))== 0.0000000000000000 ," 350770 has zero counts if not has_gamma")
863 insist(sum(this%binned(:, 228 ))== 0.0000000000000000 ," 350771 has zero counts if not has_gamma")
865 insist(sum(this%binned(:, 229 ))== 0.0000000000000000 ," 350780 has zero counts if not has_gamma")
867 insist(sum(this%binned(:, 230 ))== 0.0000000000000000 ," 350790 has zero counts if not has_gamma")
874 insist(sum(this%binned(:, 232 ))== 0.0000000000000000 ," 350801 has zero counts if not has_gamma")
876 insist(sum(this%binned(:, 233 ))== 0.0000000000000000 ," 350810 has zero counts if not has_gamma")
885 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 9.4313711645749483E-004 , new=",new,")"
893 insist(sum(this%binned(:, 237 ))== 0.0000000000000000 ," 360760 has zero counts if not has_gamma")
895 insist(sum(this%binned(:, 238 ))== 0.0000000000000000 ," 360770 has zero counts if not has_gamma")
897 insist(sum(this%binned(:, 239 ))== 0.0000000000000000 ," 360780 has zero counts if not has_gamma")
901 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 5.5476567838820726E-002 , new=",new,")"
904 insist(sum(this%binned(:, 241 ))== 0.0000000000000000 ," 360791 has zero counts if not has_gamma")
906 insist(sum(this%binned(:, 242 ))== 0.0000000000000000 ," 360800 has zero counts if not has_gamma")
908 insist(sum(this%binned(:, 243 ))== 0.0000000000000000 ," 360810 has zero counts if not has_gamma")
910 insist(sum(this%binned(:, 244 ))== 0.0000000000000000 ," 360811 has zero counts if not has_gamma")
912 insist(sum(this%binned(:, 245 ))== 0.0000000000000000 ," 360820 has zero counts if not has_gamma")
914 insist(sum(this%binned(:, 246 ))== 0.0000000000000000 ," 360830 has zero counts if not has_gamma")
916 insist(sum(this%binned(:, 247 ))== 0.0000000000000000 ," 360831 has zero counts if not has_gamma")
918 insist(sum(this%binned(:, 248 ))== 0.0000000000000000 ," 360840 has zero counts if not has_gamma")
930 insist(sum(this%binned(:, 251 ))== 0.0000000000000000 ," 360860 has zero counts if not has_gamma")
942 insist(sum(this%binned(:, 254 ))== 0.0000000000000000 ," 370820 has zero counts if not has_gamma")
944 insist(sum(this%binned(:, 255 ))== 0.0000000000000000 ," 370830 has zero counts if not has_gamma")
946 insist(sum(this%binned(:, 256 ))== 0.0000000000000000 ," 370840 has zero counts if not has_gamma")
948 insist(sum(this%binned(:, 257 ))== 0.0000000000000000 ," 370850 has zero counts if not has_gamma")
955 insist(sum(this%binned(:, 259 ))== 0.0000000000000000 ," 370861 has zero counts if not has_gamma")
959 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 6.0877798994170014E-002 , new=",new,")"
972 insist(sum(this%binned(:, 263 ))== 0.0000000000000000 ," 380820 has zero counts if not has_gamma")
974 insist(sum(this%binned(:, 264 ))== 0.0000000000000000 ," 380830 has zero counts if not has_gamma")
976 insist(sum(this%binned(:, 265 ))== 0.0000000000000000 ," 380840 has zero counts if not has_gamma")
978 insist(sum(this%binned(:, 266 ))== 0.0000000000000000 ," 380850 has zero counts if not has_gamma")
980 insist(sum(this%binned(:, 267 ))== 0.0000000000000000 ," 380851 has zero counts if not has_gamma")
982 insist(sum(this%binned(:, 268 ))== 0.0000000000000000 ," 380860 has zero counts if not has_gamma")
984 insist(sum(this%binned(:, 269 ))== 0.0000000000000000 ," 380870 has zero counts if not has_gamma")
986 insist(sum(this%binned(:, 270 ))== 0.0000000000000000 ," 380871 has zero counts if not has_gamma")
988 insist(sum(this%binned(:, 271 ))== 0.0000000000000000 ," 380880 has zero counts if not has_gamma")
1002 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.75223841715706297 , new=",new,")"
1005 insist(sum(this%binned(:, 275 ))== 0.0000000000000000 ," 380930 has zero counts if not has_gamma")
1007 insist(sum(this%binned(:, 276 ))== 0.0000000000000000 ," 390860 has zero counts if not has_gamma")
1009 insist(sum(this%binned(:, 277 ))== 0.0000000000000000 ," 390870 has zero counts if not has_gamma")
1011 insist(sum(this%binned(:, 278 ))== 0.0000000000000000 ," 390871 has zero counts if not has_gamma")
1013 insist(sum(this%binned(:, 279 ))== 0.0000000000000000 ," 390880 has zero counts if not has_gamma")
1015 insist(sum(this%binned(:, 280 ))== 0.0000000000000000 ," 390890 has zero counts if not has_gamma")
1017 insist(sum(this%binned(:, 281 ))== 0.0000000000000000 ," 390891 has zero counts if not has_gamma")
1026 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.4444724824302237E-009 , new=",new,")"
1031 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.68106379883573032 , new=",new,")"
1034 insist(sum(this%binned(:, 285 ))== 0.0000000000000000 ," 390911 has zero counts if not has_gamma")
1046 insist(sum(this%binned(:, 288 ))== 0.0000000000000000 ," 390931 has zero counts if not has_gamma")
1058 insist(sum(this%binned(:, 291 ))== 0.0000000000000000 ," 400860 has zero counts if not has_gamma")
1060 insist(sum(this%binned(:, 292 ))== 0.0000000000000000 ," 400870 has zero counts if not has_gamma")
1062 insist(sum(this%binned(:, 293 ))== 0.0000000000000000 ," 400880 has zero counts if not has_gamma")
1066 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.21884140142702282 , new=",new,")"
1069 insist(sum(this%binned(:, 295 ))== 0.0000000000000000 ," 400900 has zero counts if not has_gamma")
1071 insist(sum(this%binned(:, 296 ))== 0.0000000000000000 ," 400901 has zero counts if not has_gamma")
1073 insist(sum(this%binned(:, 297 ))== 0.0000000000000000 ," 400910 has zero counts if not has_gamma")
1075 insist(sum(this%binned(:, 298 ))== 0.0000000000000000 ," 400920 has zero counts if not has_gamma")
1079 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.0673408209577246E-003 , new=",new,")"
1082 insist(sum(this%binned(:, 300 ))== 0.0000000000000000 ," 400940 has zero counts if not has_gamma")
1086 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 9.8333829532349815E-002 , new=",new,")"
1089 insist(sum(this%binned(:, 302 ))== 0.0000000000000000 ," 400960 has zero counts if not has_gamma")
1093 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.79223635100248524 , new=",new,")"
1096 insist(sum(this%binned(:, 304 ))== 0.0000000000000000 ," 410900 has zero counts if not has_gamma")
1098 insist(sum(this%binned(:, 305 ))== 0.0000000000000000 ," 410901 has zero counts if not has_gamma")
1100 insist(sum(this%binned(:, 306 ))== 0.0000000000000000 ," 410910 has zero counts if not has_gamma")
1102 insist(sum(this%binned(:, 307 ))== 0.0000000000000000 ," 410911 has zero counts if not has_gamma")
1104 insist(sum(this%binned(:, 308 ))== 0.0000000000000000 ," 410920 has zero counts if not has_gamma")
1106 insist(sum(this%binned(:, 309 ))== 0.0000000000000000 ," 410921 has zero counts if not has_gamma")
1108 insist(sum(this%binned(:, 310 ))== 0.0000000000000000 ," 410930 has zero counts if not has_gamma")
1110 insist(sum(this%binned(:, 311 ))== 0.0000000000000000 ," 410931 has zero counts if not has_gamma")
1114 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.12969260956754219 , new=",new,")"
1119 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.0912115440728130E-005 , new=",new,")"
1124 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 2.7653486272365409E-002 , new=",new,")"
1127 insist(sum(this%binned(:, 315 ))== 0.0000000000000000 ," 410951 has zero counts if not has_gamma")
1131 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.24353126187349969 , new=",new,")"
1136 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.50117733078528348 , new=",new,")"
1139 insist(sum(this%binned(:, 318 ))== 0.0000000000000000 ," 410971 has zero counts if not has_gamma")
1146 insist(sum(this%binned(:, 320 ))== 0.0000000000000000 ," 411000 has zero counts if not has_gamma")
1148 insist(sum(this%binned(:, 321 ))== 0.0000000000000000 ," 420920 has zero counts if not has_gamma")
1150 insist(sum(this%binned(:, 322 ))== 0.0000000000000000 ," 420931 has zero counts if not has_gamma")
1152 insist(sum(this%binned(:, 323 ))== 0.0000000000000000 ," 420930 has zero counts if not has_gamma")
1154 insist(sum(this%binned(:, 324 ))== 0.0000000000000000 ," 420940 has zero counts if not has_gamma")
1156 insist(sum(this%binned(:, 325 ))== 0.0000000000000000 ," 420950 has zero counts if not has_gamma")
1158 insist(sum(this%binned(:, 326 ))== 0.0000000000000000 ," 420960 has zero counts if not has_gamma")
1160 insist(sum(this%binned(:, 327 ))== 0.0000000000000000 ," 420970 has zero counts if not has_gamma")
1162 insist(sum(this%binned(:, 328 ))== 0.0000000000000000 ," 420980 has zero counts if not has_gamma")
1166 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.40968073719541048 , new=",new,")"
1169 insist(sum(this%binned(:, 330 ))== 0.0000000000000000 ," 421000 has zero counts if not has_gamma")
1173 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.63661447291623863 , new=",new,")"
1176 insist(sum(this%binned(:, 332 ))== 0.0000000000000000 ," 430950 has zero counts if not has_gamma")
1178 insist(sum(this%binned(:, 333 ))== 0.0000000000000000 ," 430951 has zero counts if not has_gamma")
1180 insist(sum(this%binned(:, 334 ))== 0.0000000000000000 ," 430960 has zero counts if not has_gamma")
1182 insist(sum(this%binned(:, 335 ))== 0.0000000000000000 ," 430970 has zero counts if not has_gamma")
1184 insist(sum(this%binned(:, 336 ))== 0.0000000000000000 ," 430971 has zero counts if not has_gamma")
1188 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.10194629045810190 , new=",new,")"
1193 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 6.6673559170234667E-002 , new=",new,")"
1196 insist(sum(this%binned(:, 339 ))== 0.0000000000000000 ," 430991 has zero counts if not has_gamma")
1205 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.50693567226374714 , new=",new,")"
1208 insist(sum(this%binned(:, 342 ))== 0.0000000000000000 ," 440960 has zero counts if not has_gamma")
1210 insist(sum(this%binned(:, 343 ))== 0.0000000000000000 ," 440970 has zero counts if not has_gamma")
1212 insist(sum(this%binned(:, 344 ))== 0.0000000000000000 ," 440980 has zero counts if not has_gamma")
1214 insist(sum(this%binned(:, 345 ))== 0.0000000000000000 ," 440990 has zero counts if not has_gamma")
1216 insist(sum(this%binned(:, 346 ))== 0.0000000000000000 ," 441000 has zero counts if not has_gamma")
1218 insist(sum(this%binned(:, 347 ))== 0.0000000000000000 ," 441010 has zero counts if not has_gamma")
1220 insist(sum(this%binned(:, 348 ))== 0.0000000000000000 ," 441020 has zero counts if not has_gamma")
1224 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 5.6067314364812035E-002 , new=",new,")"
1227 insist(sum(this%binned(:, 350 ))== 0.0000000000000000 ," 441040 has zero counts if not has_gamma")
1231 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.41997130747915101 , new=",new,")"
1236 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 2.5935015206863810E-003 , new=",new,")"
1244 insist(sum(this%binned(:, 354 ))== 0.0000000000000000 ," 450990 has zero counts if not has_gamma")
1246 insist(sum(this%binned(:, 355 ))== 0.0000000000000000 ," 450991 has zero counts if not has_gamma")
1248 insist(sum(this%binned(:, 356 ))== 0.0000000000000000 ," 451000 has zero counts if not has_gamma")
1250 insist(sum(this%binned(:, 357 ))== 0.0000000000000000 ," 451010 has zero counts if not has_gamma")
1252 insist(sum(this%binned(:, 358 ))== 0.0000000000000000 ," 451011 has zero counts if not has_gamma")
1254 insist(sum(this%binned(:, 359 ))== 0.0000000000000000 ," 451020 has zero counts if not has_gamma")
1256 insist(sum(this%binned(:, 360 ))== 0.0000000000000000 ," 451021 has zero counts if not has_gamma")
1258 insist(sum(this%binned(:, 361 ))== 0.0000000000000000 ," 451030 has zero counts if not has_gamma")
1260 insist(sum(this%binned(:, 362 ))== 0.0000000000000000 ," 451031 has zero counts if not has_gamma")
1269 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 5.0072926284761812E-007 , new=",new,")"
1274 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.13824386501530891 , new=",new,")"
1277 insist(sum(this%binned(:, 366 ))== 0.0000000000000000 ," 451051 has zero counts if not has_gamma")
1286 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.30463977297928585 , new=",new,")"
1291 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.49020845414431330 , new=",new,")"
1294 insist(sum(this%binned(:, 370 ))== 0.0000000000000000 ," 461000 has zero counts if not has_gamma")
1296 insist(sum(this%binned(:, 371 ))== 0.0000000000000000 ," 461010 has zero counts if not has_gamma")
1298 insist(sum(this%binned(:, 372 ))== 0.0000000000000000 ," 461020 has zero counts if not has_gamma")
1300 insist(sum(this%binned(:, 373 ))== 0.0000000000000000 ," 461030 has zero counts if not has_gamma")
1302 insist(sum(this%binned(:, 374 ))== 0.0000000000000000 ," 461040 has zero counts if not has_gamma")
1304 insist(sum(this%binned(:, 375 ))== 0.0000000000000000 ," 461050 has zero counts if not has_gamma")
1306 insist(sum(this%binned(:, 376 ))== 0.0000000000000000 ," 461060 has zero counts if not has_gamma")
1310 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 2.1676466463853572E-003 , new=",new,")"
1313 insist(sum(this%binned(:, 378 ))== 0.0000000000000000 ," 461071 has zero counts if not has_gamma")
1315 insist(sum(this%binned(:, 379 ))== 0.0000000000000000 ," 461080 has zero counts if not has_gamma")
1319 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.37413798787321895 , new=",new,")"
1322 insist(sum(this%binned(:, 381 ))== 0.0000000000000000 ," 461091 has zero counts if not has_gamma")
1324 insist(sum(this%binned(:, 382 ))== 0.0000000000000000 ," 461100 has zero counts if not has_gamma")
1326 insist(sum(this%binned(:, 383 ))== 0.0000000000000000 ," 461110 has zero counts if not has_gamma")
1328 insist(sum(this%binned(:, 384 ))== 0.0000000000000000 ," 461111 has zero counts if not has_gamma")
1330 insist(sum(this%binned(:, 385 ))== 0.0000000000000000 ," 461120 has zero counts if not has_gamma")
1332 insist(sum(this%binned(:, 386 ))== 0.0000000000000000 ," 471050 has zero counts if not has_gamma")
1334 insist(sum(this%binned(:, 387 ))== 0.0000000000000000 ," 471060 has zero counts if not has_gamma")
1336 insist(sum(this%binned(:, 388 ))== 0.0000000000000000 ," 471061 has zero counts if not has_gamma")
1338 insist(sum(this%binned(:, 389 ))== 0.0000000000000000 ," 471070 has zero counts if not has_gamma")
1340 insist(sum(this%binned(:, 390 ))== 0.0000000000000000 ," 471071 has zero counts if not has_gamma")
1344 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.66592332489877126 , new=",new,")"
1347 insist(sum(this%binned(:, 392 ))== 0.0000000000000000 ," 471081 has zero counts if not has_gamma")
1349 insist(sum(this%binned(:, 393 ))== 0.0000000000000000 ," 471090 has zero counts if not has_gamma")
1351 insist(sum(this%binned(:, 394 ))== 0.0000000000000000 ," 471091 has zero counts if not has_gamma")
1360 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 5.2684700660961419E-002 , new=",new,")"
1365 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.36628860078762671 , new=",new,")"
1368 insist(sum(this%binned(:, 398 ))== 0.0000000000000000 ," 471111 has zero counts if not has_gamma")
1370 insist(sum(this%binned(:, 399 ))== 0.0000000000000000 ," 471120 has zero counts if not has_gamma")
1372 insist(sum(this%binned(:, 400 ))== 0.0000000000000000 ," 481060 has zero counts if not has_gamma")
1376 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.3637423944359017E-003 , new=",new,")"
1379 insist(sum(this%binned(:, 402 ))== 0.0000000000000000 ," 481080 has zero counts if not has_gamma")
1381 insist(sum(this%binned(:, 403 ))== 0.0000000000000000 ," 481090 has zero counts if not has_gamma")
1383 insist(sum(this%binned(:, 404 ))== 0.0000000000000000 ," 481100 has zero counts if not has_gamma")
1385 insist(sum(this%binned(:, 405 ))== 0.0000000000000000 ," 481110 has zero counts if not has_gamma")
1387 insist(sum(this%binned(:, 406 ))== 0.0000000000000000 ," 481111 has zero counts if not has_gamma")
1389 insist(sum(this%binned(:, 407 ))== 0.0000000000000000 ," 481120 has zero counts if not has_gamma")
1393 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 7.5475555081041520E-002 , new=",new,")"
1398 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.17320447776984155 , new=",new,")"
1401 insist(sum(this%binned(:, 410 ))== 0.0000000000000000 ," 481140 has zero counts if not has_gamma")
1405 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.30763816677283629 , new=",new,")"
1410 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.67684358372750064 , new=",new,")"
1413 insist(sum(this%binned(:, 413 ))== 0.0000000000000000 ," 481160 has zero counts if not has_gamma")
1415 insist(sum(this%binned(:, 414 ))== 0.0000000000000000 ," 481170 has zero counts if not has_gamma")
1417 insist(sum(this%binned(:, 415 ))== 0.0000000000000000 ," 481171 has zero counts if not has_gamma")
1419 insist(sum(this%binned(:, 416 ))== 0.0000000000000000 ," 481190 has zero counts if not has_gamma")
1421 insist(sum(this%binned(:, 417 ))== 0.0000000000000000 ," 481210 has zero counts if not has_gamma")
1423 insist(sum(this%binned(:, 418 ))== 0.0000000000000000 ," 491110 has zero counts if not has_gamma")
1425 insist(sum(this%binned(:, 419 ))== 0.0000000000000000 ," 491120 has zero counts if not has_gamma")
1427 insist(sum(this%binned(:, 420 ))== 0.0000000000000000 ," 491130 has zero counts if not has_gamma")
1429 insist(sum(this%binned(:, 421 ))== 0.0000000000000000 ," 491131 has zero counts if not has_gamma")
1433 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.88342222361784639 , new=",new,")"
1436 insist(sum(this%binned(:, 423 ))== 0.0000000000000000 ," 491141 has zero counts if not has_gamma")
1440 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.13679583884205912 , new=",new,")"
1445 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 4.0115862475793927E-004 , new=",new,")"
1455 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.30791418396150139 , new=",new,")"
1458 insist(sum(this%binned(:, 428 ))== 0.0000000000000000 ," 491170 has zero counts if not has_gamma")
1460 insist(sum(this%binned(:, 429 ))== 0.0000000000000000 ," 491171 has zero counts if not has_gamma")
1467 insist(sum(this%binned(:, 431 ))== 0.0000000000000000 ," 491190 has zero counts if not has_gamma")
1469 insist(sum(this%binned(:, 432 ))== 0.0000000000000000 ," 491191 has zero counts if not has_gamma")
1481 insist(sum(this%binned(:, 435 ))== 0.0000000000000000 ," 491210 has zero counts if not has_gamma")
1483 insist(sum(this%binned(:, 436 ))== 0.0000000000000000 ," 491211 has zero counts if not has_gamma")
1485 insist(sum(this%binned(:, 437 ))== 0.0000000000000000 ," 501120 has zero counts if not has_gamma")
1487 insist(sum(this%binned(:, 438 ))== 0.0000000000000000 ," 501130 has zero counts if not has_gamma")
1489 insist(sum(this%binned(:, 439 ))== 0.0000000000000000 ," 501131 has zero counts if not has_gamma")
1491 insist(sum(this%binned(:, 440 ))== 0.0000000000000000 ," 501140 has zero counts if not has_gamma")
1493 insist(sum(this%binned(:, 441 ))== 0.0000000000000000 ," 501150 has zero counts if not has_gamma")
1495 insist(sum(this%binned(:, 442 ))== 0.0000000000000000 ," 501160 has zero counts if not has_gamma")
1497 insist(sum(this%binned(:, 443 ))== 0.0000000000000000 ," 501170 has zero counts if not has_gamma")
1499 insist(sum(this%binned(:, 444 ))== 0.0000000000000000 ," 501171 has zero counts if not has_gamma")
1501 insist(sum(this%binned(:, 445 ))== 0.0000000000000000 ," 501180 has zero counts if not has_gamma")
1503 insist(sum(this%binned(:, 446 ))== 0.0000000000000000 ," 501190 has zero counts if not has_gamma")
1505 insist(sum(this%binned(:, 447 ))== 0.0000000000000000 ," 501191 has zero counts if not has_gamma")
1507 insist(sum(this%binned(:, 448 ))== 0.0000000000000000 ," 501200 has zero counts if not has_gamma")
1509 insist(sum(this%binned(:, 449 ))== 0.0000000000000000 ," 501210 has zero counts if not has_gamma")
1511 insist(sum(this%binned(:, 450 ))== 0.0000000000000000 ," 501211 has zero counts if not has_gamma")
1513 insist(sum(this%binned(:, 451 ))== 0.0000000000000000 ," 501220 has zero counts if not has_gamma")
1517 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.57978316783614969 , new=",new,")"
1522 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.48882191691704341 , new=",new,")"
1525 insist(sum(this%binned(:, 454 ))== 0.0000000000000000 ," 501240 has zero counts if not has_gamma")
1529 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.95998007878166880 , new=",new,")"
1534 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.90123469909512754 , new=",new,")"
1539 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.4624188393718172E-002 , new=",new,")"
1542 insist(sum(this%binned(:, 458 ))== 0.0000000000000000 ," 511180 has zero counts if not has_gamma")
1544 insist(sum(this%binned(:, 459 ))== 0.0000000000000000 ," 511190 has zero counts if not has_gamma")
1546 insist(sum(this%binned(:, 460 ))== 0.0000000000000000 ," 511200 has zero counts if not has_gamma")
1548 insist(sum(this%binned(:, 461 ))== 0.0000000000000000 ," 511201 has zero counts if not has_gamma")
1550 insist(sum(this%binned(:, 462 ))== 0.0000000000000000 ," 511210 has zero counts if not has_gamma")
1554 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.61206787508708693 , new=",new,")"
1557 insist(sum(this%binned(:, 464 ))== 0.0000000000000000 ," 511221 has zero counts if not has_gamma")
1559 insist(sum(this%binned(:, 465 ))== 0.0000000000000000 ," 511230 has zero counts if not has_gamma")
1563 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.41775164849608570 , new=",new,")"
1568 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.7386910510827434E-002 , new=",new,")"
1573 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 7.1726213357862861E-002 , new=",new,")"
1578 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.28249918769960813 , new=",new,")"
1583 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.57419892530658634 , new=",new,")"
1588 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.31553111973075260 , new=",new,")"
1591 insist(sum(this%binned(:, 472 ))== 0.0000000000000000 ," 521180 has zero counts if not has_gamma")
1593 insist(sum(this%binned(:, 473 ))== 0.0000000000000000 ," 521190 has zero counts if not has_gamma")
1595 insist(sum(this%binned(:, 474 ))== 0.0000000000000000 ," 521191 has zero counts if not has_gamma")
1597 insist(sum(this%binned(:, 475 ))== 0.0000000000000000 ," 521200 has zero counts if not has_gamma")
1599 insist(sum(this%binned(:, 476 ))== 0.0000000000000000 ," 521210 has zero counts if not has_gamma")
1601 insist(sum(this%binned(:, 477 ))== 0.0000000000000000 ," 521211 has zero counts if not has_gamma")
1603 insist(sum(this%binned(:, 478 ))== 0.0000000000000000 ," 521220 has zero counts if not has_gamma")
1605 insist(sum(this%binned(:, 479 ))== 0.0000000000000000 ," 521230 has zero counts if not has_gamma")
1607 insist(sum(this%binned(:, 480 ))== 0.0000000000000000 ," 521231 has zero counts if not has_gamma")
1609 insist(sum(this%binned(:, 481 ))== 0.0000000000000000 ," 521240 has zero counts if not has_gamma")
1611 insist(sum(this%binned(:, 482 ))== 0.0000000000000000 ," 521250 has zero counts if not has_gamma")
1613 insist(sum(this%binned(:, 483 ))== 0.0000000000000000 ," 521251 has zero counts if not has_gamma")
1615 insist(sum(this%binned(:, 484 ))== 0.0000000000000000 ," 521260 has zero counts if not has_gamma")
1619 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.21533089173324707 , new=",new,")"
1624 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.4553212944834476E-004 , new=",new,")"
1627 insist(sum(this%binned(:, 487 ))== 0.0000000000000000 ," 521280 has zero counts if not has_gamma")
1631 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.57148233870467513 , new=",new,")"
1636 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 7.7115786881851278E-002 , new=",new,")"
1639 insist(sum(this%binned(:, 490 ))== 0.0000000000000000 ," 521300 has zero counts if not has_gamma")
1643 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.79543965914349735 , new=",new,")"
1648 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.11407246174684361 , new=",new,")"
1653 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 4.2364425689839891E-002 , new=",new,")"
1656 insist(sum(this%binned(:, 494 ))== 0.0000000000000000 ," 531220 has zero counts if not has_gamma")
1658 insist(sum(this%binned(:, 495 ))== 0.0000000000000000 ," 531230 has zero counts if not has_gamma")
1660 insist(sum(this%binned(:, 496 ))== 0.0000000000000000 ," 531240 has zero counts if not has_gamma")
1662 insist(sum(this%binned(:, 497 ))== 0.0000000000000000 ," 531250 has zero counts if not has_gamma")
1666 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 6.9957713445655045E-002 , new=",new,")"
1669 insist(sum(this%binned(:, 499 ))== 0.0000000000000000 ," 531270 has zero counts if not has_gamma")
1673 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.81150481878300196 , new=",new,")"
1678 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 2.5070131073676859E-002 , new=",new,")"
1683 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.29084157692338020 , new=",new,")"
1688 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 3.1578917906401223E-002 , new=",new,")"
1693 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.16979873602472145 , new=",new,")"
1698 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.53271029128831782 , new=",new,")"
1703 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.43136643703652339 , new=",new,")"
1706 insist(sum(this%binned(:, 507 ))== 0.0000000000000000 ," 541220 has zero counts if not has_gamma")
1708 insist(sum(this%binned(:, 508 ))== 0.0000000000000000 ," 541230 has zero counts if not has_gamma")
1710 insist(sum(this%binned(:, 509 ))== 0.0000000000000000 ," 541240 has zero counts if not has_gamma")
1714 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 5.0064228739966918E-003 , new=",new,")"
1717 insist(sum(this%binned(:, 511 ))== 0.0000000000000000 ," 541251 has zero counts if not has_gamma")
1719 insist(sum(this%binned(:, 512 ))== 0.0000000000000000 ," 541260 has zero counts if not has_gamma")
1721 insist(sum(this%binned(:, 513 ))== 0.0000000000000000 ," 541270 has zero counts if not has_gamma")
1723 insist(sum(this%binned(:, 514 ))== 0.0000000000000000 ," 541271 has zero counts if not has_gamma")
1725 insist(sum(this%binned(:, 515 ))== 0.0000000000000000 ," 541280 has zero counts if not has_gamma")
1727 insist(sum(this%binned(:, 516 ))== 0.0000000000000000 ," 541290 has zero counts if not has_gamma")
1729 insist(sum(this%binned(:, 517 ))== 0.0000000000000000 ," 541291 has zero counts if not has_gamma")
1731 insist(sum(this%binned(:, 518 ))== 0.0000000000000000 ," 541300 has zero counts if not has_gamma")
1733 insist(sum(this%binned(:, 519 ))== 0.0000000000000000 ," 541310 has zero counts if not has_gamma")
1735 insist(sum(this%binned(:, 520 ))== 0.0000000000000000 ," 541311 has zero counts if not has_gamma")
1737 insist(sum(this%binned(:, 521 ))== 0.0000000000000000 ," 541320 has zero counts if not has_gamma")
1741 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.2721398470994259E-002 , new=",new,")"
1744 insist(sum(this%binned(:, 523 ))== 0.0000000000000000 ," 541331 has zero counts if not has_gamma")
1746 insist(sum(this%binned(:, 524 ))== 0.0000000000000000 ," 541340 has zero counts if not has_gamma")
1750 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.30621217934464795 , new=",new,")"
1755 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.2799521752700906E-005 , new=",new,")"
1758 insist(sum(this%binned(:, 527 ))== 0.0000000000000000 ," 541360 has zero counts if not has_gamma")
1765 insist(sum(this%binned(:, 529 ))== 0.0000000000000000 ," 551280 has zero counts if not has_gamma")
1767 insist(sum(this%binned(:, 530 ))== 0.0000000000000000 ," 551290 has zero counts if not has_gamma")
1769 insist(sum(this%binned(:, 531 ))== 0.0000000000000000 ," 551300 has zero counts if not has_gamma")
1771 insist(sum(this%binned(:, 532 ))== 0.0000000000000000 ," 551310 has zero counts if not has_gamma")
1775 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 2.1000282732913530E-003 , new=",new,")"
1778 insist(sum(this%binned(:, 534 ))== 0.0000000000000000 ," 551330 has zero counts if not has_gamma")
1782 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.14673576315679460 , new=",new,")"
1785 insist(sum(this%binned(:, 536 ))== 0.0000000000000000 ," 551341 has zero counts if not has_gamma")
1789 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 3.9443880041662993E-002 , new=",new,")"
1794 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.5233914775956710E-002 , new=",new,")"
1799 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.15841674472337336 , new=",new,")"
1807 insist(sum(this%binned(:, 541 ))== 0.0000000000000000 ," 561280 has zero counts if not has_gamma")
1809 insist(sum(this%binned(:, 542 ))== 0.0000000000000000 ," 561290 has zero counts if not has_gamma")
1811 insist(sum(this%binned(:, 543 ))== 0.0000000000000000 ," 561300 has zero counts if not has_gamma")
1815 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 4.0519245755252231E-006 , new=",new,")"
1818 insist(sum(this%binned(:, 545 ))== 0.0000000000000000 ," 561311 has zero counts if not has_gamma")
1820 insist(sum(this%binned(:, 546 ))== 0.0000000000000000 ," 561320 has zero counts if not has_gamma")
1822 insist(sum(this%binned(:, 547 ))== 0.0000000000000000 ," 561330 has zero counts if not has_gamma")
1824 insist(sum(this%binned(:, 548 ))== 0.0000000000000000 ," 561331 has zero counts if not has_gamma")
1826 insist(sum(this%binned(:, 549 ))== 0.0000000000000000 ," 561340 has zero counts if not has_gamma")
1828 insist(sum(this%binned(:, 550 ))== 0.0000000000000000 ," 561350 has zero counts if not has_gamma")
1830 insist(sum(this%binned(:, 551 ))== 0.0000000000000000 ," 561351 has zero counts if not has_gamma")
1832 insist(sum(this%binned(:, 552 ))== 0.0000000000000000 ," 561360 has zero counts if not has_gamma")
1834 insist(sum(this%binned(:, 553 ))== 0.0000000000000000 ," 561361 has zero counts if not has_gamma")
1836 insist(sum(this%binned(:, 554 ))== 0.0000000000000000 ," 561370 has zero counts if not has_gamma")
1838 insist(sum(this%binned(:, 555 ))== 0.0000000000000000 ," 561371 has zero counts if not has_gamma")
1840 insist(sum(this%binned(:, 556 ))== 0.0000000000000000 ," 561380 has zero counts if not has_gamma")
1849 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.29417625363541500 , new=",new,")"
1857 insist(sum(this%binned(:, 560 ))== 0.0000000000000000 ," 571340 has zero counts if not has_gamma")
1859 insist(sum(this%binned(:, 561 ))== 0.0000000000000000 ," 571350 has zero counts if not has_gamma")
1861 insist(sum(this%binned(:, 562 ))== 0.0000000000000000 ," 571360 has zero counts if not has_gamma")
1863 insist(sum(this%binned(:, 563 ))== 0.0000000000000000 ," 571370 has zero counts if not has_gamma")
1867 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.0724458313241581E-003 , new=",new,")"
1870 insist(sum(this%binned(:, 565 ))== 0.0000000000000000 ," 571390 has zero counts if not has_gamma")
1874 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.60218382412752902 , new=",new,")"
1882 insist(sum(this%binned(:, 568 ))== 0.0000000000000000 ," 581340 has zero counts if not has_gamma")
1884 insist(sum(this%binned(:, 569 ))== 0.0000000000000000 ," 581350 has zero counts if not has_gamma")
1886 insist(sum(this%binned(:, 570 ))== 0.0000000000000000 ," 581360 has zero counts if not has_gamma")
1888 insist(sum(this%binned(:, 571 ))== 0.0000000000000000 ," 581370 has zero counts if not has_gamma")
1890 insist(sum(this%binned(:, 572 ))== 0.0000000000000000 ," 581371 has zero counts if not has_gamma")
1892 insist(sum(this%binned(:, 573 ))== 0.0000000000000000 ," 581380 has zero counts if not has_gamma")
1894 insist(sum(this%binned(:, 574 ))== 0.0000000000000000 ," 581390 has zero counts if not has_gamma")
1896 insist(sum(this%binned(:, 575 ))== 0.0000000000000000 ," 581391 has zero counts if not has_gamma")
1898 insist(sum(this%binned(:, 576 ))== 0.0000000000000000 ," 581400 has zero counts if not has_gamma")
1902 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.13008160678741820 , new=",new,")"
1905 insist(sum(this%binned(:, 578 ))== 0.0000000000000000 ," 581420 has zero counts if not has_gamma")
1909 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.42580739165593329 , new=",new,")"
1914 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 6.3751021087238499E-002 , new=",new,")"
1919 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.68194909056210118 , new=",new,")"
1922 insist(sum(this%binned(:, 582 ))== 0.0000000000000000 ," 591400 has zero counts if not has_gamma")
1924 insist(sum(this%binned(:, 583 ))== 0.0000000000000000 ," 591410 has zero counts if not has_gamma")
1928 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.95774008341515993 , new=",new,")"
1931 insist(sum(this%binned(:, 585 ))== 0.0000000000000000 ," 591421 has zero counts if not has_gamma")
1935 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.32059556534630923 , new=",new,")"
1945 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 2.2731732977862911E-007 , new=",new,")"
1950 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.76913395252632344 , new=",new,")"
1953 insist(sum(this%binned(:, 590 ))== 0.0000000000000000 ," 601400 has zero counts if not has_gamma")
1955 insist(sum(this%binned(:, 591 ))== 0.0000000000000000 ," 601410 has zero counts if not has_gamma")
1957 insist(sum(this%binned(:, 592 ))== 0.0000000000000000 ," 601411 has zero counts if not has_gamma")
1959 insist(sum(this%binned(:, 593 ))== 0.0000000000000000 ," 601420 has zero counts if not has_gamma")
1961 insist(sum(this%binned(:, 594 ))== 0.0000000000000000 ," 601430 has zero counts if not has_gamma")
1963 insist(sum(this%binned(:, 595 ))== 0.0000000000000000 ," 601440 has zero counts if not has_gamma")
1965 insist(sum(this%binned(:, 596 ))== 0.0000000000000000 ," 601450 has zero counts if not has_gamma")
1967 insist(sum(this%binned(:, 597 ))== 0.0000000000000000 ," 601460 has zero counts if not has_gamma")
1971 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.22875016324034583 , new=",new,")"
1974 insist(sum(this%binned(:, 599 ))== 0.0000000000000000 ," 601480 has zero counts if not has_gamma")
1978 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.50091352746460194 , new=",new,")"
1981 insist(sum(this%binned(:, 601 ))== 0.0000000000000000 ," 601500 has zero counts if not has_gamma")
1985 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.68210990210782507 , new=",new,")"
1988 insist(sum(this%binned(:, 603 ))== 0.0000000000000000 ," 611430 has zero counts if not has_gamma")
1990 insist(sum(this%binned(:, 604 ))== 0.0000000000000000 ," 611440 has zero counts if not has_gamma")
1992 insist(sum(this%binned(:, 605 ))== 0.0000000000000000 ," 611450 has zero counts if not has_gamma")
1994 insist(sum(this%binned(:, 606 ))== 0.0000000000000000 ," 611460 has zero counts if not has_gamma")
1998 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 4.4805144281838238E-002 , new=",new,")"
2003 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.83675692258202616 , new=",new,")"
2008 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.12653649152724911 , new=",new,")"
2013 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.38134159373967347 , new=",new,")"
2018 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.90416439631137979 , new=",new,")"
2023 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.28534802372394646 , new=",new,")"
2026 insist(sum(this%binned(:, 613 ))== 0.0000000000000000 ," 611520 has zero counts if not has_gamma")
2028 insist(sum(this%binned(:, 614 ))== 0.0000000000000000 ," 621440 has zero counts if not has_gamma")
2030 insist(sum(this%binned(:, 615 ))== 0.0000000000000000 ," 621450 has zero counts if not has_gamma")
2032 insist(sum(this%binned(:, 616 ))== 0.0000000000000000 ," 621460 has zero counts if not has_gamma")
2034 insist(sum(this%binned(:, 617 ))== 0.0000000000000000 ," 621470 has zero counts if not has_gamma")
2036 insist(sum(this%binned(:, 618 ))== 0.0000000000000000 ," 621480 has zero counts if not has_gamma")
2038 insist(sum(this%binned(:, 619 ))== 0.0000000000000000 ," 621490 has zero counts if not has_gamma")
2040 insist(sum(this%binned(:, 620 ))== 0.0000000000000000 ," 621500 has zero counts if not has_gamma")
2044 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.5480743712219379E-003 , new=",new,")"
2047 insist(sum(this%binned(:, 622 ))== 0.0000000000000000 ," 621520 has zero counts if not has_gamma")
2051 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.21990841593176749 , new=",new,")"
2054 insist(sum(this%binned(:, 624 ))== 0.0000000000000000 ," 621540 has zero counts if not has_gamma")
2058 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.61408630543853560 , new=",new,")"
2061 insist(sum(this%binned(:, 626 ))== 0.0000000000000000 ," 631450 has zero counts if not has_gamma")
2063 insist(sum(this%binned(:, 627 ))== 0.0000000000000000 ," 631460 has zero counts if not has_gamma")
2065 insist(sum(this%binned(:, 628 ))== 0.0000000000000000 ," 631470 has zero counts if not has_gamma")
2067 insist(sum(this%binned(:, 629 ))== 0.0000000000000000 ," 631480 has zero counts if not has_gamma")
2069 insist(sum(this%binned(:, 630 ))== 0.0000000000000000 ," 631490 has zero counts if not has_gamma")
2071 insist(sum(this%binned(:, 631 ))== 0.0000000000000000 ," 631500 has zero counts if not has_gamma")
2073 insist(sum(this%binned(:, 632 ))== 0.0000000000000000 ," 631501 has zero counts if not has_gamma")
2075 insist(sum(this%binned(:, 633 ))== 0.0000000000000000 ," 631510 has zero counts if not has_gamma")
2079 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 2.4607505074389455E-002 , new=",new,")"
2084 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.40900613295608057 , new=",new,")"
2087 insist(sum(this%binned(:, 636 ))== 0.0000000000000000 ," 631530 has zero counts if not has_gamma")
2091 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.23601550649829009 , new=",new,")"
2096 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 3.0158225960669194E-002 , new=",new,")"
2101 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.43135526533026802 , new=",new,")"
2106 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.38807598956565947 , new=",new,")"
2109 insist(sum(this%binned(:, 641 ))== 0.0000000000000000 ," 641460 has zero counts if not has_gamma")
2111 insist(sum(this%binned(:, 642 ))== 0.0000000000000000 ," 641470 has zero counts if not has_gamma")
2113 insist(sum(this%binned(:, 643 ))== 0.0000000000000000 ," 641480 has zero counts if not has_gamma")
2115 insist(sum(this%binned(:, 644 ))== 0.0000000000000000 ," 641490 has zero counts if not has_gamma")
2117 insist(sum(this%binned(:, 645 ))== 0.0000000000000000 ," 641500 has zero counts if not has_gamma")
2119 insist(sum(this%binned(:, 646 ))== 0.0000000000000000 ," 641510 has zero counts if not has_gamma")
2121 insist(sum(this%binned(:, 647 ))== 0.0000000000000000 ," 641520 has zero counts if not has_gamma")
2123 insist(sum(this%binned(:, 648 ))== 0.0000000000000000 ," 641530 has zero counts if not has_gamma")
2125 insist(sum(this%binned(:, 649 ))== 0.0000000000000000 ," 641540 has zero counts if not has_gamma")
2127 insist(sum(this%binned(:, 650 ))== 0.0000000000000000 ," 641551 has zero counts if not has_gamma")
2129 insist(sum(this%binned(:, 651 ))== 0.0000000000000000 ," 641550 has zero counts if not has_gamma")
2131 insist(sum(this%binned(:, 652 ))== 0.0000000000000000 ," 641560 has zero counts if not has_gamma")
2133 insist(sum(this%binned(:, 653 ))== 0.0000000000000000 ," 641570 has zero counts if not has_gamma")
2135 insist(sum(this%binned(:, 654 ))== 0.0000000000000000 ," 641580 has zero counts if not has_gamma")
2139 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.29955126869737125 , new=",new,")"
2142 insist(sum(this%binned(:, 656 ))== 0.0000000000000000 ," 641600 has zero counts if not has_gamma")
2146 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.62205631423809715 , new=",new,")"
2151 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.33008241509130104 , new=",new,")"
2154 insist(sum(this%binned(:, 659 ))== 0.0000000000000000 ," 651520 has zero counts if not has_gamma")
2156 insist(sum(this%binned(:, 660 ))== 0.0000000000000000 ," 651530 has zero counts if not has_gamma")
2158 insist(sum(this%binned(:, 661 ))== 0.0000000000000000 ," 651540 has zero counts if not has_gamma")
2160 insist(sum(this%binned(:, 662 ))== 0.0000000000000000 ," 651541 has zero counts if not has_gamma")
2162 insist(sum(this%binned(:, 663 ))== 0.0000000000000000 ," 651550 has zero counts if not has_gamma")
2164 insist(sum(this%binned(:, 664 ))== 0.0000000000000000 ," 651560 has zero counts if not has_gamma")
2166 insist(sum(this%binned(:, 665 ))== 0.0000000000000000 ," 651561 has zero counts if not has_gamma")
2168 insist(sum(this%binned(:, 666 ))== 0.0000000000000000 ," 651570 has zero counts if not has_gamma")
2170 insist(sum(this%binned(:, 667 ))== 0.0000000000000000 ," 651580 has zero counts if not has_gamma")
2172 insist(sum(this%binned(:, 668 ))== 0.0000000000000000 ," 651590 has zero counts if not has_gamma")
2176 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.19677743286810334 , new=",new,")"
2181 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.13716866353510065 , new=",new,")"
2186 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.49386870265293037 , new=",new,")"
2189 insist(sum(this%binned(:, 672 ))== 0.0000000000000000 ," 661540 has zero counts if not has_gamma")
2191 insist(sum(this%binned(:, 673 ))== 0.0000000000000000 ," 661550 has zero counts if not has_gamma")
2193 insist(sum(this%binned(:, 674 ))== 0.0000000000000000 ," 661560 has zero counts if not has_gamma")
2195 insist(sum(this%binned(:, 675 ))== 0.0000000000000000 ," 661570 has zero counts if not has_gamma")
2197 insist(sum(this%binned(:, 676 ))== 0.0000000000000000 ," 661580 has zero counts if not has_gamma")
2199 insist(sum(this%binned(:, 677 ))== 0.0000000000000000 ," 661590 has zero counts if not has_gamma")
2201 insist(sum(this%binned(:, 678 ))== 0.0000000000000000 ," 661600 has zero counts if not has_gamma")
2203 insist(sum(this%binned(:, 679 ))== 0.0000000000000000 ," 661610 has zero counts if not has_gamma")
2205 insist(sum(this%binned(:, 680 ))== 0.0000000000000000 ," 661620 has zero counts if not has_gamma")
2207 insist(sum(this%binned(:, 681 ))== 0.0000000000000000 ," 661630 has zero counts if not has_gamma")
2209 insist(sum(this%binned(:, 682 ))== 0.0000000000000000 ," 661640 has zero counts if not has_gamma")
2213 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.47595028997703481 , new=",new,")"
2218 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.5781428807830995E-004 , new=",new,")"
2223 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 9.8151209522421748E-002 , new=",new,")"
2226 insist(sum(this%binned(:, 686 ))== 0.0000000000000000 ," 671600 has zero counts if not has_gamma")
2228 insist(sum(this%binned(:, 687 ))== 0.0000000000000000 ," 671601 has zero counts if not has_gamma")
2230 insist(sum(this%binned(:, 688 ))== 0.0000000000000000 ," 671610 has zero counts if not has_gamma")
2232 insist(sum(this%binned(:, 689 ))== 0.0000000000000000 ," 671630 has zero counts if not has_gamma")
2234 insist(sum(this%binned(:, 690 ))== 0.0000000000000000 ," 671631 has zero counts if not has_gamma")
2236 insist(sum(this%binned(:, 691 ))== 0.0000000000000000 ," 671640 has zero counts if not has_gamma")
2238 insist(sum(this%binned(:, 692 ))== 0.0000000000000000 ," 671641 has zero counts if not has_gamma")
2240 insist(sum(this%binned(:, 693 ))== 0.0000000000000000 ," 671650 has zero counts if not has_gamma")
2244 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.76381537700397173 , new=",new,")"
2249 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 6.9905578106101993E-003 , new=",new,")"
2252 insist(sum(this%binned(:, 696 ))== 0.0000000000000000 ," 681600 has zero counts if not has_gamma")
2254 insist(sum(this%binned(:, 697 ))== 0.0000000000000000 ," 681610 has zero counts if not has_gamma")
2256 insist(sum(this%binned(:, 698 ))== 0.0000000000000000 ," 681620 has zero counts if not has_gamma")
2258 insist(sum(this%binned(:, 699 ))== 0.0000000000000000 ," 681630 has zero counts if not has_gamma")
2260 insist(sum(this%binned(:, 700 ))== 0.0000000000000000 ," 681640 has zero counts if not has_gamma")
2262 insist(sum(this%binned(:, 701 ))== 0.0000000000000000 ," 681650 has zero counts if not has_gamma")
2264 insist(sum(this%binned(:, 702 ))== 0.0000000000000000 ," 681660 has zero counts if not has_gamma")
2266 insist(sum(this%binned(:, 703 ))== 0.0000000000000000 ," 681670 has zero counts if not has_gamma")
2268 insist(sum(this%binned(:, 704 ))== 0.0000000000000000 ," 681671 has zero counts if not has_gamma")
2270 insist(sum(this%binned(:, 705 ))== 0.0000000000000000 ," 681680 has zero counts if not has_gamma")
2274 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.2603648633007620E-002 , new=",new,")"
2277 insist(sum(this%binned(:, 707 ))== 0.0000000000000000 ," 681700 has zero counts if not has_gamma")
2281 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.38508704871806548 , new=",new,")"
2286 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 7.9573185132163951E-002 , new=",new,")"
2289 insist(sum(this%binned(:, 710 ))== 0.0000000000000000 ," 691650 has zero counts if not has_gamma")
2291 insist(sum(this%binned(:, 711 ))== 0.0000000000000000 ," 691660 has zero counts if not has_gamma")
2293 insist(sum(this%binned(:, 712 ))== 0.0000000000000000 ," 691670 has zero counts if not has_gamma")
2295 insist(sum(this%binned(:, 713 ))== 0.0000000000000000 ," 691680 has zero counts if not has_gamma")
2297 insist(sum(this%binned(:, 714 ))== 0.0000000000000000 ," 691690 has zero counts if not has_gamma")
2301 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.32198517535952448 , new=",new,")"
2306 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.2421613035298590E-002 , new=",new,")"
2311 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.54836751755446267 , new=",new,")"
2316 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.30832883353377594 , new=",new,")"
2319 insist(sum(this%binned(:, 719 ))== 0.0000000000000000 ," 701660 has zero counts if not has_gamma")
2321 insist(sum(this%binned(:, 720 ))== 0.0000000000000000 ," 701670 has zero counts if not has_gamma")
2323 insist(sum(this%binned(:, 721 ))== 0.0000000000000000 ," 701680 has zero counts if not has_gamma")
2325 insist(sum(this%binned(:, 722 ))== 0.0000000000000000 ," 701690 has zero counts if not has_gamma")
2327 insist(sum(this%binned(:, 723 ))== 0.0000000000000000 ," 701700 has zero counts if not has_gamma")
2329 insist(sum(this%binned(:, 724 ))== 0.0000000000000000 ," 701710 has zero counts if not has_gamma")
2331 insist(sum(this%binned(:, 725 ))== 0.0000000000000000 ," 701720 has zero counts if not has_gamma")
2333 insist(sum(this%binned(:, 726 ))== 0.0000000000000000 ," 701730 has zero counts if not has_gamma")
2335 insist(sum(this%binned(:, 727 ))== 0.0000000000000000 ," 701740 has zero counts if not has_gamma")
2339 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.11027029078394435 , new=",new,")"
2342 insist(sum(this%binned(:, 729 ))== 0.0000000000000000 ," 701751 has zero counts if not has_gamma")
2344 insist(sum(this%binned(:, 730 ))== 0.0000000000000000 ," 701760 has zero counts if not has_gamma")
2348 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.45213056628767700 , new=",new,")"
2351 insist(sum(this%binned(:, 732 ))== 0.0000000000000000 ," 711690 has zero counts if not has_gamma")
2353 insist(sum(this%binned(:, 733 ))== 0.0000000000000000 ," 711700 has zero counts if not has_gamma")
2355 insist(sum(this%binned(:, 734 ))== 0.0000000000000000 ," 711710 has zero counts if not has_gamma")
2357 insist(sum(this%binned(:, 735 ))== 0.0000000000000000 ," 711720 has zero counts if not has_gamma")
2359 insist(sum(this%binned(:, 736 ))== 0.0000000000000000 ," 711721 has zero counts if not has_gamma")
2361 insist(sum(this%binned(:, 737 ))== 0.0000000000000000 ," 711730 has zero counts if not has_gamma")
2363 insist(sum(this%binned(:, 738 ))== 0.0000000000000000 ," 711740 has zero counts if not has_gamma")
2365 insist(sum(this%binned(:, 739 ))== 0.0000000000000000 ," 711741 has zero counts if not has_gamma")
2367 insist(sum(this%binned(:, 740 ))== 0.0000000000000000 ," 711750 has zero counts if not has_gamma")
2371 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.16774423530656077 , new=",new,")"
2376 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.46985723134678636 , new=",new,")"
2381 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.11788118328851648 , new=",new,")"
2386 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.25134608291358329 , new=",new,")"
2389 insist(sum(this%binned(:, 745 ))== 0.0000000000000000 ," 721700 has zero counts if not has_gamma")
2391 insist(sum(this%binned(:, 746 ))== 0.0000000000000000 ," 721710 has zero counts if not has_gamma")
2393 insist(sum(this%binned(:, 747 ))== 0.0000000000000000 ," 721720 has zero counts if not has_gamma")
2395 insist(sum(this%binned(:, 748 ))== 0.0000000000000000 ," 721730 has zero counts if not has_gamma")
2397 insist(sum(this%binned(:, 749 ))== 0.0000000000000000 ," 721740 has zero counts if not has_gamma")
2399 insist(sum(this%binned(:, 750 ))== 0.0000000000000000 ," 721750 has zero counts if not has_gamma")
2401 insist(sum(this%binned(:, 751 ))== 0.0000000000000000 ," 721760 has zero counts if not has_gamma")
2403 insist(sum(this%binned(:, 752 ))== 0.0000000000000000 ," 721770 has zero counts if not has_gamma")
2405 insist(sum(this%binned(:, 753 ))== 0.0000000000000000 ," 721771 has zero counts if not has_gamma")
2407 insist(sum(this%binned(:, 754 ))== 0.0000000000000000 ," 721780 has zero counts if not has_gamma")
2409 insist(sum(this%binned(:, 755 ))== 0.0000000000000000 ," 721781 has zero counts if not has_gamma")
2411 insist(sum(this%binned(:, 756 ))== 0.0000000000000000 ," 721790 has zero counts if not has_gamma")
2413 insist(sum(this%binned(:, 757 ))== 0.0000000000000000 ," 721791 has zero counts if not has_gamma")
2415 insist(sum(this%binned(:, 758 ))== 0.0000000000000000 ," 721800 has zero counts if not has_gamma")
2417 insist(sum(this%binned(:, 759 ))== 0.0000000000000000 ," 721801 has zero counts if not has_gamma")
2421 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.10216851044242503 , new=",new,")"
2426 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 3.2873713067093753E-002 , new=",new,")"
2429 insist(sum(this%binned(:, 762 ))== 0.0000000000000000 ," 731770 has zero counts if not has_gamma")
2431 insist(sum(this%binned(:, 763 ))== 0.0000000000000000 ," 731780 has zero counts if not has_gamma")
2433 insist(sum(this%binned(:, 764 ))== 0.0000000000000000 ," 731790 has zero counts if not has_gamma")
2435 insist(sum(this%binned(:, 765 ))== 0.0000000000000000 ," 731801 has zero counts if not has_gamma")
2437 insist(sum(this%binned(:, 766 ))== 0.0000000000000000 ," 731800 has zero counts if not has_gamma")
2439 insist(sum(this%binned(:, 767 ))== 0.0000000000000000 ," 731810 has zero counts if not has_gamma")
2443 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.11098766851437693 , new=",new,")"
2446 insist(sum(this%binned(:, 769 ))== 0.0000000000000000 ," 731821 has zero counts if not has_gamma")
2450 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.17589090681543226 , new=",new,")"
2453 insist(sum(this%binned(:, 771 ))== 0.0000000000000000 ," 741780 has zero counts if not has_gamma")
2455 insist(sum(this%binned(:, 772 ))== 0.0000000000000000 ," 741800 has zero counts if not has_gamma")
2457 insist(sum(this%binned(:, 773 ))== 0.0000000000000000 ," 741810 has zero counts if not has_gamma")
2459 insist(sum(this%binned(:, 774 ))== 0.0000000000000000 ," 741820 has zero counts if not has_gamma")
2461 insist(sum(this%binned(:, 775 ))== 0.0000000000000000 ," 741831 has zero counts if not has_gamma")
2463 insist(sum(this%binned(:, 776 ))== 0.0000000000000000 ," 741830 has zero counts if not has_gamma")
2465 insist(sum(this%binned(:, 777 ))== 0.0000000000000000 ," 741840 has zero counts if not has_gamma")
2469 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.11064296651577800 , new=",new,")"
2472 insist(sum(this%binned(:, 779 ))== 0.0000000000000000 ," 741851 has zero counts if not has_gamma")
2474 insist(sum(this%binned(:, 780 ))== 0.0000000000000000 ," 741860 has zero counts if not has_gamma")
2478 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.26596432951907234 , new=",new,")"
2483 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.1702264111754630E-002 , new=",new,")"
2486 insist(sum(this%binned(:, 783 ))== 0.0000000000000000 ," 741890 has zero counts if not has_gamma")
2488 insist(sum(this%binned(:, 784 ))== 0.0000000000000000 ," 751810 has zero counts if not has_gamma")
2490 insist(sum(this%binned(:, 785 ))== 0.0000000000000000 ," 751820 has zero counts if not has_gamma")
2492 insist(sum(this%binned(:, 786 ))== 0.0000000000000000 ," 751821 has zero counts if not has_gamma")
2494 insist(sum(this%binned(:, 787 ))== 0.0000000000000000 ," 751830 has zero counts if not has_gamma")
2496 insist(sum(this%binned(:, 788 ))== 0.0000000000000000 ," 751840 has zero counts if not has_gamma")
2498 insist(sum(this%binned(:, 789 ))== 0.0000000000000000 ," 751841 has zero counts if not has_gamma")
2500 insist(sum(this%binned(:, 790 ))== 0.0000000000000000 ," 751850 has zero counts if not has_gamma")
2504 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.32251814027147196 , new=",new,")"
2507 insist(sum(this%binned(:, 792 ))== 0.0000000000000000 ," 751861 has zero counts if not has_gamma")
2511 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 3.6429999370878735E-014 , new=",new,")"
2516 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.88296889652508559 , new=",new,")"
2519 insist(sum(this%binned(:, 795 ))== 0.0000000000000000 ," 751881 has zero counts if not has_gamma")
2523 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.38755910558928969 , new=",new,")"
2526 insist(sum(this%binned(:, 797 ))== 0.0000000000000000 ," 761820 has zero counts if not has_gamma")
2528 insist(sum(this%binned(:, 798 ))== 0.0000000000000000 ," 761830 has zero counts if not has_gamma")
2530 insist(sum(this%binned(:, 799 ))== 0.0000000000000000 ," 761840 has zero counts if not has_gamma")
2532 insist(sum(this%binned(:, 800 ))== 0.0000000000000000 ," 761850 has zero counts if not has_gamma")
2534 insist(sum(this%binned(:, 801 ))== 0.0000000000000000 ," 761860 has zero counts if not has_gamma")
2536 insist(sum(this%binned(:, 802 ))== 0.0000000000000000 ," 761870 has zero counts if not has_gamma")
2538 insist(sum(this%binned(:, 803 ))== 0.0000000000000000 ," 761880 has zero counts if not has_gamma")
2540 insist(sum(this%binned(:, 804 ))== 0.0000000000000000 ," 761890 has zero counts if not has_gamma")
2542 insist(sum(this%binned(:, 805 ))== 0.0000000000000000 ," 761891 has zero counts if not has_gamma")
2544 insist(sum(this%binned(:, 806 ))== 0.0000000000000000 ," 761900 has zero counts if not has_gamma")
2546 insist(sum(this%binned(:, 807 ))== 0.0000000000000000 ," 761901 has zero counts if not has_gamma")
2550 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 2.2744864753734710E-002 , new=",new,")"
2553 insist(sum(this%binned(:, 809 ))== 0.0000000000000000 ," 761911 has zero counts if not has_gamma")
2555 insist(sum(this%binned(:, 810 ))== 0.0000000000000000 ," 761920 has zero counts if not has_gamma")
2557 insist(sum(this%binned(:, 811 ))== 0.0000000000000000 ," 761930 has zero counts if not has_gamma")
2561 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.0053409411032190E-002 , new=",new,")"
2564 insist(sum(this%binned(:, 813 ))== 0.0000000000000000 ," 771850 has zero counts if not has_gamma")
2566 insist(sum(this%binned(:, 814 ))== 0.0000000000000000 ," 771860 has zero counts if not has_gamma")
2568 insist(sum(this%binned(:, 815 ))== 0.0000000000000000 ," 771880 has zero counts if not has_gamma")
2570 insist(sum(this%binned(:, 816 ))== 0.0000000000000000 ," 771890 has zero counts if not has_gamma")
2572 insist(sum(this%binned(:, 817 ))== 0.0000000000000000 ," 771891 has zero counts if not has_gamma")
2574 insist(sum(this%binned(:, 818 ))== 0.0000000000000000 ," 771900 has zero counts if not has_gamma")
2576 insist(sum(this%binned(:, 819 ))== 0.0000000000000000 ," 771910 has zero counts if not has_gamma")
2578 insist(sum(this%binned(:, 820 ))== 0.0000000000000000 ," 771911 has zero counts if not has_gamma")
2582 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.15271818383563138 , new=",new,")"
2587 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.3639700666395634E-005 , new=",new,")"
2590 insist(sum(this%binned(:, 823 ))== 0.0000000000000000 ," 771930 has zero counts if not has_gamma")
2592 insist(sum(this%binned(:, 824 ))== 0.0000000000000000 ," 771931 has zero counts if not has_gamma")
2596 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.93501711795654774 , new=",new,")"
2601 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 5.2449764006308186E-002 , new=",new,")"
2604 insist(sum(this%binned(:, 827 ))== 0.0000000000000000 ," 771960 has zero counts if not has_gamma")
2606 insist(sum(this%binned(:, 828 ))== 0.0000000000000000 ," 771961 has zero counts if not has_gamma")
2608 insist(sum(this%binned(:, 829 ))== 0.0000000000000000 ," 781880 has zero counts if not has_gamma")
2610 insist(sum(this%binned(:, 830 ))== 0.0000000000000000 ," 781890 has zero counts if not has_gamma")
2612 insist(sum(this%binned(:, 831 ))== 0.0000000000000000 ," 781900 has zero counts if not has_gamma")
2614 insist(sum(this%binned(:, 832 ))== 0.0000000000000000 ," 781910 has zero counts if not has_gamma")
2616 insist(sum(this%binned(:, 833 ))== 0.0000000000000000 ," 781920 has zero counts if not has_gamma")
2618 insist(sum(this%binned(:, 834 ))== 0.0000000000000000 ," 781930 has zero counts if not has_gamma")
2620 insist(sum(this%binned(:, 835 ))== 0.0000000000000000 ," 781931 has zero counts if not has_gamma")
2622 insist(sum(this%binned(:, 836 ))== 0.0000000000000000 ," 781940 has zero counts if not has_gamma")
2624 insist(sum(this%binned(:, 837 ))== 0.0000000000000000 ," 781950 has zero counts if not has_gamma")
2626 insist(sum(this%binned(:, 838 ))== 0.0000000000000000 ," 781951 has zero counts if not has_gamma")
2628 insist(sum(this%binned(:, 839 ))== 0.0000000000000000 ," 781960 has zero counts if not has_gamma")
2632 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.18554690397005824 , new=",new,")"
2637 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 7.0945916320679387E-003 , new=",new,")"
2640 insist(sum(this%binned(:, 842 ))== 0.0000000000000000 ," 781980 has zero counts if not has_gamma")
2644 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.57670183556534194 , new=",new,")"
2647 insist(sum(this%binned(:, 844 ))== 0.0000000000000000 ," 781991 has zero counts if not has_gamma")
2649 insist(sum(this%binned(:, 845 ))== 0.0000000000000000 ," 782000 has zero counts if not has_gamma")
2651 insist(sum(this%binned(:, 846 ))== 0.0000000000000000 ," 791930 has zero counts if not has_gamma")
2653 insist(sum(this%binned(:, 847 ))== 0.0000000000000000 ," 791940 has zero counts if not has_gamma")
2655 insist(sum(this%binned(:, 848 ))== 0.0000000000000000 ," 791950 has zero counts if not has_gamma")
2657 insist(sum(this%binned(:, 849 ))== 0.0000000000000000 ," 791951 has zero counts if not has_gamma")
2659 insist(sum(this%binned(:, 850 ))== 0.0000000000000000 ," 791960 has zero counts if not has_gamma")
2661 insist(sum(this%binned(:, 851 ))== 0.0000000000000000 ," 791970 has zero counts if not has_gamma")
2663 insist(sum(this%binned(:, 852 ))== 0.0000000000000000 ," 791971 has zero counts if not has_gamma")
2667 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.31966403609022070 , new=",new,")"
2670 insist(sum(this%binned(:, 854 ))== 0.0000000000000000 ," 791981 has zero counts if not has_gamma")
2674 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 6.9247868059394008E-002 , new=",new,")"
2677 insist(sum(this%binned(:, 856 ))== 0.0000000000000000 ," 792000 has zero counts if not has_gamma")
2679 insist(sum(this%binned(:, 857 ))== 0.0000000000000000 ," 792001 has zero counts if not has_gamma")
2681 insist(sum(this%binned(:, 858 ))== 0.0000000000000000 ," 801930 has zero counts if not has_gamma")
2683 insist(sum(this%binned(:, 859 ))== 0.0000000000000000 ," 801931 has zero counts if not has_gamma")
2685 insist(sum(this%binned(:, 860 ))== 0.0000000000000000 ," 801940 has zero counts if not has_gamma")
2687 insist(sum(this%binned(:, 861 ))== 0.0000000000000000 ," 801950 has zero counts if not has_gamma")
2689 insist(sum(this%binned(:, 862 ))== 0.0000000000000000 ," 801951 has zero counts if not has_gamma")
2691 insist(sum(this%binned(:, 863 ))== 0.0000000000000000 ," 801960 has zero counts if not has_gamma")
2693 insist(sum(this%binned(:, 864 ))== 0.0000000000000000 ," 801970 has zero counts if not has_gamma")
2695 insist(sum(this%binned(:, 865 ))== 0.0000000000000000 ," 801971 has zero counts if not has_gamma")
2697 insist(sum(this%binned(:, 866 ))== 0.0000000000000000 ," 801980 has zero counts if not has_gamma")
2699 insist(sum(this%binned(:, 867 ))== 0.0000000000000000 ," 801990 has zero counts if not has_gamma")
2701 insist(sum(this%binned(:, 868 ))== 0.0000000000000000 ," 801991 has zero counts if not has_gamma")
2703 insist(sum(this%binned(:, 869 ))== 0.0000000000000000 ," 802000 has zero counts if not has_gamma")
2705 insist(sum(this%binned(:, 870 ))== 0.0000000000000000 ," 802010 has zero counts if not has_gamma")
2707 insist(sum(this%binned(:, 871 ))== 0.0000000000000000 ," 802020 has zero counts if not has_gamma")
2711 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 4.0785911570608853E-002 , new=",new,")"
2714 insist(sum(this%binned(:, 873 ))== 0.0000000000000000 ," 802040 has zero counts if not has_gamma")
2716 insist(sum(this%binned(:, 874 ))== 0.0000000000000000 ," 802050 has zero counts if not has_gamma")
2718 insist(sum(this%binned(:, 875 ))== 0.0000000000000000 ," 802060 has zero counts if not has_gamma")
2720 insist(sum(this%binned(:, 876 ))== 0.0000000000000000 ," 812000 has zero counts if not has_gamma")
2722 insist(sum(this%binned(:, 877 ))== 0.0000000000000000 ," 812010 has zero counts if not has_gamma")
2724 insist(sum(this%binned(:, 878 ))== 0.0000000000000000 ," 812020 has zero counts if not has_gamma")
2726 insist(sum(this%binned(:, 879 ))== 0.0000000000000000 ," 812030 has zero counts if not has_gamma")
2728 insist(sum(this%binned(:, 880 ))== 0.0000000000000000 ," 812040 has zero counts if not has_gamma")
2730 insist(sum(this%binned(:, 881 ))== 0.0000000000000000 ," 812050 has zero counts if not has_gamma")
2732 insist(sum(this%binned(:, 882 ))== 0.0000000000000000 ," 812060 has zero counts if not has_gamma")
2736 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.53874462950932056 , new=",new,")"
2739 insist(sum(this%binned(:, 884 ))== 0.0000000000000000 ," 822000 has zero counts if not has_gamma")
2741 insist(sum(this%binned(:, 885 ))== 0.0000000000000000 ," 822020 has zero counts if not has_gamma")
2743 insist(sum(this%binned(:, 886 ))== 0.0000000000000000 ," 822030 has zero counts if not has_gamma")
2745 insist(sum(this%binned(:, 887 ))== 0.0000000000000000 ," 822040 has zero counts if not has_gamma")
2747 insist(sum(this%binned(:, 888 ))== 0.0000000000000000 ," 822050 has zero counts if not has_gamma")
2749 insist(sum(this%binned(:, 889 ))== 0.0000000000000000 ," 822051 has zero counts if not has_gamma")
2751 insist(sum(this%binned(:, 890 ))== 0.0000000000000000 ," 822060 has zero counts if not has_gamma")
2753 insist(sum(this%binned(:, 891 ))== 0.0000000000000000 ," 822070 has zero counts if not has_gamma")
2755 insist(sum(this%binned(:, 892 ))== 0.0000000000000000 ," 822071 has zero counts if not has_gamma")
2757 insist(sum(this%binned(:, 893 ))== 0.0000000000000000 ," 822080 has zero counts if not has_gamma")
2761 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.18797484500782954 , new=",new,")"
2766 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.4747836263975245E-003 , new=",new,")"
2769 insist(sum(this%binned(:, 896 ))== 0.0000000000000000 ," 832050 has zero counts if not has_gamma")
2771 insist(sum(this%binned(:, 897 ))== 0.0000000000000000 ," 832060 has zero counts if not has_gamma")
2773 insist(sum(this%binned(:, 898 ))== 0.0000000000000000 ," 832070 has zero counts if not has_gamma")
2775 insist(sum(this%binned(:, 899 ))== 0.0000000000000000 ," 832080 has zero counts if not has_gamma")
2777 insist(sum(this%binned(:, 900 ))== 0.0000000000000000 ," 832090 has zero counts if not has_gamma")
2781 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.41111443877903026 , new=",new,")"
2784 insist(sum(this%binned(:, 902 ))== 0.0000000000000000 ," 832101 has zero counts if not has_gamma")
2788 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.2738720412561945E-006 , new=",new,")"
2791 insist(sum(this%binned(:, 904 ))== 0.0000000000000000 ," 842060 has zero counts if not has_gamma")
2793 insist(sum(this%binned(:, 905 ))== 0.0000000000000000 ," 842070 has zero counts if not has_gamma")
2795 insist(sum(this%binned(:, 906 ))== 0.0000000000000000 ," 842080 has zero counts if not has_gamma")
2797 insist(sum(this%binned(:, 907 ))== 0.0000000000000000 ," 842090 has zero counts if not has_gamma")
2799 insist(sum(this%binned(:, 908 ))== 0.0000000000000000 ," 842100 has zero counts if not has_gamma")
2801 insist(sum(this%binned(:, 909 ))== 0.0000000000000000 ," 842110 has zero counts if not has_gamma")
2803 insist(sum(this%binned(:, 910 ))== 0.0000000000000000 ," 842111 has zero counts if not has_gamma")
2805 insist(sum(this%binned(:, 911 ))== 0.0000000000000000 ," 20030 has zero counts if not has_gamma")
2807 insist(sum(this%binned(:, 912 ))== 0.0000000000000000 ," 20040 has zero counts if not has_gamma")
2809 insist(sum(this%binned(:, 913 ))== 0.0000000000000000 ," 30060 has zero counts if not has_gamma")
2811 insist(sum(this%binned(:, 914 ))== 0.0000000000000000 ," 30070 has zero counts if not has_gamma")
2813 insist(sum(this%binned(:, 915 ))== 0.0000000000000000 ," 40070 has zero counts if not has_gamma")
2815 insist(sum(this%binned(:, 916 ))== 0.0000000000000000 ," 60120 has zero counts if not has_gamma")
2817 insist(sum(this%binned(:, 917 ))== 0.0000000000000000 ," 802060 has zero counts if not has_gamma")
2819 insist(sum(this%binned(:, 918 ))== 0.0000000000000000 ," 812030 has zero counts if not has_gamma")
2821 insist(sum(this%binned(:, 919 ))== 0.0000000000000000 ," 812050 has zero counts if not has_gamma")
2823 insist(sum(this%binned(:, 920 ))== 0.0000000000000000 ," 812060 has zero counts if not has_gamma")
2827 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.53874462950932056 , new=",new,")"
2832 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.62573442505910604 , new=",new,")"
2837 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.74720925750205869 , new=",new,")"
2840 insist(sum(this%binned(:, 924 ))== 0.0000000000000000 ," 812100 has zero counts if not has_gamma")
2842 insist(sum(this%binned(:, 925 ))== 0.0000000000000000 ," 822030 has zero counts if not has_gamma")
2844 insist(sum(this%binned(:, 926 ))== 0.0000000000000000 ," 822040 has zero counts if not has_gamma")
2846 insist(sum(this%binned(:, 927 ))== 0.0000000000000000 ," 822050 has zero counts if not has_gamma")
2848 insist(sum(this%binned(:, 928 ))== 0.0000000000000000 ," 822060 has zero counts if not has_gamma")
2850 insist(sum(this%binned(:, 929 ))== 0.0000000000000000 ," 822070 has zero counts if not has_gamma")
2852 insist(sum(this%binned(:, 930 ))== 0.0000000000000000 ," 822071 has zero counts if not has_gamma")
2854 insist(sum(this%binned(:, 931 ))== 0.0000000000000000 ," 822080 has zero counts if not has_gamma")
2858 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.18797484500782954 , new=",new,")"
2863 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.4747836263975245E-003 , new=",new,")"
2868 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.49050796820090392 , new=",new,")"
2873 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.2801822919428225E-002 , new=",new,")"
2878 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.21343081639527520 , new=",new,")"
2881 insist(sum(this%binned(:, 937 ))== 0.0000000000000000 ," 832060 has zero counts if not has_gamma")
2883 insist(sum(this%binned(:, 938 ))== 0.0000000000000000 ," 832070 has zero counts if not has_gamma")
2885 insist(sum(this%binned(:, 939 ))== 0.0000000000000000 ," 832080 has zero counts if not has_gamma")
2887 insist(sum(this%binned(:, 940 ))== 0.0000000000000000 ," 832090 has zero counts if not has_gamma")
2889 insist(sum(this%binned(:, 941 ))== 0.0000000000000000 ," 832101 has zero counts if not has_gamma")
2893 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.41111443877903026 , new=",new,")"
2898 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.2738720412561945E-006 , new=",new,")"
2903 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.33996331021536363 , new=",new,")"
2906 insist(sum(this%binned(:, 945 ))== 0.0000000000000000 ," 832121 has zero counts if not has_gamma")
2910 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.45135284930860575 , new=",new,")"
2915 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.73373987625549841 , new=",new,")"
2918 insist(sum(this%binned(:, 948 ))== 0.0000000000000000 ," 842070 has zero counts if not has_gamma")
2920 insist(sum(this%binned(:, 949 ))== 0.0000000000000000 ," 842080 has zero counts if not has_gamma")
2922 insist(sum(this%binned(:, 950 ))== 0.0000000000000000 ," 842090 has zero counts if not has_gamma")
2924 insist(sum(this%binned(:, 951 ))== 0.0000000000000000 ," 842100 has zero counts if not has_gamma")
2926 insist(sum(this%binned(:, 952 ))== 0.0000000000000000 ," 842111 has zero counts if not has_gamma")
2928 insist(sum(this%binned(:, 953 ))== 0.0000000000000000 ," 842110 has zero counts if not has_gamma")
2930 insist(sum(this%binned(:, 954 ))== 0.0000000000000000 ," 842120 has zero counts if not has_gamma")
2932 insist(sum(this%binned(:, 955 ))== 0.0000000000000000 ," 842130 has zero counts if not has_gamma")
2934 insist(sum(this%binned(:, 956 ))== 0.0000000000000000 ," 842140 has zero counts if not has_gamma")
2936 insist(sum(this%binned(:, 957 ))== 0.0000000000000000 ," 842150 has zero counts if not has_gamma")
2938 insist(sum(this%binned(:, 958 ))== 0.0000000000000000 ," 842160 has zero counts if not has_gamma")
2940 insist(sum(this%binned(:, 959 ))== 0.0000000000000000 ," 842180 has zero counts if not has_gamma")
2942 insist(sum(this%binned(:, 960 ))== 0.0000000000000000 ," 852160 has zero counts if not has_gamma")
2944 insist(sum(this%binned(:, 961 ))== 0.0000000000000000 ," 852170 has zero counts if not has_gamma")
2946 insist(sum(this%binned(:, 962 ))== 0.0000000000000000 ," 852180 has zero counts if not has_gamma")
2948 insist(sum(this%binned(:, 963 ))== 0.0000000000000000 ," 862160 has zero counts if not has_gamma")
2950 insist(sum(this%binned(:, 964 ))== 0.0000000000000000 ," 862170 has zero counts if not has_gamma")
2952 insist(sum(this%binned(:, 965 ))== 0.0000000000000000 ," 862180 has zero counts if not has_gamma")
2954 insist(sum(this%binned(:, 966 ))== 0.0000000000000000 ," 862190 has zero counts if not has_gamma")
2956 insist(sum(this%binned(:, 967 ))== 0.0000000000000000 ," 862200 has zero counts if not has_gamma")
2958 insist(sum(this%binned(:, 968 ))== 0.0000000000000000 ," 862220 has zero counts if not has_gamma")
2960 insist(sum(this%binned(:, 969 ))== 0.0000000000000000 ," 872200 has zero counts if not has_gamma")
2962 insist(sum(this%binned(:, 970 ))== 0.0000000000000000 ," 872210 has zero counts if not has_gamma")
2964 insist(sum(this%binned(:, 971 ))== 0.0000000000000000 ," 872220 has zero counts if not has_gamma")
2968 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.34457590992681397 , new=",new,")"
2971 insist(sum(this%binned(:, 973 ))== 0.0000000000000000 ," 882200 has zero counts if not has_gamma")
2973 insist(sum(this%binned(:, 974 ))== 0.0000000000000000 ," 882220 has zero counts if not has_gamma")
2975 insist(sum(this%binned(:, 975 ))== 0.0000000000000000 ," 882230 has zero counts if not has_gamma")
2977 insist(sum(this%binned(:, 976 ))== 0.0000000000000000 ," 882240 has zero counts if not has_gamma")
2981 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 7.6982210094959369E-002 , new=",new,")"
2984 insist(sum(this%binned(:, 978 ))== 0.0000000000000000 ," 882260 has zero counts if not has_gamma")
2986 insist(sum(this%binned(:, 979 ))== 0.0000000000000000 ," 882270 has zero counts if not has_gamma")
2990 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.5048172072433814E-003 , new=",new,")"
2993 insist(sum(this%binned(:, 981 ))== 0.0000000000000000 ," 892240 has zero counts if not has_gamma")
2995 insist(sum(this%binned(:, 982 ))== 0.0000000000000000 ," 892250 has zero counts if not has_gamma")
2997 insist(sum(this%binned(:, 983 ))== 0.0000000000000000 ," 892260 has zero counts if not has_gamma")
3001 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 2.4139754975247558E-003 , new=",new,")"
3006 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.47638958783602731 , new=",new,")"
3009 insist(sum(this%binned(:, 986 ))== 0.0000000000000000 ," 902260 has zero counts if not has_gamma")
3011 insist(sum(this%binned(:, 987 ))== 0.0000000000000000 ," 902270 has zero counts if not has_gamma")
3013 insist(sum(this%binned(:, 988 ))== 0.0000000000000000 ," 902280 has zero counts if not has_gamma")
3015 insist(sum(this%binned(:, 989 ))== 0.0000000000000000 ," 902290 has zero counts if not has_gamma")
3017 insist(sum(this%binned(:, 990 ))== 0.0000000000000000 ," 902300 has zero counts if not has_gamma")
3021 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 6.1092506991736975E-002 , new=",new,")"
3024 insist(sum(this%binned(:, 992 ))== 0.0000000000000000 ," 902320 has zero counts if not has_gamma")
3026 insist(sum(this%binned(:, 993 ))== 0.0000000000000000 ," 902330 has zero counts if not has_gamma")
3030 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 2.9265930386982702E-002 , new=",new,")"
3033 insist(sum(this%binned(:, 995 ))== 0.0000000000000000 ," 912280 has zero counts if not has_gamma")
3035 insist(sum(this%binned(:, 996 ))== 0.0000000000000000 ," 912290 has zero counts if not has_gamma")
3037 insist(sum(this%binned(:, 997 ))== 0.0000000000000000 ," 912300 has zero counts if not has_gamma")
3039 insist(sum(this%binned(:, 998 ))== 0.0000000000000000 ," 912310 has zero counts if not has_gamma")
3043 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.1573892923436156E-002 , new=",new,")"
3048 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 5.2396245176693078E-002 , new=",new,")"
3053 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.95710188882963210 , new=",new,")"
3058 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.20125070867845496 , new=",new,")"
3063 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.51959527883799284 , new=",new,")"
3066 insist(sum(this%binned(:, 1004 ))== 0.0000000000000000 ," 922300 has zero counts if not has_gamma")
3068 insist(sum(this%binned(:, 1005 ))== 0.0000000000000000 ," 922310 has zero counts if not has_gamma")
3070 insist(sum(this%binned(:, 1006 ))== 0.0000000000000000 ," 922320 has zero counts if not has_gamma")
3072 insist(sum(this%binned(:, 1007 ))== 0.0000000000000000 ," 922330 has zero counts if not has_gamma")
3074 insist(sum(this%binned(:, 1008 ))== 0.0000000000000000 ," 922340 has zero counts if not has_gamma")
3076 insist(sum(this%binned(:, 1009 ))== 0.0000000000000000 ," 922350 has zero counts if not has_gamma")
3078 insist(sum(this%binned(:, 1010 ))== 0.0000000000000000 ," 922351 has zero counts if not has_gamma")
3080 insist(sum(this%binned(:, 1011 ))== 0.0000000000000000 ," 922360 has zero counts if not has_gamma")
3084 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 4.9673577145975742E-002 , new=",new,")"
3087 insist(sum(this%binned(:, 1013 ))== 0.0000000000000000 ," 922380 has zero counts if not has_gamma")
3091 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.42497242704025745 , new=",new,")"
3096 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.4512707180117816E-002 , new=",new,")"
3099 insist(sum(this%binned(:, 1016 ))== 0.0000000000000000 ," 922410 has zero counts if not has_gamma")
3101 insist(sum(this%binned(:, 1017 ))== 0.0000000000000000 ," 932340 has zero counts if not has_gamma")
3103 insist(sum(this%binned(:, 1018 ))== 0.0000000000000000 ," 932350 has zero counts if not has_gamma")
3107 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 3.2724980694425192E-002 , new=",new,")"
3112 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 6.6346084215321592E-004 , new=",new,")"
3115 insist(sum(this%binned(:, 1021 ))== 0.0000000000000000 ," 932370 has zero counts if not has_gamma")
3119 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.22831253906543125 , new=",new,")"
3124 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.10292849191012543 , new=",new,")"
3129 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.70032517157923735 , new=",new,")"
3134 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.27112565937590682 , new=",new,")"
3137 insist(sum(this%binned(:, 1026 ))== 0.0000000000000000 ," 932410 has zero counts if not has_gamma")
3139 insist(sum(this%binned(:, 1027 ))== 0.0000000000000000 ," 942360 has zero counts if not has_gamma")
3141 insist(sum(this%binned(:, 1028 ))== 0.0000000000000000 ," 942371 has zero counts if not has_gamma")
3143 insist(sum(this%binned(:, 1029 ))== 0.0000000000000000 ," 942370 has zero counts if not has_gamma")
3145 insist(sum(this%binned(:, 1030 ))== 0.0000000000000000 ," 942380 has zero counts if not has_gamma")
3147 insist(sum(this%binned(:, 1031 ))== 0.0000000000000000 ," 942390 has zero counts if not has_gamma")
3149 insist(sum(this%binned(:, 1032 ))== 0.0000000000000000 ," 942400 has zero counts if not has_gamma")
3153 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 5.8739440397093101E-004 , new=",new,")"
3156 insist(sum(this%binned(:, 1034 ))== 0.0000000000000000 ," 942420 has zero counts if not has_gamma")
3160 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.14868743975125570 , new=",new,")"
3163 insist(sum(this%binned(:, 1036 ))== 0.0000000000000000 ," 942440 has zero counts if not has_gamma")
3167 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.24108988684036203 , new=",new,")"
3170 insist(sum(this%binned(:, 1038 ))== 0.0000000000000000 ," 942460 has zero counts if not has_gamma")
3172 insist(sum(this%binned(:, 1039 ))== 0.0000000000000000 ," 942470 has zero counts if not has_gamma")
3174 insist(sum(this%binned(:, 1040 ))== 0.0000000000000000 ," 952390 has zero counts if not has_gamma")
3176 insist(sum(this%binned(:, 1041 ))== 0.0000000000000000 ," 952400 has zero counts if not has_gamma")
3178 insist(sum(this%binned(:, 1042 ))== 0.0000000000000000 ," 952410 has zero counts if not has_gamma")
3180 insist(sum(this%binned(:, 1043 ))== 0.0000000000000000 ," 952421 has zero counts if not has_gamma")
3184 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.12414116680770086 , new=",new,")"
3187 insist(sum(this%binned(:, 1045 ))== 0.0000000000000000 ," 952430 has zero counts if not has_gamma")
3191 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.55156508334189447 , new=",new,")"
3196 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 9.2915887708045389E-002 , new=",new,")"
3201 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.25973814116876875 , new=",new,")"
3206 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.38123717104017280 , new=",new,")"
3209 insist(sum(this%binned(:, 1050 ))== 0.0000000000000000 ," 952461 has zero counts if not has_gamma")
3211 insist(sum(this%binned(:, 1051 ))== 0.0000000000000000 ," 952470 has zero counts if not has_gamma")
3213 insist(sum(this%binned(:, 1052 ))== 0.0000000000000000 ," 962400 has zero counts if not has_gamma")
3215 insist(sum(this%binned(:, 1053 ))== 0.0000000000000000 ," 962410 has zero counts if not has_gamma")
3217 insist(sum(this%binned(:, 1054 ))== 0.0000000000000000 ," 962420 has zero counts if not has_gamma")
3219 insist(sum(this%binned(:, 1055 ))== 0.0000000000000000 ," 962430 has zero counts if not has_gamma")
3221 insist(sum(this%binned(:, 1056 ))== 0.0000000000000000 ," 962440 has zero counts if not has_gamma")
3223 insist(sum(this%binned(:, 1057 ))== 0.0000000000000000 ," 962450 has zero counts if not has_gamma")
3225 insist(sum(this%binned(:, 1058 ))== 0.0000000000000000 ," 962460 has zero counts if not has_gamma")
3227 insist(sum(this%binned(:, 1059 ))== 0.0000000000000000 ," 962470 has zero counts if not has_gamma")
3229 insist(sum(this%binned(:, 1060 ))== 0.0000000000000000 ," 962480 has zero counts if not has_gamma")
3233 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.27444291763657253 , new=",new,")"
3236 insist(sum(this%binned(:, 1062 ))== 0.0000000000000000 ," 962500 has zero counts if not has_gamma")
3238 insist(sum(this%binned(:, 1063 ))== 0.0000000000000000 ," 962510 has zero counts if not has_gamma")
3240 insist(sum(this%binned(:, 1064 ))== 0.0000000000000000 ," 972450 has zero counts if not has_gamma")
3242 insist(sum(this%binned(:, 1065 ))== 0.0000000000000000 ," 972460 has zero counts if not has_gamma")
3244 insist(sum(this%binned(:, 1066 ))== 0.0000000000000000 ," 972470 has zero counts if not has_gamma")
3246 insist(sum(this%binned(:, 1067 ))== 0.0000000000000000 ," 972480 has zero counts if not has_gamma")
3248 insist(sum(this%binned(:, 1068 ))== 0.0000000000000000 ," 972481 has zero counts if not has_gamma")
3252 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.8516828562347602E-002 , new=",new,")"
3257 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.26173507223435616 , new=",new,")"
3260 insist(sum(this%binned(:, 1071 ))== 0.0000000000000000 ," 972510 has zero counts if not has_gamma")
3262 insist(sum(this%binned(:, 1072 ))== 0.0000000000000000 ," 982460 has zero counts if not has_gamma")
3264 insist(sum(this%binned(:, 1073 ))== 0.0000000000000000 ," 982480 has zero counts if not has_gamma")
3266 insist(sum(this%binned(:, 1074 ))== 0.0000000000000000 ," 982490 has zero counts if not has_gamma")
3268 insist(sum(this%binned(:, 1075 ))== 0.0000000000000000 ," 982500 has zero counts if not has_gamma")
3270 insist(sum(this%binned(:, 1076 ))== 0.0000000000000000 ," 982510 has zero counts if not has_gamma")
3272 insist(sum(this%binned(:, 1077 ))== 0.0000000000000000 ," 982520 has zero counts if not has_gamma")
3276 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 6.1976144970685716E-002 , new=",new,")"
3279 insist(sum(this%binned(:, 1079 ))== 0.0000000000000000 ," 982540 has zero counts if not has_gamma")
3281 insist(sum(this%binned(:, 1080 ))== 0.0000000000000000 ," 982550 has zero counts if not has_gamma")
3283 insist(sum(this%binned(:, 1081 ))== 0.0000000000000000 ," 992510 has zero counts if not has_gamma")
3285 insist(sum(this%binned(:, 1082 ))== 0.0000000000000000 ," 992520 has zero counts if not has_gamma")
3287 insist(sum(this%binned(:, 1083 ))== 0.0000000000000000 ," 992530 has zero counts if not has_gamma")
3291 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.19549448853091214 , new=",new,")"
3294 insist(sum(this%binned(:, 1085 ))== 0.0000000000000000 ," 992540 has zero counts if not has_gamma")
3296 insist(sum(this%binned(:, 1086 ))== 0.0000000000000000 ," 992550 has zero counts if not has_gamma")
3300 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 5.8041311241133610E-004 , new=",new,")"
3303 insist(sum(this%binned(:, 1088 ))== 0.0000000000000000 ," 20030 has zero counts if not has_gamma")
3305 insist(sum(this%binned(:, 1089 ))== 0.0000000000000000 ," 20040 has zero counts if not has_gamma")
3307 insist(sum(this%binned(:, 1090 ))== 0.0000000000000000 ," 260650 has zero counts if not has_gamma")
3309 insist(sum(this%binned(:, 1091 ))== 0.0000000000000000 ," 270650 has zero counts if not has_gamma")
3313 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.71665109825983109 , new=",new,")"
3316 insist(sum(this%binned(:, 1093 ))== 0.0000000000000000 ," 290650 has zero counts if not has_gamma")
3318 insist(sum(this%binned(:, 1094 ))== 0.0000000000000000 ," 240660 has zero counts if not has_gamma")
3320 insist(sum(this%binned(:, 1095 ))== 0.0000000000000000 ," 250660 has zero counts if not has_gamma")
3322 insist(sum(this%binned(:, 1096 ))== 0.0000000000000000 ," 260660 has zero counts if not has_gamma")
3324 insist(sum(this%binned(:, 1097 ))== 0.0000000000000000 ," 270660 has zero counts if not has_gamma")
3328 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 4.9393710887042759E-002 , new=",new,")"
3336 insist(sum(this%binned(:, 1100 ))== 0.0000000000000000 ," 300660 has zero counts if not has_gamma")
3338 insist(sum(this%binned(:, 1101 ))== 0.0000000000000000 ," 310660 has zero counts if not has_gamma")
3340 insist(sum(this%binned(:, 1102 ))== 0.0000000000000000 ," 320660 has zero counts if not has_gamma")
3342 insist(sum(this%binned(:, 1103 ))== 0.0000000000000000 ," 240670 has zero counts if not has_gamma")
3344 insist(sum(this%binned(:, 1104 ))== 0.0000000000000000 ," 250670 has zero counts if not has_gamma")
3346 insist(sum(this%binned(:, 1105 ))== 0.0000000000000000 ," 260670 has zero counts if not has_gamma")
3348 insist(sum(this%binned(:, 1106 ))== 0.0000000000000000 ," 270670 has zero counts if not has_gamma")
3350 insist(sum(this%binned(:, 1107 ))== 0.0000000000000000 ," 280670 has zero counts if not has_gamma")
3354 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.12419624399658105 , new=",new,")"
3357 insist(sum(this%binned(:, 1109 ))== 0.0000000000000000 ," 300670 has zero counts if not has_gamma")
3359 insist(sum(this%binned(:, 1110 ))== 0.0000000000000000 ," 310670 has zero counts if not has_gamma")
3361 insist(sum(this%binned(:, 1111 ))== 0.0000000000000000 ," 320670 has zero counts if not has_gamma")
3363 insist(sum(this%binned(:, 1112 ))== 0.0000000000000000 ," 250680 has zero counts if not has_gamma")
3365 insist(sum(this%binned(:, 1113 ))== 0.0000000000000000 ," 260680 has zero counts if not has_gamma")
3367 insist(sum(this%binned(:, 1114 ))== 0.0000000000000000 ," 270680 has zero counts if not has_gamma")
3369 insist(sum(this%binned(:, 1115 ))== 0.0000000000000000 ," 280680 has zero counts if not has_gamma")
3371 insist(sum(this%binned(:, 1116 ))== 0.0000000000000000 ," 290680 has zero counts if not has_gamma")
3373 insist(sum(this%binned(:, 1117 ))== 0.0000000000000000 ," 290681 has zero counts if not has_gamma")
3375 insist(sum(this%binned(:, 1118 ))== 0.0000000000000000 ," 300680 has zero counts if not has_gamma")
3377 insist(sum(this%binned(:, 1119 ))== 0.0000000000000000 ," 310680 has zero counts if not has_gamma")
3379 insist(sum(this%binned(:, 1120 ))== 0.0000000000000000 ," 320680 has zero counts if not has_gamma")
3381 insist(sum(this%binned(:, 1121 ))== 0.0000000000000000 ," 250690 has zero counts if not has_gamma")
3383 insist(sum(this%binned(:, 1122 ))== 0.0000000000000000 ," 260690 has zero counts if not has_gamma")
3385 insist(sum(this%binned(:, 1123 ))== 0.0000000000000000 ," 270690 has zero counts if not has_gamma")
3387 insist(sum(this%binned(:, 1124 ))== 0.0000000000000000 ," 280690 has zero counts if not has_gamma")
3389 insist(sum(this%binned(:, 1125 ))== 0.0000000000000000 ," 290690 has zero counts if not has_gamma")
3393 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.32448997333713686 , new=",new,")"
3398 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 3.1907135371456243E-008 , new=",new,")"
3401 insist(sum(this%binned(:, 1128 ))== 0.0000000000000000 ," 310690 has zero counts if not has_gamma")
3403 insist(sum(this%binned(:, 1129 ))== 0.0000000000000000 ," 320690 has zero counts if not has_gamma")
3405 insist(sum(this%binned(:, 1130 ))== 0.0000000000000000 ," 330690 has zero counts if not has_gamma")
3407 insist(sum(this%binned(:, 1131 ))== 0.0000000000000000 ," 260700 has zero counts if not has_gamma")
3409 insist(sum(this%binned(:, 1132 ))== 0.0000000000000000 ," 270700 has zero counts if not has_gamma")
3411 insist(sum(this%binned(:, 1133 ))== 0.0000000000000000 ," 280700 has zero counts if not has_gamma")
3413 insist(sum(this%binned(:, 1134 ))== 0.0000000000000000 ," 290700 has zero counts if not has_gamma")
3415 insist(sum(this%binned(:, 1135 ))== 0.0000000000000000 ," 290701 has zero counts if not has_gamma")
3417 insist(sum(this%binned(:, 1136 ))== 0.0000000000000000 ," 300700 has zero counts if not has_gamma")
3421 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.72407405871151376 , new=",new,")"
3424 insist(sum(this%binned(:, 1138 ))== 0.0000000000000000 ," 320700 has zero counts if not has_gamma")
3426 insist(sum(this%binned(:, 1139 ))== 0.0000000000000000 ," 260710 has zero counts if not has_gamma")
3428 insist(sum(this%binned(:, 1140 ))== 0.0000000000000000 ," 270710 has zero counts if not has_gamma")
3430 insist(sum(this%binned(:, 1141 ))== 0.0000000000000000 ," 280710 has zero counts if not has_gamma")
3432 insist(sum(this%binned(:, 1142 ))== 0.0000000000000000 ," 290710 has zero counts if not has_gamma")
3441 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.61465169864222335 , new=",new,")"
3444 insist(sum(this%binned(:, 1145 ))== 0.0000000000000000 ," 310710 has zero counts if not has_gamma")
3446 insist(sum(this%binned(:, 1146 ))== 0.0000000000000000 ," 320710 has zero counts if not has_gamma")
3448 insist(sum(this%binned(:, 1147 ))== 0.0000000000000000 ," 320711 has zero counts if not has_gamma")
3450 insist(sum(this%binned(:, 1148 ))== 0.0000000000000000 ," 330710 has zero counts if not has_gamma")
3452 insist(sum(this%binned(:, 1149 ))== 0.0000000000000000 ," 260720 has zero counts if not has_gamma")
3454 insist(sum(this%binned(:, 1150 ))== 0.0000000000000000 ," 270720 has zero counts if not has_gamma")
3456 insist(sum(this%binned(:, 1151 ))== 0.0000000000000000 ," 280720 has zero counts if not has_gamma")
3458 insist(sum(this%binned(:, 1152 ))== 0.0000000000000000 ," 290720 has zero counts if not has_gamma")
3462 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 6.7099930788272200E-002 , new=",new,")"
3467 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.54910348605116188 , new=",new,")"
3470 insist(sum(this%binned(:, 1155 ))== 0.0000000000000000 ," 310721 has zero counts if not has_gamma")
3472 insist(sum(this%binned(:, 1156 ))== 0.0000000000000000 ," 320720 has zero counts if not has_gamma")
3474 insist(sum(this%binned(:, 1157 ))== 0.0000000000000000 ," 330720 has zero counts if not has_gamma")
3476 insist(sum(this%binned(:, 1158 ))== 0.0000000000000000 ," 340720 has zero counts if not has_gamma")
3478 insist(sum(this%binned(:, 1159 ))== 0.0000000000000000 ," 270730 has zero counts if not has_gamma")
3480 insist(sum(this%binned(:, 1160 ))== 0.0000000000000000 ," 280730 has zero counts if not has_gamma")
3482 insist(sum(this%binned(:, 1161 ))== 0.0000000000000000 ," 290730 has zero counts if not has_gamma")
3484 insist(sum(this%binned(:, 1162 ))== 0.0000000000000000 ," 300730 has zero counts if not has_gamma")
3486 insist(sum(this%binned(:, 1163 ))== 0.0000000000000000 ," 310730 has zero counts if not has_gamma")
3488 insist(sum(this%binned(:, 1164 ))== 0.0000000000000000 ," 320730 has zero counts if not has_gamma")
3490 insist(sum(this%binned(:, 1165 ))== 0.0000000000000000 ," 320731 has zero counts if not has_gamma")
3492 insist(sum(this%binned(:, 1166 ))== 0.0000000000000000 ," 330730 has zero counts if not has_gamma")
3494 insist(sum(this%binned(:, 1167 ))== 0.0000000000000000 ," 340730 has zero counts if not has_gamma")
3496 insist(sum(this%binned(:, 1168 ))== 0.0000000000000000 ," 340731 has zero counts if not has_gamma")
3498 insist(sum(this%binned(:, 1169 ))== 0.0000000000000000 ," 270740 has zero counts if not has_gamma")
3500 insist(sum(this%binned(:, 1170 ))== 0.0000000000000000 ," 280740 has zero counts if not has_gamma")
3502 insist(sum(this%binned(:, 1171 ))== 0.0000000000000000 ," 290740 has zero counts if not has_gamma")
3514 insist(sum(this%binned(:, 1174 ))== 0.0000000000000000 ," 310741 has zero counts if not has_gamma")
3516 insist(sum(this%binned(:, 1175 ))== 0.0000000000000000 ," 320740 has zero counts if not has_gamma")
3518 insist(sum(this%binned(:, 1176 ))== 0.0000000000000000 ," 330740 has zero counts if not has_gamma")
3520 insist(sum(this%binned(:, 1177 ))== 0.0000000000000000 ," 340740 has zero counts if not has_gamma")
3522 insist(sum(this%binned(:, 1178 ))== 0.0000000000000000 ," 270750 has zero counts if not has_gamma")
3524 insist(sum(this%binned(:, 1179 ))== 0.0000000000000000 ," 280750 has zero counts if not has_gamma")
3526 insist(sum(this%binned(:, 1180 ))== 0.0000000000000000 ," 290750 has zero counts if not has_gamma")
3528 insist(sum(this%binned(:, 1181 ))== 0.0000000000000000 ," 300750 has zero counts if not has_gamma")
3530 insist(sum(this%binned(:, 1182 ))== 0.0000000000000000 ," 310750 has zero counts if not has_gamma")
3534 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.44207167789574253 , new=",new,")"
3537 insist(sum(this%binned(:, 1184 ))== 0.0000000000000000 ," 320751 has zero counts if not has_gamma")
3539 insist(sum(this%binned(:, 1185 ))== 0.0000000000000000 ," 330750 has zero counts if not has_gamma")
3541 insist(sum(this%binned(:, 1186 ))== 0.0000000000000000 ," 330751 has zero counts if not has_gamma")
3543 insist(sum(this%binned(:, 1187 ))== 0.0000000000000000 ," 340750 has zero counts if not has_gamma")
3545 insist(sum(this%binned(:, 1188 ))== 0.0000000000000000 ," 350750 has zero counts if not has_gamma")
3547 insist(sum(this%binned(:, 1189 ))== 0.0000000000000000 ," 280760 has zero counts if not has_gamma")
3549 insist(sum(this%binned(:, 1190 ))== 0.0000000000000000 ," 290760 has zero counts if not has_gamma")
3551 insist(sum(this%binned(:, 1191 ))== 0.0000000000000000 ," 300760 has zero counts if not has_gamma")
3558 insist(sum(this%binned(:, 1193 ))== 0.0000000000000000 ," 320760 has zero counts if not has_gamma")
3565 insist(sum(this%binned(:, 1195 ))== 0.0000000000000000 ," 340760 has zero counts if not has_gamma")
3567 insist(sum(this%binned(:, 1196 ))== 0.0000000000000000 ," 280770 has zero counts if not has_gamma")
3569 insist(sum(this%binned(:, 1197 ))== 0.0000000000000000 ," 290770 has zero counts if not has_gamma")
3571 insist(sum(this%binned(:, 1198 ))== 0.0000000000000000 ," 300770 has zero counts if not has_gamma")
3573 insist(sum(this%binned(:, 1199 ))== 0.0000000000000000 ," 310770 has zero counts if not has_gamma")
3577 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.73257897902610103 , new=",new,")"
3587 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.21992764589933780 , new=",new,")"
3590 insist(sum(this%binned(:, 1203 ))== 0.0000000000000000 ," 340770 has zero counts if not has_gamma")
3592 insist(sum(this%binned(:, 1204 ))== 0.0000000000000000 ," 340771 has zero counts if not has_gamma")
3594 insist(sum(this%binned(:, 1205 ))== 0.0000000000000000 ," 350770 has zero counts if not has_gamma")
3596 insist(sum(this%binned(:, 1206 ))== 0.0000000000000000 ," 350771 has zero counts if not has_gamma")
3598 insist(sum(this%binned(:, 1207 ))== 0.0000000000000000 ," 360770 has zero counts if not has_gamma")
3600 insist(sum(this%binned(:, 1208 ))== 0.0000000000000000 ," 280780 has zero counts if not has_gamma")
3602 insist(sum(this%binned(:, 1209 ))== 0.0000000000000000 ," 290780 has zero counts if not has_gamma")
3604 insist(sum(this%binned(:, 1210 ))== 0.0000000000000000 ," 300780 has zero counts if not has_gamma")
3606 insist(sum(this%binned(:, 1211 ))== 0.0000000000000000 ," 310780 has zero counts if not has_gamma")
3610 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.22897752262353060 , new=",new,")"
3618 insist(sum(this%binned(:, 1214 ))== 0.0000000000000000 ," 340780 has zero counts if not has_gamma")
3620 insist(sum(this%binned(:, 1215 ))== 0.0000000000000000 ," 350780 has zero counts if not has_gamma")
3622 insist(sum(this%binned(:, 1216 ))== 0.0000000000000000 ," 360780 has zero counts if not has_gamma")
3624 insist(sum(this%binned(:, 1217 ))== 0.0000000000000000 ," 290790 has zero counts if not has_gamma")
3626 insist(sum(this%binned(:, 1218 ))== 0.0000000000000000 ," 300790 has zero counts if not has_gamma")
3628 insist(sum(this%binned(:, 1219 ))== 0.0000000000000000 ," 310790 has zero counts if not has_gamma")
3635 insist(sum(this%binned(:, 1221 ))== 0.0000000000000000 ," 320791 has zero counts if not has_gamma")
3639 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.98792303115040181 , new=",new,")"
3644 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 3.6759974434876781E-002 , new=",new,")"
3647 insist(sum(this%binned(:, 1224 ))== 0.0000000000000000 ," 340791 has zero counts if not has_gamma")
3649 insist(sum(this%binned(:, 1225 ))== 0.0000000000000000 ," 350790 has zero counts if not has_gamma")
3651 insist(sum(this%binned(:, 1226 ))== 0.0000000000000000 ," 350791 has zero counts if not has_gamma")
3655 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 5.5476567838820726E-002 , new=",new,")"
3658 insist(sum(this%binned(:, 1228 ))== 0.0000000000000000 ," 360791 has zero counts if not has_gamma")
3660 insist(sum(this%binned(:, 1229 ))== 0.0000000000000000 ," 370790 has zero counts if not has_gamma")
3662 insist(sum(this%binned(:, 1230 ))== 0.0000000000000000 ," 290800 has zero counts if not has_gamma")
3664 insist(sum(this%binned(:, 1231 ))== 0.0000000000000000 ," 300800 has zero counts if not has_gamma")
3666 insist(sum(this%binned(:, 1232 ))== 0.0000000000000000 ," 310800 has zero counts if not has_gamma")
3668 insist(sum(this%binned(:, 1233 ))== 0.0000000000000000 ," 320800 has zero counts if not has_gamma")
3675 insist(sum(this%binned(:, 1235 ))== 0.0000000000000000 ," 340800 has zero counts if not has_gamma")
3679 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.77559239318698625 , new=",new,")"
3682 insist(sum(this%binned(:, 1237 ))== 0.0000000000000000 ," 350801 has zero counts if not has_gamma")
3684 insist(sum(this%binned(:, 1238 ))== 0.0000000000000000 ," 360800 has zero counts if not has_gamma")
3686 insist(sum(this%binned(:, 1239 ))== 0.0000000000000000 ," 300810 has zero counts if not has_gamma")
3688 insist(sum(this%binned(:, 1240 ))== 0.0000000000000000 ," 310810 has zero counts if not has_gamma")
3690 insist(sum(this%binned(:, 1241 ))== 0.0000000000000000 ," 320810 has zero counts if not has_gamma")
3692 insist(sum(this%binned(:, 1242 ))== 0.0000000000000000 ," 320811 has zero counts if not has_gamma")
3694 insist(sum(this%binned(:, 1243 ))== 0.0000000000000000 ," 330810 has zero counts if not has_gamma")
3696 insist(sum(this%binned(:, 1244 ))== 0.0000000000000000 ," 340810 has zero counts if not has_gamma")
3698 insist(sum(this%binned(:, 1245 ))== 0.0000000000000000 ," 340811 has zero counts if not has_gamma")
3700 insist(sum(this%binned(:, 1246 ))== 0.0000000000000000 ," 350810 has zero counts if not has_gamma")
3702 insist(sum(this%binned(:, 1247 ))== 0.0000000000000000 ," 360810 has zero counts if not has_gamma")
3704 insist(sum(this%binned(:, 1248 ))== 0.0000000000000000 ," 360811 has zero counts if not has_gamma")
3706 insist(sum(this%binned(:, 1249 ))== 0.0000000000000000 ," 370810 has zero counts if not has_gamma")
3708 insist(sum(this%binned(:, 1250 ))== 0.0000000000000000 ," 300820 has zero counts if not has_gamma")
3710 insist(sum(this%binned(:, 1251 ))== 0.0000000000000000 ," 310820 has zero counts if not has_gamma")
3712 insist(sum(this%binned(:, 1252 ))== 0.0000000000000000 ," 320820 has zero counts if not has_gamma")
3724 insist(sum(this%binned(:, 1255 ))== 0.0000000000000000 ," 340820 has zero counts if not has_gamma")
3728 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.12020021954554848 , new=",new,")"
3733 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 9.4313711645749483E-004 , new=",new,")"
3736 insist(sum(this%binned(:, 1258 ))== 0.0000000000000000 ," 360820 has zero counts if not has_gamma")
3738 insist(sum(this%binned(:, 1259 ))== 0.0000000000000000 ," 300830 has zero counts if not has_gamma")
3740 insist(sum(this%binned(:, 1260 ))== 0.0000000000000000 ," 310830 has zero counts if not has_gamma")
3742 insist(sum(this%binned(:, 1261 ))== 0.0000000000000000 ," 320830 has zero counts if not has_gamma")
3744 insist(sum(this%binned(:, 1262 ))== 0.0000000000000000 ," 330830 has zero counts if not has_gamma")
3748 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.60150839844070703 , new=",new,")"
3758 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.32627888912931968 , new=",new,")"
3761 insist(sum(this%binned(:, 1266 ))== 0.0000000000000000 ," 360830 has zero counts if not has_gamma")
3763 insist(sum(this%binned(:, 1267 ))== 0.0000000000000000 ," 360831 has zero counts if not has_gamma")
3765 insist(sum(this%binned(:, 1268 ))== 0.0000000000000000 ," 370830 has zero counts if not has_gamma")
3767 insist(sum(this%binned(:, 1269 ))== 0.0000000000000000 ," 380830 has zero counts if not has_gamma")
3769 insist(sum(this%binned(:, 1270 ))== 0.0000000000000000 ," 310840 has zero counts if not has_gamma")
3771 insist(sum(this%binned(:, 1271 ))== 0.0000000000000000 ," 320840 has zero counts if not has_gamma")
3773 insist(sum(this%binned(:, 1272 ))== 0.0000000000000000 ," 330840 has zero counts if not has_gamma")
3775 insist(sum(this%binned(:, 1273 ))== 0.0000000000000000 ," 340840 has zero counts if not has_gamma")
3782 insist(sum(this%binned(:, 1275 ))== 0.0000000000000000 ," 350841 has zero counts if not has_gamma")
3784 insist(sum(this%binned(:, 1276 ))== 0.0000000000000000 ," 360840 has zero counts if not has_gamma")
3786 insist(sum(this%binned(:, 1277 ))== 0.0000000000000000 ," 370840 has zero counts if not has_gamma")
3788 insist(sum(this%binned(:, 1278 ))== 0.0000000000000000 ," 380840 has zero counts if not has_gamma")
3790 insist(sum(this%binned(:, 1279 ))== 0.0000000000000000 ," 310850 has zero counts if not has_gamma")
3792 insist(sum(this%binned(:, 1280 ))== 0.0000000000000000 ," 320850 has zero counts if not has_gamma")
3794 insist(sum(this%binned(:, 1281 ))== 0.0000000000000000 ," 330850 has zero counts if not has_gamma")
3796 insist(sum(this%binned(:, 1282 ))== 0.0000000000000000 ," 340850 has zero counts if not has_gamma")
3805 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.24726424389526619 , new=",new,")"
3810 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.18059291069431535 , new=",new,")"
3813 insist(sum(this%binned(:, 1286 ))== 0.0000000000000000 ," 370850 has zero counts if not has_gamma")
3815 insist(sum(this%binned(:, 1287 ))== 0.0000000000000000 ," 380850 has zero counts if not has_gamma")
3817 insist(sum(this%binned(:, 1288 ))== 0.0000000000000000 ," 380851 has zero counts if not has_gamma")
3819 insist(sum(this%binned(:, 1289 ))== 0.0000000000000000 ," 390850 has zero counts if not has_gamma")
3821 insist(sum(this%binned(:, 1290 ))== 0.0000000000000000 ," 310860 has zero counts if not has_gamma")
3823 insist(sum(this%binned(:, 1291 ))== 0.0000000000000000 ," 320860 has zero counts if not has_gamma")
3825 insist(sum(this%binned(:, 1292 ))== 0.0000000000000000 ," 330860 has zero counts if not has_gamma")
3827 insist(sum(this%binned(:, 1293 ))== 0.0000000000000000 ," 340860 has zero counts if not has_gamma")
3829 insist(sum(this%binned(:, 1294 ))== 0.0000000000000000 ," 350860 has zero counts if not has_gamma")
3831 insist(sum(this%binned(:, 1295 ))== 0.0000000000000000 ," 360860 has zero counts if not has_gamma")
3835 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.76634182647394378 , new=",new,")"
3838 insist(sum(this%binned(:, 1297 ))== 0.0000000000000000 ," 370861 has zero counts if not has_gamma")
3840 insist(sum(this%binned(:, 1298 ))== 0.0000000000000000 ," 380860 has zero counts if not has_gamma")
3842 insist(sum(this%binned(:, 1299 ))== 0.0000000000000000 ," 320870 has zero counts if not has_gamma")
3844 insist(sum(this%binned(:, 1300 ))== 0.0000000000000000 ," 330870 has zero counts if not has_gamma")
3846 insist(sum(this%binned(:, 1301 ))== 0.0000000000000000 ," 340870 has zero counts if not has_gamma")
3848 insist(sum(this%binned(:, 1302 ))== 0.0000000000000000 ," 350870 has zero counts if not has_gamma")
3857 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 6.0877798994170014E-002 , new=",new,")"
3860 insist(sum(this%binned(:, 1305 ))== 0.0000000000000000 ," 380870 has zero counts if not has_gamma")
3862 insist(sum(this%binned(:, 1306 ))== 0.0000000000000000 ," 380871 has zero counts if not has_gamma")
3864 insist(sum(this%binned(:, 1307 ))== 0.0000000000000000 ," 390870 has zero counts if not has_gamma")
3866 insist(sum(this%binned(:, 1308 ))== 0.0000000000000000 ," 390871 has zero counts if not has_gamma")
3868 insist(sum(this%binned(:, 1309 ))== 0.0000000000000000 ," 400870 has zero counts if not has_gamma")
3870 insist(sum(this%binned(:, 1310 ))== 0.0000000000000000 ," 320880 has zero counts if not has_gamma")
3872 insist(sum(this%binned(:, 1311 ))== 0.0000000000000000 ," 330880 has zero counts if not has_gamma")
3874 insist(sum(this%binned(:, 1312 ))== 0.0000000000000000 ," 340880 has zero counts if not has_gamma")
3883 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.39692490516986828 , new=",new,")"
3891 insist(sum(this%binned(:, 1316 ))== 0.0000000000000000 ," 380880 has zero counts if not has_gamma")
3893 insist(sum(this%binned(:, 1317 ))== 0.0000000000000000 ," 390880 has zero counts if not has_gamma")
3895 insist(sum(this%binned(:, 1318 ))== 0.0000000000000000 ," 400880 has zero counts if not has_gamma")
3897 insist(sum(this%binned(:, 1319 ))== 0.0000000000000000 ," 320890 has zero counts if not has_gamma")
3899 insist(sum(this%binned(:, 1320 ))== 0.0000000000000000 ," 330890 has zero counts if not has_gamma")
3901 insist(sum(this%binned(:, 1321 ))== 0.0000000000000000 ," 340890 has zero counts if not has_gamma")
3903 insist(sum(this%binned(:, 1322 ))== 0.0000000000000000 ," 350890 has zero counts if not has_gamma")
3917 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.65625788511815231 , new=",new,")"
3920 insist(sum(this%binned(:, 1326 ))== 0.0000000000000000 ," 390890 has zero counts if not has_gamma")
3922 insist(sum(this%binned(:, 1327 ))== 0.0000000000000000 ," 390891 has zero counts if not has_gamma")
3926 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.21884140142702282 , new=",new,")"
3929 insist(sum(this%binned(:, 1329 ))== 0.0000000000000000 ," 400891 has zero counts if not has_gamma")
3931 insist(sum(this%binned(:, 1330 ))== 0.0000000000000000 ," 410890 has zero counts if not has_gamma")
3933 insist(sum(this%binned(:, 1331 ))== 0.0000000000000000 ," 330900 has zero counts if not has_gamma")
3935 insist(sum(this%binned(:, 1332 ))== 0.0000000000000000 ," 340900 has zero counts if not has_gamma")
3937 insist(sum(this%binned(:, 1333 ))== 0.0000000000000000 ," 350900 has zero counts if not has_gamma")
3956 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.18447088309130777 , new=",new,")"
3966 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.4444724824302237E-009 , new=",new,")"
3969 insist(sum(this%binned(:, 1340 ))== 0.0000000000000000 ," 400900 has zero counts if not has_gamma")
3971 insist(sum(this%binned(:, 1341 ))== 0.0000000000000000 ," 400901 has zero counts if not has_gamma")
3973 insist(sum(this%binned(:, 1342 ))== 0.0000000000000000 ," 410900 has zero counts if not has_gamma")
3975 insist(sum(this%binned(:, 1343 ))== 0.0000000000000000 ," 410901 has zero counts if not has_gamma")
3977 insist(sum(this%binned(:, 1344 ))== 0.0000000000000000 ," 420900 has zero counts if not has_gamma")
3979 insist(sum(this%binned(:, 1345 ))== 0.0000000000000000 ," 330910 has zero counts if not has_gamma")
3981 insist(sum(this%binned(:, 1346 ))== 0.0000000000000000 ," 340910 has zero counts if not has_gamma")
3983 insist(sum(this%binned(:, 1347 ))== 0.0000000000000000 ," 350910 has zero counts if not has_gamma")
3985 insist(sum(this%binned(:, 1348 ))== 0.0000000000000000 ," 360910 has zero counts if not has_gamma")
3987 insist(sum(this%binned(:, 1349 ))== 0.0000000000000000 ," 370910 has zero counts if not has_gamma")
3991 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.75223841715706297 , new=",new,")"
3996 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.68106379883573032 , new=",new,")"
3999 insist(sum(this%binned(:, 1352 ))== 0.0000000000000000 ," 390911 has zero counts if not has_gamma")
4001 insist(sum(this%binned(:, 1353 ))== 0.0000000000000000 ," 400910 has zero counts if not has_gamma")
4003 insist(sum(this%binned(:, 1354 ))== 0.0000000000000000 ," 410910 has zero counts if not has_gamma")
4005 insist(sum(this%binned(:, 1355 ))== 0.0000000000000000 ," 410911 has zero counts if not has_gamma")
4007 insist(sum(this%binned(:, 1356 ))== 0.0000000000000000 ," 420910 has zero counts if not has_gamma")
4009 insist(sum(this%binned(:, 1357 ))== 0.0000000000000000 ," 330920 has zero counts if not has_gamma")
4011 insist(sum(this%binned(:, 1358 ))== 0.0000000000000000 ," 340920 has zero counts if not has_gamma")
4013 insist(sum(this%binned(:, 1359 ))== 0.0000000000000000 ," 350920 has zero counts if not has_gamma")
4027 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.18683929832809518 , new=",new,")"
4035 insist(sum(this%binned(:, 1364 ))== 0.0000000000000000 ," 400920 has zero counts if not has_gamma")
4037 insist(sum(this%binned(:, 1365 ))== 0.0000000000000000 ," 410920 has zero counts if not has_gamma")
4039 insist(sum(this%binned(:, 1366 ))== 0.0000000000000000 ," 410921 has zero counts if not has_gamma")
4041 insist(sum(this%binned(:, 1367 ))== 0.0000000000000000 ," 420920 has zero counts if not has_gamma")
4043 insist(sum(this%binned(:, 1368 ))== 0.0000000000000000 ," 340930 has zero counts if not has_gamma")
4045 insist(sum(this%binned(:, 1369 ))== 0.0000000000000000 ," 350930 has zero counts if not has_gamma")
4047 insist(sum(this%binned(:, 1370 ))== 0.0000000000000000 ," 360930 has zero counts if not has_gamma")
4049 insist(sum(this%binned(:, 1371 ))== 0.0000000000000000 ," 370930 has zero counts if not has_gamma")
4051 insist(sum(this%binned(:, 1372 ))== 0.0000000000000000 ," 380930 has zero counts if not has_gamma")
4058 insist(sum(this%binned(:, 1374 ))== 0.0000000000000000 ," 390931 has zero counts if not has_gamma")
4062 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.0673408209577246E-003 , new=",new,")"
4065 insist(sum(this%binned(:, 1376 ))== 0.0000000000000000 ," 410930 has zero counts if not has_gamma")
4067 insist(sum(this%binned(:, 1377 ))== 0.0000000000000000 ," 410931 has zero counts if not has_gamma")
4069 insist(sum(this%binned(:, 1378 ))== 0.0000000000000000 ," 420930 has zero counts if not has_gamma")
4071 insist(sum(this%binned(:, 1379 ))== 0.0000000000000000 ," 420931 has zero counts if not has_gamma")
4073 insist(sum(this%binned(:, 1380 ))== 0.0000000000000000 ," 430930 has zero counts if not has_gamma")
4075 insist(sum(this%binned(:, 1381 ))== 0.0000000000000000 ," 340940 has zero counts if not has_gamma")
4077 insist(sum(this%binned(:, 1382 ))== 0.0000000000000000 ," 350940 has zero counts if not has_gamma")
4079 insist(sum(this%binned(:, 1383 ))== 0.0000000000000000 ," 360940 has zero counts if not has_gamma")
4086 insist(sum(this%binned(:, 1385 ))== 0.0000000000000000 ," 380940 has zero counts if not has_gamma")
4093 insist(sum(this%binned(:, 1387 ))== 0.0000000000000000 ," 400940 has zero counts if not has_gamma")
4097 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.12969260956754219 , new=",new,")"
4102 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.0912115440728130E-005 , new=",new,")"
4105 insist(sum(this%binned(:, 1390 ))== 0.0000000000000000 ," 420940 has zero counts if not has_gamma")
4107 insist(sum(this%binned(:, 1391 ))== 0.0000000000000000 ," 350950 has zero counts if not has_gamma")
4109 insist(sum(this%binned(:, 1392 ))== 0.0000000000000000 ," 360950 has zero counts if not has_gamma")
4111 insist(sum(this%binned(:, 1393 ))== 0.0000000000000000 ," 370950 has zero counts if not has_gamma")
4113 insist(sum(this%binned(:, 1394 ))== 0.0000000000000000 ," 380950 has zero counts if not has_gamma")
4115 insist(sum(this%binned(:, 1395 ))== 0.0000000000000000 ," 390950 has zero counts if not has_gamma")
4119 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 9.8333829532349815E-002 , new=",new,")"
4124 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 2.7653486272365409E-002 , new=",new,")"
4127 insist(sum(this%binned(:, 1398 ))== 0.0000000000000000 ," 410951 has zero counts if not has_gamma")
4129 insist(sum(this%binned(:, 1399 ))== 0.0000000000000000 ," 420950 has zero counts if not has_gamma")
4131 insist(sum(this%binned(:, 1400 ))== 0.0000000000000000 ," 430950 has zero counts if not has_gamma")
4133 insist(sum(this%binned(:, 1401 ))== 0.0000000000000000 ," 430951 has zero counts if not has_gamma")
4135 insist(sum(this%binned(:, 1402 ))== 0.0000000000000000 ," 440950 has zero counts if not has_gamma")
4137 insist(sum(this%binned(:, 1403 ))== 0.0000000000000000 ," 350960 has zero counts if not has_gamma")
4139 insist(sum(this%binned(:, 1404 ))== 0.0000000000000000 ," 360960 has zero counts if not has_gamma")
4141 insist(sum(this%binned(:, 1405 ))== 0.0000000000000000 ," 370960 has zero counts if not has_gamma")
4143 insist(sum(this%binned(:, 1406 ))== 0.0000000000000000 ," 380960 has zero counts if not has_gamma")
4150 insist(sum(this%binned(:, 1408 ))== 0.0000000000000000 ," 390961 has zero counts if not has_gamma")
4152 insist(sum(this%binned(:, 1409 ))== 0.0000000000000000 ," 400960 has zero counts if not has_gamma")
4156 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.24353126187349969 , new=",new,")"
4159 insist(sum(this%binned(:, 1411 ))== 0.0000000000000000 ," 420960 has zero counts if not has_gamma")
4161 insist(sum(this%binned(:, 1412 ))== 0.0000000000000000 ," 430960 has zero counts if not has_gamma")
4163 insist(sum(this%binned(:, 1413 ))== 0.0000000000000000 ," 440960 has zero counts if not has_gamma")
4165 insist(sum(this%binned(:, 1414 ))== 0.0000000000000000 ," 350970 has zero counts if not has_gamma")
4167 insist(sum(this%binned(:, 1415 ))== 0.0000000000000000 ," 360970 has zero counts if not has_gamma")
4169 insist(sum(this%binned(:, 1416 ))== 0.0000000000000000 ," 370970 has zero counts if not has_gamma")
4171 insist(sum(this%binned(:, 1417 ))== 0.0000000000000000 ," 380970 has zero counts if not has_gamma")
4173 insist(sum(this%binned(:, 1418 ))== 0.0000000000000000 ," 390970 has zero counts if not has_gamma")
4175 insist(sum(this%binned(:, 1419 ))== 0.0000000000000000 ," 390971 has zero counts if not has_gamma")
4179 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.79223635100248524 , new=",new,")"
4184 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.50117733078528348 , new=",new,")"
4187 insist(sum(this%binned(:, 1422 ))== 0.0000000000000000 ," 410971 has zero counts if not has_gamma")
4189 insist(sum(this%binned(:, 1423 ))== 0.0000000000000000 ," 420970 has zero counts if not has_gamma")
4191 insist(sum(this%binned(:, 1424 ))== 0.0000000000000000 ," 430970 has zero counts if not has_gamma")
4193 insist(sum(this%binned(:, 1425 ))== 0.0000000000000000 ," 430971 has zero counts if not has_gamma")
4195 insist(sum(this%binned(:, 1426 ))== 0.0000000000000000 ," 440970 has zero counts if not has_gamma")
4197 insist(sum(this%binned(:, 1427 ))== 0.0000000000000000 ," 360980 has zero counts if not has_gamma")
4199 insist(sum(this%binned(:, 1428 ))== 0.0000000000000000 ," 370980 has zero counts if not has_gamma")
4201 insist(sum(this%binned(:, 1429 ))== 0.0000000000000000 ," 380980 has zero counts if not has_gamma")
4203 insist(sum(this%binned(:, 1430 ))== 0.0000000000000000 ," 390980 has zero counts if not has_gamma")
4205 insist(sum(this%binned(:, 1431 ))== 0.0000000000000000 ," 390981 has zero counts if not has_gamma")
4207 insist(sum(this%binned(:, 1432 ))== 0.0000000000000000 ," 400980 has zero counts if not has_gamma")
4219 insist(sum(this%binned(:, 1435 ))== 0.0000000000000000 ," 420980 has zero counts if not has_gamma")
4223 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.10194629045810190 , new=",new,")"
4226 insist(sum(this%binned(:, 1437 ))== 0.0000000000000000 ," 440980 has zero counts if not has_gamma")
4228 insist(sum(this%binned(:, 1438 ))== 0.0000000000000000 ," 360990 has zero counts if not has_gamma")
4230 insist(sum(this%binned(:, 1439 ))== 0.0000000000000000 ," 370990 has zero counts if not has_gamma")
4232 insist(sum(this%binned(:, 1440 ))== 0.0000000000000000 ," 380990 has zero counts if not has_gamma")
4234 insist(sum(this%binned(:, 1441 ))== 0.0000000000000000 ," 390990 has zero counts if not has_gamma")
4236 insist(sum(this%binned(:, 1442 ))== 0.0000000000000000 ," 400990 has zero counts if not has_gamma")
4250 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.40968073719541048 , new=",new,")"
4255 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 6.6673559170234667E-002 , new=",new,")"
4258 insist(sum(this%binned(:, 1447 ))== 0.0000000000000000 ," 430991 has zero counts if not has_gamma")
4260 insist(sum(this%binned(:, 1448 ))== 0.0000000000000000 ," 440990 has zero counts if not has_gamma")
4262 insist(sum(this%binned(:, 1449 ))== 0.0000000000000000 ," 450990 has zero counts if not has_gamma")
4264 insist(sum(this%binned(:, 1450 ))== 0.0000000000000000 ," 450991 has zero counts if not has_gamma")
4266 insist(sum(this%binned(:, 1451 ))== 0.0000000000000000 ," 460990 has zero counts if not has_gamma")
4268 insist(sum(this%binned(:, 1452 ))== 0.0000000000000000 ," 361000 has zero counts if not has_gamma")
4270 insist(sum(this%binned(:, 1453 ))== 0.0000000000000000 ," 371000 has zero counts if not has_gamma")
4272 insist(sum(this%binned(:, 1454 ))== 0.0000000000000000 ," 381000 has zero counts if not has_gamma")
4274 insist(sum(this%binned(:, 1455 ))== 0.0000000000000000 ," 391000 has zero counts if not has_gamma")
4276 insist(sum(this%binned(:, 1456 ))== 0.0000000000000000 ," 401000 has zero counts if not has_gamma")
4278 insist(sum(this%binned(:, 1457 ))== 0.0000000000000000 ," 411000 has zero counts if not has_gamma")
4280 insist(sum(this%binned(:, 1458 ))== 0.0000000000000000 ," 411001 has zero counts if not has_gamma")
4282 insist(sum(this%binned(:, 1459 ))== 0.0000000000000000 ," 421000 has zero counts if not has_gamma")
4289 insist(sum(this%binned(:, 1461 ))== 0.0000000000000000 ," 441000 has zero counts if not has_gamma")
4291 insist(sum(this%binned(:, 1462 ))== 0.0000000000000000 ," 371010 has zero counts if not has_gamma")
4293 insist(sum(this%binned(:, 1463 ))== 0.0000000000000000 ," 381010 has zero counts if not has_gamma")
4295 insist(sum(this%binned(:, 1464 ))== 0.0000000000000000 ," 391010 has zero counts if not has_gamma")
4302 insist(sum(this%binned(:, 1466 ))== 0.0000000000000000 ," 411010 has zero counts if not has_gamma")
4306 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.63661447291623863 , new=",new,")"
4311 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.50693567226374714 , new=",new,")"
4314 insist(sum(this%binned(:, 1469 ))== 0.0000000000000000 ," 441010 has zero counts if not has_gamma")
4316 insist(sum(this%binned(:, 1470 ))== 0.0000000000000000 ," 451010 has zero counts if not has_gamma")
4318 insist(sum(this%binned(:, 1471 ))== 0.0000000000000000 ," 451011 has zero counts if not has_gamma")
4320 insist(sum(this%binned(:, 1472 ))== 0.0000000000000000 ," 461010 has zero counts if not has_gamma")
4322 insist(sum(this%binned(:, 1473 ))== 0.0000000000000000 ," 371020 has zero counts if not has_gamma")
4324 insist(sum(this%binned(:, 1474 ))== 0.0000000000000000 ," 381020 has zero counts if not has_gamma")
4326 insist(sum(this%binned(:, 1475 ))== 0.0000000000000000 ," 391020 has zero counts if not has_gamma")
4328 insist(sum(this%binned(:, 1476 ))== 0.0000000000000000 ," 401020 has zero counts if not has_gamma")
4330 insist(sum(this%binned(:, 1477 ))== 0.0000000000000000 ," 411020 has zero counts if not has_gamma")
4332 insist(sum(this%binned(:, 1478 ))== 0.0000000000000000 ," 411021 has zero counts if not has_gamma")
4334 insist(sum(this%binned(:, 1479 ))== 0.0000000000000000 ," 421020 has zero counts if not has_gamma")
4346 insist(sum(this%binned(:, 1482 ))== 0.0000000000000000 ," 441020 has zero counts if not has_gamma")
4348 insist(sum(this%binned(:, 1483 ))== 0.0000000000000000 ," 451020 has zero counts if not has_gamma")
4350 insist(sum(this%binned(:, 1484 ))== 0.0000000000000000 ," 451021 has zero counts if not has_gamma")
4352 insist(sum(this%binned(:, 1485 ))== 0.0000000000000000 ," 461020 has zero counts if not has_gamma")
4354 insist(sum(this%binned(:, 1486 ))== 0.0000000000000000 ," 381030 has zero counts if not has_gamma")
4356 insist(sum(this%binned(:, 1487 ))== 0.0000000000000000 ," 391030 has zero counts if not has_gamma")
4358 insist(sum(this%binned(:, 1488 ))== 0.0000000000000000 ," 401030 has zero counts if not has_gamma")
4360 insist(sum(this%binned(:, 1489 ))== 0.0000000000000000 ," 411030 has zero counts if not has_gamma")
4362 insist(sum(this%binned(:, 1490 ))== 0.0000000000000000 ," 421030 has zero counts if not has_gamma")
4364 insist(sum(this%binned(:, 1491 ))== 0.0000000000000000 ," 431030 has zero counts if not has_gamma")
4368 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 5.6067314364812035E-002 , new=",new,")"
4371 insist(sum(this%binned(:, 1493 ))== 0.0000000000000000 ," 451030 has zero counts if not has_gamma")
4373 insist(sum(this%binned(:, 1494 ))== 0.0000000000000000 ," 451031 has zero counts if not has_gamma")
4375 insist(sum(this%binned(:, 1495 ))== 0.0000000000000000 ," 461030 has zero counts if not has_gamma")
4377 insist(sum(this%binned(:, 1496 ))== 0.0000000000000000 ," 471030 has zero counts if not has_gamma")
4379 insist(sum(this%binned(:, 1497 ))== 0.0000000000000000 ," 381040 has zero counts if not has_gamma")
4381 insist(sum(this%binned(:, 1498 ))== 0.0000000000000000 ," 391040 has zero counts if not has_gamma")
4383 insist(sum(this%binned(:, 1499 ))== 0.0000000000000000 ," 401040 has zero counts if not has_gamma")
4385 insist(sum(this%binned(:, 1500 ))== 0.0000000000000000 ," 411040 has zero counts if not has_gamma")
4387 insist(sum(this%binned(:, 1501 ))== 0.0000000000000000 ," 411041 has zero counts if not has_gamma")
4389 insist(sum(this%binned(:, 1502 ))== 0.0000000000000000 ," 421040 has zero counts if not has_gamma")
4396 insist(sum(this%binned(:, 1504 ))== 0.0000000000000000 ," 441040 has zero counts if not has_gamma")
4405 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 5.0072926284761812E-007 , new=",new,")"
4408 insist(sum(this%binned(:, 1507 ))== 0.0000000000000000 ," 461040 has zero counts if not has_gamma")
4410 insist(sum(this%binned(:, 1508 ))== 0.0000000000000000 ," 381050 has zero counts if not has_gamma")
4412 insist(sum(this%binned(:, 1509 ))== 0.0000000000000000 ," 391050 has zero counts if not has_gamma")
4414 insist(sum(this%binned(:, 1510 ))== 0.0000000000000000 ," 401050 has zero counts if not has_gamma")
4416 insist(sum(this%binned(:, 1511 ))== 0.0000000000000000 ," 411050 has zero counts if not has_gamma")
4418 insist(sum(this%binned(:, 1512 ))== 0.0000000000000000 ," 421050 has zero counts if not has_gamma")
4420 insist(sum(this%binned(:, 1513 ))== 0.0000000000000000 ," 431050 has zero counts if not has_gamma")
4424 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.41997130747915101 , new=",new,")"
4429 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.13824386501530891 , new=",new,")"
4432 insist(sum(this%binned(:, 1516 ))== 0.0000000000000000 ," 451051 has zero counts if not has_gamma")
4434 insist(sum(this%binned(:, 1517 ))== 0.0000000000000000 ," 461050 has zero counts if not has_gamma")
4436 insist(sum(this%binned(:, 1518 ))== 0.0000000000000000 ," 471050 has zero counts if not has_gamma")
4438 insist(sum(this%binned(:, 1519 ))== 0.0000000000000000 ," 471051 has zero counts if not has_gamma")
4440 insist(sum(this%binned(:, 1520 ))== 0.0000000000000000 ," 481050 has zero counts if not has_gamma")
4442 insist(sum(this%binned(:, 1521 ))== 0.0000000000000000 ," 391060 has zero counts if not has_gamma")
4444 insist(sum(this%binned(:, 1522 ))== 0.0000000000000000 ," 401060 has zero counts if not has_gamma")
4446 insist(sum(this%binned(:, 1523 ))== 0.0000000000000000 ," 411060 has zero counts if not has_gamma")
4448 insist(sum(this%binned(:, 1524 ))== 0.0000000000000000 ," 421060 has zero counts if not has_gamma")
4450 insist(sum(this%binned(:, 1525 ))== 0.0000000000000000 ," 431060 has zero counts if not has_gamma")
4454 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 2.5935015206863810E-003 , new=",new,")"
4464 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.30463977297928585 , new=",new,")"
4467 insist(sum(this%binned(:, 1529 ))== 0.0000000000000000 ," 461060 has zero counts if not has_gamma")
4469 insist(sum(this%binned(:, 1530 ))== 0.0000000000000000 ," 471060 has zero counts if not has_gamma")
4471 insist(sum(this%binned(:, 1531 ))== 0.0000000000000000 ," 471061 has zero counts if not has_gamma")
4473 insist(sum(this%binned(:, 1532 ))== 0.0000000000000000 ," 481060 has zero counts if not has_gamma")
4475 insist(sum(this%binned(:, 1533 ))== 0.0000000000000000 ," 391070 has zero counts if not has_gamma")
4477 insist(sum(this%binned(:, 1534 ))== 0.0000000000000000 ," 401070 has zero counts if not has_gamma")
4479 insist(sum(this%binned(:, 1535 ))== 0.0000000000000000 ," 411070 has zero counts if not has_gamma")
4481 insist(sum(this%binned(:, 1536 ))== 0.0000000000000000 ," 421070 has zero counts if not has_gamma")
4483 insist(sum(this%binned(:, 1537 ))== 0.0000000000000000 ," 431070 has zero counts if not has_gamma")
4492 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.49020845414431330 , new=",new,")"
4497 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 2.1676466463853572E-003 , new=",new,")"
4500 insist(sum(this%binned(:, 1541 ))== 0.0000000000000000 ," 461071 has zero counts if not has_gamma")
4502 insist(sum(this%binned(:, 1542 ))== 0.0000000000000000 ," 471070 has zero counts if not has_gamma")
4504 insist(sum(this%binned(:, 1543 ))== 0.0000000000000000 ," 471071 has zero counts if not has_gamma")
4508 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.3637423944359017E-003 , new=",new,")"
4511 insist(sum(this%binned(:, 1545 ))== 0.0000000000000000 ," 491070 has zero counts if not has_gamma")
4513 insist(sum(this%binned(:, 1546 ))== 0.0000000000000000 ," 391080 has zero counts if not has_gamma")
4515 insist(sum(this%binned(:, 1547 ))== 0.0000000000000000 ," 401080 has zero counts if not has_gamma")
4517 insist(sum(this%binned(:, 1548 ))== 0.0000000000000000 ," 411080 has zero counts if not has_gamma")
4519 insist(sum(this%binned(:, 1549 ))== 0.0000000000000000 ," 421080 has zero counts if not has_gamma")
4521 insist(sum(this%binned(:, 1550 ))== 0.0000000000000000 ," 431080 has zero counts if not has_gamma")
4523 insist(sum(this%binned(:, 1551 ))== 0.0000000000000000 ," 441080 has zero counts if not has_gamma")
4535 insist(sum(this%binned(:, 1554 ))== 0.0000000000000000 ," 461080 has zero counts if not has_gamma")
4539 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.66592332489877126 , new=",new,")"
4542 insist(sum(this%binned(:, 1556 ))== 0.0000000000000000 ," 471081 has zero counts if not has_gamma")
4544 insist(sum(this%binned(:, 1557 ))== 0.0000000000000000 ," 481080 has zero counts if not has_gamma")
4546 insist(sum(this%binned(:, 1558 ))== 0.0000000000000000 ," 401090 has zero counts if not has_gamma")
4548 insist(sum(this%binned(:, 1559 ))== 0.0000000000000000 ," 411090 has zero counts if not has_gamma")
4550 insist(sum(this%binned(:, 1560 ))== 0.0000000000000000 ," 421090 has zero counts if not has_gamma")
4552 insist(sum(this%binned(:, 1561 ))== 0.0000000000000000 ," 431090 has zero counts if not has_gamma")
4554 insist(sum(this%binned(:, 1562 ))== 0.0000000000000000 ," 441090 has zero counts if not has_gamma")
4556 insist(sum(this%binned(:, 1563 ))== 0.0000000000000000 ," 451090 has zero counts if not has_gamma")
4560 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.37413798787321895 , new=",new,")"
4563 insist(sum(this%binned(:, 1565 ))== 0.0000000000000000 ," 461091 has zero counts if not has_gamma")
4565 insist(sum(this%binned(:, 1566 ))== 0.0000000000000000 ," 471090 has zero counts if not has_gamma")
4567 insist(sum(this%binned(:, 1567 ))== 0.0000000000000000 ," 471091 has zero counts if not has_gamma")
4569 insist(sum(this%binned(:, 1568 ))== 0.0000000000000000 ," 481090 has zero counts if not has_gamma")
4571 insist(sum(this%binned(:, 1569 ))== 0.0000000000000000 ," 491090 has zero counts if not has_gamma")
4573 insist(sum(this%binned(:, 1570 ))== 0.0000000000000000 ," 401100 has zero counts if not has_gamma")
4575 insist(sum(this%binned(:, 1571 ))== 0.0000000000000000 ," 411100 has zero counts if not has_gamma")
4577 insist(sum(this%binned(:, 1572 ))== 0.0000000000000000 ," 421100 has zero counts if not has_gamma")
4579 insist(sum(this%binned(:, 1573 ))== 0.0000000000000000 ," 431100 has zero counts if not has_gamma")
4581 insist(sum(this%binned(:, 1574 ))== 0.0000000000000000 ," 441100 has zero counts if not has_gamma")
4593 insist(sum(this%binned(:, 1577 ))== 0.0000000000000000 ," 461100 has zero counts if not has_gamma")
4602 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 5.2684700660961419E-002 , new=",new,")"
4605 insist(sum(this%binned(:, 1580 ))== 0.0000000000000000 ," 481100 has zero counts if not has_gamma")
4607 insist(sum(this%binned(:, 1581 ))== 0.0000000000000000 ," 411110 has zero counts if not has_gamma")
4609 insist(sum(this%binned(:, 1582 ))== 0.0000000000000000 ," 421110 has zero counts if not has_gamma")
4611 insist(sum(this%binned(:, 1583 ))== 0.0000000000000000 ," 431110 has zero counts if not has_gamma")
4613 insist(sum(this%binned(:, 1584 ))== 0.0000000000000000 ," 441110 has zero counts if not has_gamma")
4615 insist(sum(this%binned(:, 1585 ))== 0.0000000000000000 ," 451110 has zero counts if not has_gamma")
4617 insist(sum(this%binned(:, 1586 ))== 0.0000000000000000 ," 461110 has zero counts if not has_gamma")
4619 insist(sum(this%binned(:, 1587 ))== 0.0000000000000000 ," 461111 has zero counts if not has_gamma")
4623 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.36628860078762671 , new=",new,")"
4626 insist(sum(this%binned(:, 1589 ))== 0.0000000000000000 ," 471111 has zero counts if not has_gamma")
4628 insist(sum(this%binned(:, 1590 ))== 0.0000000000000000 ," 481110 has zero counts if not has_gamma")
4630 insist(sum(this%binned(:, 1591 ))== 0.0000000000000000 ," 481111 has zero counts if not has_gamma")
4632 insist(sum(this%binned(:, 1592 ))== 0.0000000000000000 ," 491110 has zero counts if not has_gamma")
4634 insist(sum(this%binned(:, 1593 ))== 0.0000000000000000 ," 491111 has zero counts if not has_gamma")
4636 insist(sum(this%binned(:, 1594 ))== 0.0000000000000000 ," 501110 has zero counts if not has_gamma")
4638 insist(sum(this%binned(:, 1595 ))== 0.0000000000000000 ," 411120 has zero counts if not has_gamma")
4640 insist(sum(this%binned(:, 1596 ))== 0.0000000000000000 ," 421120 has zero counts if not has_gamma")
4642 insist(sum(this%binned(:, 1597 ))== 0.0000000000000000 ," 431120 has zero counts if not has_gamma")
4644 insist(sum(this%binned(:, 1598 ))== 0.0000000000000000 ," 441120 has zero counts if not has_gamma")
4646 insist(sum(this%binned(:, 1599 ))== 0.0000000000000000 ," 451120 has zero counts if not has_gamma")
4648 insist(sum(this%binned(:, 1600 ))== 0.0000000000000000 ," 461120 has zero counts if not has_gamma")
4650 insist(sum(this%binned(:, 1601 ))== 0.0000000000000000 ," 471120 has zero counts if not has_gamma")
4652 insist(sum(this%binned(:, 1602 ))== 0.0000000000000000 ," 481120 has zero counts if not has_gamma")
4654 insist(sum(this%binned(:, 1603 ))== 0.0000000000000000 ," 491120 has zero counts if not has_gamma")
4656 insist(sum(this%binned(:, 1604 ))== 0.0000000000000000 ," 491121 has zero counts if not has_gamma")
4658 insist(sum(this%binned(:, 1605 ))== 0.0000000000000000 ," 501120 has zero counts if not has_gamma")
4660 insist(sum(this%binned(:, 1606 ))== 0.0000000000000000 ," 411130 has zero counts if not has_gamma")
4662 insist(sum(this%binned(:, 1607 ))== 0.0000000000000000 ," 421130 has zero counts if not has_gamma")
4664 insist(sum(this%binned(:, 1608 ))== 0.0000000000000000 ," 431130 has zero counts if not has_gamma")
4666 insist(sum(this%binned(:, 1609 ))== 0.0000000000000000 ," 441130 has zero counts if not has_gamma")
4668 insist(sum(this%binned(:, 1610 ))== 0.0000000000000000 ," 451130 has zero counts if not has_gamma")
4670 insist(sum(this%binned(:, 1611 ))== 0.0000000000000000 ," 461130 has zero counts if not has_gamma")
4672 insist(sum(this%binned(:, 1612 ))== 0.0000000000000000 ," 471130 has zero counts if not has_gamma")
4674 insist(sum(this%binned(:, 1613 ))== 0.0000000000000000 ," 471131 has zero counts if not has_gamma")
4678 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 7.5475555081041520E-002 , new=",new,")"
4683 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.17320447776984155 , new=",new,")"
4686 insist(sum(this%binned(:, 1616 ))== 0.0000000000000000 ," 491130 has zero counts if not has_gamma")
4688 insist(sum(this%binned(:, 1617 ))== 0.0000000000000000 ," 491131 has zero counts if not has_gamma")
4690 insist(sum(this%binned(:, 1618 ))== 0.0000000000000000 ," 501130 has zero counts if not has_gamma")
4692 insist(sum(this%binned(:, 1619 ))== 0.0000000000000000 ," 501131 has zero counts if not has_gamma")
4694 insist(sum(this%binned(:, 1620 ))== 0.0000000000000000 ," 511130 has zero counts if not has_gamma")
4696 insist(sum(this%binned(:, 1621 ))== 0.0000000000000000 ," 421140 has zero counts if not has_gamma")
4698 insist(sum(this%binned(:, 1622 ))== 0.0000000000000000 ," 431140 has zero counts if not has_gamma")
4700 insist(sum(this%binned(:, 1623 ))== 0.0000000000000000 ," 441140 has zero counts if not has_gamma")
4702 insist(sum(this%binned(:, 1624 ))== 0.0000000000000000 ," 451140 has zero counts if not has_gamma")
4704 insist(sum(this%binned(:, 1625 ))== 0.0000000000000000 ," 461140 has zero counts if not has_gamma")
4711 insist(sum(this%binned(:, 1627 ))== 0.0000000000000000 ," 481140 has zero counts if not has_gamma")
4715 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.88342222361784639 , new=",new,")"
4718 insist(sum(this%binned(:, 1629 ))== 0.0000000000000000 ," 491141 has zero counts if not has_gamma")
4720 insist(sum(this%binned(:, 1630 ))== 0.0000000000000000 ," 501140 has zero counts if not has_gamma")
4722 insist(sum(this%binned(:, 1631 ))== 0.0000000000000000 ," 421150 has zero counts if not has_gamma")
4724 insist(sum(this%binned(:, 1632 ))== 0.0000000000000000 ," 431150 has zero counts if not has_gamma")
4726 insist(sum(this%binned(:, 1633 ))== 0.0000000000000000 ," 441150 has zero counts if not has_gamma")
4728 insist(sum(this%binned(:, 1634 ))== 0.0000000000000000 ," 451150 has zero counts if not has_gamma")
4730 insist(sum(this%binned(:, 1635 ))== 0.0000000000000000 ," 461150 has zero counts if not has_gamma")
4737 insist(sum(this%binned(:, 1637 ))== 0.0000000000000000 ," 471151 has zero counts if not has_gamma")
4741 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.30763816677283629 , new=",new,")"
4746 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.67684358372750064 , new=",new,")"
4751 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.13679583884205912 , new=",new,")"
4756 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 4.0115862475793927E-004 , new=",new,")"
4759 insist(sum(this%binned(:, 1642 ))== 0.0000000000000000 ," 501150 has zero counts if not has_gamma")
4761 insist(sum(this%binned(:, 1643 ))== 0.0000000000000000 ," 511150 has zero counts if not has_gamma")
4763 insist(sum(this%binned(:, 1644 ))== 0.0000000000000000 ," 521150 has zero counts if not has_gamma")
4765 insist(sum(this%binned(:, 1645 ))== 0.0000000000000000 ," 431160 has zero counts if not has_gamma")
4767 insist(sum(this%binned(:, 1646 ))== 0.0000000000000000 ," 441160 has zero counts if not has_gamma")
4769 insist(sum(this%binned(:, 1647 ))== 0.0000000000000000 ," 451160 has zero counts if not has_gamma")
4771 insist(sum(this%binned(:, 1648 ))== 0.0000000000000000 ," 461160 has zero counts if not has_gamma")
4783 insist(sum(this%binned(:, 1651 ))== 0.0000000000000000 ," 481160 has zero counts if not has_gamma")
4792 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.30791418396150139 , new=",new,")"
4795 insist(sum(this%binned(:, 1654 ))== 0.0000000000000000 ," 501160 has zero counts if not has_gamma")
4797 insist(sum(this%binned(:, 1655 ))== 0.0000000000000000 ," 431170 has zero counts if not has_gamma")
4799 insist(sum(this%binned(:, 1656 ))== 0.0000000000000000 ," 441170 has zero counts if not has_gamma")
4801 insist(sum(this%binned(:, 1657 ))== 0.0000000000000000 ," 451170 has zero counts if not has_gamma")
4803 insist(sum(this%binned(:, 1658 ))== 0.0000000000000000 ," 461170 has zero counts if not has_gamma")
4805 insist(sum(this%binned(:, 1659 ))== 0.0000000000000000 ," 471170 has zero counts if not has_gamma")
4807 insist(sum(this%binned(:, 1660 ))== 0.0000000000000000 ," 471171 has zero counts if not has_gamma")
4809 insist(sum(this%binned(:, 1661 ))== 0.0000000000000000 ," 481170 has zero counts if not has_gamma")
4811 insist(sum(this%binned(:, 1662 ))== 0.0000000000000000 ," 481171 has zero counts if not has_gamma")
4813 insist(sum(this%binned(:, 1663 ))== 0.0000000000000000 ," 491170 has zero counts if not has_gamma")
4815 insist(sum(this%binned(:, 1664 ))== 0.0000000000000000 ," 491171 has zero counts if not has_gamma")
4817 insist(sum(this%binned(:, 1665 ))== 0.0000000000000000 ," 501170 has zero counts if not has_gamma")
4819 insist(sum(this%binned(:, 1666 ))== 0.0000000000000000 ," 501171 has zero counts if not has_gamma")
4821 insist(sum(this%binned(:, 1667 ))== 0.0000000000000000 ," 511170 has zero counts if not has_gamma")
4823 insist(sum(this%binned(:, 1668 ))== 0.0000000000000000 ," 521170 has zero counts if not has_gamma")
4825 insist(sum(this%binned(:, 1669 ))== 0.0000000000000000 ," 431180 has zero counts if not has_gamma")
4827 insist(sum(this%binned(:, 1670 ))== 0.0000000000000000 ," 441180 has zero counts if not has_gamma")
4829 insist(sum(this%binned(:, 1671 ))== 0.0000000000000000 ," 451180 has zero counts if not has_gamma")
4831 insist(sum(this%binned(:, 1672 ))== 0.0000000000000000 ," 461180 has zero counts if not has_gamma")
4840 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.86840080777885975 , new=",new,")"
4843 insist(sum(this%binned(:, 1675 ))== 0.0000000000000000 ," 481180 has zero counts if not has_gamma")
4852 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.63419669953075919 , new=",new,")"
4855 insist(sum(this%binned(:, 1678 ))== 0.0000000000000000 ," 501180 has zero counts if not has_gamma")
4857 insist(sum(this%binned(:, 1679 ))== 0.0000000000000000 ," 511180 has zero counts if not has_gamma")
4859 insist(sum(this%binned(:, 1680 ))== 0.0000000000000000 ," 511181 has zero counts if not has_gamma")
4861 insist(sum(this%binned(:, 1681 ))== 0.0000000000000000 ," 521180 has zero counts if not has_gamma")
4863 insist(sum(this%binned(:, 1682 ))== 0.0000000000000000 ," 441190 has zero counts if not has_gamma")
4865 insist(sum(this%binned(:, 1683 ))== 0.0000000000000000 ," 451190 has zero counts if not has_gamma")
4867 insist(sum(this%binned(:, 1684 ))== 0.0000000000000000 ," 461190 has zero counts if not has_gamma")
4869 insist(sum(this%binned(:, 1685 ))== 0.0000000000000000 ," 471190 has zero counts if not has_gamma")
4871 insist(sum(this%binned(:, 1686 ))== 0.0000000000000000 ," 481190 has zero counts if not has_gamma")
4873 insist(sum(this%binned(:, 1687 ))== 0.0000000000000000 ," 481191 has zero counts if not has_gamma")
4875 insist(sum(this%binned(:, 1688 ))== 0.0000000000000000 ," 491190 has zero counts if not has_gamma")
4877 insist(sum(this%binned(:, 1689 ))== 0.0000000000000000 ," 491191 has zero counts if not has_gamma")
4879 insist(sum(this%binned(:, 1690 ))== 0.0000000000000000 ," 501190 has zero counts if not has_gamma")
4881 insist(sum(this%binned(:, 1691 ))== 0.0000000000000000 ," 501191 has zero counts if not has_gamma")
4883 insist(sum(this%binned(:, 1692 ))== 0.0000000000000000 ," 511190 has zero counts if not has_gamma")
4885 insist(sum(this%binned(:, 1693 ))== 0.0000000000000000 ," 521190 has zero counts if not has_gamma")
4887 insist(sum(this%binned(:, 1694 ))== 0.0000000000000000 ," 441200 has zero counts if not has_gamma")
4889 insist(sum(this%binned(:, 1695 ))== 0.0000000000000000 ," 451200 has zero counts if not has_gamma")
4891 insist(sum(this%binned(:, 1696 ))== 0.0000000000000000 ," 461200 has zero counts if not has_gamma")
4898 insist(sum(this%binned(:, 1698 ))== 0.0000000000000000 ," 471201 has zero counts if not has_gamma")
4902 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.61944403806587356 , new=",new,")"
4915 insist(sum(this%binned(:, 1702 ))== 0.0000000000000000 ," 501200 has zero counts if not has_gamma")
4917 insist(sum(this%binned(:, 1703 ))== 0.0000000000000000 ," 511200 has zero counts if not has_gamma")
4919 insist(sum(this%binned(:, 1704 ))== 0.0000000000000000 ," 511201 has zero counts if not has_gamma")
4921 insist(sum(this%binned(:, 1705 ))== 0.0000000000000000 ," 521200 has zero counts if not has_gamma")
4923 insist(sum(this%binned(:, 1706 ))== 0.0000000000000000 ," 451210 has zero counts if not has_gamma")
4925 insist(sum(this%binned(:, 1707 ))== 0.0000000000000000 ," 461210 has zero counts if not has_gamma")
4927 insist(sum(this%binned(:, 1708 ))== 0.0000000000000000 ," 471210 has zero counts if not has_gamma")
4929 insist(sum(this%binned(:, 1709 ))== 0.0000000000000000 ," 481210 has zero counts if not has_gamma")
4931 insist(sum(this%binned(:, 1710 ))== 0.0000000000000000 ," 481211 has zero counts if not has_gamma")
4933 insist(sum(this%binned(:, 1711 ))== 0.0000000000000000 ," 491210 has zero counts if not has_gamma")
4935 insist(sum(this%binned(:, 1712 ))== 0.0000000000000000 ," 491211 has zero counts if not has_gamma")
4937 insist(sum(this%binned(:, 1713 ))== 0.0000000000000000 ," 501210 has zero counts if not has_gamma")
4939 insist(sum(this%binned(:, 1714 ))== 0.0000000000000000 ," 501211 has zero counts if not has_gamma")
4941 insist(sum(this%binned(:, 1715 ))== 0.0000000000000000 ," 511210 has zero counts if not has_gamma")
4943 insist(sum(this%binned(:, 1716 ))== 0.0000000000000000 ," 521210 has zero counts if not has_gamma")
4945 insist(sum(this%binned(:, 1717 ))== 0.0000000000000000 ," 521211 has zero counts if not has_gamma")
4947 insist(sum(this%binned(:, 1718 ))== 0.0000000000000000 ," 531210 has zero counts if not has_gamma")
4949 insist(sum(this%binned(:, 1719 ))== 0.0000000000000000 ," 451220 has zero counts if not has_gamma")
4951 insist(sum(this%binned(:, 1720 ))== 0.0000000000000000 ," 461220 has zero counts if not has_gamma")
4953 insist(sum(this%binned(:, 1721 ))== 0.0000000000000000 ," 471220 has zero counts if not has_gamma")
4955 insist(sum(this%binned(:, 1722 ))== 0.0000000000000000 ," 471221 has zero counts if not has_gamma")
4957 insist(sum(this%binned(:, 1723 ))== 0.0000000000000000 ," 481220 has zero counts if not has_gamma")
4969 insist(sum(this%binned(:, 1726 ))== 0.0000000000000000 ," 501220 has zero counts if not has_gamma")
4973 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.61206787508708693 , new=",new,")"
4976 insist(sum(this%binned(:, 1728 ))== 0.0000000000000000 ," 511221 has zero counts if not has_gamma")
4978 insist(sum(this%binned(:, 1729 ))== 0.0000000000000000 ," 521220 has zero counts if not has_gamma")
4980 insist(sum(this%binned(:, 1730 ))== 0.0000000000000000 ," 461230 has zero counts if not has_gamma")
4982 insist(sum(this%binned(:, 1731 ))== 0.0000000000000000 ," 471230 has zero counts if not has_gamma")
4984 insist(sum(this%binned(:, 1732 ))== 0.0000000000000000 ," 481230 has zero counts if not has_gamma")
4986 insist(sum(this%binned(:, 1733 ))== 0.0000000000000000 ," 481231 has zero counts if not has_gamma")
4993 insist(sum(this%binned(:, 1735 ))== 0.0000000000000000 ," 491231 has zero counts if not has_gamma")
4997 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.57978316783614969 , new=",new,")"
5002 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.48882191691704341 , new=",new,")"
5005 insist(sum(this%binned(:, 1738 ))== 0.0000000000000000 ," 511230 has zero counts if not has_gamma")
5007 insist(sum(this%binned(:, 1739 ))== 0.0000000000000000 ," 521230 has zero counts if not has_gamma")
5009 insist(sum(this%binned(:, 1740 ))== 0.0000000000000000 ," 521231 has zero counts if not has_gamma")
5011 insist(sum(this%binned(:, 1741 ))== 0.0000000000000000 ," 531230 has zero counts if not has_gamma")
5013 insist(sum(this%binned(:, 1742 ))== 0.0000000000000000 ," 461240 has zero counts if not has_gamma")
5015 insist(sum(this%binned(:, 1743 ))== 0.0000000000000000 ," 471240 has zero counts if not has_gamma")
5017 insist(sum(this%binned(:, 1744 ))== 0.0000000000000000 ," 481240 has zero counts if not has_gamma")
5019 insist(sum(this%binned(:, 1745 ))== 0.0000000000000000 ," 491240 has zero counts if not has_gamma")
5021 insist(sum(this%binned(:, 1746 ))== 0.0000000000000000 ," 491241 has zero counts if not has_gamma")
5023 insist(sum(this%binned(:, 1747 ))== 0.0000000000000000 ," 501240 has zero counts if not has_gamma")
5027 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.41775164849608570 , new=",new,")"
5032 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.7386910510827434E-002 , new=",new,")"
5035 insist(sum(this%binned(:, 1750 ))== 0.0000000000000000 ," 521240 has zero counts if not has_gamma")
5037 insist(sum(this%binned(:, 1751 ))== 0.0000000000000000 ," 531240 has zero counts if not has_gamma")
5039 insist(sum(this%binned(:, 1752 ))== 0.0000000000000000 ," 541240 has zero counts if not has_gamma")
5041 insist(sum(this%binned(:, 1753 ))== 0.0000000000000000 ," 471250 has zero counts if not has_gamma")
5043 insist(sum(this%binned(:, 1754 ))== 0.0000000000000000 ," 481250 has zero counts if not has_gamma")
5045 insist(sum(this%binned(:, 1755 ))== 0.0000000000000000 ," 491250 has zero counts if not has_gamma")
5047 insist(sum(this%binned(:, 1756 ))== 0.0000000000000000 ," 491251 has zero counts if not has_gamma")
5051 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.95998007878166880 , new=",new,")"
5056 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.90123469909512754 , new=",new,")"
5061 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 7.1726213357862861E-002 , new=",new,")"
5064 insist(sum(this%binned(:, 1760 ))== 0.0000000000000000 ," 521250 has zero counts if not has_gamma")
5066 insist(sum(this%binned(:, 1761 ))== 0.0000000000000000 ," 521251 has zero counts if not has_gamma")
5068 insist(sum(this%binned(:, 1762 ))== 0.0000000000000000 ," 531250 has zero counts if not has_gamma")
5072 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 5.0064228739966918E-003 , new=",new,")"
5075 insist(sum(this%binned(:, 1764 ))== 0.0000000000000000 ," 541251 has zero counts if not has_gamma")
5077 insist(sum(this%binned(:, 1765 ))== 0.0000000000000000 ," 471260 has zero counts if not has_gamma")
5079 insist(sum(this%binned(:, 1766 ))== 0.0000000000000000 ," 481260 has zero counts if not has_gamma")
5081 insist(sum(this%binned(:, 1767 ))== 0.0000000000000000 ," 491260 has zero counts if not has_gamma")
5083 insist(sum(this%binned(:, 1768 ))== 0.0000000000000000 ," 491261 has zero counts if not has_gamma")
5087 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.4624188393718172E-002 , new=",new,")"
5092 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.28249918769960813 , new=",new,")"
5097 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.57419892530658634 , new=",new,")"
5100 insist(sum(this%binned(:, 1772 ))== 0.0000000000000000 ," 521260 has zero counts if not has_gamma")
5104 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 6.9957713445655045E-002 , new=",new,")"
5107 insist(sum(this%binned(:, 1774 ))== 0.0000000000000000 ," 541260 has zero counts if not has_gamma")
5109 insist(sum(this%binned(:, 1775 ))== 0.0000000000000000 ," 471270 has zero counts if not has_gamma")
5111 insist(sum(this%binned(:, 1776 ))== 0.0000000000000000 ," 481270 has zero counts if not has_gamma")
5113 insist(sum(this%binned(:, 1777 ))== 0.0000000000000000 ," 491270 has zero counts if not has_gamma")
5115 insist(sum(this%binned(:, 1778 ))== 0.0000000000000000 ," 491271 has zero counts if not has_gamma")
5119 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.61070897067947871 , new=",new,")"
5122 insist(sum(this%binned(:, 1780 ))== 0.0000000000000000 ," 501271 has zero counts if not has_gamma")
5126 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.31553111973075260 , new=",new,")"
5131 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.21533089173324707 , new=",new,")"
5136 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.4553212944834476E-004 , new=",new,")"
5139 insist(sum(this%binned(:, 1784 ))== 0.0000000000000000 ," 531270 has zero counts if not has_gamma")
5141 insist(sum(this%binned(:, 1785 ))== 0.0000000000000000 ," 541270 has zero counts if not has_gamma")
5143 insist(sum(this%binned(:, 1786 ))== 0.0000000000000000 ," 541271 has zero counts if not has_gamma")
5145 insist(sum(this%binned(:, 1787 ))== 0.0000000000000000 ," 551270 has zero counts if not has_gamma")
5147 insist(sum(this%binned(:, 1788 ))== 0.0000000000000000 ," 471280 has zero counts if not has_gamma")
5149 insist(sum(this%binned(:, 1789 ))== 0.0000000000000000 ," 481280 has zero counts if not has_gamma")
5151 insist(sum(this%binned(:, 1790 ))== 0.0000000000000000 ," 491280 has zero counts if not has_gamma")
5153 insist(sum(this%binned(:, 1791 ))== 0.0000000000000000 ," 491281 has zero counts if not has_gamma")
5155 insist(sum(this%binned(:, 1792 ))== 0.0000000000000000 ," 501280 has zero counts if not has_gamma")
5157 insist(sum(this%binned(:, 1793 ))== 0.0000000000000000 ," 501281 has zero counts if not has_gamma")
5161 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.41964379699613336 , new=",new,")"
5169 insist(sum(this%binned(:, 1796 ))== 0.0000000000000000 ," 521280 has zero counts if not has_gamma")
5173 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.81150481878300196 , new=",new,")"
5176 insist(sum(this%binned(:, 1798 ))== 0.0000000000000000 ," 541280 has zero counts if not has_gamma")
5178 insist(sum(this%binned(:, 1799 ))== 0.0000000000000000 ," 471290 has zero counts if not has_gamma")
5180 insist(sum(this%binned(:, 1800 ))== 0.0000000000000000 ," 481290 has zero counts if not has_gamma")
5182 insist(sum(this%binned(:, 1801 ))== 0.0000000000000000 ," 491290 has zero counts if not has_gamma")
5184 insist(sum(this%binned(:, 1802 ))== 0.0000000000000000 ," 491291 has zero counts if not has_gamma")
5186 insist(sum(this%binned(:, 1803 ))== 0.0000000000000000 ," 501290 has zero counts if not has_gamma")
5188 insist(sum(this%binned(:, 1804 ))== 0.0000000000000000 ," 501291 has zero counts if not has_gamma")
5192 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.44350162218044037 , new=",new,")"
5195 insist(sum(this%binned(:, 1806 ))== 0.0000000000000000 ," 511291 has zero counts if not has_gamma")
5199 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.57148233870467513 , new=",new,")"
5204 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 7.7115786881851278E-002 , new=",new,")"
5209 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 2.5070131073676859E-002 , new=",new,")"
5212 insist(sum(this%binned(:, 1810 ))== 0.0000000000000000 ," 541290 has zero counts if not has_gamma")
5214 insist(sum(this%binned(:, 1811 ))== 0.0000000000000000 ," 541291 has zero counts if not has_gamma")
5216 insist(sum(this%binned(:, 1812 ))== 0.0000000000000000 ," 551290 has zero counts if not has_gamma")
5218 insist(sum(this%binned(:, 1813 ))== 0.0000000000000000 ," 561290 has zero counts if not has_gamma")
5220 insist(sum(this%binned(:, 1814 ))== 0.0000000000000000 ," 471300 has zero counts if not has_gamma")
5222 insist(sum(this%binned(:, 1815 ))== 0.0000000000000000 ," 481300 has zero counts if not has_gamma")
5229 insist(sum(this%binned(:, 1817 ))== 0.0000000000000000 ," 491301 has zero counts if not has_gamma")
5233 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.36380849496143708 , new=",new,")"
5236 insist(sum(this%binned(:, 1819 ))== 0.0000000000000000 ," 501301 has zero counts if not has_gamma")
5240 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.76425034656235624 , new=",new,")"
5248 insist(sum(this%binned(:, 1822 ))== 0.0000000000000000 ," 521300 has zero counts if not has_gamma")
5252 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.29084157692338020 , new=",new,")"
5257 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 3.1578917906401223E-002 , new=",new,")"
5260 insist(sum(this%binned(:, 1825 ))== 0.0000000000000000 ," 541300 has zero counts if not has_gamma")
5262 insist(sum(this%binned(:, 1826 ))== 0.0000000000000000 ," 481310 has zero counts if not has_gamma")
5264 insist(sum(this%binned(:, 1827 ))== 0.0000000000000000 ," 491310 has zero counts if not has_gamma")
5266 insist(sum(this%binned(:, 1828 ))== 0.0000000000000000 ," 491311 has zero counts if not has_gamma")
5268 insist(sum(this%binned(:, 1829 ))== 0.0000000000000000 ," 501310 has zero counts if not has_gamma")
5270 insist(sum(this%binned(:, 1830 ))== 0.0000000000000000 ," 501311 has zero counts if not has_gamma")
5274 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.61258087053982990 , new=",new,")"
5279 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.79543965914349735 , new=",new,")"
5284 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.11407246174684361 , new=",new,")"
5289 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.16979873602472145 , new=",new,")"
5292 insist(sum(this%binned(:, 1835 ))== 0.0000000000000000 ," 541310 has zero counts if not has_gamma")
5294 insist(sum(this%binned(:, 1836 ))== 0.0000000000000000 ," 541311 has zero counts if not has_gamma")
5296 insist(sum(this%binned(:, 1837 ))== 0.0000000000000000 ," 551310 has zero counts if not has_gamma")
5300 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 4.0519245755252231E-006 , new=",new,")"
5303 insist(sum(this%binned(:, 1839 ))== 0.0000000000000000 ," 481320 has zero counts if not has_gamma")
5305 insist(sum(this%binned(:, 1840 ))== 0.0000000000000000 ," 491320 has zero counts if not has_gamma")
5309 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.84123627697318115 , new=",new,")"
5324 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 4.2364425689839891E-002 , new=",new,")"
5329 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.53271029128831782 , new=",new,")"
5332 insist(sum(this%binned(:, 1846 ))== 0.0000000000000000 ," 531321 has zero counts if not has_gamma")
5334 insist(sum(this%binned(:, 1847 ))== 0.0000000000000000 ," 541320 has zero counts if not has_gamma")
5338 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 2.1000282732913530E-003 , new=",new,")"
5341 insist(sum(this%binned(:, 1849 ))== 0.0000000000000000 ," 561320 has zero counts if not has_gamma")
5343 insist(sum(this%binned(:, 1850 ))== 0.0000000000000000 ," 491330 has zero counts if not has_gamma")
5345 insist(sum(this%binned(:, 1851 ))== 0.0000000000000000 ," 501330 has zero counts if not has_gamma")
5347 insist(sum(this%binned(:, 1852 ))== 0.0000000000000000 ," 511330 has zero counts if not has_gamma")
5351 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.94294755443960732 , new=",new,")"
5356 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.64882945766150601 , new=",new,")"
5361 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.43136643703652339 , new=",new,")"
5364 insist(sum(this%binned(:, 1856 ))== 0.0000000000000000 ," 531331 has zero counts if not has_gamma")
5368 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.2721398470994259E-002 , new=",new,")"
5371 insist(sum(this%binned(:, 1858 ))== 0.0000000000000000 ," 541331 has zero counts if not has_gamma")
5373 insist(sum(this%binned(:, 1859 ))== 0.0000000000000000 ," 551330 has zero counts if not has_gamma")
5375 insist(sum(this%binned(:, 1860 ))== 0.0000000000000000 ," 561330 has zero counts if not has_gamma")
5377 insist(sum(this%binned(:, 1861 ))== 0.0000000000000000 ," 571330 has zero counts if not has_gamma")
5379 insist(sum(this%binned(:, 1862 ))== 0.0000000000000000 ," 491340 has zero counts if not has_gamma")
5381 insist(sum(this%binned(:, 1863 ))== 0.0000000000000000 ," 501340 has zero counts if not has_gamma")
5388 insist(sum(this%binned(:, 1865 ))== 0.0000000000000000 ," 511341 has zero counts if not has_gamma")
5392 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.25847059190726424 , new=",new,")"
5397 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.69768931544394530 , new=",new,")"
5402 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 4.7926361991801342E-004 , new=",new,")"
5405 insist(sum(this%binned(:, 1869 ))== 0.0000000000000000 ," 541340 has zero counts if not has_gamma")
5407 insist(sum(this%binned(:, 1870 ))== 0.0000000000000000 ," 541341 has zero counts if not has_gamma")
5411 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.14673576315679460 , new=",new,")"
5414 insist(sum(this%binned(:, 1872 ))== 0.0000000000000000 ," 551341 has zero counts if not has_gamma")
5416 insist(sum(this%binned(:, 1873 ))== 0.0000000000000000 ," 561340 has zero counts if not has_gamma")
5418 insist(sum(this%binned(:, 1874 ))== 0.0000000000000000 ," 491350 has zero counts if not has_gamma")
5420 insist(sum(this%binned(:, 1875 ))== 0.0000000000000000 ," 501350 has zero counts if not has_gamma")
5422 insist(sum(this%binned(:, 1876 ))== 0.0000000000000000 ," 511350 has zero counts if not has_gamma")
5431 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.37838574961282295 , new=",new,")"
5436 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.30621217934464795 , new=",new,")"
5441 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.2799521752700906E-005 , new=",new,")"
5446 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 3.9443880041662993E-002 , new=",new,")"
5449 insist(sum(this%binned(:, 1882 ))== 0.0000000000000000 ," 551351 has zero counts if not has_gamma")
5451 insist(sum(this%binned(:, 1883 ))== 0.0000000000000000 ," 561350 has zero counts if not has_gamma")
5453 insist(sum(this%binned(:, 1884 ))== 0.0000000000000000 ," 561351 has zero counts if not has_gamma")
5455 insist(sum(this%binned(:, 1885 ))== 0.0000000000000000 ," 571350 has zero counts if not has_gamma")
5457 insist(sum(this%binned(:, 1886 ))== 0.0000000000000000 ," 581350 has zero counts if not has_gamma")
5459 insist(sum(this%binned(:, 1887 ))== 0.0000000000000000 ," 501360 has zero counts if not has_gamma")
5461 insist(sum(this%binned(:, 1888 ))== 0.0000000000000000 ," 511360 has zero counts if not has_gamma")
5463 insist(sum(this%binned(:, 1889 ))== 0.0000000000000000 ," 521360 has zero counts if not has_gamma")
5475 insist(sum(this%binned(:, 1892 ))== 0.0000000000000000 ," 541360 has zero counts if not has_gamma")
5479 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.5233914775956710E-002 , new=",new,")"
5482 insist(sum(this%binned(:, 1894 ))== 0.0000000000000000 ," 551361 has zero counts if not has_gamma")
5484 insist(sum(this%binned(:, 1895 ))== 0.0000000000000000 ," 561360 has zero counts if not has_gamma")
5486 insist(sum(this%binned(:, 1896 ))== 0.0000000000000000 ," 561361 has zero counts if not has_gamma")
5488 insist(sum(this%binned(:, 1897 ))== 0.0000000000000000 ," 501370 has zero counts if not has_gamma")
5490 insist(sum(this%binned(:, 1898 ))== 0.0000000000000000 ," 511370 has zero counts if not has_gamma")
5492 insist(sum(this%binned(:, 1899 ))== 0.0000000000000000 ," 521370 has zero counts if not has_gamma")
5494 insist(sum(this%binned(:, 1900 ))== 0.0000000000000000 ," 531370 has zero counts if not has_gamma")
5503 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.15841674472337336 , new=",new,")"
5506 insist(sum(this%binned(:, 1903 ))== 0.0000000000000000 ," 561370 has zero counts if not has_gamma")
5508 insist(sum(this%binned(:, 1904 ))== 0.0000000000000000 ," 561371 has zero counts if not has_gamma")
5510 insist(sum(this%binned(:, 1905 ))== 0.0000000000000000 ," 571370 has zero counts if not has_gamma")
5512 insist(sum(this%binned(:, 1906 ))== 0.0000000000000000 ," 581370 has zero counts if not has_gamma")
5514 insist(sum(this%binned(:, 1907 ))== 0.0000000000000000 ," 511380 has zero counts if not has_gamma")
5516 insist(sum(this%binned(:, 1908 ))== 0.0000000000000000 ," 521380 has zero counts if not has_gamma")
5525 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.75874852860065167 , new=",new,")"
5535 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.10407849520040637 , new=",new,")"
5538 insist(sum(this%binned(:, 1913 ))== 0.0000000000000000 ," 561380 has zero counts if not has_gamma")
5542 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.0724458313241581E-003 , new=",new,")"
5545 insist(sum(this%binned(:, 1915 ))== 0.0000000000000000 ," 581380 has zero counts if not has_gamma")
5547 insist(sum(this%binned(:, 1916 ))== 0.0000000000000000 ," 511390 has zero counts if not has_gamma")
5549 insist(sum(this%binned(:, 1917 ))== 0.0000000000000000 ," 521390 has zero counts if not has_gamma")
5551 insist(sum(this%binned(:, 1918 ))== 0.0000000000000000 ," 531390 has zero counts if not has_gamma")
5568 insist(sum(this%binned(:, 1922 ))== 0.0000000000000000 ," 571390 has zero counts if not has_gamma")
5570 insist(sum(this%binned(:, 1923 ))== 0.0000000000000000 ," 581390 has zero counts if not has_gamma")
5572 insist(sum(this%binned(:, 1924 ))== 0.0000000000000000 ," 581391 has zero counts if not has_gamma")
5574 insist(sum(this%binned(:, 1925 ))== 0.0000000000000000 ," 591390 has zero counts if not has_gamma")
5576 insist(sum(this%binned(:, 1926 ))== 0.0000000000000000 ," 521400 has zero counts if not has_gamma")
5578 insist(sum(this%binned(:, 1927 ))== 0.0000000000000000 ," 531400 has zero counts if not has_gamma")
5580 insist(sum(this%binned(:, 1928 ))== 0.0000000000000000 ," 541400 has zero counts if not has_gamma")
5582 insist(sum(this%binned(:, 1929 ))== 0.0000000000000000 ," 551400 has zero counts if not has_gamma")
5586 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.29417625363541500 , new=",new,")"
5591 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.60218382412752902 , new=",new,")"
5594 insist(sum(this%binned(:, 1932 ))== 0.0000000000000000 ," 581400 has zero counts if not has_gamma")
5596 insist(sum(this%binned(:, 1933 ))== 0.0000000000000000 ," 591400 has zero counts if not has_gamma")
5598 insist(sum(this%binned(:, 1934 ))== 0.0000000000000000 ," 601400 has zero counts if not has_gamma")
5600 insist(sum(this%binned(:, 1935 ))== 0.0000000000000000 ," 521410 has zero counts if not has_gamma")
5602 insist(sum(this%binned(:, 1936 ))== 0.0000000000000000 ," 531410 has zero counts if not has_gamma")
5604 insist(sum(this%binned(:, 1937 ))== 0.0000000000000000 ," 541410 has zero counts if not has_gamma")
5606 insist(sum(this%binned(:, 1938 ))== 0.0000000000000000 ," 551410 has zero counts if not has_gamma")
5620 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.13008160678741820 , new=",new,")"
5623 insist(sum(this%binned(:, 1942 ))== 0.0000000000000000 ," 591410 has zero counts if not has_gamma")
5625 insist(sum(this%binned(:, 1943 ))== 0.0000000000000000 ," 601410 has zero counts if not has_gamma")
5627 insist(sum(this%binned(:, 1944 ))== 0.0000000000000000 ," 601411 has zero counts if not has_gamma")
5629 insist(sum(this%binned(:, 1945 ))== 0.0000000000000000 ," 611410 has zero counts if not has_gamma")
5631 insist(sum(this%binned(:, 1946 ))== 0.0000000000000000 ," 521420 has zero counts if not has_gamma")
5633 insist(sum(this%binned(:, 1947 ))== 0.0000000000000000 ," 531420 has zero counts if not has_gamma")
5635 insist(sum(this%binned(:, 1948 ))== 0.0000000000000000 ," 541420 has zero counts if not has_gamma")
5637 insist(sum(this%binned(:, 1949 ))== 0.0000000000000000 ," 551420 has zero counts if not has_gamma")
5641 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.44279565642674240 , new=",new,")"
5649 insist(sum(this%binned(:, 1952 ))== 0.0000000000000000 ," 581420 has zero counts if not has_gamma")
5653 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.95774008341515993 , new=",new,")"
5656 insist(sum(this%binned(:, 1954 ))== 0.0000000000000000 ," 591421 has zero counts if not has_gamma")
5658 insist(sum(this%binned(:, 1955 ))== 0.0000000000000000 ," 601420 has zero counts if not has_gamma")
5660 insist(sum(this%binned(:, 1956 ))== 0.0000000000000000 ," 531430 has zero counts if not has_gamma")
5662 insist(sum(this%binned(:, 1957 ))== 0.0000000000000000 ," 541430 has zero counts if not has_gamma")
5664 insist(sum(this%binned(:, 1958 ))== 0.0000000000000000 ," 551430 has zero counts if not has_gamma")
5666 insist(sum(this%binned(:, 1959 ))== 0.0000000000000000 ," 561430 has zero counts if not has_gamma")
5668 insist(sum(this%binned(:, 1960 ))== 0.0000000000000000 ," 571430 has zero counts if not has_gamma")
5672 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.42580739165593329 , new=",new,")"
5677 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.32059556534630923 , new=",new,")"
5680 insist(sum(this%binned(:, 1963 ))== 0.0000000000000000 ," 601430 has zero counts if not has_gamma")
5682 insist(sum(this%binned(:, 1964 ))== 0.0000000000000000 ," 611430 has zero counts if not has_gamma")
5684 insist(sum(this%binned(:, 1965 ))== 0.0000000000000000 ," 621430 has zero counts if not has_gamma")
5686 insist(sum(this%binned(:, 1966 ))== 0.0000000000000000 ," 621431 has zero counts if not has_gamma")
5688 insist(sum(this%binned(:, 1967 ))== 0.0000000000000000 ," 531440 has zero counts if not has_gamma")
5690 insist(sum(this%binned(:, 1968 ))== 0.0000000000000000 ," 541440 has zero counts if not has_gamma")
5692 insist(sum(this%binned(:, 1969 ))== 0.0000000000000000 ," 551440 has zero counts if not has_gamma")
5694 insist(sum(this%binned(:, 1970 ))== 0.0000000000000000 ," 561440 has zero counts if not has_gamma")
5696 insist(sum(this%binned(:, 1971 ))== 0.0000000000000000 ," 571440 has zero counts if not has_gamma")
5700 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 6.3751021087238499E-002 , new=",new,")"
5710 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 2.2731732977862911E-007 , new=",new,")"
5713 insist(sum(this%binned(:, 1975 ))== 0.0000000000000000 ," 601440 has zero counts if not has_gamma")
5715 insist(sum(this%binned(:, 1976 ))== 0.0000000000000000 ," 611440 has zero counts if not has_gamma")
5717 insist(sum(this%binned(:, 1977 ))== 0.0000000000000000 ," 621440 has zero counts if not has_gamma")
5719 insist(sum(this%binned(:, 1978 ))== 0.0000000000000000 ," 541450 has zero counts if not has_gamma")
5721 insist(sum(this%binned(:, 1979 ))== 0.0000000000000000 ," 551450 has zero counts if not has_gamma")
5723 insist(sum(this%binned(:, 1980 ))== 0.0000000000000000 ," 561450 has zero counts if not has_gamma")
5725 insist(sum(this%binned(:, 1981 ))== 0.0000000000000000 ," 571450 has zero counts if not has_gamma")
5729 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.68194909056210118 , new=",new,")"
5734 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.76913395252632344 , new=",new,")"
5737 insist(sum(this%binned(:, 1984 ))== 0.0000000000000000 ," 601450 has zero counts if not has_gamma")
5739 insist(sum(this%binned(:, 1985 ))== 0.0000000000000000 ," 611450 has zero counts if not has_gamma")
5741 insist(sum(this%binned(:, 1986 ))== 0.0000000000000000 ," 621450 has zero counts if not has_gamma")
5743 insist(sum(this%binned(:, 1987 ))== 0.0000000000000000 ," 541460 has zero counts if not has_gamma")
5745 insist(sum(this%binned(:, 1988 ))== 0.0000000000000000 ," 551460 has zero counts if not has_gamma")
5747 insist(sum(this%binned(:, 1989 ))== 0.0000000000000000 ," 561460 has zero counts if not has_gamma")
5749 insist(sum(this%binned(:, 1990 ))== 0.0000000000000000 ," 571460 has zero counts if not has_gamma")
5751 insist(sum(this%binned(:, 1991 ))== 0.0000000000000000 ," 571461 has zero counts if not has_gamma")
5755 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.23345326440065151 , new=",new,")"
5758 insist(sum(this%binned(:, 1993 ))== 0.0000000000000000 ," 591460 has zero counts if not has_gamma")
5760 insist(sum(this%binned(:, 1994 ))== 0.0000000000000000 ," 601460 has zero counts if not has_gamma")
5762 insist(sum(this%binned(:, 1995 ))== 0.0000000000000000 ," 611460 has zero counts if not has_gamma")
5764 insist(sum(this%binned(:, 1996 ))== 0.0000000000000000 ," 621460 has zero counts if not has_gamma")
5766 insist(sum(this%binned(:, 1997 ))== 0.0000000000000000 ," 541470 has zero counts if not has_gamma")
5768 insist(sum(this%binned(:, 1998 ))== 0.0000000000000000 ," 551470 has zero counts if not has_gamma")
5770 insist(sum(this%binned(:, 1999 ))== 0.0000000000000000 ," 561470 has zero counts if not has_gamma")
5772 insist(sum(this%binned(:, 2000 ))== 0.0000000000000000 ," 571470 has zero counts if not has_gamma")
5774 insist(sum(this%binned(:, 2001 ))== 0.0000000000000000 ," 581470 has zero counts if not has_gamma")
5776 insist(sum(this%binned(:, 2002 ))== 0.0000000000000000 ," 591470 has zero counts if not has_gamma")
5780 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.22875016324034583 , new=",new,")"
5785 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 4.4805144281838238E-002 , new=",new,")"
5788 insist(sum(this%binned(:, 2005 ))== 0.0000000000000000 ," 621470 has zero counts if not has_gamma")
5790 insist(sum(this%binned(:, 2006 ))== 0.0000000000000000 ," 631470 has zero counts if not has_gamma")
5792 insist(sum(this%binned(:, 2007 ))== 0.0000000000000000 ," 641470 has zero counts if not has_gamma")
5794 insist(sum(this%binned(:, 2008 ))== 0.0000000000000000 ," 551480 has zero counts if not has_gamma")
5796 insist(sum(this%binned(:, 2009 ))== 0.0000000000000000 ," 561480 has zero counts if not has_gamma")
5798 insist(sum(this%binned(:, 2010 ))== 0.0000000000000000 ," 571480 has zero counts if not has_gamma")
5800 insist(sum(this%binned(:, 2011 ))== 0.0000000000000000 ," 581480 has zero counts if not has_gamma")
5807 insist(sum(this%binned(:, 2013 ))== 0.0000000000000000 ," 591481 has zero counts if not has_gamma")
5809 insist(sum(this%binned(:, 2014 ))== 0.0000000000000000 ," 601480 has zero counts if not has_gamma")
5813 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.83675692258202616 , new=",new,")"
5818 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.12653649152724911 , new=",new,")"
5821 insist(sum(this%binned(:, 2017 ))== 0.0000000000000000 ," 621480 has zero counts if not has_gamma")
5823 insist(sum(this%binned(:, 2018 ))== 0.0000000000000000 ," 551490 has zero counts if not has_gamma")
5825 insist(sum(this%binned(:, 2019 ))== 0.0000000000000000 ," 561490 has zero counts if not has_gamma")
5827 insist(sum(this%binned(:, 2020 ))== 0.0000000000000000 ," 571490 has zero counts if not has_gamma")
5829 insist(sum(this%binned(:, 2021 ))== 0.0000000000000000 ," 581490 has zero counts if not has_gamma")
5838 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.50091352746460194 , new=",new,")"
5843 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.38134159373967347 , new=",new,")"
5846 insist(sum(this%binned(:, 2025 ))== 0.0000000000000000 ," 621490 has zero counts if not has_gamma")
5848 insist(sum(this%binned(:, 2026 ))== 0.0000000000000000 ," 631490 has zero counts if not has_gamma")
5850 insist(sum(this%binned(:, 2027 ))== 0.0000000000000000 ," 641490 has zero counts if not has_gamma")
5852 insist(sum(this%binned(:, 2028 ))== 0.0000000000000000 ," 551500 has zero counts if not has_gamma")
5854 insist(sum(this%binned(:, 2029 ))== 0.0000000000000000 ," 561500 has zero counts if not has_gamma")
5856 insist(sum(this%binned(:, 2030 ))== 0.0000000000000000 ," 571500 has zero counts if not has_gamma")
5858 insist(sum(this%binned(:, 2031 ))== 0.0000000000000000 ," 581500 has zero counts if not has_gamma")
5860 insist(sum(this%binned(:, 2032 ))== 0.0000000000000000 ," 591500 has zero counts if not has_gamma")
5862 insist(sum(this%binned(:, 2033 ))== 0.0000000000000000 ," 601500 has zero counts if not has_gamma")
5866 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.90416439631137979 , new=",new,")"
5869 insist(sum(this%binned(:, 2035 ))== 0.0000000000000000 ," 621500 has zero counts if not has_gamma")
5871 insist(sum(this%binned(:, 2036 ))== 0.0000000000000000 ," 551510 has zero counts if not has_gamma")
5873 insist(sum(this%binned(:, 2037 ))== 0.0000000000000000 ," 561510 has zero counts if not has_gamma")
5875 insist(sum(this%binned(:, 2038 ))== 0.0000000000000000 ," 571510 has zero counts if not has_gamma")
5877 insist(sum(this%binned(:, 2039 ))== 0.0000000000000000 ," 581510 has zero counts if not has_gamma")
5879 insist(sum(this%binned(:, 2040 ))== 0.0000000000000000 ," 591510 has zero counts if not has_gamma")
5883 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.68210990210782507 , new=",new,")"
5888 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.28534802372394646 , new=",new,")"
5893 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.5480743712219379E-003 , new=",new,")"
5896 insist(sum(this%binned(:, 2044 ))== 0.0000000000000000 ," 631510 has zero counts if not has_gamma")
5898 insist(sum(this%binned(:, 2045 ))== 0.0000000000000000 ," 641510 has zero counts if not has_gamma")
5900 insist(sum(this%binned(:, 2046 ))== 0.0000000000000000 ," 651510 has zero counts if not has_gamma")
5902 insist(sum(this%binned(:, 2047 ))== 0.0000000000000000 ," 561520 has zero counts if not has_gamma")
5904 insist(sum(this%binned(:, 2048 ))== 0.0000000000000000 ," 571520 has zero counts if not has_gamma")
5906 insist(sum(this%binned(:, 2049 ))== 0.0000000000000000 ," 581520 has zero counts if not has_gamma")
5908 insist(sum(this%binned(:, 2050 ))== 0.0000000000000000 ," 591520 has zero counts if not has_gamma")
5910 insist(sum(this%binned(:, 2051 ))== 0.0000000000000000 ," 601520 has zero counts if not has_gamma")
5912 insist(sum(this%binned(:, 2052 ))== 0.0000000000000000 ," 611520 has zero counts if not has_gamma")
5914 insist(sum(this%binned(:, 2053 ))== 0.0000000000000000 ," 611521 has zero counts if not has_gamma")
5916 insist(sum(this%binned(:, 2054 ))== 0.0000000000000000 ," 621520 has zero counts if not has_gamma")
5920 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 2.4607505074389455E-002 , new=",new,")"
5925 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.40900613295608057 , new=",new,")"
5928 insist(sum(this%binned(:, 2057 ))== 0.0000000000000000 ," 641520 has zero counts if not has_gamma")
5930 insist(sum(this%binned(:, 2058 ))== 0.0000000000000000 ," 561530 has zero counts if not has_gamma")
5932 insist(sum(this%binned(:, 2059 ))== 0.0000000000000000 ," 571530 has zero counts if not has_gamma")
5934 insist(sum(this%binned(:, 2060 ))== 0.0000000000000000 ," 581530 has zero counts if not has_gamma")
5936 insist(sum(this%binned(:, 2061 ))== 0.0000000000000000 ," 591530 has zero counts if not has_gamma")
5938 insist(sum(this%binned(:, 2062 ))== 0.0000000000000000 ," 601530 has zero counts if not has_gamma")
5942 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.72474542147517440 , new=",new,")"
5947 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.21990841593176749 , new=",new,")"
5950 insist(sum(this%binned(:, 2065 ))== 0.0000000000000000 ," 631530 has zero counts if not has_gamma")
5952 insist(sum(this%binned(:, 2066 ))== 0.0000000000000000 ," 641530 has zero counts if not has_gamma")
5954 insist(sum(this%binned(:, 2067 ))== 0.0000000000000000 ," 651530 has zero counts if not has_gamma")
5956 insist(sum(this%binned(:, 2068 ))== 0.0000000000000000 ," 571540 has zero counts if not has_gamma")
5958 insist(sum(this%binned(:, 2069 ))== 0.0000000000000000 ," 581540 has zero counts if not has_gamma")
5960 insist(sum(this%binned(:, 2070 ))== 0.0000000000000000 ," 591540 has zero counts if not has_gamma")
5962 insist(sum(this%binned(:, 2071 ))== 0.0000000000000000 ," 601540 has zero counts if not has_gamma")
5964 insist(sum(this%binned(:, 2072 ))== 0.0000000000000000 ," 611540 has zero counts if not has_gamma")
5966 insist(sum(this%binned(:, 2073 ))== 0.0000000000000000 ," 611541 has zero counts if not has_gamma")
5968 insist(sum(this%binned(:, 2074 ))== 0.0000000000000000 ," 621540 has zero counts if not has_gamma")
5972 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.23601550649829009 , new=",new,")"
5975 insist(sum(this%binned(:, 2076 ))== 0.0000000000000000 ," 631541 has zero counts if not has_gamma")
5977 insist(sum(this%binned(:, 2077 ))== 0.0000000000000000 ," 641540 has zero counts if not has_gamma")
5979 insist(sum(this%binned(:, 2078 ))== 0.0000000000000000 ," 571550 has zero counts if not has_gamma")
5981 insist(sum(this%binned(:, 2079 ))== 0.0000000000000000 ," 581550 has zero counts if not has_gamma")
5983 insist(sum(this%binned(:, 2080 ))== 0.0000000000000000 ," 591550 has zero counts if not has_gamma")
5985 insist(sum(this%binned(:, 2081 ))== 0.0000000000000000 ," 601550 has zero counts if not has_gamma")
5987 insist(sum(this%binned(:, 2082 ))== 0.0000000000000000 ," 611550 has zero counts if not has_gamma")
5991 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.61408630543853560 , new=",new,")"
5996 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 3.0158225960669194E-002 , new=",new,")"
5999 insist(sum(this%binned(:, 2085 ))== 0.0000000000000000 ," 641550 has zero counts if not has_gamma")
6001 insist(sum(this%binned(:, 2086 ))== 0.0000000000000000 ," 641551 has zero counts if not has_gamma")
6003 insist(sum(this%binned(:, 2087 ))== 0.0000000000000000 ," 651550 has zero counts if not has_gamma")
6005 insist(sum(this%binned(:, 2088 ))== 0.0000000000000000 ," 661550 has zero counts if not has_gamma")
6007 insist(sum(this%binned(:, 2089 ))== 0.0000000000000000 ," 581560 has zero counts if not has_gamma")
6009 insist(sum(this%binned(:, 2090 ))== 0.0000000000000000 ," 591560 has zero counts if not has_gamma")
6011 insist(sum(this%binned(:, 2091 ))== 0.0000000000000000 ," 601560 has zero counts if not has_gamma")
6013 insist(sum(this%binned(:, 2092 ))== 0.0000000000000000 ," 611560 has zero counts if not has_gamma")
6017 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.17136474786326161 , new=",new,")"
6022 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.43135526533026802 , new=",new,")"
6025 insist(sum(this%binned(:, 2095 ))== 0.0000000000000000 ," 641560 has zero counts if not has_gamma")
6027 insist(sum(this%binned(:, 2096 ))== 0.0000000000000000 ," 651560 has zero counts if not has_gamma")
6029 insist(sum(this%binned(:, 2097 ))== 0.0000000000000000 ," 651561 has zero counts if not has_gamma")
6031 insist(sum(this%binned(:, 2098 ))== 0.0000000000000000 ," 661560 has zero counts if not has_gamma")
6033 insist(sum(this%binned(:, 2099 ))== 0.0000000000000000 ," 581570 has zero counts if not has_gamma")
6035 insist(sum(this%binned(:, 2100 ))== 0.0000000000000000 ," 591570 has zero counts if not has_gamma")
6037 insist(sum(this%binned(:, 2101 ))== 0.0000000000000000 ," 601570 has zero counts if not has_gamma")
6039 insist(sum(this%binned(:, 2102 ))== 0.0000000000000000 ," 611570 has zero counts if not has_gamma")
6041 insist(sum(this%binned(:, 2103 ))== 0.0000000000000000 ," 621570 has zero counts if not has_gamma")
6045 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.38807598956565947 , new=",new,")"
6048 insist(sum(this%binned(:, 2105 ))== 0.0000000000000000 ," 641570 has zero counts if not has_gamma")
6050 insist(sum(this%binned(:, 2106 ))== 0.0000000000000000 ," 651570 has zero counts if not has_gamma")
6052 insist(sum(this%binned(:, 2107 ))== 0.0000000000000000 ," 661570 has zero counts if not has_gamma")
6054 insist(sum(this%binned(:, 2108 ))== 0.0000000000000000 ," 591580 has zero counts if not has_gamma")
6056 insist(sum(this%binned(:, 2109 ))== 0.0000000000000000 ," 601580 has zero counts if not has_gamma")
6058 insist(sum(this%binned(:, 2110 ))== 0.0000000000000000 ," 611580 has zero counts if not has_gamma")
6060 insist(sum(this%binned(:, 2111 ))== 0.0000000000000000 ," 621580 has zero counts if not has_gamma")
6062 insist(sum(this%binned(:, 2112 ))== 0.0000000000000000 ," 631580 has zero counts if not has_gamma")
6064 insist(sum(this%binned(:, 2113 ))== 0.0000000000000000 ," 641580 has zero counts if not has_gamma")
6066 insist(sum(this%binned(:, 2114 ))== 0.0000000000000000 ," 651580 has zero counts if not has_gamma")
6068 insist(sum(this%binned(:, 2115 ))== 0.0000000000000000 ," 651581 has zero counts if not has_gamma")
6070 insist(sum(this%binned(:, 2116 ))== 0.0000000000000000 ," 661580 has zero counts if not has_gamma")
6072 insist(sum(this%binned(:, 2117 ))== 0.0000000000000000 ," 591590 has zero counts if not has_gamma")
6074 insist(sum(this%binned(:, 2118 ))== 0.0000000000000000 ," 601590 has zero counts if not has_gamma")
6076 insist(sum(this%binned(:, 2119 ))== 0.0000000000000000 ," 611590 has zero counts if not has_gamma")
6078 insist(sum(this%binned(:, 2120 ))== 0.0000000000000000 ," 621590 has zero counts if not has_gamma")
6087 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.29955126869737125 , new=",new,")"
6090 insist(sum(this%binned(:, 2123 ))== 0.0000000000000000 ," 651590 has zero counts if not has_gamma")
6092 insist(sum(this%binned(:, 2124 ))== 0.0000000000000000 ," 661590 has zero counts if not has_gamma")
6094 insist(sum(this%binned(:, 2125 ))== 0.0000000000000000 ," 671590 has zero counts if not has_gamma")
6096 insist(sum(this%binned(:, 2126 ))== 0.0000000000000000 ," 671591 has zero counts if not has_gamma")
6098 insist(sum(this%binned(:, 2127 ))== 0.0000000000000000 ," 601600 has zero counts if not has_gamma")
6100 insist(sum(this%binned(:, 2128 ))== 0.0000000000000000 ," 611600 has zero counts if not has_gamma")
6102 insist(sum(this%binned(:, 2129 ))== 0.0000000000000000 ," 621600 has zero counts if not has_gamma")
6104 insist(sum(this%binned(:, 2130 ))== 0.0000000000000000 ," 631600 has zero counts if not has_gamma")
6106 insist(sum(this%binned(:, 2131 ))== 0.0000000000000000 ," 641600 has zero counts if not has_gamma")
6110 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.19677743286810334 , new=",new,")"
6113 insist(sum(this%binned(:, 2133 ))== 0.0000000000000000 ," 661600 has zero counts if not has_gamma")
6115 insist(sum(this%binned(:, 2134 ))== 0.0000000000000000 ," 601610 has zero counts if not has_gamma")
6117 insist(sum(this%binned(:, 2135 ))== 0.0000000000000000 ," 611610 has zero counts if not has_gamma")
6119 insist(sum(this%binned(:, 2136 ))== 0.0000000000000000 ," 621610 has zero counts if not has_gamma")
6121 insist(sum(this%binned(:, 2137 ))== 0.0000000000000000 ," 631610 has zero counts if not has_gamma")
6125 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.62205631423809715 , new=",new,")"
6130 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.13716866353510065 , new=",new,")"
6133 insist(sum(this%binned(:, 2140 ))== 0.0000000000000000 ," 661610 has zero counts if not has_gamma")
6135 insist(sum(this%binned(:, 2141 ))== 0.0000000000000000 ," 671610 has zero counts if not has_gamma")
6137 insist(sum(this%binned(:, 2142 ))== 0.0000000000000000 ," 671611 has zero counts if not has_gamma")
6139 insist(sum(this%binned(:, 2143 ))== 0.0000000000000000 ," 681610 has zero counts if not has_gamma")
6141 insist(sum(this%binned(:, 2144 ))== 0.0000000000000000 ," 611620 has zero counts if not has_gamma")
6143 insist(sum(this%binned(:, 2145 ))== 0.0000000000000000 ," 621620 has zero counts if not has_gamma")
6145 insist(sum(this%binned(:, 2146 ))== 0.0000000000000000 ," 631620 has zero counts if not has_gamma")
6149 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.33008241509130104 , new=",new,")"
6154 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.49386870265293037 , new=",new,")"
6157 insist(sum(this%binned(:, 2149 ))== 0.0000000000000000 ," 661620 has zero counts if not has_gamma")
6159 insist(sum(this%binned(:, 2150 ))== 0.0000000000000000 ," 671620 has zero counts if not has_gamma")
6161 insist(sum(this%binned(:, 2151 ))== 0.0000000000000000 ," 671621 has zero counts if not has_gamma")
6163 insist(sum(this%binned(:, 2152 ))== 0.0000000000000000 ," 681620 has zero counts if not has_gamma")
6165 insist(sum(this%binned(:, 2153 ))== 0.0000000000000000 ," 611630 has zero counts if not has_gamma")
6167 insist(sum(this%binned(:, 2154 ))== 0.0000000000000000 ," 621630 has zero counts if not has_gamma")
6169 insist(sum(this%binned(:, 2155 ))== 0.0000000000000000 ," 631630 has zero counts if not has_gamma")
6171 insist(sum(this%binned(:, 2156 ))== 0.0000000000000000 ," 641630 has zero counts if not has_gamma")
6173 insist(sum(this%binned(:, 2157 ))== 0.0000000000000000 ," 651630 has zero counts if not has_gamma")
6175 insist(sum(this%binned(:, 2158 ))== 0.0000000000000000 ," 661630 has zero counts if not has_gamma")
6177 insist(sum(this%binned(:, 2159 ))== 0.0000000000000000 ," 671630 has zero counts if not has_gamma")
6179 insist(sum(this%binned(:, 2160 ))== 0.0000000000000000 ," 671631 has zero counts if not has_gamma")
6181 insist(sum(this%binned(:, 2161 ))== 0.0000000000000000 ," 681630 has zero counts if not has_gamma")
6183 insist(sum(this%binned(:, 2162 ))== 0.0000000000000000 ," 621640 has zero counts if not has_gamma")
6185 insist(sum(this%binned(:, 2163 ))== 0.0000000000000000 ," 631640 has zero counts if not has_gamma")
6187 insist(sum(this%binned(:, 2164 ))== 0.0000000000000000 ," 641640 has zero counts if not has_gamma")
6191 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.90820332488854727 , new=",new,")"
6194 insist(sum(this%binned(:, 2166 ))== 0.0000000000000000 ," 661640 has zero counts if not has_gamma")
6196 insist(sum(this%binned(:, 2167 ))== 0.0000000000000000 ," 671640 has zero counts if not has_gamma")
6198 insist(sum(this%binned(:, 2168 ))== 0.0000000000000000 ," 671641 has zero counts if not has_gamma")
6200 insist(sum(this%binned(:, 2169 ))== 0.0000000000000000 ," 681640 has zero counts if not has_gamma")
6202 insist(sum(this%binned(:, 2170 ))== 0.0000000000000000 ," 621650 has zero counts if not has_gamma")
6204 insist(sum(this%binned(:, 2171 ))== 0.0000000000000000 ," 631650 has zero counts if not has_gamma")
6206 insist(sum(this%binned(:, 2172 ))== 0.0000000000000000 ," 641650 has zero counts if not has_gamma")
6208 insist(sum(this%binned(:, 2173 ))== 0.0000000000000000 ," 651650 has zero counts if not has_gamma")
6212 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.47595028997703481 , new=",new,")"
6217 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.5781428807830995E-004 , new=",new,")"
6220 insist(sum(this%binned(:, 2176 ))== 0.0000000000000000 ," 671650 has zero counts if not has_gamma")
6222 insist(sum(this%binned(:, 2177 ))== 0.0000000000000000 ," 681650 has zero counts if not has_gamma")
6224 insist(sum(this%binned(:, 2178 ))== 0.0000000000000000 ," 691650 has zero counts if not has_gamma")
6226 insist(sum(this%binned(:, 2179 ))== 0.0000000000000000 ," 631660 has zero counts if not has_gamma")
6228 insist(sum(this%binned(:, 2180 ))== 0.0000000000000000 ," 641660 has zero counts if not has_gamma")
6230 insist(sum(this%binned(:, 2181 ))== 0.0000000000000000 ," 651660 has zero counts if not has_gamma")
6234 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 9.8151209522421748E-002 , new=",new,")"
6239 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.76381537700397173 , new=",new,")"
6244 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 6.9905578106101993E-003 , new=",new,")"
6247 insist(sum(this%binned(:, 2185 ))== 0.0000000000000000 ," 681660 has zero counts if not has_gamma")
6249 insist(sum(this%binned(:, 2186 ))== 0.0000000000000000 ," 691660 has zero counts if not has_gamma")
6251 insist(sum(this%binned(:, 2187 ))== 0.0000000000000000 ," 701660 has zero counts if not has_gamma")
6253 insist(sum(this%binned(:, 2188 ))== 0.0000000000000000 ," 631670 has zero counts if not has_gamma")
6255 insist(sum(this%binned(:, 2189 ))== 0.0000000000000000 ," 641670 has zero counts if not has_gamma")
6257 insist(sum(this%binned(:, 2190 ))== 0.0000000000000000 ," 651670 has zero counts if not has_gamma")
6259 insist(sum(this%binned(:, 2191 ))== 0.0000000000000000 ," 661670 has zero counts if not has_gamma")
6261 insist(sum(this%binned(:, 2192 ))== 0.0000000000000000 ," 671670 has zero counts if not has_gamma")
6263 insist(sum(this%binned(:, 2193 ))== 0.0000000000000000 ," 681670 has zero counts if not has_gamma")
6265 insist(sum(this%binned(:, 2194 ))== 0.0000000000000000 ," 681671 has zero counts if not has_gamma")
6267 insist(sum(this%binned(:, 2195 ))== 0.0000000000000000 ," 691670 has zero counts if not has_gamma")
6269 insist(sum(this%binned(:, 2196 ))== 0.0000000000000000 ," 701670 has zero counts if not has_gamma")
6271 insist(sum(this%binned(:, 2197 ))== 0.0000000000000000 ," 641680 has zero counts if not has_gamma")
6273 insist(sum(this%binned(:, 2198 ))== 0.0000000000000000 ," 651680 has zero counts if not has_gamma")
6275 insist(sum(this%binned(:, 2199 ))== 0.0000000000000000 ," 661680 has zero counts if not has_gamma")
6277 insist(sum(this%binned(:, 2200 ))== 0.0000000000000000 ," 671680 has zero counts if not has_gamma")
6279 insist(sum(this%binned(:, 2201 ))== 0.0000000000000000 ," 681680 has zero counts if not has_gamma")
6281 insist(sum(this%binned(:, 2202 ))== 0.0000000000000000 ," 691680 has zero counts if not has_gamma")
6283 insist(sum(this%binned(:, 2203 ))== 0.0000000000000000 ," 701680 has zero counts if not has_gamma")
6285 insist(sum(this%binned(:, 2204 ))== 0.0000000000000000 ," 641690 has zero counts if not has_gamma")
6287 insist(sum(this%binned(:, 2205 ))== 0.0000000000000000 ," 651690 has zero counts if not has_gamma")
6289 insist(sum(this%binned(:, 2206 ))== 0.0000000000000000 ," 661690 has zero counts if not has_gamma")
6291 insist(sum(this%binned(:, 2207 ))== 0.0000000000000000 ," 671690 has zero counts if not has_gamma")
6295 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 8.2603648633007620E-002 , new=",new,")"
6298 insist(sum(this%binned(:, 2209 ))== 0.0000000000000000 ," 691690 has zero counts if not has_gamma")
6300 insist(sum(this%binned(:, 2210 ))== 0.0000000000000000 ," 701690 has zero counts if not has_gamma")
6302 insist(sum(this%binned(:, 2211 ))== 0.0000000000000000 ," 701691 has zero counts if not has_gamma")
6304 insist(sum(this%binned(:, 2212 ))== 0.0000000000000000 ," 711690 has zero counts if not has_gamma")
6306 insist(sum(this%binned(:, 2213 ))== 0.0000000000000000 ," 711691 has zero counts if not has_gamma")
6308 insist(sum(this%binned(:, 2214 ))== 0.0000000000000000 ," 651700 has zero counts if not has_gamma")
6310 insist(sum(this%binned(:, 2215 ))== 0.0000000000000000 ," 661700 has zero counts if not has_gamma")
6312 insist(sum(this%binned(:, 2216 ))== 0.0000000000000000 ," 671700 has zero counts if not has_gamma")
6314 insist(sum(this%binned(:, 2217 ))== 0.0000000000000000 ," 671701 has zero counts if not has_gamma")
6316 insist(sum(this%binned(:, 2218 ))== 0.0000000000000000 ," 681700 has zero counts if not has_gamma")
6320 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.32198517535952448 , new=",new,")"
6323 insist(sum(this%binned(:, 2220 ))== 0.0000000000000000 ," 701700 has zero counts if not has_gamma")
6325 insist(sum(this%binned(:, 2221 ))== 0.0000000000000000 ," 651710 has zero counts if not has_gamma")
6327 insist(sum(this%binned(:, 2222 ))== 0.0000000000000000 ," 661710 has zero counts if not has_gamma")
6329 insist(sum(this%binned(:, 2223 ))== 0.0000000000000000 ," 671710 has zero counts if not has_gamma")
6333 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.38508704871806548 , new=",new,")"
6338 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 1.2421613035298590E-002 , new=",new,")"
6341 insist(sum(this%binned(:, 2226 ))== 0.0000000000000000 ," 701710 has zero counts if not has_gamma")
6343 insist(sum(this%binned(:, 2227 ))== 0.0000000000000000 ," 711710 has zero counts if not has_gamma")
6345 insist(sum(this%binned(:, 2228 ))== 0.0000000000000000 ," 711711 has zero counts if not has_gamma")
6347 insist(sum(this%binned(:, 2229 ))== 0.0000000000000000 ," 721710 has zero counts if not has_gamma")
6349 insist(sum(this%binned(:, 2230 ))== 0.0000000000000000 ," 661720 has zero counts if not has_gamma")
6351 insist(sum(this%binned(:, 2231 ))== 0.0000000000000000 ," 671720 has zero counts if not has_gamma")
6355 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 7.9573185132163951E-002 , new=",new,")"
6360 write(msg,"(a,es18.6,a)")" total counts are conserved (ref= 0.54836751755446267 , new=",new,")"
6363 insist(sum(this%binned(:, 2234 ))== 0.0000000000000000 ," 701720 has zero counts if not has_gamma")
6365 insist(sum(this%binned(:, 2235 ))== 0.0000000000000000 ," 711720 has zero counts if not has_gamma")
6367 insist(sum(this%binned(:, 2236 ))== 0.0000000000000000 ," 711721 has zero counts if not has_gamma")