double f55( double mass ) /* differential dsigma-1/dM========>dsigma-1/(M*dlnM ) */
{ /* dsigma-1/dm =-1/sigma^2*dsigma/dm */
double h = 0.000001;
//return (( f44( exp(log(mass) + h) ) - f44( exp(log(mass) ) )) / h) / mass;
;return (( f33( exp(log(mass) + h) ) - f33( exp(log(mass) ) )) / h) * (-1.0 / pow( f33( mass ),2 ));
return ((( f33( exp(log(mass) + h) ) - f33( exp(log(mass) ) )) / h) / mass) * (-1.0 / pow( f33( mass ),2 ));
}
int main(void) {
f0();
f1();
f2();
f3();
f4();
f5();
f6();
f7();
f8();
f9();
f10();
f11();
f12();
f13();
f14();
f15();
f16();
f17();
f18();
f19();
f20();
f21();
f22();
f23();
f24();
f25();
f26();
f27();
f28();
f29();
f30();
f31();
f32();
f33();
f34();
f35();
f36();
f37();
f38();
f39();
f40();
f41();
f42();
f43();
f44();
f45();
f46();
f47();
f48();
f49();
f50();
f51();
f52();
return 0;
}
int main(void) {
f0();
f1();
f2();
f3();
f4();
f5();
f6();
f7();
f8();
f9();
f10();
f11();
f12();
f13();
f14();
f15();
f16();
f17();
f18();
f19();
f20();
f21();
f22();
f23();
f24();
f25();
f26();
f27();
f28();
f29();
f30();
f31();
f32();
f33();
f34();
f35();
f36();
f37();
return 0;
}
// primitive types, by value in, out
void
test1()
{
printf("\ntesting primitive type functions: fxx(0) ...\n");
f11(0);
f12(0);
f13(0);
f14(0);
f15(0);
f16(0);
f17(0);
f18(0);
f19(0);
f20(0);
f21(0);
f22(0);
f23(0);
f24(0);
f25(0);
f31(0);
f32(0);
f33(0);
f34(0);
f35(0);
f36(0);
f37(0);
f38(0);
f39(0);
f40(0);
f41(0);
f42(0);
f43(0);
f44(0);
f45(0);
}
int main(void) {
f0();
f1();
f2();
f3();
f4();
f5();
f6();
f7();
f8();
f9();
f10();
f11();
f12();
f13();
f14();
f15();
f16();
f17();
f18();
f19();
f20();
f21();
f22();
f23();
f24();
f25();
f26();
f27();
f28();
f29();
f30();
f31();
f32();
f33();
f34();
f35();
f36();
f37();
f38();
f39();
f40();
f41();
f42();
f43();
f44();
f45();
f46();
f47();
f48();
f49();
f50();
f51();
f52();
f53();
f54();
f55();
f56();
f57();
f58();
f59();
f60();
f61();
f62();
f63();
f64();
f65();
f66();
f67();
f68();
f69();
f70();
f71();
f72();
f73();
f74();
f75();
f76();
f77();
f78();
f79();
f80();
f81();
f82();
f83();
f84();
f85();
f86();
f87();
f88();
f89();
f90();
f91();
f92();
f93();
f94();
f95();
f96();
f97();
f98();
f99();
f100();
f101();
f102();
f103();
f104();
f105();
f106();
f107();
f108();
f109();
f110();
f111();
f112();
f113();
f114();
f115();
f116();
f117();
f118();
f119();
f120();
f121();
f122();
f123();
f124();
f125();
f126();
f127();
f128();
f129();
f130();
f131();
f132();
f133();
f134();
f135();
f136();
f137();
f138();
f139();
f140();
f141();
f142();
f143();
f144();
f145();
f146();
f147();
f148();
f149();
f150();
f151();
f152();
f153();
f154();
f155();
f156();
f157();
f158();
f159();
f160();
f161();
f162();
f163();
f164();
f165();
f166();
f167();
f168();
f169();
f170();
f171();
f172();
f173();
f174();
f175();
f176();
f177();
f178();
f179();
f180();
f181();
f182();
f183();
f184();
f185();
f186();
f187();
f188();
f189();
f190();
f191();
f192();
f193();
f194();
f195();
f196();
f197();
f198();
f199();
return 0;
}
int main(void) {
f0();
f1();
f2();
f3();
f4();
f5();
f6();
f7();
f8();
f9();
f10();
f11();
f12();
f13();
f14();
f15();
f16();
f17();
f18();
f19();
f20();
f21();
f22();
f23();
f24();
f25();
f26();
f27();
f28();
f29();
f30();
f31();
f32();
f33();
f34();
f35();
f36();
f37();
f38();
f39();
f40();
f41();
f42();
f43();
f44();
f45();
f46();
f47();
f48();
f49();
f50();
f51();
f52();
f53();
f54();
f55();
f56();
f57();
f58();
f59();
f60();
f61();
f62();
f63();
f64();
f65();
f66();
f67();
f68();
f69();
f70();
f71();
f72();
f73();
f74();
f75();
f76();
f77();
f78();
f79();
f80();
f81();
f82();
f83();
f84();
f85();
f86();
f87();
f88();
f89();
f90();
f91();
f92();
f93();
f94();
f95();
f96();
f97();
f98();
f99();
f100();
f101();
f102();
f103();
f104();
f105();
f106();
f107();
f108();
f109();
f110();
f111();
f112();
f113();
f114();
f115();
f116();
f117();
f118();
f119();
f120();
f121();
f122();
f123();
return 0;
}
int f32(){
totalSum += 32;
return f33() + 32;
}
void CompareWithNikos_23Jun08(){
/*
{
gROOT->ProcessLine(".L ~/tdrstyle.C");
setTDRStyle();
TFile f2("histos.root");
TH1F *eff = (TH1F*)f2.Get("gsfelectron_gsfelectronIso_Eta_reReco");
//
TFile f("./Nikos-reReco/sc_eta_iso.root");
// common for all sequences ...........................
TH1F *base = (TH1F*)f.Get("base");
TH1F *destination = (TH1F*)f.Get("destination");
TGraphAsymmErrors *g1 = new TGraphAsymmErrors();
g1->BayesDivide(destination, base, "");
g1->GetYaxis()->SetRangeUser(0.5, 1.05);
eff->SetLineColor(2);
eff->SetMarkerStyle(22);
eff->SetMarkerSize(1.4);
eff->SetMarkerColor(2);
//
g1->GetYaxis()->SetTitle("Efficiency");
g1->GetXaxis()->SetTitle("#eta");
g1->Draw("APE");
eff->Draw("esame");
gPad->Update();
*/
TString stringname1 = "supercluster_gsfelectron";
TString stringname2 = "gsfelectron_gsfelectronIso";
TString stringname3 = "gsfelectronIso_gsfelectronIsoLooseID";
TString stringname4 = "gsfelectronIsoLooseID_HLT";
// TFile myhistos("histos_sbs.root");
// TFile myhistos("histos_truth.root");
TFile myhistos("histos.root");
TString name11 = stringname1 + TString("_Eta_reReco");
TString name12 = stringname1 + TString("_Eta_Reco");
TString name13 = stringname1 + TString("_Phi_reReco");
TString name14 = stringname1 + TString("_Phi_Reco");
TString name15 = stringname1 + TString("_Pt_reReco");
TString name16 = stringname1 + TString("_Pt_Reco");
TString name17 = stringname1 + TString("_Pt_Eta_reReco");
TString name18 = stringname1 + TString("_Pt_Eta_Reco");
TString name19 = stringname1 + TString("_Phi_Eta_reReco");
TString name10 = stringname1 + TString("_Phi_Eta_Reco");
TH1F* h11 = (TH1F*) myhistos.Get(name11);
TH1F* h12 = (TH1F*) myhistos.Get(name12);
TH1F* h13 = (TH1F*) myhistos.Get(name13);
TH1F* h14 = (TH1F*) myhistos.Get(name14);
TH1F* h15 = (TH1F*) myhistos.Get(name15);
TH1F* h16 = (TH1F*) myhistos.Get(name16);
TH2F* h17 = (TH2F*) myhistos.Get(name17);
TH2F* h18 = (TH2F*) myhistos.Get(name18);
TH2F* h19 = (TH2F*) myhistos.Get(name19);
TH2F* h10 = (TH2F*) myhistos.Get(name10);
TFile f11("./Nikos-reReco/sc_eta_recoEle.root");
TFile f12("./Nikos-Reco/sc_eta_recoEle.root");
TFile f13("./Nikos-reReco/sc_phi_recoEle.root");
TFile f14("./Nikos-Reco/sc_phi_recoEle.root");
TFile f15("./Nikos-reReco/sc_et_recoEle.root");
TFile f16("./Nikos-Reco/sc_et_recoEle.root");
TH1F* base11 = (TH1F*)f11.Get("base");
TH1F* destination11 = (TH1F*)f11.Get("destination");
TH1F* base12 = (TH1F*)f12.Get("base");
TH1F* destination12 = (TH1F*)f12.Get("destination");
TH1F* base13 = (TH1F*)f13.Get("base");
TH1F* destination13 = (TH1F*)f13.Get("destination");
TH1F* base14 = (TH1F*)f14.Get("base");
TH1F* destination14 = (TH1F*)f14.Get("destination");
TH1F* base15 = (TH1F*)f15.Get("base");
TH1F* destination15 = (TH1F*)f15.Get("destination");
TH1F* base16 = (TH1F*)f16.Get("base");
TH1F* destination16 = (TH1F*)f16.Get("destination");
// Process(*h11, *base11, *destination11);
// Process(*h12, *base12, *destination12);
// Process(*h13, *base13, *destination13);
// Process(*h14, *base14, *destination14);
// Process(*h15, *base15, *destination15);
// Process(*h16, *base16, *destination16);
TString name21 = stringname2 + TString("_Eta_reReco");
TString name22 = stringname2 + TString("_Eta_Reco");
TString name23 = stringname2 + TString("_Phi_reReco");
TString name24 = stringname2 + TString("_Phi_Reco");
TString name25 = stringname2 + TString("_Pt_reReco");
TString name26 = stringname2 + TString("_Pt_Reco");
TString name27 = stringname2 + TString("_Pt_Eta_reReco");
TString name28 = stringname2 + TString("_Pt_Eta_Reco");
TString name29 = stringname2 + TString("_Phi_Eta_reReco");
TString name20 = stringname2 + TString("_Phi_Eta_Reco");
TH1F* h21 = (TH1F*) myhistos.Get(name21);
TH1F* h22 = (TH1F*) myhistos.Get(name22);
TH1F* h23 = (TH1F*) myhistos.Get(name23);
TH1F* h24 = (TH1F*) myhistos.Get(name24);
TH1F* h25 = (TH1F*) myhistos.Get(name25);
TH1F* h26 = (TH1F*) myhistos.Get(name26);
TH2F* h27 = (TH2F*) myhistos.Get(name27);
TH2F* h28 = (TH2F*) myhistos.Get(name28);
TH2F* h29 = (TH2F*) myhistos.Get(name29);
TH2F* h20 = (TH2F*) myhistos.Get(name20);
TFile f21("./Nikos-reReco/sc_eta_iso.root");
TFile f22("./Nikos-Reco/sc_eta_iso.root");
TFile f23("./Nikos-reReco/sc_phi_iso.root");
TFile f24("./Nikos-Reco/sc_phi_iso.root");
TFile f25("./Nikos-reReco/sc_et_iso.root");
TFile f26("./Nikos-Reco/sc_et_iso.root");
TH1F* base21 = (TH1F*)f21.Get("base");
TH1F* destination21 = (TH1F*)f21.Get("destination");
TH1F* base22 = (TH1F*)f22.Get("base");
TH1F* destination22 = (TH1F*)f22.Get("destination");
TH1F* base23 = (TH1F*)f23.Get("base");
TH1F* destination23 = (TH1F*)f23.Get("destination");
TH1F* base24 = (TH1F*)f24.Get("base");
TH1F* destination24 = (TH1F*)f24.Get("destination");
TH1F* base25 = (TH1F*)f25.Get("base");
TH1F* destination25 = (TH1F*)f25.Get("destination");
TH1F* base26 = (TH1F*)f26.Get("base");
TH1F* destination26 = (TH1F*)f26.Get("destination");
// Process(*h21, *base21, *destination21);
// Process(*h22, *base22, *destination22);
// Process(*h23, *base23, *destination23);
// Process(*h24, *base24, *destination24);
// Process(*h25, *base25, *destination25);
// Process(*h26, *base26, *destination26);
TString name31 = stringname3 + TString("_Eta_reReco");
TString name32 = stringname3 + TString("_Eta_Reco");
TString name33 = stringname3 + TString("_Phi_reReco");
TString name34 = stringname3 + TString("_Phi_Reco");
TString name35 = stringname3 + TString("_Pt_reReco");
TString name36 = stringname3 + TString("_Pt_Reco");
TString name37 = stringname3 + TString("_Pt_Eta_reReco");
TString name38 = stringname3 + TString("_Pt_Eta_Reco");
TString name39 = stringname3 + TString("_Phi_Eta_reReco");
TString name30 = stringname3 + TString("_Phi_Eta_Reco");
TH1F* h31 = (TH1F*) myhistos.Get(name31);
TH1F* h32 = (TH1F*) myhistos.Get(name32);
TH1F* h33 = (TH1F*) myhistos.Get(name33);
TH1F* h34 = (TH1F*) myhistos.Get(name34);
TH1F* h35 = (TH1F*) myhistos.Get(name35);
TH1F* h36 = (TH1F*) myhistos.Get(name36);
TH2F* h37 = (TH2F*) myhistos.Get(name37);
TH2F* h38 = (TH2F*) myhistos.Get(name38);
TH2F* h39 = (TH2F*) myhistos.Get(name39);
TH2F* h30 = (TH2F*) myhistos.Get(name30);
TFile f31("./Nikos-reReco/sc_eta_id_robust.root");
TFile f32("./Nikos-Reco/sc_eta_id_robust.root");
TFile f33("./Nikos-reReco/sc_phi_id_robust.root");
TFile f34("./Nikos-Reco/sc_phi_id_robust.root");
TFile f35("./Nikos-reReco/sc_et_id_robust.root");
TFile f36("./Nikos-Reco/sc_et_id_robust.root");
TH1F* base31 = (TH1F*)f31.Get("base");
TH1F* destination31 = (TH1F*)f31.Get("destination");
TH1F* base32 = (TH1F*)f32.Get("base");
TH1F* destination32 = (TH1F*)f32.Get("destination");
TH1F* base33 = (TH1F*)f33.Get("base");
TH1F* destination33 = (TH1F*)f33.Get("destination");
TH1F* base34 = (TH1F*)f34.Get("base");
TH1F* destination34 = (TH1F*)f34.Get("destination");
TH1F* base35 = (TH1F*)f35.Get("base");
TH1F* destination35 = (TH1F*)f35.Get("destination");
TH1F* base36 = (TH1F*)f36.Get("base");
TH1F* destination36 = (TH1F*)f36.Get("destination");
TString name41 = stringname4 + TString("_Eta_reReco");
TString name42 = stringname4 + TString("_Eta_Reco");
TString name43 = stringname4 + TString("_Phi_reReco");
TString name44 = stringname4 + TString("_Phi_Reco");
TString name45 = stringname4 + TString("_Pt_reReco");
TString name46 = stringname4 + TString("_Pt_Reco");
TString name47 = stringname4 + TString("_Pt_Eta_reReco");
TString name48 = stringname4 + TString("_Pt_Eta_Reco");
TString name49 = stringname4 + TString("_Phi_Eta_reReco");
TString name40 = stringname4 + TString("_Phi_Eta_Reco");
TH1F* h41 = (TH1F*) myhistos.Get(name41);
TH1F* h42 = (TH1F*) myhistos.Get(name42);
TH1F* h43 = (TH1F*) myhistos.Get(name43);
TH1F* h44 = (TH1F*) myhistos.Get(name44);
TH1F* h45 = (TH1F*) myhistos.Get(name45);
TH1F* h46 = (TH1F*) myhistos.Get(name46);
TH2F* h47 = (TH2F*) myhistos.Get(name47);
TH2F* h48 = (TH2F*) myhistos.Get(name48);
TH2F* h49 = (TH2F*) myhistos.Get(name49);
TH2F* h40 = (TH2F*) myhistos.Get(name40);
TFile f41("./Nikos-reReco/sc_eta_trigger_idRobust.root");
TFile f42("./Nikos-Reco/sc_eta_trigger_idRobust.root");
TFile f43("./Nikos-reReco/sc_phi_trigger_idRobust.root");
TFile f44("./Nikos-Reco/sc_phi_trigger_idRobust.root");
TFile f45("./Nikos-reReco/sc_et_trigger_idRobust.root");
TFile f46("./Nikos-Reco/sc_et_trigger_idRobust.root");
TH1F* base41 = (TH1F*)f41.Get("base");
TH1F* destination41 = (TH1F*)f41.Get("destination");
TH1F* base42 = (TH1F*)f42.Get("base");
TH1F* destination42 = (TH1F*)f42.Get("destination");
TH1F* base43 = (TH1F*)f43.Get("base");
TH1F* destination43 = (TH1F*)f43.Get("destination");
TH1F* base44 = (TH1F*)f44.Get("base");
TH1F* destination44 = (TH1F*)f44.Get("destination");
TH1F* base45 = (TH1F*)f45.Get("base");
TH1F* destination45 = (TH1F*)f45.Get("destination");
TH1F* base46 = (TH1F*)f46.Get("base");
TH1F* destination46 = (TH1F*)f46.Get("destination");
// Make 1D comparison plots
makeplots1D( *h11, *base11, *destination11, name11);
makeplots1D( *h12, *base12, *destination12, name12);
makeplots1D( *h13, *base13, *destination13, name13);
makeplots1D( *h14, *base14, *destination14, name14);
makeplots1D( *h15, *base15, *destination15, name15);
makeplots1D( *h16, *base16, *destination16, name16);
makeplots1D( *h21, *base21, *destination21, name21);
makeplots1D( *h22, *base22, *destination22, name22);
makeplots1D( *h23, *base23, *destination23, name23);
makeplots1D( *h24, *base24, *destination24, name24);
makeplots1D( *h25, *base25, *destination25, name25);
makeplots1D( *h26, *base26, *destination26, name26);
makeplots1D( *h31, *base31, *destination31, name31);
makeplots1D( *h32, *base32, *destination32, name32);
makeplots1D( *h33, *base33, *destination33, name33);
makeplots1D( *h34, *base34, *destination34, name34);
makeplots1D( *h35, *base35, *destination35, name35);
makeplots1D( *h36, *base36, *destination36, name36);
makeplots1D( *h41, *base41, *destination41, name41);
makeplots1D( *h42, *base42, *destination42, name42);
makeplots1D( *h43, *base43, *destination43, name43);
makeplots1D( *h44, *base44, *destination44, name44);
makeplots1D( *h45, *base45, *destination45, name45);
makeplots1D( *h46, *base46, *destination46, name46);
// Make 2D efficiency plots
makeplots2D( *h17, name17);
makeplots2D( *h18, name18);
makeplots2D( *h19, name19);
makeplots2D( *h10, name10);
makeplots2D( *h27, name27);
makeplots2D( *h28, name28);
makeplots2D( *h29, name29);
makeplots2D( *h20, name20);
makeplots2D( *h37, name37);
makeplots2D( *h38, name38);
makeplots2D( *h39, name39);
makeplots2D( *h30, name30);
makeplots2D( *h47, name47);
makeplots2D( *h48, name48);
makeplots2D( *h49, name49);
makeplots2D( *h40, name40);
}
double f44(double mass) /*1/sigmaM*/
{
return 1.0 / f33( mass );
}