Fork me on GitHub

Ignore:
Timestamp:
Jan 24, 2018, 12:29:43 PM (7 years ago)
Author:
Ulrike Schnoor <schnooru@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
1f4c2c3, 270c0ed
Parents:
3037765
Message:

added BTagging of VLC jets for the 25 combinations of R and N

File:
1 edited

Legend:

Unmodified
Added
Removed

  • cards/delphes_card_CLICdet.tcl

    r3037765 rcd60feb  
    8080    GenMissingET
    8181
    82     JetEnergyScale
    83 
    84     JetFlavorAssociation
    85     BTagging
    86    
    87     TauTagging
     82   
     83
     84    JetFlavorAssociation_R05N2
     85        JetFlavorAssociation_R05N3
     86        JetFlavorAssociation_R05N4
     87        JetFlavorAssociation_R05N5
     88        JetFlavorAssociation_R05N6
     89       
     90        JetFlavorAssociation_R07N2
     91        JetFlavorAssociation_R07N3
     92        JetFlavorAssociation_R07N4
     93        JetFlavorAssociation_R07N5
     94        JetFlavorAssociation_R07N6
     95
     96        JetFlavorAssociation_R10N2
     97        JetFlavorAssociation_R10N3
     98        JetFlavorAssociation_R10N4
     99        JetFlavorAssociation_R10N5
     100        JetFlavorAssociation_R10N6
     101
     102        JetFlavorAssociation_R12N2
     103        JetFlavorAssociation_R12N3
     104        JetFlavorAssociation_R12N4
     105        JetFlavorAssociation_R12N5
     106        JetFlavorAssociation_R12N6
     107
     108        JetFlavorAssociation_R15N2
     109        JetFlavorAssociation_R15N3
     110        JetFlavorAssociation_R15N4
     111        JetFlavorAssociation_R15N5
     112        JetFlavorAssociation_R15N6
     113
     114        BTaggingWP50_R05N2
     115        BTaggingWP70_R05N2
     116        BTaggingWP90_R05N2
     117        BTaggingWP50_R05N3
     118        BTaggingWP70_R05N3
     119        BTaggingWP90_R05N3
     120        BTaggingWP50_R05N4
     121        BTaggingWP70_R05N4
     122        BTaggingWP90_R05N4
     123        BTaggingWP50_R05N5
     124        BTaggingWP70_R05N5
     125        BTaggingWP90_R05N5
     126        BTaggingWP50_R05N6
     127        BTaggingWP70_R05N6
     128        BTaggingWP90_R05N6
     129        BTaggingWP50_R07N2
     130        BTaggingWP70_R07N2
     131        BTaggingWP90_R07N2
     132        BTaggingWP50_R07N3
     133        BTaggingWP70_R07N3
     134        BTaggingWP90_R07N3
     135        BTaggingWP50_R07N4
     136        BTaggingWP70_R07N4
     137        BTaggingWP90_R07N4
     138        BTaggingWP50_R07N5
     139        BTaggingWP70_R07N5
     140        BTaggingWP90_R07N5
     141        BTaggingWP50_R07N6
     142        BTaggingWP70_R07N6
     143        BTaggingWP90_R07N6
     144        BTaggingWP50_R10N2
     145        BTaggingWP70_R10N2
     146        BTaggingWP90_R10N2
     147        BTaggingWP50_R10N3
     148        BTaggingWP70_R10N3
     149        BTaggingWP90_R10N3
     150        BTaggingWP50_R10N4
     151        BTaggingWP70_R10N4
     152        BTaggingWP90_R10N4
     153        BTaggingWP50_R10N5
     154        BTaggingWP70_R10N5
     155        BTaggingWP90_R10N5
     156        BTaggingWP50_R10N6
     157        BTaggingWP70_R10N6
     158        BTaggingWP90_R10N6
     159        BTaggingWP50_R12N2
     160        BTaggingWP70_R12N2
     161        BTaggingWP90_R12N2
     162        BTaggingWP50_R12N3
     163        BTaggingWP70_R12N3
     164        BTaggingWP90_R12N3
     165        BTaggingWP50_R12N4
     166        BTaggingWP70_R12N4
     167        BTaggingWP90_R12N4
     168        BTaggingWP50_R12N5
     169        BTaggingWP70_R12N5
     170        BTaggingWP90_R12N5
     171        BTaggingWP50_R12N6
     172        BTaggingWP70_R12N6
     173        BTaggingWP90_R12N6
     174        BTaggingWP50_R15N2
     175        BTaggingWP70_R15N2
     176        BTaggingWP90_R15N2
     177        BTaggingWP50_R15N3
     178        BTaggingWP70_R15N3
     179        BTaggingWP90_R15N3
     180        BTaggingWP50_R15N4
     181        BTaggingWP70_R15N4
     182        BTaggingWP90_R15N4
     183        BTaggingWP50_R15N5
     184        BTaggingWP70_R15N5
     185        BTaggingWP90_R15N5
     186        BTaggingWP50_R15N6
     187        BTaggingWP70_R15N6
     188        BTaggingWP90_R15N6
     189
     190   
     191    TauTagging_R05N2
    88192
    89193    ScalarHT
    90 
    91194
    92195
    … …  
    712815}
    713816
    714 ############
    715 # Jet finder VLC
    716 ############
    717 #R05 N2
    718 module FastJetFinder FastJetFinderVLC_R05_N2 {
    719     #  set InputArray Calorimeter/towers
    720     set InputArray EFlowFilter/eflow
    721 
    722     set OutputArray VLCjetsR05N2
    723 
    724     # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt, 7 anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness, 9 Valencia
    725     set NJets 2
    726     set ExclusiveClustering true
    727     set JetAlgorithm 9
    728     set ParameterR 0.5
    729     set Beta 1.0
    730     set Gamma 1.0
    731 
    732     set JetPTMin 20.0
    733 }
    734 #R05 N3
    735 module FastJetFinder FastJetFinderVLC_R05_N3 {
    736     set InputArray EFlowFilter/eflow
    737     set OutputArray VLCjetsR05N3
    738     set NJets 3
    739     set ExclusiveClustering true
    740     set JetAlgorithm 9
    741     set ParameterR 0.5
    742     set Beta 1.0
    743     set Gamma 1.0
    744     set JetPTMin 20.0
    745 }
    746 #R05 N4
    747 module FastJetFinder FastJetFinderVLC_R05_N4 {
    748     set InputArray EFlowFilter/eflow
    749     set OutputArray VLCjetsR05N4
    750     set NJets 4
    751     set ExclusiveClustering true
    752     set JetAlgorithm 9
    753     set ParameterR 0.5
    754     set Beta 1.0
    755     set Gamma 1.0
    756     set JetPTMin 20.0
    757 }
    758 #R05 N5
    759 module FastJetFinder FastJetFinderVLC_R05_N5 {
    760     set InputArray EFlowFilter/eflow
    761     set OutputArray VLCjetsR05N5
    762     set NJets 5
    763     set ExclusiveClustering true
    764     set JetAlgorithm 9
    765     set ParameterR 0.5
    766     set Beta 1.0
    767     set Gamma 1.0
    768     set JetPTMin 20.0
    769 }
    770 #R05 N6
    771 module FastJetFinder FastJetFinderVLC_R05_N6 {
    772     set InputArray EFlowFilter/eflow
    773     set OutputArray VLCjetsR05N6
    774     set NJets 6
    775     set ExclusiveClustering true
    776     set JetAlgorithm 9
    777     set ParameterR 0.5
    778     set Beta 1.0
    779     set Gamma 1.0
    780     set JetPTMin 20.0
    781 }
    782 #R07 N2
    783 module FastJetFinder FastJetFinderVLC_R07_N2 {
    784     set InputArray EFlowFilter/eflow
    785     set OutputArray VLCjetsR07N2
    786     set NJets 2
    787     set ExclusiveClustering true
    788     set JetAlgorithm 9
    789     set ParameterR 0.7
    790     set Beta 1.0
    791     set Gamma 1.0
    792     set JetPTMin 20.0
    793 }
    794 #R07 N3
    795 module FastJetFinder FastJetFinderVLC_R07_N3 {
    796     set InputArray EFlowFilter/eflow
    797     set OutputArray VLCjetsR07N3
    798     set NJets 3
    799     set ExclusiveClustering true
    800     set JetAlgorithm 9
    801     set ParameterR 0.7
    802     set Beta 1.0
    803     set Gamma 1.0
    804     set JetPTMin 20.0
    805 }
    806 #R07 N4
    807 module FastJetFinder FastJetFinderVLC_R07_N4 {
    808     set InputArray EFlowFilter/eflow
    809     set OutputArray VLCjetsR07N4
    810     set NJets 4
    811     set ExclusiveClustering true
    812     set JetAlgorithm 9
    813     set ParameterR 0.7
    814     set Beta 1.0
    815     set Gamma 1.0
    816     set JetPTMin 20.0
    817 }
    818 #R07 N5
    819 module FastJetFinder FastJetFinderVLC_R07_N5 {
    820     set InputArray EFlowFilter/eflow
    821     set OutputArray VLCjetsR07N5
    822     set NJets 5
    823     set ExclusiveClustering true
    824     set JetAlgorithm 9
    825     set ParameterR 0.7
    826     set Beta 1.0
    827     set Gamma 1.0
    828     set JetPTMin 20.0
    829 }
    830 #R07 N6
    831 module FastJetFinder FastJetFinderVLC_R07_N6 {
    832     set InputArray EFlowFilter/eflow
    833     set OutputArray VLCjetsR07N6
    834     set NJets 6
    835     set ExclusiveClustering true
    836     set JetAlgorithm 9
    837     set ParameterR 0.7
    838     set Beta 1.0
    839     set Gamma 1.0
    840     set JetPTMin 20.0
    841 }
    842 
    843 #R10N2
    844 module FastJetFinder FastJetFinderVLC_R10_N2 {
    845     set InputArray EFlowFilter/eflow
    846     set OutputArray VLCjetsR10N2
    847     set NJets 2
    848     set ExclusiveClustering true
    849     set JetAlgorithm 9
    850     set ParameterR 1.0
    851     set Beta 1.0
    852     set Gamma 1.0
    853     set JetPTMin 20.0
    854 }
    855 #R10 N3
    856 module FastJetFinder FastJetFinderVLC_R10_N3 {
    857     set InputArray EFlowFilter/eflow
    858     set OutputArray VLCjetsR10N3
    859     set NJets 3
    860     set ExclusiveClustering true
    861     set JetAlgorithm 9
    862     set ParameterR 1.0
    863     set Beta 1.0
    864     set Gamma 1.0
    865     set JetPTMin 20.0
    866 }
    867 #R10 N4
    868 module FastJetFinder FastJetFinderVLC_R10_N4 {
    869     set InputArray EFlowFilter/eflow
    870     set OutputArray VLCjetsR10N4
    871     set NJets 4
    872     set ExclusiveClustering true
    873     set JetAlgorithm 9
    874     set ParameterR 1.0
    875     set Beta 1.0
    876     set Gamma 1.0
    877     set JetPTMin 20.0
    878 }
    879 #R10 N5
    880 module FastJetFinder FastJetFinderVLC_R10_N5 {
    881     set InputArray EFlowFilter/eflow
    882     set OutputArray VLCjetsR10N5
    883     set NJets 5
    884     set ExclusiveClustering true
    885     set JetAlgorithm 9
    886     set ParameterR 1.0
    887     set Beta 1.0
    888     set Gamma 1.0
    889     set JetPTMin 20.0
    890 }
    891 #R10 N6
    892 module FastJetFinder FastJetFinderVLC_R10_N6 {
    893     set InputArray EFlowFilter/eflow
    894     set OutputArray VLCjetsR10N6
    895     set NJets 6
    896     set ExclusiveClustering true
    897     set JetAlgorithm 9
    898     set ParameterR 1.0
    899     set Beta 1.0
    900     set Gamma 1.0
    901     set JetPTMin 20.0
    902 }
    903 
    904 #R12 N2
    905 module FastJetFinder FastJetFinderVLC_R12_N2 {
    906     set InputArray EFlowFilter/eflow
    907     set OutputArray VLCjetsR12N2
    908     set NJets 2
    909     set ExclusiveClustering true
    910     set JetAlgorithm 9
    911     set ParameterR 1.2
    912     set Beta 1.0
    913     set Gamma 1.0
    914     set JetPTMin 20.0
    915 }
    916 #R12 N3
    917 module FastJetFinder FastJetFinderVLC_R12_N3 {
    918     set InputArray EFlowFilter/eflow
    919     set OutputArray VLCjetsR12N3
    920     set NJets 3
    921     set ExclusiveClustering true
    922     set JetAlgorithm 9
    923     set ParameterR 1.2
    924     set Beta 1.0
    925     set Gamma 1.0
    926     set JetPTMin 20.0
    927 }
    928 #R12 N4
    929 module FastJetFinder FastJetFinderVLC_R12_N4 {
    930     set InputArray EFlowFilter/eflow
    931     set OutputArray VLCjetsR12N4
    932     set NJets 4
    933     set ExclusiveClustering true
    934     set JetAlgorithm 9
    935     set ParameterR 1.2
    936     set Beta 1.0
    937     set Gamma 1.0
    938     set JetPTMin 20.0
    939 }
    940 #R12 N5
    941 module FastJetFinder FastJetFinderVLC_R12_N5 {
    942     set InputArray EFlowFilter/eflow
    943     set OutputArray VLCjetsR12N5
    944     set NJets 5
    945     set ExclusiveClustering true
    946     set JetAlgorithm 9
    947     set ParameterR 1.2
    948     set Beta 1.0
    949     set Gamma 1.0
    950     set JetPTMin 20.0
    951 }
    952 #R12 N6
    953 module FastJetFinder FastJetFinderVLC_R12_N6 {
    954     set InputArray EFlowFilter/eflow
    955     set OutputArray VLCjetsR12N6
    956     set NJets 6
    957     set ExclusiveClustering true
    958     set JetAlgorithm 9
    959     set ParameterR 1.2
    960     set Beta 1.0
    961     set Gamma 1.0
    962     set JetPTMin 20.0
    963 }
    964 
    965 
    966 #R15 N2
    967 module FastJetFinder FastJetFinderVLC_R15_N2 {
    968     set InputArray EFlowFilter/eflow
    969     set OutputArray VLCjetsR15N2
    970     set NJets 2
    971     set ExclusiveClustering true
    972     set JetAlgorithm 9
    973     set ParameterR 1.5
    974     set Beta 1.0
    975     set Gamma 1.0
    976     set JetPTMin 20.0
    977 }
    978 #R15 N3
    979 module FastJetFinder FastJetFinderVLC_R15_N3 {
    980     set InputArray EFlowFilter/eflow
    981     set OutputArray VLCjetsR15N3
    982     set NJets 3
    983     set ExclusiveClustering true
    984     set JetAlgorithm 9
    985     set ParameterR 1.5
    986     set Beta 1.0
    987     set Gamma 1.0
    988     set JetPTMin 20.0
    989 }
    990 #R15 N4
    991 module FastJetFinder FastJetFinderVLC_R15_N4 {
    992     set InputArray EFlowFilter/eflow
    993     set OutputArray VLCjetsR15N4
    994     set NJets 4
    995     set ExclusiveClustering true
    996     set JetAlgorithm 9
    997     set ParameterR 1.5
    998     set Beta 1.0
    999     set Gamma 1.0
    1000     set JetPTMin 20.0
    1001 }
    1002 #R15 N5
    1003 module FastJetFinder FastJetFinderVLC_R15_N5 {
    1004     set InputArray EFlowFilter/eflow
    1005     set OutputArray VLCjetsR15N5
    1006     set NJets 5
    1007     set ExclusiveClustering true
    1008     set JetAlgorithm 9
    1009     set ParameterR 1.5
    1010     set Beta 1.0
    1011     set Gamma 1.0
    1012     set JetPTMin 20.0
    1013 }
    1014 #R15 N6
    1015 module FastJetFinder FastJetFinderVLC_R15_N6 {
    1016     set InputArray EFlowFilter/eflow
    1017     set OutputArray VLCjetsR15N6
    1018     set NJets 6
    1019     set ExclusiveClustering true
    1020     set JetAlgorithm 9
    1021     set ParameterR 1.5
    1022     set Beta 1.0
    1023     set Gamma 1.0
    1024     set JetPTMin 20.0
    1025 }
    1026 ##################
    1027 # Jet Energy Scale
    1028 ##################
    1029 
    1030 module EnergyScale JetEnergyScale {
    1031     set InputArray FastJetFinderVLC_R10_N4/VLCjetsR10N4
    1032     set OutputArray jets
    1033 
    1034     # scale formula for jets
    1035     set ScaleFormula {1.00}
    1036 }
     817source CLICdet_JetReco.tcl
     818
     819
     820###################
     821## Jet Energy Scale
     822###################
     823#
     824#module EnergyScale JetEnergyScale {
     825#    set InputArray FastJetFinderVLC_R10_N4/VLCjetsR10N4
     826#    set OutputArray jets
     827#
     828#    # scale formula for jets
     829#    set ScaleFormula {1.00}
     830#}
    1037831
    1038832
    … …  
    1041835########################
    1042836
    1043 module JetFlavorAssociation JetFlavorAssociation {
    1044 
    1045     set PartonInputArray Delphes/partons
    1046     set ParticleInputArray Delphes/allParticles
    1047     set ParticleLHEFInputArray Delphes/allParticlesLHEF
    1048     set JetInputArray JetEnergyScale/jets
    1049 
    1050     set DeltaR 0.5
    1051     set PartonPTMin 1.0
    1052     set PartonEtaMax 2.5
    1053 
    1054 }
    1055 
     837source CLICdet_JetFlavorAssociation.tcl
    1056838
    1057839###########
    1058840# b-tagging
    1059841###########
    1060 
    1061 module BTagging BTagging {
    1062     set JetInputArray JetEnergyScale/jets
    1063 
    1064     set BitNumber 0
    1065 
    1066     # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
    1067     # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
    1068     # gluon's PDG code has the lowest priority
    1069 
    1070     # based on arXiv:1211.4462 (CMS)
    1071    
    1072     # default efficiency formula (misidentification rate)
    1073     add EfficiencyFormula {0} {0.01+0.000038*pt}
    1074 
    1075     # efficiency formula for c-jets (misidentification rate)
    1076     add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
    1077 
    1078     # efficiency formula for b-jets
    1079     add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
    1080 }
     842# based on CLICdp-Note-2014-002   
     843
     844source CLICdet_BTagging.tcl
     845
    1081846
    1082847#############
    … …  
    1085850
    1086851
    1087 module TauTagging TauTagging {
     852module TauTagging TauTagging_R05N2 {
    1088853    set ParticleInputArray Delphes/allParticles
    1089854    set PartonInputArray Delphes/partons
    1090     set JetInputArray JetEnergyScale/jets
     855    set JetInputArray FastJetFinderVLC_R05_N2/VLCjetsR05N2
    1091856
    1092857    set DeltaR 0.5
    … …  
    1160925    add Branch FastJetFinderVLC_R15_N6/VLCjetsR15N6 VLCjetR15N6 Jet
    1161926
    1162    
    1163     add Branch GenMissingET/momentum GenMissingET MissingET
     927#       add Branch BTaggingWP50/jets BTagged Jet
     928#       add Branch FastJetFinderVLC_R05_N2/VLCjetsR05N2 BTagged Jet
     929
     930        add Branch GenMissingET/momentum GenMissingET MissingET
    1164931
    1165932    add Branch TrackMerger/tracks Track Track
Note: See TracChangeset for help on using the changeset viewer.