| 51 |
|
h1["trkDzPV"+suffix1] = new TH1D("trkDzPV"+suffix1, "Dz of rec tracks "+suffix1, 300,-50,50); |
| 52 |
|
h1["nrecPV"+suffix1] = new TH1D("nrecPV"+suffix1, "Num of rec pvtx"+suffix1, 50,0,50); |
| 53 |
|
} |
| 54 |
– |
else if (type == "pixVtx") { |
| 55 |
– |
// Book histograms about pixelVertices |
| 56 |
– |
h1["trkdz_pxlpvtxdz"] = new TH1D("trkdz_pxlpvtxdz", "(Track dz - pixelpvtx dz) in cm",300,-0.5,0.5); |
| 57 |
– |
h1["trkdz_pxlpvtxdz_pxlpvtxdzerr"] = new TH1D("trkdz_pxlpvtxdz_pxlpvtxdzerr", "|Track dz - pixelpvtx dz| / pxlpvtxdzErr",300,0,100); |
| 58 |
– |
h1["trkdz_pxlpvtxdz_trkdzerr"] = new TH1D("trkdz_pxlpvtxdz_trkdzerr", "|Track dz - pixelpvtx dz| / trkdzErr",300,0,50); |
| 59 |
– |
h1["trkdzErr_pxlpvtx"] = new TH1D("trkdzErr_pxlpvtxdz", "Track dzErr of leading pixelpvtx ",300,0,0.5); |
| 60 |
– |
h1["trkdzErr_pvtx"] = new TH1D("trkdzErr_pvtx", "Track dzErr of the leading pvtx ",300,0,0.5); |
| 61 |
– |
h1["dzErr_pxlpvtx"] = new TH1D("dzErr_pxlpvtx", "zError of the leading pvtx ",300,0,0.5); |
| 62 |
– |
|
| 63 |
– |
// Compare offlinePrimaryVertices with pixelVertices |
| 64 |
– |
h1["nrecPV_minus_nrecPxlPV"] = new TH1D("nrecPV_minus_nrecPxlPV", "nrecPV_minus_nrecPxlPV",21,-10.5,10.5); |
| 65 |
– |
h1["recxPV_minus_recxPxlPV"] = new TH1D("recxPV_minus_recxPxlPV", "recxPV_minus_recxPxlPV",300,-0.02,0.02); |
| 66 |
– |
h1["recyPV_minus_recyPxlPV"] = new TH1D("recyPV_minus_recyPxlPV", "recyPV_minus_recyPxlPV",300,-0.02,0.02); |
| 67 |
– |
h1["reczPV_minus_reczPxlPV"] = new TH1D("reczPV_minus_reczPxlPV", "reczPV_minus_reczPxlPV",300,-0.1,0.1); |
| 68 |
– |
} |
| 54 |
|
else if (type == "generator") { |
| 55 |
|
h1["genPart_T"] = new TH1D("genPart_T","t component of gen particles",300,-0.5,0.5); |
| 56 |
|
h1["genPart_T"]->GetXaxis()->SetTitle("t (nanosecond)"); |
