---

DisplayCanvas.cxx

Go to the documentation of this file.

#include "bfemDisplay/EDControl.h"
#include "bfemDisplay/DisplayCanvas.h"

// Options for AnnotatePad() method.

#define kDrawTitle 1
#define kSetRange 2
#define kShowAxes 4


extern EDControl *m_evtCtrl;   // Copy of the event control class object.

// 


ClassImp(DisplayCanvas)

DisplayCanvas::DisplayCanvas() {




   m_TowerHeight     = 38.81;
   m_TowerWidth      = 16.5; 

   frame_x           = m_TowerWidth;
   frame_y           = m_TowerWidth;
   frame_z           = m_TowerHeight;

   m_eventTextList = 0;   // List of TPaveTexts for one event.


   // Data trees.

   fTree = 0;
   rTree = 0;

   m_event = 0;
   m_recon = 0;

   // CAL, ACD, XGT pedestals.

   m_ped = m_evtCtrl->getCALPeds();

   m_acdped = m_evtCtrl->getACDPeds();

   // m_xgtped = m_evtCtrl->getXGTPeds();
    

   // Get pointer to display canvas.
   m_canvas = m_evtCtrl->getCanvas();


   // Make a top-level reference node ( old "ttt" node; seems necessary in order to properly set color
   // of  subsequent TBRIKs ).
   TBRIK* top = new TBRIK( "top", "top", "void", 0.001, 0.001, 0.001 );

   TNode* m_topNode = new TNode( "topNode", "topNode", "top" );

   // Set up pads
   //   setDrawingPads();

   // Create the reference frame.
   TBRIK *frame = new TBRIK( "frame", "frame", "void", frame_x, frame_y, frame_z );
   frame->SetLineColor( 10 );    // Do not draw - color is white.

   m_frameNode = new TNode( "frameNode", "frameNode", "frame" );
  
   // Create frame representing the tower.
   TBRIK *tower = new TBRIK( "tower", "tower", "void", m_TowerWidth, m_TowerWidth, m_TowerHeight );
   tower->SetLineColor( 10 );  // Do not draw - color is white.

   m_towerNode = new TNode( "towerNode", "towerNode", "tower" );
   m_towerNode->cd();

   // Create node and shape structure for the tracker.
   m_Tkr = new TkrDisplay( m_towerNode );

   // Create node and shape structure for the CAL.
   m_CAL = new CalDisplay( m_towerNode, m_ped );

   // Create node and shape structure for the CAL.
   m_Acd = new AcdDisplay( m_towerNode, m_acdped );




   m_currentevent = 0;

   m_scanForward = kTRUE;


}





