loon: AlignmentHistograms Class Reference

Constructor & Destructor Documentation

AlignmentHistograms::AlignmentHistograms	(	TFile *	rootfile,
		const VldContext &	vld
	)

Definition at line 48 of file AlignmentHistograms.cxx.

References f2hChargevsResidualHitStrips, f2hNumberOfTrackHitsPerModule, f2hNumberOfTrackHitsPerStrip, f2hPlanePlexChargeRatio, f2hPlanePlexHits, f2hResidualvsPlane, fChargeLimit, fChargesInAllHitStrip, fChargevsResidualHitStrips, fNChargeBins, fnplanes, fNResidualBins, fnstrips, fPlanePlexDirectory, fpMeanStripResponsevsResidual, fpResidualvsCosZUview, fpResidualvsCosZVview, fpResidualvsPlane, fpResidualvsStripChargeUview, fpResidualvsStripChargeVview, fpResidualvsTPosSigmaUview, fpResidualvsTPosSigmaVview, fpResidualvsTrackChargeRatioUview, fpResidualvsTrackChargeRatioVview, fpResidualvsTrackChargeUview, fpResidualvsTrackChargeVview, fpResidualvsTrackNStripsUview, fpResidualvsTrackNStripsVview, fResidualAllStrips, fResidualLimit, fResidualSummaryDirectory, fResidualUStrips, fResidualVStrips, fRootFile, fTrackDirectory, UgliGeomHandle::GetScintPlnHandleVector(), MSG, and SetAxisLabels().

