00001 void FillExtHisto(char *filename="ExtValidNTuple.root"){ 00002 00003 TFile *fin = new TFile(filename,"READ"); 00004 TTree *T = (TTree *)fin->Get("ExtTrks"); 00005 00006 double depth, distance, deltaPhi, chi2; 00007 double maxHits, numLayers, numHits; 00008 double xPos, yPos, zPos; 00009 00010 double tof1, tof1path, tof1PosSigX; 00011 double tof1PosX, tof1PosY, tof1PosZ; 00012 double emTof1_11, emTof1_33, emTof1_55; 00013 double emTof1_66, emTof1_23, emTof1_35; 00014 00015 T->SetBranchAddress("tof1", &tof1); 00016 T->SetBranchAddress("tof1Path", &tof1path); 00017 T->SetBranchAddress("tof1PosSigmaAlongX", &tof1PosSigX); 00018 T->SetBranchAddress("xTof1Pos", &tof1PosX); 00019 T->SetBranchAddress("yTof1Pos", &tof1PosY); 00020 T->SetBranchAddress("zTof1Pos", &tof1PosZ); 00021 T->SetBranchAddress("emTof1_11", &emTof1_11); 00022 T->SetBranchAddress("emTof1_33", &emTof1_33); 00023 T->SetBranchAddress("emTof1_55", &emTof1_55); 00024 T->SetBranchAddress("emTof1_66", &emTof1_66); 00025 T->SetBranchAddress("emTof1_23", &emTof1_23); 00026 T->SetBranchAddress("emTof1_35", &emTof1_35); 00027 00028 TFile *fout = new TFile("ExtValidHist.root","recreate"); 00029 00030 TH1F *H_tof1 = new TH1F("H_tof1","tof1",100,0,1000); 00031 TH1F *H_tof1path = new TH1F("H_tof1path","tof1path",100,0,300); 00032 TH1F *H_tof1PosSigX = new TH1F("H_tof1PosSigX","tof1PosSigX",100,0,100); 00033 TH1F *H_tof1PosX = new TH1F("H_tof1PosX","tof1PosX",100,-200,200); 00034 TH1F *H_tof1PosY = new TH1F("H_tof1PosY","tof1PosY",100,-200,200); 00035 TH1F *H_tof1PosZ = new TH1F("H_tof1PosZ","tof1PosZ",100,-400,400); 00036 TH1F *H_emTof1_11 = new TH1F("H_emTof1_11","emTof1_11",100,-100,100); 00037 TH1F *H_emTof1_33 = new TH1F("H_emTof1_33","emTof1_33",100,-100,100); 00038 TH1F *H_emTof1_55 = new TH1F("H_emTof1_55","emTof1_55",100,-100,100); 00039 TH1F *H_emTof1_66 = new TH1F("H_emTof1_66","emTof1_66",100,-100,100); 00040 TH1F *H_emTof1_23 = new TH1F("H_emTof1_23","emTof1_23",100,-100,100); 00041 TH1F *H_emTof1_35 = new TH1F("H_emTof1_35","emTof1_35",100,-100,100); 00042 00043 00044 for(int i = 0; i < T->GetEntries(); i++){ 00045 00046 T->GetEntry(i); 00047 00048 H_tof1->Fill(tof1); 00049 H_tof1path->Fill(tof1path); 00050 H_tof1PosSigX->Fill(tof1PosSigX); 00051 H_tof1PosX->Fill(tof1PosX); 00052 H_tof1PosY->Fill(tof1PosY); 00053 H_tof1PosZ->Fill(tof1PosZ); 00054 H_emTof1_11->Fill(emTof1_11); 00055 H_emTof1_33->Fill(emTof1_33); 00056 H_emTof1_55->Fill(emTof1_55); 00057 H_emTof1_66->Fill(emTof1_66); 00058 H_emTof1_23->Fill(emTof1_23); 00059 H_emTof1_35->Fill(emTof1_35); 00060 00061 } 00062 00063 H_tof1->Write(); 00064 H_tof1path->Write(); 00065 H_tof1PosSigX->Write(); 00066 H_tof1PosX->Write(); 00067 H_tof1PosY->Write(); 00068 H_tof1PosZ->Write(); 00069 H_emTof1_11->Write(); 00070 H_emTof1_33->Write(); 00071 H_emTof1_55->Write(); 00072 H_emTof1_66->Write(); 00073 H_emTof1_23->Write(); 00074 H_emTof1_35->Write(); 00075 00076 fin->Close(); 00077 fout->Close(); 00078 00079 }