00001 void FillValidRecMdcTrackHisto(char *filename="ValidRecMdcTrackNTuple_new.root"){
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;
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;
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
00040
00041
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
00091
00092
00093
00094
00095
00096
00097
00098
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 }