本文整理汇总了C++中TLorentzVector::DrEtaPhi方法的典型用法代码示例。如果您正苦于以下问题:C++ TLorentzVector::DrEtaPhi方法的具体用法?C++ TLorentzVector::DrEtaPhi怎么用?C++ TLorentzVector::DrEtaPhi使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类TLorentzVector
的用法示例。
在下文中一共展示了TLorentzVector::DrEtaPhi方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: fill
//.........这里部分代码省略.........
double max4Pt = 0.0;
double maxPt = 0.0;
double M3first = 0.0;
TLorentzVector maxPtsystem;
TLorentzVector phovec;
phovec.SetPtEtaPhiM(0.00001,0.0,0.0,0.0);
if( selEvent->Photons.size() > 0 ){
int phoi = selEvent->Photons[0];
phovec.SetPtEtaPhiM(tree->phoEt_->at(phoi), tree->phoEta_->at(phoi), tree->phoPhi_->at(phoi), 0.0);
}
for(int jet1I=0; jet1I < selEvent->Jets.size()-2; jet1I++){
jetI = selEvent->Jets[jet1I];
j1.SetPtEtaPhiM(tree->jetPt_->at(jetI), tree->jetEta_->at(jetI), tree->jetPhi_->at(jetI), 0.0);
for(int jet2I=jet1I+1; jet2I < selEvent->Jets.size()-1; jet2I++){
jetI = selEvent->Jets[jet2I];
j2.SetPtEtaPhiM(tree->jetPt_->at(jetI), tree->jetEta_->at(jetI), tree->jetPhi_->at(jetI), 0.0);
for(int jet3I=jet2I+1; jet3I < selEvent->Jets.size(); jet3I++){
jetI = selEvent->Jets[jet3I];
j3.SetPtEtaPhiM(tree->jetPt_->at(jetI), tree->jetEta_->at(jetI), tree->jetPhi_->at(jetI), 0.0);
double m3 = (j1+j2+j3).M();
double totalPt = (j1+j2+j3).Pt();
if(jet1I==0 && jet2I==1 && jet3I==2) M3first = m3;
if(m3 < minM3) minM3 = m3;
if(minPt > totalPt){
minPt = totalPt;
M3minPt = m3;
}
if(maxPt < totalPt){
maxPt = totalPt;
M3maxPt = m3;
maxPtsystem = (j1+j2+j3);
}
if( phovec.DrEtaPhi(j1) < 0.3 )
j1 = j1 - phovec;
if( phovec.DrEtaPhi(j2) < 0.3 )
j2 = j2 - phovec;
if( phovec.DrEtaPhi(j3) < 0.3)
j3 = j3 - phovec;
double m4 = (phovec+j1+j2+j3).M();
double total4Pt = (phovec+j1+j2+j3).Pt();
hists["M3phoMulti"]->Fill(m4, weight);
if(max4Pt < total4Pt ){
max4Pt=total4Pt;
M4maxPt=m4;
}
}
}
}
double toppt=0.0;
double antitoppt=0.0;
for(int mcInd=0; mcInd<tree->nMC_; ++mcInd){
if(tree->mcPID->at(mcInd)==6) toppt = tree->mcPt->at(mcInd);
if(tree->mcPID->at(mcInd)==-6) antitoppt = tree->mcPt->at(mcInd);
}
double maxtoppt = std::max(toppt,antitoppt);
if( maxtoppt < 30 ) hists["M3_0_30"]->Fill(M3maxPt, weight);
else if( maxtoppt < 100) hists["M3_30_100"]->Fill(M3maxPt, weight);
else if( maxtoppt < 200) hists["M3_100_200"]->Fill(M3maxPt, weight);
else if( maxtoppt < 300) hists["M3_200_300"]->Fill(M3maxPt, weight);
else hists["M3_300_up"]->Fill(M3maxPt, weight);
hists["M3first"]->Fill(M3first, weight);
hists["M3"]->Fill(M3maxPt, weight);
hists["M3minPt"]->Fill(M3minPt, weight);
if( selEvent->Photons.size() > 0 ) {
hists["M3pho"]->Fill(M4maxPt, weight);
hists["dRpho3j"]->Fill(phovec.DrEtaPhi(maxPtsystem), weight);
}
hists["minM3"]->Fill(minM3, weight);
hists2d["MTW_M3"]->Fill( MTW, M3maxPt, weight);
}
if( selEvent->Jets.size() > 0 ){
int ind = selEvent->Jets[0];
hists["jet1Pt"]->Fill( tree->jetPt_->at(ind), weight );
hists["jet1Eta"]->Fill( tree->jetEta_->at(ind), weight );
}
if( selEvent->Jets.size() > 1 ){
int ind = selEvent->Jets[1];
hists["jet2Pt"]->Fill( tree->jetPt_->at(ind), weight );
hists["jet2Eta"]->Fill( tree->jetEta_->at(ind), weight );
}
if( selEvent->Jets.size() > 2 ){
int ind = selEvent->Jets[2];
hists["jet3Pt"]->Fill( tree->jetPt_->at(ind), weight );
hists["jet3Eta"]->Fill( tree->jetEta_->at(ind), weight );
}
if( selEvent->Jets.size() > 3 ){
int ind = selEvent->Jets[3];
hists["jet4Pt"]->Fill( tree->jetPt_->at(ind), weight );
hists["jet4Eta"]->Fill( tree->jetEta_->at(ind), weight );
}
}