00049    :fRootFile(rootfile),
00050     fnplanes(0),
00051     fnstrips(0),
00052     fNResidualBins(400),
00053     fResidualLimit(0.2),
00054     fResidualStripBinWidth(0.001),
00055     fResidualMdlBinWidth(0.001),
00056     fNResidualBinsStrip(60),
00057     fResidualLimitStrip(0.12),
00058     fNChargeBins(1000),
00059     fChargeLimit(10000.0),
00060     fNTrackChargeBins(2000),
00061     fTrackChargeLimit(100000.0),
00062     fNHitsBins(100),
00063     fTrackResidualRMSBins(1000),
00064     fTrackResidualRMSLimit(0.2),
00065     fNCosBins(600),
00066     //Histograms from here down
00067     fhCharge2dTracksU(0),
00068     fhCharge2dTracksV(0),
00069     fhCharge2dTracksDiffNorm(0),
00070     fhNStrip2dTracksU(0),
00071     fhNStrip2dTracksV(0),
00072     fhNStrip2dTracksDiffNorm(0),
00073     fhSigmaOfTPosVview(0),
00074     fhSigmaOfTPosUview(0),
00075     fhSigmaOfTPosDiffNorm(0),
00076     fChargesInAllHitStrip(0),
00077     fChargevsResidualHitStrips(0),  
00078     fResidualAllStrips(0), 
00079     fResidualVStrips(0), 
00080     fResidualUStrips(0), 
00081     fTrackResidualRMSVview(0),
00082     fTrackResidualRMSUview(0),
00083     fChargeRatioOf2dTracks(0),
00084     fhTrackHoughCosU(0),
00085     fhTrackHoughCosV(0),
00086     fhTrackHoughCosZ(0),
00087     fhTrackChargeOverTotalChargeRatio(0),
00088     //2D histograms from here down        
00089     f2hNumberOfTrackHitsPerStrip(0), 
00090     f2hNumberOfTrackHitsPerModule(0), 
00091     f2hChargevsResidualHitStrips(0),  
00092     f2hResidualvsPlane(0),
00093     f2hMeanChargePlaneStrip(0),
00094     f2hMeanPathCorrChargePlaneStrip(0),
00095     //Track TProfile histograms from here down
00096     fpTPosSigmavsCosZVview(0),
00097     fpTPosSigmavsCosZUview(0),
00098     fpTPosSigmavsTrackChargeVview(0),
00099     fpTPosSigmavsTrackChargeUview(0),
00100     fpTPosSigmavsTrackNStripsVview(0),
00101     fpTPosSigmavsTrackNStripsUview(0),
00102     fpTPosSigmavsTrackChargeRatioVview(0),
00103     fpTPosSigmavsTrackChargeRatioUview(0),
00104     fp3dTrackChargevsTrackChargeRatio(0),
00105     fp3dTrackChargevsTrackCosz(0),
00106     fp3dTrackNStripsvsCosz(0),
00107     fp3dTrackChargevsNStrips(0),
00108     fp3dTrackChargeRatiovsCosz(0),
00109     //Strip TProfile histograms
00110     fpResidualvsPlane(0),
00111     fpResidualvsCosZVview(0),
00112     fpResidualvsCosZUview(0),
00113     fpResidualvsTPosSigmaVview(0),
00114     fpResidualvsTPosSigmaUview(0),
00115     fpResidualvsStripChargeVview(0),
00116     fpResidualvsStripChargeUview(0),
00117     fpResidualvsTrackChargeVview(0),
00118     fpResidualvsTrackChargeUview(0),
00119     fpResidualvsTrackNStripsVview(0),
00120     fpResidualvsTrackNStripsUview(0),
00121     fpResidualvsTrackChargeRatioVview(0),
00122     fpResidualvsTrackChargeRatioUview(0),
00123     fpMeanStripResponsevsResidual(0)
00124 {
00125 
00126    MSG("Align", Msg::kInfo) << "Constructor AlignmentHistograms() " << endl;
00127 
00128    if(!fRootFile)
00129    {
00130       MSG("Align", Msg::kFatal) << "ROOT file pointer is zero. Abort()" <<endl;
00131       abort();
00132    }
00133 
00134    fRootFile -> cd();   
00135 
00136    fPlanePlexDirectory       = new TDirectory("PlanePlex","PlanePlex");
00137    fResidualSummaryDirectory = new TDirectory("ResidualSummary","ResidualSummary");   
00138    fTrackDirectory           = new TDirectory("TrackData","TrackData");   
00139 
00140    UgliGeomHandle ugh(vld);
00141    vector<UgliScintPlnHandle> uph  = ugh.GetScintPlnHandleVector();
00142    for(vector<UgliScintPlnHandle>::iterator it = uph.begin(); it != uph.end(); ++it)
00143    {
00144       const int plane = it -> GetPlaneNumber();
00145       if(plane > fnplanes)
00146          fnplanes = plane;
00147 
00148       vector<UgliStripHandle> ushv = it -> GetStripHandleVector();
00149       for(vector<UgliStripHandle>::iterator sit = ushv.begin(); sit != ushv.end(); ++sit)
00150       {
00151          const int strip = (sit -> GetSEId()).GetStrip();
00152          if(strip > fnstrips)
00153             fnstrips = strip;
00154       }
00155       
00156       const int nstrip = it -> NumberOfStrips();
00157       char name[100], title[100];
00158       
00159       sprintf(name,"PlexChargeRatioPlane%03d", plane);
00160       sprintf(title,"cand strip charge over trackv strip charge plane %03d", plane);
00161       TH2D *h_ratio = new TH2D(name, title, nstrip, 0, nstrip, nstrip, 0, nstrip);
00162       h_ratio -> SetDirectory(fPlanePlexDirectory);
00163       SetAxisLabels(h_ratio, "track strip #", "candidate strip #");
00164       f2hPlanePlexChargeRatio[plane] = h_ratio;
00165 
00166       sprintf(name,"PlexHitsPlane%03d", plane);
00167       sprintf(title,"Number of entries in each bin plane %03d", plane);
00168       TH2D *h_hits = new TH2D(name, title, nstrip, 0, nstrip, nstrip, 0, nstrip);
00169       h_hits -> SetDirectory(fPlanePlexDirectory);
00170       SetAxisLabels(h_hits, "track strip #", "candidate strip #");
00171       f2hPlanePlexHits[plane] = h_hits; 
00172    }
00173    
00174    fChargesInAllHitStrip = new TH1D("ChargesInAllHitStrip",
00175                                     "Charge in all hit strips",
00176                                     fNChargeBins, 0, fChargeLimit); 
00177    SetAxisLabels(fChargesInAllHitStrip, "Raw ADC", 0);
00178    
00179 
00180    fChargevsResidualHitStrips = new TH1D("ChargevsResidualHitStrips2h", "Charge of hit strip vs fited residual",
00181                                          2*fNResidualBins , -fResidualLimit, fResidualLimit);
00182    SetAxisLabels(fChargevsResidualHitStrips, "residual (m)", "Raw ADC");
00183    
00184    f2hChargevsResidualHitStrips = new TH2D("ChargevsResidualHitStrips", "Charge of hit strip vs fited residual",
00185                                            2*fNResidualBins , -fResidualLimit, fResidualLimit,
00186                                            fNChargeBins, 0, fChargeLimit);
00187    SetAxisLabels(f2hChargevsResidualHitStrips, "residual (m)", "Raw ADC");
00188    
00189    fResidualSummaryDirectory -> cd();
00190    
00191    fResidualAllStrips = new TH1D("ResidualAllStrips", "Residuals of all strips",
00192                                  2*fNResidualBins , -fResidualLimit, fResidualLimit);   
00193    SetAxisLabels(fResidualAllStrips,"residual (m)",0);
00194    
00195    fResidualUStrips = new TH1D("ResidualUStrips", "Residuals of U strips",
00196                                2*fNResidualBins , -fResidualLimit, fResidualLimit);   
00197    SetAxisLabels(fResidualUStrips, "residual (m)",0);
00198    
00199    fResidualVStrips = new TH1D("ResidualVStrips",
00200                                "Residuals of V strips",
00201                                2*fNResidualBins , -fResidualLimit, fResidualLimit);   
00202    fResidualVStrips -> GetXaxis() -> SetTitle("residual (m)");
00203    fResidualVStrips -> GetXaxis() -> CenterTitle();
00204    
00205    f2hResidualvsPlane = new TH2D("ResidualvsPlane",
00206                                  "Residual vs plane number, all events",
00207                                  fnplanes, 0, fnplanes,
00208                                  2*fNResidualBins , -fResidualLimit, fResidualLimit); 
00209    SetAxisLabels(f2hResidualvsPlane, "residual (m)","plane number");
00210    
00211    fRootFile -> cd();   
00212    
00213    f2hNumberOfTrackHitsPerStrip = new TH2D("NumberOfTrackHitsPerStrip", "Strip vs plane for track hits",
00214                                            fnplanes, 0, fnplanes, fnstrips, 0, fnstrips);
00215    SetAxisLabels(f2hNumberOfTrackHitsPerStrip,"plane number", "strip number");
00216    
00217    
00218    f2hNumberOfTrackHitsPerModule = new TH2D("NumberOfTrackHitsPerModule",
00219                                           "Module vs plane for track hits",
00220                                           fnplanes, 0, fnplanes,
00221                                           5, 0, 5);
00222    SetAxisLabels(f2hNumberOfTrackHitsPerModule, "plane number", "module number");
00223 
00224    fResidualSummaryDirectory -> cd();   
00225    
00226    //____________________________________________________________________________________________
00227    fpMeanStripResponsevsResidual = new TProfile("MeanStripResponsevsResidual", 
00228                                                 "Response of scintilator vs residual, path length corrected",
00229                                                 2*fNResidualBins , -fResidualLimit, fResidualLimit);   
00230    SetAxisLabels(fpMeanStripResponsevsResidual, "residual (m)", "charge (ADC)");
00231 
00232    //_______________________________________________________________________________
00233    fpResidualvsPlane = new TProfile("ProfileResidualvsPlane", "Residual vs plane number", 
00234                                     fnplanes, 0.0, fnplanes, "s");
00235    SetAxisLabels(fpResidualvsPlane, "residual (m)");
00236    
00237    //_______________________________________________________________________________
00238    fpResidualvsCosZVview = new TProfile("ResidualvsCosZVview", 
00239                                         "Residual vs cos with Z axis V planes", 
00240                                         100, -1.0, 1.0);
00241    SetAxisLabels(fpResidualvsCosZVview, "residual (m)");
00242    
00243    fpResidualvsCosZUview = new TProfile("ResidualvsCosZUview", 
00244                                         "Residual vs cos with Z axis U plane", 
00245                                         100, -1.0, 1.0);
00246    SetAxisLabels(fpResidualvsCosZUview, "residual (m)");
00247    
00248    //____________________________________________________________________________________________
00249    fpResidualvsTPosSigmaVview = new TProfile("ResidualvsTPosSigmaVview", 
00250                                              "Residual vs tpos error V planes", 
00251                                              100, 0, 0.05);
00252    SetAxisLabels(fpResidualvsTPosSigmaVview, "residual (m)","error in tpos (m)");
00253 
00254    fpResidualvsTPosSigmaUview = new TProfile("ResidualvsTPosSigmaUview", 
00255                                              "Residual vs tpos error U planes", 
00256                                              100, 0, 0.05);
00257    SetAxisLabels(fpResidualvsTPosSigmaUview, "residual (m)", "error in tpos (m)");
00258 
00259    //____________________________________________________________________________________________
00260    fpResidualvsStripChargeVview = new TProfile("ResidualvsStripChargeVview", 
00261                                                "Residual vs strip charge V planes", 
00262                                                200, 0, 10000);
00263    SetAxisLabels(fpResidualvsStripChargeVview, "residual (m)", "strip charge (ADC)");
00264 
00265    fpResidualvsStripChargeUview = new TProfile("ResidualvsStripChargeUview", 
00266                                                "Residual vs strip charge U planes", 
00267                                                200, 0, 10000);
00268    SetAxisLabels(fpResidualvsStripChargeUview, "residual (m)", "strip charge (ADC)");
00269 
00270    //____________________________________________________________________________________________
00271    fpResidualvsTrackChargeVview = new TProfile("ResidualvsTrackChargeVview", 
00272                                                "Residual vs track charge V planes", 
00273                                                100, 0, 60000);
00274    SetAxisLabels(fpResidualvsTrackChargeVview, "residual (m)", "track charge (ADC)");
00275    
00276    fpResidualvsTrackChargeUview = new TProfile("ResidualvsTrackChargeUview", 
00277                                                "Residual vs track charge U planes", 
00278                                                100, 0, 60000);
00279    SetAxisLabels(fpResidualvsTrackChargeUview, "residual (m)", "track charge (ADC)");
00280 
00281    //____________________________________________________________________________________________
00282    fpResidualvsTrackNStripsVview = new TProfile("ResidualvsTrackNStripsVview", 
00283                                                 "Residual vs # of strips in track V planes", 
00284                                                 100, 0.0, 100.0);
00285    SetAxisLabels(fpResidualvsTrackNStripsVview, "residual (m)", "#strips");
00286    
00287    fpResidualvsTrackNStripsUview = new TProfile("ResidualvsTrackNStripsUview", 
00288                                                 "Residual vs # of strips in track U planes", 
00289                                                 100, 0.0, 100.0);
00290    SetAxisLabels(fpResidualvsTrackNStripsUview, "residual (m)", "#strips");
00291 
00292    //____________________________________________________________________________________________
00293    fpResidualvsTrackChargeRatioVview = new TProfile("ResidualvsTrackChargeRatioVview", 
00294                                                     "Error in tpos vs # track/cand charge V planes", 
00295                                                     100, 0.0, 1.0);
00296    SetAxisLabels(fpResidualvsTrackChargeRatioVview, "error in tpos (m)");
00297 
00298    fpResidualvsTrackChargeRatioUview = new TProfile("ResidualvsTrackChargeRatioUview", 
00299                                                     "Error in tpos vs # track/cand charge U planes", 
00300                                                     100, 0.0, 1.0);
00301    SetAxisLabels(fpResidualvsTrackChargeRatioUview, "error in tpos (m)");
00302 
00303    fRootFile -> cd();
00304 
00305 }

