#include <FWTrackResidualDetailView.h>

Inheritance diagram for FWTrackResidualDetailView:

Public Member Functions
	FWTrackResidualDetailView ()

virtual	~FWTrackResidualDetailView ()

Public Member Functions inherited from FWDetailViewCanvas< reco::Track >
	FWDetailViewCanvas ()

virtual void	init (TEveWindowSlot *)

virtual	~FWDetailViewCanvas ()

Public Member Functions inherited from FWDetailView< reco::Track >
	FWDetailView ()

Public Member Functions inherited from FWDetailViewBase
void	build (const FWModelId &)

const fireworks::Context &	context () const

const FWEventItem *	item ()

virtual void	setBackgroundColor (Color_t col)

void	setItem (const FWEventItem *x)

virtual	~FWDetailViewBase ()

Private Member Functions
virtual void	build (const FWModelId &id, const reco::Track *)

	FWTrackResidualDetailView (const FWTrackResidualDetailView &)

double	getSignedResidual (const FWGeometry *geom, unsigned int id, double resX)

const FWTrackResidualDetailView &	operator= (const FWTrackResidualDetailView &)

void	prepareData (const FWModelId &id, const reco::Track *)

void	printDebug ()

virtual void	setTextInfo (const FWModelId &id, const reco::Track *)

Private Attributes
int	hittype [64]

int	m_det [64]

int	m_detector [64]

Color_t	m_invalidCol

Int_t	m_invalidFill

int	m_ndet

int	m_nhits

Color_t	m_resXCol

Int_t	m_resXFill

Color_t	m_resYCol

Int_t	m_resYFill

Color_t	m_stereoCol

Int_t	m_stereoFill

float	res [2][64]

int	stereo [64]

int	substruct [64]

int	subsubstruct [64]

Static Private Attributes
static const char *	m_det_tracker_str [] ={"PXB","PXF","TIB","TID","TOB","TEC"}

Additional Inherited Members
Static Public Member Functions inherited from FWDetailView< reco::Track >
static std::string	classRegisterTypeName ()

static std::string	classTypeName ()

Static Public Member Functions inherited from FWDetailViewBase
static void	drawCanvasBox (Double_t *pos, Color_t fillCol, Int_t fillType=0, bool bg=kTRUE)

static void	drawCanvasDot (Float_t x, Float_t y, Float_t r, Color_t)

Protected Member Functions inherited from FWDetailViewBase
	FWDetailViewBase (const std::type_info &)

Protected Attributes inherited from FWDetailViewCanvas< reco::Track >
TGCompositeFrame *	m_guiFrame

TCanvas *	m_infoCanvas

TCanvas *	m_viewCanvas

Detailed Description

Definition at line 46 of file FWTrackResidualDetailView.h.

Constructor & Destructor Documentation

FWTrackResidualDetailView::FWTrackResidualDetailView ( )

Definition at line 28 of file FWTrackResidualDetailView.cc.

References hittype, m_det, m_detector, res, stereo, substruct, and subsubstruct.

                                                      :
    m_ndet(0),
    m_nhits(0),
    m_resXFill(3007),
    m_resXCol(kGreen-9),
    m_resYFill(3006),
    m_resYCol(kWhite),
    m_stereoFill(3004),
    m_stereoCol(kCyan-9),
    m_invalidFill(3001),
    m_invalidCol(kRed)
 {
    memset(m_det,        0, sizeof(m_det));
    memset(res,          0, sizeof(res));
    memset(hittype,      0, sizeof(hittype));
    memset(stereo,       0, sizeof(stereo));
    memset(substruct,    0, sizeof(substruct));
    memset(subsubstruct, 0, sizeof(subsubstruct));
    memset(m_detector,   0, sizeof(m_detector));
 }

FWTrackResidualDetailView::~FWTrackResidualDetailView ( )

virtual

Definition at line 49 of file FWTrackResidualDetailView.cc.

50 {

51 }

FWTrackResidualDetailView::FWTrackResidualDetailView ( const FWTrackResidualDetailView & )

private

Member Function Documentation

void FWTrackResidualDetailView::build	(	const FWModelId &	id,
		const reco::Track *	track
	)

privatevirtual

Implements FWDetailView< reco::Track >.

Definition at line 102 of file FWTrackResidualDetailView.cc.

