- fInvariantMass->Fill(anp4comb.M());
- fInvmassPt->Fill(anp4comb.M(), anp4comb.Vect().Pt());
- if ((mumpdgpos==resonancePDGcode)&&(mumpdgneg==resonancePDGcode)&&(mumlabelpos==mumlabelneg)
- &&(pdgpos==fdaughter1pdg)&&(pdgneg==(-fdaughter2pdg))&&(TMath::Abs(trackpos->GetLabel())>=0)&&(TMath::Abs(trackneg->GetLabel())>=0)&&(mumlabelpos>=0) &&(mumlabelneg>=0)) {
- fOpeningAngle->Fill(openingAngle);
- fInvMassTrue->Fill(p4comb.M());
- fInvmassPtTrue->Fill(anp4comb.M(), anp4comb.Vect().Pt());
- if((TMath::Abs(anp4neg.Vect().Eta())<1.1)&&(TMath::Abs(anp4pos.Vect().Eta())<1.1)&&(anp4comb.Vect().Eta()<1.1)) {
- fRecPt->Fill(anp4comb.Vect().Pt());
- fRecEta->Fill(anp4comb.Vect().Eta());
- fRecEtaPt->Fill(anp4comb.Vect().Pt(), anp4comb.Vect().Eta());
+
+ if(anp4comb.Vect().Pt()<=0.25) continue;
+
+ if((TMath::Abs(anp4neg.Vect().Eta())<0.9)&&(TMath::Abs(anp4pos.Vect().Eta())<0.9)&&(anp4comb.Vect().Eta()<0.9)) {
+ fInvariantMass->Fill(anp4comb.M());
+ fInvmassPt->Fill(anp4comb.M(), anp4comb.Vect().Pt());
+ if ((mumpdgpos==resonancePDGcode)&&(mumpdgneg==resonancePDGcode)&&(mumlabelpos==mumlabelneg)
+ &&(pdgpos==fdaughter1pdg)&&(pdgneg==(-fdaughter2pdg))&&(TMath::Abs(trackpos->GetLabel())>=0)&&(TMath::Abs(trackneg->GetLabel())>=0)&&(mumlabelpos>=0) &&(mumlabelneg>=0)) {
+ fOpeningAngle->Fill(openingAngle);
+ fInvMassTrue->Fill(anp4comb.M());
+ fInvmassPtTrue->Fill(anp4comb.M(), anp4comb.Vect().Pt());
+
+ fRecPt->Fill(anp4comb.Vect().Pt());
+ fRecEta->Fill(anp4comb.Vect().Eta());
+ fRecEtaPt->Fill(anp4comb.Vect().Pt(), anp4comb.Vect().Eta());