Fork me on GitHub

source: svn/trunk/examples/delphes_card_CMS_PileUp.tcl@ 1190

Last change on this file since 1190 was 1139, checked in by Pavel Demin, 11 years ago

fix tau-tagging input collection
File size: 20.5 KB
Line 
1#######################################
2# Order of execution of various modules
3#######################################
4
5set ExecutionPath {
6
7 Weighter
8
9 PileUpMerger
10 ParticlePropagator
11
12 ChargedHadronTrackingEfficiency
13 ElectronTrackingEfficiency
14 MuonTrackingEfficiency
15
16 ChargedHadronMomentumSmearing
17 ElectronEnergySmearing
18 MuonMomentumSmearing
19
20 TrackMerger
21 Calorimeter
22 TrackPileUpSubtractor
23 EFlowMerger
24
25 GenJetFinder
26
27 Rho
28 FastJetFinder
29 JetPileUpSubtractor
30
31 ConstituentFilter
32
33 PhotonEfficiency
34 PhotonIsolation
35
36 ElectronEfficiency
37 ElectronIsolation
38
39 MuonEfficiency
40 MuonIsolation
41
42 MissingET
43
44 BTagging
45 TauTagging
46
47 UniqueObjectFinder
48
49 ScalarHT
50
51 TreeWriter
52}
53
54##########
55# Weighter
56##########
57
58module Weighter Weighter {
59 set InputArray Delphes/allParticles
60
61 set OutputArray weight
62
63 # add Weight {PID1 PID2} {weight}
64 # default weight
65 add Weight {0} {1.0}
66
67 add Weight {23} {0.5}
68 add Weight {24} {0.5}
69 add Weight {-24} {0.5}
70 add Weight {23 24} {0.25}
71 add Weight {23 -24} {0.25}
72}
73
74###############
75# PileUp Merger
76###############
77
78module PileUpMerger PileUpMerger {
79 set InputArray Delphes/stableParticles
80
81 set OutputArray stableParticles
82
83 # pre-generated minbias input file
84 set PileUpFile MinBias.pileup
85
86 # average expected pile up
87 set MeanPileUp 50
88 # spread in the beam direction in m (assumes gaussian)
89 set ZVertexSpread 0.07
90}
91
92#################################
93# Propagate particles in cylinder
94#################################
95
96module ParticlePropagator ParticlePropagator {
97 set InputArray PileUpMerger/stableParticles
98
99 set OutputArray stableParticles
100 set ChargedHadronOutputArray chargedHadrons
101 set ElectronOutputArray electrons
102 set MuonOutputArray muons
103
104 # radius of the magnetic field coverage, in m
105 set Radius 1.29
106 # half-length of the magnetic field coverage, in m
107 set HalfLength 3.00
108
109 # magnetic field
110 set Bz 3.8
111}
112
113####################################
114# Charged hadron tracking efficiency
115####################################
116
117module Efficiency ChargedHadronTrackingEfficiency {
118 set InputArray ParticlePropagator/chargedHadrons
119 set OutputArray chargedHadrons
120
121 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
122
123 # tracking efficiency formula for charged hadrons
124 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
125 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
126 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + \
127 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + \
128 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) + \
129 (abs(eta) > 2.5) * (0.00)}
130}
131
132##############################
133# Electron tracking efficiency
134##############################
135
136module Efficiency ElectronTrackingEfficiency {
137 set InputArray ParticlePropagator/electrons
138 set OutputArray electrons
139
140 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
141
142 # tracking efficiency formula for electrons
143 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
144 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + \
145 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + \
146 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + \
147 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + \
148 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + \
149 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + \
150 (abs(eta) > 2.5) * (0.00)}
151}
152
153##########################
154# Muon tracking efficiency
155##########################
156
157module Efficiency MuonTrackingEfficiency {
158 set InputArray ParticlePropagator/muons
159 set OutputArray muons
160
161 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
162
163 # tracking efficiency formula for muons
164 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
165 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \
166 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \
167 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
168 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + \
169 (abs(eta) > 2.5) * (0.00)}
170}
171
172########################################
173# Momentum resolution for charged tracks
174########################################
175
176module MomentumSmearing ChargedHadronMomentumSmearing {
177 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
178 set OutputArray chargedHadrons
179
180 # set ResolutionFormula {resolution formula as a function of eta and pt}
181
182 # resolution formula for charged hadrons
183 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
184 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \
185 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
186 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
187 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
188 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \
189 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
190 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)}
191}
192
193#################################
194# Energy resolution for electrons
195#################################
196
197module EnergySmearing ElectronEnergySmearing {
198 set InputArray ElectronTrackingEfficiency/electrons
199 set OutputArray electrons
200
201 # set ResolutionFormula {resolution formula as a function of eta and energy}
202
203 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.5e1) * (energy*0.015) + \
204 (abs(eta) <= 2.5) * (energy > 2.5e1) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
205 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
206 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
207
208}
209
210###############################
211# Momentum resolution for muons
212###############################
213
214module MomentumSmearing MuonMomentumSmearing {
215 set InputArray MuonTrackingEfficiency/muons
216 set OutputArray muons
217
218 # set ResolutionFormula {resolution formula as a function of eta and pt}
219
220 # resolution formula for muons
221 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
222 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \
223 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
224 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
225 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.04) + \
226 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.03) + \
227 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
228 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)}
229}
230
231##############
232# Track merger
233##############
234
235module Merger TrackMerger {
236# add InputArray InputArray
237 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
238 add InputArray ElectronEnergySmearing/electrons
239 set OutputArray tracks
240}
241
242#############
243# Calorimeter
244#############
245
246module Calorimeter Calorimeter {
247 set ParticleInputArray ParticlePropagator/stableParticles
248 set TrackInputArray TrackMerger/tracks
249
250 set TowerOutputArray towers
251 set PhotonOutputArray photons
252
253 set EFlowTrackOutputArray eflowTracks
254 set EFlowTowerOutputArray eflowTowers
255
256 set pi [expr {acos(-1)}]
257
258 # lists of the edges of each tower in eta and phi
259 # each list starts with the lower edge of the first tower
260 # the list ends with the higher edged of the last tower
261
262 # 5 degrees towers
263 set PhiBins {}
264 for {set i -36} {$i <= 36} {incr i} {
265 add PhiBins [expr {$i * $pi/36.0}]
266 }
267 foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} {
268 add EtaPhiBins $eta $PhiBins
269 }
270
271 # 10 degrees towers
272 set PhiBins {}
273 for {set i -18} {$i <= 18} {incr i} {
274 add PhiBins [expr {$i * $pi/18.0}]
275 }
276 foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} {
277 add EtaPhiBins $eta $PhiBins
278 }
279
280 # 20 degrees towers
281 set PhiBins {}
282 for {set i -9} {$i <= 9} {incr i} {
283 add PhiBins [expr {$i * $pi/9.0}]
284 }
285 foreach eta {-5 -4.7 -4.525 4.7 5} {
286 add EtaPhiBins $eta $PhiBins
287 }
288
289 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
290 add EnergyFraction {0} {0.0 1.0}
291 # energy fractions for e, gamma and pi0
292 add EnergyFraction {11} {1.0 0.0}
293 add EnergyFraction {22} {1.0 0.0}
294 add EnergyFraction {111} {1.0 0.0}
295 # energy fractions for muon, neutrinos and neutralinos
296 add EnergyFraction {12} {0.0 0.0}
297 add EnergyFraction {13} {0.0 0.0}
298 add EnergyFraction {14} {0.0 0.0}
299 add EnergyFraction {16} {0.0 0.0}
300 add EnergyFraction {1000022} {0.0 0.0}
301 add EnergyFraction {1000023} {0.0 0.0}
302 add EnergyFraction {1000025} {0.0 0.0}
303 add EnergyFraction {1000035} {0.0 0.0}
304 add EnergyFraction {1000045} {0.0 0.0}
305 # energy fractions for K0short and Lambda
306 add EnergyFraction {310} {0.3 0.7}
307 add EnergyFraction {3122} {0.3 0.7}
308
309 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
310 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
311 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
312
313 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
314 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + \
315 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
316}
317
318##########################
319# Track pile-up subtractor
320##########################
321
322module TrackPileUpSubtractor TrackPileUpSubtractor {
323# add InputArray InputArray OutputArray
324 add InputArray Calorimeter/eflowTracks eflowTracks
325 add InputArray ElectronEnergySmearing/electrons electrons
326 add InputArray MuonMomentumSmearing/muons muons
327
328 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
329 # Z vertex resolution in m
330 set ZVertexResolution 0.0001
331}
332
333####################
334# Energy flow merger
335####################
336
337module Merger EFlowMerger {
338# add InputArray InputArray
339 add InputArray TrackPileUpSubtractor/eflowTracks
340 add InputArray Calorimeter/eflowTowers
341 add InputArray MuonMomentumSmearing/muons
342 set OutputArray eflow
343}
344
345#############
346# Rho pile-up
347#############
348
349module FastJetFinder Rho {
350# set InputArray Calorimeter/towers
351 set InputArray EFlowMerger/eflow
352
353 set ComputeRho true
354 set RhoOutputArray rho
355
356 # area algorithm: 0 Do not compute area, 1 Active area explicit ghosts, 2 One ghost passive area, 3 Passive area, 4 Voronoi, 5 Active area
357 set AreaAlgorithm 5
358
359 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
360 set JetAlgorithm 4
361 set ParameterR 0.6
362 set GhostEtaMax 5.0
363 set RhoEtaMax 5.0
364
365 set JetPTMin 0.0
366}
367
368#####################
369# MC truth jet finder
370#####################
371
372module FastJetFinder GenJetFinder {
373 set InputArray Delphes/stableParticles
374
375 set OutputArray jets
376
377 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
378 set JetAlgorithm 6
379 set ParameterR 0.5
380
381 set JetPTMin 20.0
382}
383
384############
385# Jet finder
386############
387
388module FastJetFinder FastJetFinder {
389# set InputArray Calorimeter/towers
390 set InputArray EFlowMerger/eflow
391
392 set OutputArray jets
393
394 # area algorithm: 0 Do not compute area, 1 Active area explicit ghosts, 2 One ghost passive area, 3 Passive area, 4 Voronoi, 5 Active area
395 set AreaAlgorithm 5
396
397 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
398 set JetAlgorithm 6
399 set ParameterR 0.5
400
401 set JetPTMin 20.0
402}
403
404####################
405# Constituent filter
406####################
407
408module ConstituentFilter ConstituentFilter {
409
410# add JetInputArray InputArray
411 add JetInputArray GenJetFinder/jets
412 add JetInputArray FastJetFinder/jets
413
414# add ConstituentInputArray InputArray OutputArray
415 add ConstituentInputArray Delphes/stableParticles stableParticles
416 add ConstituentInputArray TrackPileUpSubtractor/eflowTracks eflowTracks
417 add ConstituentInputArray Calorimeter/eflowTowers eflowTowers
418 add ConstituentInputArray MuonMomentumSmearing/muons muons
419}
420
421###########################
422# Jet Pile-Up Subtraction
423###########################
424
425module JetPileUpSubtractor JetPileUpSubtractor {
426 set JetInputArray FastJetFinder/jets
427 set RhoInputArray Rho/rho
428
429 set OutputArray jets
430
431 set JetPTMin 20.0
432}
433
434###################
435# Photon efficiency
436###################
437
438module Efficiency PhotonEfficiency {
439 set InputArray Calorimeter/photons
440 set OutputArray photons
441
442 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
443
444 # efficiency formula for photons
445 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
446 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
447 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
448 (abs(eta) > 2.5) * (0.00)}
449}
450
451##################
452# Photon isolation
453##################
454
455module Isolation PhotonIsolation {
456 set CandidateInputArray PhotonEfficiency/photons
457 set IsolationInputArray EFlowMerger/eflow
458 set RhoInputArray Rho/rho
459
460 set OutputArray photons
461
462 set DeltaRMax 0.5
463
464 set PTMin 0.5
465
466 set PTRatioMax 0.1
467}
468
469#####################
470# Electron efficiency
471#####################
472
473module Efficiency ElectronEfficiency {
474 set InputArray TrackPileUpSubtractor/electrons
475 set OutputArray electrons
476
477 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
478
479 # efficiency formula for electrons
480 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
481 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
482 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
483 (abs(eta) > 2.5) * (0.00)}
484}
485
486####################
487# Electron isolation
488####################
489
490module Isolation ElectronIsolation {
491 set CandidateInputArray ElectronEfficiency/electrons
492 set IsolationInputArray EFlowMerger/eflow
493 set RhoInputArray Rho/rho
494
495 set OutputArray electrons
496
497 set DeltaRMax 0.5
498
499 set PTMin 0.5
500
501 set PTRatioMax 0.1
502}
503
504#################
505# Muon efficiency
506#################
507
508module Efficiency MuonEfficiency {
509 set InputArray TrackPileUpSubtractor/muons
510 set OutputArray muons
511
512 # set EfficiencyFormula {efficiency as a function of eta and pt}
513
514 # efficiency formula for muons
515 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
516 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
517 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0) * (0.85) + \
518 (abs(eta) > 2.4) * (0.00)}
519}
520
521################
522# Muon isolation
523################
524
525module Isolation MuonIsolation {
526 set CandidateInputArray MuonEfficiency/muons
527 set IsolationInputArray EFlowMerger/eflow
528 set RhoInputArray Rho/rho
529
530 set OutputArray muons
531
532 set DeltaRMax 0.5
533
534 set PTMin 0.5
535
536 set PTRatioMax 0.1
537}
538
539###################
540# Missing ET merger
541###################
542
543module Merger MissingET {
544# add InputArray InputArray
545 add InputArray EFlowMerger/eflow
546 set MomentumOutputArray momentum
547}
548
549##################
550# Scalar HT merger
551##################
552
553module Merger ScalarHT {
554# add InputArray InputArray
555 add InputArray UniqueObjectFinder/jets
556 add InputArray UniqueObjectFinder/electrons
557 add InputArray UniqueObjectFinder/photons
558 add InputArray MuonIsolation/muons
559 set EnergyOutputArray energy
560}
561
562
563###########
564# b-tagging
565###########
566
567module BTagging BTagging {
568 set PartonInputArray Delphes/partons
569 set JetInputArray JetPileUpSubtractor/jets
570
571 set BitNumber 0
572
573 set DeltaR 0.5
574
575 set PartonPTMin 1.0
576
577 set PartonEtaMax 2.5
578
579 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
580 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
581 # gluon's PDG code has the lowest priority
582
583 # https://twiki.cern.ch/twiki/bin/view/CMSPublic/PhysicsResultsBTV
584 # default efficiency formula (misidentification rate)
585 add EfficiencyFormula {0} {0.001}
586
587 # efficiency formula for c-jets (misidentification rate)
588 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \
589 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
590 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
591 (abs(eta) > 2.5) * (0.000)}
592
593 # efficiency formula for b-jets
594 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \
595 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
596 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
597 (abs(eta) > 2.5) * (0.000)}
598}
599
600module TauTagging TauTagging {
601 set ParticleInputArray Delphes/allParticles
602 set PartonInputArray Delphes/partons
603 set JetInputArray JetPileUpSubtractor/jets
604
605 set DeltaR 0.5
606
607 set TauPTMin 1.0
608
609 set TauEtaMax 2.5
610
611 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
612
613 # default efficiency formula (misidentification rate)
614 add EfficiencyFormula {0} {0.001}
615 # efficiency formula for tau-jets
616 add EfficiencyFormula {15} {0.4}
617}
618
619#####################################################
620# Find uniquely identified photons/electrons/tau/jets
621#####################################################
622
623module UniqueObjectFinder UniqueObjectFinder {
624# earlier arrays take precedence over later ones
625# add InputArray InputArray OutputArray
626 add InputArray PhotonIsolation/photons photons
627 add InputArray ElectronIsolation/electrons electrons
628 add InputArray JetPileUpSubtractor/jets jets
629}
630
631##################
632# ROOT tree writer
633##################
634
635module TreeWriter TreeWriter {
636# add Branch InputArray BranchName BranchClass
637 add Branch Delphes/allParticles Particle GenParticle
638 add Branch TrackMerger/tracks Track Track
639 add Branch Calorimeter/towers Tower Tower
640 add Branch ConstituentFilter/eflowTracks EFlowTrack Track
641 add Branch ConstituentFilter/eflowTowers EFlowTower Tower
642 add Branch ConstituentFilter/muons EFlowMuon Muon
643 add Branch GenJetFinder/jets GenJet Jet
644 add Branch UniqueObjectFinder/jets Jet Jet
645 add Branch UniqueObjectFinder/electrons Electron Electron
646 add Branch UniqueObjectFinder/photons Photon Photon
647 add Branch MuonIsolation/muons Muon Muon
648 add Branch MissingET/momentum MissingET MissingET
649 add Branch ScalarHT/energy ScalarHT ScalarHT
650 add Branch Rho/rho Rho Rho
651 add Branch Weighter/weight Weight Weight
652}
653
Note: See TracBrowser for help on using the repository browser.