本文整理汇总了C++中HiForest::SelectJetAlgo方法的典型用法代码示例。如果您正苦于以下问题：C++ HiForest::SelectJetAlgo方法的具体用法？C++ HiForest::SelectJetAlgo怎么用？C++ HiForest::SelectJetAlgo使用的例子？那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类HiForest的用法示例。
在下文中一共展示了HiForest::SelectJetAlgo方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: IndResponse
int IndResponse(double kPt=80,const char *ksp="pbpb")
{
timer.Start();
//! Load Lib
gSystem->Load("hiForest_h.so");
//! Load the hiforest input root file
HiForest *c = 0;
if(strcmp("pp",ksp)==0)c = new HiForest(Form("/net/hisrv0001/home/icali/hadoop/Pythia/Z2/ppDijet_merged/pp276Dijet%0.0f_merged.root",kPt),"pp2012",1,1);
else {
if(kPt==30)c = new HiForest("/d102/yjlee/hiForest2MC/Pythia30_HydjetDrum_mix01_HiForest2_v19.root","pbpb2012",0,1);
else if(kPt==50)c = new HiForest("/d102/yjlee/hiForest2MC/Pythia50_HydjetDrum_mix01_HiForest2_v19.root","pbpb2012",0,1);
else if(kPt==80)c = new HiForest("/d102/yjlee/hiForest2MC/Pythia80_HydjetDrum_mix01_HiForest2_v20.root","pbpb2012",0,1);
else if(kPt==120)c = new HiForest("/d102/yjlee/hiForest2MC/Pythia120_HydjetDrum_mix01_HiForest2_v21_ivan.root","pbpb2012",0,1);
else if(kPt==170)c = new HiForest("/d102/yjlee/hiForest2MC/Pythia170_HydjetDrum_mix01_HiForest2_v19.root","pbpb2012",0,1);
else if(kPt==200)c = new HiForest("/d102/yjlee/hiForest2MC/Pythia200_HydjetDrum_mix01_HiForest2_v21_ivan.root","pbpb2012",0,1);
else if(kPt==250)c = new HiForest("/d102/yjlee/hiForest2MC/Pythia250_HydjetDrum_mix01_HiForest2_v21_ivan.root","pbpb2012",0,1);
//c = new HiForest(Form("/net/hisrv0001/home/icali/hadoop/Hydjet1.8/Z2/Dijet_merged/Dijet%0.0f_merged.root",kPt),"pbpb2012",0,1);
}
double xsection=0;
double xup=0;
double xsub=0;
double maxpthat=9999;
if(kPt==15){
maxpthat=30;
xup=1.566e-01;
xsub=1.079e-02;
}
else if(kPt==30){
maxpthat=50;
xup=1.079e-02;
xsub=1.021e-03;
}
else if(kPt==50){
maxpthat=80;
xup=1.021e-03;
xsub=9.913e-05;
}
else if(kPt==80){
maxpthat=120;
xup=9.913e-05;
xsub=1.128e-05;
}
else if(kPt==120){
maxpthat=170;
xup=1.128e-05;
xsub=1.470e-06;
}
else if(kPt==170){
maxpthat=200;
xup=1.470e-06;
xsub=5.310e-07;
}
else if(kPt==200){
maxpthat=250;
xup=5.310e-07;
xsub=1.192e-07;
}
else if(kPt==250){
maxpthat=300;
xup=1.192e-07;
xsub=3.176e-08;
}
else if(kPt==300){
maxpthat=9999;
xup=3.176e-08;
xsub=0;
}
xsection = xup-xsub;
std::cout<<std::endl;
std::cout<<std::endl;
//! Don't want to loop over trees which is not used in the analysis
//! event trees
//c->hasEvtTree=0;
//! jet trees
//! Switch on only the jet trees which you require
const char *cjets[7] = {"icPu5","ak2PF","ak3PF","ak4PF","akPu2PF","akPu3PF","akPu4PF"};
c->SelectJetAlgo(cjets,7);
//! photon tree
//c->hasPhotonTree=0;
//! Select only the jet branches which you are going to use
//! This increases the speed for running over trees with lot of branches.
//! This is currently for only jet algos and evtTree uniformly applied to all the
//! Event Tree
const char *evlist[]={"hiBin","vz","hiHF"};
const int kevbr = sizeof(evlist)/sizeof(const char *);
//.........这里部分代码省略.........