2052 double pt1 = recoP1.pt();
2053 double eta1 = recoP1.eta();
2054 double pt2 = recoP2.pt();
2055 double eta2 = recoP2.eta();
2057 double diffPt1 = (pt1 - genP1.pt())/genP1.pt();
2058 double diffPt2 = (pt2 - genP2.pt())/genP2.pt();
2060 double genTheta1 = genP1.theta();
2061 double genTheta2 = genP2.theta();
2062 double recoTheta1 = recoP1.theta();
2063 double recoTheta2 = recoP2.theta();
2065 double genCotgTheta1 = TMath::Cos(genTheta1)/(TMath::Sin(genTheta1));
2066 double genCotgTheta2 = TMath::Cos(genTheta2)/(TMath::Sin(genTheta2));
2067 double recoCotgTheta1 = TMath::Cos(recoTheta1)/(TMath::Sin(recoTheta1));
2068 double recoCotgTheta2 = TMath::Cos(recoTheta2)/(TMath::Sin(recoTheta2));
2072 double diffCotgTheta1 = recoCotgTheta1 - genCotgTheta1;
2073 double diffCotgTheta2 = recoCotgTheta2 - genCotgTheta2;
2081 mapHisto_[name_+
"Pt"]->Fill(pt1, eta1, diffPt1, diffPt1);
2082 mapHisto_[name_+
"Pt"]->Fill(pt2, eta2, diffPt2, diffPt2);
2083 mapHisto_[name_+
"CotgTheta"]->Fill(pt1, eta1, diffCotgTheta1, diffCotgTheta1);
2084 mapHisto_[name_+
"CotgTheta"]->Fill(pt2, eta2, diffCotgTheta2, diffCotgTheta2);
2085 mapHisto_[name_+
"Phi"]->Fill(pt1, eta1, diffPhi1, diffPhi1);
2086 mapHisto_[name_+
"Phi"]->Fill(pt2, eta2, diffPhi2, diffPhi2);
2089 mapHisto_[name_+
"Pt-CotgTheta"]->Fill(pt1, eta1, diffPt1, diffCotgTheta1 );
2090 mapHisto_[name_+
"Pt-CotgTheta"]->Fill(pt2, eta2, diffPt2, diffCotgTheta2 );
2091 mapHisto_[name_+
"Pt-Phi"]->Fill(pt1, eta1, diffPt1, diffPhi1);
2092 mapHisto_[name_+
"Pt-Phi"]->Fill(pt2, eta2, diffPt2, diffPhi2);
2093 mapHisto_[name_+
"CotgTheta-Phi"]->Fill(pt1, eta1, diffCotgTheta1, diffPhi1);
2094 mapHisto_[name_+
"CotgTheta-Phi"]->Fill(pt2, eta2, diffCotgTheta2, diffPhi2);
2099 mapHisto_[name_+
"Pt1-Pt2"]->Fill(pt1, eta1, diffPt1, diffPt2);
2100 mapHisto_[name_+
"Pt1-Pt2"]->Fill(pt2, eta2, diffPt1, diffPt2);
2101 mapHisto_[name_+
"CotgTheta1-CotgTheta2"]->Fill(pt1, eta1, diffCotgTheta1, diffCotgTheta2);
2102 mapHisto_[name_+
"CotgTheta1-CotgTheta2"]->Fill(pt2, eta2, diffCotgTheta1, diffCotgTheta2);
2103 mapHisto_[name_+
"Phi1-Phi2"]->Fill(pt1, eta1, diffPhi1, diffPhi2);
2104 mapHisto_[name_+
"Phi1-Phi2"]->Fill(pt2, eta2, diffPhi1, diffPhi2);
2109 mapHisto_[name_+
"Pt12-CotgTheta21"]->Fill(pt1, eta1, diffPt1, diffCotgTheta2);
2110 mapHisto_[name_+
"Pt12-CotgTheta21"]->Fill(pt2, eta2, diffPt2, diffCotgTheta1);
2111 mapHisto_[name_+
"Pt12-Phi21"]->Fill(pt1, eta1, diffPt1, diffPhi2);
2112 mapHisto_[name_+
"Pt12-Phi21"]->Fill(pt2, eta2, diffPt2, diffPhi1);
2113 mapHisto_[name_+
"CotgTheta12-Phi21"]->Fill(pt1, eta1, diffCotgTheta1, diffPhi2);
2114 mapHisto_[name_+
"CotgTheta12-Phi21"]->Fill(pt2, eta2, diffCotgTheta2, diffPhi1);
static double deltaPhiNoFabs(const double &phi1, const double &phi2)
Without fabs at the end, used to have a symmetric distribution for the resolution fits and variance c...
std::map< TString, HCovarianceVSxy * > mapHisto_