Draft: [QEE] Wcharm line
NB: Closure
CLOSED. Superseded by: !3217 (merged)
See: discussion/reasoning !3152 (comment 7774158)
Development of asymmetric dijet lines targeting W->Charm production, but also general use for SVTag + Inclusive line for other future physics.
Specific line rates:
MinBias efficiency:
LAZY_AND: Hlt2QEE_DiJetIncSVTag_pT10M40dR20Full #=25000 Sum=31 Eff=|(0.1240000 +- 0.0222572)%|
Signal MC efficiency:
LAZY_AND: Hlt2QEE_DiJetIncSVTag_pT10M40dR20Full #=5843 Sum=278 Eff=|( 4.757830 +- 0.278485)%|
pT
, di-jet mass, and dR
tuning:
WcharmInfo.txt
With ALL Changes:
Ran over 25000 events.
Fullstream filesize = 20 MB. Fullstream Events fired = 157
Fullstream Rate = 1.14 MHz * 157 / 25,000 = 7.159 kHz
Fullstream BW = 1.14 MHz * 20 MB / 25,000 = 912 MB/s
Turbo filesize = 8.1 MB. 1758 events passed (overlap 58)
**Note: ** Rates here will have non-negligible overlap with di-jet lines, and likely cannot be reasonably measured independently
With ONLY PV matching on all di-jet lines:
Ran over 5000 events.
Fullstream filesize = 3.1 MB. Fullstream Events fired = 27
Fullstream Rate = 1.14 MHz * 27 / 5,000 = 6.156 kHz
Fullstream BW = 1.14 MHz * 3.1 MB / 5,000 = 706.8 MB/s
Turbo filesize = 1.7 MB. 346 events passed (overlap 1)
Without Changes:
Ran over 5000 events.
Fullstream filesize = 3.6 MB. Fullstream Events fired = 28
Fullstream Rate = 1.14 MHz * 28 / 5,000 = 6.384 kHz
Fullstream BW = 1.14 MHz * 3.6 MB / 5,000 = 843.6 MB/s
Turbo filesize = 1.7 MB. 346 events passed (overlap 1)
Merge request reports
Activity
added 2024data-selections qee labels
assigned to @ngrieser
added RTA label
mentioned in issue #718
mentioned in merge request lhcb-simulation/mc-requests!478 (merged)
mentioned in merge request !3151 (merged)
207 220 calo_clusters=True, 208 221 persistreco=persistreco) 209 222 223 ############ Asymmetric dijet lines ################ 224 @register_line_builder(all_lines) 225 @configurable 226 def DiJetIncSVTag_line(name='Hlt2QEE_DiJetIncSVTag', 227 prescale=1, 228 persistreco=True): 229 jets = make_dijets( 230 tagpair=('SV', 'Inc'), prod_pt_min=30 * GeV, min_dphi=0.0) changed this line in version 3 of the diff
For completeness: the question was
Is this additional line doing what I intend physically (One SVtag jet, one inclusive)? It "works" but I am not fully certain that it's behaving exactly how I'd hope.
make_dijets('SV',None): jet1 = make_jets(pt_min=prod_pt_min, tags='SV') jet2 = make_jets(pt_min=prod_pt_min, tags=None) # ... physics related code return ParticleCombiner( Inputs=[jet1, jet2], DecayDescriptor="CLUSjet -> CELLjet CELLjet", CombinationCut=combination_code, CompositeCut=F.ALL, AllowDiffInputsForSameIDChildren = True, ParticleCombiner="ParticleAdder")
This format you have, summarised by me above, looks like it is generating a DiJet of ['SV' jet, None Jet] in your new line. Which if I understand your meaning is exactly what you want to do.
AllowDiffInputsForSameIdChildren=True
is needed, because they are bothCELLjet
, but are different containers. As I understand this wasn't being done up until now, i.e. the jets were always 'symmetric'.I'm not sure if there's a way you can verify with Data as I'm not a jet expert, but if there is I'd recommend trying it locally and convincing yourself this looks physically correct. But to me this looks like the correct way to perform 'asymmetric' Jet combinations.
Edited by Luke Grazette
added 1 commit
- 00afbe76 - Add allowing different inputs to support asymmetric dijet combinations.
@rjhunter since this will have a lot of overlap with inclusive di-jet lines (waiting on Allen!1372 (merged)), would you want them all here to get a more reasonable first-order bw estimate, or wait for it all in !3151 (merged)
added 36 commits
- 00afbe76...c01cc4be - 26 earlier commits
- 3dac80ad - Update References for: Allen!1436 (merged) based on lhcb-master-mr/11117
- 4941e648 - set pv_tracks true for persistreco
- 1407f576 - Update References for: !3103 (merged) based on lhcb-master-mr/11140
- e06b010b - Update References for: Allen!1446 (merged) based on lhcb-master-mr/11152
- e0abec0b - decreports filter without functors
- ffeeb2d3 - Update References for: LHCb!4454 (merged), !3054 (merged) based on lhcb-master-mr/11156
- a7632c4d - Update BandQ inclusive line data processing flow, turbo PersistReco -> full...
- 66ea35bf - adapt to removal of additionalInfo for calo and global pid
- bb41a9ab - Update References for: LHCb!4338 (merged), Rec!3740 (merged), !3149 (merged), DaVinci!1050 (merged),...
- a0475889 - Add dijet mass cut to suppress non-resonant contribution.
Toggle commit list