AlignmentHistograms::~AlignmentHistograms ( ) [virtual]

Definition at line 615 of file AlignmentHistograms.cxx.

References Fill(), fNMdlHits, fNStripHits, fRootFile, PlexScintMdlId::GetEncoded(), PlexStripEndId::GetEncoded(), PlexScintMdlId::GetModule(), PlexPlaneId::GetPlane(), PlexStripEndId::GetStrip(), and MSG.

00616 {
00617    MSG("Align", Msg::kInfo) << "Destructor ~AlignmentHistograms()" << endl;
00618 
00619    TTree *stree = new TTree("StripHits","Number of track hits per strip");
00620    stree -> SetDirectory(fRootFile);
00621    Int_t plane, strip, module;
00622    UInt_t plexseid, plexmdlid, nhits;
00623    stree -> Branch("nhits", &nhits, "nhits/i");
00624    stree -> Branch("plane", &plane, "plane/I");
00625    stree -> Branch("strip", &strip, "strip/I");
00626    stree -> Branch("plexseid", &plexseid, "plexseid/i");
00627 
00628    for(map<PlexStripEndId, unsigned int>::const_iterator sit = fNStripHits.begin();
00629        sit != fNStripHits.end(); ++sit)
00630    {
00631       PlexStripEndId plexid = sit -> first;
00632       plane    = plexid.GetPlane();
00633       strip    = plexid.GetStrip();
00634       plexseid = plexid.GetEncoded();
00635       nhits    = sit -> second;
00636       stree -> Fill();
00637    }
00638 
00639    TTree *mtree = new TTree("MdlHits","Number of track hits per module");
00640    stree -> SetDirectory(fRootFile);
00641    mtree -> Branch("nhits", &nhits, "nhits/i");
00642    mtree -> Branch("plane", &plane, "plane/I");
00643    mtree -> Branch("module", &module, "module/I");
00644    mtree -> Branch("plexmdlid", &plexmdlid, "plexmdlid/i"); 
00645 
00646    for(map<PlexScintMdlId, unsigned int>::const_iterator mit = fNMdlHits.begin(); 
00647        mit != fNMdlHits.end(); ++mit)
00648    {
00649       PlexScintMdlId plexid = mit -> first;
00650       plane     = plexid.GetPlane();
00651       module    = plexid.GetModule();
00652       plexmdlid = plexid.GetEncoded();
00653       nhits     = mit -> second;
00654       mtree -> Fill();
00655    }
00656    MSG("Align", Msg::kInfo) << "Destructor ~AlignmentHistograms().. Done" << endl;
00657 }

