L1CADebug.h

Go to the documentation of this file.

#ifndef L1CADebug_h
#define L1CADebug_h 1
 
#define DEBUG_CA
 
#include "CbmL1.h"
 
TH1F  
  *h_pick_res_x[20][3][3], *h_pick_res_y[20][3][3],
  *h_pick_pull_x[20][3][3], *h_pick_pull_y[20][3][3],
  *h_nduplets[20][3], 
  *h_res_x[20][3][3], *h_res_y[20][3][3],
  *h_res_tx[20][3][3],*h_res_ty[20][3][3], *h_res_qp[20][3][3],
  *h_pull_x[20][3][3], *h_pull_y[20][3][3], 
  *h_pull_tx[20][3][3], *h_pull_ty[20][3][3], *h_pull_qp[20][3][3]
  ;
//TH2F *h_dyvsy[20], *h_dxvsx[20];
 
void Pulls( int i, int j, int k, 
      double *mc, L1TrackPar &T, fvec qp0, L1FieldRegion &fld )
{
  L1TrackPar tmp_T = T;
  fvec z = mc[5];
  L1Extrapolate( tmp_T, z, qp0, fld );
  double dx = mc[0] - tmp_T.x[0];
  double dy = mc[1] - tmp_T.y[0];
  double dtx = mc[2] - tmp_T.tx[0];
  double dty = mc[3] - tmp_T.ty[0];
  double dqp = mc[4] - tmp_T.qp[0];
  double sx = sqrt(tmp_T.C00[0]);
  double sy = sqrt(tmp_T.C11[0]);
  double stx = sqrt(tmp_T.C22[0]);
  double sty = sqrt(tmp_T.C33[0]);
  double sqp = sqrt(tmp_T.C44[0]);
  h_res_x  [i][j][k]->Fill(dx);
  h_res_y  [i][j][k]->Fill(dy);       
  h_res_tx [i][j][k]->Fill(dtx);
  h_res_ty [i][j][k]->Fill(dty);        
  h_res_qp [i][j][k]->Fill(dqp/mc[4]);        
  h_pull_x [i][j][k]->Fill(dx/sx);
  h_pull_y [i][j][k]->Fill(dy/sy);        
  h_pull_tx[i][j][k]->Fill(dtx/stx);
  h_pull_ty[i][j][k]->Fill(dty/sty);        
  h_pull_qp[i][j][k]->Fill(dqp/sqp);        
}
 
void InitL1Histo(){
  static bool first_call = 1;
  if( first_call ){
    TDirectory *curdir = gDirectory;
    CbmL1::Instance()->histodir->cd();
    L1Algo *algo = CbmL1::Instance()->algo;
    gDirectory->mkdir("L1ALGO");
    gDirectory->cd("L1ALGO");
    for( int ist=0; ist<algo->NStations; ist++ ){      
      char stname[225];//, namedir[225], title[225];
      sprintf(stname,"%i",ist);
      gDirectory->mkdir(stname);  
      gDirectory->cd(stname);
      {
  //h_dyvsy[ist] = new TH2F("h_dyvsy","h_dyvsy",1000,0,100, 1000, 0.,.5);
  //h_dxvsx[ist] = new TH2F("h_dxvsx","h_dxvsx",1000,0,100, 1000, 0.,2.5);
  const char* stepname[3]={"D0","Ref","Sec"};
  for( int istep=0; istep<3; istep++){
    gDirectory->mkdir(stepname[istep]);  
    gDirectory->cd(stepname[istep]);
 
    h_nduplets[ist][istep] = new TH1F("h_nduplets","h_nduplets",50,0,5000);
 
    const char* lmr[3]={"left","middle","right"};
    for( int ilmr=0; ilmr<3; ilmr++){
      gDirectory->mkdir(lmr[ilmr]);  
      gDirectory->cd(lmr[ilmr]);
      {
        h_pick_res_x[ist][ilmr][istep] = new TH1F("h_pick_res_x","h_pick_res_x",100,-2,3);
        h_pick_res_y[ist][ilmr][istep] = new TH1F("h_pick_res_y","h_pick_res_y",100,-1,1);
        h_pick_pull_x[ist][ilmr][istep] = new TH1F("h_pick_pull_x","h_pick_pull_x",100,-10,10);
        h_pick_pull_y[ist][ilmr][istep] = new TH1F("h_pick_pull_y","h_pick_pull_y",100,-10,10);
        h_res_x[ist][ilmr][istep] = new TH1F("h_res_x","h_res_x",100,-.01,.01);
        h_res_y[ist][ilmr][istep] = new TH1F("h_res_y","h_res_y",100,-.1,.1);
        h_res_tx[ist][ilmr][istep] = new TH1F("h_res_tx","h_res_tx",100,-.01,.01);
        h_res_ty[ist][ilmr][istep] = new TH1F("h_res_ty","h_res_ty",100,-.01,.01);
        h_res_qp[ist][ilmr][istep] = new TH1F("h_res_qp","h_res_qp",100,-.1,.1);
        h_pull_x[ist][ilmr][istep] = new TH1F("h_pull_x","h_pull_x",100,-10,10);
        h_pull_y[ist][ilmr][istep] = new TH1F("h_pull_y","h_pull_y",100,-10,10);
        h_pull_tx[ist][ilmr][istep] = new TH1F("h_pull_tx","h_pull_tx",100,-10,10);
        h_pull_ty[ist][ilmr][istep] = new TH1F("h_pull_ty","h_pull_ty",100,-10,10);
        h_pull_qp[ist][ilmr][istep] = new TH1F("h_pull_qp","h_pull_qp",100,-10,10);
        gDirectory->cd("..");
      }
    } gDirectory->cd("..");
  }
      } gDirectory->cd("..");
    } 
    gDirectory->cd("..");    
    curdir->cd();
    first_call = 0;
  }
}
 
 
#endif