ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/Jeng/PVStudy/src/PVHistograms.cc
(Generate patch)

Comparing UserCode/Jeng/PVStudy/src/PVHistograms.cc (file contents):
Revision 1.2 by jengbou, Fri Dec 4 20:41:58 2009 UTC vs.
Revision 1.4 by yygao, Sat Jan 30 18:25:53 2010 UTC

# Line 33 | Line 33 | void PVHistograms::Init(TString type, TS
33    if (type == "pvTrk") {
34      // Histograms of track properties
35      h1["nTrk"+suffix1]     = new TH1D("nTrk"+suffix1, "Num of rec tracks"+suffix1, 300,0,300);
36 <    h1["trkPt"+suffix1]    = new TH1D("trkPt"+suffix1, "Pt of rec tracks "+suffix1, 100,0,100);
36 >    h1["trkPt"+suffix1]    = new TH1D("trkPt"+suffix1, "Pt of rec tracks "+suffix1, 800,0,100);
37      h1["trkEta"+suffix1]   = new TH1D("trkEta"+suffix1, "#eta of rec tracks "+suffix1, 100,-3,3);
38      h1["trkPhi"+suffix1]   = new TH1D("trkPhi"+suffix1, "#phi of rec tracks "+suffix1, 100,-3.2,3.2);
39      h1["trkDxy"+suffix1]   = new TH1D("trkDxy"+suffix1, "Dxy of rec tracks "+suffix1, 100,-0.5,0.5);
40 <    h1["trkDz"+suffix1]    = new TH1D("trkDz"+suffix1, "Dz of rec tracks "+suffix1, 100,-50,50);
40 >    h1["trkDxyCorr"+suffix1] = new TH1D("trkDxyCorr"+suffix1, TString("BS Corrected Dxy of rec tracks "+suffix1),100,-0.5,0.5);
41 >    h1["trkDz"+suffix1]    = new TH1D("trkDz"+suffix1, "Dz of rec tracks "+suffix1, 300,-50,50);
42  
43      h1["nTrkPV"+suffix1]   = new TH1D("nTrkPV"+suffix1, "Num of rec tracks in PV"+suffix1, 300,0,300);
44 <    h1["trkPtPV"+suffix1]  = new TH1D("trkPtPV"+suffix1, "Pt of rec tracks in "+suffix1 ,100,0,100);
44 >    h1["nHWTrkPV"+suffix1]   = new TH1D("nHWTrkPV"+suffix1, "Num of rec tracks with over 0.5 weight in PV"+suffix1, 300,0,300);
45 >    h1["ndofPV"+suffix1]   = new TH1D("ndofPV"+suffix1, "PV ndof"+suffix1, 300,0,100);
46 >    h1["trkPtPV"+suffix1]  = new TH1D("trkPtPV"+suffix1, "Pt of rec tracks in "+suffix1 ,800,0,100);
47      h1["trkEtaPV"+suffix1] = new TH1D("trkEtaPV"+suffix1, "#eta of rec tracks in PV"+suffix1, 100,-3,3);
48      h1["trkPhiPV"+suffix1] = new TH1D("trkPhiPV"+suffix1, "#phi of rec tracks in PV"+suffix1, 100,-3.2,3.2);
49 <    h1["trkDxyPV"+suffix1] = new TH1D("trkDxyPV"+suffix1, "Dxy of rec tracks in PV"+suffix1, 100,-5,5);
50 <    h1["trkDzPV"+suffix1]  = new TH1D("trkDzPV"+suffix1, "Dz of rec tracks "+suffix1, 100,-50,50);
49 >    h1["trkDxyPV"+suffix1] = new TH1D("trkDxyPV"+suffix1, "Dxy of rec tracks in PV"+suffix1, 100,-0.5,0.5);
50 >    h1["trkDxyCorrPV"+suffix1] = new TH1D("trkDxyCorrPV"+suffix1, "BS corrected Dxy of rec tracks in PV"+suffix1, 100,-0.5,0.5);
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    }
50 –  else if (type == "pixVtx") {
51 –    // Book histograms about pixelVertices
52 –    h1["trkdz_pxlpvtxdz"]              = new TH1D("trkdz_pxlpvtxdz", "(Track dz - pixelpvtx dz) in cm",300,-0.5,0.5);
53 –    h1["trkdz_pxlpvtxdz_pxlpvtxdzerr"] = new TH1D("trkdz_pxlpvtxdz_pxlpvtxdzerr", "|Track dz - pixelpvtx dz| / pxlpvtxdzErr",300,0,100);  
54 –    h1["trkdz_pxlpvtxdz_trkdzerr"]     = new TH1D("trkdz_pxlpvtxdz_trkdzerr", "|Track dz - pixelpvtx dz| / trkdzErr",300,0,50);
55 –    h1["trkdzErr_pxlpvtx"]             = new TH1D("trkdzErr_pxlpvtxdz", "Track dzErr of leading pixelpvtx ",300,0,0.5);
56 –    h1["trkdzErr_pvtx"]                = new TH1D("trkdzErr_pvtx", "Track dzErr of the leading pvtx ",300,0,0.5);
57 –    h1["dzErr_pxlpvtx"]                = new TH1D("dzErr_pxlpvtx", "zError of the leading pvtx ",300,0,0.5);
58 –
59 –    // Compare offlinePrimaryVertices with pixelVertices
60 –    h1["nrecPV_minus_nrecPxlPV"]       = new TH1D("nrecPV_minus_nrecPxlPV", "nrecPV_minus_nrecPxlPV",21,-10.5,10.5);
61 –    h1["recxPV_minus_recxPxlPV"]       = new TH1D("recxPV_minus_recxPxlPV", "recxPV_minus_recxPxlPV",300,-0.02,0.02);
62 –    h1["recyPV_minus_recyPxlPV"]       = new TH1D("recyPV_minus_recyPxlPV", "recyPV_minus_recyPxlPV",300,-0.02,0.02);
63 –    h1["reczPV_minus_reczPxlPV"]       = new TH1D("reczPV_minus_reczPxlPV", "reczPV_minus_reczPxlPV",300,-0.1,0.1);
64 –  }
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)");
# Line 72 | Line 61 | void PVHistograms::Init(TString type, TS
61    else if (type == "misc") {
62      h1["nrecPVDiff"]     = new TH1D("nrecPVDiff","nrecPV1-nRecPV2",21,-10.5,10.5);
63      h1["nTrkPVDiff"]     = new TH1D("nTrkPVDiff","nTrkPV1-nTrkPV2",41,-20.5,20.5);
64 <    h1["nTrkPVRelDiff"]  = new TH1D("nTrkPVRelDiff","(nTrkPV1-nTrkPV2)/(nTrkPV1+nTrkPV2)",100,-1,1);
64 >    h1["nTrkPVRelDiff"]  = new TH1D("nTrkPVRelDiff","(nTrkPV1-nTrkPV2)/(nTrkPV1+nTrkPV2)",100,-1,1);  
65 >    h1["ndofPVDiff"]     = new TH1D("ndofPVDiff","ndofPV1-ndofPV2",100,-10,10);
66 >    h1["ndofPVRelDiff"]  = new TH1D("ndofPVRelDiff","(ndofPV1-ndofPV2)/(ndofPV1+ndofPV2)",100,-1,1);  
67 >    h1["twovtxzsign"]     = new TH1D("twovtxzsign", "two vtx z signficance",300,-20,20);
68      // Histograms on comparing the multi-vertices
69      // Difference in reconstructed vtx position
70      h1["min_xsep"]       = new TH1D("min_xsep", "min x diff of primary and secondary pvtx",300,0,0.1);
# Line 84 | Line 76 | void PVHistograms::Init(TString type, TS
76      // Difference in reconstructed vtx position
77      h1["min_ntrksep"]    = new TH1D("min_ntrksep", "min nTrk diff of primary and secondary pvtx",201,-50.5,150.5);
78      h1["min_sumpt2sep"]  = new TH1D("min_sumpt2sep", "min sumpt2 diff of primary and secondary pvtx",300,0,10000);
79 +
80    }
81    else if (type == "summary") {
82 <    h1["deltax"+suffix1] = new TH1D("deltax"+suffix1, "x-residual pvtx"+suffix1, 800,-0.04,0.04);
83 <    h1["deltay"+suffix1] = new TH1D("deltay"+suffix1, "y-residual pvtx"+suffix1, 800,-0.04,0.04);
84 <    h1["deltaz"+suffix1] = new TH1D("deltaz"+suffix1, "z-residual pvtx"+suffix1, 800,-0.04,0.04);
85 <    h1["pullx"+suffix1]  = new TH1D("pullx"+suffix1, "x-pull pvtx"+suffix1, 800,-10,10);
86 <    h1["pully"+suffix1]  = new TH1D("pully"+suffix1, "y-pull pvtx"+suffix1, 800,-10,10);
87 <    h1["pullz"+suffix1]  = new TH1D("pullz"+suffix1, "z-pull pvtx"+suffix1, 800,-10,10);
88 <    h1["errPVx"+suffix1] = new TH1D("errPVx"+suffix1, "X"+suffix1+" vertex error", 100,0.,0.02);
89 <    h1["errPVy"+suffix1] = new TH1D("errPVy"+suffix1, "Y"+suffix1+" vertex error", 100,0.,0.02);
90 <    h1["errPVz"+suffix1] = new TH1D("errPVz"+suffix1, "Z"+suffix1+" vertex error", 100,0.,0.02);
82 >    h1["deltax"+suffix1] = new TH1D("deltax"+suffix1, "x-residual pvtx"+suffix1, 800,-0.1,0.1);
83 >    h1["deltay"+suffix1] = new TH1D("deltay"+suffix1, "y-residual pvtx"+suffix1, 800,-0.1,0.1);
84 >    h1["deltaz"+suffix1] = new TH1D("deltaz"+suffix1, "z-residual pvtx"+suffix1, 800,-0.1,0.1);
85 >    h1["pullx"+suffix1]  = new TH1D("pullx"+suffix1, "x-pull pvtx"+suffix1, 800,-5,5);
86 >    h1["pully"+suffix1]  = new TH1D("pully"+suffix1, "y-pull pvtx"+suffix1, 800,-5,5);
87 >    h1["pullz"+suffix1]  = new TH1D("pullz"+suffix1, "z-pull pvtx"+suffix1, 800,-5,5);
88 >    h1["errPVx"+suffix1] = new TH1D("errPVx"+suffix1, "X"+suffix1+" vertex error", 200,0.,0.1);
89 >    h1["errPVy"+suffix1] = new TH1D("errPVy"+suffix1, "Y"+suffix1+" vertex error", 200,0.,0.1);
90 >    h1["errPVz"+suffix1] = new TH1D("errPVz"+suffix1, "Z"+suffix1+" vertex error", 200,0.,0.1);
91    }
92    else if (type == "others") {
93 <    h1["deltax"+suffix1] = new TH1D("deltax"+suffix1, "x-residual of pvtx"+methodName, 100,-0.02,0.02);
93 >    h1["deltax"+suffix1] = new TH1D("deltax"+suffix1, "x-residual of pvtx"+methodName, 200,-0.1,0.1);
94      h1["deltax"+suffix1]->GetXaxis()->SetTitle("cm");
95 <    h1["deltay"+suffix1] = new TH1D("deltay"+suffix1, "y-residual of pvtx"+methodName, 100,-0.02,0.02);
95 >    h1["deltay"+suffix1] = new TH1D("deltay"+suffix1, "y-residual of pvtx"+methodName, 200,-0.1,0.1);
96      h1["deltay"+suffix1]->GetXaxis()->SetTitle("cm");
97 <    h1["deltaz"+suffix1] = new TH1D("deltaz"+suffix1, "z-residual of pvtx"+methodName, 100,-0.02,0.02);
97 >    h1["deltaz"+suffix1] = new TH1D("deltaz"+suffix1, "z-residual of pvtx"+methodName, 200,-0.1,0.1);
98      h1["deltaz"+suffix1]->GetXaxis()->SetTitle("cm");
99 <    h1["pullx"+suffix1]  = new TH1D("pullx"+suffix1, "x-pull of pvtx"+methodName, 100,-10.,10.);
100 <    h1["pully"+suffix1]  = new TH1D("pully"+suffix1, "y-pull of pvtx"+methodName, 100,-10.,10.);
101 <    h1["pullz"+suffix1]  = new TH1D("pullz"+suffix1, "z-pull of pvtx"+methodName, 100,-10.,10.);
102 <    h1["errPVx"+suffix1] = new TH1D("errPVx"+suffix1, "X"+methodName+" vertex error", 100,0.,0.02);
103 <    h1["errPVy"+suffix1] = new TH1D("errPVy"+suffix1, "Y"+methodName+" vertex error", 100,0.,0.02);
104 <    h1["errPVz"+suffix1] = new TH1D("errPVz"+suffix1, "Z"+methodName+" vertex error", 100,0.,0.02);
99 >    h1["pullx"+suffix1]  = new TH1D("pullx"+suffix1, "x-pull of pvtx"+methodName, 200,-5.,5.);
100 >    h1["pully"+suffix1]  = new TH1D("pully"+suffix1, "y-pull of pvtx"+methodName, 200,-5.,5.);
101 >    h1["pullz"+suffix1]  = new TH1D("pullz"+suffix1, "z-pull of pvtx"+methodName, 200,-5.,5.);
102 >    h1["errPVx"+suffix1] = new TH1D("errPVx"+suffix1, "X"+methodName+" vertex error", 200,0.,0.1);
103 >    h1["errPVy"+suffix1] = new TH1D("errPVy"+suffix1, "Y"+methodName+" vertex error", 200,0.,0.1);
104 >    h1["errPVz"+suffix1] = new TH1D("errPVz"+suffix1, "Z"+methodName+" vertex error", 200,0.,0.1);
105    }
106  
107   }
# Line 119 | Line 112 | void PVHistograms::InitA(TString type, T
112    if (suffix2 != "") suffix1 += "_" + suffix2;
113  
114    if (type == "analysis") {
115 <    h2["resx"+suffix1]= new TH2D("resx"+suffix1, "x-resolution vs number of tracks in pvtx"+methodName, nTrkMax_-nTrkMin_,nTrkMin_,nTrkMax_,400,0.,200);
115 >    h2["resx"+suffix1]= new TH2D("resx"+suffix1, "x-resolution vs number of tracks in pvtx"+methodName, nTrkMax_-nTrkMin_,nTrkMin_,nTrkMax_,400,0.,500);
116      h2["resx"+suffix1]->SetMarkerStyle(21);
117      h2["resx"+suffix1]->SetMarkerColor(4);
118      h2["resx"+suffix1]->GetXaxis()->SetTitle("Num of tracks");
119      h2["resx"+suffix1]->GetYaxis()->SetTitle("#mum");
120 <    h2["resy"+suffix1]= new TH2D("resy"+suffix1, "y-resolution vs number of tracks in pvtx"+methodName, nTrkMax_-nTrkMin_,nTrkMin_,nTrkMax_,400,0.,200);
120 >    h2["resy"+suffix1]= new TH2D("resy"+suffix1, "y-resolution vs number of tracks in pvtx"+methodName, nTrkMax_-nTrkMin_,nTrkMin_,nTrkMax_,400,0.,500);
121      h2["resy"+suffix1]->SetMarkerStyle(21);
122      h2["resy"+suffix1]->SetMarkerColor(4);
123      h2["resy"+suffix1]->GetXaxis()->SetTitle("Num of tracks");
124      h2["resy"+suffix1]->GetYaxis()->SetTitle("#mum");
125 <    h2["resz"+suffix1]= new TH2D("resz"+suffix1, "z-resolution vs number of tracks in pvtx"+methodName, nTrkMax_-nTrkMin_,nTrkMin_,nTrkMax_,400,0.,200);
125 >    h2["resz"+suffix1]= new TH2D("resz"+suffix1, "z-resolution vs number of tracks in pvtx"+methodName, nTrkMax_-nTrkMin_,nTrkMin_,nTrkMax_,400,0.,500);
126      h2["resz"+suffix1]->SetMarkerStyle(21);
127      h2["resz"+suffix1]->SetMarkerColor(4);
128      h2["resz"+suffix1]->GetXaxis()->SetTitle("Num of tracks");

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines