diff --git a/NNLOJET_ATLAS_Z0J_13TEV_PT_LL/ATLAS_Z0J_13TEV_PT_LL.yaml b/NNLOJET_ATLAS_Z0J_13TEV_PT_LL/ATLAS_Z0J_13TEV_PT_LL.yaml new file mode 100644 index 00000000..178caed6 --- /dev/null +++ b/NNLOJET_ATLAS_Z0J_13TEV_PT_LL/ATLAS_Z0J_13TEV_PT_LL.yaml @@ -0,0 +1,217 @@ +# arXiv number: 1912.02844.pdf, hepdata entry: https://www.hepdata.net/record/ins1768911 (tables: [5]) + +runname: ATLAS_Z0J_13TEV_PT_LL + +process: + proc: ZJ + sqrts: 13000.0 + jet: antikt[0.5] + +pdf: NNPDF40_nnlo_as_01180 + +techcut: 1e-07 + +histograms: + - name: pT + observable: ptz + bins: + - 30.0 + - 33.0 + - 36.0 + - 39.0 + - 42.0 + - 45.0 + - 48.0 + - 51.0 + - 54.0 + - 57.0 + - 61.0 + - 65.0 + - 70.0 + - 75.0 + - 80.0 + - 85.0 + - 95.0 + - 105.0 + - 125.0 + - 150.0 + - 175.0 + - 200.0 + - 250.0 + - 300.0 + - 350.0 + - 400.0 + - 470.0 + - 550.0 + - 650.0 + - 900.0 + - 2500.0 + +multi_channel: 1 + +channels: + LO: LO + R: R + V: V + RRa_1: "163 177 191 205 206 245 246 285 301" # channel: [10,0] + RRa_2: "164 178 192 207 208 247 248 286 302" # channel: [20,0] + RRa_3: "165 179 193 217 218 257 258 287 303" # channel: [10,-10] + RRa_4: "166 180 194 219 220 259 260 288 304" # channel: [20,-20] + RRa_5: "167 181 195 221 222 261 262 291 307" # channel: [0,10] + RRa_6: "168 182 196 223 224 263 264 292 308" # channel: [0,20] + RRa_7: "169 183 197 225 226 265 266 293 309" # channel: [-10,10] + RRa_8: "170 184 198 227 228 267 268 294 310" # channel: [-20,20] + RRa_9: "173 187 201 229 230 269 270 295 311" # channel: [0,-10] + RRa_10: "174 188 202 231 232 271 272 296 312" # channel: [0,-20] + RRa_11: "175 189 203 233 234 273 274 297 313" # channel: [-10,0] + RRa_12: "176 190 204 235 236 275 276 298 314" # channel: [-20,0] + RRa_13: "171 172 185 186 199 200" # channel: [0,0] + RRa_14: "209 249" # channel: [11,-11] + RRa_15: "210 250" # channel: [11,-21] + RRa_16: "211 251" # channel: [21,-11] + RRa_17: "212 252" # channel: [21,-21] + RRa_18: "213 253" # channel: [11,11] + RRa_19: "214 254" # channel: [11,21] + RRa_20: "215 255" # channel: [21,11] + RRa_21: "216 256" # channel: [21,21] + RRa_22: "237 277" # channel: [-11,-11] + RRa_23: "238 278" # channel: [-11,-21] + RRa_24: "239 279" # channel: [-21,-11] + RRa_25: "240 280" # channel: [-21,-21] + RRa_26: "241 281" # channel: [-11,11] + RRa_27: "242 282" # channel: [-11,21] + RRa_28: "243 283" # channel: [-21,11] + RRa_29: "244 284" # channel: [-21,21] + RRa_30: "289 305" # channel: [10,10] + RRa_31: "290 306" # channel: [20,20] + RRa_32: "299 315" # channel: [-10,-10] + RRa_33: "300 316" # channel: [-20,-20] + RRb_1: "163 177 191 205 206 245 246 285 301" # channel: [10,0] + RRb_2: "164 178 192 207 208 247 248 286 302" # channel: [20,0] + RRb_3: "165 179 193 217 218 257 258 287 303" # channel: [10,-10] + RRb_4: "166 180 194 219 220 259 260 288 304" # channel: [20,-20] + RRb_5: "167 181 195 221 222 261 262 291 307" # channel: [0,10] + RRb_6: "168 182 196 223 224 263 264 292 308" # channel: [0,20] + RRb_7: "169 183 197 225 226 265 266 293 309" # channel: [-10,10] + RRb_8: "170 184 198 227 228 267 268 294 310" # channel: [-20,20] + RRb_9: "173 187 201 229 230 269 270 295 311" # channel: [0,-10] + RRb_10: "174 188 202 231 232 271 272 296 312" # channel: [0,-20] + RRb_11: "175 189 203 233 234 273 274 297 313" # channel: [-10,0] + RRb_12: "176 190 204 235 236 275 276 298 314" # channel: [-20,0] + RRb_13: "171 172 185 186 199 200" # channel: [0,0] + RRb_14: "209 249" # channel: [11,-11] + RRb_15: "210 250" # channel: [11,-21] + RRb_16: "211 251" # channel: [21,-11] + RRb_17: "212 252" # channel: [21,-21] + RRb_18: "213 253" # channel: [11,11] + RRb_19: "214 254" # channel: [11,21] + RRb_20: "215 255" # channel: [21,11] + RRb_21: "216 256" # channel: [21,21] + RRb_22: "237 277" # channel: [-11,-11] + RRb_23: "238 278" # channel: [-11,-21] + RRb_24: "239 279" # channel: [-21,-11] + RRb_25: "240 280" # channel: [-21,-21] + RRb_26: "241 281" # channel: [-11,11] + RRb_27: "242 282" # channel: [-11,21] + RRb_28: "243 283" # channel: [-21,11] + RRb_29: "244 284" # channel: [-21,21] + RRb_30: "289 305" # channel: [10,10] + RRb_31: "290 306" # channel: [20,20] + RRb_32: "299 315" # channel: [-10,-10] + RRb_33: "300 316" # channel: [-20,-20] + RV_1: "317 347 377 391 405 435 481 482 521 522 561 562 581 597 613" # channel: [10,0] + RV_2: "318 348 378 392 406 436 483 484 523 524 563 564 582 598 614" # channel: [20,0] + RV_3: "319 349 379 393 407 437 457 458 497 498 537 538 569 585 601" # channel: [10,-10] + RV_4: "320 350 380 394 408 438 459 460 499 500 539 540 570 586 602" # channel: [20,-20] + RV_5: "321 351 381 395 409 439 473 474 513 514 553 554 577 593 609" # channel: [0,10] + RV_6: "322 352 382 396 410 440 475 476 515 516 555 556 578 594 610" # channel: [0,20] + RV_7: "323 353 383 397 411 441 461 462 501 502 541 542 573 589 605" # channel: [-10,10] + RV_8: "324 354 384 398 412 442 463 464 503 504 543 544 574 590 606" # channel: [-20,20] + RV_9: "327 357 387 401 415 445 477 478 517 518 557 558 579 595 611" # channel: [0,-10] + RV_10: "328 358 388 402 416 446 479 480 519 520 559 560 580 596 612" # channel: [0,-20] + RV_11: "329 359 389 403 417 447 485 486 525 526 565 566 583 599 615" # channel: [-10,0] + RV_12: "330 360 390 404 418 448 487 488 527 528 567 568 584 600 616" # channel: [-20,0] + RV_13: "325 326 355 356 385 386 399 400 413 414 443 444" # channel: [0,0] + RV_14: "335 365 423 449 489 529" # channel: [11,-11] + RV_15: "336 366 424 450 490 530" # channel: [11,-21] + RV_16: "337 367 425 451 491 531" # channel: [21,-11] + RV_17: "338 368 426 452 492 532" # channel: [21,-21] + RV_18: "331 361 419 453 493 533" # channel: [11,11] + RV_19: "332 362 420 454 494 534" # channel: [11,21] + RV_20: "333 363 421 455 495 535" # channel: [21,11] + RV_21: "334 364 422 456 496 536" # channel: [21,21] + RV_22: "343 373 431 465 505 545" # channel: [-11,-11] + RV_23: "344 374 432 466 506 546" # channel: [-11,-21] + RV_24: "345 375 433 467 507 547" # channel: [-21,-11] + RV_25: "346 376 434 468 508 548" # channel: [-21,-21] + RV_26: "339 369 427 469 509 549" # channel: [-11,11] + RV_27: "340 370 428 470 510 550" # channel: [-11,21] + RV_28: "341 371 429 471 511 551" # channel: [-21,11] + RV_29: "342 372 430 472 512 552" # channel: [-21,21] + RV_30: "571 587 603" # channel: [10,10] + RV_31: "572 588 604" # channel: [20,20] + RV_32: "575 591 607" # channel: [-10,-10] + RV_33: "576 592 608" # channel: [-20,-20] + VV_1: "617 651 685 719 731 733 769 781 783 819 831 833" # channel: [10,0] + VV_2: "618 652 686 720 732 734 770 782 784 820 832 834" # channel: [20,0] + VV_3: "619 653 687 721 771 821" # channel: [10,-10] + VV_4: "620 654 688 722 772 822" # channel: [20,-20] + VV_5: "621 655 689 723 739 741 773 789 791 823 839 841" # channel: [0,10] + VV_6: "622 656 690 724 740 742 774 790 792 824 840 842" # channel: [0,20] + VV_7: "623 657 691 725 775 825" # channel: [-10,10] + VV_8: "624 658 692 726 776 826" # channel: [-20,20] + VV_9: "625 659 693 727 743 745 777 793 795 827 843 845" # channel: [0,-10] + VV_10: "626 660 694 728 744 746 778 794 796 828 844 846" # channel: [0,-20] + VV_11: "627 661 695 729 735 737 779 785 787 829 835 837" # channel: [-10,0] + VV_12: "628 662 696 730 736 738 780 786 788 830 836 838" # channel: [-20,0] + VV_13: "629 630 663 664 697 698 747 748 797 798 847 848" # channel: [0,0] + VV_14: "635 669 703 753 803 853" # channel: [11,-11] + VV_15: "636 670 704 754 804 854" # channel: [11,-21] + VV_16: "637 671 705 755 805 855" # channel: [21,-11] + VV_17: "638 672 706 756 806 856" # channel: [21,-21] + VV_18: "631 665 699 749 799 849" # channel: [11,11] + VV_19: "632 666 700 750 800 850" # channel: [11,21] + VV_20: "633 667 701 751 801 851" # channel: [21,11] + VV_21: "634 668 702 752 802 852" # channel: [21,21] + VV_22: "645 679 713 763 813 863" # channel: [-11,-11] + VV_23: "646 680 714 764 814 864" # channel: [-11,-21] + VV_24: "647 681 715 765 815 865" # channel: [-21,-11] + VV_25: "648 682 716 766 816 866" # channel: [-21,-21] + VV_26: "641 675 709 759 809 859" # channel: [-11,11] + VV_27: "642 676 710 760 810 860" # channel: [-11,21] + VV_28: "643 677 711 761 811 861" # channel: [-21,11] + VV_29: "644 678 712 762 812 862" # channel: [-21,21] + VV_30: "639 673 707 757 807 857" # channel: [10,10] + VV_31: "640 674 708 758 808 858" # channel: [20,20] + VV_32: "649 683 717 767 817 867" # channel: [-10,-10] + VV_33: "650 684 718 768 818 868" # channel: [-20,-20] + +parameters: + MASS[Z]: 91.1876 + MASS[W]: 80.379 + WIDTH[Z]: 2.4952 + WIDTH[W]: 2.085 + SCHEME[alpha]: Gmu + GF: 1.1663787d-5 + +selectors: + - observable: ptz + min: 30.0 + max: 2500.0 + - observable: mll + min: 66.0 + max: 116.0 + - observable: abs_ylp + min: 0.0 + max: 2.5 + - observable: abs_ylm + min: 0.0 + max: 2.5 + - observable: ptlp + min: 27.0 + - observable: ptlm + min: 27.0 + +scales: + mur: etz + muf: etz diff --git a/NNLOJET_ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT-lowpt.yaml b/NNLOJET_ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT-lowpt.yaml new file mode 100644 index 00000000..2a836b08 --- /dev/null +++ b/NNLOJET_ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT-lowpt.yaml @@ -0,0 +1,245 @@ +# arXiv number: 1512.02192, hepdata entry: https://www.hepdata.net/record/ins1408516 +# Mass distribution (tables: [35, 36, 37, 38, 40]) +# Rapidity distribution (tables: [29, 30, 31, 32, 33, 34]) +# The fine-grained pt distribution: table 27 + +runname: ATLAS_Z0J_8TEV_PT + +process: + proc: ZJ + sqrts: 8000.0 + +pdf: NNPDF40_nnlo_as_01180 + +techcut: 1e-07 + +histograms: + - name: pT_mll_46_66 + observable: ptz + bins: &lowpt + - 2.0 + - 4.0 + - 6.0 + - 8.0 + - 10.0 + - 13.0 + - 16.0 + - 20.0 + - 25.0 + - 30.0 + extra_selectors: + - observable: mll + min: 46.0 + max: 66.0 + + - name: pT_mll_116_150 + observable: ptz + bins: *lowpt + extra_selectors: + - observable: mll + min: 116.0 + max: 150.0 + + - name: pt_mll_66_116_y_0.0_0.4 + observable: ptz + bins: *lowpt + extra_selectors: + - observable: abs_yz + min: 0.0 + max: 0.4 + - ¢ral_mll + observable: mll + min: 66.0 + max: 116.0 + + - name: pt_mll_66_116_y_0.4_0.8 + observable: ptz + bins: *lowpt + extra_selectors: + - observable: abs_yz + min: 0.4 + max: 0.8 + - *central_mll + + - name: pt_mll_66_116_y_0.8_1.2 + observable: ptz + bins: *lowpt + extra_selectors: + - observable: abs_yz + min: 0.8 + max: 1.2 + - *central_mll + + - name: pt_mll_66_116_y_1.2_1.6 + observable: ptz + bins: *lowpt + extra_selectors: + - observable: abs_yz + min: 1.2 + max: 1.6 + - *central_mll + + - name: pt_mll_66_116_y_1.6_2.0 + observable: ptz + bins: *lowpt + extra_selectors: + - observable: abs_yz + min: 1.6 + max: 2.0 + - *central_mll + + - name: pt_mll_66_116_y_2.0_2.4 + observable: ptz + bins: *lowpt + extra_selectors: + - observable: abs_yz + min: 2.0 + max: 2.4 + - *central_mll + + - name: pt_mll_66_116_y_0.0_2.4 + observable: ptz + extra_selectors: + - observable: abs_yz + min: 0.0 + max: 2.4 + - *central_mll + bins: *lowpt + +multi_channel: 0 + +parameters: + MASS[Z]: 91.1876 + MASS[W]: 80.379 + WIDTH[Z]: 2.4952 + WIDTH[W]: 2.085 + SCHEME[alpha]: Gmu + GF: 1.1663787d-5 + +selectors: + - observable: ptz + min: 2.0 + max: 30.0 + - observable: abs_yz + min: 0.0 + max: 2.4 + - observable: mll + min: 46.0 + max: 150.0 + - observable: abs_ylp + max: 2.4 + - observable: abs_ylm + max: 2.4 + - observable: ptlp + min: 20.0 + - observable: ptlm + min: 20.0 + +scales: + mur: etz + muf: etz + +channels: + LO: LO + R: R + V: V + RRa_1: "163 177 191 205 206 245 246 285 301" # channel: [10,0] + RRa_2: "164 178 192 207 208 247 248 286 302" # channel: [20,0] + RRa_3: "165 179 193 217 218 257 258 287 303" # channel: [10,-10] + RRa_4: "166 180 194 219 220 259 260 288 304" # channel: [20,-20] + RRa_5: "167 181 195 221 222 261 262 291 307" # channel: [0,10] + RRa_6: "168 182 196 223 224 263 264 292 308" # channel: [0,20] + RRa_7: "169 183 197 225 226 265 266 293 309" # channel: [-10,10] + RRa_8: "170 184 198 227 228 267 268 294 310" # channel: [-20,20] + RRa_9: "173 187 201 229 230 269 270 295 311" # channel: [0,-10] + RRa_10: "174 188 202 231 232 271 272 296 312" # channel: [0,-20] + RRa_11: "175 189 203 233 234 273 274 297 313" # channel: [-10,0] + RRa_12: "176 190 204 235 236 275 276 298 314" # channel: [-20,0] + RRa_13: "171 172 185 186 199 200" # channel: [0,0] + RRa_14: "209 249" # channel: [11,-11] + RRa_15: "210 250" # channel: [11,-21] + RRa_16: "211 251" # channel: [21,-11] + RRa_17: "212 252" # channel: [21,-21] + RRa_18: "213 253" # channel: [11,11] + RRa_19: "214 254" # channel: [11,21] + RRa_20: "215 255" # channel: [21,11] + RRa_21: "216 256" # channel: [21,21] + RRa_22: "237 277" # channel: [-11,-11] + RRa_23: "238 278" # channel: [-11,-21] + RRa_24: "239 279" # channel: [-21,-11] + RRa_25: "240 280" # channel: [-21,-21] + RRa_26: "241 281" # channel: [-11,11] + RRa_27: "242 282" # channel: [-11,21] + RRa_28: "243 283" # channel: [-21,11] + RRa_29: "244 284" # channel: [-21,21] + RRa_30: "289 305" # channel: [10,10] + RRa_31: "290 306" # channel: [20,20] + RRa_32: "299 315" # channel: [-10,-10] + RRa_33: "300 316" # channel: [-20,-20] + RRb_1: "163 177 191 205 206 245 246 285 301" # channel: [10,0] + RRb_2: "164 178 192 207 208 247 248 286 302" # channel: [20,0] + RRb_3: "165 179 193 217 218 257 258 287 303" # channel: [10,-10] + RRb_4: "166 180 194 219 220 259 260 288 304" # channel: [20,-20] + RRb_5: "167 181 195 221 222 261 262 291 307" # channel: [0,10] + RRb_6: "168 182 196 223 224 263 264 292 308" # channel: [0,20] + RRb_7: "169 183 197 225 226 265 266 293 309" # channel: [-10,10] + RRb_8: "170 184 198 227 228 267 268 294 310" # channel: [-20,20] + RRb_9: "173 187 201 229 230 269 270 295 311" # channel: [0,-10] + RRb_10: "174 188 202 231 232 271 272 296 312" # channel: [0,-20] + RRb_11: "175 189 203 233 234 273 274 297 313" # channel: [-10,0] + RRb_12: "176 190 204 235 236 275 276 298 314" # channel: [-20,0] + RRb_13: "171 172 185 186 199 200" # channel: [0,0] + RRb_14: "209 249" # channel: [11,-11] + RRb_15: "210 250" # channel: [11,-21] + RRb_16: "211 251" # channel: [21,-11] + RRb_17: "212 252" # channel: [21,-21] + RRb_18: "213 253" # channel: [11,11] + RRb_19: "214 254" # channel: [11,21] + RRb_20: "215 255" # channel: [21,11] + RRb_21: "216 256" # channel: [21,21] + RRb_22: "237 277" # channel: [-11,-11] + RRb_23: "238 278" # channel: [-11,-21] + RRb_24: "239 279" # channel: [-21,-11] + RRb_25: "240 280" # channel: [-21,-21] + RRb_26: "241 281" # channel: [-11,11] + RRb_27: "242 282" # channel: [-11,21] + RRb_28: "243 283" # channel: [-21,11] + RRb_29: "244 284" # channel: [-21,21] + RRb_30: "289 305" # channel: [10,10] + RRb_31: "290 306" # channel: [20,20] + RRb_32: "299 315" # channel: [-10,-10] + RRb_33: "300 316" # channel: [-20,-20] + RV_1: "317 347 377 391 405 435 481 482 521 522 561 562 581 597 613" # channel: [10,0] + RV_2: "318 348 378 392 406 436 483 484 523 524 563 564 582 598 614" # channel: [20,0] + RV_3: "319 349 379 393 407 437 457 458 497 498 537 538 569 585 601" # channel: [10,-10] + RV_4: "320 350 380 394 408 438 459 460 499 500 539 540 570 586 602" # channel: [20,-20] + RV_5: "321 351 381 395 409 439 473 474 513 514 553 554 577 593 609" # channel: [0,10] + RV_6: "322 352 382 396 410 440 475 476 515 516 555 556 578 594 610" # channel: [0,20] + RV_7: "323 353 383 397 411 441 461 462 501 502 541 542 573 589 605" # channel: [-10,10] + RV_8: "324 354 384 398 412 442 463 464 503 504 543 544 574 590 606" # channel: [-20,20] + RV_9: "327 357 387 401 415 445 477 478 517 518 557 558 579 595 611" # channel: [0,-10] + RV_10: "328 358 388 402 416 446 479 480 519 520 559 560 580 596 612" # channel: [0,-20] + RV_11: "329 359 389 403 417 447 485 486 525 526 565 566 583 599 615" # channel: [-10,0] + RV_12: "330 360 390 404 418 448 487 488 527 528 567 568 584 600 616" # channel: [-20,0] + RV_13: "325 326 355 356 385 386 399 400 413 414 443 444" # channel: [0,0] + RV_14: "335 365 423 449 489 529" # channel: [11,-11] + RV_15: "336 366 424 450 490 530" # channel: [11,-21] + RV_16: "337 367 425 451 491 531" # channel: [21,-11] + RV_17: "338 368 426 452 492 532" # channel: [21,-21] + RV_18: "331 361 419 453 493 533" # channel: [11,11] + RV_19: "332 362 420 454 494 534" # channel: [11,21] + RV_20: "333 363 421 455 495 535" # channel: [21,11] + RV_21: "334 364 422 456 496 536" # channel: [21,21] + RV_22: "343 373 431 465 505 545" # channel: [-11,-11] + RV_23: "344 374 432 466 506 546" # channel: [-11,-21] + RV_24: "345 375 433 467 507 547" # channel: [-21,-11] + RV_25: "346 376 434 468 508 548" # channel: [-21,-21] + RV_26: "339 369 427 469 509 549" # channel: [-11,11] + RV_27: "340 370 428 470 510 550" # channel: [-11,21] + RV_28: "341 371 429 471 511 551" # channel: [-21,11] + RV_29: "342 372 430 472 512 552" # channel: [-21,21] + RV_30: "571 587 603" # channel: [10,10] + RV_31: "572 588 604" # channel: [20,20] + RV_32: "575 591 607" # channel: [-10,-10] + RV_33: "576 592 608" # channel: [-20,-20] + VV: VV diff --git a/NNLOJET_ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT.yaml b/NNLOJET_ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT.yaml new file mode 100644 index 00000000..6287b702 --- /dev/null +++ b/NNLOJET_ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT.yaml @@ -0,0 +1,315 @@ +# arXiv number: 1512.02192, hepdata entry: https://www.hepdata.net/record/ins1408516 +# Mass distribution (tables: [35, 36, 37, 38, 40]) +# Rapidity distribution (tables: [29, 30, 31, 32, 33, 34]) +# The fine-grained pt distribution: table 27 +# +# NOTE: this includes both the PT-Y and PT-MLL datasets in one single run +# the low-pt region is separated in order not to skew the statistics +# The low-pt .yaml runcard is provided for reference +# The bins at low mass are better run separately as well to avoid needing extreme statistics + +runname: ATLAS_Z0J_8TEV_PT + +process: + proc: ZJ + sqrts: 8000.0 + +pdf: NNPDF40_nnlo_as_01180 + +techcut: 1e-07 + +histograms: + - name: pT_mll_12_20 + observable: ptz + bins: &lomass + - 45.0 + - 55.0 + - 65.0 + - 75.0 + - 85.0 + - 105.0 + - 150.0 + - 200.0 + - 900.0 + extra_selectors: + - observable: mll + min: 12.0 + max: 20.0 + + - name: pT_mll_20_30 + observable: ptz + bins: *lowmass + extra_selectors: + - observable: mll + min: 20.0 + max: 30.0 + + - name: pT_mll_30_46 + observable: ptz + bins: *lowmass + extra_selectors: + - observable: mll + min: 30.0 + max: 46.0 + + - name: pT_mll_46_66 + observable: ptz + bins: &hipt + - 30.0 + - 37.0 + - 45.0 + - 55.0 + - 65.0 + - 75.0 + - 85.0 + - 105.0 + - 150.0 + - 200.0 + - 900.0 + extra_selectors: + - observable: mll + min: 46.0 + max: 66.0 + + - name: pT_mll_116_150 + observable: ptz + bins: *hipt + extra_selectors: + - observable: mll + min: 116.0 + max: 150.0 + + - name: pt_mll_66_116_y_0.0_0.4 + observable: ptz + bins: *hipt + extra_selectors: + - observable: abs_yz + min: 0.0 + max: 0.4 + - ¢ral_mll + observable: mll + min: 66.0 + max: 116.0 + + - name: pt_mll_66_116_y_0.4_0.8 + observable: ptz + bins: *hipt + extra_selectors: + - observable: abs_yz + min: 0.4 + max: 0.8 + - *central_mll + + - name: pt_mll_66_116_y_0.8_1.2 + observable: ptz + bins: *hipt + extra_selectors: + - observable: abs_yz + min: 0.8 + max: 1.2 + - *central_mll + + - name: pt_mll_66_116_y_1.2_1.6 + observable: ptz + bins: *hipt + extra_selectors: + - observable: abs_yz + min: 1.2 + max: 1.6 + - *central_mll + + - name: pt_mll_66_116_y_1.6_2.0 + observable: ptz + bins: *hipt + extra_selectors: + - observable: abs_yz + min: 1.6 + max: 2.0 + - *central_mll + + - name: pt_mll_66_116_y_2.0_2.4 + observable: ptz + bins: *hipt + extra_selectors: + - observable: abs_yz + min: 2.0 + max: 2.4 + - *central_mll + + - name: pt_mll_66_116_y_0.0_2.4 + observable: ptz + extra_selectors: + - observable: abs_yz + min: 0.0 + max: 2.4 + - *central_mll + bins: + - 30.0 + - 33.0 + - 36.0 + - 39.0 + - 42.0 + - 45.0 + - 48.0 + - 51.0 + - 54.0 + - 57.0 + - 61.0 + - 65.0 + - 70.0 + - 75.0 + - 80.0 + - 85.0 + - 95.0 + - 105.0 + - 125.0 + - 150.0 + - 175.0 + - 200.0 + - 250.0 + - 300.0 + - 350.0 + - 400.0 + - 470.0 + - 550.0 + +multi_channel: 0 + +parameters: + MASS[Z]: 91.1876 + MASS[W]: 80.379 + WIDTH[Z]: 2.4952 + WIDTH[W]: 2.085 + SCHEME[alpha]: Gmu + GF: 1.1663787d-5 + +selectors: + - observable: ptz + min: 30.0 + #min: 45.0 + max: 900.0 + - observable: abs_yz + min: 0.0 + max: 2.4 + - observable: mll + min: 12.0 + #min: 46.0 + max: 150.0 + #max: 46.0 + - observable: abs_ylp + max: 2.4 + - observable: abs_ylm + max: 2.4 + - observable: ptlp + min: 20.0 + - observable: ptlm + min: 20.0 + +scales: + mur: etz + muf: etz + +channels: + LO: LO + R: R + V: V + RRa_1: "163 177 191 205 206 245 246 285 301" # channel: [10,0] + RRa_2: "164 178 192 207 208 247 248 286 302" # channel: [20,0] + RRa_3: "165 179 193 217 218 257 258 287 303" # channel: [10,-10] + RRa_4: "166 180 194 219 220 259 260 288 304" # channel: [20,-20] + RRa_5: "167 181 195 221 222 261 262 291 307" # channel: [0,10] + RRa_6: "168 182 196 223 224 263 264 292 308" # channel: [0,20] + RRa_7: "169 183 197 225 226 265 266 293 309" # channel: [-10,10] + RRa_8: "170 184 198 227 228 267 268 294 310" # channel: [-20,20] + RRa_9: "173 187 201 229 230 269 270 295 311" # channel: [0,-10] + RRa_10: "174 188 202 231 232 271 272 296 312" # channel: [0,-20] + RRa_11: "175 189 203 233 234 273 274 297 313" # channel: [-10,0] + RRa_12: "176 190 204 235 236 275 276 298 314" # channel: [-20,0] + RRa_13: "171 172 185 186 199 200" # channel: [0,0] + RRa_14: "209 249" # channel: [11,-11] + RRa_15: "210 250" # channel: [11,-21] + RRa_16: "211 251" # channel: [21,-11] + RRa_17: "212 252" # channel: [21,-21] + RRa_18: "213 253" # channel: [11,11] + RRa_19: "214 254" # channel: [11,21] + RRa_20: "215 255" # channel: [21,11] + RRa_21: "216 256" # channel: [21,21] + RRa_22: "237 277" # channel: [-11,-11] + RRa_23: "238 278" # channel: [-11,-21] + RRa_24: "239 279" # channel: [-21,-11] + RRa_25: "240 280" # channel: [-21,-21] + RRa_26: "241 281" # channel: [-11,11] + RRa_27: "242 282" # channel: [-11,21] + RRa_28: "243 283" # channel: [-21,11] + RRa_29: "244 284" # channel: [-21,21] + RRa_30: "289 305" # channel: [10,10] + RRa_31: "290 306" # channel: [20,20] + RRa_32: "299 315" # channel: [-10,-10] + RRa_33: "300 316" # channel: [-20,-20] + RRb_1: "163 177 191 205 206 245 246 285 301" # channel: [10,0] + RRb_2: "164 178 192 207 208 247 248 286 302" # channel: [20,0] + RRb_3: "165 179 193 217 218 257 258 287 303" # channel: [10,-10] + RRb_4: "166 180 194 219 220 259 260 288 304" # channel: [20,-20] + RRb_5: "167 181 195 221 222 261 262 291 307" # channel: [0,10] + RRb_6: "168 182 196 223 224 263 264 292 308" # channel: [0,20] + RRb_7: "169 183 197 225 226 265 266 293 309" # channel: [-10,10] + RRb_8: "170 184 198 227 228 267 268 294 310" # channel: [-20,20] + RRb_9: "173 187 201 229 230 269 270 295 311" # channel: [0,-10] + RRb_10: "174 188 202 231 232 271 272 296 312" # channel: [0,-20] + RRb_11: "175 189 203 233 234 273 274 297 313" # channel: [-10,0] + RRb_12: "176 190 204 235 236 275 276 298 314" # channel: [-20,0] + RRb_13: "171 172 185 186 199 200" # channel: [0,0] + RRb_14: "209 249" # channel: [11,-11] + RRb_15: "210 250" # channel: [11,-21] + RRb_16: "211 251" # channel: [21,-11] + RRb_17: "212 252" # channel: [21,-21] + RRb_18: "213 253" # channel: [11,11] + RRb_19: "214 254" # channel: [11,21] + RRb_20: "215 255" # channel: [21,11] + RRb_21: "216 256" # channel: [21,21] + RRb_22: "237 277" # channel: [-11,-11] + RRb_23: "238 278" # channel: [-11,-21] + RRb_24: "239 279" # channel: [-21,-11] + RRb_25: "240 280" # channel: [-21,-21] + RRb_26: "241 281" # channel: [-11,11] + RRb_27: "242 282" # channel: [-11,21] + RRb_28: "243 283" # channel: [-21,11] + RRb_29: "244 284" # channel: [-21,21] + RRb_30: "289 305" # channel: [10,10] + RRb_31: "290 306" # channel: [20,20] + RRb_32: "299 315" # channel: [-10,-10] + RRb_33: "300 316" # channel: [-20,-20] + RV_1: "317 347 377 391 405 435 481 482 521 522 561 562 581 597 613" # channel: [10,0] + RV_2: "318 348 378 392 406 436 483 484 523 524 563 564 582 598 614" # channel: [20,0] + RV_3: "319 349 379 393 407 437 457 458 497 498 537 538 569 585 601" # channel: [10,-10] + RV_4: "320 350 380 394 408 438 459 460 499 500 539 540 570 586 602" # channel: [20,-20] + RV_5: "321 351 381 395 409 439 473 474 513 514 553 554 577 593 609" # channel: [0,10] + RV_6: "322 352 382 396 410 440 475 476 515 516 555 556 578 594 610" # channel: [0,20] + RV_7: "323 353 383 397 411 441 461 462 501 502 541 542 573 589 605" # channel: [-10,10] + RV_8: "324 354 384 398 412 442 463 464 503 504 543 544 574 590 606" # channel: [-20,20] + RV_9: "327 357 387 401 415 445 477 478 517 518 557 558 579 595 611" # channel: [0,-10] + RV_10: "328 358 388 402 416 446 479 480 519 520 559 560 580 596 612" # channel: [0,-20] + RV_11: "329 359 389 403 417 447 485 486 525 526 565 566 583 599 615" # channel: [-10,0] + RV_12: "330 360 390 404 418 448 487 488 527 528 567 568 584 600 616" # channel: [-20,0] + RV_13: "325 326 355 356 385 386 399 400 413 414 443 444" # channel: [0,0] + RV_14: "335 365 423 449 489 529" # channel: [11,-11] + RV_15: "336 366 424 450 490 530" # channel: [11,-21] + RV_16: "337 367 425 451 491 531" # channel: [21,-11] + RV_17: "338 368 426 452 492 532" # channel: [21,-21] + RV_18: "331 361 419 453 493 533" # channel: [11,11] + RV_19: "332 362 420 454 494 534" # channel: [11,21] + RV_20: "333 363 421 455 495 535" # channel: [21,11] + RV_21: "334 364 422 456 496 536" # channel: [21,21] + RV_22: "343 373 431 465 505 545" # channel: [-11,-11] + RV_23: "344 374 432 466 506 546" # channel: [-11,-21] + RV_24: "345 375 433 467 507 547" # channel: [-21,-11] + RV_25: "346 376 434 468 508 548" # channel: [-21,-21] + RV_26: "339 369 427 469 509 549" # channel: [-11,11] + RV_27: "340 370 428 470 510 550" # channel: [-11,21] + RV_28: "341 371 429 471 511 551" # channel: [-21,11] + RV_29: "342 372 430 472 512 552" # channel: [-21,21] + RV_30: "571 587 603" # channel: [10,10] + RV_31: "572 588 604" # channel: [20,20] + RV_32: "575 591 607" # channel: [-10,-10] + RV_33: "576 592 608" # channel: [-20,-20] + VV: VV diff --git a/NNLOJET_ATLAS_Z0J_8TEV_PT/metadata.txt b/NNLOJET_ATLAS_Z0J_8TEV_PT/metadata.txt new file mode 100644 index 00000000..e5c4b3a4 --- /dev/null +++ b/NNLOJET_ATLAS_Z0J_8TEV_PT/metadata.txt @@ -0,0 +1 @@ +hepdata="https://www.hepdata.net/record/ins1408516" diff --git a/NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh b/NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh new file mode 100755 index 00000000..4044d9d7 --- /dev/null +++ b/NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh @@ -0,0 +1,146 @@ +#!/usr/bin/bash + +# +# This script finalizes the grids and makes them compatible with +# the previous set of NNPDF grids for the PT-Y distribution. +# For the PT-MLL distribution the grids span a different phase space. +# The old NNPDF grid were called ATLASZPT8TEV* +# +# It uses an external $PINECARD environment variable which should point to the PINECARD folder +# +# For this script to be able to finalize the grid correctly (with all metadata) +# it needs access to the Final/ folder in which the information about the grid and .dat combination is stored +# +# Note: the lowPT region might mess up with the statistics and it is convenient +# to run it separated from the region that is actually used in the PDF fit. +# If the lowPT region has not been run separately set the following flag to false +# +# Note 2: this script only works with pineappl >= 1 + +if [[ -z $PINECARD ]] +then + # If the external PINECARD variable is not set, set it to the cwd + PINECARD="." +fi + +if [[ -z $FAILFAST ]] +then + FAILFAST=true +fi + +SEPARATED_LOWPT=true +pinecard=$(ls ${PINECARD}/*.yaml) +# This script assumes that there is a .run runcard available in the run folder +nnlojet_runcard=$(ls *.run) +results_pdf="NNPDF40_nnlo_as_01180" +tmp_name="merged.lz4" + +merge_and_finalize() { + # usage: merge_and_finalize $high_pt_grid + # + # Takes a high pt grid, autogenerates the naming for the low pt one + # then merges them together and finally puts in all the metadata that it can find + high_pt_grid=${1} + + if [[ ! -f ${high_pt_grid} ]] + then + echo ${high_pt_grid} not found + if [[ ${FAILFAST} == true ]] + then + echo "Fail and exit" + exit -1 + fi + return 0 + fi + + rm -f ${tmp_name} + + + + prefix=${high_pt_grid%.NNLO*} + suffix=${high_pt_grid#*.NNLO} + low_pt_grid=${prefix}-LOWPT.NNLO${suffix} + pineappl merge ${tmp_name} ${low_pt_grid} ${high_pt_grid} + + results=Final/NNLO${suffix%.pineappl*}.dat + final_grid=${prefix}${suffix} + + pineappl write ${tmp_name} --set-key-value nnlojet_runcard "$(cat $nnlojet_runcard)" --set-key-value pinecard "$(cat $pinecard)" --set-key-value results "$(cat $results)" --set-key-value results_pdf ${results_pdf} ${final_grid} + echo Written ${final_grid} + + rm -f ${tmp_name} +} + +### NNLOJET name - NNPDF name +## PT-MLL distribution +# In this case the NNPDF grids already start with a pT cut (for mll > 46 GeV, pt > 30 GeV) +# As a result, if we want to consider bins with lower pT, we need to use a `legacy_theory` implementation anyway. +# Therefore, while they are keep here for reference, the old names are not being used + +# Bins that have changed +bin_1=("ATLAS_Z0J_8TEV_PT.NNLO.pT_mll_46_66.pineappl.lz4" "ATLASZPT8TEVMDIST-ATLASZPT8TEV-MLLBIN4_ptZ.pineappl.lz4") +bin_2=("ATLAS_Z0J_8TEV_PT.NNLO.pT_mll_116_150.pineappl.lz4" "ATLASZPT8TEVMDIST-ATLASZPT8TEV-MLLBIN6_ptZ.pineappl.lz4") +# +for bin in bin_1 bin_2 +do + rm -f ${tmp_name} + declare -n bin_ref=$bin + + high_pt_grid=${bin_ref[0]} + final_name=${high_pt_grid/.NNLO./.} + + merge_and_finalize ${bin_ref[0]} + + mv ${final_name} ${tmp_name} + bins=$(pineappl read ${tmp_name} --bins | awk 'END{print $1}') + pineappl write ${tmp_name} --delete-bins 1-${bins} --scale 0.0 --set-bins 0,2 empty_bin_${tmp_name} + pineappl merge ${final_name} empty_bin_${tmp_name} ${tmp_name} + rm -f ${tmp_name} empty_bin_${tmp_name} +done + +# Bins that are kept the same (but we keep the new names anyway) +bin_1=("ATLAS_Z0J_8TEV_PT.NNLO.pT_mll_12_20.pineappl.lz4" "ATLASZPT8TEVMDIST-ATLASZPT8TEV-MLLBIN1_ptZ.pineappl.lz4") +bin_2=("ATLAS_Z0J_8TEV_PT.NNLO.pT_mll_20_30.pineappl.lz4" "ATLASZPT8TEVMDIST-ATLASZPT8TEV-MLLBIN2_ptZ.pineappl.lz4") +bin_3=("ATLAS_Z0J_8TEV_PT.NNLO.pT_mll_30_46.pineappl.lz4" "ATLASZPT8TEVMDIST-ATLASZPT8TEV-MLLBIN3_ptZ.pineappl.lz4") +for bin in bin_1 bin_2 bin_3 +do + rm -f ${tmp_name} + declare -n high_pt_grid=${bin} + + # For these three there's no need to do any merging, just add the metadata + prefix=${high_pt_grid%.NNLO*} + suffix=${high_pt_grid#*.NNLO} + results=Final/NNLO${suffix%.pineappl*}.dat + final_grid=${prefix}${suffix} + + pineappl write ${high_pt_grid} --set-key-value nnlojet_runcard "$(cat $nnlojet_runcard)" --set-key-value pinecard "$(cat $pinecard)" --set-key-value results "$(cat $results)" --set-key-value results_pdf ${results_pdf} ${final_grid} + + echo Written ${final_grid} +done + + +## PT-Y distribution +# In this case the old names will used: ATLASZPT8TEVYDIST-ATLASZPT8TEVYDIST-BIN${i}_ptZ.pineappl.lz4 +bin_1="ATLAS_Z0J_8TEV_PT.NNLO.pt_mll_66_116_y_0.0_0.4.pineappl.lz4" +bin_2="ATLAS_Z0J_8TEV_PT.NNLO.pt_mll_66_116_y_0.4_0.8.pineappl.lz4" +bin_3="ATLAS_Z0J_8TEV_PT.NNLO.pt_mll_66_116_y_0.8_1.2.pineappl.lz4" +bin_4="ATLAS_Z0J_8TEV_PT.NNLO.pt_mll_66_116_y_1.2_1.6.pineappl.lz4" +bin_5="ATLAS_Z0J_8TEV_PT.NNLO.pt_mll_66_116_y_1.6_2.0.pineappl.lz4" +bin_6="ATLAS_Z0J_8TEV_PT.NNLO.pt_mll_66_116_y_2.0_2.4.pineappl.lz4" + +for i in {1..6} +do + declare -n high_pt_grid="bin_$i" + + merge_and_finalize ${high_pt_grid} + + intermediate_name=${high_pt_grid/.NNLO./.} + final_name=ATLASZPT8TEVYDIST-ATLASZPT8TEVYDIST-BIN${i}_ptZ.pineappl.lz4 + + # Add a 0 pt bin at the beginning + rm -f ${tmp_name} + bins=$(pineappl read ${intermediate_name} --bins | awk 'END{print $1}') + pineappl write ${intermediate_name} --delete-bins 1-${bins} --scale 0.0 --set-bins 0,2 ${tmp_name} + pineappl merge ${final_name} ${tmp_name} ${intermediate_name} + rm -f ${tmp_name} +done diff --git a/NNLOJET_CMS_Z0J_8TEV_PT-Y/CMS_Z0J_8TEV_PT-Y.yaml b/NNLOJET_CMS_Z0J_8TEV_PT-Y/CMS_Z0J_8TEV_PT-Y.yaml new file mode 100644 index 00000000..31054336 --- /dev/null +++ b/NNLOJET_CMS_Z0J_8TEV_PT-Y/CMS_Z0J_8TEV_PT-Y.yaml @@ -0,0 +1,274 @@ +# arXiv number: 1504.03511, hepdata entry: https://doi.org/None (tables: [2]) + +runname: CMS_Z0J_8TEV_PT-Y + +process: + proc: ZJ + sqrts: 8000.0 + jet: antikt[0.5] + +pdf: NNPDF40_nnlo_as_01180 + +techcut: 1e-07 + +histograms: + - name: pT_bin_0 + observable: ptz + bins: + - 20.0 + - 40.0 + - 60.0 + - 80.0 + - 100.0 + - 120.0 + - 140.0 + - 170.0 + - 200.0 + - 1000.0 + extra_selectors: + - observable: abs_yz + min: 0.0 + max: 0.4 + + - name: pT_bin_1 + observable: ptz + bins: + - 20.0 + - 40.0 + - 60.0 + - 80.0 + - 100.0 + - 120.0 + - 140.0 + - 170.0 + - 200.0 + - 1000.0 + extra_selectors: + - observable: abs_yz + min: 0.4 + max: 0.8 + + - name: pT_bin_2 + observable: ptz + bins: + - 20.0 + - 40.0 + - 60.0 + - 80.0 + - 100.0 + - 120.0 + - 140.0 + - 170.0 + - 200.0 + - 1000.0 + extra_selectors: + - observable: abs_yz + min: 0.8 + max: 1.2 + + - name: pT_bin_3 + observable: ptz + bins: + - 20.0 + - 40.0 + - 60.0 + - 80.0 + - 100.0 + - 120.0 + - 140.0 + - 170.0 + - 200.0 + - 1000.0 + extra_selectors: + - observable: abs_yz + min: 1.2 + max: 1.6 + + - name: pT_bin_4 + observable: ptz + bins: + - 20.0 + - 40.0 + - 60.0 + - 80.0 + - 100.0 + - 120.0 + - 140.0 + - 170.0 + - 200.0 + - 1000.0 + extra_selectors: + - observable: abs_yz + min: 1.6 + max: 2.0 + +multi_channel: 0 + + +parameters: + MASS[Z]: 91.1876 + MASS[W]: 80.379 + WIDTH[Z]: 2.4952 + WIDTH[W]: 2.085 + SCHEME[alpha]: Gmu + GF: 1.1663787d-5 + +selectors: + - observable: abs_yz + min: 0.0 + max: 2.0 + - observable: ptz + min: 20.0 + max: 1000.0 + - observable: ptl1 + min: 25.0 + - observable: ptl2 + min: 10.0 + - observable: abs_yl1 + max: 2.1 + - observable: abs_yl2 + max: 2.4 + - observable: mll + min: 81.0 + max: 101.0 + +scales: + mur: etz + muf: etz + +channels: + LO: LO + R: R + V: V + RRa_1: "163 177 191 205 206 245 246 285 301" # channel: [10,0] + RRa_2: "164 178 192 207 208 247 248 286 302" # channel: [20,0] + RRa_3: "165 179 193 217 218 257 258 287 303" # channel: [10,-10] + RRa_4: "166 180 194 219 220 259 260 288 304" # channel: [20,-20] + RRa_5: "167 181 195 221 222 261 262 291 307" # channel: [0,10] + RRa_6: "168 182 196 223 224 263 264 292 308" # channel: [0,20] + RRa_7: "169 183 197 225 226 265 266 293 309" # channel: [-10,10] + RRa_8: "170 184 198 227 228 267 268 294 310" # channel: [-20,20] + RRa_9: "173 187 201 229 230 269 270 295 311" # channel: [0,-10] + RRa_10: "174 188 202 231 232 271 272 296 312" # channel: [0,-20] + RRa_11: "175 189 203 233 234 273 274 297 313" # channel: [-10,0] + RRa_12: "176 190 204 235 236 275 276 298 314" # channel: [-20,0] + RRa_13: "171 172 185 186 199 200" # channel: [0,0] + RRa_14: "209 249" # channel: [11,-11] + RRa_15: "210 250" # channel: [11,-21] + RRa_16: "211 251" # channel: [21,-11] + RRa_17: "212 252" # channel: [21,-21] + RRa_18: "213 253" # channel: [11,11] + RRa_19: "214 254" # channel: [11,21] + RRa_20: "215 255" # channel: [21,11] + RRa_21: "216 256" # channel: [21,21] + RRa_22: "237 277" # channel: [-11,-11] + RRa_23: "238 278" # channel: [-11,-21] + RRa_24: "239 279" # channel: [-21,-11] + RRa_25: "240 280" # channel: [-21,-21] + RRa_26: "241 281" # channel: [-11,11] + RRa_27: "242 282" # channel: [-11,21] + RRa_28: "243 283" # channel: [-21,11] + RRa_29: "244 284" # channel: [-21,21] + RRa_30: "289 305" # channel: [10,10] + RRa_31: "290 306" # channel: [20,20] + RRa_32: "299 315" # channel: [-10,-10] + RRa_33: "300 316" # channel: [-20,-20] + RRb_1: "163 177 191 205 206 245 246 285 301" # channel: [10,0] + RRb_2: "164 178 192 207 208 247 248 286 302" # channel: [20,0] + RRb_3: "165 179 193 217 218 257 258 287 303" # channel: [10,-10] + RRb_4: "166 180 194 219 220 259 260 288 304" # channel: [20,-20] + RRb_5: "167 181 195 221 222 261 262 291 307" # channel: [0,10] + RRb_6: "168 182 196 223 224 263 264 292 308" # channel: [0,20] + RRb_7: "169 183 197 225 226 265 266 293 309" # channel: [-10,10] + RRb_8: "170 184 198 227 228 267 268 294 310" # channel: [-20,20] + RRb_9: "173 187 201 229 230 269 270 295 311" # channel: [0,-10] + RRb_10: "174 188 202 231 232 271 272 296 312" # channel: [0,-20] + RRb_11: "175 189 203 233 234 273 274 297 313" # channel: [-10,0] + RRb_12: "176 190 204 235 236 275 276 298 314" # channel: [-20,0] + RRb_13: "171 172 185 186 199 200" # channel: [0,0] + RRb_14: "209 249" # channel: [11,-11] + RRb_15: "210 250" # channel: [11,-21] + RRb_16: "211 251" # channel: [21,-11] + RRb_17: "212 252" # channel: [21,-21] + RRb_18: "213 253" # channel: [11,11] + RRb_19: "214 254" # channel: [11,21] + RRb_20: "215 255" # channel: [21,11] + RRb_21: "216 256" # channel: [21,21] + RRb_22: "237 277" # channel: [-11,-11] + RRb_23: "238 278" # channel: [-11,-21] + RRb_24: "239 279" # channel: [-21,-11] + RRb_25: "240 280" # channel: [-21,-21] + RRb_26: "241 281" # channel: [-11,11] + RRb_27: "242 282" # channel: [-11,21] + RRb_28: "243 283" # channel: [-21,11] + RRb_29: "244 284" # channel: [-21,21] + RRb_30: "289 305" # channel: [10,10] + RRb_31: "290 306" # channel: [20,20] + RRb_32: "299 315" # channel: [-10,-10] + RRb_33: "300 316" # channel: [-20,-20] + RV_1: "317 347 377 391 405 435 481 482 521 522 561 562 581 597 613" # channel: [10,0] + RV_2: "318 348 378 392 406 436 483 484 523 524 563 564 582 598 614" # channel: [20,0] + RV_3: "319 349 379 393 407 437 457 458 497 498 537 538 569 585 601" # channel: [10,-10] + RV_4: "320 350 380 394 408 438 459 460 499 500 539 540 570 586 602" # channel: [20,-20] + RV_5: "321 351 381 395 409 439 473 474 513 514 553 554 577 593 609" # channel: [0,10] + RV_6: "322 352 382 396 410 440 475 476 515 516 555 556 578 594 610" # channel: [0,20] + RV_7: "323 353 383 397 411 441 461 462 501 502 541 542 573 589 605" # channel: [-10,10] + RV_8: "324 354 384 398 412 442 463 464 503 504 543 544 574 590 606" # channel: [-20,20] + RV_9: "327 357 387 401 415 445 477 478 517 518 557 558 579 595 611" # channel: [0,-10] + RV_10: "328 358 388 402 416 446 479 480 519 520 559 560 580 596 612" # channel: [0,-20] + RV_11: "329 359 389 403 417 447 485 486 525 526 565 566 583 599 615" # channel: [-10,0] + RV_12: "330 360 390 404 418 448 487 488 527 528 567 568 584 600 616" # channel: [-20,0] + RV_13: "325 326 355 356 385 386 399 400 413 414 443 444" # channel: [0,0] + RV_14: "335 365 423 449 489 529" # channel: [11,-11] + RV_15: "336 366 424 450 490 530" # channel: [11,-21] + RV_16: "337 367 425 451 491 531" # channel: [21,-11] + RV_17: "338 368 426 452 492 532" # channel: [21,-21] + RV_18: "331 361 419 453 493 533" # channel: [11,11] + RV_19: "332 362 420 454 494 534" # channel: [11,21] + RV_20: "333 363 421 455 495 535" # channel: [21,11] + RV_21: "334 364 422 456 496 536" # channel: [21,21] + RV_22: "343 373 431 465 505 545" # channel: [-11,-11] + RV_23: "344 374 432 466 506 546" # channel: [-11,-21] + RV_24: "345 375 433 467 507 547" # channel: [-21,-11] + RV_25: "346 376 434 468 508 548" # channel: [-21,-21] + RV_26: "339 369 427 469 509 549" # channel: [-11,11] + RV_27: "340 370 428 470 510 550" # channel: [-11,21] + RV_28: "341 371 429 471 511 551" # channel: [-21,11] + RV_29: "342 372 430 472 512 552" # channel: [-21,21] + RV_30: "571 587 603" # channel: [10,10] + RV_31: "572 588 604" # channel: [20,20] + RV_32: "575 591 607" # channel: [-10,-10] + RV_33: "576 592 608" # channel: [-20,-20] + VV_1: "617 651 685 719 731 733 769 781 783 819 831 833" # channel: [10,0] + VV_2: "618 652 686 720 732 734 770 782 784 820 832 834" # channel: [20,0] + VV_3: "619 653 687 721 771 821" # channel: [10,-10] + VV_4: "620 654 688 722 772 822" # channel: [20,-20] + VV_5: "621 655 689 723 739 741 773 789 791 823 839 841" # channel: [0,10] + VV_6: "622 656 690 724 740 742 774 790 792 824 840 842" # channel: [0,20] + VV_7: "623 657 691 725 775 825" # channel: [-10,10] + VV_8: "624 658 692 726 776 826" # channel: [-20,20] + VV_9: "625 659 693 727 743 745 777 793 795 827 843 845" # channel: [0,-10] + VV_10: "626 660 694 728 744 746 778 794 796 828 844 846" # channel: [0,-20] + VV_11: "627 661 695 729 735 737 779 785 787 829 835 837" # channel: [-10,0] + VV_12: "628 662 696 730 736 738 780 786 788 830 836 838" # channel: [-20,0] + VV_13: "629 630 663 664 697 698 747 748 797 798 847 848" # channel: [0,0] + VV_14: "635 669 703 753 803 853" # channel: [11,-11] + VV_15: "636 670 704 754 804 854" # channel: [11,-21] + VV_16: "637 671 705 755 805 855" # channel: [21,-11] + VV_17: "638 672 706 756 806 856" # channel: [21,-21] + VV_18: "631 665 699 749 799 849" # channel: [11,11] + VV_19: "632 666 700 750 800 850" # channel: [11,21] + VV_20: "633 667 701 751 801 851" # channel: [21,11] + VV_21: "634 668 702 752 802 852" # channel: [21,21] + VV_22: "645 679 713 763 813 863" # channel: [-11,-11] + VV_23: "646 680 714 764 814 864" # channel: [-11,-21] + VV_24: "647 681 715 765 815 865" # channel: [-21,-11] + VV_25: "648 682 716 766 816 866" # channel: [-21,-21] + VV_26: "641 675 709 759 809 859" # channel: [-11,11] + VV_27: "642 676 710 760 810 860" # channel: [-11,21] + VV_28: "643 677 711 761 811 861" # channel: [-21,11] + VV_29: "644 678 712 762 812 862" # channel: [-21,21] + VV_30: "639 673 707 757 807 857" # channel: [10,10] + VV_31: "640 674 708 758 808 858" # channel: [20,20] + VV_32: "649 683 717 767 817 867" # channel: [-10,-10] + VV_33: "650 684 718 768 818 868" # channel: [-20,-20] diff --git a/NNLOJET_CMS_Z0J_8TEV_PT-Y/metadata.txt b/NNLOJET_CMS_Z0J_8TEV_PT-Y/metadata.txt new file mode 100644 index 00000000..fbc55cee --- /dev/null +++ b/NNLOJET_CMS_Z0J_8TEV_PT-Y/metadata.txt @@ -0,0 +1 @@ +hepdata="https://www.hepdata.net/record/ins1359450" diff --git a/NNLOJET_CMS_Z0J_8TEV_PT-Y/postrun.sh b/NNLOJET_CMS_Z0J_8TEV_PT-Y/postrun.sh new file mode 100755 index 00000000..38e3ef00 --- /dev/null +++ b/NNLOJET_CMS_Z0J_8TEV_PT-Y/postrun.sh @@ -0,0 +1,88 @@ +#!/usr/bin/bash + +# +# This script finalizes the grids and makes them compatible with +# the previous set of NNPDF grids +# (which were called CMSZDIFF12-CMSZDIFF12-BINX) +# +# It uses an external $PINECARD environment variable which should point to the PINECARD folder +# +# For this script to be able to finalize the grid correctly (with all metadata) +# it needs access to the folder (remote or local) in which the grids were combined. +# Please modify the remote_path directory below accordingly. +# +# +# Note: the lowPT region might mess up with the statistics and it is convenient +# to run it separated from the region that is actually used in the PDF fit. +# If the lowPT region has not been run separately set the following flag to false +# +# Note 2: this script only works with pineappl >= 1 + +if [[ -z $PINECARD ]] +then + # If the external PINECARD variable is not set, set it to the cwd + PINECARD="." +fi + +if [[ -z $FAILFAST ]] +then + FAILFAST=true +fi + + +remote_path="" +get_results() { + bin=${1} + data_file=combined/Final/NNLO.pT_bin_${i}.dat + + if [[ -z ${remote_path} ]] + then + cp ${data_file} results.dat + else + scp ${remote_path}/${data_file} results.dat + fi +} + +# This script assumes that there is a .run runcard available +nnlojet_runcard=$(ls *.run) + +SEPARATED_LOWPT=true +low_pt_grid=CMS_Z0J_8TEV_PT-Y-LOWPT.NNLO.abs_yz.pineappl.lz4 + +pinecard=$(ls ${PINECARD}/*.yaml) +results_pdf=NNPDF40_nnlo_as_01180 + +results=results.dat + +for i in {0..4} +do + input_grid=CMS_Z0J_8TEV_PT-Y.NNLO.pT_bin_${i}.pineappl.lz4 + + # Get the results as given by NNLOJET + get_results ${1} + + if [ "$SEPARATED_LOWPT" = true ] + then + # Prepare the combination of bins to remove + remove_bins=$(seq 0 4 | tr '\n' ',' | sed "s/$i,//" | sed 's/,$//') + + # Use as output name the name from the old NNPDF grids + # output_name=CMS_Z0J_8TEV_PT-Y-LOWPT.NNLO.pT_bin_${i}.pineappl.lz4 + output_name=CMSZDIFF12-CMSZDIFF12-BIN$(( i+1 ))_ptZ.pineappl.lz4 + + # Make the bins compatible + rm -f tmp.lz4 + pineappl write $low_pt_grid --delete-bins "${remove_bins}" --set-bins "2,20" --set-key-value x1_label "pT_bin_${i}" tmp.lz4 + + # Merge! + rm -f merged_tmp.lz4 + pineappl merge merged_tmp.lz4 tmp.lz4 $input_grid + rm -f tmp.lz4 + + # And now the input grid is merged_tmp.lz4 + input_grid=merged_tmp.lz4 + fi + + # And, finally, set up the metadata + pineappl write ${input_grid} --set-key-value nnlojet_runcard "$(cat $nnlojet_runcard)" --set-key-value pinecard "$(cat $pinecard)" --set-key-value results "$(cat $results)" --set-key-value results_pdf ${results_pdf} ${output_name} +done