00001 void FillMdcDedxHisto(char *filename="/ihepbatch/besd13/xcao/6.3.3/valid/rec/NTuple/MdcDedxValidNTuple_e_1Gev.root"){
00002
00003 TFile *fin = new TFile(filename,"READ");
00004 TTree *Mdc_dedx = (TTree *)fin->Get("Mdc_dedx");
00005
00006 double dedx, parttype, chidedx, chidedxE, chidedxMu,chidedxPi,chidedxK,chidedxP;
00007
00008 Mdc_dedx->SetBranchAddress("dedx", &dedx);
00009 Mdc_dedx->SetBranchAddress("parttype",&parttype);
00010 Mdc_dedx->SetBranchAddress("chidedx",&chidedx);
00011 Mdc_dedx->SetBranchAddress("chidedxE",&chidedxE);
00012 Mdc_dedx->SetBranchAddress("chidedxMu",&chidedxMu);
00013 Mdc_dedx->SetBranchAddress("chidedxPi",&chidedxPi);
00014 Mdc_dedx->SetBranchAddress("chidedxK",&chidedxK);
00015 Mdc_dedx->SetBranchAddress("chidedxP",&chidedxP);
00016
00017 TFile *fout = new TFile("/ihepbatch/besd13/xcao/6.3.3/valid/rec/Histgram_10/MdcDedxValidHist_e_1Gev.root","recreate");
00018
00019 TH1F *Dedx = new TH1F("Dedx","phase shape",200,0,4000);
00020 TH1F *PartType = new TH1F("PartType","particle type", 200, 0,5);
00021 TH1F *ChiDedx = new TH1F("ChiDedx","chi value of most probability particle",200, -10, 10);
00022 TH1F *ChiDedxE = new TH1F("ChiDedxE","chi value of Electron ",200, -10, 10);
00023 TH1F *ChiDedxMu = new TH1F("ChiDedxMu","chi value of Muon",200, -10, 10);
00024 TH1F *ChiDedxPi = new TH1F("ChiDedxPi","chi value of Pion",200, -10, 10);
00025 TH1F *ChiDedxK = new TH1F("ChiDedxK","chi value of kaon",200, -10, 10);
00026 TH1F *ChiDedxP = new TH1F("ChiDedxP","chi value of Proton",200, -10, 10);
00027
00028
00029 for(int i = 0; i < Mdc_dedx->GetEntries(); i++){
00030
00031 Mdc_dedx->GetEntry(i);
00032 Dedx->Fill(dedx);
00033 PartType->Fill(parttype);
00034 ChiDedx->Fill(chidedx);
00035 ChiDedxE->Fill(chidedxE);
00036 ChiDedxMu->Fill(chidedxMu);
00037 ChiDedxPi->Fill(chidedxPi);
00038 ChiDedxK->Fill(chidedxK);
00039 ChiDedxP->Fill(chidedxP);
00040
00041 }
00042
00043 Dedx->Write();
00044 PartType->Write();
00045 ChiDedx->Write();
00046 ChiDedxE->Write();
00047 ChiDedxMu->Write();
00048 ChiDedxPi->Write();
00049 ChiDedxK->Write();
00050 ChiDedxP->Write();
00051
00052 fin->Close();
00053 fout->Close();
00054
00055 }
