{ TH1F *hPtJets3 = new TH1F("hPtJets3", "hPtJets3", 50, 0, 150); TH1F *hEtaJets3 = new TH1F("hEtaJets3", "hEtaJets3", 50, -5.5, 5.5); TH1F *hPtJets5 = new TH1F("hPtJets5", "hPtJets5", 50, 0, 150); TH1F *hEtaJets5 = new TH1F("hEtaJets5", "hEtaJets5", 50, -5.5, 5.5); TH1F *hPtJets7 = new TH1F("hPtJets7", "hPtJets7", 50, 0, 150); TH1F *hEtaJets7 = new TH1F("hEtaJets7", "hEtaJets7", 50, -5.5, 5.5); TH1F *hPtJets3_eta1 = new TH1F("hPtJets3_eta1", "hPtJets3_eta1", 50, 0, 150); TH1F *hPtJets3_eta2 = new TH1F("hPtJets3_eta2", "hPtJets3_eta2", 50, 0, 150); TH1F *hPtJets5_eta1 = new TH1F("hPtJets5_eta1", "hPtJets5_eta1", 50, 0, 150); TH1F *hPtJets5_eta2 = new TH1F("hPtJets5_eta2", "hPtJets5_eta2", 50, 0, 150); TH1F *hPtJets7_eta1 = new TH1F("hPtJets7_eta1", "hPtJets7_eta1", 50, 0, 150); TH1F *hPtJets7_eta2 = new TH1F("hPtJets7_eta2", "hPtJets7_eta2", 50, 0, 150); TH1F *x1_3 = new TH1F("x1_3", "x1_3", 50, -6, 0); TH1F *x1_5 = new TH1F("x1_5", "x1_5", 50, -6, 0); TH1F *x1_7 = new TH1F("x1_7", "x1_7", 50, -6, 0); TH1F *x2_3 = new TH1F("x2_3", "x2_3", 50, -6, 0); TH1F *x2_5 = new TH1F("x2_5", "x2_5", 50, -6, 0); TH1F *x2_7 = new TH1F("x2_7", "x2_7", 50, -6, 0); TH1F *xf_3 = new TH1F("xf_3", "xf_3", 50, -8, 0); TH1F *xf_5 = new TH1F("xf_5", "xf_5", 50, -8, 0); TH1F *xf_7 = new TH1F("xf_7", "xf_7", 50, -8, 0); TH1F *Q2_3 = new TH1F("Q2_3", "Q2_3", 50, 0, 500); TH1F *Q2_5 = new TH1F("Q2_5", "Q2_5", 50, 0, 500); TH1F *Q2_7 = new TH1F("Q2_7", "Q2_7", 50, 0, 500); /* TH1F *x1_3_eta1 = new TH1F("x1_3_eta1", "x1_3_eta1", 50, -6, 0); TH1F *x1_3_eta2 = new TH1F("x1_3_eta2", "x1_3_eta2", 50, -6, 0); TH1F *x1_5_eta1 = new TH1F("x1_5_eta1", "x1_5_eta1", 50, -6, 0); TH1F *x1_5_eta2 = new TH1F("x1_5_eta2", "x1_5_eta2", 50, -6, 0); TH1F *x1_7_eta1 = new TH1F("x1_7_eta1", "x1_7_eta1", 50, -6, 0); TH1F *x1_7_eta2 = new TH1F("x1_7_eta2", "x1_7_eta2", 50, -6, 0); TH1F *x2_3_eta1 = new TH1F("x2_3_eta1", "x2_3_eta1", 50, -6, 0); TH1F *x2_3_eta2 = new TH1F("x2_3_eta2", "x2_3_eta2", 50, -6, 0); TH1F *x2_5_eta1 = new TH1F("x2_5_eta1", "x2_5_eta1", 50, -6, 0); TH1F *x2_5_eta2 = new TH1F("x2_5_eta2", "x2_5_eta2", 50, -6, 0); TH1F *x2_7_eta1 = new TH1F("x2_7_eta1", "x2_7_eta1", 50, -6, 0); TH1F *x2_7_eta2 = new TH1F("x2_7_eta2", "x2_7_eta2", 50, -6, 0); TH1F *Q2_3_eta1 = new TH1F("Q2_3_eta1", "Q2_3_eta1", 50, 0, 500); TH1F *Q2_3_eta2 = new TH1F("Q2_3_eta2", "Q2_3_eta2", 50, 0, 500); TH1F *Q2_5_eta1 = new TH1F("Q2_5_eta1", "Q2_5_eta1", 50, 0, 500); TH1F *Q2_5_eta2 = new TH1F("Q2_5_eta2", "Q2_5_eta2", 50, 0, 500); TH1F *Q2_7_eta1 = new TH1F("Q2_7_eta1", "Q2_7_eta1", 50, 0, 500); TH1F *Q2_7_eta2 = new TH1F("Q2_7_eta2", "Q2_7_eta2", 50, 0, 500); */ /* TH1F *x1Fx_Q100_3 = new TH1F("x1Fx_Q100_3", "x1*f(x1) for Q = 100 +-#epsilon, cteq3d", 50, -6, 0); TH1F *x2Fx_Q100_3 = new TH1F("x2Fx_Q100_3", "x2*f(x2) for Q = 100 +-#epsilon, cteq3d", 50, -6, 0); TH1F *x1Fx_Q100_5 = new TH1F("x1Fx_Q100_5", "x1*f(x1) for Q = 100 +-#epsilon, grv94", 50, -6, 0); TH1F *x2Fx_Q100_5 = new TH1F("x2Fx_Q100_5", "x2*f(x2) for Q = 100 +-#epsilon, grv94", 50, -6, 0); TH1F *x1Fx_Q100_7 = new TH1F("x1Fx_Q100_7", "x1*f(x1) for Q = 100 +-#epsilon, cteq5l", 50, -6, 0); TH1F *x2Fx_Q100_7 = new TH1F("x2Fx_Q100_7", "x2*f(x2) for Q = 100 +-#epsilon, cteq5l", 50, -6, 0); */ //----------------------------------------------------------------------- Int_t njet2; Float_t pxj2[50], pyj2[50], pzj2[50], ptj2[50], etaj2[50]; Float_t ptJet[50]; Float_t x1, x2, xf, Q2; Float_t sigma[4], coef; Float_t LUMI = 1.; //Luminosity needed nb^-1 Int_t nevt; //number of events Int_t nfiles; //number of files //----------------------------------------------------------------------- sigma[1] = 67340000; sigma[2] = 80360000; sigma[3] = 89270000; nfiles = 100; nevt = nfiles*200000; char filename[100]; for (Int_t file_num = 1; file_num <= nfiles; file_num++) { sprintf (filename, "rfio:/castor/cern.ch/user/a/andrey/batch/pp14TeV_pdf5_%d.root", file_num); cout <<" file_num "<< file_num << endl; TFile *f1_1 = TFile::Open(filename); TTree *h107 = (TTree*)f1_1->Get("h101"); TTree *h105 = (TTree*)f1_1->Get("h102"); TTree *h103 = (TTree*)f1_1->Get("h103"); h103->SetBranchAddress("njet2",&njet2); h103->SetBranchAddress("etaj2",&etaj2); h103->SetBranchAddress("pxj2",&pxj2); h103->SetBranchAddress("pyj2",&pyj2); h103->SetBranchAddress("pzj2",&pzj2); h103->SetBranchAddress("ptj2",&ptj2); h103->SetBranchAddress("x1",&x1); h103->SetBranchAddress("x2",&x2); h103->SetBranchAddress("xf",&xf); h103->SetBranchAddress("Q2",&Q2); coef = sigma[1]*LUMI/nevt; cout << "coef1 for 'CTEQ3D' ->> "<GetEntries(); for (Int_t i=0; iGetEntry(i); // if(log10(x1)<-3. || log10(x2)<-3.) // { for (Int_t j=0; j Fill(ptj2[j], coef); hEtaJets3 -> Fill(etaj2[j], coef); if (2 < fabs(etaj2[j]) && fabs(etaj2[j]) < 4) {hPtJets3_eta1 -> Fill(ptj2[j], coef);} if (4 < fabs(etaj2[j]) && fabs(etaj2[j]) < 5) {hPtJets3_eta2 -> Fill(ptj2[j], coef);} } // if (fabs(Q2-100) < 10) {x1Fx_Q100_3 -> Fill(log10(x1), x1); x2Fx_Q100_3 -> Fill(log10(x2), x2);} x1_3 -> Fill(log10(x1)); x2_3 -> Fill(log10(x2)); xf_3 -> Fill(log10(xf)); Q2_3 -> Fill(Q2); // } } h105->SetBranchAddress("njet2",&njet2); h105->SetBranchAddress("etaj2",&etaj2); h105->SetBranchAddress("pxj2",&pxj2); h105->SetBranchAddress("pyj2",&pyj2); h105->SetBranchAddress("pzj2",&pzj2); h105->SetBranchAddress("ptj2",&ptj2); h105->SetBranchAddress("x1",&x1); h105->SetBranchAddress("x2",&x2); h105->SetBranchAddress("xf",&xf); h105->SetBranchAddress("Q2",&Q2); coef = sigma[2]*LUMI/nevt; cout << "coef1 for 'GRV94M' ->> "<GetEntries(); for (Int_t i=0; iGetEntry(i); // if(log10(x1)<-3. || log10(x2)<-3.) // { for (Int_t j=0; j Fill(ptj2[j], coef); hEtaJets5 -> Fill(etaj2[j], coef); if (2 < fabs(etaj2[j]) && fabs(etaj2[j]) < 4) {hPtJets5_eta1 -> Fill(ptj2[j], coef);} if (4 < fabs(etaj2[j]) && fabs(etaj2[j]) < 5) {hPtJets5_eta2 -> Fill(ptj2[j], coef);} } // if (fabs(Q2-100) < 10) {x1Fx_Q100_5 -> Fill(log10(x1),x1); x2Fx_Q100_5 -> Fill(log10(x2),x2);} x1_5 -> Fill(log10(x1)); x2_5 -> Fill(log10(x2)); xf_5 -> Fill(log10(xf)); Q2_5 -> Fill(Q2); // } } h107->SetBranchAddress("njet2",&njet2); h107->SetBranchAddress("etaj2",&etaj2); h107->SetBranchAddress("pxj2",&pxj2); h107->SetBranchAddress("pyj2",&pyj2); h107->SetBranchAddress("pzj2",&pzj2); h107->SetBranchAddress("ptj2",&ptj2); h107->SetBranchAddress("x1",&x1); h107->SetBranchAddress("x2",&x2); h107->SetBranchAddress("xf",&xf); h107->SetBranchAddress("Q2",&Q2); coef = sigma[3]*LUMI/nevt; cout << "coef1 for 'CTEQ5L' ->> "<GetEntries(); for (Int_t i=0; iGetEntry(i); // if(log10(x1)<-3. || log10(x2)<-3.) // { for (Int_t j=0; j Fill(ptj2[j], coef); hEtaJets7 -> Fill(etaj2[j], coef); if (2 < fabs(etaj2[j]) && fabs(etaj2[j]) < 4) {hPtJets7_eta1 -> Fill(ptj2[j], coef);} if (4 < fabs(etaj2[j]) && fabs(etaj2[j]) < 5) {hPtJets7_eta2 -> Fill(ptj2[j], coef);} } // if (fabs(Q2-100) < 10) {x1Fx_Q100_7 -> Fill(log10(x1),x1); x2Fx_Q100_7 -> Fill(log10(x2),x2);} x1_7 -> Fill(log10(x1)); x2_7 -> Fill(log10(x2)); xf_7 -> Fill(log10(xf)); Q2_7 -> Fill(Q2); // } } } TFile ftemp("temp6.root", "recreate"); hPtJets3 -> Sumw2(); hPtJets3_eta1 -> Sumw2(); hPtJets3_eta2 -> Sumw2(); hPtJets3 -> Write(); hPtJets3_eta1 -> Write(); hPtJets3_eta2 -> Write(); hPtJets5 -> Sumw2(); hPtJets5_eta1 -> Sumw2(); hPtJets5_eta2 -> Sumw2(); hPtJets5 -> Write(); hPtJets5_eta1 -> Write(); hPtJets5_eta2 -> Write(); hPtJets7 -> Sumw2(); hPtJets7_eta1 -> Sumw2(); hPtJets7_eta2 -> Sumw2(); hPtJets7 -> Write(); hPtJets7_eta1 -> Write(); hPtJets7_eta2 -> Write(); /* x1Fx_Q100_3 -> Write(); x2Fx_Q100_3 -> Write(); x1Fx_Q100_5 -> Write(); x2Fx_Q100_5 -> Write(); x1Fx_Q100_7 -> Write(); x2Fx_Q100_7 -> Write(); */ x1_3 -> Write(); x1_5 -> Write(); x1_7 -> Write(); x2_3 -> Write(); x2_5 -> Write(); x2_7 -> Write(); xf_3 -> Write(); xf_5 -> Write(); xf_7 -> Write(); Q2_3 -> Write(); Q2_5 -> Write(); Q2_7 -> Write(); /* x1_3_eta1 -> Write(); x1_3_eta2 -> Write(); x1_5_eta1 -> Write(); x1_5_eta2 -> Write(); x1_7_eta1 -> Write(); x1_7_eta2 -> Write(); */ //ftemp.Close(); }