/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Validation/ValidRecMdcKalTrack/ValidRecMdcKalTrack-00-00-07/share/FillMdcKalTrackHist.C File Reference

#include <cmath>

Go to the source code of this file.

Functions
void	FillMdcKalTrackHist (void)

---

Function Documentation

void FillMdcKalTrackHist

(

void

)

Definition at line 4 of file FillMdcKalTrackHist.C.

References check_raw_filter::filename, and genRecEmupikp::i.

                              {

  const char* const filename = "MdcKalTrackValidNTuple.root";

  TFile *fin = new TFile(filename);
  
  TTree *r1 = (TTree *)fin->Get("n101");
  //TTree *r3 = (TTree *)fin->Get("n103");
  
  Double_t zhelixe[5];
  Double_t zhelixmu[5];
  Double_t zhelix[5];
  Double_t zhelixk[5];
  Double_t zhelixp[5];

  Double_t fhelixe[5];
  Double_t fhelixmu[5];
  Double_t fhelix[5];
  Double_t fhelixk[5];
  Double_t fhelixp[5];
   
  Double_t helix_rec[5];
  Double_t stat[5][2];
 
  r1->SetBranchAddress("zhelixe", &zhelixe);
  r1->SetBranchAddress("zhelixmu", &zhelixmu);
  r1->SetBranchAddress("zhelix", &zhelix);
  r1->SetBranchAddress("zhelixk", &zhelixk);
  r1->SetBranchAddress("zhelixp", &zhelixp);
 
  r1->SetBranchAddress("fhelixe", &fhelixe);
  r1->SetBranchAddress("fhelixmu", &fhelixmu);
  r1->SetBranchAddress("fhelix", &fhelix);
  r1->SetBranchAddress("fhelixk", &fhelixk);
  r1->SetBranchAddress("fhelixp", &fhelixp);
  
  r1->SetBranchAddress("stat", &stat);
  //r3->SetBranchAddress("trkhelix", &helix_rec);
 
 
  TFile *fout = new TFile("MdcKalTrackValidHist.root","recreate");

  TH1F *zdr_e  = new TH1F("zdr_e", "zdr_e of zhelixe",100,-1,1);
  TH1F *zdr_mu = new TH1F("zdr_mu","zdr_mu of zhelixe",100,-1,1);
  TH1F *zdr_pi = new TH1F("zdr_pi","zdr_pi of zhelixe",100,-1,1);
  TH1F *zdr_k  = new TH1F("zdr_k", "zdr_k of zhelixe",100,-1,1);
  TH1F *zdr_p  = new TH1F("zdr_p", "zdr_p of zhelixe",100,-1,1);
 
  TH1F *zdz_e  = new TH1F("zdz_e", "zdz_e of zhelixe",100,-10,10);
  TH1F *zdz_mu = new TH1F("zdz_mu","zdz_mu of zhelixe",100,-10,10);
  TH1F *zdz_pi = new TH1F("zdz_pi","zdz_pi of zhelixe",100,-10,10);
  TH1F *zdz_k  = new TH1F("zdz_k", "zdz_k of zhelixe",100,-10,10);
  TH1F *zdz_p  = new TH1F("zdz_p", "zdz_p of zhelixe",100,-10,10);
        

  TH1F *zpt_e  = new TH1F("zpt_e", "zpt of e",100,0,1.5);
  TH1F *zpt_mu = new TH1F("zpt_mu","zpt of mu",100,0,1.5);
  TH1F *zpt_pi = new TH1F("zpt_pi","zpt of pi",100,0,1.5);
  TH1F *zpt_k  = new TH1F("zpt_k", "zpt of k",100,0,1.5);
  TH1F *zpt_p  = new TH1F("zpt_p", "zpt of p",100,0,1.5);

  TH1F *zptot_e  = new TH1F("zptot_e", "zptot of e",100,0,2);
  TH1F *zptot_mu = new TH1F("zptot_mu","zptot of mu",100,0,2);
  TH1F *zptot_pi = new TH1F("zptot_pi","zptot of pi",100,0,2);
  TH1F *zptot_k  = new TH1F("zptot_k", "zptot of k",100,0,2);
  TH1F *zptot_p  = new TH1F("zptot_p", "zptot of p",100,0,2);

  TH1F *dr_rec   = new TH1F("dr_rec", "zpt of e",100,-1,1);
  TH1F *dz_rec   = new TH1F("dz_rec", "zpt of e",100,-10,10);
  TH1F *pt_rec   = new TH1F("pt_rec", "zpt of e",100,0,1.5);
  TH1F *ptot_rec = new TH1F("ptot_rec", "zpt of e",100,0,2);


  for(Int_t i=0; i<r1->GetEntries(); i++){

    r1->GetEntry(i);
 
    if(stat[0][0]==0) {
      zdr_e->Fill(zhelixe[0]);
      zdz_e->Fill(zhelixe[3]);
      zpt_e->Fill(1./fabs(zhelixe[2]));
      zptot_e->Fill(sqrt(1+zhelixe[4]*zhelixe[4])/fabs(zhelixe[2]));
    }
 

    if(stat[1][0]==0) {
      zdr_mu->Fill(zhelixmu[0]);
      zdz_mu->Fill(zhelixmu[3]);
      zpt_mu->Fill(1./fabs(zhelixmu[2]));
      zptot_mu->Fill(sqrt(1+zhelixmu[4]*zhelixmu[4])/fabs(zhelixmu[2]));
    }

    if(stat[2][0]==0) {
      zdr_pi->Fill(zhelix[0]);
      zdz_pi->Fill(zhelix[3]);
      zpt_pi->Fill(1./fabs(zhelix[2]));
      zptot_pi->Fill(sqrt(1+zhelix[4]*zhelix[4])/fabs(zhelix[2]));
    }

    if(stat[3][0]==0) {
      zdr_k->Fill(zhelixk[0]);
      zdz_k->Fill(zhelixk[3]);
      zpt_k->Fill(1./fabs(zhelixk[2]));
      zptot_k->Fill(sqrt(1+zhelixk[4]*zhelixk[4])/fabs(zhelixk[2]));
    }

    if(stat[4][0]==0) {
      zdr_p->Fill(zhelixp[0]);
      zdz_p->Fill(zhelixp[3]);
      zpt_p->Fill(1./fabs(zhelixp[2]));
      zptot_p->Fill(sqrt(1+zhelixp[4]*zhelixp[4])/fabs(zhelixp[2]));
    }
  } 
 
/*
  for(Int_t j=0; j<r3->GetEntries(); j++){

    dr_rec->Fill(helix_rec[0]);
    dz_rec->Fill(helix_rec[3]);
    pt_rec->Fill(1./fabs(helix_rec[2]));
    ptot_rec->Fill(sqrt(1+helix_rec[4]*helix_rec[4])/fabs(helix_rec[2]));
  }
 
  */

  zdr_e->Write();
  zdr_mu->Write();
  zdr_pi->Write();
  zdr_k->Write();
  zdr_p->Write();
 
  zdz_e->Write();
  zdz_mu->Write();
  zdz_pi->Write();
  zdz_k->Write();
  zdz_p->Write();

  zpt_e->Write();
  zpt_mu->Write();
  zpt_pi->Write();
  zdr_k->Write();
  zdr_p->Write();

  zptot_e->Write();
  zptot_mu->Write();
  zptot_pi->Write();
  zptot_k->Write();
  zptot_p->Write();
 
  fin->Close();
  fout->Close();
}

---