|
|---|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Contacts:
|
Subjet analysis with KT algorithm
jobOptions to select
ClusterMaker.Strategy = "Kt"; //Sets clustering stratgy to KT
ClusterMaker.KtRParameter = 1.0; //Sets R-Parameter of Kt algorithm
ClusterMaker.KtAngle = "deltaR"; //Sets distance scheme
ClusterMaker.KtRecomScheme = "E"; //Sets recombination scheme
ClusterMaker.KtYCut = "0" //Sets number of jets merged when y-cut of cluster is set
To run the subjet anaylsis you need to set KtYCut to non-zero. This will then calculate the y-cut value at the point when the clusters components are merged into the KtYCut number of subjets. This value is then stored by the cluster.
ClusterKtStratgy makes Atlfast::KtClusters which honour the IKtCluster and ICluster interfaces. These clusters are however stored by ClusterMaker as an ICellCollection. This is a collection of ICells from which you cannot retrieve the y-cut value. Thus you must retrieve the origanal object using for example TypeVisitor that can be found in AtlfastUtils.
Once you have the KtCluster you can retrieve the ycut value with KtCluster::yCut()
Note: JetSmearer was parameterised for the cone algorithm, so should not be used with the Kt algorithm. Instead, JetMaker.DoSmearing should be set to false and either the shower parameterisation OR the simplified CellMaker.DoSmearing should be used.