References create_public_lumi_plots::color, diffTreeTool::diff, FWDetailViewBase::drawCanvasBox(), reco::Track::extra(), fwLog, h, hittype, edm::Ref< C, T, F >::isAvailable(), j, gen::k, fwlog::kError, ConfigFiles::l, groupFilesInBlocks::lines, m_det, m_det_tracker_str, m_invalidCol, m_invalidFill, m_ndet, m_nhits, m_resXCol, m_resXFill, m_resYCol, m_resYFill, m_stereoCol, m_stereoFill, FWDetailViewCanvas< reco::Track >::m_viewCanvas, prepareData(), RecoTauCleanerPlugins::pt, res, setTextInfo(), stereo, substruct, subsubstruct, and create_public_lumi_plots::width.

 {
    if (!track->extra().isAvailable()) {
       fwLog(fwlog::kError) << " no track extra information is available.\n";
      m_viewCanvas->cd();
      TLatex* latex = new TLatex();
      latex->SetTextSize(0.1);
      latex->DrawLatex(0.1, 0.5, "No track extra information");
      return;
    }
    prepareData(id, track);
 
    // draw histogram
    m_viewCanvas->cd();
    m_viewCanvas->SetHighLightColor(-1);
    TH2F* h_res = new TH2F("h_resx","h_resx",10,-5.5,5.5,m_ndet,0,m_ndet);
    TPad* padX = new TPad("pad1","pad1", 0.2, 0., 0.8, 0.99);
    padX->SetBorderMode(0);
    padX->SetLeftMargin(0.2);
    padX->Draw();
    padX->cd();
    padX->SetFrameLineWidth(0);
    padX->Modified();
    h_res->SetDirectory(0);
    h_res->SetStats(kFALSE);
    h_res->SetTitle("");
    h_res->SetXTitle("residual");
    h_res->GetXaxis()->SetTickLength(0);
    h_res->GetYaxis()->SetTickLength(0);
    h_res->GetXaxis()->SetNdivisions(20);
    h_res->GetYaxis()->SetLabelSize(0.06);
    h_res->Draw();
    padX->SetGridy();
 
    float larray[9]={0.0, 0.5, 1.0, 1.5, 2.0, 2.5, 3.5, 4.5, 5.5};
    float larray2[8];
    for(int l=0; l<8; l++) {
       float diff2=(larray[l+1]-larray[l])/2;
       larray2[l]=larray[l]+diff2;
       //  printf("(%.1f,%.1f),",larray[i],larray[i+1]);
    }
 
    int resi[2][64];
    for(int l=0; l<m_nhits; l++) {
       for(int k=0; k<8; k++) {
          if(fabs(res[0][l])==larray2[k])
             resi[0][l]=k;
          if(fabs(res[1][l])==larray2[k])
             resi[1][l]=k;
       }
    }
 
    TLine* lines[17];
    for(int l=0; l<17; l++) {
       int ix=l%9;
       int sign=1;
       sign = (l>8) ? -1 : 1;
       lines[l] = new TLine(sign*larray[ix],0,sign*larray[ix],m_ndet);
       if(l!=9)
          lines[l]->SetLineStyle(3);
       padX->cd();
       lines[l]->Draw();
    }
 
    float width=0.25;
    int filltype;
    Color_t color;
    Double_t box[4];
    padX->cd();
 
    for(int h=0; h<2; h++) {
       float height1=0;
       for(int j=0; j<m_ndet; j++) {
          // take only X res and Y pixel residals
          if (strcmp(m_det_tracker_str[substruct[m_det[j]]-1], "PXB") && h)
             continue;
 
          char det_str2[256];
          snprintf(det_str2, 255, "%s/%i",m_det_tracker_str[substruct[m_det[j]]-1],subsubstruct[m_det[j]]);
          h_res->GetYaxis()->SetBinLabel(j+1, det_str2);
 
          int diff=m_det[j+1]-m_det[j];
          int k=0;
          width=1.0/diff;
 
          for(int l=m_det[j]; l<(m_det[j]+diff); l++) {
             //      g->SetPoint(l,resx[l],j+0.5);
             //  printf("%i, %f %f %f\n",l,resx[l],sign*larray[resxi[l]],sign*larray[resxi[l]+1]);
             int sign = (res[h][l]<0) ? -1 : 1;
             box[0] = (hittype[l]==0) ? sign*larray[resi[h][l]] : -5.5;
             box[2] = (hittype[l]==0) ? sign*larray[resi[h][l]+1] : 5.5;
             box[1] = height1+width*k;
             box[3] = height1+width*(k+1);
 
             if(stereo[l]==1) {
                color    = m_stereoCol;
                filltype = m_stereoFill;
             }
             else if(hittype[l]!=0) {
                color    = m_invalidCol;
                filltype = m_invalidFill;
             }
             else {
                filltype = h ? m_resYFill : m_resXFill;
                color    = h ? m_resYCol  : m_resXCol;
             }
 
             drawCanvasBox(box, color, filltype, h<1);
             k++;
          }
          height1 +=1;
       }
    }
    
    //  title
    const char* res_str= "residuals in Si detector local x-y coord.";
    TPaveText *pt = new TPaveText(0.0,0.91, 1,0.99,"blNDC");
    pt->SetBorderSize(0);
    pt->SetFillColor(kWhite);
    pt->AddText(res_str);
    pt->Draw();
    
    m_viewCanvas->cd();
    m_viewCanvas->SetEditable(kFALSE);
 
    setTextInfo(id, track);
 }