00310 { 00311 PlexStripEndId plexid(astrip.plexseid); 00312 00313 fChargesInAllHitStrip -> Fill(astrip.charge); 00314 f2hNumberOfTrackHitsPerStrip -> Fill(plexid.GetPlane(), plexid.GetStrip(), 1.0); 00315 f2hNumberOfTrackHitsPerModule -> Fill(plexid.GetPlane(), plexid.GetScintMdlId().GetModule(), 1.0); 00316 fChargevsResidualHitStrips -> Fill(astrip.residual, astrip.charge); 00317 f2hChargevsResidualHitStrips -> Fill(astrip.residual, astrip.charge, 1.0); 00318 fResidualAllStrips -> Fill(astrip.residual); 00319 f2hResidualvsPlane -> Fill(plexid.GetPlane(), astrip.residual); 00320 fpResidualvsPlane -> Fill(plexid.GetPlane(), astrip.residual); 00321 fpMeanStripResponsevsResidual -> Fill(astrip.residual, astrip.charge*fabs(ntprec->hcosz), 1.0); 00322 00323 double sum = ntprec->vcharge + ntprec->ucharge + ntprec->vcandcharge + ntprec->ucandcharge; 00324 double ratio = -1.0; 00325 if(sum>1.0) 00326 ratio = (ntprec->vcharge + ntprec->ucharge)/sum; 00327 00328 if(plexid.GetPlaneView() == PlaneView::kV) 00329 { 00330 fResidualVStrips -> Fill(astrip.residual); 00331 fpResidualvsStripChargeVview -> Fill(astrip.charge, astrip.residual); 00332 fpResidualvsTPosSigmaVview -> Fill(ntprec->vsigmaoftpos, astrip.residual); 00333 fpResidualvsCosZVview -> Fill(ntprec->hcosz, astrip.residual); 00334 fpResidualvsTrackChargeVview -> Fill(ntprec->vcharge, astrip.residual); 00335 fpResidualvsTrackNStripsVview -> Fill(ntprec->ntrackvstrip, astrip.residual); 00336 fpResidualvsTrackChargeRatioVview -> Fill(ratio, astrip.residual); 00337 } else 00338 { 00339 fResidualUStrips -> Fill(astrip.residual); 00340 fpResidualvsStripChargeUview -> Fill(astrip.charge, astrip.residual); 00341 fpResidualvsTPosSigmaUview -> Fill(ntprec->usigmaoftpos, astrip.residual); 00342 fpResidualvsCosZUview -> Fill(ntprec->hcosz, astrip.residual); 00343 fpResidualvsTrackChargeUview -> Fill(ntprec->ucharge, astrip.residual); 00344 fpResidualvsTrackNStripsUview -> Fill(ntprec->ntrackustrip, astrip.residual); 00345 fpResidualvsTrackChargeRatioUview -> Fill(ratio, astrip.residual); 00346 } 00347 00348 FillTrackInfo(ntprec); 00349 FillHits(plexid); 00350 }

