本文整理汇总了C++中Jet::getWeightedDepthNoEE方法的典型用法代码示例。如果您正苦于以下问题:C++ Jet::getWeightedDepthNoEE方法的具体用法?C++ Jet::getWeightedDepthNoEE怎么用?C++ Jet::getWeightedDepthNoEE使用的例子?那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类Jet
的用法示例。
在下文中一共展示了Jet::getWeightedDepthNoEE方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: fillPlots
void fillPlots ( CollectionPtr & jets,
CollectionPtr & core_jets,
CollectionPtr & partons,
CollectionPtr & hgceeClusters,
CollectionPtr & hgchefClusters,
CollectionPtr & hgchebClusters,
likelihoodGetter & l,
int index,
std::vector<TH1F*> quark_th1, std::vector<TH1F*> gluon_th1,
std::vector<TH2F*> quark_th2, std::vector<TH2F*> gluon_th2 ){
int n_jets = jets -> GetSize();
std::vector<double> likelihood_variables;
for (int i = 0; i < n_jets; ++i){
Jet jet = jets -> GetConstituent<Jet>(i);
CollectionPtr all_core_jets = core_jets;
CollectionPtr matched_core_jets = core_jets -> SkimByRequireDRMatch<PFCoreJet, PFPrunedJet>( jet, 0.3 ) ;
rechit_type my_type = NO_TYPE;
int my_index = -1;
jet.getLeadClusterTypeAndIndex(my_type, my_index);
double clusterLength = -1.;
double clusterVolume = -1.;
if ( my_type == NO_TYPE ) continue;
if ( my_type == HGCEE ){
HGCEECluster c(*hgceeClusters, my_index);
clusterLength = c.getLength();
clusterVolume = c.getVolume();
}
else if ( my_type == HGCHEF ){
HGCHEFCluster c(*hgchefClusters, my_index);
clusterLength = c.getLength();
clusterVolume = c.getVolume();
}
else {
HGCHEBCluster c(*hgchebClusters, my_index);
clusterLength = c.getLength();
clusterVolume = c.getVolume();
}
double core_over_total_pt = -1.;
double lead_core_jet_dr = -1.;
double n_core_jets = matched_core_jets -> GetSize();
if ( n_core_jets > 0 ){
PFCoreJet lead_core_jet = matched_core_jets -> GetLeadPtObject<PFCoreJet> ();
double core_pt = lead_core_jet.Pt();
core_over_total_pt = core_pt / jet.Pt();
lead_core_jet_dr = jet.DeltaR( &lead_core_jet );
}
int flavor = getFlavor ( jet , partons );
if ( flavor == -1 ) continue;
bool isQuarkJet ( flavor == 1 );
bool isGluonJet ( flavor == 0 );
likelihood_variables.clear();
double weighted_depth = jet.getWeightedDepth();
double weighted_depth_noEE = jet.getWeightedDepthNoEE();
double profile[n_radii];
jet.getProfile(n_radii, radii, profile);
double profile50 = jet.getProfileRadius(n_radii, radii, profile, 0.5);
double profile90 = jet.getProfileRadius(n_radii, radii, profile, 0.9);
double profile95 = jet.getProfileRadius(n_radii, radii, profile, 0.95);
double profile99 = jet.getProfileRadius(n_radii, radii, profile, 0.99);
likelihood_variables.push_back ( jet.TrimmedNSubjets(index) );
likelihood_variables.push_back ( jet.getWidth() );
likelihood_variables.push_back ( jet.getNPFCandidates() );
likelihood_variables.push_back ( jet.getPTD() );
likelihood_variables.push_back ( jet.NSubJettiness() );
likelihood_variables.push_back ( jet.TrimmedMass(index) * 1000. );
likelihood_variables.push_back ( profile50 );
likelihood_variables.push_back ( core_over_total_pt );
likelihood_variables.push_back ( lead_core_jet_dr );
double ee_energy = jet.getEEEnergy();
double ee015_energy = jet.getEEEnergy(0, 15);
double ee1631_energy = jet.getEEEnergy(16, 31);
double hef_energy = jet.getHEFEnergy();
double heb_energy = jet.getHEBEnergy();
double all_energy = ee_energy + hef_energy + heb_energy;
double hef11_energy = jet.getHEFEnergy(1, 1);
double hef22_energy = jet.getHEFEnergy(2, 2);
double hef34_energy = jet.getHEFEnergy(3, 4);
double hef56_energy = jet.getHEFEnergy(5, 6);
double hef712_energy = jet.getHEFEnergy(7,12);
double likelihood = l.getLikelihood ( likelihood_variables );
//.........这里部分代码省略.........