double FWTrackResidualDetailView::getSignedResidual	(	const FWGeometry *	geom,
		unsigned int	id,
		double	resX
	)

private

Definition at line 231 of file FWTrackResidualDetailView.cc.

References diffTwoXMLs::g1, diffTwoXMLs::g2, FWGeometry::getMatrix(), and m.

Referenced by prepareData().

 {
    double local1[3] = { 0, 0, 0 };
    double local2[3] = { resX, 0, 0 };
    double global1[3], global2[3];
    const TGeoMatrix *m = geom->getMatrix(id);
    assert(m != 0);
    m->LocalToMaster(local1, global1);
    m->LocalToMaster(local2, global2);
    TVector3 g1 = global1;
    TVector3 g2 = global2;
    if (g2.DeltaPhi(g1) > 0)
       return resX;
    else return -resX;
 }

const FWTrackResidualDetailView& FWTrackResidualDetailView::operator= ( const FWTrackResidualDetailView & )

private

void FWTrackResidualDetailView::prepareData	(	const FWModelId &	id,
		const reco::Track *	track
	)

private

Definition at line 54 of file FWTrackResidualDetailView.cc.

References relativeConstraints::geom, reco::HitPattern::getHitPattern(), getSignedResidual(), reco::TrackBase::hitPattern(), hittype, i, j, gen::k, m_det, m_detector, m_ndet, m_nhits, reco::HitPattern::numberOfHits(), reco::Track::recHitsBegin(), res, reco::Track::residuals(), reco::TrackResiduals::residualX(), reco::TrackResiduals::residualY(), stereo, substruct, and subsubstruct.

Referenced by build().

 {
    HitPattern hitpat = track->hitPattern();
    TrackResiduals residuals = track->residuals();
 
    const FWGeometry *geom = id.item()->getGeom();
    assert(geom != 0);
    m_nhits=hitpat.numberOfHits();
    for (int i = 0; i < m_nhits; ++i) {
       //    printf("there are %d hits in the pattern, %d in the vector, this is %u\n",
       //           m_nhits, track->recHitsEnd() - track->recHitsBegin(), (*(track->recHitsBegin() + i))->geographicalId().rawId());
       hittype[i] = 0x3 & hitpat.getHitPattern(i);
       stereo[i] = 0x1 & hitpat.getHitPattern(i) >> 2;
       subsubstruct[i] = 0xf & hitpat.getHitPattern(i) >> 3;
       substruct[i] = 0x7 & hitpat.getHitPattern(i) >> 7;
       m_detector[i] = 0x01 & hitpat.getHitPattern(i) >> 10;
       if ((*(track->recHitsBegin() + i))->isValid()) {
          res[0][i] = getSignedResidual(geom,
                                        (*(track->recHitsBegin() + i))->geographicalId().rawId(),
                                        residuals.residualX(i, hitpat));
       } else {
          res[0][i] = 0;
       }
       res[1][i] = residuals.residualY(i, hitpat);
       // printf("%s, %i\n",m_det_tracker_str[substruct[i]-1],subsubstruct[i]);
    }
 
    m_det[0]=0;
    for(int j=0; j < m_nhits-1;) {
       int k=j+1;
       for(; k<m_nhits ; k++) {
          if(substruct[j]==substruct[k]  && subsubstruct[j]==subsubstruct[k]) {
             if(k==(m_nhits-1)) j=k;
          }
          else {
             m_ndet++;
             j=k;
             m_det[m_ndet]=j;
             break;
          }
       }
    }
    m_ndet++;
    m_det[m_ndet]=m_nhits;
    // printDebug();
 }

