Go to the source code of this file.
Functions | |
void | FillValidRecMdcTrackHisto (char *filename="ValidRecMdcTrackNTuple_new.root") |
void FillValidRecMdcTrackHisto | ( | char * | filename = "ValidRecMdcTrackNTuple_new.root" |
) |
Definition at line 1 of file FillValidRecMdcTrackHisto.c.
References eventNo, check_raw_filter::filename, genRecEmupikp::i, m_doca, m_layer, m_wire, mdc, and runNo.
00001 { 00002 TFile *fin = new TFile(filename,"READ"); 00003 TTree *mdc = (TTree *)fin->Get("mdcTrack"); 00004 00005 TFile *fin1 = new TFile(filename,"READ"); 00006 TTree *mdc1 = (TTree *)fin1->Get("mdcHit"); 00007 00008 int nhits; 00009 float Nhits,Nster; 00010 float eventNo,Cpa, runNo, Id, Phi0;// kappa; 00011 float Tanl, Dr, Dz, Chisq, FiTerm; 00012 float Hits_ddl,Hits_id; 00013 float Hits_ddr, doca, Hits_mdc_id, layer, wire; 00014 float Hits_tdc, Hits_adc;// energy, nhits, path1; 00015 00016 00017 mdc->SetBranchAddress("eventNo", &eventNo); 00018 mdc->SetBranchAddress("runNo",&runNo); 00019 mdc->SetBranchAddress("nhits",&nhits); 00020 mdc->SetBranchAddress("Id", &Id); 00021 mdc->SetBranchAddress("Phi0", &Phi0); 00022 mdc->SetBranchAddress("Tanl",&Tanl); 00023 mdc->SetBranchAddress("Cpa",&Cpa); 00024 mdc->SetBranchAddress("Dr", &Dr); 00025 mdc->SetBranchAddress("Dz", &Dz); 00026 mdc->SetBranchAddress("Chisq", &Chisq); 00027 mdc->SetBranchAddress("FiTerm", &FiTerm); 00028 mdc1->SetBranchAddress("Hits_id", &Hits_id); 00029 mdc1->SetBranchAddress("Hits_ddl", &Hits_ddl); 00030 mdc1->SetBranchAddress("Hits_ddr", &Hits_ddr); 00031 mdc->SetBranchAddress("Nhits", &Nhits); 00032 mdc->SetBranchAddress("Nster", &Nster); 00033 mdc1->SetBranchAddress("doca", &doca); 00034 mdc1->SetBranchAddress("Hits_mdc_id", &Hits_mdc_id); 00035 mdc1->SetBranchAddress("layer", &layer); 00036 mdc1->SetBranchAddress("wire", &wire); 00037 mdc1->SetBranchAddress("Hits_tdc", &Hits_tdc); 00038 mdc1->SetBranchAddress("Hits_adc", &Hits_adc); 00039 // mdc->SetBranchAddress("energy", &energy); 00040 // mdc->SetBranchAddress("nhits", &nhits); 00041 // mdc->SetBranchAddress("path1", &path1); 00042 00043 00044 00045 TH1F *m_nhits = new TH1F("nhits","nhits",100,0,5000); 00046 TH1F *m_eventNo = new TH1F("eventNo","eventNo", 100,0,10000); 00047 TH1F *m_runNo = new TH1F("runNo","runNo",100, 0, 100000); 00048 TH1F *m_Id = new TH1F("Id","track_id",100,0,100); 00049 TH1F *m_Phi0 = new TH1F("Phi0","phi0",100,0,3.14); 00050 TH1F *m_Tanl = new TH1F("Tanl","tanl",100,-300,300); 00051 TH1F *m_Dr = new TH1F("Dr","dr",100,-300,300); 00052 TH1F *m_Dz = new TH1F("Dz","dz",100,-300,300); 00053 TH1F *m_Chisq = new TH1F("Chisq","chisq_of_fit",100,-300,3000); 00054 TH1F *m_FiTerm = new TH1F("FiTerm","phi_terminal",100,0,6.3); 00055 TH1F *m_Nhits = new TH1F("Nhits","number_of_hits",100,0,500); 00056 TH1F *m_Nster = new TH1F("Nster","number_of_stereo_hits",100,0,5000); 00057 TH1F *m_Cpa = new TH1F("Cpa","cpa",100,0,5000); 00058 00059 TFile *fout1 = new TFile("ValidRecMdcTrackHist.root","recreate"); 00060 00061 TH1F *m_Hits_id = new TH1F("Hits_id","hits_id",100,0,300); 00062 TH1F *m_Hits_ddl = new TH1F("Hits_ddl","hits_ddl",100,-300,300); 00063 TH1F *m_Hits_ddr = new TH1F("Hits_ddr","hits_ddr",100,-300,300); 00064 TH1F *m_doca = new TH1F("doca","doca",100,-300,300); 00065 TH1F *m_Hits_mdc_id = new TH1F("Hits_mdc_id","hits_mdc_id",100,-300,300); 00066 TH1F *m_layer = new TH1F("layer","layer",100,-300,300); 00067 TH1F *m_wire = new TH1F("wire","wire",100,-300,300); 00068 TH1F *m_Hits_tdc = new TH1F("Hits_tdc","hits_tdc",100,-300,300); 00069 TH1F *m_Hits_adc = new TH1F("Hits_adc","hits_adc",100,-300,300); 00070 00071 00072 for(int i = 0; i < mdc->GetEntries(); i++){ 00073 mdc->GetEntry(i); 00074 00075 m_nhits->Fill(nhits); 00076 m_eventNo->Fill(eventNo); 00077 m_runNo->Fill(runNo); 00078 m_Id->Fill(Id); 00079 m_Phi0->Fill(Phi0); 00080 m_Cpa->Fill(Cpa); 00081 m_Tanl->Fill(Tanl); 00082 00083 m_Dr->Fill(Dr); 00084 m_Dz->Fill(Dz); 00085 m_Chisq->Fill(Chisq); 00086 m_FiTerm->Fill(FiTerm); 00087 m_Nhits->Fill(Nhits); 00088 m_Nster->Fill(Nster); 00089 /* 00090 m_Hits_id->Fill(Hits_id); 00091 m_Hits_ddl->Fill(Hits_ddl); 00092 m_Hits_ddr->Fill(Hits_ddr); 00093 m_doca->Fill(doca); 00094 m_Hits_mdc_id->Fill(Hits_mdc_id); 00095 m_layer->Fill(layer); 00096 m_wire->Fill(wire); 00097 m_Hits_tdc->Fill(Hits_tdc); 00098 m_Hits_adc->Fill(Hits_adc); 00099 */ 00100 00101 } 00102 for(int i = 0; i < mdc1->GetEntries(); i++){ 00103 mdc1->GetEntry(i); 00104 00105 m_Hits_id->Fill(Hits_id); 00106 m_Hits_ddl->Fill(Hits_ddl); 00107 m_Hits_ddr->Fill(Hits_ddr); 00108 m_doca->Fill(doca); 00109 m_Hits_mdc_id->Fill(Hits_mdc_id); 00110 m_layer->Fill(layer); 00111 m_wire->Fill(wire); 00112 m_Hits_tdc->Fill(Hits_tdc); 00113 m_Hits_adc->Fill(Hits_adc); 00114 } 00115 m_nhits->Write(); 00116 m_eventNo->Write(); 00117 m_runNo->Write(); 00118 m_Id->Write(); 00119 m_Phi0->Write(); 00120 m_Cpa->Write(); 00121 00122 m_Tanl->Write(); 00123 m_Dr->Write(); 00124 m_Dz->Write(); 00125 m_Chisq->Write(); 00126 m_FiTerm->Write(); 00127 m_Nhits->Write(); 00128 m_Nster->Write(); 00129 m_Hits_id->Write(); 00130 m_Hits_ddl->Write(); 00131 m_Hits_ddr->Write(); 00132 m_doca->Write(); 00133 m_Hits_mdc_id->Write(); 00134 m_layer->Write(); 00135 m_wire->Write(); 00136 m_Hits_tdc->Write(); 00137 m_Hits_adc->Write(); 00138 00139 fin->Close(); 00140 fin1->Close(); 00141 fout1->Close(); 00142 }