/* ====================================
* Plot Magnetic Permeability vs Field
====================================*/
void makePlot_uvB()
{
const bool ideal_cloak = true;
const bool plot_FM_Fm602 = true;
const bool plot_FM_Fm625 = false;
const bool plot_FM_Fm612 = true;
const bool plot_FM_Fm651 = true;
TCanvas *c_uvB = new TCanvas();
TH1 *h_uvB = c_uvB->DrawFrame(0, 1.0, 60, 4.5);
h_uvB->SetTitle("Permeability of Epoxy/Steel Ferromagnets");
h_uvB->GetXaxis()->SetTitle("B_{0} (mT)");
h_uvB->GetXaxis()->SetTitleOffset(1.5);
h_uvB->GetYaxis()->SetTitle("#mu_{r}");
h_uvB->GetYaxis()->SetTitleOffset(1.0);
/*Create Legend*/
TLegend *leg_uvB = new TLegend(0.55,0.72,0.8,0.87);
leg_uvB->SetNColumns(1);
leg_uvB->Draw();
vector<double> u, Fm;
// u.push_back(1.0);
// Fm.push_back(0.0);
// u.push_back(550.);
// Fm.push_back(1.0);
if(plot_FM_Fm651)
{
/*
* Permeability Plot for New FM = 0.602 Ferromagnet
*/
/*Run calibration*/
const TString calib_fm651 = "../Data/Calib_Data/DataFile_160916_211714.txt";
const TString fmscan_fm651 = "../Data/FMScan_Data/DataFile_160916_212729_Part2.txt";
TF1 *calib_fit1 = Calibrate(calib_fm651);
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "../Data/Calib_Data/fm503_di.txt";
const TString do_file = "../Data/Calib_Data/fm503_do.txt";
double R_fm651 = ratio(di_file, do_file);
double R_sig_fm651 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm651 = plot_uvB(fmscan_fm651, calib_fit1, R_fm651, R_sig_fm651);
g_fm651->Draw("LP");
g_fm651->SetLineColor(kViolet);
g_fm651->SetMarkerColor(kViolet);
TF1 *fit_fm651 = new TF1("fit_fm651", "[0]/x + [1]", 0, 60);
fit_fm651->SetLineStyle(2);
fit_fm651->SetLineColor(kViolet);
g_fm651->Fit("fit_fm651");
leg_uvB->AddEntry( g_fm651 , "F_{m} = 0.651" , "lp");
u.push_back( fit_fm651->Eval(50.) );
Fm.push_back( 0.651 );
}
if(plot_FM_Fm602)
{
/*
* Permeability Plot for New FM = 0.602 Ferromagnet
*/
/*Run calibration*/
const TString calib_fm602 = "../Data/Calib_Data/DataFile_160727_162618.txt";
const TString fmscan_fm602 = "../Data/FMScan_Data/DataFile_160727_163533.txt";
TF1 *calib_fit1 = Calibrate(calib_fm602);
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "../Data/Calib_Data/fm612_di.txt";
const TString do_file = "../Data/Calib_Data/fm612_do.txt";
double R_fm602 = ratio(di_file, do_file);
double R_sig_fm602 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm602 = plot_uvB(fmscan_fm602, calib_fit1, R_fm602, R_sig_fm602);
g_fm602->Draw("LP");
g_fm602->SetLineColor(kRed);
g_fm602->SetMarkerColor(kRed);
TF1 *fit_fm602 = new TF1("fit_fm602", "[0]/x + [1]", 0, 60);
fit_fm602->SetLineStyle(2);
fit_fm602->SetLineColor(kRed);
g_fm602->Fit("fit_fm602");
leg_uvB->AddEntry( g_fm602 , "F_{m} = 0.602" , "lp");
u.push_back( fit_fm602->Eval(50.) );
Fm.push_back( 0.602 );
}
if(plot_FM_Fm612)
{
/*
* Permeability Plot for New FM = 0.612 Ferromagnet
*/
/*Run calibration*/
const TString calib_fm612 = "../Data/Calib_Data/DataFile_160805_142716.txt";
const TString fmscan_fm612 = "../Data/FMScan_Data/DataFile_160805_143508.txt";
TF1 *calib_fit1 = Calibrate(calib_fm612);
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "../Data/Calib_Data/fm612_di.txt";
const TString do_file = "../Data/Calib_Data/fm612_do.txt";
double R_fm612 = ratio(di_file, do_file);
double R_sig_fm612 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm612 = plot_uvB(fmscan_fm612, calib_fit1, R_fm612, R_sig_fm612);
g_fm612->Draw("LP");
g_fm612->SetLineColor(kBlue+2);
g_fm612->SetMarkerColor(kBlue+2);
TF1 *fit_fm612 = new TF1("fit_fm612", "[0]/x + [1]", 0, 60);
fit_fm612->SetLineStyle(2);
fit_fm612->SetLineColor(kBlue+2);
g_fm612->Fit("fit_fm612");
// cout << fit_fm612->Eval(500) << endl;
leg_uvB->AddEntry( g_fm612 , "F_{m} = 0.612" , "lp");
u.push_back( fit_fm612->Eval(50.) );
Fm.push_back( 0.612 );
}
if(plot_FM_Fm625)
{
/*
* Permeability Plot for New FM = 0.625 Ferromagnet
*/
/*Run calibration*/
const TString calib_fm625 = "../Data/Calib_Data/DataFile_160801_121348.txt";
const TString fmscan_fm625 = "../Data/FMScan_Data/DataFile_160801_122800.txt";
TF1 *calib_fit1 = Calibrate(calib_fm625);
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "../Data/Calib_Data/fm612_di.txt";
const TString do_file = "../Data/Calib_Data/fm612_do.txt";
double R_fm625 = ratio(di_file, do_file);
double R_sig_fm625 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm625 = plot_uvB(fmscan_fm625, calib_fit1, R_fm625, R_sig_fm625);
g_fm625->Draw("LP");
g_fm625->SetLineColor(kGreen+2);
g_fm625->SetMarkerColor(kGreen+2);
// TF1 *fit_fm625 = new TF1("fit_fm625", "[0]/x + [1]", 0, 60);
// fit_fm625->SetLineStyle(2);
// fit_fm625->SetLineColor(kBlue+2);
// g_fm625->Fit("fit_fm625");
leg_uvB->AddEntry( g_fm625 , "F_{m} = 0.625" , "lp");
}
if(ideal_cloak)
{
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "../Data/Calib_Data/fm503_di.txt";
const TString do_file = "../Data/Calib_Data/fm503_do.txt";
double R_fv30 = ratio(di_file, do_file);
// cout << R_fv30 << endl;
double R_sig_fv30 = r_sig(di_file, do_file);
// cout << R_sig_fv30 << endl;
/*Calculate theoretical permeability of ferromagnet*/
double u_cloak = (1 + pow(R_fv30,2)) / (1 - pow(R_fv30,2));
cout << "Desired Permeability: " << u_cloak << endl;
/*Draw Line for Permeability of Ideal Cloak*/
TLine *l_ucloak = new TLine(0.0, u_cloak, 60.0, u_cloak);
l_ucloak->SetLineStyle(2);
// l_ucloak->SetLineColor(kRed);
l_ucloak->Draw();
leg_uvB->AddEntry( l_ucloak , "Ideal Cloak" , "l");
}
/*
* Plot of Permeability vs Fractional Mass of Steel Powder in Ferromagnet
*/
TGraph *g_uvFm = new TGraph(Fm.size(), &(Fm[0]), &(u[0]) );
TCanvas *c_uvFm = new TCanvas();
g_uvFm->Draw("AP");
// g_uvFm->Fit("expo");
g_uvFm->SetTitle("#mu_{r} vs F_{m}");
g_uvFm->GetXaxis()->SetTitle("F_{m}");
g_uvFm->GetYaxis()->SetTitle("#mu_{r}");
/*Save plots to png and eps files*/
/*
c_uvB->Print("../../Plots/uvB/EpoxySteel_uvB_160916.png");
c_uvB->Print("../../Plots/uvB/EpoxySteel_uvB_160916.eps");
*/
return;
}
int makePlot_uvB_MRI()
{
TCanvas *c_uvB = new TCanvas();
TH1 *h_uvB = c_uvB->DrawFrame(0, 1.0, 550, 5.0);
// h_uvB->SetTitle("Permeability of Epoxy/Steel Ferromagnets");
h_uvB->SetTitle("");
h_uvB->GetXaxis()->SetTitle("B_{0} (mT)");
h_uvB->GetXaxis()->SetTitleOffset(1.5);
h_uvB->GetYaxis()->SetTitle("#mu_{r}");
h_uvB->GetYaxis()->SetTitleOffset(1.0);
/*Create Legend*/
TLegend *leg_uvB = new TLegend(0.45,0.65,0.83,0.85);
leg_uvB->SetNColumns(1);
leg_uvB->Draw();
const bool plot_554ramp = true;
const bool plot_699ramp2 = true;
const bool plot_699ramp3 = true;
const bool plot_699ramp4 = true;
const bool plot_699ramp5 = true;
const bool plot_618ramp = true;
const bool plot_554at200 = true;
const bool plot_618at200 = true;
const bool plot_673at200 = true;
const bool plot_699at200 = true;
const bool plot_554at300 = true;
const bool plot_618at300 = true;
const bool plot_673at300 = true;
const bool plot_699at300 = true;
const bool plot_554at400 = true;
const bool plot_618at400 = true;
const bool plot_673at400 = true;
const bool plot_699at400 = true;
const bool plot_554at500 = true;
const bool plot_618at500 = true;
const bool plot_673at500 = true;
const bool plot_699at500 = true;
const bool ideal_cloak = true;
const bool exp_cloak = true;
// vector<double> u, Fm, u_err, Fm_err;
// /*Include Point for Fm = 0.0, No Steel Powder*/
// u.push_back(1.0);
// Fm.push_back(0.0);
// u_err.push_back(0.03);
// Fm_err.push_back(0.0);
// /*Include Point for Fm = 1.0, Only Steel Powder*/
// // u.push_back(1000.);
// // Fm.push_back(1.0);
// // u_err.push_back(0.0);
// // Fm_err.push_back(0.0);
if(plot_554at200)
{
/*
* Permeability Plot for Fm = 0.544 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_20-31-45.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm554_di.txt";
const TString do_file = "diameter-files/fm554_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_point(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
// g_fm554->SetLineColor(kRed);
g_fm554->SetMarkerColor(kBlue);
g_fm554->SetMarkerStyle(22);
g_fm554->SetMarkerSize(1.5);
leg_uvB->AddEntry( g_fm554 , "F_{m} = 0.554" , "p");
}
if(plot_554at300)
{
/*
* Permeability Plot for Fm = 0.544 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_20-58-06.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm554_di.txt";
const TString do_file = "diameter-files/fm554_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_point(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kBlue);
g_fm554->SetMarkerStyle(22);
g_fm554->SetMarkerSize(1.5);
}
if(plot_554at400)
{
/*
* Permeability Plot for Fm = 0.544 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_21-18-15.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm554_di.txt";
const TString do_file = "diameter-files/fm554_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_point(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kBlue);
g_fm554->SetMarkerStyle(22);
g_fm554->SetMarkerSize(1.5);
}
if(plot_554at500)
{
/*
* Permeability Plot for Fm = 0.544 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_21-37-52.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm554_di.txt";
const TString do_file = "diameter-files/fm554_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_point(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kBlue);
g_fm554->SetMarkerStyle(22);
g_fm554->SetMarkerSize(1.5);
}
if(plot_618at200)
{
/*
* Permeability Plot for Fm = 0.618 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_20-35-17.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm618_di.txt";
const TString do_file = "diameter-files/fm618_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_point(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kGreen+2);
g_fm554->SetMarkerStyle(22);
g_fm554->SetMarkerSize(1.5);
leg_uvB->AddEntry( g_fm554 , "F_{m} = 0.618" , "p");
}
if(plot_618at300)
{
/*
* Permeability Plot for Fm = 0.618 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_21-00-57.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm618_di.txt";
const TString do_file = "diameter-files/fm618_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_point(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kGreen+2);
g_fm554->SetMarkerStyle(22);
g_fm554->SetMarkerSize(1.5);
}
if(plot_618at400)
{
/*
* Permeability Plot for Fm = 0.618 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_21-20-51.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm618_di.txt";
const TString do_file = "diameter-files/fm618_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_point(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kGreen+2);
g_fm554->SetMarkerStyle(22);
g_fm554->SetMarkerSize(1.5);
}
if(plot_618at500)
{
/*
* Permeability Plot for Fm = 0.618 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_21-40-20.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm618_di.txt";
const TString do_file = "diameter-files/fm618_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_point(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kGreen+2);
g_fm554->SetMarkerStyle(22);
g_fm554->SetMarkerSize(1.5);
}
if(plot_673at200)
{
/*
* Permeability Plot for Fm = 0.673 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_20-37-24.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm673_di.txt";
const TString do_file = "diameter-files/fm673_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_point(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kMagenta);
g_fm554->SetMarkerStyle(22);
g_fm554->SetMarkerSize(1.5);
leg_uvB->AddEntry( g_fm554 , "F_{m} = 0.673" , "p");
}
if(plot_673at300)
{
/*
* Permeability Plot for Fm = 0.618 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_21-03-21.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm673_di.txt";
const TString do_file = "diameter-files/fm673_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_point(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kMagenta);
g_fm554->SetMarkerStyle(22);
g_fm554->SetMarkerSize(1.5);
}
if(plot_673at400)
{
/*
* Permeability Plot for Fm = 0.618 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_21-22-51.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm673_di.txt";
const TString do_file = "diameter-files/fm673_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_point(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kMagenta);
g_fm554->SetMarkerStyle(22);
g_fm554->SetMarkerSize(1.5);
}
if(plot_673at500)
{
/*
* Permeability Plot for Fm = 0.618 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_21-42-13.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm673_di.txt";
const TString do_file = "diameter-files/fm673_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_point(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kMagenta);
g_fm554->SetMarkerStyle(22);
g_fm554->SetMarkerSize(1.5);
}
if(plot_699at200)
{
/*
* Permeability Plot for Fm = 0.699 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_20-29-24.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm699_di.txt";
const TString do_file = "diameter-files/fm699_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_point(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kRed);
g_fm554->SetMarkerStyle(22);
g_fm554->SetMarkerSize(1.5);
leg_uvB->AddEntry( g_fm554 , "F_{m} = 0.699" , "p");
}
if(plot_699at300)
{
/*
* Permeability Plot for Fm = 0.699 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_20-54-13.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm699_di.txt";
const TString do_file = "diameter-files/fm699_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_point(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kRed);
g_fm554->SetMarkerStyle(22);
g_fm554->SetMarkerSize(1.5);
}
if(plot_699at400)
{
/*
* Permeability Plot for Fm = 0.699 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_21-15-59.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm699_di.txt";
const TString do_file = "diameter-files/fm699_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_point(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kRed);
g_fm554->SetMarkerStyle(22);
g_fm554->SetMarkerSize(1.5);
}
if(plot_699at500)
{
/*
* Permeability Plot for Fm = 0.699 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_21-36-02.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm699_di.txt";
const TString do_file = "diameter-files/fm699_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_point(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kRed);
g_fm554->SetMarkerStyle(22);
g_fm554->SetMarkerSize(1.5);
}
if(plot_554ramp)
{
/*
* Permeability Plot for Fm = 0.554 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_19-56-40.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm554_di.txt";
const TString do_file = "diameter-files/fm554_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_ramp(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
// g_fm554->SetLineColor(kRed);
g_fm554->SetMarkerColor(kBlue);
// leg_uvB->AddEntry( g_fm554 , "F_{m} = 0.554 Ramp" , "p");
}
if(plot_699ramp2)
{
/*
* Permeability Plot for Fm = 0.699 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_20-14-15.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm699_di.txt";
const TString do_file = "diameter-files/fm699_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_ramp(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kRed);
// leg_uvB->AddEntry( g_fm554 , "F_{m} = 0.699 Ramp" , "p");
}
if(plot_699ramp3)
{
/*
* Permeability Plot for Fm = 0.699 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_20-41-31.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm699_di.txt";
const TString do_file = "diameter-files/fm699_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_ramp(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kRed);
// leg_uvB->AddEntry( g_fm554 , "F_{m} = 0.699 Ramp" , "p");
}
if(plot_699ramp4)
{
/*
* Permeability Plot for Fm = 0.699 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_21-06-35.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm699_di.txt";
const TString do_file = "diameter-files/fm699_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_ramp(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kRed);
// leg_uvB->AddEntry( g_fm554 , "F_{m} = 0.699 Ramp" , "p");
}
if(plot_699ramp5)
{
/*
* Permeability Plot for Fm = 0.699 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_21-25-26.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm699_di.txt";
const TString do_file = "diameter-files/fm699_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_ramp(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kRed);
// leg_uvB->AddEntry( g_fm554 , "F_{m} = 0.699 Ramp" , "p");
}
if(plot_618ramp)
{
/*
* Permeability Plot for Fm = 0.618 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_21-45-58.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm699_di.txt";
const TString do_file = "diameter-files/fm699_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_ramp(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kGreen+2);
// leg_uvB->AddEntry( g_fm554 , "F_{m} = 0.699 Ramp" , "p");
}
if(plot_673ramp)
{
/*
* Permeability Plot for Fm = 0.673 Ferromagnet
*/
const TString fmscan_fm554 = "data/DATA_MegaVIEW/DataFile_2016-12-08_21-45-58.txt";
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm673_di.txt";
const TString do_file = "diameter-files/fm673_do.txt";
double R_fm554 = ratio(di_file, do_file);
double R_sig_fm554 = r_sig(di_file, do_file);
/*Plot u vs B for FM*/
TGraphErrors *g_fm554 = plot_ramp(fmscan_fm554, R_fm554, R_sig_fm554);
g_fm554->Draw("P");
g_fm554->SetMarkerColor(kGreen+2);
// leg_uvB->AddEntry( g_fm554 , "F_{m} = 0.699 Ramp" , "p");
}
if(ideal_cloak)
{
/*Calculate Radius Ratio and Uncertainty*/
const TString di_file = "diameter-files/fm554_di.txt";
const TString do_file = "diameter-files/fm554_do.txt";
double R_fv30 = ratio(di_file, do_file);
// cout << R_fv30 << endl;
double R_sig_fv30 = r_sig(di_file, do_file);
// cout << R_sig_fv30 << endl;
/*Calculate theoretical permeability of ferromagnet*/
double u_cloak = (1. + pow(R_fv30,2.)) / (1. - pow(R_fv30,2.));
cout << "Desired Permeability: " << u_cloak << endl;
/*Draw Line for Permeability of Ideal Cloak*/
TLine *l_ucloak = new TLine(0.0, u_cloak, 500.0, u_cloak);
l_ucloak->SetLineStyle(2);
l_ucloak->Draw();
leg_uvB->AddEntry( l_ucloak , "Ideal Cloak - Theory" , "l");
TLine *l_ucloak2 = new TLine(0.0, 2.0706, 500.0, 2.0706);
l_ucloak2->SetLineStyle(2);
l_ucloak2->SetLineColor(kBlue);
l_ucloak2->Draw();
leg_uvB->AddEntry( l_ucloak2 , "Ideal Cloak - Experiment" , "l");
}
// /*
// * Plot of Permeability vs Fractional Mass of Steel Powder in Ferromagnet
// */
// TGraphErrors *g_uvFm = new TGraphErrors(Fm.size(), &(Fm[0]), &(u[0]), &(Fm_err[0]), &(u_err[0]) );
// //TGraph *g_uvFm = new TGraph(Fm.size(), &(Fm[0]), &(u[0]) );
// TCanvas *c_uvFm = new TCanvas();
// g_uvFm->Draw("AP");
// TF1 *fit = new TF1("fit", "[0]/TMath::Tan([1]*x + [2]) + [3]", 0.0, 1.0);
// fit->SetParameter(0, 2.0);
// fit->SetParameter(1, 1.0);
// fit->SetParameter(2, 1.0);
// fit->SetParameter(3, 1.0);
// g_uvFm->Fit("fit");
// // g_uvFm->Fit("pol3");
// g_uvFm->SetTitle("#mu_{r} vs F_{m}");
// g_uvFm->GetXaxis()->SetTitle("F_{m}");
// g_uvFm->GetYaxis()->SetTitle("#mu_{r}");
// cout << g_uvFm->GetFunction("fit")->GetChisquare() << endl;
// cout << g_uvFm->GetFunction("fit")->GetNDF() << endl;
// /*Save plots to png and eps files*/
// // c_uvB->Print("../../Plots/uvB/EpoxySteel_uvB_161208.png");
// // c_uvFm->Print("../../Plots/uvB/EpoxySteel_uvFm_161208.png");
return 0;
}