void DisplayCanvas::Draw( Bool_t dataType ) {


  // Refresh data tree pointers.

  fTree = m_evtCtrl->getRawTree();
  rTree = m_evtCtrl->getReconTree();

  m_event = m_evtCtrl->getEvent();
  m_recon = m_evtCtrl->getRecon();

  if ( m_event ) m_event->Clean();


  if ( fTree ) {
    
    if ( ( m_currentevent < 0 ) || ( m_currentevent > m_evtCtrl->getEntries() ) ) return;

    fTree->SetBranchStatus("*",0);  // disable all branches
    fTree->SetBranchStatus("m_Acd*",1);  // activate branchname
    fTree->SetBranchStatus("m_Xgt*",1);  // activate branchname
    fTree->SetBranchStatus("m_Tkr*",1);  // activate branchname
    fTree->SetBranchStatus("m_Cal*",1);  // activate branchname
    fTree->SetBranchStatus("m_eventId", 1); // activate branchname
    fTree->SetBranchStatus("m_run", 1); // activate branchname
    fTree->GetEvent(m_currentevent);

  }

  if ( rTree ) {
    TBranch* TKR = rTree->GetBranch("m_tkr");
    TKR->GetEvent( m_currentevent );
    TBranch* CALrec = rTree->GetBranch("m_cal");
    CALrec->GetEvent( m_currentevent );
  }

  m_canvas->Clear();
  setDrawingPads();


  if ( fTree ) {
    m_Tkr->DrawEvents( m_event, m_recon );
    m_CAL->DrawEvents( m_event, m_recon, dataType );
    m_Acd->DrawEvents( m_event );
  }


  // Add some info at top of canvas.

  addEnergyScale();


  //
  // Draw each of the various views in different pads on the canvas.
  //


  // Delete old paves for previous event.
  if ( m_eventTextList ) m_eventTextList->Delete();
  m_eventTextList = new TObjArray();


   // ------------------------------------- 3D display

   m_3DPad->cd();

   m_Tkr->setVisibility( kHideNone );  // All hits visible.
   m_CAL->setVisibility( kHideNone );

   m_towerNode->DrawClone();

   // N.B. Need DrawClone, otherwise visibility pointers for this pad are changed when re-drawing tower
   // on next pad.

   AnnotatePad( m_3DPad, "3D", kSetRange );


   DrawEventInfo( m_3DPad );   // Show some info about this event.


   // ------------------------------------- Front - ZX projection.

   m_zxpad->cd();

   m_Tkr->setVisibility( kHideY );  // Hide hits on strips parallel to X-axis ( Y measured coord. )
   m_CAL->setVisibility( kHideY );

   m_towerNode->DrawClone();

   AnnotatePad( m_zxpad, "Z-X", kDrawTitle | kSetRange );    // Add title and axes.

  // ------------------------------------- Side - ZY projection.

   m_zypad->cd();

   m_Tkr->setVisibility( kHideX );  // Hide hits on strips parallel to Y-axis ( X measured coord. )
   m_CAL->setVisibility( kHideX );

   m_towerNode->DrawClone();

   AnnotatePad( m_zypad, "Z-Y", kDrawTitle | kSetRange );

  // ------------------------------------- Top - XY projection.

   m_xyTopPad->cd();

   m_Tkr->setVisibility( kHideAll );
   m_CAL->setVisibility( kHideAll );

   m_towerNode->DrawClone();

   AnnotatePad( m_xyTopPad, "X-Y",  kDrawTitle | kSetRange );    // Add title and axes.

  // ------------------------------------- Top - YX projection.

   m_yxTopPad->cd();

   m_Tkr->setVisibility( kHideAll );
   m_CAL->setVisibility( kHideAll );

   m_towerNode->Draw();  // Last can and MUST be normal draw ( otherwise seg. violation! )

   AnnotatePad( m_yxTopPad, "Y-X",  kDrawTitle | kSetRange );    // Add title and axes.


  

   printf("\nFinished displaying event %i\n",m_currentevent);
  

   if ( m_scanForward ) m_currentevent++;

   else m_currentevent--;

  
   //   gObjectTable->Print();

   m_3DPad->cd();

   //   m_canvas->SetCanvasSize( 800, 840 );

   m_canvas->Update();


}




void DisplayCanvas::setDrawingPads() {



   m_infoPad = new TPad("mainpad","GLAST Display",0.0,0.95,1.0,1.0);
   m_infoPad->Draw();



   m_3DPad = new TPad("3DPad","GLAST Display",0.0,0.0,0.33,0.95);
   m_3DPad->Draw();
   m_3DPad->SetTheta(10.0); 
   m_3DPad->SetPhi(30.0);



   m_zxpad = new TPad("zxpad","GLAST Display",0.33,0.0,0.66,0.64);
   m_zxpad->Draw();
   m_zxpad->SetTheta(0.); 
   m_zxpad->SetPhi(0.);

   

   m_zypad = new TPad("zypad","GLAST Display",0.66,0.0,1.0,0.64);
   m_zypad->Draw();
   m_zypad->SetPhi(-90.); 
   m_zypad->SetTheta(0.);



   m_xyTopPad = new TPad("xyTopPad","GLAST Display",0.33,0.64,0.66,0.95);
   m_xyTopPad->Draw();
   m_xyTopPad->SetTheta(90.); 
   m_xyTopPad->SetPhi(0.);



   m_yxTopPad = new TPad("yxTopPad","GLAST Display",0.66,0.64,1.0,0.95);
   m_yxTopPad->Draw();
   m_yxTopPad->SetTheta(90.); 
   m_yxTopPad->SetPhi(-90.);

}


void DisplayCanvas::AnnotatePad( TPad *pad, Char_t *title, UShort_t opt ) {

  // Add a title and possibly axes to a pad.

  pad->cd();

  TView *view;
  TPaveText *text;


  if ( opt & kDrawTitle ) {

    Bool_t isTopView = !strcmp( title, "X-Y" ) || !strcmp( title, "Y-X" );

    if ( isTopView ) text = new TPaveText( 0.80, 0.90, 0.95, 0.98, "NDC" );
    else text = new TPaveText( 0.80, 0.95, 0.95, 0.99, "NDC" );

    text->SetFillColor( 42 );

    text->AddText( title );
    text->SetAllWith( title, "Size", 0.07 );

    text->Draw();

    m_eventTextList->Add( text );

  }

  view = pad->GetView();

  if ( opt & kSetRange ) {

    if ( !strcmp( title, "3D" ) ) view->SetRange( -23.0, -23.0, -35.0, 23.0, 23.0, 55.0 ); // Range for 3D view.
    else view->SetRange( -23.0, -23.0, -20.0, 23.0, 23.0, 50.0 );  // Range for all other views.

  }

  if ( opt & kShowAxes ) {  // This option causes crash for some reason!

    view->ShowAxis();

  }

}




