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