00355 { 00356 map<PlexStripEndId, unsigned int>::iterator sit = fNStripHits.find(plexid); 00357 if(sit == fNStripHits.end()) 00358 fNStripHits[plexid] = 1; 00359 else 00360 sit->second += 1; 00361 00362 PlexScintMdlId plexmdlid = plexid.GetScintMdlId(); 00363 map<PlexScintMdlId, unsigned int>::iterator mit = fNMdlHits.find(plexmdlid); 00364 if(mit == fNMdlHits.end()) 00365 fNMdlHits[plexmdlid] = 1; 00366 else 00367 mit->second += 1; 00368 }

00565 { 00566 00567 //for(unsigned int i = 0; i < fTrackVStrip.size(); ++i) 00568 for(vector<AlignmentStrip>::const_iterator it = track_strips.begin(); 00569 it != track_strips.end(); ++it) 00570 { 00571 const AlignmentStrip &track_strip = *it; 00572 00573 if(track_strip.charge < 0.1) continue; 00574 00575 for(vector<AlignmentStrip>::const_iterator cit = cand_strips.begin(); 00576 cit != cand_strips.end(); ++cit) 00577 { 00578 const AlignmentStrip &cand_strip = *cit; 00579 if(track_strip.plane == cand_strip.plane) 00580 { 00581 const double ratio = cand_strip.charge/track_strip.charge; 00582 map<int, TH2D *>::iterator fit_r = f2hPlanePlexChargeRatio.find(track_strip.plane); 00583 if(fit_r != f2hPlanePlexChargeRatio.end()) 00584 fit_r -> second -> Fill(track_strip.strip, cand_strip.strip, ratio); 00585 map<int, TH2D *>::iterator fit_h = f2hPlanePlexHits.find(track_strip.plane); 00586 if(fit_h != f2hPlanePlexHits.end()) 00587 fit_h -> second -> Fill(track_strip.strip, cand_strip.strip, ratio); 00588 } 00589 } 00590 } 00591 }

00373 { 00374 if(!fhCharge2dTracksV){ 00375 fTrackDirectory -> cd(); 00376 00377 fp3dTrackChargevsTrackChargeRatio = new TProfile("TrackChargevsTrackChargeRatio", 00378 "3d track charge vs track charge/record charge", 00379 200, 0.0, 1.0); 00380 00381 fp3dTrackChargeRatiovsCosz = new TProfile("TrackChargeRatiovsCosz", 00382 "3d track charge/record charge vs Hough CosZ", 00383 200, -1.0, 1.0); 00384 00385 fp3dTrackChargevsTrackCosz = new TProfile("TrackChargevsTrackCosz", 00386 "3d track charge vs Hough CosZ", 00387 200, -1.0, 1.0); 00388 00389 fp3dTrackChargevsNStrips = new TProfile("TrackChargevsNStrips", 00390 "3d track charge vs # of strips", 00391 100, 0.0, 100.0); 00392 00393 fp3dTrackNStripsvsCosz = new TProfile("TrackNStripsvsCosz", 00394 "3d track # of strip vs Hough CosZ", 00395 200, -1.0, 1.0); 00396 00397 fpTPosSigmavsCosZVview = new TProfile("TPosSigmavsCosZVview", 00398 "Error in tpos vs cosz V planes", 00399 100, -1.0, 1.0); 00400 SetAxisLabels(fpTPosSigmavsCosZVview, "error in tpos (m)"); 00401 00402 fpTPosSigmavsCosZUview = new TProfile("TPosSigmavsCosZUview", 00403 "Error in tpos vs cosz U planes", 00404 100, -1.0, 1.0); 00405 SetAxisLabels(fpTPosSigmavsCosZUview, "error in tpos (m)"); 00406 00407 fpTPosSigmavsTrackChargeVview = new TProfile("TPosSigmavsTrackChargeVview", 00408 "Error in tpos vs track charge V planes", 00409 1000, 0.0, 100000); 00410 SetAxisLabels(fpTPosSigmavsTrackChargeVview, "error in tpos (m)", "2d track charge (ADC)"); 00411 00412 fpTPosSigmavsTrackChargeUview = new TProfile("TPosSigmavsTrackChargeUview", 00413 "Error in tpos vs # track charge U planes", 00414 1000, 0.0, 100000); 00415 SetAxisLabels(fpTPosSigmavsTrackChargeUview, "error in tpos (m)", "2d track charge (ADC)"); 00416 00417 fpTPosSigmavsTrackNStripsVview = new TProfile("TPosSigmavsTrackNStripsVview", 00418 "Error in tpos vs # of track V strips", 00419 100, 0.0, 100); 00420 SetAxisLabels(fpTPosSigmavsTrackNStripsVview, "error in tpos (m)", "#strips"); 00421 00422 fpTPosSigmavsTrackNStripsUview = new TProfile("TPosSigmavsTrackNStripsUview", 00423 "Error in tpos vs # of track U strips", 00424 100, 0.0, 100); 00425 SetAxisLabels(fpTPosSigmavsTrackNStripsUview, "error in tpos (m)", "#strips"); 00426 00427 fpTPosSigmavsTrackChargeRatioVview = new TProfile("TPosSigmavsTrackChargeRatioVview", 00428 "Error in tpos vs # track/record charge V planes", 00429 100, 0.0, 1.0); 00430 SetAxisLabels(fpTPosSigmavsTrackChargeRatioVview, "error in tpos (m)"); 00431 00432 fpTPosSigmavsTrackChargeRatioUview = new TProfile("TPosSigmavsTrackChargeRatioUview", 00433 "Error in tpos vs # track/record charge U planes", 00434 100, 0.0, 1.0); 00435 SetAxisLabels(fpTPosSigmavsTrackChargeRatioUview, "error in tpos (m)"); 00436 00437 fhTrackChargeOverTotalChargeRatio = new TH1D("TrackChargeOverTotalChargeRatio", 00438 "Ratio of 3d track charge over total record charge", 00439 200, 0.0, 1.0); 00440 00441 fhCharge2dTracksV = new TH1D("Charge2dTracksV", 00442 "Charge in 2d V view tracks", 00443 fNTrackChargeBins, 0, fTrackChargeLimit); 00444 SetAxisLabels(fhCharge2dTracksV, "Raw ADC"); 00445 00446 fhCharge2dTracksU = new TH1D("Charge2dTracksU", 00447 "Charge in 2d U view tracks", 00448 fNTrackChargeBins, 0, fTrackChargeLimit); 00449 SetAxisLabels(fhCharge2dTracksU, "Raw ADC"); 00450 00451 fhCharge2dTracksDiffNorm = new TH1D("Charge2dTracksDiffNorm", 00452 "(vcharge - ucharge)/(vcharge + ucharge)", 00453 100, -1.0, 1.0); 00454 00455 fhNStrip2dTracksV = new TH1D("NStrip2dTracksV", "Number of strips in 2d V view tracks", 80, 0, 80); 00456 SetAxisLabels(fhNStrip2dTracksV, "Raw ADC"); 00457 00458 fhNStrip2dTracksU = new TH1D("NStrip2dTracksU", "Number of strips in 2d U view tracks", 80, 0, 80); 00459 SetAxisLabels(fhNStrip2dTracksU, "Raw ADC"); 00460 00461 fhNStrip2dTracksDiffNorm = new TH1D("NStrip2dTracksDiffNorm", 00462 "(vnstrip - unstrip)/(vnstrip + unstrip)", 00463 100, -1.0, 1.0); 00464 00465 fhTrackHoughCosU = new TH1D("TrackHoughCosU", "Residual rms of 2d tracks in V view", 00466 fNCosBins, -1, 1); 00467 fhTrackHoughCosV = new TH1D("TrackHoughCosV", 00468 "Residual rms of 2d tracks in V view", 00469 fNCosBins, -1, 1); 00470 fhTrackHoughCosZ = new TH1D("TrackHoughCosZ", 00471 "Residual rms of 2d tracks in V view", 00472 fNCosBins, -1, 1); 00473 fhSigmaOfTPosVview = new TH1D("SigmaOfTPosVview", 00474 "Error in tpos from fit in V view", 00475 fTrackResidualRMSBins, 0, fTrackResidualRMSLimit); 00476 SetAxisLabels(fhSigmaOfTPosVview, "tpos uncertainty (m)"); 00477 00478 fhSigmaOfTPosUview = new TH1D("SigmaOfTPosUview", 00479 "Error in tpos from fit in U view", 00480 fTrackResidualRMSBins, 0, fTrackResidualRMSLimit); 00481 SetAxisLabels(fhSigmaOfTPosUview, "tpos uncertainty (m)"); 00482 00483 fTrackResidualRMSVview = new TH1D ("TrackResidualRMSVview", 00484 "Roughly RMS of residuals along V view tracks", 00485 fTrackResidualRMSBins, 0, fTrackResidualRMSLimit); 00486 SetAxisLabels(fTrackResidualRMSVview, "track rms residual (m)"); 00487 00488 fhSigmaOfTPosDiffNorm = new TH1D("SigmaOfTPosDiffNorm", 00489 "(vsigma - usigma)/(vsigma + usigma)", 00490 100, -1.0, 1.0); 00491 00492 fTrackResidualRMSUview = new TH1D ("TrackResidualRMSUview", 00493 "Roughly RMS of residuals along U view tracks", 00494 fTrackResidualRMSBins, 0, fTrackResidualRMSLimit); 00495 SetAxisLabels(fTrackResidualRMSUview, "track rms residual (m)"); 00496 00497 fChargeRatioOf2dTracks = new TH1D ("ChargeRatioOf2dTracks", 00498 "Ratio of track charge, V over U", 00499 100, 0, 5); 00500 fRootFile->cd(); 00501 } 00502 00503 00504 const Double_t &ucharge = ntprec -> ucharge; 00505 const Double_t &vcharge = ntprec -> vcharge; 00506 const Double_t &vcandcharge = ntprec -> vcandcharge; 00507 const Double_t &ucandcharge = ntprec -> ucandcharge; 00508 const Double_t &vnstrip = ntprec -> ntrackvstrip; 00509 const Double_t &unstrip = ntprec -> ntrackustrip; 00510 const Double_t &vsigmaoftpos = ntprec -> vsigmaoftpos; 00511 const Double_t &usigmaoftpos = ntprec -> usigmaoftpos; 00512 00513 fhCharge2dTracksV -> Fill(vcharge); 00514 fhCharge2dTracksU -> Fill(ucharge); 00515 fhNStrip2dTracksV -> Fill(vnstrip); 00516 fhNStrip2dTracksU -> Fill(unstrip); 00517 fhSigmaOfTPosVview -> Fill(vsigmaoftpos); 00518 fhSigmaOfTPosUview -> Fill(usigmaoftpos); 00519 fTrackResidualRMSVview -> Fill(ntprec -> vtrackrms); 00520 fTrackResidualRMSUview -> Fill(ntprec -> utrackrms); 00521 fhTrackHoughCosV -> Fill(ntprec -> hcosv); 00522 fhTrackHoughCosU -> Fill(ntprec -> hcosu); 00523 fhTrackHoughCosZ -> Fill(ntprec -> hcosz); 00524 00525 if(vcharge + ucharge > 0.0) 00526 fhCharge2dTracksDiffNorm -> Fill((vcharge-ucharge)/(vcharge+ucharge)); 00527 00528 if(vnstrip + unstrip > 0.0) 00529 fhNStrip2dTracksDiffNorm -> Fill((vnstrip-unstrip)/(vnstrip+unstrip)); 00530 00531 if(vsigmaoftpos + usigmaoftpos > 0.0) 00532 fhSigmaOfTPosDiffNorm -> Fill((vsigmaoftpos-usigmaoftpos)/(vsigmaoftpos+usigmaoftpos)); 00533 00534 if(ucharge > 1.0) 00535 fChargeRatioOf2dTracks -> Fill(vcharge/ucharge); 00536 00537 const double sum = vcharge + ucharge + vcandcharge + ucandcharge; 00538 if(sum < 1.0){ 00539 MSG("Align", Msg::kError) << "Total record charge is less than 1.0."<< endl; 00540 return; 00541 } 00542 00543 const double ratio = (vcharge+ucharge)/sum; 00544 fhTrackChargeOverTotalChargeRatio -> Fill(ratio); 00545 00546 //Fill TProfile histograms 00547 fpTPosSigmavsCosZVview -> Fill(ntprec->hcosz, vsigmaoftpos); 00548 fpTPosSigmavsCosZUview -> Fill(ntprec->hcosz, usigmaoftpos); 00549 fpTPosSigmavsTrackChargeVview -> Fill(ntprec->vcharge, vsigmaoftpos); 00550 fpTPosSigmavsTrackChargeUview -> Fill(ntprec->ucharge, usigmaoftpos); 00551 fpTPosSigmavsTrackNStripsVview -> Fill(vnstrip, vsigmaoftpos); 00552 fpTPosSigmavsTrackNStripsUview -> Fill(unstrip, usigmaoftpos); 00553 fpTPosSigmavsTrackChargeRatioVview -> Fill(ratio, vsigmaoftpos); 00554 fpTPosSigmavsTrackChargeRatioUview -> Fill(ratio, usigmaoftpos); 00555 fp3dTrackChargevsTrackChargeRatio -> Fill(ratio, vcharge+ucharge); 00556 fp3dTrackChargevsTrackCosz -> Fill(ntprec->hcosz, vcharge+ucharge); 00557 fp3dTrackNStripsvsCosz -> Fill(ntprec->hcosz, vnstrip+unstrip); 00558 fp3dTrackChargevsNStrips -> Fill(vnstrip+unstrip, vcharge+ucharge); 00559 fp3dTrackChargeRatiovsCosz -> Fill(ntprec->hcosz, ratio); 00560 }

00595 { 00596 MSG("Align", Msg::kVerbose) << "AlignmentHistograms::SetAxisLabels()" << endl; 00597 00598 if(!h) return; 00599 00600 if(xlabel) 00601 { 00602 h -> GetXaxis() -> SetTitle(xlabel); 00603 h -> GetXaxis() -> CenterTitle(); 00604 } 00605 00606 if(ylabel) 00607 { 00608 h -> GetYaxis() -> SetTitle(ylabel); 00609 h -> GetYaxis() -> CenterTitle(); 00610 } 00611 }

Member Data Documentation

Generated on Thu Nov 1 11:56:36 2007 for loon by

1.3.9.1

AlignmentHistograms Class Reference

Public Member Functions

Private Member Functions

Private Attributes

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation


Public Member Functions
	AlignmentHistograms (TFile *rootfile, const VldContext &vld)
virtual	~AlignmentHistograms ()
void	Fill (const AlignmentStrip &astrip, const NtpAlignmentRecord *ntprec)
void	FillPlanePlex (const vector< AlignmentStrip > &track_strips, const vector< AlignmentStrip > &cand_strips)
Private Member Functions
void	FillTrackInfo (const NtpAlignmentRecord *ntprec)
void	FillHits (PlexStripEndId &plexid)
void	SetAxisLabels (TH1 h, const char xlabel=0, const char *ylabel=0)
Private Attributes
map< int, TH2D * >	f2hPlanePlexChargeRatio
map< int, TH2D * >	f2hPlanePlexHits
map< PlexStripEndId, unsigned int >	fNStripHits
map< PlexScintMdlId, unsigned int >	fNMdlHits
TFile *	fRootFile
TDirectory *	fPlanePlexDirectory
TDirectory *	fTrackDirectory
TDirectory *	fResidualSummaryDirectory
int	fnplanes
int	fnstrips
const int	fNResidualBins
const Double_t	fResidualLimit
const Double_t	fResidualStripBinWidth
const Double_t	fResidualMdlBinWidth
const int	fNResidualBinsStrip
const Double_t	fResidualLimitStrip
const int	fNChargeBins
const Double_t	fChargeLimit
const int	fNTrackChargeBins
const Double_t	fTrackChargeLimit
const int	fNHitsBins
const int	fTrackResidualRMSBins
const Double_t	fTrackResidualRMSLimit
const int	fNCosBins
TH1D *	fhCharge2dTracksU
TH1D *	fhCharge2dTracksV
TH1D *	fhCharge2dTracksDiffNorm
TH1D *	fhNStrip2dTracksU
TH1D *	fhNStrip2dTracksV
TH1D *	fhNStrip2dTracksDiffNorm
TH1D *	fhSigmaOfTPosVview
TH1D *	fhSigmaOfTPosUview
TH1D *	fhSigmaOfTPosDiffNorm
TH1D *	fChargesInAllHitStrip
TH1D *	fChargevsResidualHitStrips
TH1D *	fResidualAllStrips
TH1D *	fResidualVStrips
TH1D *	fResidualUStrips
TH1D *	fTrackResidualRMSVview
TH1D *	fTrackResidualRMSUview
TH1D *	fChargeRatioOf2dTracks
TH1D *	fhTrackHoughCosU
TH1D *	fhTrackHoughCosV
TH1D *	fhTrackHoughCosZ
TH1D *	fhPlexPlaneEntries
TH1D *	fhTrackChargeOverTotalChargeRatio
TH2D *	f2hNumberOfTrackHitsPerStrip
TH2D *	f2hNumberOfTrackHitsPerModule
TH2D *	f2hChargevsResidualHitStrips
TH2D *	f2hResidualvsPlane
TH2D *	f2hMeanChargePlaneStrip
TH2D *	f2hMeanPathCorrChargePlaneStrip
TProfile *	fpTPosSigmavsCosZVview
TProfile *	fpTPosSigmavsCosZUview
TProfile *	fpTPosSigmavsTrackChargeVview
TProfile *	fpTPosSigmavsTrackChargeUview
TProfile *	fpTPosSigmavsTrackNStripsVview
TProfile *	fpTPosSigmavsTrackNStripsUview
TProfile *	fpTPosSigmavsTrackChargeRatioVview
TProfile *	fpTPosSigmavsTrackChargeRatioUview
TProfile *	fp3dTrackChargevsTrackChargeRatio
TProfile *	fp3dTrackChargevsTrackCosz
TProfile *	fp3dTrackNStripsvsCosz
TProfile *	fp3dTrackChargevsNStrips
TProfile *	fp3dTrackChargeRatiovsCosz
TProfile *	fpResidualvsPlane
TProfile *	fpResidualvsCosZVview
TProfile *	fpResidualvsCosZUview
TProfile *	fpResidualvsTPosSigmaVview
TProfile *	fpResidualvsTPosSigmaUview
TProfile *	fpResidualvsStripChargeVview
TProfile *	fpResidualvsStripChargeUview
TProfile *	fpResidualvsTrackChargeVview
TProfile *	fpResidualvsTrackChargeUview
TProfile *	fpResidualvsTrackNStripsVview
TProfile *	fpResidualvsTrackNStripsUview
TProfile *	fpResidualvsTrackChargeRatioVview
TProfile *	fpResidualvsTrackChargeRatioUview
TProfile *	fpMeanStripResponsevsResidual