void DisplayCanvas::DrawEventInfo( TPad *pad ) {

  pad->cd();

  Char_t evt_num[ 50 ];


  // Allert if there was a gamma convertion
  if (m_CAL->m_Gammas) {

    TPaveText *gammas = new TPaveText( 0., 0.12, 0.9, 0.15, "NDC" );

    gammas->SetTextColor(2);
    gammas->SetBorderSize( 0 );
    gammas->SetTextSize(0.05);

    sprintf( evt_num, "Gamma convertion!!!");

    gammas->AddText( evt_num );

    gammas->Draw();

    m_eventTextList->Add( gammas );  // Add to list of TPaveTexts for event.

  }

  // Info on reconstructed energy in CAL.

  Float_t eTot=0, eCorr=0, eLeak=0;

  if ( rTree ) {

    eTot = m_CAL->m_eTot;
    eCorr = m_CAL->m_eCorr;
    eLeak = m_CAL->m_eLeak;

  

    TPaveText *ecal = new TPaveText( 0., 0., 0.9, 0.1, "NDC" );

    ecal->SetFillColor( 0 );
    ecal->SetBorderSize( 0 );
    ecal->SetTextSize(0.05);
    ecal->SetTextAlign(12);

    sprintf( evt_num, "CAL: Etot(GeV): %4.2f\n", eTot);

    ecal->AddText( evt_num );

    ecal->Draw();

    sprintf( evt_num, "      Ecorr(GeV): %4.2f\n", eCorr);
 
    ecal->AddText( evt_num );

    ecal->Draw();

    sprintf( evt_num, "      Eleak(GeV): %4.2f\n", eLeak);
 
    ecal->AddText( evt_num );
    
    ecal->Draw();

    m_eventTextList->Add( ecal );  // Add to list of TPaveTexts for event.

  }

  // Event identifier and run number

  TPaveText *text = new TPaveText( 0.0, 0.95, 0.85, 0.98, "NDC" );

  text->SetFillColor( 0 );
  text->SetBorderSize( 0 );
  text->SetTextSize(0.05);
  text->SetTextAlign(12);

  Int_t runID;
  Int_t rawEventID;


  if ( m_event ) {

    runID = m_event->getRunId();
    rawEventID = m_event->getEventId();

  }

  else { runID = 0; rawEventID = 0; }

  sprintf( evt_num, " Event Id: %d  --  Run. %d ", rawEventID, runID );

  text->AddText( evt_num );

  text->Draw();

  m_eventTextList->Add( text );  // Add to list of TPaveTexts for event.

}


void DisplayCanvas::setThresholds( Float_t thresh[ 3 ][ 4 ] ) {

  m_cal_lev[ 0 ] = thresh[ 0 ][ 0 ];
  m_cal_lev[ 1 ] = thresh[ 0 ][ 1 ];
  m_cal_lev[ 2 ] = thresh[ 0 ][ 2 ];
  m_cal_lev[ 3 ] = thresh[ 0 ][ 3 ];

  m_CAL->setThresholds( m_cal_lev );

  m_acd_lev[ 0 ] = thresh[ 1 ][ 0 ];
  m_acd_lev[ 1 ] = thresh[ 1 ][ 1 ];
  m_acd_lev[ 2 ] = thresh[ 1 ][ 2 ];
  m_acd_lev[ 3 ] = thresh[ 1 ][ 3 ];

  m_xgt_lev[ 0 ] = thresh[ 2 ][ 0 ];
  m_xgt_lev[ 1 ] = thresh[ 2 ][ 1 ];
  m_xgt_lev[ 2 ] = thresh[ 2 ][ 2 ];
  m_xgt_lev[ 3 ] = thresh[ 2 ][ 3 ];

  m_Acd->setThresholds( m_acd_lev, m_xgt_lev );

}



