{ #include #include char str[100]; char str2[100]; for (int r=0; r<20; r++){ std::string s1; if (r<10) s1 = "listPeriod0"; else s1 = "listPeriod"; std::stringstream out; out << r; s1 += out.str(); char* select1; select1 = s1.c_str(); std::string s2; if (r<10) s2 = "./Period0"; else s2 = "./Period"; s2 += out.str(); s2 += ".root"; char* select2; select2 = s2.c_str(); int start, end; TFile *f2 = TFile::Open(select2); int numEntries = bctpc->GetEntriesFast()-1; UInt_t startYear, startMonth, startDay; UInt_t startHour, startMinute, startSecond; UInt_t endYear, endMonth, endDay; UInt_t endHour, endMinute, endSecond; TTimeStamp s,e; bctpc->GetEntry(0); s= bctpc->GetLeaf("timestamp")->GetValue(); int startD = s.GetDate(1,0,&startYear,&startMonth,&startDay); int startT = s.GetTime(1,0,&startHour,&startMinute,&startSecond); bctpc->GetEntry(numEntries); e= bctpc->GetLeaf("timestamp")->GetValue(); int endD = e.GetDate(1,0,&endYear,&endMonth,&endDay); int endT = e.GetTime(1,0,&endHour,&endMinute,&endSecond); f2.Close(); fstream file_op(select1,ios::in); TH1F *hsum_[72]= new TH1F("hsum_0","hsum_0",100,0,2000); int sumEntr = 0; while(!file_op.eof()) { int entries = 0; file_op.getline(str,100); if (strlen(str)==0) break; TFile *f = TFile::Open(str); CosmicCalib_Mod_0->SetDirectory(0); TH1F *h0 = new TH1F(*CosmicCalib_Mod_0); TH1F *h1 = new TH1F(*CosmicCalib_Mod_1); TH1F *h2 = new TH1F(*CosmicCalib_Mod_2); TH1F *h3 = new TH1F(*CosmicCalib_Mod_3); TH1F *h4 = new TH1F(*CosmicCalib_Mod_4); TH1F *h5 = new TH1F(*CosmicCalib_Mod_5); TH1F *h6 = new TH1F(*CosmicCalib_Mod_6); TH1F *h7 = new TH1F(*CosmicCalib_Mod_7); TH1F *h8 = new TH1F(*CosmicCalib_Mod_8); TH1F *h9 = new TH1F(*CosmicCalib_Mod_9); TH1F *h10 = new TH1F(*CosmicCalib_Mod_10); TH1F *h11 = new TH1F(*CosmicCalib_Mod_11); TH1F *h12 = new TH1F(*CosmicCalib_Mod_12); TH1F *h13 = new TH1F(*CosmicCalib_Mod_13); TH1F *h14 = new TH1F(*CosmicCalib_Mod_14); TH1F *h15 = new TH1F(*CosmicCalib_Mod_15); TH1F *h16 = new TH1F(*CosmicCalib_Mod_16); TH1F *h17 = new TH1F(*CosmicCalib_Mod_17); TH1F *h18 = new TH1F(*CosmicCalib_Mod_18); TH1F *h19 = new TH1F(*CosmicCalib_Mod_19); TH1F *h20 = new TH1F(*CosmicCalib_Mod_20); TH1F *h21 = new TH1F(*CosmicCalib_Mod_21); TH1F *h22 = new TH1F(*CosmicCalib_Mod_22); TH1F *h23 = new TH1F(*CosmicCalib_Mod_23); TH1F *h24 = new TH1F(*CosmicCalib_Mod_24); TH1F *h25 = new TH1F(*CosmicCalib_Mod_25); TH1F *h26 = new TH1F(*CosmicCalib_Mod_26); TH1F *h27 = new TH1F(*CosmicCalib_Mod_27); TH1F *h28 = new TH1F(*CosmicCalib_Mod_28); TH1F *h29 = new TH1F(*CosmicCalib_Mod_29); TH1F *h30 = new TH1F(*CosmicCalib_Mod_30); TH1F *h31 = new TH1F(*CosmicCalib_Mod_31); TH1F *h32 = new TH1F(*CosmicCalib_Mod_32); TH1F *h33 = new TH1F(*CosmicCalib_Mod_33); TH1F *h34 = new TH1F(*CosmicCalib_Mod_34); TH1F *h35 = new TH1F(*CosmicCalib_Mod_35); TH1F *h36 = new TH1F(*CosmicCalib_Mod_36); TH1F *h37 = new TH1F(*CosmicCalib_Mod_37); TH1F *h38 = new TH1F(*CosmicCalib_Mod_38); TH1F *h39 = new TH1F(*CosmicCalib_Mod_39); TH1F *h40 = new TH1F(*CosmicCalib_Mod_40); TH1F *h41 = new TH1F(*CosmicCalib_Mod_41); TH1F *h42 = new TH1F(*CosmicCalib_Mod_42); TH1F *h43 = new TH1F(*CosmicCalib_Mod_43); TH1F *h44 = new TH1F(*CosmicCalib_Mod_44); TH1F *h45 = new TH1F(*CosmicCalib_Mod_45); TH1F *h46 = new TH1F(*CosmicCalib_Mod_46); TH1F *h47 = new TH1F(*CosmicCalib_Mod_47); TH1F *h48 = new TH1F(*CosmicCalib_Mod_48); TH1F *h49 = new TH1F(*CosmicCalib_Mod_49); TH1F *h50 = new TH1F(*CosmicCalib_Mod_50); TH1F *h51 = new TH1F(*CosmicCalib_Mod_51); TH1F *h52 = new TH1F(*CosmicCalib_Mod_52); TH1F *h53 = new TH1F(*CosmicCalib_Mod_53); TH1F *h54 = new TH1F(*CosmicCalib_Mod_54); TH1F *h55 = new TH1F(*CosmicCalib_Mod_55); TH1F *h56 = new TH1F(*CosmicCalib_Mod_56); TH1F *h57 = new TH1F(*CosmicCalib_Mod_57); TH1F *h58 = new TH1F(*CosmicCalib_Mod_58); TH1F *h59 = new TH1F(*CosmicCalib_Mod_59); TH1F *h60 = new TH1F(*CosmicCalib_Mod_60); TH1F *h61 = new TH1F(*CosmicCalib_Mod_61); TH1F *h62 = new TH1F(*CosmicCalib_Mod_62); TH1F *h63 = new TH1F(*CosmicCalib_Mod_63); TH1F *h64 = new TH1F(*CosmicCalib_Mod_64); TH1F *h65 = new TH1F(*CosmicCalib_Mod_65); TH1F *h66 = new TH1F(*CosmicCalib_Mod_66); TH1F *h67 = new TH1F(*CosmicCalib_Mod_67); TH1F *h68 = new TH1F(*CosmicCalib_Mod_68); TH1F *h69 = new TH1F(*CosmicCalib_Mod_69); TH1F *h70 = new TH1F(*CosmicCalib_Mod_70); TH1F *h71 = new TH1F(*CosmicCalib_Mod_71); h0->Sumw2(); hsum_0->Add(h0); f->Close(); } hsum_0->Fit("landau","q","",50,1000); hsum_0->Fit("landau","q","",hsum_0->GetFunction("landau")->GetParameter(1)-1.5*hsum_0->GetFunction("landau")->GetParameter(2),hsum_0->GetFunction("landau")->GetParameter(1)+2*hsum_0->GetFunction("landau")->GetParameter(2)); TDatime present; present.Set(); char * nowDate; nowDate = present.AsSQLString(); cout << " BEGIN_TABLE MICROMEGAS_GAIN_TABLE "<< "'" << startYear <<"-"<< setfill('0') << setw(2)<< startMonth << "-" << setfill('0') << setw(2)<< startDay <<" " << setfill('0') << setw(2)<< startHour <<":" << setfill('0') << setw(2)<< startMinute <<":" << setfill('0') << setw(2)<GetFunction("landau")->GetParameter(1)/364.384<<" , "<< hsum_0->GetFunction("landau")->GetParError(1)/364.384<< endl; endl; file_op.close(); } }