Double_t cosp=d->CosPointingAngle();
Double_t sumD02=d->Getd0Prong(0)*d->Getd0Prong(0)+d->Getd0Prong(1)*d->Getd0Prong(1)+d->Getd0Prong(2)*d->Getd0Prong(2);
Double_t dca=d->GetDCA();
+ Double_t sigvert=d->GetSigmaVert();
Double_t ptmax=0;
for(Int_t i=0;i<3;i++){
if(d->PtProng(i)>ptmax)ptmax=d->PtProng(i);
fCosPHistLS[indexcut]->Fill(cosp);
fDLenHistLS[indexcut]->Fill(dlen);
fSumd02HistLS[indexcut]->Fill(sumD02);
+ fSigVertHistLS[indexcut]->Fill(sigvert);
fPtMaxHistLS[indexcut]->Fill(ptmax);
fDCAHistLS[indexcut]->Fill(dca);
Double_t dlen=d->DecayLength();
Double_t cosp=d->CosPointingAngle();
Double_t sumD02=d->Getd0Prong(0)*d->Getd0Prong(0)+d->Getd0Prong(1)*d->Getd0Prong(1)+d->Getd0Prong(2)*d->Getd0Prong(2);
- Double_t dca=d->GetDCA();
-Double_t ptmax=0;
+ Double_t dca=d->GetDCA();
+ Double_t sigvert=d->GetSigmaVert();
+ Double_t ptmax=0;
for(Int_t i=0;i<3;i++){
if(d->PtProng(i)>ptmax)ptmax=d->PtProng(i);
}
fCosPHist[index]->Fill(cosp);
fDLenHist[index]->Fill(dlen);
fSumd02Hist[index]->Fill(sumD02);
+ fSigVertHist[index]->Fill(sigvert);
fPtMaxHist[index]->Fill(ptmax);
fDCAHist[index]->Fill(dca);
fCosPHist[index]->Fill(cosp);
fDLenHist[index]->Fill(dlen);
fSumd02Hist[index]->Fill(sumD02);
+ fSigVertHist[index]->Fill(sigvert);
fPtMaxHist[index]->Fill(ptmax);
fDCAHist[index]->Fill(dca);
if(passTightCuts==1){
fCosPHist[index]->Fill(cosp);
fDLenHist[index]->Fill(dlen);
fSumd02Hist[index]->Fill(sumD02);
+ fSigVertHist[index]->Fill(sigvert);
fPtMaxHist[index]->Fill(ptmax);
fDCAHist[index]->Fill(dca);
if(passTightCuts==1){