void drawEventRPhi(double vxlb, double vylb, double vxl0, double vyl0, double pmu1, double phimu1, double pmu2, double phimu2, double ppr, double phipr, double ppi, double phipi, int colors = 0) { Color_t colMu1(1), colMu2(1), colPr(1), colPi(1), colL0(1), colLb(1); switch(colors) { case 0: colLb=1; colL0=2; colMu1=2; colMu2=2; colPr=3; colPi=4; break; case 1: colLb=11; colL0=50; colMu1=50; colMu2=50; colPr=8; colPi=9; break; } // draw the vertices TMarker *m; const double xlb = vxlb; const double ylb = vylb; m = new TMarker(xlb,ylb,7); m->SetMarkerColor(colLb); m->Draw(); const double xl0 = vxl0; const double yl0 = vyl0; m = new TMarker(xl0,yl0,7); m->SetMarkerColor(colL0); m->Draw(); // draw the l0 flight line TLine *l; l = new TLine(vxlb,vylb,vxl0,vyl0); l->SetLineColor(colL0); l->SetLineStyle(2); l->Draw(); // draw the muons TArrow *a; const double xmu1 = xlb + scalemu * pmu1 * TMath::Cos(phimu1); const double ymu1 = ylb + scalemu * pmu1 * TMath::Sin(phimu1); a = new TArrow(xlb,ylb,xmu1,ymu1,.01,">"); a->SetLineColor(colMu1); a->Draw(); const double xmu2 = xlb + scalemu * pmu2 * TMath::Cos(phimu2); const double ymu2 = ylb + scalemu * pmu2 * TMath::Sin(phimu2); a = new TArrow(xlb,ylb,xmu2,ymu2,.01,">"); a->SetLineColor(colMu2); a->Draw(); // draw the p and pi const double xpr = xl0 + scalepr * ppr * TMath::Cos(phipr); const double ypr = yl0 + scalepr * ppr * TMath::Sin(phipr); a = new TArrow(xl0,yl0,xpr,ypr,.01,">"); a->SetLineColor(colPr); a->Draw(); const double xpi = xl0 + scalepi * ppi * TMath::Cos(phipi); const double ypi = yl0 + scalepi * ppi * TMath::Sin(phipi); a = new TArrow(xl0,yl0,xpi,ypi,.01,">"); a->SetLineColor(colPi); a->Draw(); }
void DrawLabel(TString txt, double x, double y, int mstyle, int col, double msize) { TMarker *m = new TMarker(x,y,mstyle); m->SetNDC(); m->SetMarkerSize(msize); m->SetMarkerColor(col); TLine *l = new TLine(); l->SetLineWidth(2); l->SetLineColor(col); m->Draw(); tex->SetTextSize(0.04); tex->SetTextAlign(12); tex->SetTextColor(col); tex->DrawLatex(x+0.025,y,txt); tex->SetTextColor(kBlack); }
void myMarkerText(Double_t x,Double_t y,Int_t color,Int_t mstyle, const char *text,Float_t msize) { Double_t tsize=0.06; TMarker *marker = new TMarker(x-(0.4*tsize),y,8); marker->SetMarkerColor(color); marker->SetNDC(); marker->SetMarkerStyle(mstyle); marker->SetMarkerSize(msize); marker->Draw(); TLatex l; l.SetTextAlign(12); //l.SetTextSize(tsize); l.SetNDC(); l.DrawLatex(x,y,text); }
void myMarkerText(Double_t x,Double_t y,Int_t color,Int_t mstyle,char *text) { // printf("**myMarker: text= %s\ m ",text); Double_t tsize=0.06; TMarker *marker = new TMarker(x-(0.4*tsize),y,8); marker->SetMarkerColor(color); marker->SetNDC(); marker->SetMarkerStyle(mstyle); marker->SetMarkerSize(2.0); marker->Draw(); TLatex l; l.SetTextAlign(12); //l.SetTextSize(tsize); l.SetNDC(); l.DrawLatex(x,y,text); }
void showHistogram1d(TH1* histogram, const std::string& xAxisTitle, Float_t* genX, const std::string& outputFileName) { TCanvas* canvas = new TCanvas("canvas", "canvas", 800, 600); canvas->SetFillColor(10); canvas->SetBorderSize(2); canvas->SetTopMargin(0.10); canvas->SetLeftMargin(0.16); canvas->SetRightMargin(0.14); canvas->SetBottomMargin(0.12); TAxis* xAxis = histogram->GetXaxis(); xAxis->SetTitle(xAxisTitle.data()); xAxis->SetTitleOffset(1.15); TAxis* yAxis = histogram->GetYaxis(); yAxis->SetTitle("Sampling Points"); yAxis->SetTitleOffset(1.60); histogram->SetLineColor(1); histogram->SetLineWidth(2); histogram->SetMarkerColor(1); histogram->SetMarkerStyle(20); histogram->Draw("e1p"); TMarker* genMarker = 0; if ( genX ) { genMarker = new TMarker(*genX, 0.10, 34); genMarker->SetMarkerColor(1); genMarker->SetMarkerSize(2); genMarker->Draw(); } canvas->Update(); size_t idx = outputFileName.find_last_of('.'); std::string outputFileName_plot = std::string(outputFileName, 0, idx); if ( idx != std::string::npos ) canvas->Print(std::string(outputFileName_plot).append(std::string(outputFileName, idx)).data()); canvas->Print(std::string(outputFileName_plot).append(".png").data()); canvas->Print(std::string(outputFileName_plot).append(".pdf").data()); canvas->Print(std::string(outputFileName_plot).append(".root").data()); delete genMarker; delete canvas; }
// Draw on a given canvas the full set of working points void overlayWorkingPoints(TCanvas *c1, TTree *signalTree, TTree *backgroundTree, const TString *cutFileNames, int markerColor, int markerStyle, TLegend *leg, const TString legendText){ // Now loop over working points for(int iwp = 0; iwp<nWP; iwp++){ // Load the working point from a ROOT file TFile *cutFile = new TFile(cutFileNames[iwp]); if( !cutFile ) assert(0); VarCut *cutObject = (VarCut*)cutFile->Get("cuts"); if( !cutObject ) assert(0); // Compute the efficiencies float effSignal, effBackground; findEfficiencies(signalTree, backgroundTree, effSignal, effBackground, cutObject); printf("Computed eff for cut from %s, effS= %.3f effB= %.3f\n", cutFileNames[iwp].Data(), effSignal, effBackground); // Make a marker and draw it. TMarker *marker = new TMarker(effSignal, 1.0-effBackground, 20); marker->SetMarkerSize(2); marker->SetMarkerColor(markerColor); marker->SetMarkerStyle(markerStyle); marker->Draw("same"); // Add marker to the legend only once. Do not draw the legend here, // it is drawn in the main function later if( iwp == 0 ){ if( !leg ) assert(0); leg->AddEntry(marker, legendText, "p"); } c1->Update(); cutFile->Close(); } }
void drawTransition(){ TBox *bTrans = new TBox(0., 0., 0.97, 0.97); if(iColor)bTrans->SetFillColor(kYellow-9); if (! iColor)bTrans->SetFillStyle(4017); bTrans->Draw(); TEllipse *eTrans1 = new TEllipse(0., 0., 0.9, 0.7, 70., 90., 0.); eTrans1->SetNoEdges(); eTrans1->SetFillColor(17); eTrans1->SetLineWidth(2); eTrans1->SetLineStyle(7); eTrans1->Draw(); TEllipse *eTrans2 = new TEllipse(0., 0., 0.9, 0.7, 0., 70., 0.); eTrans2->SetNoEdges(); eTrans2->SetFillColor(17); eTrans2->SetLineWidth(2); eTrans2->Draw(); TMarker *mCrit = new TMarker(0.29, 0.6625, 20); mCrit->SetMarkerSize(1.4); mCrit->Draw(); }
void long_Ay_nu_05() { gROOT->SetStyle("HALLA"); TCanvas *cn = new TCanvas("cn","cn",540,360); cn->Draw(); cn->UseCurrentStyle(); TH1F *frm = new TH1F("frm","",100,0.,10.); frm->GetXaxis()->SetTitle("#nu (GeV)"); frm->GetYaxis()->SetTitle("Ay for Q2=0.456 (GeV/c)2"); frm->SetMinimum(0); // frm->SetMinimum(0); // frm->SetMaximum(1.0); frm->SetMaximum(0.35); frm->UseCurrentStyle(); frm->Draw(); frm->SetAxisRange(0.120,0.350,"X"); // frm->SetAxisRange(0.5,1.1,"X"); // TF1* galster = new TF1("galster","x/(4.*0.938*.938)*1.91/(1.+x/.71)^2/(1.+5.6*x/(4.*.938*.938))",0.,4.); // galster->SetLineColor(6); // galster->SetLineStyle(3); // galster->SetLineWidth(2); TF1 *genf = new TF1("genf",genff,1.,10.,1); genf->SetLineColor(2); genf->SetLineStyle(2); genf->SetParameter(0,1.); // match to Madey point just below 1.5 // genf->SetParameter(0,.0411/genf->Eval(1.45)); genf->SetParameter(0,-0.558645); // TF1 *bbba05 = new TF1("BBBA05",gen_bbba05,0.,10.,0); // bbba05->SetLineColor(7); // bbba05->SetLineStyle(3); TMultiGraph* mgrDta = new TMultiGraph("Data","G_{E}^{n}"); TLegend *legDta = new TLegend(.54,.6,.875,.90,"","brNDC"); TMultiGraph* wgr = mgrDta; TLegend *wlg = legDta; // the data legDta->SetBorderSize(0); // turn off border legDta->SetFillStyle(0); datafile_t *f = datafiles; TGraph* gr=0; while ( f && f->filename ) { gr=OneGraph(f); if (gr) { if (f->lnpt) { mgrDta->Add(gr,f->lnpt); legDta->AddEntry(gr,f->label,f->lnpt); } else if (gr->GetMarkerStyle()>=20) { mgrDta->Add(gr,"p"); legDta->AddEntry(gr,f->label,"p"); } else { mgrDta->Add(gr,"l"); legDta->AddEntry(gr,f->label,"l"); } } f++; } mgrDta->Draw("p"); // legDta->Draw(); TF1 *theFit = new TF1("theFit","pol0"); gr->Fit(theFit); theFit->Draw("same"); TMultiGraph* mgrThry = new TMultiGraph("Theory","G_{E}^{n}"); TLegend *legThry = new TLegend(.54,.3,.875,.6,"","brNDC"); wgr = mgrThry; wlg = legThry; // the theory wlg->SetBorderSize(0); // turn off border wlg->SetFillStyle(0); f = theoryfiles1; gr=0; while ( f && f->filename ) { gr=OneGraph(f); if (gr) { TGraphAsymmErrors *egr = dynamic_cast<TGraphAsymmErrors*>(gr); if (egr && egr->GetN()>1 && egr->GetEYhigh() && egr->GetEYhigh()[1]>0) { gr = toerror_band(egr); gr->SetFillStyle(3000+f->style); } if (f->lnpt) { wgr->Add(gr,f->lnpt); wlg->AddEntry(gr,f->label,f->lnpt); } else if (gr->GetMarkerStyle()>=20) { wgr->Add(gr,"p"); wlg->AddEntry(gr,f->label,"p"); } else { wgr->Add(gr,"l"); wlg->AddEntry(gr,f->label,"l"); } } f++; } // genf->Draw("same"); mgrThry->Draw("c"); // galster->Draw("same"); // bbba05->Draw("same"); // legThry->AddEntry(genf,"F_{2}/F_{1} #propto ln^{2}(Q^{2}/#Lambda^{2})/Q^{2}","l"); // legThry->AddEntry(galster,"Galster fit","l"); // legThry->AddEntry(bbba05,"BBBA05","l"); // legThry->Draw(); // legDta->Draw(); // draw a line at 1 cn->Modified(); cn->Update(); cn->SaveAs(Form("%s.eps",psfile)); cn->SaveAs(Form("%s.root",psfile)); gSystem->Exec(Form("./replace_symbols.pl %s.eps",psfile)); return; // LEAVING HERE // now an overlay, hopefully matching dimensions // remove everything but the graph cn->Update(); TList *clist = cn->GetListOfPrimitives(); TFrame* frame = cn->GetFrame(); for (int i=0; i<clist->GetSize(); ) { if (clist->At(i) != frame) { clist->RemoveAt(i); } else i++; } // draw markers in the corners TMarker *mkr = new TMarker(frame->GetX1(),frame->GetY1(),2); mkr->Draw(); mkr = new TMarker(frame->GetX2(),frame->GetY1(),2); mkr->Draw(); mkr = new TMarker(frame->GetX1(),frame->GetY2(),2); mkr->Draw(); mkr = new TMarker(frame->GetX2(),frame->GetY2(),2); mkr->Draw(); frame->SetLineColor(10); cn->Update(); datafile_t miller = { "figure_input/Miller/lattice.GEn.rtf","Miller", "[0]","[1]","[1]-[3]","[2]-[1]",0,0,1,3,"F" }; gr = OneGraph(&miller); TGraphAsymmErrors* egr = dynamic_cast<TGraphAsymmErrors*>(gr); if (egr && egr->GetEYhigh() && egr->GetEYhigh()[egr->GetN()/2]>0) { gr = toerror_band(egr); gr->SetLineStyle(1); gr->SetFillColor(gr->GetLineColor()); gr->SetFillStyle(3000+miller.style); } gr->Draw("F"); cn->Update(); cn->SaveAs("gen_Miller_Overlay.eps"); cn->SaveAs("gen_Miller_Overlay.root"); }
void q4gmn_jan2009() { gROOT->SetStyle("HALLA"); TCanvas *cn = new TCanvas("cn"); cn->Draw(); // TH1F *frm = new TH1F("frm","",100,0.,5.5); TH1F *frm = new TH1F("frm","",100,0.,20.); frm->GetXaxis()->SetTitle("Q^{2} [GeV^{2}]"); frm->GetXaxis()->CenterTitle(); frm->GetYaxis()->SetTitle("Q^{4} G_{M}^{n}/#mu_{n}"); frm->GetYaxis()->CenterTitle(); frm->SetMinimum(0.0); frm->SetMaximum(0.6); frm->UseCurrentStyle(); frm->Draw(); TMultiGraph* mgrDta = new TMultiGraph("Data","G_{M}^{n}/#mu_{n}G_{D}"); TLegend *legDta = new TLegend(0.5029,.2076,.9095,.3877,"","brNDC"); TMultiGraph* wgr = mgrDta; TLegend *wlg = legDta; // the data legDta->SetBorderSize(0); // turn off border legDta->SetFillStyle(0); datafile_t *f = datafiles; TGraph* gr=0; TGraph* ogr=0; while ( f && f->filename ) { ogr=OneGraph(f); if (ogr) { gr = fromGMnGDtoQ4GMn(ogr); if (f->lnpt) { mgrDta->Add(gr,f->lnpt); legDta->AddEntry(gr,f->label,f->lnpt); } else if (gr->GetMarkerStyle()>=20) { mgrDta->Add(gr,"p"); legDta->AddEntry(gr,f->label,"p"); } else { mgrDta->Add(gr,"l"); legDta->AddEntry(gr,f->label,"l"); } } f++; } mgrDta->Draw("p"); legDta->Draw(); // don't draw the data legend TMultiGraph* mgrThry = new TMultiGraph("Theory",""); TLegend *legThry = new TLegend(0.1868,.1949,.5287,.3432,"","brNDC"); wgr = mgrThry; wlg = legThry; // the theory wlg->SetBorderSize(0); // turn off border wlg->SetFillStyle(0); f = theoryfiles1; ogr=0; while ( f && f->filename ) { ogr=OneGraph(f); if (ogr) { gr = fromGMnGDtoQ4GMn(ogr); TGraphAsymmErrors *egr = dynamic_cast<TGraphAsymmErrors*>(gr); if (egr && egr->GetEYhigh() && egr->GetEYhigh()[0]>0) { gr = toerror_band(egr); gr->SetFillStyle(3000+f->style); } gr->SetLineWidth(2); if (f->lnpt) { wgr->Add(gr,f->lnpt); wlg->AddEntry(gr,f->label,f->lnpt); } else if (gr->GetMarkerStyle()>=20) { wgr->Add(gr,"p"); wlg->AddEntry(gr,f->label,"p"); } else { wgr->Add(gr,"l"); wlg->AddEntry(gr,f->label,"l"); } } f++; } mgrThry->Draw("c"); legThry->Draw(); // draw a line at 1 cn->Modified(); cn->Update(); TFrame* frame = gPad->GetFrame(); TLine *line1 = new TLine(frame->GetX1(),1.,frame->GetX2(),1.); line1->SetLineStyle(1); line1->Draw(); cn->Update(); cn->SaveAs(Form("%s.eps",psfile)); cn->SaveAs(Form("%s.root",psfile)); gSystem->Exec(Form("./replace_symbols.pl %s.eps",psfile)); return; // remove everything but the graph cn->Update(); TList *clist = cn->GetListOfPrimitives(); frame = cn->GetFrame(); for (int i=0; i<clist->GetSize(); ) { if (clist->At(i) != frame) { clist->RemoveAt(i); } else i++; } // draw markers in the corners TMarker *mkr = new TMarker(frame->GetX1(),frame->GetY1(),2); mkr->Draw(); mkr = new TMarker(frame->GetX2(),frame->GetY1(),2); mkr->Draw(); mkr = new TMarker(frame->GetX1(),frame->GetY2(),2); mkr->Draw(); mkr = new TMarker(frame->GetX2(),frame->GetY2(),2); mkr->Draw(); frame->SetLineColor(10); cn->Update(); { datafile_t miller = { "figure_input/Miller/lattice.GMn.rtf","Miller", "[0]","[1]/(-1.91)*(1+[0]/.71)*(1+[0]/.71)", "([3]-[1])/(-1.91)*(1+[0]/.71)*(1+[0]/.71)", "([1]-[2])/(-1.91)*(1+[0]/.71)*(1+[0]/.71)", 0,0,1,3,"F" }; gr = OneGraph(&miller); TGraphAsymmErrors* egr = dynamic_cast<TGraphAsymmErrors*>(gr); if (egr && egr->GetEYhigh() && egr->GetEYhigh()[egr->GetN()/2]>0) { gr = toerror_band(egr); gr->SetLineStyle(1); gr->SetFillColor(gr->GetLineColor()); gr->SetFillStyle(3000+miller.style); } gr->Draw("F"); cn->Update(); cn->SaveAs("gmngd_Miller_Overlay.eps"); cn->SaveAs("gmngd_Miller_Overlay.root"); clist->Remove(gr); } { datafile_t miller_spl = { "figure_input/Miller/lattice.GMn.spl3","Miller", "[0]","[1]/(-1.91)*(1+[0]/.71)*(1+[0]/.71)", "([3]-[1])/(-1.91)*(1+[0]/.71)*(1+[0]/.71)", "([1]-[2])/(-1.91)*(1+[0]/.71)*(1+[0]/.71)", 0,0,1,3,"F" }; gr = OneGraph(&miller_spl); TGraphAsymmErrors* egr = dynamic_cast<TGraphAsymmErrors*>(gr); if (egr && egr->GetEYhigh() && egr->GetEYhigh()[egr->GetN()/2]>0) { gr = toerror_band(egr); gr->SetLineStyle(1); gr->SetFillColor(gr->GetLineColor()); gr->SetFillStyle(3000+miller_spl.style); } gr->Draw("F"); cn->Update(); cn->SaveAs("gmngd_Miller_Overlay_spl.eps"); cn->SaveAs("gmngd_Miller_Overlay_spl.root"); clist->Remove(gr); } }
void drawEventZR(double vrlb, double vzlb, double vrl0, double vzl0, double pmu1, double etamu1, double pmu2, double etamu2, double ppr, double etapr, double ppi, double etapi, int colors = 0) { Color_t colMu1(1), colMu2(1), colPr(1), colPi(1), colL0(1), colLb(1); switch(colors) { case 0: colLb=1; colL0=2; colMu1=2; colMu2=2; colPr=3; colPi=4; break; case 1: colLb=11; colL0=50; colMu1=50; colMu2=50; colPr=8; colPi=9; break; } cout << etamu1 << " " << etamu2 << " " << etapr << " " << etapi << endl; const double thetamu1 = 2*TMath::ATan(TMath::Exp(-etamu1)); const double thetamu2 = 2*TMath::ATan(TMath::Exp(-etamu2)); const double thetapr = 2*TMath::ATan(TMath::Exp(-etapr)); const double thetapi = 2*TMath::ATan(TMath::Exp(-etapi)); //const double thetamu1 = 2*TMath::ATan(TMath::Exp(-TMath::Abs(etamu1)))*sign(etamu1); //const double thetamu2 = 2*TMath::ATan(TMath::Exp(-TMath::Abs(etamu2)))*sign(etamu2); //const double thetapr = 2*TMath::ATan(TMath::Exp(-TMath::Abs(etapr)))*sign(etapr); //const double thetapi = 2*TMath::ATan(TMath::Exp(-TMath::Abs(etapi)))*sign(etapi); // draw the vertices TMarker *m; const double xlb = vzlb; const double ylb = vrlb; m = new TMarker(xlb,ylb,7); m->SetMarkerColor(colLb); m->Draw(); const double xl0 = vzl0; const double yl0 = vrl0; m = new TMarker(xl0,yl0,7); m->SetMarkerColor(colL0); m->Draw(); // draw the l0 flight line TLine *l; l = new TLine(vzlb,vrlb,vzl0,vrl0); l->SetLineColor(colL0); l->SetLineStyle(2); l->Draw(); // draw the muons TArrow *a; const double xmu1 = xlb + scalemu * pmu1 * TMath::Cos(thetamu1); const double ymu1 = ylb + scalemu * pmu1 * TMath::Sin(thetamu1); a = new TArrow(xlb,ylb,xmu1,ymu1,.01,">"); a->SetLineColor(colMu1); a->Draw(); const double xmu2 = xlb + scalemu * pmu2 * TMath::Cos(thetamu2); const double ymu2 = ylb + scalemu * pmu2 * TMath::Sin(thetamu2); a = new TArrow(xlb,ylb,xmu2,ymu2,.01,">"); a->SetLineColor(colMu2); a->Draw(); // draw the p and pi const double xpr = xl0 + scalepr * ppr * TMath::Cos(thetapr); const double ypr = yl0 + scalepr * ppr * TMath::Sin(thetapr); a = new TArrow(xl0,yl0,xpr,ypr,.01,">"); a->SetLineColor(colPr); a->Draw(); const double xpi = xl0 + scalepi * ppi * TMath::Cos(thetapi); const double ypi = yl0 + scalepi * ppi * TMath::Sin(thetapi); a = new TArrow(xl0,yl0,xpi,ypi,.01,">"); a->SetLineColor(colPi); a->Draw(); }
void plotAnitaEventMap(Adu5Pat *patPtr,double phi,double theta){ double sourceLon,sourceLat,headLon,headLat,phi10Lon,phi10Lat,phi6Lon,phi6Lat,phi14Lon,phi14Lat,actualLat,actualLon,actual2Lat,actual2Lon; float xEvent,yEvent,xAnita,yAnita,anitaLat,anitaLon,anitaAlt,xHead,yHead,x10,y10,x6,y6,x14,y14,yActual,xActual,yActual2,xActual2; anitaLat = patPtr->latitude; anitaLon = patPtr->longitude; anitaAlt = patPtr->altitude; UsefulAdu5Pat usefulPat(patPtr); std::cout << "source " << std::endl; //int sourceLoc = usefulPat.getSourceLonAndLatAltZero((180-phi)/180.*PI,(theta)/180.*PI,sourceLon,sourceLat); int sourceLoc = usefulPat.getSourceLonAndLatAltZero((phi)/180.*PI,(theta)/180.*PI,sourceLon,sourceLat); std::cout << std::endl << "heading " << std::endl; int headLoc = usefulPat.getSourceLonAndLatAltZero(0./180.*PI,10./180.*PI,headLon,headLat); std::cout << std::endl << "phi 10 " << std::endl; int headLoc10 = usefulPat.getSourceLonAndLatAltZero(180./180.*PI,10./180.*PI,phi10Lon,phi10Lat); std::cout << std::endl << "phi 14 " << std::endl; int headLoc14 = usefulPat.getSourceLonAndLatAltZero(270./180.*PI,10./180.*PI,phi14Lon,phi14Lat); std::cout << std::endl << "phi 6 " << std::endl; int headLoc6 = usefulPat.getSourceLonAndLatAltZero(90./180.*PI,10./180.*PI,phi6Lon,phi6Lat); std::cout << std::endl << "actual 14.5 " << std::endl; int actualLoc = usefulPat.getSourceLonAndLatAltZero(phi/180.*PI,12.5/180.*PI,actualLon,actualLat); std::cout << std::endl << "actual 4.5 " << std::endl; int actualLoc2 = usefulPat.getSourceLonAndLatAltZero(phi/180.*PI,7.5/180.*PI,actual2Lon,actual2Lat); //int sourceLoc = usefulPat.getSourceLonAndLatAltZero((phi)/180.*PI,(TMath::PiOver2()-theta)/180.*PI,sourceLon,sourceLat); TImage *map = TImage::Open("/home/mottram/work/eventCorrelator/macros/antarcticaIceMap.png"); std::cout << "sourceLoc " << sourceLoc << " phi " << phi << " theta " << theta << " lon " << sourceLon << " lat " << sourceLat << std::endl; gStyle->SetMarkerColor(kBlack); //gStyle->SetMarkerSize(2); gStyle->SetTextSize(0.02); TMarker *anitaPos = new TMarker(xAnita,yAnita,23); getRelXYFromLatLong(anitaLat,anitaLon,xAnita,yAnita); getRelXYFromLatLong(static_cast<float>(sourceLat),static_cast<float>(sourceLon),xEvent,yEvent); getRelXYFromLatLong(static_cast<float>(headLat),static_cast<float>(headLon),xHead,yHead); getRelXYFromLatLong(static_cast<float>(phi10Lat),static_cast<float>(phi10Lon),x10,y10); getRelXYFromLatLong(static_cast<float>(phi14Lat),static_cast<float>(phi14Lon),x14,y14); getRelXYFromLatLong(static_cast<float>(phi6Lat),static_cast<float>(phi6Lon),x6,y6); getRelXYFromLatLong(static_cast<float>(actualLat),static_cast<float>(actualLon),xActual,yActual); getRelXYFromLatLong(static_cast<float>(actual2Lat),static_cast<float>(actual2Lon),xActual2,yActual2); TCanvas *canMap=(TCanvas*)gROOT->FindObject("canMap"); if(!canMap) canMap = new TCanvas("canMap","canMap",(int)xSize,(int)ySize); canMap->Clear(); canMap->SetLogz(); canMap->SetTopMargin(0); canMap->SetBottomMargin(0); canMap->SetLeftMargin(0); canMap->SetRightMargin(0); map->Draw(""); TMarker *headingPos = new TMarker(xHead,yHead,29); TMarker *heading14Pos = new TMarker(x14,y14,29); TMarker *heading10Pos = new TMarker(x10,y10,29); TMarker *heading6Pos = new TMarker(x6,y6,29); TMarker *actualPos = new TMarker(xActual,yActual,29); TMarker *actual2Pos = new TMarker(xActual2,yActual2,29); headingPos->SetMarkerColor(kRed); heading14Pos->SetMarkerColor(kGray); heading10Pos->SetMarkerColor(kGray+2); heading6Pos->SetMarkerColor(kViolet); actualPos->SetMarkerColor(kRed+2);//12.5 theta actual2Pos->SetMarkerColor(kBlue+2);//7.5 theta headingPos->Draw(""); heading14Pos->Draw(""); heading10Pos->Draw(""); heading6Pos->Draw(""); actualPos->Draw(""); actual2Pos->Draw(""); anitaPos->DrawMarker(xAnita,yAnita); TLatex *positionLabel=0; char label[FILENAME_MAX]; if(sourceLoc==0){ if(anitaAlt<0){ sprintf(label,"Could not get event position, ANITA below 0 altitude!"); } else if(theta>0){ sprintf(label,"Pointing upwards! Cannot locate source at ground position"); } else{ sprintf(label,"Unkown error, cannot position source at 0 altitude"); } positionLabel = new TLatex(); positionLabel->DrawLatex(0.05,0.95,label); return; } TMarker *eventPos = new TMarker(xEvent,yEvent,29); eventPos->SetMarkerColor(kBlack); eventPos->Draw(""); sprintf(label,"ANITA location: lat %f; long %f; alt %f, x %f, y %f",anitaLat,anitaLon,anitaAlt,xAnita,yAnita); positionLabel = new TLatex(); positionLabel->DrawLatex(0.05,0.97,label); sprintf(label,"Event location: lat %f; long %f, x %f, y %f",sourceLat,sourceLon,xEvent,yEvent); positionLabel->DrawLatex(0.05,0.94,label); }
void PrintVertexRosenbrock() { // open file std::ifstream data ("Vertex.txt", std::ios::in); double buf; std::vector<double> x1; std::vector<double> x2; std::vector<double> x3; std::vector<double> y1; std::vector<double> y2; std::vector<double> y3; // loop on file while (1) { data >> buf; if (data.eof()) {break;} x1.push_back(buf); data >> buf; y1.push_back(buf); data >> buf; x2.push_back(buf); data >> buf; y2.push_back(buf); data >> buf; x3.push_back(buf); data >> buf; y3.push_back(buf); } // functions //TF2* f = new TF2 ("himmelblau", "(x*x + y - 11)*(x*x + y - 11) + (x + y*y - 7)*(x + y*y - 7)", 0, 10, 0, 10); // himmelblau gStyle->SetOptTitle(0); TF2* f = new TF2 ("rosenbrock", "100*(y-x*x)*(y-x*x) + (1-x)*(1-x)", -1.5, 2.2, -0.5, 2.2); // rosenbrock f->SetNpx(1000); f->SetNpy(1000); f->GetXaxis()->SetTitle("x"); f->GetYaxis()->SetTitle("y"); TCanvas* c1 = new TCanvas; c1->SetLogz(); rosenbrock->Draw("COLZ"); // loop on vertex std::vector<TLine*> vl1; std::vector<TLine*> vl2; std::vector<TLine*> vl3; TLine* l1; TLine* l2; TLine* l3; for (int i = 0; i < x1.size(); i++) { l1 = new TLine(0,0,0,0); l2 = new TLine(0,0,0,0); l3 = new TLine(0,0,0,0); l1->SetLineWidth(1.5); l2->SetLineWidth(1.5); l3->SetLineWidth(1.5); l1->SetLineColor(kBlue -10 + i); l2->SetLineColor(kBlue -10 + i); l3->SetLineColor(kBlue -10 + i); l1->SetX1(x1.at(i)); l1->SetY1(y1.at(i)); l1->SetX2(x2.at(i)); l1->SetY2(y2.at(i)); l2->SetX1(x2.at(i)); l2->SetY1(y2.at(i)); l2->SetX2(x3.at(i)); l2->SetY2(y3.at(i)); l3->SetX1(x3.at(i)); l3->SetY1(y3.at(i)); l3->SetX2(x1.at(i)); l3->SetY2(y1.at(i)); l1->Draw("same"); l2->Draw("same"); l3->Draw("same"); } TMarker* min = new TMarker (1,1,0); min ->SetMarkerColor (kBlue); min->SetMarkerStyle(20); min->SetMarkerSize (1.0); min->Draw(); }
void /home/users/m/k/mkomm/Diff13/analysis/unfolding/result/nominal/multi_top_y_nol() { //=========Macro generated from canvas: cvMulti0.984295447508/ //========= (Mon May 15 15:02:17 2017) by ROOT version6.02/05 TCanvas *cvMulti0.984295447508 = new TCanvas("cvMulti0.984295447508", "",0,0,800,650); gStyle->SetOptFit(1); gStyle->SetOptStat(0); gStyle->SetOptTitle(0); cvMulti0.984295447508->SetHighLightColor(2); cvMulti0.984295447508->Range(-0.5419355,-0.2962025,3.329032,1.982278); cvMulti0.984295447508->SetFillColor(0); cvMulti0.984295447508->SetBorderMode(0); cvMulti0.984295447508->SetBorderSize(2); cvMulti0.984295447508->SetTickx(1); cvMulti0.984295447508->SetTicky(1); cvMulti0.984295447508->SetLeftMargin(0.14); cvMulti0.984295447508->SetRightMargin(0.24); cvMulti0.984295447508->SetTopMargin(0.08); cvMulti0.984295447508->SetBottomMargin(0.13); cvMulti0.984295447508->SetFrameFillStyle(0); cvMulti0.984295447508->SetFrameLineWidth(0); cvMulti0.984295447508->SetFrameBorderMode(0); cvMulti0.984295447508->SetFrameBorderSize(0); cvMulti0.984295447508->SetFrameFillStyle(0); cvMulti0.984295447508->SetFrameLineWidth(0); cvMulti0.984295447508->SetFrameBorderMode(0); cvMulti0.984295447508->SetFrameBorderSize(0); TH2F *axis0.09336527505115 = new TH2F("axis0.09336527505115","",50,0,2.4,50,0,1.8); axis0.09336527505115->SetLineStyle(0); axis0.09336527505115->SetMarkerStyle(20); axis0.09336527505115->SetMarkerSize(0.16); axis0.09336527505115->GetXaxis()->SetTitle("top quark |y|"); axis0.09336527505115->GetXaxis()->SetNdivisions(1005); axis0.09336527505115->GetXaxis()->SetLabelFont(43); axis0.09336527505115->GetXaxis()->SetLabelOffset(0.0077); axis0.09336527505115->GetXaxis()->SetLabelSize(32); axis0.09336527505115->GetXaxis()->SetTitleSize(35); axis0.09336527505115->GetXaxis()->SetTickLength(0.04032258); axis0.09336527505115->GetXaxis()->SetTitleOffset(1.15); axis0.09336527505115->GetXaxis()->SetTitleFont(43); axis0.09336527505115->GetYaxis()->SetTitle(" Scale factor"); axis0.09336527505115->GetYaxis()->SetNdivisions(512); axis0.09336527505115->GetYaxis()->SetLabelFont(43); axis0.09336527505115->GetYaxis()->SetLabelOffset(0.0077); axis0.09336527505115->GetYaxis()->SetLabelSize(32); axis0.09336527505115->GetYaxis()->SetTitleSize(35); axis0.09336527505115->GetYaxis()->SetTickLength(0.03164557); axis0.09336527505115->GetYaxis()->SetTitleOffset(1.3); axis0.09336527505115->GetYaxis()->SetTitleFont(43); axis0.09336527505115->GetZaxis()->SetLabelFont(43); axis0.09336527505115->GetZaxis()->SetLabelOffset(0.0077); axis0.09336527505115->GetZaxis()->SetLabelSize(32); axis0.09336527505115->GetZaxis()->SetTitleSize(34); axis0.09336527505115->GetZaxis()->SetTitleOffset(1.3); axis0.09336527505115->GetZaxis()->SetTitleFont(43); axis0.09336527505115->Draw("AXIS"); TF1 *tf13 = new TF1("tf1","1",0,2.4); tf13->SetFillColor(19); tf13->SetFillStyle(0); tf13->SetMarkerStyle(20); tf13->SetMarkerSize(0.16); tf13->SetLineWidth(1); tf13->GetXaxis()->SetNdivisions(505); tf13->GetXaxis()->SetLabelFont(43); tf13->GetXaxis()->SetLabelOffset(0.0077); tf13->GetXaxis()->SetLabelSize(32); tf13->GetXaxis()->SetTitleSize(34); tf13->GetXaxis()->SetTickLength(0.05); tf13->GetXaxis()->SetTitleOffset(1.15); tf13->GetXaxis()->SetTitleFont(43); tf13->GetYaxis()->SetNdivisions(512); tf13->GetYaxis()->SetLabelFont(43); tf13->GetYaxis()->SetLabelOffset(0.0077); tf13->GetYaxis()->SetLabelSize(32); tf13->GetYaxis()->SetTitleSize(34); tf13->GetYaxis()->SetTitleOffset(1.3); tf13->GetYaxis()->SetTitleFont(43); tf13->Draw("LSame"); TMarker *marker = new TMarker(0.0892,0.9606917,20); Int_t ci; // for color index setting TColor *color; // for color definition with alpha ci = TColor::GetColor("#f9910c"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); TLine *line = new TLine(0.0892,0.9359605,0.0892,0.9854228); ci = TColor::GetColor("#f9910c"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(0.2,0,0.2,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(0.0892,0.9606917,0.3142,0.9554859); ci = TColor::GetColor("#f9910c"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(0.3142,0.9554859,20); ci = TColor::GetColor("#f9910c"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(0.3142,0.9323344,0.3142,0.9786374); ci = TColor::GetColor("#f9910c"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(0.45,0,0.45,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(0.3142,0.9554859,0.5642,0.9562376); ci = TColor::GetColor("#f9910c"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(0.5642,0.9562376,20); ci = TColor::GetColor("#f9910c"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(0.5642,0.9291835,0.5642,0.9832918); ci = TColor::GetColor("#f9910c"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(0.7,0,0.7,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(0.5642,0.9562376,0.8142,0.940729); ci = TColor::GetColor("#f9910c"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(0.8142,0.940729,20); ci = TColor::GetColor("#f9910c"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(0.8142,0.9117791,0.8142,0.9696789); ci = TColor::GetColor("#f9910c"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(0.95,0,0.95,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(0.8142,0.940729,1.0642,1.028565); ci = TColor::GetColor("#f9910c"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(1.0642,1.028565,20); ci = TColor::GetColor("#f9910c"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(1.0642,0.9987074,1.0642,1.058422); ci = TColor::GetColor("#f9910c"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(1.2,0,1.2,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(1.0642,1.028565,1.3392,1.024151); ci = TColor::GetColor("#f9910c"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(1.3392,1.024151,20); ci = TColor::GetColor("#f9910c"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(1.3392,0.9916412,1.3392,1.05666); ci = TColor::GetColor("#f9910c"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(1.5,0,1.5,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(1.3392,1.024151,1.5892,0.9990193); ci = TColor::GetColor("#f9910c"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(1.5892,0.9990193,20); ci = TColor::GetColor("#f9910c"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(1.5892,0.9494996,1.5892,1.048539); ci = TColor::GetColor("#f9910c"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(1.7,0,1.7,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(1.5892,0.9990193,2.0392,0.9212695); ci = TColor::GetColor("#f9910c"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(2.0392,0.9212695,20); ci = TColor::GetColor("#f9910c"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(2.0392,0.8767255,2.0392,0.9658136); ci = TColor::GetColor("#f9910c"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(2.4,0,2.4,1.8); line->SetLineStyle(2); line->Draw(); marker = new TMarker(0.1036,1.197582,20); ci = TColor::GetColor("#0a9128"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(0.1036,1.070044,0.1036,1.325121); ci = TColor::GetColor("#0a9128"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(0.2,0,0.2,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(0.1036,1.197582,0.3286,1.07941); ci = TColor::GetColor("#0a9128"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(0.3286,1.07941,20); ci = TColor::GetColor("#0a9128"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(0.3286,0.9690778,0.3286,1.189742); ci = TColor::GetColor("#0a9128"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(0.45,0,0.45,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(0.3286,1.07941,0.5786,1.238685); ci = TColor::GetColor("#0a9128"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(0.5786,1.238685,20); ci = TColor::GetColor("#0a9128"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(0.5786,1.103282,0.5786,1.374089); ci = TColor::GetColor("#0a9128"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(0.7,0,0.7,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(0.5786,1.238685,0.8286,1.070291); ci = TColor::GetColor("#0a9128"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(0.8286,1.070291,20); ci = TColor::GetColor("#0a9128"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(0.8286,0.9410973,0.8286,1.199485); ci = TColor::GetColor("#0a9128"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(0.95,0,0.95,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(0.8286,1.070291,1.0786,0.7963992); ci = TColor::GetColor("#0a9128"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(1.0786,0.7963992,20); ci = TColor::GetColor("#0a9128"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(1.0786,0.6813629,1.0786,0.9114354); ci = TColor::GetColor("#0a9128"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(1.2,0,1.2,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(1.0786,0.7963992,1.3536,0.9233731); ci = TColor::GetColor("#0a9128"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(1.3536,0.9233731,20); ci = TColor::GetColor("#0a9128"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(1.3536,0.7879375,1.3536,1.058809); ci = TColor::GetColor("#0a9128"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(1.5,0,1.5,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(1.3536,0.9233731,1.6036,1.066334); ci = TColor::GetColor("#0a9128"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(1.6036,1.066334,20); ci = TColor::GetColor("#0a9128"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(1.6036,0.8539173,1.6036,1.278751); ci = TColor::GetColor("#0a9128"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(1.7,0,1.7,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(1.6036,1.066334,2.0536,1.283384); ci = TColor::GetColor("#0a9128"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(2.0536,1.283384,20); ci = TColor::GetColor("#0a9128"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(2.0536,1.067281,2.0536,1.499488); ci = TColor::GetColor("#0a9128"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(2.4,0,2.4,1.8); line->SetLineStyle(2); line->Draw(); marker = new TMarker(0.118,0.1968883,20); ci = TColor::GetColor("#999999"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(0.118,0.1722209,0.118,0.2215558); ci = TColor::GetColor("#999999"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(0.2,0,0.2,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(0.118,0.1968883,0.343,0.2348717); ci = TColor::GetColor("#999999"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(0.343,0.2348717,20); ci = TColor::GetColor("#999999"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(0.343,0.2130683,0.343,0.2566752); ci = TColor::GetColor("#999999"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(0.45,0,0.45,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(0.343,0.2348717,0.593,0.2059685); ci = TColor::GetColor("#999999"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(0.593,0.2059685,20); ci = TColor::GetColor("#999999"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(0.593,0.1802414,0.593,0.2316957); ci = TColor::GetColor("#999999"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(0.7,0,0.7,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(0.593,0.2059685,0.843,0.2345659); ci = TColor::GetColor("#999999"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(0.843,0.2345659,20); ci = TColor::GetColor("#999999"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(0.843,0.208081,0.843,0.2610507); ci = TColor::GetColor("#999999"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(0.95,0,0.95,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(0.843,0.2345659,1.093,0.3094085); ci = TColor::GetColor("#999999"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(1.093,0.3094085,20); ci = TColor::GetColor("#999999"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(1.093,0.2799161,1.093,0.3389009); ci = TColor::GetColor("#999999"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(1.2,0,1.2,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(1.093,0.3094085,1.368,0.347837); ci = TColor::GetColor("#999999"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(1.368,0.347837,20); ci = TColor::GetColor("#999999"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(1.368,0.2999248,1.368,0.3957492); ci = TColor::GetColor("#999999"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(1.5,0,1.5,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(1.368,0.347837,1.618,0.5309299); ci = TColor::GetColor("#999999"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(1.618,0.5309299,20); ci = TColor::GetColor("#999999"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(1.618,0.4081874,1.618,0.6536723); ci = TColor::GetColor("#999999"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(1.7,0,1.7,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(1.618,0.5309299,2.068,0.2949066); ci = TColor::GetColor("#999999"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(2.068,0.2949066,20); ci = TColor::GetColor("#999999"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(2.068,0,2.068,0.7323252); ci = TColor::GetColor("#999999"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(2.4,0,2.4,1.8); line->SetLineStyle(2); line->Draw(); marker = new TMarker(0.1252,0.9784018,20); ci = TColor::GetColor("#ff0505"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(0.1252,0.8521009,0.1252,1.104703); ci = TColor::GetColor("#ff0505"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(0.2,0,0.2,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(0.1252,0.9784018,0.3502,0.9447267); ci = TColor::GetColor("#ff0505"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(0.3502,0.9447267,20); ci = TColor::GetColor("#ff0505"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(0.3502,0.8411522,0.3502,1.048301); ci = TColor::GetColor("#ff0505"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(0.45,0,0.45,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(0.3502,0.9447267,0.6002,0.7547158); ci = TColor::GetColor("#ff0505"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(0.6002,0.7547158,20); ci = TColor::GetColor("#ff0505"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(0.6002,0.6180429,0.6002,0.8913887); ci = TColor::GetColor("#ff0505"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(0.7,0,0.7,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(0.6002,0.7547158,0.8502,1.12315); ci = TColor::GetColor("#ff0505"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(0.8502,1.12315,20); ci = TColor::GetColor("#ff0505"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(0.8502,0.9840732,0.8502,1.262226); ci = TColor::GetColor("#ff0505"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(0.95,0,0.95,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(0.8502,1.12315,1.1002,1.107947); ci = TColor::GetColor("#ff0505"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(1.1002,1.107947,20); ci = TColor::GetColor("#ff0505"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(1.1002,0.9620515,1.1002,1.253843); ci = TColor::GetColor("#ff0505"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(1.2,0,1.2,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(1.1002,1.107947,1.3752,0.6245485); ci = TColor::GetColor("#ff0505"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(1.3752,0.6245485,20); ci = TColor::GetColor("#ff0505"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(1.3752,0.462954,1.3752,0.786143); ci = TColor::GetColor("#ff0505"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(1.5,0,1.5,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(1.3752,0.6245485,1.6252,0.8379717); ci = TColor::GetColor("#ff0505"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(1.6252,0.8379717,20); ci = TColor::GetColor("#ff0505"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(1.6252,0.5506004,1.6252,1.125343); ci = TColor::GetColor("#ff0505"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(1.7,0,1.7,1.8); line->SetLineStyle(2); line->Draw(); line = new TLine(1.6252,0.8379717,2.0752,0.3802364); ci = TColor::GetColor("#ff0505"); line->SetLineColor(ci); line->SetLineWidth(3); line->Draw(); marker = new TMarker(2.0752,0.3802364,20); ci = TColor::GetColor("#ff0505"); marker->SetMarkerColor(ci); marker->SetMarkerStyle(20); marker->SetMarkerSize(1.2); marker->Draw(); line = new TLine(2.0752,0.06910222,2.0752,0.6913706); ci = TColor::GetColor("#ff0505"); line->SetLineColor(ci); line->SetLineWidth(2); line->Draw(); line = new TLine(2.4,0,2.4,1.8); line->SetLineStyle(2); line->Draw(); TH2F *axis0.0933652750511_copy6 = new TH2F("axis0.0933652750511_copy6","",50,0,2.4,50,0,1.8); axis0.0933652750511_copy6->SetDirectory(0); axis0.0933652750511_copy6->SetLineStyle(0); axis0.0933652750511_copy6->SetMarkerStyle(20); axis0.0933652750511_copy6->SetMarkerSize(0.16); axis0.0933652750511_copy6->GetXaxis()->SetTitle("top quark |y|"); axis0.0933652750511_copy6->GetXaxis()->SetNdivisions(1005); axis0.0933652750511_copy6->GetXaxis()->SetLabelFont(43); axis0.0933652750511_copy6->GetXaxis()->SetLabelOffset(0.0077); axis0.0933652750511_copy6->GetXaxis()->SetLabelSize(32); axis0.0933652750511_copy6->GetXaxis()->SetTitleSize(35); axis0.0933652750511_copy6->GetXaxis()->SetTickLength(0.04032258); axis0.0933652750511_copy6->GetXaxis()->SetTitleOffset(1.15); axis0.0933652750511_copy6->GetXaxis()->SetTitleFont(43); axis0.0933652750511_copy6->GetYaxis()->SetTitle(" Scale factor"); axis0.0933652750511_copy6->GetYaxis()->SetNdivisions(512); axis0.0933652750511_copy6->GetYaxis()->SetLabelFont(43); axis0.0933652750511_copy6->GetYaxis()->SetLabelOffset(0.0077); axis0.0933652750511_copy6->GetYaxis()->SetLabelSize(32); axis0.0933652750511_copy6->GetYaxis()->SetTitleSize(35); axis0.0933652750511_copy6->GetYaxis()->SetTickLength(0.03164557); axis0.0933652750511_copy6->GetYaxis()->SetTitleOffset(1.3); axis0.0933652750511_copy6->GetYaxis()->SetTitleFont(43); axis0.0933652750511_copy6->GetZaxis()->SetLabelFont(43); axis0.0933652750511_copy6->GetZaxis()->SetLabelOffset(0.0077); axis0.0933652750511_copy6->GetZaxis()->SetLabelSize(32); axis0.0933652750511_copy6->GetZaxis()->SetTitleSize(34); axis0.0933652750511_copy6->GetZaxis()->SetTitleOffset(1.3); axis0.0933652750511_copy6->GetZaxis()->SetTitleFont(43); axis0.0933652750511_copy6->Draw("sameaxis"); TPaveText *pt = new TPaveText(0.14,0.94,0.14,0.94,"brNDC"); pt->SetBorderSize(0); pt->SetFillColor(0); pt->SetTextAlign(11); pt->SetTextFont(63); pt->SetTextSize(34); TText *AText = pt->AddText("CMS"); pt->Draw(); pt = new TPaveText(0.24,0.94,0.24,0.94,"brNDC"); pt->SetBorderSize(0); pt->SetFillColor(0); pt->SetTextAlign(11); pt->SetTextFont(53); pt->SetTextSize(34); AText = pt->AddText("Preliminary"); pt->Draw(); cvMulti0.984295447508->Modified(); cvMulti0.984295447508->cd(); cvMulti0.984295447508->SetSelected(cvMulti0.984295447508); }
void drawMeasurement(int i, double m[5], char label[2][100], int aux[5], double vstep, TH2F* histo, TCanvas* canvas) { double lowY = (i+1)*vstep; double uppY = (i+2)*vstep; //double lowX = histo->GetBinLowEdge(1); //double uppX = histo->GetBinLowEdge(histo->GetNbinsX()) + // histo->GetBinWidth(histo->GetNbinsX()); double lowX = 0.25; double uppX = 2.80; double widthX = uppX - lowX; // y-range of the histogram is [0...1] double startX = lowX + 0.04*widthX; TPaveText* text = new TPaveText(startX, lowY, startX, uppY, "BR"); text->SetTextAlign(12); text->SetFillColor(aux[1]); text->SetTextColor(aux[0]); text->SetLineColor(1); text->SetBorderSize(0); TText* t0 = text->AddText(" "); t0->SetTextSize(0.08); t0->SetTextFont(aux[2]); TText* t1 = text->AddText(label[0]); t1->SetTextSize(0.08); t1->SetTextFont(aux[2]); TText* t2 = text->AddText(label[1]); t2->SetTextSize(0.08); t2->SetTextFont(aux[2]); text->Draw(); double ypos = 0.5*(lowY+uppY); double mean = m[0]; double nErr1 = m[1]; double pErr1 = m[2]; double nErr2 = sqrt(m[1]*m[1]+m[3]*m[3]); double pErr2 = sqrt(m[2]*m[2]+m[4]*m[4]); // draw TGraphAsymmErrors 1 (stat only) |---*---| TMarker* measurement = new TMarker(mean, ypos, aux[4]); measurement->SetMarkerColor(aux[0]); measurement->SetMarkerStyle(aux[4]); measurement->SetMarkerSize(1.5); //measurement->SetMarkerSize(1.75); measurement->Draw(); double vsizeErr1 = 0.09*vstep; TLine* l1 = new TLine(mean, ypos, mean-nErr1, ypos); l1->SetLineWidth(aux[3]); l1->SetLineColor(aux[0]); l1->Draw(); TLine* l2 = new TLine(mean, ypos, mean+pErr1, ypos); l2->SetLineWidth(aux[3]); l2->SetLineColor(aux[0]); l2->Draw(); TLine* l3 = new TLine(mean-nErr1, ypos-vsizeErr1, mean-nErr1, ypos+vsizeErr1); l3->SetLineWidth(aux[3]); l3->SetLineColor(aux[0]); l3->Draw(); TLine* l4 = new TLine(mean+pErr1, ypos-vsizeErr1, mean+pErr1, ypos+vsizeErr1); l4->SetLineWidth(aux[3]); l4->SetLineColor(aux[0]); l4->Draw(); // overlay TGraphAsymmErrors 2 (stat+syst) |----*-----| double vsizeErr2 = 0.12*vstep; TLine* l5 = new TLine(mean, ypos, mean-nErr2, ypos); l5->SetLineWidth(aux[3]); l5->SetLineColor(aux[0]); l5->Draw(); TLine* l6 = new TLine(mean, ypos, mean+pErr2, ypos); l6->SetLineWidth(aux[3]); l6->SetLineColor(aux[0]); l6->Draw(); TLine* l7 = new TLine(mean-nErr2, ypos-vsizeErr2, mean-nErr2, ypos+vsizeErr2); l7->SetLineWidth(aux[3]); l7->SetLineColor(aux[0]); l7->Draw(); TLine* l8 = new TLine(mean+pErr2, ypos-vsizeErr2, mean+pErr2, ypos+vsizeErr2); l8->SetLineWidth(aux[3]); l8->SetLineColor(aux[0]); l8->Draw(); // draw measurement label "XXX+/-YY+/-ZZ" TPaveText* num = new TPaveText(uppX-0.32*widthX, lowY, uppX-0.02*widthX, uppY, "BR"); num->SetTextAlign(12); num->SetFillColor(aux[1]); num->SetTextColor(aux[0]); num->SetLineColor(aux[0]); num->SetBorderSize(0); TString str; char s[100]; sprintf(s, "%4.2f#color[%d]{X}", m[0], aux[1]); str +=s; if (m[1]==m[2]) { // sym. stat. errors sprintf(s, "#pm %4.2f", m[1]); str +=s; } else { sprintf(s, "^{+%4.2f}", m[2]); str +=s; sprintf(s, "_{-#color[%d]{|}%4.2f}", aux[1], m[1]); str +=s; } str += " (stat)"; if (m[3]!=0.0 || m[4]!=0.0) { if (m[3]==m[4]) { // sym. syst. errors sprintf(s, "#color[%d]{X}#pm% 4.2f", aux[1], m[3]); str +=s; } else { sprintf(s, "#color[%d]{X}", aux[1]); str +=s; sprintf(s, "^{+%4.2f}", m[4]); str +=s; sprintf(s, "_{-#color[%d]{|}%4.2f}", aux[1], m[3]); str +=s; } } str += " (syst)"; TText* n0 = num->AddText(str); n0->SetTextFont(aux[2]); num->Draw(); return; }
void geugd_riordan() { gROOT->SetStyle("HALLA"); TCanvas *cn = new TCanvas("cn"); cn->Draw(); cn->UseCurrentStyle(); TH1F *frm = new TH1F("frm","",100,0.,12.); frm->GetXaxis()->SetTitle("Q^{2} [GeV^{2}]"); frm->GetXaxis()->CenterTitle(); frm->GetYaxis()->SetTitle("G_{E}^{u}/G_{D}"); frm->GetYaxis()->CenterTitle(); frm->SetMinimum(.40); frm->SetMaximum(2.4); //frm->SetMaximum(0.3); frm->UseCurrentStyle(); frm->Draw(); frm->SetAxisRange(0.10,12.,"X"); TF1* galster = new TF1("galster", "x/(4.*0.938*.938)*1.91/(1.+x/.71)^2/(1.+5.6*x/(4.*.938*.938))", 0.,4.); galster->SetLineColor(kBlack); galster->SetLineStyle(kBlack); galster->SetLineWidth(2); TF1* gen0 = new TF1("gen0", f1dugen0, 0.,12.); gen0->SetLineColor(kBlack); gen0->SetLineStyle(kBlack); gen0->SetLineWidth(1); TF1 *genf = new TF1("genf",genff,1.,10.,1); genf->SetLineColor(kBlue); genf->SetLineStyle(1); genf->SetParameter(0,1.); // match to Madey point just below 1.5 // genf->SetParameter(0,.0411/genf->Eval(1.45)); // genf->SetParameter(0,-0.558645); genf->SetParameter(0,-0.182645); TF1 *roberts_curve = new TF1("roberts",roberts_gen,0.035,12.344,1); roberts_curve->SetLineColor(kRed); roberts_curve->SetLineStyle(9); TF1 *ourfit = new TF1("ourfit",gen_ourfit,0.,10.,0); ourfit->SetLineColor(kBlue); ourfit->SetLineStyle(0); /* TF1 *bbba05 = new TF1("BBBA05",gen_bbba05,0.,10.,0); bbba05->SetLineColor(kGreen); bbba05->SetLineStyle(3); */ // TF1 *lomon = new TF1("Lomon",Lomon_GEn,0.,10.,0); // lomon->SetLineColor(7); // lomon->SetLineStyle(4); TMultiGraph* mgrDta = new TMultiGraph("Data","G_{E}^{n}"); //TLegend *legDta = new TLegend(.3448,.6123,.6810,.9110,"","brNDC"); TLegend *legDta = new TLegend(.6020,.4004,.9382,.9089,"","brNDC"); TMultiGraph* wgr = mgrDta; TLegend *wlg = legDta; // the data legDta->SetBorderSize(0); // turn off border legDta->SetFillStyle(0); datafile_t *f = datafiles; TGraph* gr=0; TGraph* ogr=0; while ( f && f->filename ) { ogr=OneGraph(f); if (ogr) { gr = fromGEntransform(ogr); gr->SetLineStyle(0); if (f->lnpt) { mgrDta->Add(gr,f->lnpt); if( f->label[0] != 'x' ) legDta->AddEntry(gr,f->label,f->lnpt); } mgrDta->Add(gr,"p"); if( f->label[0] != 'x' ) legDta->AddEntry(gr,f->label,"p"); /* else if (gr->GetMarkerStyle()>=20) { mgrDta->Add(gr,"p"); if( f->label[0] != 'x' ) legDta->AddEntry(gr,f->label,"p"); } else { mgrDta->Add(gr,"l"); if( f->label[0] != 'x' ) legDta->AddEntry(gr,f->label,"l"); } */ } f++; } mgrDta->Draw("p"); legDta->Draw(); TMultiGraph* mgrThry = new TMultiGraph("Theory","G_{E}^{n}"); // TLegend *legThry = new TLegend(.546,.6208,.8822,.9195,"","brNDC"); // TLegend *legThry = new TLegend(.2055,.7013,.6020,.8893,"","brNDC"); TLegend *legThry = new TLegend(.4267,.8453,.6236,.8962,"","brNDC"); wgr = mgrThry; wlg = legThry; // the theory wlg->SetBorderSize(0); // turn off border wlg->SetFillStyle(0); f = theoryfiles1; gr=0; Bool_t isfill; while ( f && f->filename ) { gr=OneGraph(f); gr->SetLineWidth(2); isfill = kFALSE; if (gr) { TGraphAsymmErrors *egr = dynamic_cast<TGraphAsymmErrors*>(gr); if (egr && egr->GetN()>1 && egr->GetEYhigh() && egr->GetEYhigh()[1]>0) { gr = toerror_band(egr); gr->SetFillStyle(3000+f->style); gr->SetFillStyle(1); gr->SetFillColor(gr->GetLineColor()); isfill = kTRUE; } if (f->lnpt) { wgr->Add(gr,f->lnpt); if( f->label[0] != 'x' ) wlg->AddEntry(gr,f->label,f->lnpt); } else if (gr->GetMarkerStyle()>=20) { wgr->Add(gr,"p"); if( f->label[0] != 'x' ) wlg->AddEntry(gr,f->label,"p"); } else { if( isfill ){ gr->SetLineStyle(0); wgr->Add(gr,"f"); // wgr->Add( new TGraph(*egr),"l"); if( f->label[0] != 'x' ) wlg->AddEntry(gr,f->label,"lf"); } else { wgr->Add(gr,"l"); if( f->label[0] != 'x' ) wlg->AddEntry(gr,f->label,"l"); } } } f++; } TPaveLabel *prelim = new TPaveLabel( 0.3204, 0.1886, 0.9310, 0.4643, "PRELIMINARY", "NDC" ); prelim->SetTextAngle(0); TColor *pink = new TColor(300, 1.0, 0.7, 0.7, "pink"); prelim->SetTextColor(300); prelim->SetFillStyle(0); prelim->SetBorderSize(0); #ifdef FAKE_SCHIAVILLA TGraph *schiagraph = draw_schiavilla(); schiagraph->Draw("F"); schiagraph = draw_schiavilla(1); schiagraph->Draw("C"); legThry->AddEntry(schiagraph, "d(e,e'd) T_{20} - Schiavilla & Sick", "LF"); #endif #ifdef PRELIMINARY prelim->Draw("same"); #endif TPaveLabel *nofsilab = new TPaveLabel( 0.4569, 0.1674, 0.7514, 0.2415, "No FSI Corrections", "NDC" ); nofsilab->SetFillStyle(0); nofsilab->SetBorderSize(0); nofsilab->SetTextColor(kRed); #ifdef NOFSI // nofsilab->Draw("same"); #endif #ifdef PQCD genf->Draw("same"); #endif mgrThry->Draw("c"); #ifdef GALSTER galster->Draw("same"); #endif // gen0->Draw("same"); #ifdef NEW_ROBERTS roberts_curve->Draw("same"); #endif // ourfit->Draw("same"); // bbba05->Draw("same"); //lomon->Draw("same"); // bandi->Draw("same"); #ifdef PQCD legThry->AddEntry(genf,"F_{2}/F_{1} #propto ln^{2}(Q^{2}/#Lambda^{2})/Q^{2}","l"); #endif #ifdef GALSTER legThry->AddEntry(galster,"Galster fit (1971)","l"); #endif // legThry->AddEntry(gen0,"G_{E}^{n} = 0","l"); #ifdef NEW_ROBERTS legThry->AddEntry(roberts_curve,"q(qq) Faddeev - I. Cloet, ANL","l"); #endif // legThry->AddEntry(bbba05,"BBBA05","l"); // legThry->AddEntry(lomon, "Lomon", "l"); // legThry->AddEntry(ourfit, "Our Fit", "l"); legThry->Draw(); legDta->Draw(); mgrDta->Draw("p"); legDta->Draw(); TEllipse *cir1 = new TEllipse(1.31, 0.0, 0.17, 0.0065); TEllipse *cir2 = new TEllipse(2.4, 0.0, 0.17, 0.0065); cir1->SetFillStyle(0); cir2->SetFillStyle(0); #ifdef PLOTKIN1 cir1->Draw("same"); // cir2->Draw("same"); // TPaveLabel *exacc = new TPaveLabel( 0.3779, 0.2691, 0.5963, 0.3114, "Expected Accuracy", "NDC"); TPaveLabel *exacc = new TPaveLabel( 0.2716, 0.2733, 0.4899, 0.3136, "Expected Accuracy", "NDC"); exacc->SetFillStyle(0); exacc->SetBorderSize(0); exacc->Draw("same"); #endif // draw a line at 1 cn->Modified(); cn->Update(); cn->SaveAs(Form("%s.pdf",psfile)); cn->SaveAs(Form("%s.eps",psfile)); // cn->SaveAs(Form("%s.root",psfile)); // gSystem->Exec(Form("./replace_symbols.pl %s.eps",psfile)); // cn->SaveAs("bogdan_gen.eps"); return; // LEAVING HERE // now an overlay, hopefully matching dimensions // remove everything but the graph cn->Update(); TList *clist = cn->GetListOfPrimitives(); TFrame* frame = cn->GetFrame(); for (int i=0; i<clist->GetSize(); ) { if (clist->At(i) != frame) { clist->RemoveAt(i); } else i++; } // draw markers in the corners TMarker *mkr = new TMarker(frame->GetX1(),frame->GetY1(),2); mkr->Draw(); mkr = new TMarker(frame->GetX2(),frame->GetY1(),2); mkr->Draw(); mkr = new TMarker(frame->GetX1(),frame->GetY2(),2); mkr->Draw(); mkr = new TMarker(frame->GetX2(),frame->GetY2(),2); mkr->Draw(); frame->SetLineColor(10); cn->Update(); datafile_t miller = { "figure_input/Miller/lattice.GEn.rtf","Miller", "[0]","[1]","[1]-[3]","[2]-[1]",0,0,1,3,"F" }; gr = OneGraph(&miller); TGraphAsymmErrors* egr = dynamic_cast<TGraphAsymmErrors*>(gr); if (egr && egr->GetEYhigh() && egr->GetEYhigh()[egr->GetN()/2]>0) { gr = toerror_band(egr); gr->SetLineStyle(1); gr->SetFillColor(gr->GetLineColor()); gr->SetFillStyle(3000+miller.style); } gr->Draw("F"); cn->Update(); cn->SaveAs("gen_Miller_Overlay.pdf"); cn->SaveAs("gen_Miller_Overlay.root"); }
void showHistogram2d(TH2* histogram, const std::string& xAxisTitle, const std::string& yAxisTitle, int zAxisNormOption, double zMin, double zMax, Float_t* genX, Float_t* genY, const std::string& outputFileName) { TCanvas* canvas = new TCanvas("canvas", "canvas", 900, 800); canvas->SetFillColor(10); canvas->SetBorderSize(2); canvas->SetTopMargin(0.10); canvas->SetLeftMargin(0.12); canvas->SetRightMargin(0.14); canvas->SetBottomMargin(0.12); histogram->SetTitle(""); histogram->SetStats(false); int numBinsX = histogram->GetNbinsX(); int numBinsY = histogram->GetNbinsY(); if ( zAxisNormOption == kNormByQuantiles ) { std::vector<double> binContents; for ( int iBinX = 1; iBinX <= numBinsX; ++iBinX ) { for ( int iBinY = 1; iBinY <= numBinsY; ++iBinY ) { binContents.push_back(histogram->GetBinContent(iBinX, iBinY)); } } std::sort(binContents.begin(), binContents.end()); histogram->SetMinimum(binContents[TMath::Nint(0.05*binContents.size())]); histogram->SetMaximum(binContents[TMath::Nint(0.95*binContents.size())]); } else if ( zAxisNormOption == kNormByNegLogMax ) { double maxBinContent = 0.; for ( int iBinX = 1; iBinX <= numBinsX; ++iBinX ) { for ( int iBinY = 1; iBinY <= numBinsY; ++iBinY ) { double binContent = histogram->GetBinContent(iBinX, iBinY); if ( binContent > maxBinContent ) { std::cout << "binX = " << iBinX << " (x = " << histogram->GetXaxis()->GetBinCenter(iBinX) << ")," << " binY = " << iBinY << " (y = " << histogram->GetYaxis()->GetBinCenter(iBinY) << "): maxBinContent = " << maxBinContent << std::endl; maxBinContent = binContent; } } } double logMaxBinContent = TMath::Log(maxBinContent); for ( int iBinX = 1; iBinX <= numBinsX; ++iBinX ) { for ( int iBinY = 1; iBinY <= numBinsY; ++iBinY ) { double binContent = histogram->GetBinContent(iBinX, iBinY); if ( binContent > 0. ) { histogram->SetBinContent(iBinX, iBinY, -TMath::Log(binContent) + logMaxBinContent); } else { histogram->SetBinContent(iBinX, iBinY, -1.); } } } histogram->SetMinimum(0.); histogram->SetMaximum(zMax); } else if ( zAxisNormOption == kNormByValue ) { histogram->SetMinimum(zMin); histogram->SetMaximum(zMax); } else assert(0); TAxis* xAxis = histogram->GetXaxis(); xAxis->SetTitle(xAxisTitle.data()); xAxis->SetTitleOffset(1.15); TAxis* yAxis = histogram->GetYaxis(); yAxis->SetTitle(yAxisTitle.data()); yAxis->SetTitleOffset(1.30); gStyle->SetPalette(1,0); histogram->Draw("COLZ"); TMarker* genMarker = 0; if ( genX && genY ) { genMarker = new TMarker(*genX, *genY, 34); genMarker->SetMarkerColor(1); genMarker->SetMarkerSize(2); genMarker->Draw(); } canvas->Update(); size_t idx = outputFileName.find_last_of('.'); std::string outputFileName_plot = std::string(outputFileName, 0, idx); if ( idx != std::string::npos ) canvas->Print(std::string(outputFileName_plot).append(std::string(outputFileName, idx)).data()); canvas->Print(std::string(outputFileName_plot).append(".png").data()); canvas->Print(std::string(outputFileName_plot).append(".pdf").data()); canvas->Print(std::string(outputFileName_plot).append(".root").data()); delete genMarker; delete canvas; }
void kees_gen() { gROOT->SetStyle("HALLA"); TCanvas *cn = new TCanvas("cn"); cn->Draw(); cn->UseCurrentStyle(); TH1F *frm = new TH1F("frm","",100,0.,10.); frm->GetXaxis()->SetTitle("Q^{2} [GeV^{2}]"); frm->GetYaxis()->SetTitle("G_{E}^{n}"); frm->SetMinimum(-.02); frm->SetMaximum(0.1); frm->UseCurrentStyle(); frm->Draw(); frm->SetAxisRange(0.,5.,"X"); TF1 *genf = new TF1("genf",genff,1.,10.,1); genf->SetLineColor(2); genf->SetLineStyle(2); genf->SetParameter(0,1.); genf->SetParameter(1,.3); genf->SetParameter(0,-0.632); // match to Madey point just below 1.5 // genf->SetParameter(0,.0411/genf->Eval(1.45)); TMultiGraph* mgrDta = new TMultiGraph("Data","G_{E}^{n}"); TLegend *legDta = new TLegend(.54,.6,.875,.90,"","brNDC"); TMultiGraph* wgr = mgrDta; TLegend *wlg = legDta; // the data legDta->SetBorderSize(0); // turn off border legDta->SetFillStyle(0); datafile_t *f = datafiles; TGraph* gr=0; while ( f && f->filename ) { gr=OneGraph(f); if (gr) { if (f->lnpt) { mgrDta->Add(gr,f->lnpt); legDta->AddEntry(gr,f->label,f->lnpt); } else if (gr->GetMarkerStyle()>=20) { mgrDta->Add(gr,"p"); legDta->AddEntry(gr,f->label,"p"); } else { mgrDta->Add(gr,"l"); legDta->AddEntry(gr,f->label,"l"); } } f++; } mgrDta->Draw("p"); // legDta->Draw(); don't draw the data legend TMultiGraph* mgrThry = new TMultiGraph("Theory","G_{E}^{n}"); TLegend *legThry = new TLegend(.54,.6,.875,.9,"","brNDC"); wgr = mgrThry; wlg = legThry; // the theory wlg->SetBorderSize(0); // turn off border wlg->SetFillStyle(0); f = theoryfiles1; gr=0; while ( f && f->filename ) { gr=OneGraph(f); if (gr) { TGraphAsymmErrors *egr = dynamic_cast<TGraphAsymmErrors*>(gr); if (egr && egr->GetN()>1 && egr->GetEYhigh() && egr->GetEYhigh()[1]>0) { gr = toerror_band(egr); gr->SetFillStyle(3000+f->style); } if (f->lnpt) { wgr->Add(gr,f->lnpt); wlg->AddEntry(gr,f->label,f->lnpt); } else if (gr->GetMarkerStyle()>=20) { wgr->Add(gr,"p"); wlg->AddEntry(gr,f->label,"p"); } else { wgr->Add(gr,"l"); wlg->AddEntry(gr,f->label,"l"); } } f++; } genf->Draw("same"); mgrThry->Draw("c"); legThry->AddEntry(genf,"F_{2}/F_{1} #propto ln^{2}(Q^{2}/#Lambda^{2})/Q^{2}","l"); legThry->Draw(); // draw a line at 1 cn->Modified(); cn->Update(); cn->SaveAs(Form("%s.eps",psfile)); cn->SaveAs(Form("%s.root",psfile)); gSystem->Exec(Form("./replace_symbols.pl %s.eps",psfile)); // now an overlay, hopefully matching dimensions // remove everything but the graph cn->Update(); TList *clist = cn->GetListOfPrimitives(); TFrame* frame = cn->GetFrame(); for (int i=0; i<clist->GetSize(); ) { if (clist->At(i) != frame) { clist->RemoveAt(i); } else i++; } // draw markers in the corners TMarker *mkr = new TMarker(frame->GetX1(),frame->GetY1(),2); mkr->Draw(); mkr = new TMarker(frame->GetX2(),frame->GetY1(),2); mkr->Draw(); mkr = new TMarker(frame->GetX1(),frame->GetY2(),2); mkr->Draw(); mkr = new TMarker(frame->GetX2(),frame->GetY2(),2); mkr->Draw(); frame->SetLineColor(10); cn->Update(); datafile_t miller = { "figure_input/Miller/lattice.GEn.rtf","Miller", "[0]","[1]","[1]-[3]","[2]-[1]",0,0,1,3,"F" }; gr = OneGraph(&miller); TGraphAsymmErrors* egr = dynamic_cast<TGraphAsymmErrors*>(gr); if (egr && egr->GetEYhigh() && egr->GetEYhigh()[egr->GetN()/2]>0) { gr = toerror_band(egr); gr->SetLineStyle(1); gr->SetFillColor(gr->GetLineColor()); gr->SetFillStyle(3000+miller.style); } gr->Draw("F"); cn->Update(); cn->SaveAs("gen_Miller_Overlay.eps"); cn->SaveAs("gen_Miller_Overlay.root"); }
// draw the same thing but after reco void genPlots02(std::string fullPath, int nOverlay = 500, bool custBinning = false) { const int fVerbose(1); setTDRStyle(); gStyle->SetOptStat(112211); gStyle->SetPalette(1); // Canvas c = new TCanvas("c2","c2",1000,600); const unsigned int nPadX = 1; const unsigned int nPadY = 1; c->Divide(nPadX,nPadY); const unsigned int nPads=nPadX*nPadY; for(unsigned int i=1; i<=nPads; i++) { TPad* pad= (TPad*)c->cd(i); pad->SetTopMargin(0.10); pad->SetRightMargin(0.20); pad->SetLeftMargin(0.15); } // Open file TFile *f = TFile::Open(fullPath.c_str()); if (f==0) { cout << "File " << fullPath << " not found -- exiting" << endl; return; } if(fVerbose>0) cout << "Succesfully opened file " << fullPath << endl; // Get TTree TTree* t = (TTree*) f->Get("events"); if(fVerbose>0) cout << "Got TTree with " << t->GetEntries() << " entries" << endl; // Do a cut, if needed //t->Draw(">>lst","chi2lb>.1&&mlb>5.61&&mlb<5.63"); //t->Draw(">>lst","chi2lb>.1&&isSig==1"); t->Draw(">>lst","(rid1m&4)==4&&(rid2m&4)==4&&mjp>2.895&&mjp<3.295&&prob1m>0.1&&prob2m>0.1&&ptjp>2&&probjp>0.005&&ml0>1.101&&ml0<1.129&&probpr>0.02&&probpi>0.02&&rptpr>rptpi&&ptl0>3&&rptpr>1&&rptpi>0.5&&probl0>0.02&&alphal0<0.3&&d3l0>1&&d3l0/d3El0>10&&problb>0.001&&alphalb<0.3"); TEventList *lst; lst = (TEventList*)gDirectory->Get("lst"); t->SetEventList(lst); if(fVerbose>0) cout << "Got TTree with " << t->GetEntries() << " entries" << endl; // Do plots c->cd(1); //doPlot2d(t,"hrzL0vtx", "vrl0:TMath::Abs(vzl0)",30,0,300,30,0,120,"Tit","|z|","r","cm","cm"); if (custBinning) { double newbinsX[]={0,2,4,6,8,10,12,14,16,18,20,22,24,26,28,30,32,34,36,38,40,42,44,46,48,50}; const int newbinsX_size = sizeof(newbinsX)/sizeof(double); std::vector<double> binvecX(newbinsX,newbinsX+newbinsX_size); //double newbinsY[]={0,1,2,3,4,5,6,7,8,9,10}; double newbinsY[]={0,0.5,1,2,4,8,16,32}; const int newbinsY_size = sizeof(newbinsY)/sizeof(double); std::vector<double> binvecY(newbinsY,newbinsY+newbinsY_size); doPlot2d(t,"hrzL0vtxreco", "vrl0:TMath::Abs(vzl0)",binvecX, binvecY,"#Lambda vertices","|z|","r","cm","cm"); } else { doPlot2d(t,"hrzL0vtxreco", "vrl0:TMath::Abs(vzl0)",30,0,50,30,0,30,"#Lambda vertices","|z|","r","cm","cm"); } // add tracker TPad* pad; pad = (TPad*)c->cd(1); pad->Modified(); pad->Update(); repositionPalette("hrzL0vtxreco"); pad->Update(); pad->SetLogz(); drawTracker(pad); if (nOverlay<=0) return; int maxN = nOverlay; if (maxN > t->GetEntries()) maxN = t->GetEntries(); double vrl0,vzl0,ppr,ppi,etapr,etapi; t->SetBranchAddress("vrl0",&vrl0); t->SetBranchAddress("vzl0",&vzl0); t->SetBranchAddress("ppr",&ppr); t->SetBranchAddress("etapr",&etapr); t->SetBranchAddress("ppi",&ppi); t->SetBranchAddress("etapi",&etapi); double scalepr = 4; double scalepi = 8; { // reference indicator const double x1pr = 0; const double y1pr = -3; const double x2pr = scalepr; const double y2pr = y1pr; const double versatz = 14; const double x1pi = x1pr+versatz; const double y1pi = -3; const double x2pi = x2pr+versatz+scalepi; const double y2pi = y1pi; TArrow *a; a = new TArrow(x1pr,y1pr,x2pr,y2pr,.01,">"); a->SetLineColor(24); a->Draw(); TLatex tl; tl.SetTextSize(20); tl.SetTextFont(4); tl.DrawLatex(x1pr,y2pr-1.2,"p(p) / 1 GeV"); a = new TArrow(x1pi,y1pi,x2pi,y2pi,.01,">"); a->SetLineColor(20); a->Draw(); tl.SetTextSize(20); tl.SetTextFont(4); tl.DrawLatex(x1pi,y2pi-1.2,"p(#pi) / 1 GeV"); } for (int i = 0; i!=maxN; i++) { t->GetEntry(i); const double thetapr = 2*TMath::ATan(TMath::Exp(-TMath::Abs(etapr))); const double thetapi = 2*TMath::ATan(TMath::Exp(-TMath::Abs(etapi))); const double x1=TMath::Abs(vzl0); const double y1=vrl0; const double x2pr=x1+scalepr*ppr*TMath::Cos(thetapr); const double y2pr=y1+scalepr*ppr*TMath::Sin(thetapr); const double x2pi=x1+scalepi*ppi*TMath::Cos(thetapi); const double y2pi=y1+scalepi*ppi*TMath::Sin(thetapi); TArrow *a; a = new TArrow(x1,y1,x2pr,y2pr,.01,">"); a->SetLineColor(24); a->Draw(); a = new TArrow(x1,y1,x2pi,y2pi,.01,">"); a->SetLineColor(20); a->Draw(); TMarker *m = new TMarker(x1,y1,7); m->SetMarkerColor(28); m->Draw(); } }
// Working at 2D scan by first copying the original ProbScan int MethodDatasetsProbScan::scan2d() { if ( arg->debug ) cout << "MethodDatasetsProbScan::scan2d() : starting ..." << endl; nScansDone++; sanityChecks(); if ( startPars ) delete startPars; // Define whether the 2d contours in hCL are "1D sigma" (ndof=1) or "2D sigma" (ndof=2). // Titus: Change this to 2, since there is no reason to do wrong hCL contours. int ndof = 2; // Set up storage for fit results of this particular // scan. This is used for the drag start parameters. // We cannot use the curveResults2d member because that // only holds better results. vector<vector<RooSlimFitResult*> > mycurveResults2d; for ( int i=0; i<nPoints2dx; i++ ){ vector<RooSlimFitResult*> tmp; for ( int j=0; j<nPoints2dy; j++ ) tmp.push_back(0); mycurveResults2d.push_back(tmp); } // ///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// // // Titus: Saving is done via saveSolutions2d() in the combination, but maybe we need it inline for now \todo: implement a saving function // // Define outputfile // system("mkdir -p root"); // TString probResName = Form("root/scan1dDatasetsProb_" + this->pdf->getName() + "_%ip" + "_" + scanVar1 + ".root", arg->npoints1d); // TFile* outputFile = new TFile(probResName, "RECREATE"); // // Set up toy root tree // this->probScanTree = new ToyTree(this->pdf, arg); // this->probScanTree->init(); // this->probScanTree->nrun = -999; //\todo: why does this branch even exist in the output tree of the prob scan? // ///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// // store start parameters so we can reset them later startPars = new RooDataSet("startPars", "startPars", *w->set(parsName)); startPars->add(*w->set(parsName)); // // start scan from global minimum (not always a good idea as we need to set from other places as well) // setParameters(w, parsName, globalMin); // Define scan parameters and scan range: RooRealVar *par1 = w->var(scanVar1); RooRealVar *par2 = w->var(scanVar2); // Set limit to all parameters. this->loadParameterLimits(); // fix scan parameters par1->setConstant(true); par2->setConstant(true); // Report on the smallest new minimum we come across while scanning. // Sometimes the scan doesn't find the minimum // that was found before. Warn if this happens. double bestMinOld = chi2minGlobal; double bestMinFoundInScan = 100.; // for the status bar int nSteps = 0; float nTotalSteps = nPoints2dx*nPoints2dy; float printFreq = nTotalSteps>100 && !arg->probforce ? 100 : nTotalSteps; ///< number of messages // initialize some control plots gStyle->SetOptTitle(1); TCanvas *cDbg = newNoWarnTCanvas(getUniqueRootName(), Form("DeltaChi2 for 2D scan %i",nScansDone)); cDbg->SetMargin(0.1,0.15,0.1,0.1); float hChi2min2dMin = hChi2min2d->GetMinimum(); bool firstScanDone = hChi2min2dMin<1e5; TH2F *hDbgChi2min2d = histHardCopy(hChi2min2d, firstScanDone); hDbgChi2min2d->SetTitle(Form("#Delta#chi^{2} for scan %i, %s",nScansDone,title.Data())); if ( firstScanDone ) hDbgChi2min2d->GetZaxis()->SetRangeUser(hChi2min2dMin,hChi2min2dMin+25); hDbgChi2min2d->GetXaxis()->SetTitle(par1->GetTitle()); hDbgChi2min2d->GetYaxis()->SetTitle(par2->GetTitle()); hDbgChi2min2d->GetZaxis()->SetTitle("#Delta#chi^{2}"); TH2F *hDbgStart = histHardCopy(hChi2min2d, false); // start coordinates //Titus: start at the global minimum // don't allow the under/overflow bins int iStart = min(hCL2d->GetXaxis()->FindBin(par1->getVal()), hCL2d->GetNbinsX()); int jStart = min(hCL2d->GetYaxis()->FindBin(par2->getVal()), hCL2d->GetNbinsY()); iStart = max(iStart, 1); jStart = max(jStart, 1); hDbgStart->SetBinContent(iStart, jStart, 500.); TMarker *startpointmark = new TMarker(par1->getVal(),par2->getVal(),3); // timer TStopwatch tFit; TStopwatch tSlimResult; TStopwatch tScan; TStopwatch tMemory; // set up the scan spiral int X = 2*nPoints2dx; int Y = 2*nPoints2dy; int x,y,dx,dy; x = y = dx = 0; dy = -1; int t = std::max(X,Y); int maxI = t*t; for ( int spiralstep=0; spiralstep<maxI; spiralstep++ ) { if ((-X/2 <= x) && (x <= X/2) && (-Y/2 <= y) && (y <= Y/2)) { int i = x+iStart; int j = y+jStart; if ( i>0 && i<=nPoints2dx && j>0 && j<=nPoints2dy ) { tScan.Start(false); // status bar if (((int)nSteps % (int)(nTotalSteps/printFreq)) == 0){ cout << Form("MethodDatasetsProbScan::scan2d() : scanning %3.0f%%", (float)nSteps/(float)nTotalSteps*100.) << " \r" << flush; } nSteps++; // status histogram if ( spiralstep>0 ) hDbgStart->SetBinContent(i, j, 500./*firstScan ? 1. : hChi2min2dMin+36*/); // set start parameters from inner turn of the spiral int xStartPars, yStartPars; computeInnerTurnCoords(iStart, jStart, i, j, xStartPars, yStartPars, 1); RooSlimFitResult *rStartPars = mycurveResults2d[xStartPars-1][yStartPars-1]; if ( rStartPars ) setParameters(w, parsName, rStartPars); // memory management: tMemory.Start(false); // delete old, inner fit results, that we don't need for start parameters anymore // for this we take the second-inner-most turn. int iOld, jOld; bool innerTurnExists = computeInnerTurnCoords(iStart, jStart, i, j, iOld, jOld, 2); if ( innerTurnExists ){ deleteIfNotInCurveResults2d(mycurveResults2d[iOld-1][jOld-1]); mycurveResults2d[iOld-1][jOld-1] = 0; } tMemory.Stop(); // alternative choice for start parameters: always from what we found at function call // setParameters(w, parsName, startPars->get(0)); // set scan point float scanvalue1 = hCL2d->GetXaxis()->GetBinCenter(i); float scanvalue2 = hCL2d->GetYaxis()->GetBinCenter(j); par1->setVal(scanvalue1); par2->setVal(scanvalue2); // fit! tFit.Start(false); RooFitResult *fr; // if ( !arg->probforce ) fr = fitToMinBringBackAngles(w->pdf(pdfName), false, -1); // else fr = fitToMinForce(w, combiner->getPdfName()); fr = this->loadAndFit(this->pdf); //Titus: change fitting strategy to the one from the datasets \todo: should be possible to use the fittominforce etc methods // double chi2minScan = 2 * fr->minNll(); //Titus: take 2*minNll vs. minNll? Where is the squared in the main gammacombo? double chi2minScan = 2 * pdf->getMinNll(); tFit.Stop(); tSlimResult.Start(false); RooSlimFitResult *r = new RooSlimFitResult(fr); // try to save memory by using the slim fit result tSlimResult.Stop(); delete fr; allResults.push_back(r); bestMinFoundInScan = TMath::Min((double)chi2minScan, (double)bestMinFoundInScan); mycurveResults2d[i-1][j-1] = r; // If we find a new global minumum, this means that all // previous 1-CL values are too high. We'll save the new possible solution, adjust the global // minimum, return a status code, and stop. // if ( chi2minScan > -500 && chi2minScan<chi2minGlobal ){ //Titus: the hard coded minimum chi2 to avoid ridiculous minima (e.g. at boundaries) only sensible when using the Utils::fitToMin, since the chi2 of the best fit with that fitting method is nominally 0. if ( chi2minScan<chi2minGlobal ){ // warn only if there was a significant improvement if ( arg->debug || chi2minScan<chi2minGlobal-1e-2 ){ if ( arg->verbose ) cout << "MethodDatasetsProbScan::scan2d() : WARNING : '" << title << "' new global minimum found! chi2minGlobal=" << chi2minGlobal << " chi2minScan=" << chi2minScan << endl; } chi2minGlobal = chi2minScan; // recompute previous 1-CL values for ( int k=1; k<=hCL2d->GetNbinsX(); k++ ) for ( int l=1; l<=hCL2d->GetNbinsY(); l++ ){ hCL2d->SetBinContent(k, l, TMath::Prob(hChi2min2d->GetBinContent(k,l)-chi2minGlobal, ndof)); hCLs2d->SetBinContent(k, l, TMath::Prob(hChi2min2d->GetBinContent(k,l)-chi2minBkg, ndof)); } } double deltaChi2 = chi2minScan - chi2minGlobal; double oneMinusCL = TMath::Prob(deltaChi2, ndof); // if ( arg->debug ) { // cout << "chi2minScan: " << chi2minScan << endl; // cout << "chi2minGlobal: " << chi2minGlobal << endl; // cout << "deltaChi2: " << deltaChi2 << endl; // cout << "ndof: " << ndof << endl; // cout << "oneMinusCL: " << oneMinusCL << endl << endl; // } // Save the 1-CL value. But only if better than before! if ( hCL2d->GetBinContent(i, j) < oneMinusCL ){ hCL2d->SetBinContent(i, j, oneMinusCL); double cls_pval = chi2minScan > chi2minBkg ? chi2minScan - chi2minBkg : 0.; hCLs2d->SetBinContent(i, j, TMath::Prob(cls_pval, ndof)); hChi2min2d->SetBinContent(i, j, chi2minScan); hDbgChi2min2d->SetBinContent(i, j, chi2minScan); curveResults2d[i-1][j-1] = r; } // draw/update histograms - doing only every 10th update saves // a lot of time for small combinations if ( ( arg->interactive && ((int)nSteps % 10 == 0) ) || nSteps==nTotalSteps ){ hDbgChi2min2d->Draw("colz"); hDbgStart->Draw("boxsame"); startpointmark->Draw(); cDbg->Update(); } tScan.Stop(); } } // spiral stuff: if( (x == y) || ((x < 0) && (x == -y)) || ((x > 0) && (x == 1-y))) { t = dx; dx = -dy; dy = t; } x += dx; y += dy; } cout << "MethodDatasetsProbScan::scan2d() : scan done. " << endl; if ( arg->debug ){ cout << "MethodDatasetsProbScan::scan2d() : full scan time: "; tScan.Print(); cout << "MethodDatasetsProbScan::scan2d() : - fitting: "; tFit.Print(); cout << "MethodDatasetsProbScan::scan2d() : - create RooSlimFitResults: "; tSlimResult.Print(); cout << "MethodDatasetsProbScan::scan2d() : - memory management: "; tMemory.Print(); } setParameters(w, parsName, startPars->get(0)); saveSolutions2d(); if ( arg->debug ) printLocalMinima(); // confirmSolutions(); //Titus: Leave this out for now, since using it requires compatibility to Utils:fitToMin(), which is not achieved yet // clean all fit results that didn't make it into the final result for ( int i=0; i<allResults.size(); i++ ){ deleteIfNotInCurveResults2d(allResults[i]); } if ( bestMinFoundInScan-bestMinOld > 0.1 ) { cout << "MethodDatasetsProbScan::scan2d() : WARNING: Scan didn't find minimum that was found before!" << endl; cout << "MethodDatasetsProbScan::scan2d() : Are you using too strict parameter limits?" << endl; cout << "MethodDatasetsProbScan::scan2d() : min chi2 found in scan: " << bestMinFoundInScan << ", old min chi2: " << bestMinOld << endl; return 1; } return 0; }
void plotFit(TString filename = "fp-d", TString pltmd = "tph") { // CHECK FOR RIGHT INPUT //////////////////////////////////////////////// string strpltmd = pltmd; if( strpltmd.compare("tph") != 0 && strpltmd.compare("s2b") != 0 && strpltmd.compare("mmp") != 0 ) {error(4);}; // OPEN THE ROOT FILE ////////////////////////////////////////////////// gROOT->Reset(); gROOT->SetStyle("Plain"); gStyle->SetTitleBorderSize(0); gStyle->SetPalette(1); TCanvas *MyC = new TCanvas("MyC","Plot of the GAPP fit",200,10,700,500); // Still to do: Automate the frame boundaries. string strfile = filename, rootname = strfile + ".root"; TFile *rootfile = TFile::Open(rootname.c_str()); if(rootfile == NULL) error(1); TTree *tree = (TTree*)rootfile->Get(strfile.c_str()); if(tree == NULL) error(2); TBranch *fits2bbranch = (TBranch*)tree->GetBranch("fits2b"); TBranch *fittphbranch = (TBranch*)tree->GetBranch("fittph"); TBranch *fitxbranch = (TBranch*)tree->GetBranch("fitx"); if( (fits2bbranch == NULL) || (fittphbranch == NULL) || (fitxbranch == NULL) ) error(3); Float_t fits2b, fittph, fitx; tree->SetBranchAddress("fits2b",&fits2b); tree->SetBranchAddress("fittph",&fittph); tree->SetBranchAddress("fitx", &fitx); // GET GRID ///////////////////////////////////////////////////////////// Int_t Npoints = (Int_t)tree->GetEntries(); Int_t point, zSteps = 0, ySteps = 0; Float_t fitxMin = 100000, fitxMax = -1.0, s2bMin = 100000, tphMin = 100000; Float_t s2bMax = -1.0, tphMax = -1.0; for(point=0; point<Npoints; point++) { tree->GetEntry(point); if( fits2b > s2bMax ) {zSteps++; s2bMax = fits2b;} if( fittph > tphMax ) {ySteps++; tphMax = fittph;} }; const int s2bSteps = zSteps, tphSteps = ySteps, mmpSteps = Npoints; Float_t s2bValues[s2bSteps], tphValues[tphSteps]; s2bMax = -1.0, tphMax = -1.0; int s2bStep = 0, tphStep = 0; for(point=0; point<Npoints; point++) { tree->GetEntry(point); if( fits2b > s2bMax ) {s2bValues[s2bStep] = fits2b; s2bStep++; s2bMax = fits2b;} if( fittph > tphMax ) {tphValues[tphStep] = fittph; tphStep++; tphMax = fittph;} }; // PREPARE PLOT ///////////////////////////////////////////////////////// string plottitle = "Model: " + strfile + " | Plot: "; if( strpltmd.compare("tph") == 0 ) { plottitle += "tan^{2}(#phi) over x for fixed sin^{2}(2#beta)"; } if( strpltmd.compare("s2b") == 0 ) { plottitle += "sin^{2}(2#beta) over x for fixed tan^{2}(#phi)"; } if( strpltmd.compare("mmp") == 0 ) { plottitle += "Masses of the new heavy gauge bosons"; } // PLOT DATA /////////////////////////////////////////////////////////// if (strpltmd.compare("tph") == 0) { tphMax = -1.0; // TGraph *tphplots[s2bSteps]; TGraph *tphplots[s2bSteps-30]; // for(s2bStep=0; s2bStep<s2bSteps; s2bStep++) { for(s2bStep=0; s2bStep<s2bSteps-30; s2bStep++) { Float_t tphArray[tphSteps], fitxArray[tphSteps]; tphStep = 0; for(point=0; point<Npoints; point++) { tree->GetEntry(point); if(fits2b == s2bValues[s2bStep]) { // tphArray[tphStep] = fittph; tphArray[tphStep] = sqrt(1.0/(1.0+fittph)); if (fittph < tphMin) tphMin = fittph; if (fittph > tphMax) tphMax = fittph; fitxArray[tphStep] = fitx; if (fitx < fitxMin) fitxMin = fitx; if (fitx > fitxMax) fitxMax = fitx; tphStep++; TMarker *m = new TMarker(fitxArray[tphStep],tphArray[tphStep],20); m->SetMarkerSize(2); m->SetMarkerColor(31+tphStep); m->Draw(); } } if (s2bStep == 0) { TH1F* frame = MyC->DrawFrame(0.0,0.0,1.1*fitxMax,1.1); // TH1F* frame = MyC->DrawFrame(0.7*fitxMin,0.7*tphMin,1.1*fitxMax,1.1*tphMax); TAxis *xaxis = frame->GetXaxis(); TAxis *yaxis = frame->GetYaxis(); xaxis->SetTitle("x = u^{2}/v^{2}"); xaxis->CenterTitle(); xaxis->SetTitleOffset(1.); xaxis->SetDecimals(); xaxis->SetLabelSize(0.03); xaxis->SetLabelOffset(0.01); yaxis->SetTitle("tan^{2}(#phi)"); yaxis->CenterTitle(); yaxis->SetTitleOffset(1.); yaxis->SetDecimals(); yaxis->SetLabelSize(0.03); yaxis->SetLabelOffset(0.01); frame->SetTitle(plottitle.c_str()); } tphplots[s2bStep] = new TGraph(tphSteps,fitxArray,tphArray); tphplots[s2bStep]->SetMarkerStyle(20); tphplots[s2bStep]->SetMarkerSize(0.4); tphplots[s2bStep]->Draw("CP"); } } else if (strpltmd.compare("s2b") == 0) { s2bMax = -1.0; TGraph *s2bplots[tphSteps-100]; for(tphStep=0; tphStep<tphSteps-100; tphStep++) { Float_t s2bArray[s2bSteps], fitxArray[s2bSteps]; s2bStep = 0; for(point=0; point<Npoints; point++) { tree->GetEntry(point); if(fittph == tphValues[tphStep+20]) { s2bArray[s2bStep] = fits2b; if (fits2b < s2bMin) s2bMin = fits2b; if (fits2b > s2bMax) s2bMax = fits2b; fitxArray[s2bStep] = fitx; if (fitx < fitxMin) fitxMin = fitx; if (fitx > fitxMax) fitxMax = fitx; s2bStep++; } } if (tphStep == 0) { TH1F* frame = MyC->DrawFrame(0.8*fitxMin,0.95*s2bMin,1.2*fitxMax,1.05*s2bMax); TAxis *xaxis = frame->GetXaxis(); TAxis *yaxis = frame->GetYaxis(); xaxis->SetTitle("x = u^{2}/v^{2}"); xaxis->CenterTitle(); xaxis->SetTitleOffset(1.); xaxis->SetDecimals(); xaxis->SetLabelSize(0.03); xaxis->SetLabelOffset(0.01); yaxis->SetTitle("sin^{2}(2#beta)"); yaxis->CenterTitle(); yaxis->SetTitleOffset(1.25); yaxis->SetDecimals(); yaxis->SetLabelSize(0.03); yaxis->SetLabelOffset(0.01); frame->SetTitle(plottitle.c_str()); } s2bplots[tphStep] = new TGraph(s2bSteps,fitxArray,s2bArray); s2bplots[tphStep]->SetMarkerStyle(20); s2bplots[tphStep]->SetMarkerSize(0.4); s2bplots[tphStep]->Draw("C"); } } else if (strpltmd.compare("mmp") == 0) { Float_t mzpArray[mmpSteps], mwpArray[mmpSteps]; Float_t mzpMin = 100000, mzpMax = -1.0, mwpMin = 100000, mwpMax = -1.0; Float_t fitsph, fitcph; Float_t Cz1, Cz2, Cz3, Cw1, Cw2, Cw3, Cw4, C1, C2; string mdl(strfile,0,2); if ( (mdl.compare("lr") == 0) || (mdl.compare("lp") == 0) || (mdl.compare("hp") == 0) || (mdl.compare("fp") == 0) ) { string Higgs(strfile,3,1); if (Higgs.compare("d") == 0) { Cz1 = 11.95349795785275; Cz2 = 30.63269990028513; Cz3 = 42.58619785813789; Cw1 = 21.29309892906894; Cw2 = 9.339600971216193; Cw3 = 30.63269990028513; Cw4 = 42.58619785813789; } else if (Higgs.compare("t") == 0) { Cz1 = 5.976748978926375; Cz2 = 30.63269990028513; Cz3 = 85.17239571627579; Cw1 = 15.05649464522066; Cw2 = 3.302047590161717; Cw3 = 21.66058982554409; Cw4 = 60.22597858088265; } else {error(6);} for(point=0; point<Npoints; point++) { tree->GetEntry(point); fitsph = fittph / (1.0 + fittph); fitcph = 1.0 - fitsph; if (fitsph != 0.0) { mzpArray[point] = (0.001/sqrt(fitsph*fitcph*fitx)) * (Cz1*fitcph*fitcph + Cz2*fits2b + Cz3*fitx); if (mzpArray[point] < mzpMin) mzpMin = mzpArray[point]; if (mzpArray[point] > mzpMax) mzpMax = mzpArray[point]; mwpArray[point] = (0.001/sqrt(fitsph*fitx)) * (Cw1 - Cw2*fitcph*fitcph + Cw3*fits2b + Cw4*fitx); if (mwpArray[point] < mwpMin) mwpMin = mwpArray[point]; if (mwpArray[point] > mwpMax) mwpMax = mwpArray[point]; } else { mzpArray[point] = 0.0; mwpArray[point] = 0.0; } } } else if ( (mdl.compare("uu") == 0) || (mdl.compare("nu") == 0) ) { C1 = 94.0397928463607 C2 = 77.1253849720165 for(point=0; point<Npoints; point++) { tree->GetEntry(point); fitsph = fittph / (1.0+fittph); fitcph = 1.0 - fitsph; if (fitsph != 0.0) { mzpArray[point] = (0.001/sqrt(fitsph*fitcph*fitx)) * (C1*fitsph*fitsph + C2*fitx); if (mzpArray[point] < mzpMin) mzpMin = mzpArray[point]; if (mzpArray[point] > mzpMax) mzpMax = mzpArray[point]; mwpArray[point] = (0.001/sqrt(fitsph*fitcph*fitx)) * (C1*fitsph*fitsph + C2*fitx); if (mwpArray[point] < mwpMin) mwpMin = mwpArray[point]; if (mwpArray[point] > mwpMax) mwpMax = mwpArray[point]; } else { mzpArray[point] = 0.0; mwpArray[point] = 0.0; } } }
void q4gep_jan2009() { gROOT->SetStyle("HALLA"); TCanvas *cn = new TCanvas("cn"); cn->Draw(); TH1F *frm = new TH1F("frm","",100,0.,15.); #ifndef DIFF_LAB frm->GetXaxis()->SetTitle("Q^{2} [GeV^{2}]"); #endif #ifdef DIFF_LAB frm->GetXaxis()->SetTitle("Q^{2} in GeV^{2}"); #endif frm->GetXaxis()->CenterTitle(); #ifndef DIFF_LAB frm->GetYaxis()->SetTitle("Q^{4} G_{E}^{p}"); #endif #ifdef DIFF_LAB frm->GetYaxis()->SetTitle("#muG_{Ep}/G_{Mp}"); #endif frm->GetYaxis()->CenterTitle(); frm->SetMinimum(-0.1); frm->SetMaximum(0.4); frm->UseCurrentStyle(); frm->Draw(); TF1 *ff = new TF1("ff",fff,1.,10.,1.); ff->SetLineColor(2); ff->SetLineStyle(2); ff->SetParameter(0,1.); ff->SetParameter(1,.3); // match to Gayou point at 5.54 GeV^2 ff->SetParameter(0,.1035); TMultiGraph* mgrDta = new TMultiGraph("Data","#frac{#mu_{p}G_{E}^{p}}{G_{M}^{p}}"); //TLegend *legDta = new TLegend(.58,.565,.875,.9,"","brNDC"); TLegend *legDta = new TLegend(.648,.682,.945,.900,"","brNDC"); TMultiGraph* wgr = mgrDta; TLegend *wlg = legDta; // the data legDta->SetBorderSize(0); // turn off border legDta->SetFillStyle(0); datafile_t *f = datafiles; TGraph* gr=0; TGraph* ogr=0; while ( f && f->filename ) { ogr=OneGraph(f); if (ogr) { gr = fromGEpGMptoQ4GEp(ogr); gr->SetLineWidth(2); if (f->lnpt) { mgrDta->Add(gr,f->lnpt); legDta->AddEntry(gr,f->label,f->lnpt); } else if (gr->GetMarkerStyle()>=20) { mgrDta->Add(gr,"p"); legDta->AddEntry(gr,f->label,"p"); } else { mgrDta->Add(gr,"l"); legDta->AddEntry(gr,f->label,"l"); } } f++; } mgrDta->Draw("p"); legDta->Draw();// don't draw the data legend TMultiGraph* mgrThry = new TMultiGraph("Theory","#frac{#mu_{p}G_{E}^{p}}{G_{M}^{p}}"); // TLegend *legThry = new TLegend(.585,.6,.9,.9,"","brNDC"); TLegend *legThry = new TLegend(.335,.689,.649,.911,"","brNDC"); wgr = mgrThry; wlg = legThry; // the theory wlg->SetBorderSize(0); // turn off border wlg->SetFillStyle(0); f = theoryfiles1; ogr=0; while ( f && f->filename ) { ogr=OneGraph(f); if (ogr) { gr = fromGEpGMptoQ4GEp(ogr); TGraphAsymmErrors *egr = dynamic_cast<TGraphAsymmErrors*>(gr); if (egr && egr->GetEYhigh() && egr->GetEYhigh()[0]>0) { gr = toerror_band(egr); gr->SetFillStyle(3000+f->style); } gr->SetLineWidth(2); if (f->lnpt) { wgr->Add(gr,f->lnpt); if( f->label[0] != 'x' ){ wlg->AddEntry(gr,f->label,f->lnpt); } } else if (gr->GetMarkerStyle()>=20) { wgr->Add(gr,"p"); if( f->label[0] != 'x' ){ wlg->AddEntry(gr,f->label,"p"); } } else { wgr->Add(gr,"l"); if( f->label[0] != 'x' ){ wlg->AddEntry(gr,f->label,"l"); } } } f++; } mgrThry->Draw("c"); #ifdef PQCD ff->Draw("same"); legThry->AddEntry(ff,"F2/F1 #propto ln^{2}(Q^{2}/#Lambda^{2})/Q^{2}","l"); #endif// PQCD legThry->Draw(); cn->Modified(); cn->Update(); TFrame* frame = gPad->GetFrame(); // draw a line at 1 TLine *line1 = new TLine(frame->GetX1(),1.,frame->GetX2(),1.); line1->SetLineStyle(1); line1->Draw(); return; cn->Update(); cn->SaveAs(Form("%s.eps",psfile)); cn->SaveAs(Form("%s.root",psfile)); gSystem->Exec(Form("./replace_symbols.pl %s.eps",psfile)); // now an overlay, hopefully matching dimensions // remove everything except the frame cn->Update(); TList *clist = cn->GetListOfPrimitives(); frame = cn->GetFrame(); for (int i=0; i<clist->GetSize(); ) { if (clist->At(i) != frame) { clist->RemoveAt(i); } else i++; } // draw markers in the corners TMarker *mkr = new TMarker(frame->GetX1(),frame->GetY1(),2); mkr->Draw(); mkr = new TMarker(frame->GetX2(),frame->GetY1(),2); mkr->Draw(); mkr = new TMarker(frame->GetX1(),frame->GetY2(),2); mkr->Draw(); mkr = new TMarker(frame->GetX2(),frame->GetY2(),2); mkr->Draw(); frame->SetLineColor(10); { datafile_t goeckeler = { "figure_input/Goeckeler/goeckeler.dat","Goeckeler", "[0]","[1]","[1]-[2]","[3]-[1]",0,0,1,4,"F" }; gr = OneGraph(&goeckeler); TGraphAsymmErrors *egr = dynamic_cast<TGraphAsymmErrors*>(gr); if (egr && egr->GetEYhigh() && egr->GetEYhigh()[egr->GetN()/2]>0) { gr = toerror_band(egr); gr->SetLineStyle(1); gr->SetFillColor(gr->GetLineColor()); gr->SetFillStyle(3000+goeckeler.style); } gr->Draw("F"); cn->Update(); cn->SaveAs("mup_gep_gmp_goeckeler_Overlay.eps"); cn->SaveAs("mup_gep_gmp_goeckeler_Overlay.root"); // remove Goeckeler curve from plot clist->Remove(gr); } { datafile_t miller = { "figure_input/Miller/lattice.GEp_GMp.rtf","Miller", "[0]","[1]","[1]-[3]","[2]-[1]",0,0,1,3,"F" }; gr = OneGraph(&miller); TGraphAsymmErrors* egr = dynamic_cast<TGraphAsymmErrors*>(gr); if (egr && egr->GetEYhigh() && egr->GetEYhigh()[egr->GetN()/2]>0) { gr = toerror_band(egr); gr->SetLineStyle(1); gr->SetFillColor(gr->GetLineColor()); gr->SetFillStyle(3000+miller.style); } gr->Draw("F"); cn->Update(); cn->SaveAs("mup_gep_gmp_Miller_Overlay.eps"); cn->SaveAs("mup_gep_gmp_Miller_Overlay.root"); clist->Remove(gr); } { datafile_t zanotti = { "figure_input/Zanotti/zanotti_gepgmp.dat","Zanotti", "[0]","[1]","[1]-[2]","[3]-[1]",0,0,1,6,"F" }; gr = OneGraph(&zanotti); TGraphAsymmErrors* egr = dynamic_cast<TGraphAsymmErrors*>(gr); if (egr && egr->GetEYhigh() && egr->GetEYhigh()[egr->GetN()/2]>0) { gr = toerror_band(egr); gr->SetLineStyle(1); gr->SetFillColor(gr->GetLineColor()); gr->SetFillStyle(3000+zanotti.style); } gr->Draw("F"); cn->Update(); cn->SaveAs("mup_gep_gmp_Zanotti_Overlay.eps"); cn->SaveAs("mup_gep_gmp_Zanotti_Overlay.root"); clist->Remove(gr); } }