From c584bf0d730b17a5eed19088cc027d069900a0f3 Mon Sep 17 00:00:00 2001 From: juacrumar Date: Thu, 20 Mar 2025 13:29:46 +0100 Subject: [PATCH 1/5] first batch of ZpT runcards --- .../ATLAS_Z0J_13TEV_PT_LL.yaml | 217 +++++++++++++ ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT.yaml | 307 ++++++++++++++++++ CMS_Z0J_8TEV_PT-Y/CMS_Z0J_8TEV_PT-Y.yaml | 274 ++++++++++++++++ CMS_Z0J_8TEV_PT-Y/finalize.sh | 77 +++++ CMS_Z0J_8TEV_PT-Y/metadata.txt | 1 + 5 files changed, 876 insertions(+) create mode 100644 ATLAS_Z0J_13TEV_PT_LL/ATLAS_Z0J_13TEV_PT_LL.yaml create mode 100644 ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT.yaml create mode 100644 CMS_Z0J_8TEV_PT-Y/CMS_Z0J_8TEV_PT-Y.yaml create mode 100644 CMS_Z0J_8TEV_PT-Y/finalize.sh create mode 100644 CMS_Z0J_8TEV_PT-Y/metadata.txt diff --git a/ATLAS_Z0J_13TEV_PT_LL/ATLAS_Z0J_13TEV_PT_LL.yaml b/ATLAS_Z0J_13TEV_PT_LL/ATLAS_Z0J_13TEV_PT_LL.yaml new file mode 100644 index 00000000..178caed6 --- /dev/null +++ b/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/ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT.yaml b/ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT.yaml new file mode 100644 index 00000000..373f5693 --- /dev/null +++ b/ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT.yaml @@ -0,0 +1,307 @@ +# 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: &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: 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: 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: 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: 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: 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: yz + min: 2.0 + max: 2.4 + - *central_mll + + - name: pt_mll_66_116_y_0.0_2.4 + observable: ptz + extra_selectors: + - observable: 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 + max: 900.0 + - observable: 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 +# 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/CMS_Z0J_8TEV_PT-Y/CMS_Z0J_8TEV_PT-Y.yaml b/CMS_Z0J_8TEV_PT-Y/CMS_Z0J_8TEV_PT-Y.yaml new file mode 100644 index 00000000..31054336 --- /dev/null +++ b/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/CMS_Z0J_8TEV_PT-Y/finalize.sh b/CMS_Z0J_8TEV_PT-Y/finalize.sh new file mode 100644 index 00000000..f0f8ef6f --- /dev/null +++ b/CMS_Z0J_8TEV_PT-Y/finalize.sh @@ -0,0 +1,77 @@ +#!/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) +# +# 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 + + +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 +} + +# In principle, if there was no local modification we can autogenerate the LO runcard +# this is enough for NNLOJET dokan to run the calculation again +nnlojet_runcard=LO/runcard_production.run + +SEPARATED_LOWPT=true +low_pt_grid=CMS_Z0J_8TEV_PT-Y-LOWPT.NNLO.abs_yz.pineappl.lz4 + +pinecard=$(ls *.yaml) +results_pdf=NNPDF40_nnlo_as_01180 + +nnlojet_runcard=runcard.run +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 diff --git a/CMS_Z0J_8TEV_PT-Y/metadata.txt b/CMS_Z0J_8TEV_PT-Y/metadata.txt new file mode 100644 index 00000000..fbc55cee --- /dev/null +++ b/CMS_Z0J_8TEV_PT-Y/metadata.txt @@ -0,0 +1 @@ +hepdata="https://www.hepdata.net/record/ins1359450" From f04aac7324451618f54836f6fda4a19bf7484430 Mon Sep 17 00:00:00 2001 From: juacrumar Date: Mon, 31 Mar 2025 11:58:24 +0200 Subject: [PATCH 2/5] finalize pinecards for CMS and ATLAS Z0J 8 TeV with postrun.sh scripts --- .../ATLAS_Z0J_13TEV_PT_LL.yaml | 0 .../ATLAS_Z0J_8TEV_PT-lowpt.yaml | 245 ++++++++++++++++++ .../ATLAS_Z0J_8TEV_PT.yaml | 53 +--- NNLOJET_ATLAS_Z0J_8TEV_PT/metadata.txt | 1 + NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh | 94 +++++++ .../CMS_Z0J_8TEV_PT-Y.yaml | 0 .../metadata.txt | 0 .../postrun.sh | 21 +- 8 files changed, 368 insertions(+), 46 deletions(-) rename {ATLAS_Z0J_13TEV_PT_LL => NNLOJET_ATLAS_Z0J_13TEV_PT_LL}/ATLAS_Z0J_13TEV_PT_LL.yaml (100%) create mode 100644 NNLOJET_ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT-lowpt.yaml rename {ATLAS_Z0J_8TEV_PT => NNLOJET_ATLAS_Z0J_8TEV_PT}/ATLAS_Z0J_8TEV_PT.yaml (78%) create mode 100644 NNLOJET_ATLAS_Z0J_8TEV_PT/metadata.txt create mode 100755 NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh rename {CMS_Z0J_8TEV_PT-Y => NNLOJET_CMS_Z0J_8TEV_PT-Y}/CMS_Z0J_8TEV_PT-Y.yaml (100%) rename {CMS_Z0J_8TEV_PT-Y => NNLOJET_CMS_Z0J_8TEV_PT-Y}/metadata.txt (100%) rename CMS_Z0J_8TEV_PT-Y/finalize.sh => NNLOJET_CMS_Z0J_8TEV_PT-Y/postrun.sh (85%) mode change 100644 => 100755 diff --git a/ATLAS_Z0J_13TEV_PT_LL/ATLAS_Z0J_13TEV_PT_LL.yaml b/NNLOJET_ATLAS_Z0J_13TEV_PT_LL/ATLAS_Z0J_13TEV_PT_LL.yaml similarity index 100% rename from ATLAS_Z0J_13TEV_PT_LL/ATLAS_Z0J_13TEV_PT_LL.yaml rename to NNLOJET_ATLAS_Z0J_13TEV_PT_LL/ATLAS_Z0J_13TEV_PT_LL.yaml 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/ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT.yaml b/NNLOJET_ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT.yaml similarity index 78% rename from ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT.yaml rename to NNLOJET_ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT.yaml index 373f5693..ce756396 100644 --- a/ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT.yaml +++ b/NNLOJET_ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT.yaml @@ -2,6 +2,10 @@ # 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 runname: ATLAS_Z0J_8TEV_PT @@ -45,7 +49,7 @@ histograms: observable: ptz bins: *hipt extra_selectors: - - observable: yz + - observable: abs_yz min: 0.0 max: 0.4 - ¢ral_mll @@ -57,7 +61,7 @@ histograms: observable: ptz bins: *hipt extra_selectors: - - observable: yz + - observable: abs_yz min: 0.4 max: 0.8 - *central_mll @@ -66,7 +70,7 @@ histograms: observable: ptz bins: *hipt extra_selectors: - - observable: yz + - observable: abs_yz min: 0.8 max: 1.2 - *central_mll @@ -75,7 +79,7 @@ histograms: observable: ptz bins: *hipt extra_selectors: - - observable: yz + - observable: abs_yz min: 1.2 max: 1.6 - *central_mll @@ -84,7 +88,7 @@ histograms: observable: ptz bins: *hipt extra_selectors: - - observable: yz + - observable: abs_yz min: 1.6 max: 2.0 - *central_mll @@ -93,7 +97,7 @@ histograms: observable: ptz bins: *hipt extra_selectors: - - observable: yz + - observable: abs_yz min: 2.0 max: 2.4 - *central_mll @@ -101,7 +105,7 @@ histograms: - name: pt_mll_66_116_y_0.0_2.4 observable: ptz extra_selectors: - - observable: yz + - observable: abs_yz min: 0.0 max: 2.4 - *central_mll @@ -149,7 +153,7 @@ selectors: - observable: ptz min: 30.0 max: 900.0 - - observable: yz + - observable: abs_yz min: 0.0 max: 2.4 - observable: mll @@ -272,36 +276,3 @@ channels: RV_32: "575 591 607" # channel: [-10,-10] RV_33: "576 592 608" # channel: [-20,-20] VV: VV -# 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_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..851fb247 --- /dev/null +++ b/NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh @@ -0,0 +1,94 @@ +#!/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" + +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} + tmp_name="merged.lz4" + 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} + + 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 +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 +# + merge_and_finalize ${bin_ref[0]} +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 + mv ${intermediate_name} ${final_name} +done diff --git a/CMS_Z0J_8TEV_PT-Y/CMS_Z0J_8TEV_PT-Y.yaml b/NNLOJET_CMS_Z0J_8TEV_PT-Y/CMS_Z0J_8TEV_PT-Y.yaml similarity index 100% rename from CMS_Z0J_8TEV_PT-Y/CMS_Z0J_8TEV_PT-Y.yaml rename to NNLOJET_CMS_Z0J_8TEV_PT-Y/CMS_Z0J_8TEV_PT-Y.yaml diff --git a/CMS_Z0J_8TEV_PT-Y/metadata.txt b/NNLOJET_CMS_Z0J_8TEV_PT-Y/metadata.txt similarity index 100% rename from CMS_Z0J_8TEV_PT-Y/metadata.txt rename to NNLOJET_CMS_Z0J_8TEV_PT-Y/metadata.txt diff --git a/CMS_Z0J_8TEV_PT-Y/finalize.sh b/NNLOJET_CMS_Z0J_8TEV_PT-Y/postrun.sh old mode 100644 new mode 100755 similarity index 85% rename from CMS_Z0J_8TEV_PT-Y/finalize.sh rename to NNLOJET_CMS_Z0J_8TEV_PT-Y/postrun.sh index f0f8ef6f..38e3ef00 --- a/CMS_Z0J_8TEV_PT-Y/finalize.sh +++ b/NNLOJET_CMS_Z0J_8TEV_PT-Y/postrun.sh @@ -5,6 +5,8 @@ # 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. @@ -16,6 +18,17 @@ # # 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() { @@ -30,17 +43,15 @@ get_results() { fi } -# In principle, if there was no local modification we can autogenerate the LO runcard -# this is enough for NNLOJET dokan to run the calculation again -nnlojet_runcard=LO/runcard_production.run +# 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 *.yaml) +pinecard=$(ls ${PINECARD}/*.yaml) results_pdf=NNPDF40_nnlo_as_01180 -nnlojet_runcard=runcard.run results=results.dat for i in {0..4} From a7e71531f437c71b00f20680d18697dbf1fbb1ec Mon Sep 17 00:00:00 2001 From: juacrumar Date: Tue, 1 Apr 2025 18:12:57 +0200 Subject: [PATCH 3/5] add a bin0 for easier comparison with older grids --- NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) diff --git a/NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh b/NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh index 851fb247..41f36eae 100755 --- a/NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh +++ b/NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh @@ -90,5 +90,11 @@ do intermediate_name=${high_pt_grid/.NNLO./.} final_name=ATLASZPT8TEVYDIST-ATLASZPT8TEVYDIST-BIN${i}_ptZ.pineappl.lz4 - mv ${intermediate_name} ${final_name} + + # 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 From e05b0ac8fe47bf35120dba34b399fd0042eb7856 Mon Sep 17 00:00:00 2001 From: juacrumar Date: Mon, 28 Apr 2025 12:06:48 +0200 Subject: [PATCH 4/5] update pinecard to add the low mass bins --- .../ATLAS_Z0J_8TEV_PT.yaml | 39 ++++++++++++++++++- NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh | 39 ++++++++++++++++++- 2 files changed, 76 insertions(+), 2 deletions(-) diff --git a/NNLOJET_ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT.yaml b/NNLOJET_ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT.yaml index ce756396..6287b702 100644 --- a/NNLOJET_ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT.yaml +++ b/NNLOJET_ATLAS_Z0J_8TEV_PT/ATLAS_Z0J_8TEV_PT.yaml @@ -6,6 +6,7 @@ # 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 @@ -18,6 +19,39 @@ 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 @@ -152,13 +186,16 @@ parameters: selectors: - observable: ptz min: 30.0 + #min: 45.0 max: 900.0 - observable: abs_yz min: 0.0 max: 2.4 - observable: mll - min: 46.0 + min: 12.0 + #min: 46.0 max: 150.0 + #max: 46.0 - observable: abs_ylp max: 2.4 - observable: abs_ylm diff --git a/NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh b/NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh index 41f36eae..e664d39e 100755 --- a/NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh +++ b/NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh @@ -33,6 +33,7 @@ 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 @@ -40,10 +41,22 @@ merge_and_finalize() { # 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} - tmp_name="merged.lz4" + + 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} @@ -53,6 +66,7 @@ merge_and_finalize() { 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} } @@ -62,6 +76,8 @@ merge_and_finalize() { # 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") # @@ -73,6 +89,27 @@ do merge_and_finalize ${bin_ref[0]} 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" From d29598840886500e8a84d55fa3a6b13c39b3fac5 Mon Sep 17 00:00:00 2001 From: juacrumar Date: Tue, 13 May 2025 10:49:13 +0200 Subject: [PATCH 5/5] update PT-M script --- NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh | 11 ++++++++++- 1 file changed, 10 insertions(+), 1 deletion(-) diff --git a/NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh b/NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh index e664d39e..4044d9d7 100755 --- a/NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh +++ b/NNLOJET_ATLAS_Z0J_8TEV_PT/postrun.sh @@ -85,8 +85,17 @@ 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)