void FWTrackResidualDetailView::printDebug ( )

private

Definition at line 248 of file FWTrackResidualDetailView.cc.

References gather_cfg::cout, i, m_det, m_det_tracker_str, m_ndet, substruct, and subsubstruct.

 {
    for(int i=0; i<m_ndet; i++)
    {
       std::cout <<"<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<" << std::endl;
       std::cout << "idx " << i << " det[idx] " <<  m_det[i] << std::endl;
       std::cout << "m_det idx " << m_det[i] <<   std::endl;
       std::cout << "m_det_tracker_str idx " << substruct[m_det[i]]-1  << std::endl;
       printf("m_det[idx] %i m_det_tracker_str %s substruct %i\n",m_det[i], m_det_tracker_str[substruct[m_det[i]]-1], subsubstruct[m_det[i]]);
    }
 }

void FWTrackResidualDetailView::setTextInfo	(	const FWModelId &	id,
		const reco::Track *
	)

privatevirtual

Implements FWDetailView< reco::Track >.

Definition at line 261 of file FWTrackResidualDetailView.cc.

References FWDetailViewBase::drawCanvasBox(), hittype, i, FWDetailViewCanvas< reco::Track >::m_infoCanvas, m_invalidCol, m_invalidFill, m_ndet, m_nhits, m_resXCol, m_resXFill, m_resYCol, m_resYFill, m_stereoCol, m_stereoFill, substruct, and detailsBasic3DVector::y.

Referenced by build().

 {
    m_infoCanvas->cd();
 
    char mytext[256];
    Double_t fontsize = 0.07;
    TLatex* latex = new TLatex();
    latex->SetTextSize(fontsize);
    latex->Draw();
    Double_t x0 = 0.02;
    Double_t y = 0.95;
 
    // summary
    int nvalid=0;
    int npix=0;
    int nstrip=0;
    for(int i=0; i<m_nhits; i++)
    {
       if(hittype[i]==0) nvalid++;
       if(substruct[i]<3) npix++;
       else nstrip++;
    }
 
    latex->SetTextSize(fontsize);
    Double_t boxH = 0.25*fontsize;
    
    double yStep = 0.04;
 
    latex->DrawLatex(x0, y, "Residual:");
    y-= yStep;
    latex->DrawLatex(x0, y, "sgn(#hat{X}#bullet#hat{#phi}) #times #frac{X_{hit} - X_{traj}}{#sqrt{#sigma^{2}_{hit} + #sigma^{2}_{traj}}}" );
    y-= 2.5*yStep;
    snprintf(mytext, 255, "layers hit: %i", m_ndet);
    latex->DrawLatex(x0, y, mytext);
    y -= yStep;
    snprintf(mytext, 255,"valid Si hits: %i", nvalid);
    latex->DrawLatex(x0, y, mytext);
    y -= yStep;
    snprintf(mytext, 255,"total Si hits: %i", m_nhits);
    latex->DrawLatex(x0, y, mytext);
    y -= yStep;
    snprintf(mytext, 255,"valid Si pixel hits: %i", npix);
    latex->DrawLatex(x0, y, mytext);
    y -= yStep;
    snprintf(mytext, 255, "valid Si strip hits: %i", nstrip);
    latex->DrawLatex(x0, y, mytext);
    
 
    Double_t pos[4];
    pos[0] = 0.4;
    pos[2] = 0.55;
 
    y -= yStep*2;
    latex->DrawLatex(x0, y, "X hit");
    pos[1] = y; pos[3] = pos[1] + boxH;
    drawCanvasBox(pos, m_resXCol, m_resXFill);
 
    y -=  yStep;
    latex->DrawLatex(x0, y, "Y hit");
    pos[1] = y; pos[3] = pos[1] + boxH;
    drawCanvasBox(pos, m_resYCol, m_resYFill, 0);
 
    y -= yStep;
    latex->DrawLatex(x0, y, "stereo hit");
    pos[1] = y; pos[3] = pos[1] + boxH;
    drawCanvasBox(pos, m_stereoCol, m_stereoFill);
 
    y -= yStep;
    latex->DrawLatex(x0, y, "invalid hit");
    pos[1] = y; pos[3] = pos[1] + boxH;
    drawCanvasBox(pos, m_invalidCol, m_invalidFill);
 }