void DisplayCanvas::addEnergyScale(){

   // Add energy scale for CAL and ACDs 

  Int_t colors[ 4 ] = { 3, 4, 2 ,5 }; // green,blue,red,yellow

  Float_t xOffset;

  Float_t x[ 5 ];
  Float_t y[ 5 ];

  y[ 0 ] = y[ 1 ] = y[ 4 ] = 0.3;
  y[ 2 ] = y[ 3 ] = 0.7;

  m_infoPad->cd();

  TPaveText *CALtext = new TPaveText( 0.6, 0.5, 0.69, 0.9, "NDC" );
  CALtext->SetTextSize(0.35);
  CALtext->SetFillColor( 0 );
  CALtext->SetBorderSize( 0 );
  CALtext->AddText( "XGT (MIP): " );
  CALtext->Draw();
   
  TPaveText *CALtext1 = new TPaveText( 0.6, 0.1, 0.69, 0.5, "NDC" );
  CALtext1->SetTextSize(0.35); 
  CALtext1->SetFillColor( 0 );
  CALtext1->SetBorderSize( 0 );
  CALtext1->SetTextAlign(32);
  CALtext1->AddText( "ACD (MIP) " );
  CALtext1->Draw();
  
  Char_t text[ 10 ];
  Float_t mip;

  for ( Int_t color=0; color<4; ++color ) {

    xOffset = 0.7 + 0.065*color;
    
    x[ 0 ] = x[ 3 ] = x[ 4 ] = xOffset;
    x[ 1 ] = x[ 2 ] = xOffset + 0.02;

    TPolyLine *box = new TPolyLine( 5, x, y );

    box->SetLineColor( colors[ color ] );

    box->Draw();

    TPaveText *CALtext3 = new TPaveText( xOffset+0.025, 0.5, xOffset+0.06, 0.9, "NDC" );
    CALtext3->SetTextSize(0.35);
    CALtext3->SetFillColor( 0 );
    CALtext3->SetBorderSize( 0 );
    sprintf( text, "%5.1f ", m_xgt_lev[ color ] );
    CALtext3->AddText( text );
    CALtext3->Draw();
    TPaveText *CALtext4 = new TPaveText( xOffset+0.025, 0.1, xOffset+0.06, 0.5, "NDC" );
    CALtext4->SetTextSize(0.35);
    CALtext4->SetFillColor( 0 );
    CALtext4->SetBorderSize( 0 );
    mip = m_acd_lev[ color ];
    sprintf( text, "%5.1f ", mip );
    CALtext4->AddText( text );
    CALtext4->Draw();
  }


  // Add energy scale for XGTs (in MIP)

  TPaveText *CALtext5 = new TPaveText( 0.0, 0.3, 0.1, 0.8, "NDC" );
  CALtext5->SetTextSize(0.35); 
  CALtext5->SetFillColor( 0 );
  CALtext5->SetBorderSize( 0 );
  CALtext5->AddText( "CAL (MIP) " );
  CALtext5->Draw();


  for ( Int_t color=0; color<4; ++color ) {

    xOffset = 0.1 + 0.065*color;
    
    x[ 0 ] = x[ 3 ] = x[ 4 ] = xOffset;
    x[ 1 ] = x[ 2 ] = xOffset + 0.02;

    TPolyLine *box = new TPolyLine( 5, x, y );

    box->SetLineColor( colors[ color ] );

    box->Draw();

    TPaveText *CALtext4 = new TPaveText( xOffset+0.025, 0.3, xOffset+0.06, 0.8, "NDC" );
    CALtext4->SetTextSize(0.35);
    CALtext4->SetFillColor( 0 );
    CALtext4->SetBorderSize( 0 );
    mip = m_cal_lev[ color ];
    sprintf( text, "%5.1f ", mip );
    CALtext4->AddText( text );
    CALtext4->Draw();

  }

}


void DisplayCanvas::drawAxes(){



   TBRIK* x_axis = new TBRIK("x_axis","x_axis","void",frame_x/2.0,m_axis_width,m_axis_width);

   x_axis->SetLineColor(5);

   TBRIK* y_axis = new TBRIK("y_axis","y_axis","void",m_axis_width,frame_y/2.0,m_axis_width);

   y_axis->SetLineColor(6);

   TBRIK* z_axis = new TBRIK("z_axis","z_axis","void",m_axis_width,m_axis_width,frame_z/2.0);

   z_axis->SetLineColor(6);

}




Char_t* DisplayCanvas::getText(TString color, Float_t energy){

   Char_t lev[40];



   m_TextLine = new TString(color);

   m_TextLine->Append(" > ");

   sprintf(lev,"%6.4f",energy);

   m_TextLine->Append(lev);

   m_TextLine->Append(" GeV"); 

   return (Char_t*)m_TextLine->Data();



}




Char_t* DisplayCanvas::getTextInt(TString str, Int_t number){

   Char_t lev[40];



   m_TextLine = new TString(str);

   sprintf(lev,"%i",number);

   m_TextLine->Append(lev);

   return (Char_t*)m_TextLine->Data();



}






/*void DisplayCanvas::setCalDisplay(Float_t cal1, Float_t cal2, Float_t cal3){

   if ( cal1 < cal2 || cal2 < cal3) {

      printf("invalid input numbers\n");

      return;

   } else {

      cal_lev1 = (Float_t)cal1;

      cal_lev2 = (Float_t)cal2;

      cal_lev3 = (Float_t)cal3;

      return;

   }

}*/

---