FillValidRecMdcTrackHisto.c

Go to the documentation of this file.

void FillValidRecMdcTrackHisto(char *filename="ValidRecMdcTrackNTuple_new.root"){
  TFile *fin = new TFile(filename,"READ");
  TTree *mdc = (TTree *)fin->Get("mdcTrack");
 
  TFile *fin1 = new TFile(filename,"READ");
  TTree *mdc1 = (TTree *)fin1->Get("mdcHit");
 
  int   nhits;
  float Nhits,Nster;
  float eventNo,Cpa, runNo, Id, Phi0;// kappa;
  float Tanl, Dr, Dz, Chisq, FiTerm;
  float Hits_ddl,Hits_id;
  float Hits_ddr, doca, Hits_mdc_id, layer, wire;
  float Hits_tdc, Hits_adc;// energy, nhits, path1;
 
 
  mdc->SetBranchAddress("eventNo",   &eventNo);
  mdc->SetBranchAddress("runNo",&runNo);
  mdc->SetBranchAddress("nhits",&nhits);
  mdc->SetBranchAddress("Id",    &Id);
  mdc->SetBranchAddress("Phi0", &Phi0);
  mdc->SetBranchAddress("Tanl",&Tanl);
  mdc->SetBranchAddress("Cpa",&Cpa);
  mdc->SetBranchAddress("Dr", &Dr);
  mdc->SetBranchAddress("Dz",    &Dz);
  mdc->SetBranchAddress("Chisq",    &Chisq);
  mdc->SetBranchAddress("FiTerm",    &FiTerm);
  mdc1->SetBranchAddress("Hits_id",    &Hits_id);
  mdc1->SetBranchAddress("Hits_ddl",    &Hits_ddl);
  mdc1->SetBranchAddress("Hits_ddr",    &Hits_ddr);
  mdc->SetBranchAddress("Nhits",    &Nhits);
  mdc->SetBranchAddress("Nster",    &Nster);
  mdc1->SetBranchAddress("doca",    &doca);
  mdc1->SetBranchAddress("Hits_mdc_id",    &Hits_mdc_id);
  mdc1->SetBranchAddress("layer",    &layer);
  mdc1->SetBranchAddress("wire",    &wire);
  mdc1->SetBranchAddress("Hits_tdc",    &Hits_tdc);
  mdc1->SetBranchAddress("Hits_adc",    &Hits_adc);
//  mdc->SetBranchAddress("energy",    &energy);
//  mdc->SetBranchAddress("nhits",    &nhits);
//  mdc->SetBranchAddress("path1",    &path1);
 
 
 
  TH1F *m_nhits                = new TH1F("nhits","nhits",100,0,5000);
  TH1F *m_eventNo              = new TH1F("eventNo","eventNo", 100,0,10000);
  TH1F *m_runNo                = new TH1F("runNo","runNo",100, 0, 100000);
  TH1F *m_Id                   = new TH1F("Id","track_id",100,0,100);
  TH1F *m_Phi0                 = new TH1F("Phi0","phi0",100,0,3.14);
  TH1F *m_Tanl                 = new TH1F("Tanl","tanl",100,-300,300);
  TH1F *m_Dr                   = new TH1F("Dr","dr",100,-300,300);
  TH1F *m_Dz                   = new TH1F("Dz","dz",100,-300,300);
  TH1F *m_Chisq                = new TH1F("Chisq","chisq_of_fit",100,-300,3000);
  TH1F *m_FiTerm               = new TH1F("FiTerm","phi_terminal",100,0,6.3);
  TH1F *m_Nhits                = new TH1F("Nhits","number_of_hits",100,0,500);
  TH1F *m_Nster                = new TH1F("Nster","number_of_stereo_hits",100,0,5000);
  TH1F *m_Cpa                  = new TH1F("Cpa","cpa",100,0,5000);
 
  TFile *fout1   = new TFile("ValidRecMdcTrackHist.root","recreate");
 
  TH1F *m_Hits_id              = new TH1F("Hits_id","hits_id",100,0,300);
  TH1F *m_Hits_ddl             = new TH1F("Hits_ddl","hits_ddl",100,-300,300);
  TH1F *m_Hits_ddr             = new TH1F("Hits_ddr","hits_ddr",100,-300,300);
  TH1F *m_doca                 = new TH1F("doca","doca",100,-300,300);
  TH1F *m_Hits_mdc_id          = new TH1F("Hits_mdc_id","hits_mdc_id",100,-300,300);
  TH1F *m_layer                = new TH1F("layer","layer",100,-300,300);  
  TH1F *m_wire                 = new TH1F("wire","wire",100,-300,300);
  TH1F *m_Hits_tdc             = new TH1F("Hits_tdc","hits_tdc",100,-300,300);
  TH1F *m_Hits_adc             = new TH1F("Hits_adc","hits_adc",100,-300,300);
 
 
  for(int i = 0; i < mdc->GetEntries(); i++){
    mdc->GetEntry(i);
 
    m_nhits->Fill(nhits);
    m_eventNo->Fill(eventNo);
    m_runNo->Fill(runNo);
    m_Id->Fill(Id);
    m_Phi0->Fill(Phi0);
    m_Cpa->Fill(Cpa);
    m_Tanl->Fill(Tanl);
 
    m_Dr->Fill(Dr);
    m_Dz->Fill(Dz);
    m_Chisq->Fill(Chisq);
    m_FiTerm->Fill(FiTerm);
    m_Nhits->Fill(Nhits);
    m_Nster->Fill(Nster);
/*
    m_Hits_id->Fill(Hits_id);
    m_Hits_ddl->Fill(Hits_ddl);
    m_Hits_ddr->Fill(Hits_ddr);
    m_doca->Fill(doca);
    m_Hits_mdc_id->Fill(Hits_mdc_id);
    m_layer->Fill(layer);
    m_wire->Fill(wire);
    m_Hits_tdc->Fill(Hits_tdc);
    m_Hits_adc->Fill(Hits_adc);
*/   
 
  }
  for(int i = 0; i < mdc1->GetEntries(); i++){
    mdc1->GetEntry(i);
 
    m_Hits_id->Fill(Hits_id);
    m_Hits_ddl->Fill(Hits_ddl);
    m_Hits_ddr->Fill(Hits_ddr);
    m_doca->Fill(doca);
    m_Hits_mdc_id->Fill(Hits_mdc_id);
    m_layer->Fill(layer);
    m_wire->Fill(wire);
    m_Hits_tdc->Fill(Hits_tdc);
    m_Hits_adc->Fill(Hits_adc);
  }
  m_nhits->Write();
  m_eventNo->Write();
  m_runNo->Write();
  m_Id->Write();
  m_Phi0->Write();
  m_Cpa->Write();
 
  m_Tanl->Write();
  m_Dr->Write();
  m_Dz->Write();
  m_Chisq->Write();
  m_FiTerm->Write();
  m_Nhits->Write();
  m_Nster->Write();
  m_Hits_id->Write();
  m_Hits_ddl->Write();
  m_Hits_ddr->Write();
  m_doca->Write();
  m_Hits_mdc_id->Write();
  m_layer->Write();
  m_wire->Write();
  m_Hits_tdc->Write();
  m_Hits_adc->Write();
 
  fin->Close();
   fin1->Close();
  fout1->Close();
}