// Resize any histograms (they all inherit from h1) stat box. Uses NDC coordinates (e.g. 0.65, 0.5)
bool resizeStatBox( TCanvas *c1, TH1* h1, float x1, float y1, float x2=-1, float y2=-1,
TString newname = "", int color = -1, TString oldname = "", int iDebug = 0 )
{
if( !c1 || !h1 ) {cerr<<"Null pointer given to resizeStatBox!"<<endl; return 0;}
c1->Update();
if( oldname.Length() <= 0 ) oldname = "stats";
if( iDebug > 20 ) print_list_of_functions( *h1 );
TPaveStats *stats = (TPaveStats*) h1->FindObject( oldname );
if( iDebug ) cout<<"Found old name (\""<<oldname<<"\"? "<<( stats != 0 )<<endl;
if( (!stats) && newname.Length() > 0 ){ // maybe it was already renamed
stats = (TPaveStats*) h1->FindObject( newname );
if( iDebug ) cout<<"Found new name (\""<<newname<<"\"? "<<(stats != 0)<<endl;
}
if( !stats ) {cerr<<"Can't find stat box"<<endl; return 0;}
stats->SetX1NDC( x1 );
stats->SetY1NDC( y1 );
if( x2 >= 0 ) stats->SetX2NDC( x2 );
if( y2 >= 0 ) stats->SetY2NDC( y2 );
if( newname.Length() > 0 ) {
stats->SetName( newname );
if( iDebug ) cout<<"SetName to "<<newname<<endl;
}
if( color != -1 ) stats->SetTextColor (color);
stats->Draw(); // maybe, just maybe, this will finally make them appear every time, even with draw "same"
return 1;
}
void repositionStatbox(std::string name)
{
// Reposition and resize statbox
TH2F *h = (TH2F*)gDirectory->GetList()->FindObject(name.c_str());
TPaveStats *st = (TPaveStats*)h->GetListOfFunctions()->FindObject("stats");
st->SetX1NDC(0.63);
st->SetY1NDC(0.72);
st->SetX2NDC(0.99);
st->SetY2NDC(0.99);
st->Draw();
}
// überladen: Statbox-Größen manuell eingeben
void drawStatBox(int& step, TH1D* histo, int color = -1, double statboxHeight = 0.1, double statboxSpacing = 0.15){
TPaveStats* statBox = dynamic_cast<TPaveStats*>( histo->GetListOfFunctions()->FindObject("stats") );
if(color == -1) color = step+1;
statBox->SetX1NDC(0.80);
statBox->SetX2NDC(0.99);
statBox->SetY2NDC(0.95-step*statboxSpacing);
statBox->SetY1NDC(0.95-step*statboxSpacing-statboxHeight);
statBox->SetTextColor(color);
statBox->Draw();
step++;
}
// ---------------------------------------------------------------------------------------------
// The SLAC version
bool relocateStatBox (float x1, float y1, float x2=-1, float y2=-1, TString newname = "")
{
TPaveStats *st = (TPaveStats*)gPad->GetPrimitive("stats");
if (0 == st) {cerr<<"Can't find stat box in pad"<<endl; return 0;}
if (newname.Length() > 0) st->SetName(newname);
st->SetX1NDC(x1);
if (x2 != -1) st->SetX2NDC(x2);
st->SetY1NDC(y1);
if (y2 != -1) st->SetY2NDC(y2);
st->Draw();
return 1;
}
// Statistik-Boxen an den Rand der Canvas malen, automatische Größenanpassung
void drawStatBox(TH1D* histo, int& step, int nSteps, double FrameSize, int color = -1, int style = 0){
TPaveStats* statBox = dynamic_cast<TPaveStats*>( histo->GetListOfFunctions()->FindObject("stats") );
double statboxSpacing = FrameSize / nSteps; // gleichmaessiges Aufteilen der Statboxes ueber den Rand
double statboxHeight = 0.1 * statboxSpacing * (9.+1./nSteps);; // 1/10 Abstand, (9/10+Abstand/nSteps) Statbox
if(color == -1) color = step+1;
statBox->SetX1NDC(0.80);
statBox->SetX2NDC(0.99);
statBox->SetY2NDC(0.95-step*statboxSpacing);
statBox->SetY1NDC(0.95-step*statboxSpacing-statboxHeight);
statBox->SetFillColor(color);
statBox->SetFillStyle(style);
statBox->Draw();
step++;
}
int hw1(){
//define landau formula
TF1* myf = new TF1("mylandau",landauf,-3,10,3);
myf->SetParameters(1.0,2,1);
//get ready to initial hist
/*
const int num = 3;
string histvName[num] = {"ld_1e2","ld_1e3","ld_1e4"};
string histvTitle[num] = {"Landau 100 entris","Landau 1000 entris","Landau 10000 entris"};
int histEntries[num] = {100,1000,10000};
//new hist and fill them
TCanvas* canvas = new TCanvas("myc","HW1",800,600);
TMultiGraph* mg = new TMultiGraph();
for(int i = 0;i<num;i++){
TH1F* tmp=new TH1F(histvName[i].c_str(),histvTitle[i].c_str(),100,-3,10);
tmp->FillRandom("mylandau",histEntries[i]);
}
gDirectory->Get(histvName[2])->Draw("e");
gDirectory->Get(histvName[1])->Draw("esame");
gDirectory->Get(histvName[0])->Draw("esame");
*/
TCanvas* canvas = new TCanvas("myc","HW1",800,600);
TH1F* h1 = new TH1F("h1","Landau 100 entries",15,-3,10);
TH1F* h2 = new TH1F("h2","Landau 1000 entries",60,-3,10);
TH1F* h3 = new TH1F("h3","Landau 10000 entries",200,-3,10);
h1->FillRandom("mylandau",100);
h2->FillRandom("mylandau",1000);
h3->FillRandom("mylandau",10000);
h3->GetXaxis()->SetTitle("x");
h3->GetXaxis()->CenterTitle(true);
h3->GetYaxis()->SetTitle("Entries");
h3->GetYaxis()->CenterTitle(true);
h3->SetStats(kFALSE);
TF1* fit3 = new TF1("fit3",landauf,-3,10,3);
fit3->SetParameters(1.0,2,1);
h3->SetMarkerColor(kRed);
h3->SetLineColor(kRed);
fit3->SetLineColor(kRed+2);
h3->Fit("fit3","q","e");
TF1* fit2 = new TF1("fit2",landauf,-3,10,3);
fit2->SetParameters(1.0,2,1);
h2->SetMarkerColor(kBlue);
h2->SetLineColor(kBlue);
fit2->SetLineColor(kBlue+2);
h2->Fit("fit2","q","esame");
TF1* fit1 = new TF1("fit1",landauf,-3,10,3);
fit1->SetParameters(1.0,2,1);
h1->SetMarkerColor(kGreen);
h1->SetLineColor(kGreen);
fit1->SetLineColor(kGreen+2);
h1->Fit("fit1","q","esame");
double h1M = h1->GetMean();
double h1R = h1->GetRMS();
double h1S = h1->GetSkewness();
double h2M = h2->GetMean();
double h2R = h2->GetRMS();
double h2S = h2->GetSkewness();
double h3M = h3->GetMean();
double h3R = h3->GetRMS();
double h3S = h3->GetSkewness();
cout<<setprecision(4);
cout<<"+----+--------+--------+--------+"<<endl;
cout<<"|hist| Mean | RMS |skewness|"<<endl;
cout<<"+----+--------+--------+--------+"<<endl;
cout<<"+ 1 |"<<setw(8)<<h1M<<"|"<<setw(8)<<h1R<<"|"<<setw(8)<<h1S<<"|"<<endl;
cout<<"+----+--------+--------+--------+"<<endl;
cout<<"+ 2 |"<<setw(8)<<h2M<<"|"<<setw(8)<<h2R<<"|"<<setw(8)<<h2S<<"|"<<endl;
cout<<"+----+--------+--------+--------+"<<endl;
cout<<"+ 3 |"<<setw(8)<<h3M<<"|"<<setw(8)<<h3R<<"|"<<setw(8)<<h3S<<"|"<<endl;
cout<<"+----+--------+--------+--------+"<<endl;
TPaveStats *ptstats = new TPaveStats(0.5747126,0.6334746,0.9353448,0.934322,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetLineWidth(2);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(132);
ptstats->AddText(" hist Mean RMS skewness ");
ostringstream os;
cout<<setprecision(4);
os<<" 1 "<<setw(8)<<h1M<<" "<<setw(8)<<h1R<<" "<<setw(8)<<h1S<<" ";
ptstats->AddText(os.str().c_str());
os.str(string());
os<<" 2 "<<setw(8)<<h2M<<" "<<setw(8)<<h2R<<" "<<setw(8)<<h2S<<" ";
ptstats->AddText(os.str().c_str());
os.str(string());
os<<" 3 "<<setw(8)<<h3M<<" "<<setw(8)<<h3R<<" "<<setw(8)<<h3S<<" ";
ptstats->AddText(os.str().c_str());
ptstats->SetOptStat(0);
ptstats->SetOptFit(111);
ptstats->Draw();
h1->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(h1);
TLegend *leg = new TLegend(0.1582915,0.7237762,0.3580402,0.9230769,NULL,"brNDC");
leg->AddEntry(h1,"100 entries","pe");
leg->AddEntry(h2,"1000 entries","pe");
leg->AddEntry(h3,"10000 entries","pe");
leg->SetFillColor(0);
leg->Draw();
//export png
TImage *img = TImage::Create();
img->FromPad(canvas);
img->WriteImage("hw1.png");
//delete img;
//return
return 0;
}
void c()
{
//=========Macro generated from canvas: c/
//========= (Wed Jul 15 12:00:01 2015) by ROOT version6.04/00
TCanvas *c = new TCanvas("c", "",241,92,600,502);
c->SetHighLightColor(2);
c->Range(-162.5,793.0968,162.5,5327.128);
c->SetFillColor(0);
c->SetBorderMode(0);
c->SetBorderSize(2);
c->SetFrameBorderMode(0);
c->SetFrameBorderMode(0);
TH1D *hist__1__1 = new TH1D("hist__1__1","Difference in PMT-constructed position and mc-truth start position",100,-130,130);
hist__1__1->SetBinContent(0,12746);
hist__1__1->SetBinContent(1,1411);
hist__1__1->SetBinContent(2,1544);
hist__1__1->SetBinContent(3,1661);
hist__1__1->SetBinContent(4,1683);
hist__1__1->SetBinContent(5,1784);
hist__1__1->SetBinContent(6,2533);
hist__1__1->SetBinContent(7,2615);
hist__1__1->SetBinContent(8,2743);
hist__1__1->SetBinContent(9,2797);
hist__1__1->SetBinContent(10,2647);
hist__1__1->SetBinContent(11,2958);
hist__1__1->SetBinContent(12,3082);
hist__1__1->SetBinContent(13,2971);
hist__1__1->SetBinContent(14,3232);
hist__1__1->SetBinContent(15,3016);
hist__1__1->SetBinContent(16,3336);
hist__1__1->SetBinContent(17,3712);
hist__1__1->SetBinContent(18,3809);
hist__1__1->SetBinContent(19,3902);
hist__1__1->SetBinContent(20,3742);
hist__1__1->SetBinContent(21,3718);
hist__1__1->SetBinContent(22,3825);
hist__1__1->SetBinContent(23,4008);
hist__1__1->SetBinContent(24,3934);
hist__1__1->SetBinContent(25,4167);
hist__1__1->SetBinContent(26,4107);
hist__1__1->SetBinContent(27,4153);
hist__1__1->SetBinContent(28,4438);
hist__1__1->SetBinContent(29,4538);
hist__1__1->SetBinContent(30,4540);
hist__1__1->SetBinContent(31,4452);
hist__1__1->SetBinContent(32,4157);
hist__1__1->SetBinContent(33,4423);
hist__1__1->SetBinContent(34,4408);
hist__1__1->SetBinContent(35,4611);
hist__1__1->SetBinContent(36,4621);
hist__1__1->SetBinContent(37,4339);
hist__1__1->SetBinContent(38,4458);
hist__1__1->SetBinContent(39,4609);
hist__1__1->SetBinContent(40,4519);
hist__1__1->SetBinContent(41,4430);
hist__1__1->SetBinContent(42,4648);
hist__1__1->SetBinContent(43,4358);
hist__1__1->SetBinContent(44,4503);
hist__1__1->SetBinContent(45,4561);
hist__1__1->SetBinContent(46,4632);
hist__1__1->SetBinContent(47,4566);
hist__1__1->SetBinContent(48,4420);
hist__1__1->SetBinContent(49,4395);
hist__1__1->SetBinContent(50,4628);
hist__1__1->SetBinContent(51,4563);
hist__1__1->SetBinContent(52,4471);
hist__1__1->SetBinContent(53,4539);
hist__1__1->SetBinContent(54,4464);
hist__1__1->SetBinContent(55,4581);
hist__1__1->SetBinContent(56,4531);
hist__1__1->SetBinContent(57,4598);
hist__1__1->SetBinContent(58,4456);
hist__1__1->SetBinContent(59,4341);
hist__1__1->SetBinContent(60,4654);
hist__1__1->SetBinContent(61,4672);
hist__1__1->SetBinContent(62,4511);
hist__1__1->SetBinContent(63,4446);
hist__1__1->SetBinContent(64,4356);
hist__1__1->SetBinContent(65,4380);
hist__1__1->SetBinContent(66,4383);
hist__1__1->SetBinContent(67,4469);
hist__1__1->SetBinContent(68,4606);
hist__1__1->SetBinContent(69,4563);
hist__1__1->SetBinContent(70,4413);
hist__1__1->SetBinContent(71,4701);
hist__1__1->SetBinContent(72,4643);
hist__1__1->SetBinContent(73,4427);
hist__1__1->SetBinContent(74,3933);
hist__1__1->SetBinContent(75,3939);
hist__1__1->SetBinContent(76,4093);
hist__1__1->SetBinContent(77,3913);
hist__1__1->SetBinContent(78,4026);
hist__1__1->SetBinContent(79,3843);
hist__1__1->SetBinContent(80,3641);
hist__1__1->SetBinContent(81,3916);
hist__1__1->SetBinContent(82,4042);
hist__1__1->SetBinContent(83,3798);
hist__1__1->SetBinContent(84,3601);
hist__1__1->SetBinContent(85,3249);
hist__1__1->SetBinContent(86,3146);
hist__1__1->SetBinContent(87,3199);
hist__1__1->SetBinContent(88,3035);
hist__1__1->SetBinContent(89,3021);
hist__1__1->SetBinContent(90,2861);
hist__1__1->SetBinContent(91,2848);
hist__1__1->SetBinContent(92,2819);
hist__1__1->SetBinContent(93,2744);
hist__1__1->SetBinContent(94,2542);
hist__1__1->SetBinContent(95,2198);
hist__1__1->SetBinContent(96,1717);
hist__1__1->SetBinContent(97,1743);
hist__1__1->SetBinContent(98,1585);
hist__1__1->SetBinContent(99,1504);
hist__1__1->SetBinContent(100,1475);
hist__1__1->SetBinContent(101,12137);
hist__1__1->SetEntries(397426);
TPaveStats *ptstats = new TPaveStats(0.78,0.775,0.98,0.935,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
TText *AText = ptstats->AddText("hist__1");
AText->SetTextSize(0.0368);
AText = ptstats->AddText("Entries = 397426 ");
AText = ptstats->AddText("Mean = -0.08411");
AText = ptstats->AddText("Std Dev = 66.16");
ptstats->SetOptStat(1111);
ptstats->SetOptFit(0);
ptstats->Draw();
hist__1__1->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(hist__1__1);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = 926;
color = new TColor(ci, 0, 0, 1, " ", 0.273);
hist__1__1->SetFillColor(ci);
hist__1__1->SetLineColor(4);
hist__1__1->SetLineWidth(2);
hist__1__1->GetXaxis()->SetTitle("#DeltaY [cm]");
hist__1__1->GetXaxis()->SetLabelFont(42);
hist__1__1->GetXaxis()->SetLabelSize(0.05);
hist__1__1->GetXaxis()->SetTitleSize(0.05);
hist__1__1->GetXaxis()->SetTitleOffset(0.93);
hist__1__1->GetXaxis()->SetTitleFont(22);
hist__1__1->GetYaxis()->SetLabelFont(42);
hist__1__1->GetYaxis()->SetLabelSize(0.05);
hist__1__1->GetYaxis()->SetTitleSize(0.035);
hist__1__1->GetYaxis()->SetTitleFont(42);
hist__1__1->GetZaxis()->SetLabelFont(42);
hist__1__1->GetZaxis()->SetLabelSize(0.035);
hist__1__1->GetZaxis()->SetTitleSize(0.035);
hist__1__1->GetZaxis()->SetTitleFont(42);
hist__1__1->Draw("");
TPaveText *pt = new TPaveText(0.15,0.9368947,0.85,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
AText = pt->AddText("Difference in PMT-constructed position and mc-truth start position");
pt->Draw();
c->Modified();
c->cd();
c->SetSelected(c);
}
void plot_PbPb_fit_ratio_eff()
{
//=========Macro generated from canvas: Fit/Fit on Efficiency Ratio
//========= (Mon Aug 29 17:39:14 2016) by ROOT version6.02/10
TCanvas *Fit = new TCanvas("Fit", "Fit on Efficiency Ratio",0,0,800,500);
gStyle->SetOptFit(1);
Fit->Range(1.625,-0.3125,55.375,2.8125);
Fit->SetFillColor(0);
Fit->SetBorderMode(0);
Fit->SetBorderSize(2);
Fit->SetFrameBorderMode(0);
Fit->SetFrameBorderMode(0);
Double_t xAxis1[6] = {7, 10, 15, 20, 30, 50};
TH1D *PbPb_fit1 = new TH1D("PbPb_fit1","",5, xAxis1);
PbPb_fit1->SetBinContent(1,0.8304669);
PbPb_fit1->SetBinContent(2,0.8392248);
PbPb_fit1->SetBinContent(3,0.8551938);
PbPb_fit1->SetBinContent(4,0.8559738);
PbPb_fit1->SetBinContent(5,0.8749088);
PbPb_fit1->SetBinError(1,0.06688803);
PbPb_fit1->SetBinError(2,0.03626477);
PbPb_fit1->SetBinError(3,0.03301156);
PbPb_fit1->SetBinError(4,0.02217037);
PbPb_fit1->SetBinError(5,0.01963207);
PbPb_fit1->SetMinimum(0);
PbPb_fit1->SetMaximum(2.5);
PbPb_fit1->SetEntries(2335.214);
TF1 *myfit1 = new TF1("myfit","pow(10,[0]+[1]*x+x*x*[2])+pow(10,[3]*x*x+[4]*x*x*x+[5])",2,100);
myfit1->SetFillColor(19);
myfit1->SetFillStyle(0);
myfit1->SetMarkerStyle(20);
myfit1->SetLineColor(2);
myfit1->SetLineWidth(2);
myfit1->SetChisquare(5.134234e-08);
myfit1->SetNDF(0);
myfit1->GetXaxis()->SetLabelFont(42);
myfit1->GetXaxis()->SetLabelSize(0.035);
myfit1->GetXaxis()->SetTitleSize(0.035);
myfit1->GetXaxis()->SetTitleFont(42);
myfit1->GetYaxis()->SetLabelFont(42);
myfit1->GetYaxis()->SetLabelSize(0.035);
myfit1->GetYaxis()->SetTitleSize(0.035);
myfit1->GetYaxis()->SetTitleFont(42);
myfit1->SetParameter(0,-0.1325891);
myfit1->SetParError(0,0.07303413);
myfit1->SetParLimits(0,0,0);
myfit1->SetParameter(1,-0.004821085);
myfit1->SetParError(1,0.00252747);
myfit1->SetParLimits(1,0,0);
myfit1->SetParameter(2,0.0001629026);
myfit1->SetParError(2,6.139007e-05);
myfit1->SetParLimits(2,0,0);
myfit1->SetParameter(3,0.001400688);
myfit1->SetParError(3,0.0009400359);
myfit1->SetParLimits(3,0,0);
myfit1->SetParameter(4,-4.985927e-05);
myfit1->SetParError(4,3.268586e-05);
myfit1->SetParLimits(4,0,0);
myfit1->SetParameter(5,-0.9198558);
myfit1->SetParError(5,0.3015603);
myfit1->SetParLimits(5,0,0);
PbPb_fit1->GetListOfFunctions()->Add(myfit1);
TPaveStats *ptstats = new TPaveStats(0.558,0.6,0.9,0.9,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
TText *AText = ptstats->AddText("PbPb_fit");
AText->SetTextSize(0.01971428);
AText = ptstats->AddText("Entries = 2335 ");
AText = ptstats->AddText("Mean = 20.89");
AText = ptstats->AddText("RMS = 11.16");
AText = ptstats->AddText("Underflow = 0");
AText = ptstats->AddText("Overflow = 0");
AText = ptstats->AddText("Integral = 4.256");
AText = ptstats->AddText("#chi^{2} / ndf = 5.134e-08 / 0");
AText = ptstats->AddText("p0 = -0.1326 #pm 0.0730 ");
AText = ptstats->AddText("p1 = -0.004821 #pm 0.002527 ");
AText = ptstats->AddText("p2 = 0.0001629 #pm 0.0000614 ");
AText = ptstats->AddText("p3 = 0.001401 #pm 0.000940 ");
AText = ptstats->AddText("p4 = -4.986e-05 #pm 3.269e-05 ");
AText = ptstats->AddText("p5 = -0.9199 #pm 0.3016 ");
ptstats->SetOptStat(1111111);
ptstats->SetOptFit(111);
ptstats->Draw();
PbPb_fit1->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(PbPb_fit1);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#000099");
PbPb_fit1->SetLineColor(ci);
PbPb_fit1->SetMarkerStyle(20);
PbPb_fit1->GetXaxis()->SetTitle("Gen p_{T}(GeV)");
PbPb_fit1->GetXaxis()->CenterTitle(true);
PbPb_fit1->GetXaxis()->SetLabelFont(42);
PbPb_fit1->GetXaxis()->SetLabelSize(0.035);
PbPb_fit1->GetXaxis()->SetTitleSize(0.035);
PbPb_fit1->GetXaxis()->SetTitleFont(42);
PbPb_fit1->GetYaxis()->SetTitle("Efficiency p_{T} weight / no weight");
PbPb_fit1->GetYaxis()->CenterTitle(true);
PbPb_fit1->GetYaxis()->SetLabelFont(42);
PbPb_fit1->GetYaxis()->SetLabelSize(0.035);
PbPb_fit1->GetYaxis()->SetTitleSize(0.035);
PbPb_fit1->GetYaxis()->SetTitleOffset(1.2);
PbPb_fit1->GetYaxis()->SetTitleFont(42);
PbPb_fit1->GetZaxis()->SetLabelFont(42);
PbPb_fit1->GetZaxis()->SetLabelSize(0.035);
PbPb_fit1->GetZaxis()->SetTitleSize(0.035);
PbPb_fit1->GetZaxis()->SetTitleFont(42);
PbPb_fit1->Draw("");
tex = new TLatex(15,2,"B+ PbPb #sqrt{s}= 5.02 TeV");
tex->SetLineWidth(2);
tex->Draw();
Fit->Modified();
Fit->cd();
Fit->SetSelected(Fit);
}
void QCDHT_cutflow_weighted()
{
//=========Macro generated from canvas: c1/c1
//========= (Wed Feb 24 17:48:49 2016) by ROOT version6.02/05
TCanvas *c1 = new TCanvas("c1", "c1",65,52,700,500);
c1->Range(-1.25,-2351156,16.25,2.11604e+07);
c1->SetFillColor(0);
c1->SetBorderMode(0);
c1->SetBorderSize(2);
c1->SetFrameBorderMode(0);
c1->SetFrameBorderMode(0);
TH1D *cutflow1 = new TH1D("cutflow1","cut flow",14,0.5,14.5);
cutflow1->SetBinContent(1,1.791357e+07);
cutflow1->SetBinContent(2,1.791357e+07);
cutflow1->SetBinContent(3,1.482016e+07);
cutflow1->SetBinContent(4,1.397007e+07);
cutflow1->SetBinContent(5,8159947);
cutflow1->SetBinContent(6,1368453);
cutflow1->SetBinContent(7,12035.83);
cutflow1->SetBinContent(8,10653.79);
cutflow1->SetBinContent(9,4100.401);
cutflow1->SetBinContent(10,4323.004);
cutflow1->SetBinContent(11,1820.803);
cutflow1->SetBinContent(12,374.2397);
cutflow1->SetBinContent(13,35.33991);
cutflow1->SetBinContent(14,214.6597);
cutflow1->SetEntries(1.137194e+08);
TPaveStats *ptstats = new TPaveStats(0.78,0.775,0.98,0.935,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
TText *AText = ptstats->AddText("cutflow");
AText->SetTextSize(0.0368);
AText = ptstats->AddText("Entries = 1.137194e+08");
AText = ptstats->AddText("Mean = 0");
AText = ptstats->AddText("RMS = 0");
ptstats->SetOptStat(1111);
ptstats->SetOptFit(0);
ptstats->Draw();
cutflow1->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(cutflow1);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#000099");
cutflow1->SetLineColor(ci);
cutflow1->GetXaxis()->SetBinLabel(1,"All");
cutflow1->GetXaxis()->SetBinLabel(2,"Trigger");
cutflow1->GetXaxis()->SetBinLabel(3,"p_{T}");
cutflow1->GetXaxis()->SetBinLabel(4,"|#eta|");
cutflow1->GetXaxis()->SetBinLabel(5,"|#Delta#eta|");
cutflow1->GetXaxis()->SetBinLabel(6,"M(jet_{0},jet(1))");
cutflow1->GetXaxis()->SetBinLabel(7,"M(jets)");
cutflow1->GetXaxis()->SetBinLabel(8,"#tau_{21}");
cutflow1->GetXaxis()->SetBinLabel(9,"0b");
cutflow1->GetXaxis()->SetBinLabel(10,"1b");
cutflow1->GetXaxis()->SetBinLabel(11,"2b");
cutflow1->GetXaxis()->SetBinLabel(12,"3b");
cutflow1->GetXaxis()->SetBinLabel(13,"4b");
cutflow1->GetXaxis()->SetBinLabel(14,"3b+HPHP");
cutflow1->GetXaxis()->SetLabelFont(42);
cutflow1->GetXaxis()->SetLabelSize(0.035);
cutflow1->GetXaxis()->SetTitleSize(0.035);
cutflow1->GetXaxis()->SetTitleFont(42);
cutflow1->GetYaxis()->SetLabelFont(42);
cutflow1->GetYaxis()->SetLabelSize(0.035);
cutflow1->GetYaxis()->SetTitleSize(0.035);
cutflow1->GetYaxis()->SetTitleFont(42);
cutflow1->GetZaxis()->SetLabelFont(42);
cutflow1->GetZaxis()->SetLabelSize(0.035);
cutflow1->GetZaxis()->SetTitleSize(0.035);
cutflow1->GetZaxis()->SetTitleFont(42);
cutflow1->Draw("");
TPaveText *pt = new TPaveText(0.4246264,0.94,0.5753736,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
AText = pt->AddText("cut flow");
pt->Draw();
c1->Modified();
c1->cd();
c1->SetSelected(c1);
}
void x_dist_cuts()
{
//=========Macro generated from canvas: c/
//========= (Fri Jul 24 12:08:10 2015) by ROOT version6.04/00
TCanvas *c = new TCanvas("c", "",0,45,600,500);
c->SetHighLightColor(2);
c->Range(-187.5,-2144.494,187.5,19300.44);
c->SetFillColor(0);
c->SetBorderMode(0);
c->SetBorderSize(2);
c->SetFrameBorderMode(0);
c->SetFrameBorderMode(0);
TH1D *hist__1 = new TH1D("hist__1","Worry Later",100,-150,150);
hist__1->SetBinContent(0,37);
hist__1->SetBinContent(1,11);
hist__1->SetBinContent(2,6);
hist__1->SetBinContent(3,9);
hist__1->SetBinContent(4,11);
hist__1->SetBinContent(5,17);
hist__1->SetBinContent(6,29);
hist__1->SetBinContent(7,31);
hist__1->SetBinContent(8,40);
hist__1->SetBinContent(9,65);
hist__1->SetBinContent(10,62);
hist__1->SetBinContent(11,79);
hist__1->SetBinContent(12,132);
hist__1->SetBinContent(13,148);
hist__1->SetBinContent(14,184);
hist__1->SetBinContent(15,187);
hist__1->SetBinContent(16,257);
hist__1->SetBinContent(17,316);
hist__1->SetBinContent(18,354);
hist__1->SetBinContent(19,449);
hist__1->SetBinContent(20,495);
hist__1->SetBinContent(21,568);
hist__1->SetBinContent(22,607);
hist__1->SetBinContent(23,624);
hist__1->SetBinContent(24,686);
hist__1->SetBinContent(25,704);
hist__1->SetBinContent(26,739);
hist__1->SetBinContent(27,738);
hist__1->SetBinContent(28,743);
hist__1->SetBinContent(29,732);
hist__1->SetBinContent(30,730);
hist__1->SetBinContent(31,706);
hist__1->SetBinContent(32,773);
hist__1->SetBinContent(33,716);
hist__1->SetBinContent(34,788);
hist__1->SetBinContent(35,760);
hist__1->SetBinContent(36,900);
hist__1->SetBinContent(37,945);
hist__1->SetBinContent(38,1039);
hist__1->SetBinContent(39,1194);
hist__1->SetBinContent(40,1368);
hist__1->SetBinContent(41,1682);
hist__1->SetBinContent(42,1928);
hist__1->SetBinContent(43,2500);
hist__1->SetBinContent(44,3260);
hist__1->SetBinContent(45,4389);
hist__1->SetBinContent(46,5759);
hist__1->SetBinContent(47,7620);
hist__1->SetBinContent(48,9955);
hist__1->SetBinContent(49,13012);
hist__1->SetBinContent(50,16339);
hist__1->SetBinContent(51,16048);
hist__1->SetBinContent(52,13155);
hist__1->SetBinContent(53,10318);
hist__1->SetBinContent(54,7953);
hist__1->SetBinContent(55,6143);
hist__1->SetBinContent(56,4420);
hist__1->SetBinContent(57,3422);
hist__1->SetBinContent(58,2687);
hist__1->SetBinContent(59,2027);
hist__1->SetBinContent(60,1643);
hist__1->SetBinContent(61,1359);
hist__1->SetBinContent(62,1189);
hist__1->SetBinContent(63,1057);
hist__1->SetBinContent(64,929);
hist__1->SetBinContent(65,840);
hist__1->SetBinContent(66,792);
hist__1->SetBinContent(67,759);
hist__1->SetBinContent(68,719);
hist__1->SetBinContent(69,673);
hist__1->SetBinContent(70,702);
hist__1->SetBinContent(71,769);
hist__1->SetBinContent(72,693);
hist__1->SetBinContent(73,728);
hist__1->SetBinContent(74,719);
hist__1->SetBinContent(75,698);
hist__1->SetBinContent(76,729);
hist__1->SetBinContent(77,695);
hist__1->SetBinContent(78,732);
hist__1->SetBinContent(79,692);
hist__1->SetBinContent(80,620);
hist__1->SetBinContent(81,592);
hist__1->SetBinContent(82,525);
hist__1->SetBinContent(83,463);
hist__1->SetBinContent(84,407);
hist__1->SetBinContent(85,357);
hist__1->SetBinContent(86,300);
hist__1->SetBinContent(87,229);
hist__1->SetBinContent(88,192);
hist__1->SetBinContent(89,154);
hist__1->SetBinContent(90,113);
hist__1->SetBinContent(91,123);
hist__1->SetBinContent(92,73);
hist__1->SetBinContent(93,55);
hist__1->SetBinContent(94,41);
hist__1->SetBinContent(95,27);
hist__1->SetBinContent(96,27);
hist__1->SetBinContent(97,18);
hist__1->SetBinContent(98,14);
hist__1->SetBinContent(99,8);
hist__1->SetBinContent(100,11);
hist__1->SetBinContent(101,55);
hist__1->SetEntries(173117);
TPaveStats *ptstats = new TPaveStats(0.78,0.775,0.98,0.935,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
TText *AText = ptstats->AddText("hist");
AText->SetTextSize(0.0368);
AText = ptstats->AddText("Entries = 173117 ");
AText = ptstats->AddText("Mean = 0.6933");
AText = ptstats->AddText("Std Dev = 34.33");
ptstats->SetOptStat(1111);
ptstats->SetOptFit(0);
ptstats->Draw();
hist__1->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(hist__1);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#000099");
hist__1->SetLineColor(ci);
hist__1->GetXaxis()->SetTitle("#DeltaZ [cm]");
hist__1->GetXaxis()->SetLabelFont(42);
hist__1->GetXaxis()->SetLabelSize(0.035);
hist__1->GetXaxis()->SetTitleSize(0.035);
hist__1->GetXaxis()->SetTitleFont(42);
hist__1->GetYaxis()->SetLabelFont(42);
hist__1->GetYaxis()->SetLabelSize(0.035);
hist__1->GetYaxis()->SetTitleSize(0.035);
hist__1->GetYaxis()->SetTitleFont(42);
hist__1->GetZaxis()->SetLabelFont(42);
hist__1->GetZaxis()->SetLabelSize(0.035);
hist__1->GetZaxis()->SetTitleSize(0.035);
hist__1->GetZaxis()->SetTitleFont(42);
hist__1->Draw("");
TPaveText *pt = new TPaveText(0.3779599,0.9368947,0.6220401,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
AText = pt->AddText("Worry Later");
pt->Draw();
c->Modified();
c->cd();
c->SetSelected(c);
}
int CaloL3(const char *ctype="Calo")
{
bool iSave=false;
double minx=1.0;
LoadStyle();
gStyle->SetOptStat(1);
gStyle->SetOptFit(1);
//string outdir="txtfiles/";
string outdir="junk/";
cout << " # of bins in pt : " << nbins << endl;
//return 0;
TVirtualFitter::SetDefaultFitter("Minuit");
TH1F *hrsp[knj][nbins];
TGraphErrors *grsp[knj], *gcor[knj];
TH1F *hMean[knj], *hRMS[knj];
//TFile *fin = new TFile("jra_hiF_akcalo_ppReco_757p1_dijet.root","r");
//TFile *fin = new TFile("jra_hiF_akcalo_ppReco_757p1_dijet_fulleta_fullpt.root","r");
TFile *fin = new TFile("jra_hiF_akcalo_ppReco_757p1_HcalRespCorrs_v4_00_mc_dijet.root","r");
double refptmin=0;
double refptmax=1000;
TLegend *l0 = getLegend(0.6349664,0.3265565,0.8491835,0.7128335);
l0->SetHeader(Form("Anti-k_{T}, %s",ctype));
TLegend *l1 = getLegend(0.1500646,0.5052258,0.3653091,0.7490942);
l1->SetHeader("");
TF1 *fgaus=0;
TH1F *hrefpt=0, *hjetpt=0;
int njmin=0;
int njmax=knj;
for(int i=njmin; i<njmax; i++){
grsp[i] = new TGraphErrors(nbins);
gcor[i] = new TGraphErrors(nbins);
//grsp[i]->SetName(Form("L3RspVsRefPt_%d",i));
grsp[i]->SetName(Form("L3RspVsRefPt_%s",calgo[i]));
grsp[i]->SetTitle(Form("L3RspVsRefPt %s",calgo[i]));
grsp[i]->SetLineColor(1);
grsp[i]->SetMarkerColor(1);
grsp[i]->SetMarkerStyle(20);
grsp[i]->SetMarkerSize(1.2);
//gcor[i]->SetName(Form("L3CorVsJetPt_%d",i));
gcor[i]->SetName(Form("L3CorVsJetPt_%s",calgo[i]));
gcor[i]->SetTitle(Form("L3CorVsJetPt %s",calgo[i]));
gcor[i]->SetLineColor(1);
gcor[i]->SetMarkerColor(1);
gcor[i]->SetMarkerStyle(20);
gcor[i]->SetMarkerSize(1.2);
hMean [i] = new TH1F(Form("hMean_%d",i),Form("%s%s",calgo[i],ctype),nbins,ptbins);
hMean [i]->SetMarkerColor(1);
hMean [i]->SetMarkerStyle(20);
hMean [i]->SetLineColor(1);
hMean [i]->GetXaxis()->SetMoreLogLabels();
hMean [i]->GetXaxis()->SetNoExponent();
hRMS [i] = new TH1F(Form("hRMS_%d",i),Form("%s%s",calgo[i],ctype),nbins,ptbins);
hRMS [i]->SetMarkerColor(1);
hRMS [i]->SetMarkerStyle(20);
hRMS [i]->SetLineColor(1);
std::ostringstream strs;
//strs << "ak" << i+1 << "CaloJetAnalyzer/";
strs << calgo[i] << "CaloJetAnalyzer/";
cout << " Running for : " << strs.str().c_str() << endl;
for(int j=0; j<nbins; j++){
if(ptbins[j]<minx)continue;
//if(j<2)continue;
std::ostringstream spt;
spt << ptbins[j] << "to" << ptbins[j+1];
//cout << " ptbins : " << ptbins[j] << endl;
std::string sone (strs.str()+"RefPt_Barrel_RefPt"+spt.str());
std::string stwo (strs.str()+"JetPt_Barrel_RefPt"+spt.str());
std::string sthree(strs.str()+"RelRsp_Barrel_RefPt"+spt.str());
hrefpt = (TH1F*)fin->Get(sone.c_str());
hjetpt = (TH1F*)fin->Get(stwo.c_str());
hrsp[i][j] = (TH1F*)fin->Get(sthree.c_str());
//hrsp[i][j]->Rebin(2);
assert(hrefpt->GetEntries()>0 && hjetpt->GetEntries()>0);
// if( ptbins[j] < 30.0 ){
// SetHistRange(hrsp[i][j]);
// }
double refpt =hrefpt->GetMean();
double erefpt =hrefpt->GetMeanError();
double jetpt =hjetpt->GetMean();
double ejetpt =hjetpt->GetMeanError();
hrsp[i][j]->SetLineColor(1);
hrsp[i][j]->SetMarkerColor(1);
hrsp[i][j]->SetMarkerStyle(20);
if(j==0)l0->AddEntry(hrsp[i][j],Form("R = %0.1f",(i+1)*0.1),"p");
//if(i==0)cout <<" refpt : " << refpt << " rawpt : " << jetpt << endl;
//if(j<10)hrsp[i][j]->Rebin(2);
//hrsp[i][j]->Rebin(5);
double norm = hrsp[i][j]->GetMaximumStored();
double rms = hrsp[i][j]->GetRMS();
double peak = hrsp[i][j]->GetMean();
//double peak = (GetPeak(hrsp[i][j]) + hrsp[i][j]->GetMean())/2.;
double epeak = hrsp[i][j]->GetMeanError();
// double tmppeak = GetPeak(hrsp[i][j]);
// fgaus = new TF1("fgaus","gaus", 0.07, 1.2);
// fgaus->SetParameters(norm, peak, 0.15);
// //fgaus->SetParLimits(1, peak - 1.5*rms, peak + 1.5*rms);
// int fitstatus = hrsp[i][j]->Fit(fgaus,"QR");
// peak = (fgaus==0) ? hrsp[i][j]->GetMean() : fgaus->GetParameter(1);
// epeak = (fgaus==0) ? hrsp[i][j]->GetMeanError() : fgaus->GetParError(1);
// if(i==0)fit_dscb(hrsp[i][j], 1.03, minx, 5, Form("%s%s",calgo[i],ctype), 0.17, 1.70);
// else if(i==5)fit_dscb(hrsp[i][j], 1.00, minx, 5, Form("%s%s",calgo[i],ctype), 0.05, 1.45);
// else fit_dscb(hrsp[i][j], 1.03, minx, 5, Form("%s%s",calgo[i],ctype), 0.10, 1.50);
if(i==0)fit_dscb(hrsp[i][j], 1.03, minx, 5, Form("%s%s",calgo[i],ctype), 0.17, 1.70);
else if(i !=5 )fit_dscb(hrsp[i][j], 1.03, minx, 5, Form("%s%s",calgo[i],ctype), 0.10, 1.50);
//! pp reco
// if( i==0 )fit_dscb(hrsp[i][j], 0.85, minx, 5, Form("%s%s",calgo[i],ctype));
// else if(i==1)fit_dscb(hrsp[i][j], 1.03, minx, 5, Form("%s%s",calgo[i],ctype));
// else if(i==2)fit_dscb(hrsp[i][j], 1.00, minx, 5, Form("%s%s",calgo[i],ctype));
// else if(i==3)fit_dscb(hrsp[i][j], 1.05, minx, 5, Form("%s%s",calgo[i],ctype));
// else if(i==4)fit_dscb(hrsp[i][j], 1.05, minx, 5, Form("%s%s",calgo[i],ctype));
// else fit_dscb(hrsp[i][j], 1.00, minx, 5, Form("%s%s",calgo[i],ctype));
//! HI reco
// if( i==0 )fit_dscb(hrsp[i][j], 0.75, minx, 5, Form("ak%d%s",i+1,ctype));
// else if(i==1)fit_dscb(hrsp[i][j], 0.95, minx, 5, Form("ak%d%s",i+1,ctype));
// else if(i==2)fit_dscb(hrsp[i][j], 1.45, minx, 5, Form("ak%d%s",i+1,ctype));
// else if(i==3 || i==4)fit_dscb(hrsp[i][j], 1.50, minx, 4, Form("ak%d%s",i+1,ctype));
// else fit_dscb(hrsp[i][j], 1.30, minx, 5, Form("ak%d%s",i+1,ctype));
TF1* frelrsp = (TF1*)hrsp[i][j]->GetListOfFunctions()->Last();
peak = (frelrsp==0) ? hrsp[i][j]->GetMean() : frelrsp->GetParameter(1);
epeak = (frelrsp==0) ? hrsp[i][j]->GetMeanError() : frelrsp->GetParError(1);
//cout << " peak : " << peak << " mean : " << hrsp[i][j]->GetMean() << " " << frelrsp->GetParameter(1) << endl;
double cor = 1.0/peak;
double ecor = cor*cor*epeak;
//if(i==0)cout <<ptbins[j] << " ratio of " << refpt/jetpt << " " << cor << " " << ftemp->Eval(jetpt) << " " << 1/ ftemp->Eval(jetpt) << endl;
grsp[i]->SetPoint (j, refpt, peak);
grsp[i]->SetPointError(j, erefpt, epeak);
gcor[i]->SetPoint (j, jetpt, cor);
gcor[i]->SetPointError(j, ejetpt, ecor);
}
}
if( gPad )gPad->Close();
//hrsp[0][0]->Draw("p");
//return 0;
int ipad=0;
//! 0 - 20 GeV
TCanvas *c98[knj];
for(int nj=njmin;nj<njmax;nj++){
c98[nj] = new TCanvas(Form("c99_%d",nj),Form("Fine %s Fitting plots",calgo[nj]),100,102,1399,942);
c98[nj]->Divide(6,3,0,0);
//c98[nj]->Divide(10,5,0,0);
//c98[nj]->Divide(6,5,0,0);
//c98[nj]->Divide(7,4,0,0);
//c98[nj]->Divide(6,4,0,0);
ipad=0;
for(int ip=0;ip<nbins;ip++){
if(ptbins[ip]<minx)continue;
c98[nj]->cd(++ipad);
gPad->SetLeftMargin(0.15);
gPad->SetRightMargin(0.01);
gPad->SetBottomMargin(0.15);
gPad->SetLogy();
//if(ipad%10==0)gPad->SetRightMargin(0.02);
//if(ipad%(ip+1)==0)gPad->SetLeftMargin(0.15);
//hrsp[nj][ip]->SetMaximum(hrsp[nj][ip]->GetMaximum() + 0.25*hrsp[nj][ip]->GetMaximum());
hrsp[nj][ip]->SetMaximum(hrsp[nj][ip]->Integral()*2e-00);
hrsp[nj][ip]->SetMinimum(hrsp[nj][ip]->Integral()*1e-07);
hrsp[nj][ip]->SetTitle(0);
hrsp[nj][ip]->GetXaxis()->SetTitle("<reco jet p_{T} / gen jet p_{T}>");
hrsp[nj][ip]->GetXaxis()->SetTitleFont(42);
hrsp[nj][ip]->GetXaxis()->SetLabelFont(42);
hrsp[nj][ip]->GetXaxis()->SetLabelSize(0.08);
hrsp[nj][ip]->GetXaxis()->SetTitleSize(0.07);
hrsp[nj][ip]->GetXaxis()->SetNdivisions(507);
hrsp[nj][ip]->GetYaxis()->SetTitle("");
hrsp[nj][ip]->GetYaxis()->SetTitleFont(42);
hrsp[nj][ip]->GetYaxis()->SetLabelFont(42);
hrsp[nj][ip]->GetYaxis()->SetLabelSize(0.08);
hrsp[nj][ip]->GetYaxis()->SetNdivisions(507);
hrsp[nj][ip]->SetMarkerStyle(20);
hrsp[nj][ip]->SetMarkerColor(1);
hrsp[nj][ip]->SetLineColor(1);
hrsp[nj][ip]->SetMarkerSize(1.1);
hrsp[nj][ip]->Draw("hist");
TF1 *fdscb = (TF1*)hrsp[nj][ip]->GetFunction("fdscb");
if(fdscb){
fdscb->SetLineWidth(1);
fdscb->SetLineColor(2);
fdscb->SetLineColor(2);
fdscb->Draw("lsame");
}
TF1 *fgaus = (TF1*)hrsp[nj][ip]->GetFunction("fgaus");
if(fgaus){
fgaus->SetLineWidth(1);
fgaus->SetLineColor(4);
fgaus->Draw("lsame");
}
// TF1 *fdg = (TF1*)hrsp[nj][ip]->GetFunction("fdgaus");
// if(fdg){
// TF1 *fg1 = (TF1*)hrsp[nj][ip]->GetFunction("fg1");
// TF1 *fg2 = (TF1*)hrsp[nj][ip]->GetFunction("fg2");
// fdg->Draw("lsame");
// fg1->Draw("lsame");
// fg2->Draw("lsame");
// }
std::ostringstream strs;
strs << ptbins[ip] << "< p_{T} (GeV/c) <" << ptbins[ip+1];
std::string spt = strs.str();
if(ipad==1){drawText2(Form("%s%s",calgo[nj],ctype),0.28,0.90,18);
drawText2(spt.c_str(),0.22,0.80,15);
} else drawText2(spt.c_str(),0.17,0.80,15);
}
//if(nj!=0)c98[nj]->Close();
//if( strcmp(calgo[nj],"ak4")!=0 )c98[nj]->Close();
}
//return 0;
cout << endl;
cout << endl;
cout << "Fitting the L3 Response " << endl;
//cout << endl;
//cout << endl;
TF1 *fitrsp=0, *fitcor=0;
std::string fcn_rsp="";
std::string fcn_cor="";
//! Org
// fcn_rsp="[0]-[1]/(pow(log10(x),[2])+[3])";
// fcn_cor="[0]+[1]/(pow(log10(x),[2])+[3])";
//! Using this one
//fcn_rsp="[0]-[1]/(pow(log10(x),[2])+[3])+[4]/x";
//fcn_cor="[0]+[1]/(pow(log10(x),[2])+[3])-[4]/x";
fcn_rsp="[0]-[1]/(pow(log10(x),[2])+[3])+[4]/x+[5]/x/x";
fcn_cor="[0]+[1]/(pow(log10(x),[2])+[3])-[4]/x-[5]/x/x";
// double *rawpt = 0;
// double *genpt = 0;
TCanvas *c11 = new TCanvas("c11","L3 Response",1217,707/*1531,683*/);
c11->Divide(3,2);
TCanvas *c12 = new TCanvas("c12","L3 Corrections",1217,707/*1531,683*/);
c12->Divide(3,2);
for(int i=njmin; i<njmax; i++){
c11->cd(i+1);
fitrsp = new TF1("fitrsp",fcn_rsp.c_str(), minx, grsp[i]->GetX()[grsp[i]->GetN()-1]);
fitrsp->SetParameter(0,1.0);
fitrsp->SetParameter(1,1.0);
fitrsp->SetParameter(2,1.0);
fitrsp->SetParameter(3,1.0);
fitrsp->SetParameter(4,1.0);
fitrsp->SetParameter(5,1.0);
grsp[i]->Fit(fitrsp,"QR");
//if( i<knj )grsp[i]->Fit(fitrsp,"QR","",grsp[i]->GetX()[1],grsp[i]->GetX()[grsp[i]->GetN()-1]);
//else grsp[i]->Fit(fitrsp,"QR");
grsp[i]->SetMaximum(1.2);
grsp[i]->SetMinimum(0.2);
grsp[i]->Draw("ap");
gPad->SetLogx();
gPad->SetLeftMargin(0.15);
gPad->SetBottomMargin(0.15);
gPad->Update();
grsp[i]->GetXaxis()->SetTitle("< Ref p_{T} > (GeV/c)");
grsp[i]->GetXaxis()->SetTitleFont(42);
grsp[i]->GetXaxis()->SetLabelFont(42);
grsp[i]->GetXaxis()->SetLabelSize(0.06);
grsp[i]->GetXaxis()->SetTitleSize(0.06);
grsp[i]->GetXaxis()->SetNdivisions(507);
grsp[i]->GetXaxis()->SetNoExponent();
grsp[i]->GetXaxis()->SetMoreLogLabels();
grsp[i]->GetYaxis()->SetTitle("L3 Response");
grsp[i]->GetYaxis()->SetTitleFont(42);
grsp[i]->GetYaxis()->SetLabelFont(42);
grsp[i]->GetYaxis()->SetLabelSize(0.06);
grsp[i]->GetYaxis()->SetTitleSize(0.06);
grsp[i]->GetYaxis()->SetNdivisions(507);
grsp[i]->GetYaxis()->SetRangeUser(0,1.160);
//grsp[i]->GetXaxis()->SetRangeUser(15,600);
drawText2(Form("%s%s",calgo[i],ctype),0.20,0.78,21);
TPaveStats *strsp = (TPaveStats*)grsp[i]->FindObject("stats");
strsp->SetX1NDC(0.46);
strsp->SetY1NDC(0.18);
strsp->SetX2NDC(0.87);
strsp->SetY2NDC(0.59);
strsp->Draw();
//delete fitrsp;
//continue;
c12->cd(i+1);
//cout << endl;
//cout << endl;
// cout << "Fitting the L3 Corrections " << endl;
// cout << endl;
fitcor = new TF1("fitcor",fcn_cor.c_str(), minx, gcor[i]->GetX()[gcor[i]->GetN()-1]);
fitcor->SetParameter(0,1.0);
fitcor->SetParameter(1,1.0);
fitcor->SetParameter(2,1.0);
fitcor->SetParameter(3,1.0);
fitcor->SetParameter(4,1.0);
fitcor->SetParameter(5,1.0);
gcor[i]->Fit(fitcor,"QR");
//if(i<knj)gcor[i]->Fit(fitcor,"QR","",gcor[i]->GetX()[1],gcor[i]->GetX()[gcor[i]->GetN()-1]);
//else gcor[i]->Fit(fitcor,"QR");
string algn=Form("%d",i+1);
string era="JEC_pp_PYTHIA_TuneCUETP8M1_5020GeV_patch3";
//string era="JEC_pp_HIReco_PYTHIA_TuneCUETP8M1_5020GeV_patch3";
string txtfilename = outdir+era+"_L3Absolute_AK"+algn+"Calo.txt";
ofstream outf(txtfilename.c_str());
//outf.setf(ios::left);
//cout<<"{1 JetEta 1 JetPt "<<fcn_cor<<" Correction L3Absolute}"<<endl;
outf<<"{1 JetEta 1 JetPt "<<fcn_cor<<" Correction L3Absolute}"<<endl;
// cout<<setw(12)<<-3.000 // eta_min
// <<setw(12)<<+3.000 // eta_max
// <<setw(12)<<fitcor->GetNpar()+2 // number of parameters + 2
// <<setw(12)<<7.0 // minimum pT
// <<setw(11)<<1000.0; // maximum pT
outf<<setw(12)<<-3.000 // eta_min
<<setw(12)<<+3.000 // eta_max
<<setw(12)<<fitcor->GetNpar()+2 // number of parameters + 2
<<setw(12)<<7.0 // minimum pT
<<setw(12)<<1000.0; // maximum pT
for(int p=0; p<fitcor->GetNpar(); p++){
//cout<<setw(12)<<fitcor->GetParameter(p); // p0-p4
outf<<setw(12)<<fitcor->GetParameter(p); // p0-p4
}
outf.close();
//cout<<endl;
// if(i==0)cout <<Form("if(i==%d)fl3cor->SetParameters(",i);
// else cout <<Form("else if(i==%d)fl3cor->SetParameters(",i);
// for(int p=0; p<fitcor->GetNpar(); p++){
// if(p==fitcor->GetNpar()-1)cout<<setw(12)<<fitcor->GetParameter(p); // p0-p4
// else cout<<setw(12)<<fitcor->GetParameter(p)<<","; // p0-p4
// }
// cout <<");"<<endl;
// rawpt = gcor[i]->GetX();
// genpt = grsp[i]->GetX();
// for(int ix=0; ix<gcor[i]->GetN(); ix++){
// double cor_fac = fitcor->Eval(rawpt[ix]);
// double corrpt = cor_fac*rawpt[ix];
// //std::cout << " genpt : " << genpt[ix] << " rawpt : " << rawpt[ix] << " corrpt : " << corrpt << " cor_fac : " << cor_fac << std::endl;
// std::cout << genpt[ix] << " : " << (fabs(corrpt - genpt[ix])/genpt[ix])*100. << std::endl;
// }
// gcor[i]->GetXaxis()->SetRangeUser(15,600);
gcor[i]->SetMaximum(3.35);
gcor[i]->SetMinimum(0.89);
gcor[i]->Draw("ap");
gPad->SetLogx();
gPad->SetLeftMargin(0.15);
gPad->SetBottomMargin(0.15);
gPad->Update();
gcor[i]->GetXaxis()->SetTitle("< raw jet p_{T} > (GeV/c)");
gcor[i]->GetXaxis()->SetTitleFont(42);
gcor[i]->GetXaxis()->SetLabelFont(42);
gcor[i]->GetXaxis()->SetLabelSize(0.06);
gcor[i]->GetXaxis()->SetTitleSize(0.06);
gcor[i]->GetXaxis()->SetNdivisions(507);
gcor[i]->GetXaxis()->SetNoExponent();
gcor[i]->GetXaxis()->SetMoreLogLabels();
gcor[i]->GetYaxis()->SetTitle("L3 Correction");
gcor[i]->GetYaxis()->SetTitleFont(42);
gcor[i]->GetYaxis()->SetLabelFont(42);
gcor[i]->GetYaxis()->SetLabelSize(0.06);
gcor[i]->GetYaxis()->SetTitleSize(0.06);
gcor[i]->GetYaxis()->SetNdivisions(507);
gcor[i]->GetXaxis()->SetMoreLogLabels();
gcor[i]->GetXaxis()->SetNoExponent();
drawText2(Form("%s%s",calgo[i],ctype),0.20,0.78,21);
TPaveStats *stcor = (TPaveStats*)gcor[i]->FindObject("stats");
stcor->SetX1NDC(0.48);
stcor->SetY1NDC(0.53);
stcor->SetX2NDC(0.89);
stcor->SetY2NDC(0.89);
stcor->Draw();
// cout << endl;
//delete fitcor;
}
if(iSave){
c12->SaveAs("CorrectionPlots/L3Absolute_Corrections_ppReco_CaloJets_757p1.pdf");
}
TFile *fout = new TFile("l3calo_input.root","RECREATE");
for(int nj=njmin;nj<njmax;nj++){
fout->mkdir(Form("%sCaloJetAnalyzer",calgo[nj]),Form("%sCaloJetAnalyzer",calgo[nj]));
fout->cd(Form("%sCaloJetAnalyzer",calgo[nj]));
grsp[nj]->Write();
gcor[nj]->Write();
fout->cd("../");
}
fout->Close();
return 0;
}
void BSA_InCoherent_t()
{
//=========Macro generated from canvas: c17/Incoherent channe
//========= (Tue Aug 21 12:00:26 2018) by ROOT version6.09/01
TCanvas *c17 = new TCanvas("c17", "Incoherent channe",0,0,1200,400);
gStyle->SetOptFit(1);
c17->Range(0,0,1,1);
c17->SetFillColor(0);
c17->SetBorderMode(0);
c17->SetBorderSize(2);
c17->SetLeftMargin(0.025);
c17->SetRightMargin(0.02325581);
c17->SetFrameBorderMode(0);
// ------------>Primitives in pad: c17_1
TPad *c17_1 = new TPad("c17_1", "Incoherent channe_1",0,0,0.262936,0.9);
c17_1->Draw();
c17_1->cd();
c17_1->Range(-32.86883,-0.4888889,353.6809,0.4);
c17_1->SetFillColor(0);
c17_1->SetBorderMode(0);
c17_1->SetBorderSize(2);
c17_1->SetGridx();
c17_1->SetGridy();
c17_1->SetRightMargin(0);
c17_1->SetTopMargin(0);
c17_1->SetFrameBorderMode(0);
c17_1->SetFrameBorderMode(0);
TMultiGraph *multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("0.05<-t<0.21 [GeV^{2}/c^{2}]");
Double_t Graph_fx1154[9] = {
21.59954,
59.03178,
96.06391,
141.7703,
179.2242,
222.8541,
259.5793,
303.38,
337.8675};
Double_t Graph_fy1154[9] = {
0.1020231,
0.03114899,
0.2309077,
0.1341245,
0.01684177,
-0.1267998,
-0.2824572,
-0.2260703,
-0.0186068};
Double_t Graph_fex1154[9] = {
0,
0,
0,
0,
0,
0,
0,
0,
0};
Double_t Graph_fey1154[9] = {
0.04242532,
0.05165487,
0.05483741,
0.07646672,
0.1006626,
0.09089962,
0.05905575,
0.04462244,
0.05521722};
TGraphErrors *gre = new TGraphErrors(9,Graph_fx1154,Graph_fy1154,Graph_fex1154,Graph_fey1154);
gre->SetName("Graph");
gre->SetTitle("Graph");
gre->SetFillColor(1);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#00ff00");
gre->SetMarkerColor(ci);
gre->SetMarkerStyle(21);
TH1F *Graph_Graph1154 = new TH1F("Graph_Graph1154","Graph",100,0,369.4943);
Graph_Graph1154->SetMinimum(-0.4042387);
Graph_Graph1154->SetMaximum(0.3484709);
Graph_Graph1154->SetDirectory(0);
Graph_Graph1154->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1154->SetLineColor(ci);
Graph_Graph1154->GetXaxis()->SetLabelFont(42);
Graph_Graph1154->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1154->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1154->GetXaxis()->SetTitleFont(42);
Graph_Graph1154->GetYaxis()->SetLabelFont(42);
Graph_Graph1154->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1154->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1154->GetYaxis()->SetTitleFont(42);
Graph_Graph1154->GetZaxis()->SetLabelFont(42);
Graph_Graph1154->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1154->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1154->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1154);
TF1 *myfit1155 = new TF1("myfit","[0]*sin(x*3.1416/180.0)/(1 + [1]*cos(x*3.1416/180.0))",0,369.4943);
myfit1155->SetFillColor(19);
myfit1155->SetFillStyle(0);
myfit1155->SetLineColor(2);
myfit1155->SetLineWidth(2);
myfit1155->SetChisquare(12.03234);
myfit1155->SetNDF(7);
myfit1155->GetXaxis()->SetLabelFont(42);
myfit1155->GetXaxis()->SetLabelSize(0.035);
myfit1155->GetXaxis()->SetTitleSize(0.035);
myfit1155->GetXaxis()->SetTitleFont(42);
myfit1155->GetYaxis()->SetLabelFont(42);
myfit1155->GetYaxis()->SetLabelSize(0.035);
myfit1155->GetYaxis()->SetTitleSize(0.035);
myfit1155->GetYaxis()->SetTitleFont(42);
myfit1155->SetParameter(0,0.2160608);
myfit1155->SetParError(0,0.02670852);
myfit1155->SetParLimits(0,0,0);
myfit1155->SetParameter(1,0.2378828);
myfit1155->SetParError(1,0.2292802);
myfit1155->SetParLimits(1,0,0);
gre->GetListOfFunctions()->Add(myfit1155);
TPaveStats *ptstats = new TPaveStats(0.27,0.74,0.9,0.9,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
TText *AText = ptstats->AddText("#chi^{2} / ndf = 12.03 / 7");
AText = ptstats->AddText("p0 = 0.2161 #pm 0.02671 ");
AText = ptstats->AddText("p1 = 0.2379 #pm 0.2293 ");
ptstats->SetOptStat(0);
ptstats->SetOptFit(111);
ptstats->Draw();
gre->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(gre->GetListOfFunctions());
multigraph->Add(gre,"");
multigraph->Draw("AP");
multigraph->GetXaxis()->SetTitle("#phi [deg.]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetTitle("A_{LU}");
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
gPad->Modified();
multigraph->GetYaxis()->SetRangeUser(-0.35, 0.35);
TPaveText *pt = new TPaveText(0.2292451,0.9268996,0.7707549,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
AText = pt->AddText("0.05<-t<0.21 [GeV^{2}/c^{2}]");
pt->Draw();
c17_1->Modified();
c17->cd();
// ------------>Primitives in pad: c17_2
TPad *c17_2 = new TPad("c17_2", "Incoherent channe_2",0.262936,0,0.5008721,0.9);
c17_2->Draw();
c17_2->cd();
c17_2->Range(5.895345,-0.4888889,352.6198,0.4);
c17_2->SetFillColor(0);
c17_2->SetBorderMode(0);
c17_2->SetBorderSize(2);
c17_2->SetGridx();
c17_2->SetGridy();
c17_2->SetLeftMargin(0);
c17_2->SetRightMargin(0);
c17_2->SetTopMargin(0);
c17_2->SetFrameBorderMode(0);
c17_2->SetFrameBorderMode(0);
multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("0.21<-t<0.36 [GeV^{2}/c^{2}]");
Double_t Graph_fx1157[9] = {
21.65555,
61.41754,
96.1904,
139.625,
182.9173,
222.6445,
259.1526,
302.149,
336.8596};
Double_t Graph_fy1157[9] = {
0.04372413,
0.1019325,
0.1737435,
0.1708402,
0.2316287,
-0.1435087,
-0.1227365,
-0.1969296,
-0.07714485};
Double_t Graph_fex1157[9] = {
0,
0,
0,
0,
0,
0,
0,
0,
0};
Double_t Graph_fey1157[9] = {
0.0517677,
0.05885852,
0.0646232,
0.09789667,
0.1185262,
0.09557657,
0.06953456,
0.04921871,
0.06282415};
gre = new TGraphErrors(9,Graph_fx1157,Graph_fy1157,Graph_fex1157,Graph_fey1157);
gre->SetName("Graph");
gre->SetTitle("Graph");
gre->SetFillColor(1);
ci = TColor::GetColor("#00ff00");
gre->SetMarkerColor(ci);
gre->SetMarkerStyle(21);
TH1F *Graph_Graph1157 = new TH1F("Graph_Graph1157","Graph",100,0,368.38);
Graph_Graph1157->SetMinimum(-0.3057786);
Graph_Graph1157->SetMaximum(0.4097851);
Graph_Graph1157->SetDirectory(0);
Graph_Graph1157->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1157->SetLineColor(ci);
Graph_Graph1157->GetXaxis()->SetLabelFont(42);
Graph_Graph1157->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1157->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1157->GetXaxis()->SetTitleFont(42);
Graph_Graph1157->GetYaxis()->SetLabelFont(42);
Graph_Graph1157->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1157->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1157->GetYaxis()->SetTitleFont(42);
Graph_Graph1157->GetZaxis()->SetLabelFont(42);
Graph_Graph1157->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1157->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1157->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1157);
TF1 *myfit1158 = new TF1("myfit","[0]*sin(x*3.1416/180.0)/(1 + [1]*cos(x*3.1416/180.0))",0,368.38);
myfit1158->SetFillColor(19);
myfit1158->SetFillStyle(0);
myfit1158->SetLineColor(2);
myfit1158->SetLineWidth(2);
myfit1158->SetChisquare(6.955655);
myfit1158->SetNDF(7);
myfit1158->GetXaxis()->SetLabelFont(42);
myfit1158->GetXaxis()->SetLabelSize(0.035);
myfit1158->GetXaxis()->SetTitleSize(0.035);
myfit1158->GetXaxis()->SetTitleFont(42);
myfit1158->GetYaxis()->SetLabelFont(42);
myfit1158->GetYaxis()->SetLabelSize(0.035);
myfit1158->GetYaxis()->SetTitleSize(0.035);
myfit1158->GetYaxis()->SetTitleFont(42);
myfit1158->SetParameter(0,0.1721707);
myfit1158->SetParError(0,0.0324291);
myfit1158->SetParLimits(0,0,0);
myfit1158->SetParameter(1,0.001708824);
myfit1158->SetParError(1,0.4101565);
myfit1158->SetParLimits(1,0,0);
gre->GetListOfFunctions()->Add(myfit1158);
ptstats = new TPaveStats(0.27,0.74,0.9,0.9,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
AText = ptstats->AddText("#chi^{2} / ndf = 6.956 / 7");
AText = ptstats->AddText("p0 = 0.1722 #pm 0.03243 ");
AText = ptstats->AddText("p1 = 0.001709 #pm 0.4102 ");
ptstats->SetOptStat(0);
ptstats->SetOptFit(111);
ptstats->Draw();
gre->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(gre->GetListOfFunctions());
multigraph->Add(gre,"");
multigraph->Draw("AP");
multigraph->GetXaxis()->SetTitle("#phi [deg.]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
gPad->Modified();
multigraph->GetYaxis()->SetRangeUser(-0.35, 0.35);
pt = new TPaveText(0.2194177,0.9304839,0.7805823,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
AText = pt->AddText("0.21<-t<0.36 [GeV^{2}/c^{2}]");
pt->Draw();
c17_2->Modified();
c17->cd();
// ------------>Primitives in pad: c17_3
TPad *c17_3 = new TPad("c17_3", "Incoherent channe_3",0.5008721,0,0.7388081,0.9);
c17_3->Draw();
c17_3->cd();
c17_3->Range(7.121854,-0.4888889,349.9594,0.4);
c17_3->SetFillColor(0);
c17_3->SetBorderMode(0);
c17_3->SetBorderSize(2);
c17_3->SetGridx();
c17_3->SetGridy();
c17_3->SetLeftMargin(0);
c17_3->SetRightMargin(0);
c17_3->SetTopMargin(0);
c17_3->SetFrameBorderMode(0);
c17_3->SetFrameBorderMode(0);
multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("0.36<-t<0.67 [GeV^{2}/c^{2}]");
Double_t Graph_fx1160[9] = {
22.70538,
60.97743,
94.23011,
138.1892,
183.0975,
222.79,
259.8951,
302.6808,
334.3759};
Double_t Graph_fy1160[9] = {
0.1761498,
0.1945204,
0.1034809,
-0.09380414,
0.2967375,
-0.1094475,
-0.1583282,
-0.1819771,
-0.09543453};
Double_t Graph_fex1160[9] = {
0,
0,
0,
0,
0,
0,
0,
0,
0};
Double_t Graph_fey1160[9] = {
0.06358895,
0.0584754,
0.07622601,
0.115077,
0.1268505,
0.1278475,
0.07389284,
0.04701977,
0.07415904};
gre = new TGraphErrors(9,Graph_fx1160,Graph_fy1160,Graph_fex1160,Graph_fey1160);
gre->SetName("Graph");
gre->SetTitle("Graph");
gre->SetFillColor(1);
ci = TColor::GetColor("#00ff00");
gre->SetMarkerColor(ci);
gre->SetMarkerStyle(21);
TH1F *Graph_Graph1160 = new TH1F("Graph_Graph1160","Graph",100,0,365.5429);
Graph_Graph1160->SetMinimum(-0.3033832);
Graph_Graph1160->SetMaximum(0.4896763);
Graph_Graph1160->SetDirectory(0);
Graph_Graph1160->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1160->SetLineColor(ci);
Graph_Graph1160->GetXaxis()->SetLabelFont(42);
Graph_Graph1160->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1160->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1160->GetXaxis()->SetTitleFont(42);
Graph_Graph1160->GetYaxis()->SetLabelFont(42);
Graph_Graph1160->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1160->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1160->GetYaxis()->SetTitleFont(42);
Graph_Graph1160->GetZaxis()->SetLabelFont(42);
Graph_Graph1160->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1160->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1160->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1160);
TF1 *myfit1161 = new TF1("myfit","[0]*sin(x*3.1416/180.0)/(1 + [1]*cos(x*3.1416/180.0))",0,365.5429);
myfit1161->SetFillColor(19);
myfit1161->SetFillStyle(0);
myfit1161->SetLineColor(2);
myfit1161->SetLineWidth(2);
myfit1161->SetChisquare(8.926216);
myfit1161->SetNDF(7);
myfit1161->GetXaxis()->SetLabelFont(42);
myfit1161->GetXaxis()->SetLabelSize(0.035);
myfit1161->GetXaxis()->SetTitleSize(0.035);
myfit1161->GetXaxis()->SetTitleFont(42);
myfit1161->GetYaxis()->SetLabelFont(42);
myfit1161->GetYaxis()->SetLabelSize(0.035);
myfit1161->GetYaxis()->SetTitleSize(0.035);
myfit1161->GetYaxis()->SetTitleFont(42);
myfit1161->SetParameter(0,0.1371093);
myfit1161->SetParError(0,0.03383483);
myfit1161->SetParLimits(0,0,0);
myfit1161->SetParameter(1,-0.6800102);
myfit1161->SetParError(1,0.1832383);
myfit1161->SetParLimits(1,0,0);
gre->GetListOfFunctions()->Add(myfit1161);
ptstats = new TPaveStats(0.27,0.74,0.9,0.9,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
AText = ptstats->AddText("#chi^{2} / ndf = 8.926 / 7");
AText = ptstats->AddText("p0 = 0.1371 #pm 0.03383 ");
AText = ptstats->AddText("p1 = -0.68 #pm 0.1832 ");
ptstats->SetOptStat(0);
ptstats->SetOptFit(111);
ptstats->Draw();
gre->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(gre->GetListOfFunctions());
multigraph->Add(gre,"");
multigraph->Draw("AP");
multigraph->GetXaxis()->SetTitle("#phi [deg.]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
gPad->Modified();
multigraph->GetYaxis()->SetRangeUser(-0.35, 0.35);
pt = new TPaveText(0.2194177,0.9304839,0.7805823,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
AText = pt->AddText("0.36<-t<0.67 [GeV^{2}/c^{2}]");
pt->Draw();
c17_3->Modified();
c17->cd();
// ------------>Primitives in pad: c17_4
TPad *c17_4 = new TPad("c17_4", "Incoherent channe_4",0.7388081,0,0.9767442,0.9);
c17_4->Draw();
c17_4->cd();
c17_4->Range(7.951117,-0.4888889,350.3641,0.4);
c17_4->SetFillColor(0);
c17_4->SetBorderMode(0);
c17_4->SetBorderSize(2);
c17_4->SetGridx();
c17_4->SetGridy();
c17_4->SetLeftMargin(0);
c17_4->SetRightMargin(0);
c17_4->SetTopMargin(0);
c17_4->SetFrameBorderMode(0);
c17_4->SetFrameBorderMode(0);
multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("0.67<-t<0.90 [GeV^{2}/c^{2}]");
Double_t Graph_fx1163[9] = {
23.51534,
59.3759,
90.55341,
131.0631,
176.9476,
229.0094,
260.8643,
304.7922,
334.7999};
Double_t Graph_fy1163[9] = {
0.0123906,
0.195403,
0.18963,
-0.1677207,
-0.8232306,
-0.7760588,
-0.2579708,
-0.2579435,
-0.08614626};
Double_t Graph_fex1163[9] = {
0,
0,
0,
0,
0,
0,
0,
0,
0};
Double_t Graph_fey1163[9] = {
0.09527872,
0.1168373,
0.1651642,
0.3137499,
0.3909388,
0.4596559,
0.1665094,
0.09339681,
0.09961735};
gre = new TGraphErrors(9,Graph_fx1163,Graph_fy1163,Graph_fex1163,Graph_fey1163);
gre->SetName("Graph");
gre->SetTitle("Graph");
gre->SetFillColor(1);
ci = TColor::GetColor("#00ff00");
gre->SetMarkerColor(ci);
gre->SetMarkerStyle(21);
TH1F *Graph_Graph1163 = new TH1F("Graph_Graph1163","Graph",100,0,365.9283);
Graph_Graph1163->SetMinimum(-1.394766);
Graph_Graph1163->SetMaximum(0.5138451);
Graph_Graph1163->SetDirectory(0);
Graph_Graph1163->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1163->SetLineColor(ci);
Graph_Graph1163->GetXaxis()->SetLabelFont(42);
Graph_Graph1163->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1163->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1163->GetXaxis()->SetTitleFont(42);
Graph_Graph1163->GetYaxis()->SetLabelFont(42);
Graph_Graph1163->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1163->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1163->GetYaxis()->SetTitleFont(42);
Graph_Graph1163->GetZaxis()->SetLabelFont(42);
Graph_Graph1163->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1163->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1163->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1163);
TF1 *myfit1164 = new TF1("myfit","[0]*sin(x*3.1416/180.0)/(1 + [1]*cos(x*3.1416/180.0))",0,365.9283);
myfit1164->SetFillColor(19);
myfit1164->SetFillStyle(0);
myfit1164->SetLineColor(2);
myfit1164->SetLineWidth(2);
myfit1164->SetChisquare(8.810638);
myfit1164->SetNDF(7);
myfit1164->GetXaxis()->SetLabelFont(42);
myfit1164->GetXaxis()->SetLabelSize(0.035);
myfit1164->GetXaxis()->SetTitleSize(0.035);
myfit1164->GetXaxis()->SetTitleFont(42);
myfit1164->GetYaxis()->SetLabelFont(42);
myfit1164->GetYaxis()->SetLabelSize(0.035);
myfit1164->GetYaxis()->SetTitleSize(0.035);
myfit1164->GetYaxis()->SetTitleFont(42);
myfit1164->SetParameter(0,0.2352442);
myfit1164->SetParError(0,0.0807391);
myfit1164->SetParLimits(0,0,0);
myfit1164->SetParameter(1,0.003684647);
myfit1164->SetParError(1,0.5673774);
myfit1164->SetParLimits(1,0,0);
gre->GetListOfFunctions()->Add(myfit1164);
ptstats = new TPaveStats(0.27,0.74,0.9,0.9,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
AText = ptstats->AddText("#chi^{2} / ndf = 8.811 / 7");
AText = ptstats->AddText("p0 = 0.2352 #pm 0.08074 ");
AText = ptstats->AddText("p1 = 0.003685 #pm 0.5674 ");
ptstats->SetOptStat(0);
ptstats->SetOptFit(111);
ptstats->Draw();
gre->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(gre->GetListOfFunctions());
multigraph->Add(gre,"");
multigraph->Draw("AP");
multigraph->GetXaxis()->SetTitle("#phi [deg.]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
gPad->Modified();
multigraph->GetYaxis()->SetRangeUser(-0.35, 0.35);
pt = new TPaveText(0.2194177,0.9304839,0.7805823,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
AText = pt->AddText("0.67<-t<0.90 [GeV^{2}/c^{2}]");
pt->Draw();
c17_4->Modified();
c17->cd();
c17->Modified();
c17->cd();
c17->SetSelected(c17);
}
void BSA_InCoherent_xB()
{
//=========Macro generated from canvas: c16/Incoherent channe
//========= (Tue Aug 21 12:00:25 2018) by ROOT version6.09/01
TCanvas *c16 = new TCanvas("c16", "Incoherent channe",0,0,1200,400);
gStyle->SetOptFit(1);
c16->Range(0,0,1,1);
c16->SetFillColor(0);
c16->SetBorderMode(0);
c16->SetBorderSize(2);
c16->SetLeftMargin(0.025);
c16->SetRightMargin(0.02325581);
c16->SetFrameBorderMode(0);
// ------------>Primitives in pad: c16_1
TPad *c16_1 = new TPad("c16_1", "Incoherent channe_1",0,0,0.262936,0.9);
c16_1->Draw();
c16_1->cd();
c16_1->Range(-29.34641,-0.4888889,350.2581,0.4);
c16_1->SetFillColor(0);
c16_1->SetBorderMode(0);
c16_1->SetBorderSize(2);
c16_1->SetGridx();
c16_1->SetGridy();
c16_1->SetRightMargin(0);
c16_1->SetTopMargin(0);
c16_1->SetFrameBorderMode(0);
c16_1->SetFrameBorderMode(0);
TMultiGraph *multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("0.12<x_{B}<0.21 ");
Double_t Graph_fx1134[9] = {
24.14331,
62.70082,
96.14623,
141.7289,
180.8857,
220.5862,
258.7856,
303.3336,
334.7288};
Double_t Graph_fy1134[9] = {
0.1330996,
-0.03883743,
0.1388413,
0.2120595,
0.03749243,
-0.09129921,
-0.1273289,
-0.1643618,
0.03103783};
Double_t Graph_fex1134[9] = {
0,
0,
0,
0,
0,
0,
0,
0,
0};
Double_t Graph_fey1134[9] = {
0.07077824,
0.06445297,
0.05939535,
0.06833722,
0.08559024,
0.08401363,
0.06136237,
0.05845869,
0.08718727};
TGraphErrors *gre = new TGraphErrors(9,Graph_fx1134,Graph_fy1134,Graph_fex1134,Graph_fey1134);
gre->SetName("Graph");
gre->SetTitle("Graph");
gre->SetFillColor(1);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#00ff00");
gre->SetMarkerColor(ci);
gre->SetMarkerStyle(21);
TH1F *Graph_Graph1134 = new TH1F("Graph_Graph1134","Graph",100,0,365.7873);
Graph_Graph1134->SetMinimum(-0.2731422);
Graph_Graph1134->SetMaximum(0.3307184);
Graph_Graph1134->SetDirectory(0);
Graph_Graph1134->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1134->SetLineColor(ci);
Graph_Graph1134->GetXaxis()->SetLabelFont(42);
Graph_Graph1134->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1134->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1134->GetXaxis()->SetTitleFont(42);
Graph_Graph1134->GetYaxis()->SetLabelFont(42);
Graph_Graph1134->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1134->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1134->GetYaxis()->SetTitleFont(42);
Graph_Graph1134->GetZaxis()->SetLabelFont(42);
Graph_Graph1134->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1134->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1134->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1134);
TF1 *myfit1135 = new TF1("myfit","[0]*sin(x*3.1416/180.0)/(1 + [1]*cos(x*3.1416/180.0))",0,365.7873);
myfit1135->SetFillColor(19);
myfit1135->SetFillStyle(0);
myfit1135->SetLineColor(2);
myfit1135->SetLineWidth(2);
myfit1135->SetChisquare(10.13595);
myfit1135->SetNDF(7);
myfit1135->GetXaxis()->SetLabelFont(42);
myfit1135->GetXaxis()->SetLabelSize(0.035);
myfit1135->GetXaxis()->SetTitleSize(0.035);
myfit1135->GetXaxis()->SetTitleFont(42);
myfit1135->GetYaxis()->SetLabelFont(42);
myfit1135->GetYaxis()->SetLabelSize(0.035);
myfit1135->GetYaxis()->SetTitleSize(0.035);
myfit1135->GetYaxis()->SetTitleFont(42);
myfit1135->SetParameter(0,0.1205343);
myfit1135->SetParError(0,0.03387604);
myfit1135->SetParLimits(0,0,0);
myfit1135->SetParameter(1,0.6858633);
myfit1135->SetParError(1,0.2646153);
myfit1135->SetParLimits(1,0,0);
gre->GetListOfFunctions()->Add(myfit1135);
TPaveStats *ptstats = new TPaveStats(0.27,0.74,0.9,0.9,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
TText *AText = ptstats->AddText("#chi^{2} / ndf = 10.14 / 7");
AText = ptstats->AddText("p0 = 0.1205 #pm 0.03388 ");
AText = ptstats->AddText("p1 = 0.6859 #pm 0.2646 ");
ptstats->SetOptStat(0);
ptstats->SetOptFit(111);
ptstats->Draw();
gre->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(gre->GetListOfFunctions());
multigraph->Add(gre,"");
multigraph->Draw("AP");
multigraph->GetXaxis()->SetTitle("#phi [deg.]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetTitle("A_{LU}");
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
gPad->Modified();
multigraph->GetYaxis()->SetRangeUser(-0.35, 0.35);
c16_1->Modified();
c16->cd();
// ------------>Primitives in pad: c16_2
TPad *c16_2 = new TPad("c16_2", "Incoherent channe_2",0.262936,0,0.5008721,0.9);
c16_2->Draw();
c16_2->cd();
c16_2->Range(5.85529,-0.4888889,352.0975,0.4);
c16_2->SetFillColor(0);
c16_2->SetBorderMode(0);
c16_2->SetBorderSize(2);
c16_2->SetGridx();
c16_2->SetGridy();
c16_2->SetLeftMargin(0);
c16_2->SetRightMargin(0);
c16_2->SetTopMargin(0);
c16_2->SetFrameBorderMode(0);
c16_2->SetFrameBorderMode(0);
multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("0.21<x_{B}<0.31 ");
Double_t Graph_fx1137[9] = {
21.59357,
60.08662,
95.00395,
140.4897,
181.0443,
224.9978,
260.1264,
302.891,
336.3592};
Double_t Graph_fy1137[9] = {
0.0984356,
0.1891659,
0.2419886,
0.08897874,
0.2077051,
-0.3071227,
-0.1933684,
-0.2875592,
-0.08265727};
Double_t Graph_fex1137[9] = {
0,
0,
0,
0,
0,
0,
0,
0,
0};
Double_t Graph_fey1137[9] = {
0.04446748,
0.058618,
0.06898028,
0.1082398,
0.135301,
0.1152745,
0.0716065,
0.04680018,
0.05547391};
gre = new TGraphErrors(9,Graph_fx1137,Graph_fy1137,Graph_fex1137,Graph_fey1137);
gre->SetName("Graph");
gre->SetTitle("Graph");
gre->SetFillColor(1);
ci = TColor::GetColor("#00ff00");
gre->SetMarkerColor(ci);
gre->SetMarkerStyle(21);
TH1F *Graph_Graph1137 = new TH1F("Graph_Graph1137","Graph",100,0,367.8358);
Graph_Graph1137->SetMinimum(-0.4989376);
Graph_Graph1137->SetMaximum(0.4195464);
Graph_Graph1137->SetDirectory(0);
Graph_Graph1137->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1137->SetLineColor(ci);
Graph_Graph1137->GetXaxis()->SetLabelFont(42);
Graph_Graph1137->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1137->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1137->GetXaxis()->SetTitleFont(42);
Graph_Graph1137->GetYaxis()->SetLabelFont(42);
Graph_Graph1137->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1137->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1137->GetYaxis()->SetTitleFont(42);
Graph_Graph1137->GetZaxis()->SetLabelFont(42);
Graph_Graph1137->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1137->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1137->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1137);
TF1 *myfit1138 = new TF1("myfit","[0]*sin(x*3.1416/180.0)/(1 + [1]*cos(x*3.1416/180.0))",0,367.8358);
myfit1138->SetFillColor(19);
myfit1138->SetFillStyle(0);
myfit1138->SetLineColor(2);
myfit1138->SetLineWidth(2);
myfit1138->SetChisquare(7.336971);
myfit1138->SetNDF(7);
myfit1138->GetXaxis()->SetLabelFont(42);
myfit1138->GetXaxis()->SetLabelSize(0.035);
myfit1138->GetXaxis()->SetTitleSize(0.035);
myfit1138->GetXaxis()->SetTitleFont(42);
myfit1138->GetYaxis()->SetLabelFont(42);
myfit1138->GetYaxis()->SetLabelSize(0.035);
myfit1138->GetYaxis()->SetTitleSize(0.035);
myfit1138->GetYaxis()->SetTitleFont(42);
myfit1138->SetParameter(0,0.2508538);
myfit1138->SetParError(0,0.03583101);
myfit1138->SetParLimits(0,0,0);
myfit1138->SetParameter(1,-0.1422748);
myfit1138->SetParError(1,0.2426713);
myfit1138->SetParLimits(1,0,0);
gre->GetListOfFunctions()->Add(myfit1138);
ptstats = new TPaveStats(0.27,0.74,0.9,0.9,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
AText = ptstats->AddText("#chi^{2} / ndf = 7.337 / 7");
AText = ptstats->AddText("p0 = 0.2509 #pm 0.03583 ");
AText = ptstats->AddText("p1 = -0.1423 #pm 0.2427 ");
ptstats->SetOptStat(0);
ptstats->SetOptFit(111);
ptstats->Draw();
gre->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(gre->GetListOfFunctions());
multigraph->Add(gre,"");
multigraph->Draw("AP");
multigraph->GetXaxis()->SetTitle("#phi [deg.]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
gPad->Modified();
multigraph->GetYaxis()->SetRangeUser(-0.35, 0.35);
c16_2->Modified();
c16->cd();
// ------------>Primitives in pad: c16_3
TPad *c16_3 = new TPad("c16_3", "Incoherent channe_3",0.5008721,0,0.7388081,0.9);
c16_3->Draw();
c16_3->cd();
c16_3->Range(4.159976,-0.4888889,357.6921,0.4);
c16_3->SetFillColor(0);
c16_3->SetBorderMode(0);
c16_3->SetBorderSize(2);
c16_3->SetGridx();
c16_3->SetGridy();
c16_3->SetLeftMargin(0);
c16_3->SetRightMargin(0);
c16_3->SetTopMargin(0);
c16_3->SetFrameBorderMode(0);
c16_3->SetFrameBorderMode(0);
multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("0.31<x_{B}<0.39 ");
Double_t Graph_fx1140[9] = {
20.22962,
57.37981,
95.86069,
135.3486,
182.9227,
229.5531,
263.4275,
301.2495,
341.6225};
Double_t Graph_fy1140[9] = {
0.09662004,
0.06965272,
0.269744,
-0.170935,
-0.1509905,
0.6309487,
0.06270791,
-0.2582676,
-0.152689};
Double_t Graph_fex1140[9] = {
0,
0,
0,
0,
0,
0,
0,
0,
0};
Double_t Graph_fey1140[9] = {
0.05297617,
0.09107484,
0.121636,
0.2287053,
0.3066164,
0.3879093,
0.1400804,
0.07623737,
0.06739512};
gre = new TGraphErrors(9,Graph_fx1140,Graph_fy1140,Graph_fex1140,Graph_fey1140);
gre->SetName("Graph");
gre->SetTitle("Graph");
gre->SetFillColor(1);
ci = TColor::GetColor("#00ff00");
gre->SetMarkerColor(ci);
gre->SetMarkerStyle(21);
TH1F *Graph_Graph1140 = new TH1F("Graph_Graph1140","Graph",100,0,373.7618);
Graph_Graph1140->SetMinimum(-0.6052533);
Graph_Graph1140->SetMaximum(1.166504);
Graph_Graph1140->SetDirectory(0);
Graph_Graph1140->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1140->SetLineColor(ci);
Graph_Graph1140->GetXaxis()->SetLabelFont(42);
Graph_Graph1140->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1140->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1140->GetXaxis()->SetTitleFont(42);
Graph_Graph1140->GetYaxis()->SetLabelFont(42);
Graph_Graph1140->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1140->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1140->GetYaxis()->SetTitleFont(42);
Graph_Graph1140->GetZaxis()->SetLabelFont(42);
Graph_Graph1140->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1140->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1140->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1140);
TF1 *myfit1141 = new TF1("myfit","[0]*sin(x*3.1416/180.0)/(1 + [1]*cos(x*3.1416/180.0))",0,373.7618);
myfit1141->SetFillColor(19);
myfit1141->SetFillStyle(0);
myfit1141->SetLineColor(2);
myfit1141->SetLineWidth(2);
myfit1141->SetChisquare(10.84482);
myfit1141->SetNDF(7);
myfit1141->GetXaxis()->SetLabelFont(42);
myfit1141->GetXaxis()->SetLabelSize(0.035);
myfit1141->GetXaxis()->SetTitleSize(0.035);
myfit1141->GetXaxis()->SetTitleFont(42);
myfit1141->GetYaxis()->SetLabelFont(42);
myfit1141->GetYaxis()->SetLabelSize(0.035);
myfit1141->GetYaxis()->SetTitleSize(0.035);
myfit1141->GetYaxis()->SetTitleFont(42);
myfit1141->SetParameter(0,0.124049);
myfit1141->SetParError(0,0.05069177);
myfit1141->SetParLimits(0,0,0);
myfit1141->SetParameter(1,-0.6977989);
myfit1141->SetParError(1,0.2049165);
myfit1141->SetParLimits(1,0,0);
gre->GetListOfFunctions()->Add(myfit1141);
ptstats = new TPaveStats(0.27,0.74,0.9,0.9,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
AText = ptstats->AddText("#chi^{2} / ndf = 10.84 / 7");
AText = ptstats->AddText("p0 = 0.124 #pm 0.05069 ");
AText = ptstats->AddText("p1 = -0.6978 #pm 0.2049 ");
ptstats->SetOptStat(0);
ptstats->SetOptFit(111);
ptstats->Draw();
gre->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(gre->GetListOfFunctions());
multigraph->Add(gre,"");
multigraph->Draw("AP");
multigraph->GetXaxis()->SetTitle("#phi [deg.]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
gPad->Modified();
multigraph->GetYaxis()->SetRangeUser(-0.35, 0.35);
c16_3->Modified();
c16->cd();
// ------------>Primitives in pad: c16_4
TPad *c16_4 = new TPad("c16_4", "Incoherent channe_4",0.7388081,0,0.9767442,0.9);
c16_4->Draw();
c16_4->cd();
c16_4->Range(-17.15326,-0.4888889,360.2184,0.4);
c16_4->SetFillColor(0);
c16_4->SetBorderMode(0);
c16_4->SetBorderSize(2);
c16_4->SetGridx();
c16_4->SetGridy();
c16_4->SetLeftMargin(0);
c16_4->SetRightMargin(0);
c16_4->SetTopMargin(0);
c16_4->SetFrameBorderMode(0);
c16_4->SetFrameBorderMode(0);
multigraph = new TMultiGraph();
multigraph->SetName("");
multigraph->SetTitle("0.39<x_{B}<0.60 ");
Double_t Graph_fx1143[9] = {
19.13311,
55.36321,
90.36243,
127.5684,
0,
232.109,
268.7816,
302.6639,
343.0651};
Double_t Graph_fy1143[9] = {
0.1307559,
0.1996612,
0.006608881,
-2.830736,
0,
-2.034047,
0.1276631,
-0.2950426,
-0.003784223};
Double_t Graph_fex1143[9] = {
0,
0,
0,
0,
0,
0,
0,
0,
0};
Double_t Graph_fey1143[9] = {
0.05218982,
0.09871893,
0.1820512,
4.457049,
0,
1.502801,
0.207584,
0.0825667,
0.06276702};
gre = new TGraphErrors(9,Graph_fx1143,Graph_fy1143,Graph_fex1143,Graph_fey1143);
gre->SetName("Graph");
gre->SetTitle("Graph");
gre->SetFillColor(1);
ci = TColor::GetColor("#00ff00");
gre->SetMarkerColor(ci);
gre->SetMarkerStyle(21);
TH1F *Graph_Graph1143 = new TH1F("Graph_Graph1143","Graph",100,0,377.3716);
Graph_Graph1143->SetMinimum(-8.179195);
Graph_Graph1143->SetMaximum(2.517722);
Graph_Graph1143->SetDirectory(0);
Graph_Graph1143->SetStats(0);
ci = TColor::GetColor("#000099");
Graph_Graph1143->SetLineColor(ci);
Graph_Graph1143->GetXaxis()->SetLabelFont(42);
Graph_Graph1143->GetXaxis()->SetLabelSize(0.035);
Graph_Graph1143->GetXaxis()->SetTitleSize(0.035);
Graph_Graph1143->GetXaxis()->SetTitleFont(42);
Graph_Graph1143->GetYaxis()->SetLabelFont(42);
Graph_Graph1143->GetYaxis()->SetLabelSize(0.035);
Graph_Graph1143->GetYaxis()->SetTitleSize(0.035);
Graph_Graph1143->GetYaxis()->SetTitleFont(42);
Graph_Graph1143->GetZaxis()->SetLabelFont(42);
Graph_Graph1143->GetZaxis()->SetLabelSize(0.035);
Graph_Graph1143->GetZaxis()->SetTitleSize(0.035);
Graph_Graph1143->GetZaxis()->SetTitleFont(42);
gre->SetHistogram(Graph_Graph1143);
TF1 *myfit1144 = new TF1("myfit","[0]*sin(x*3.1416/180.0)/(1 + [1]*cos(x*3.1416/180.0))",0,377.3716);
myfit1144->SetFillColor(19);
myfit1144->SetFillStyle(0);
myfit1144->SetLineColor(2);
myfit1144->SetLineWidth(2);
myfit1144->SetChisquare(8.6938);
myfit1144->SetNDF(7);
myfit1144->GetXaxis()->SetLabelFont(42);
myfit1144->GetXaxis()->SetLabelSize(0.035);
myfit1144->GetXaxis()->SetTitleSize(0.035);
myfit1144->GetXaxis()->SetTitleFont(42);
myfit1144->GetYaxis()->SetLabelFont(42);
myfit1144->GetYaxis()->SetLabelSize(0.035);
myfit1144->GetYaxis()->SetTitleSize(0.035);
myfit1144->GetYaxis()->SetTitleFont(42);
myfit1144->SetParameter(0,0.1663497);
myfit1144->SetParError(0,0.06722025);
myfit1144->SetParLimits(0,0,0);
myfit1144->SetParameter(1,-0.5093054);
myfit1144->SetParError(1,0.2952021);
myfit1144->SetParLimits(1,0,0);
gre->GetListOfFunctions()->Add(myfit1144);
ptstats = new TPaveStats(0.27,0.74,0.9,0.9,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
AText = ptstats->AddText("#chi^{2} / ndf = 8.694 / 7");
AText = ptstats->AddText("p0 = 0.1663 #pm 0.06722 ");
AText = ptstats->AddText("p1 = -0.5093 #pm 0.2952 ");
ptstats->SetOptStat(0);
ptstats->SetOptFit(111);
ptstats->Draw();
gre->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(gre->GetListOfFunctions());
multigraph->Add(gre,"");
multigraph->Draw("AP");
multigraph->GetXaxis()->SetTitle("#phi [deg.]");
multigraph->GetXaxis()->SetLabelFont(42);
multigraph->GetXaxis()->SetLabelSize(0.035);
multigraph->GetXaxis()->SetTitleSize(0.035);
multigraph->GetXaxis()->SetTitleFont(42);
multigraph->GetYaxis()->SetLabelFont(42);
multigraph->GetYaxis()->SetLabelSize(0.035);
multigraph->GetYaxis()->SetTitleSize(0.035);
multigraph->GetYaxis()->SetTitleFont(42);
gPad->Modified();
multigraph->GetYaxis()->SetRangeUser(-0.35, 0.35);
c16_4->Modified();
c16->cd();
c16->Modified();
c16->cd();
c16->SetSelected(c16);
}
void yz-plane-length()
{
//=========Macro generated from canvas: c/
//========= (Tue Jul 28 19:55:12 2015) by ROOT version6.04/00
TCanvas *c = new TCanvas("c", "",0,45,600,500);
c->SetHighLightColor(2);
c->Range(-187.5,-2928.319,187.5,26354.87);
c->SetFillColor(0);
c->SetBorderMode(0);
c->SetBorderSize(2);
c->SetFrameBorderMode(0);
c->SetFrameBorderMode(0);
TH1D *hist__1 = new TH1D("hist__1","Worry Later",100,-150,150);
hist__1->SetBinContent(0,14);
hist__1->SetBinContent(1,1);
hist__1->SetBinContent(3,1);
hist__1->SetBinContent(4,2);
hist__1->SetBinContent(5,4);
hist__1->SetBinContent(6,4);
hist__1->SetBinContent(7,9);
hist__1->SetBinContent(8,20);
hist__1->SetBinContent(9,20);
hist__1->SetBinContent(10,19);
hist__1->SetBinContent(11,21);
hist__1->SetBinContent(12,41);
hist__1->SetBinContent(13,42);
hist__1->SetBinContent(14,72);
hist__1->SetBinContent(15,82);
hist__1->SetBinContent(16,124);
hist__1->SetBinContent(17,144);
hist__1->SetBinContent(18,204);
hist__1->SetBinContent(19,254);
hist__1->SetBinContent(20,298);
hist__1->SetBinContent(21,353);
hist__1->SetBinContent(22,426);
hist__1->SetBinContent(23,534);
hist__1->SetBinContent(24,636);
hist__1->SetBinContent(25,718);
hist__1->SetBinContent(26,800);
hist__1->SetBinContent(27,817);
hist__1->SetBinContent(28,856);
hist__1->SetBinContent(29,901);
hist__1->SetBinContent(30,928);
hist__1->SetBinContent(31,930);
hist__1->SetBinContent(32,1032);
hist__1->SetBinContent(33,1025);
hist__1->SetBinContent(34,1121);
hist__1->SetBinContent(35,1070);
hist__1->SetBinContent(36,1200);
hist__1->SetBinContent(37,1355);
hist__1->SetBinContent(38,1453);
hist__1->SetBinContent(39,1633);
hist__1->SetBinContent(40,1926);
hist__1->SetBinContent(41,2337);
hist__1->SetBinContent(42,2676);
hist__1->SetBinContent(43,3487);
hist__1->SetBinContent(44,4580);
hist__1->SetBinContent(45,5833);
hist__1->SetBinContent(46,7742);
hist__1->SetBinContent(47,10511);
hist__1->SetBinContent(48,13797);
hist__1->SetBinContent(49,17765);
hist__1->SetBinContent(50,22311);
hist__1->SetBinContent(51,22160);
hist__1->SetBinContent(52,17962);
hist__1->SetBinContent(53,14130);
hist__1->SetBinContent(54,10947);
hist__1->SetBinContent(55,8303);
hist__1->SetBinContent(56,6081);
hist__1->SetBinContent(57,4643);
hist__1->SetBinContent(58,3602);
hist__1->SetBinContent(59,2825);
hist__1->SetBinContent(60,2235);
hist__1->SetBinContent(61,1955);
hist__1->SetBinContent(62,1716);
hist__1->SetBinContent(63,1461);
hist__1->SetBinContent(64,1266);
hist__1->SetBinContent(65,1191);
hist__1->SetBinContent(66,1098);
hist__1->SetBinContent(67,1054);
hist__1->SetBinContent(68,999);
hist__1->SetBinContent(69,957);
hist__1->SetBinContent(70,982);
hist__1->SetBinContent(71,980);
hist__1->SetBinContent(72,932);
hist__1->SetBinContent(73,910);
hist__1->SetBinContent(74,873);
hist__1->SetBinContent(75,812);
hist__1->SetBinContent(76,785);
hist__1->SetBinContent(77,748);
hist__1->SetBinContent(78,649);
hist__1->SetBinContent(79,593);
hist__1->SetBinContent(80,459);
hist__1->SetBinContent(81,400);
hist__1->SetBinContent(82,321);
hist__1->SetBinContent(83,244);
hist__1->SetBinContent(84,198);
hist__1->SetBinContent(85,180);
hist__1->SetBinContent(86,147);
hist__1->SetBinContent(87,93);
hist__1->SetBinContent(88,73);
hist__1->SetBinContent(89,50);
hist__1->SetBinContent(90,43);
hist__1->SetBinContent(91,37);
hist__1->SetBinContent(92,22);
hist__1->SetBinContent(93,15);
hist__1->SetBinContent(94,10);
hist__1->SetBinContent(95,7);
hist__1->SetBinContent(96,7);
hist__1->SetBinContent(97,3);
hist__1->SetBinContent(98,2);
hist__1->SetBinContent(99,3);
hist__1->SetBinContent(100,3);
hist__1->SetBinContent(101,28);
hist__1->SetEntries(227323);
TPaveStats *ptstats = new TPaveStats(0.78,0.775,0.98,0.935,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
TText *AText = ptstats->AddText("hist");
AText->SetTextSize(0.0368);
AText = ptstats->AddText("Entries = 227323 ");
AText = ptstats->AddText("Mean = 0.5797");
AText = ptstats->AddText("Std Dev = 28.66");
ptstats->SetOptStat(1111);
ptstats->SetOptFit(0);
ptstats->Draw();
hist__1->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(hist__1);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#000099");
hist__1->SetLineColor(ci);
hist__1->GetXaxis()->SetTitle("#DeltaZ [cm]");
hist__1->GetXaxis()->SetLabelFont(42);
hist__1->GetXaxis()->SetLabelSize(0.035);
hist__1->GetXaxis()->SetTitleSize(0.035);
hist__1->GetXaxis()->SetTitleFont(42);
hist__1->GetYaxis()->SetLabelFont(42);
hist__1->GetYaxis()->SetLabelSize(0.035);
hist__1->GetYaxis()->SetTitleSize(0.035);
hist__1->GetYaxis()->SetTitleFont(42);
hist__1->GetZaxis()->SetLabelFont(42);
hist__1->GetZaxis()->SetLabelSize(0.035);
hist__1->GetZaxis()->SetTitleSize(0.035);
hist__1->GetZaxis()->SetTitleFont(42);
hist__1->Draw("");
TPaveText *pt = new TPaveText(0.3779599,0.9368947,0.6220401,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
AText = pt->AddText("Worry Later");
pt->Draw();
c->Modified();
c->cd();
c->SetSelected(c);
}
void DrawAllHistos(TString filename,Bool_t isMC,TString format){
// setTDRStyle("logredblue");
std::vector<TH1F*> h1vec;
std::vector<TH2F*> h2vec;
h1vec.clear();
h2vec.clear();
TFile *file = TFile::Open(filename);
// to get the names of the conditions
TObjArray *conditions_from_name = filename.Tokenize("_");
TString sa = conditions_from_name->At(1)->GetName();
TString sb = conditions_from_name->At(3)->GetName();
sb.ReplaceAll(".root","");
file->cd();
TDirectory *stardir = gDirectory;
TObject *thesourcedir;
TIter nextdir(stardir->GetListOfKeys());
while((thesourcedir=nextdir())){
TString dirName = thesourcedir->GetName();
stardir->cd(dirName);
TDirectory *current_sourcedir = gDirectory;
TH1::AddDirectory(kFALSE);
std::cout << "*************************" <<std::endl;
std::cout << "Reading Directory: " << dirName <<std::endl;
TObject *obj;
TIter next(current_sourcedir->GetListOfKeys());
while ((obj=next())) {
TString objName =obj->GetName();
if(objName.Contains("pfx")) continue;
if (objName.Contains("h2") && !objName.Contains("pfx")) {
//std::cout << "Reading: " << obj->GetName() <<std::endl;
TH2F* h2 = (TH2F*)file->Get(dirName+"/"+objName);
h2->SetName(objName+dirName);
h2vec.push_back(h2);
TCanvas *c = new TCanvas(h2->GetName(),h2->GetName(),800,600);
c->cd();
gPad->SetTopMargin(0.08);
gPad->SetRightMargin(0.15);
h2->SetStats(kFALSE);
h2->Draw("colz");
c->Draw();
c->cd();
TProfile *hpfx_tmp = (TProfile*) h2->ProfileX("_pfx",1,-1,"o");
hpfx_tmp->SetStats(kFALSE);
//hpfx_tmp->SetMarkerColor(kBlack);
hpfx_tmp->SetMarkerColor(kRed);
hpfx_tmp->SetMarkerSize(1.2);
hpfx_tmp->SetMarkerStyle(20);
hpfx_tmp->Draw("psame");
c->Draw();
cmsPrel(60.,sa,sb, isMC);
TString canvName = h2->GetName();
c->cd()->SetLogz();
c->SaveAs(canvName+"."+format);
} else {
TH1F* h1 = (TH1F*)file->Get(dirName+"/"+objName);
h1->SetName(objName+dirName);
h1vec.push_back(h1);
TCanvas *c = new TCanvas(h1->GetName(),h1->GetName(),600,600);
c->cd()->SetLogy();
h1->SetMarkerColor(kBlack);
h1->SetMarkerStyle(20);
h1->SetLineWidth(1.5);
h1->SetFillColor(393);
//h1->SetFillStyle(3005);
h1->Draw("hist");
h1->Draw("e1same");
c->Draw();
TObject *statObj;
TPaveStats *stats;
statObj = h1->GetListOfFunctions()->FindObject("stats");
stats= static_cast<TPaveStats*>(statObj);
stats->SetFillColor(10);
stats->SetLineWidth(1);
stats->SetShadowColor(0);
stats->SetTextFont(42);
stats->SetTextSize(0.025);
//stats->SetLineColor(LineColors);
//stats->SetTextColor(LineColors);
stats->SetX1NDC(0.75);
stats->SetY1NDC(0.72);
stats->SetX2NDC(0.97);
stats->SetY2NDC(0.92);
stats->Draw("same");
cmsPrel(60.,sa,sb,isMC);
TString canvName = h1->GetName();
c->SaveAs(canvName+"."+format);
}
}
}
}
void JetEtaBC(void){
std::cout << "Hellowww!!" << std::endl;
TFile *file1;
file1 = TFile::Open("PF_second.root");
TH1F* f1 = (TH1F*)file1->FindObjectAny("JetEtaBC");
TFile *file2;
file2 = TFile::Open("PAT_second.root");
TH1F* f2 = (TH1F*)file2->FindObjectAny("JetEtaBC");
TFile *file3;
file3 = TFile::Open("PF_IC5_Jets_Corrected.root");
TH1F* f3 = (TH1F*)file3->FindObjectAny("JetEtaBC");
TFile *file4;
file4 = TFile::Open("PF_ak5_tot.root");
TH1F* f4 = (TH1F*)file4->FindObjectAny("JetEtaBC");
TCanvas *c1 = new TCanvas("canvas","Plotting Canvas",150,10,990,660);
c1->Range(-0.1486535,-622.7141,1.350808,5500.641);
c1->SetFillColor(0);
c1->SetBorderSize(2);
c1->SetLogy();
c1->SetTickx(1);
c1->SetTicky(1);
f1->SetStats(kFALSE);
f4->GetXaxis()->SetTitle("#eta");
f4->GetXaxis()->SetLabelFont(42);
f4->GetXaxis()->SetTitleFont(42);
f4->GetYaxis()->SetTitle("Events");
f4->GetYaxis()->SetLabelFont(42);
f4->GetYaxis()->SetTitleOffset(1.19);
f4->GetYaxis()->SetTitleFont(42);
f4->SetLineWidth(4);
f1->SetLineColor(1);
f2->SetStats(kFALSE);
f2->SetLineColor(2);
f2->SetLineWidth(3);
f3->SetStats(kFALSE);
f3->SetLineColor(3);
f3->SetLineWidth(2);
f4->SetStats(kFALSE);
f4->SetLineColor(4);
f1->SetLineWidth(1);
int entries = f1->GetEntries();
double mean = f1->GetMean();
double RMS = f1->GetRMS();
int entries1 = f2->GetEntries();
double mean1 = f2->GetMean();
double RMS1 = f2->GetRMS();
int entries2 = f3->GetEntries();
double mean2 = f3->GetMean();
double RMS2 = f3->GetRMS();
int entries3 = f4->GetEntries();
double mean3 = f4->GetMean();
double RMS3 = f4->GetRMS();
TPaveStats *ptstats = new TPaveStats(0.8314943,0.7488136,0.9926437,0.8877119,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
ptstats->SetTextColor(1);
TPaveStats *ptstats1 = new TPaveStats(0.8314943,0.5988136,0.9926437,0.7377119,"brNDC");
ptstats1->SetName("stats");
ptstats1->SetBorderSize(1);
ptstats1->SetFillColor(0);
ptstats1->SetTextAlign(12);
ptstats1->SetTextFont(42);
ptstats1->SetTextColor(2);
TPaveStats *ptstats2 = new TPaveStats(0.8314943,0.4488136,0.9926437,0.5877119,"brNDC");
ptstats2->SetName("stats");
ptstats2->SetBorderSize(1);
ptstats2->SetFillColor(0);
ptstats2->SetTextAlign(12);
ptstats2->SetTextFont(42);
ptstats2->SetTextColor(3);
TPaveStats *ptstats3 = new TPaveStats(0.8314943,0.2988136,0.9926437,0.4377119,"brNDC");
ptstats3->SetName("stats");
ptstats3->SetBorderSize(1);
ptstats3->SetFillColor(0);
ptstats3->SetTextAlign(12);
ptstats3->SetTextFont(42);
ptstats3->SetTextColor(4);
#include <sstream>
std::string s,s1,s2;
std::stringstream out,out1,out2;
out << "Entries = " << entries << " ";
out1 << "Mean = " << mean << " ";
out2 << "RMS = " << RMS << " ";
s = out.str();
s1 = out1.str();
s2 = out2.str();
std::string s_1,s1_1,s2_1;
std::stringstream out_1,out1_1,out2_1;
out_1 << "Entries = " << entries1 << " ";
out1_1 << "Mean = " << mean1 << " ";
out2_1 << "RMS = " << RMS1 << " ";
s_1 = out_1.str();
s1_1 = out1_1.str();
s2_1 = out2_1.str();
std::string s_2,s1_2,s2_2;
std::stringstream out_2,out1_2,out2_2;
out_2 << "Entries = " << entries2 << " ";
out1_2 << "Mean = " << mean2 << " ";
out2_2 << "RMS = " << RMS2 << " ";
s_2 = out_2.str();
s1_2 = out1_2.str();
s2_2 = out2_2.str();
std::string s_3,s1_3,s2_3;
std::stringstream out_3,out1_3,out2_3;
out_3 << "Entries = " << entries3 << " ";
out1_3 << "Mean = " << mean3 << " ";
out2_3 << "RMS = " << RMS3 << " ";
s_3 = out_3.str();
s1_3 = out1_3.str();
s2_3 = out2_3.str();
TText *text = ptstats->AddText(s.c_str());
text = ptstats->AddText(s1.c_str());
text = ptstats->AddText(s2.c_str());
ptstats->SetOptStat(1110);
ptstats->SetOptFit(0);
TText *text1 = ptstats1->AddText(s_1.c_str());
text1 = ptstats1->AddText(s1_1.c_str());
text1 = ptstats1->AddText(s2_1.c_str());
ptstats1->SetOptStat(1110);
ptstats1->SetOptFit(0);
TText *text2 = ptstats2->AddText(s_2.c_str());
text2 = ptstats2->AddText(s1_2.c_str());
text2 = ptstats2->AddText(s2_2.c_str());
ptstats2->SetOptStat(1110);
ptstats2->SetOptFit(0);
TText *text3 = ptstats3->AddText(s_3.c_str());
text3 = ptstats3->AddText(s1_3.c_str());
text3 = ptstats3->AddText(s2_3.c_str());
ptstats3->SetOptStat(1110);
ptstats3->SetOptFit(0);
TPaveText *pt = new TPaveText(0.3066092,0.8061864,0.7255747,0.87055085,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetTextFont(42);
pt->SetFillColor(0);
pt->SetTextColor(1);
text = pt->AddText("PF Number of Jets before cuts");
text = pt->AddText("PF Jets Pt > 10 GeV Before Cuts");
TPaveText *pt1 = new TPaveText(0.3066092,0.7261864,0.7055747,0.77055085,"blNDC");
pt1->SetName("title");
pt1->SetBorderSize(0);
pt1->SetTextFont(42);
pt1->SetFillColor(0);
pt1->SetTextColor(2);
text1 = pt1->AddText("PAT Jets Eta Before Cuts");
TPaveText *pt2 = new TPaveText(0.3066092,0.6561864,0.7055747,0.70055085,"blNDC");
pt2->SetName("title");
pt2->SetBorderSize(0);
pt2->SetTextFont(42);
pt2->SetFillColor(0);
pt2->SetTextColor(3);
text2 = pt2->AddText("PF IC5 Pt > 0 GeV Jets Eta Before Cuts");
TPaveText *pt3 = new TPaveText(0.3266092,0.5861864,0.7055747,0.63055085,"blNDC");
pt3->SetName("title");
pt3->SetBorderSize(0);
pt3->SetTextFont(42);
pt3->SetFillColor(0);
pt3->SetTextColor(4);
text3 = pt3->AddText("PF ak5 Eta Before Cuts");
c1->Modified();
f4->Draw();
f2->Draw("SAME");
f3->Draw("SAME");
f1->Draw("SAME");
pt->Draw("SAME");
ptstats->Draw("SAME");
pt1->Draw("SAME");
ptstats1->Draw("SAME");
pt2->Draw("SAME");
ptstats2->Draw("SAME");
pt3->Draw("SAME");
ptstats3->Draw("SAME");
c1->SaveAs("JetEtaBeforeCutsComparisonAllIhaverightnow.gif");
return 0;
};
int main(int argc, char** argv){
//=========Macro generated from canvas: Canvas_1/Canvas_1
//========= (Wed Aug 27 23:28:50 2014) by ROOT version5.34/04
TCanvas *Canvas_1 = new TCanvas("Canvas_1", "Canvas_1",258,73,1243,811);
Canvas_1->Range(-0.0625,-0.0625,0.0625,0.0625);
Canvas_1->SetFillColor(0);
Canvas_1->SetBorderMode(0);
Canvas_1->SetBorderSize(2);
Canvas_1->SetFrameBorderMode(0);
Canvas_1->SetFrameBorderMode(0);
TFile f(argv[1],"r");
TH2F *h1 = (TH2F*)f.Get("h2dir/h4");
/*
TPaletteAxis *palette = new TPaletteAxis(0.050625,-0.05,0.05625,0.05,h1);
palette->SetLabelColor(1);
palette->SetLabelFont(42);
palette->SetLabelOffset(0.005);
palette->SetLabelSize(0.035);
palette->SetTitleOffset(1);
palette->SetTitleSize(0.035);
palette->SetFillColor(100);
palette->SetFillStyle(1001);
h1->GetListOfFunctions()->Add(palette,"br");
*/
TPaveStats *ptstats = new TPaveStats(0.700565,0.6608696,0.9007264,0.9006211,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
TText *text = ptstats->AddText("h1");
text->SetTextSize(0.03676191);
text = ptstats->AddText("Entries = 101 ");
text = ptstats->AddText("Mean x = -0.0003512");
text = ptstats->AddText("Mean y = -0.0003318");
text = ptstats->AddText("RMS x = 0.01432");
text = ptstats->AddText("RMS y = 0.01428");
ptstats->SetOptStat(1111);
ptstats->SetOptFit(0);
ptstats->Draw();
h1->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(h1);
Int_t ci; // for color index setting
ci = TColor::GetColor("#000099");
h1->SetLineColor(ci);
h1->GetXaxis()->SetTitle("mm");
h1->GetXaxis()->SetLabelFont(42);
h1->GetXaxis()->SetLabelSize(0.035);
h1->GetXaxis()->SetTitleSize(0.035);
h1->GetXaxis()->SetTitleFont(42);
h1->GetYaxis()->SetTitle("mm");
h1->GetYaxis()->SetLabelFont(42);
h1->GetYaxis()->SetLabelSize(0.035);
h1->GetYaxis()->SetTitleSize(0.035);
h1->GetYaxis()->SetTitleFont(42);
h1->GetZaxis()->SetLabelFont(42);
h1->GetZaxis()->SetLabelSize(0.035);
h1->GetZaxis()->SetTitleSize(0.035);
h1->GetZaxis()->SetTitleFont(42);
h1->Draw("COLZ ");
TPaveText *pt = new TPaveText(0.1594915,0.9368634,0.8405085,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
text = pt->AddText("2D Energy deposit distribution at z=-25um");
pt->Draw();
Canvas_1->Modified();
Canvas_1->cd();
Canvas_1->SetSelected(Canvas_1);
std::cout<<argv[1]<<std::endl;
Canvas_1->SaveAs(TString(argv[1])+"35#.png");
}
void initial_z_diff()
{
//=========Macro generated from canvas: c/
//========= (Fri Jul 24 10:47:18 2015) by ROOT version6.04/00
TCanvas *c = new TCanvas("c", "",0,45,600,500);
c->SetHighLightColor(2);
c->Range(-187.5,-3568.557,187.5,32117.01);
c->SetFillColor(0);
c->SetBorderMode(0);
c->SetBorderSize(2);
c->SetFrameBorderMode(0);
c->SetFrameBorderMode(0);
TH1D *hist__1 = new TH1D("hist__1","Difference in PMT-constructed position and mc-truth start position",100,-150,150);
hist__1->SetBinContent(0,49);
hist__1->SetBinContent(1,12);
hist__1->SetBinContent(2,7);
hist__1->SetBinContent(3,9);
hist__1->SetBinContent(4,12);
hist__1->SetBinContent(5,21);
hist__1->SetBinContent(6,30);
hist__1->SetBinContent(7,36);
hist__1->SetBinContent(8,47);
hist__1->SetBinContent(9,67);
hist__1->SetBinContent(10,70);
hist__1->SetBinContent(11,90);
hist__1->SetBinContent(12,140);
hist__1->SetBinContent(13,162);
hist__1->SetBinContent(14,208);
hist__1->SetBinContent(15,208);
hist__1->SetBinContent(16,317);
hist__1->SetBinContent(17,367);
hist__1->SetBinContent(18,460);
hist__1->SetBinContent(19,581);
hist__1->SetBinContent(20,685);
hist__1->SetBinContent(21,838);
hist__1->SetBinContent(22,963);
hist__1->SetBinContent(23,1079);
hist__1->SetBinContent(24,1231);
hist__1->SetBinContent(25,1281);
hist__1->SetBinContent(26,1264);
hist__1->SetBinContent(27,1318);
hist__1->SetBinContent(28,1318);
hist__1->SetBinContent(29,1265);
hist__1->SetBinContent(30,1282);
hist__1->SetBinContent(31,1235);
hist__1->SetBinContent(32,1315);
hist__1->SetBinContent(33,1290);
hist__1->SetBinContent(34,1399);
hist__1->SetBinContent(35,1363);
hist__1->SetBinContent(36,1528);
hist__1->SetBinContent(37,1700);
hist__1->SetBinContent(38,1822);
hist__1->SetBinContent(39,2053);
hist__1->SetBinContent(40,2441);
hist__1->SetBinContent(41,2903);
hist__1->SetBinContent(42,3359);
hist__1->SetBinContent(43,4326);
hist__1->SetBinContent(44,5618);
hist__1->SetBinContent(45,7284);
hist__1->SetBinContent(46,9575);
hist__1->SetBinContent(47,12857);
hist__1->SetBinContent(48,16804);
hist__1->SetBinContent(49,21562);
hist__1->SetBinContent(50,27189);
hist__1->SetBinContent(51,26947);
hist__1->SetBinContent(52,21828);
hist__1->SetBinContent(53,17254);
hist__1->SetBinContent(54,13388);
hist__1->SetBinContent(55,10171);
hist__1->SetBinContent(56,7452);
hist__1->SetBinContent(57,5752);
hist__1->SetBinContent(58,4452);
hist__1->SetBinContent(59,3522);
hist__1->SetBinContent(60,2826);
hist__1->SetBinContent(61,2449);
hist__1->SetBinContent(62,2170);
hist__1->SetBinContent(63,1855);
hist__1->SetBinContent(64,1635);
hist__1->SetBinContent(65,1515);
hist__1->SetBinContent(66,1393);
hist__1->SetBinContent(67,1344);
hist__1->SetBinContent(68,1279);
hist__1->SetBinContent(69,1232);
hist__1->SetBinContent(70,1248);
hist__1->SetBinContent(71,1277);
hist__1->SetBinContent(72,1226);
hist__1->SetBinContent(73,1247);
hist__1->SetBinContent(74,1229);
hist__1->SetBinContent(75,1187);
hist__1->SetBinContent(76,1265);
hist__1->SetBinContent(77,1253);
hist__1->SetBinContent(78,1269);
hist__1->SetBinContent(79,1210);
hist__1->SetBinContent(80,1055);
hist__1->SetBinContent(81,914);
hist__1->SetBinContent(82,761);
hist__1->SetBinContent(83,632);
hist__1->SetBinContent(84,515);
hist__1->SetBinContent(85,450);
hist__1->SetBinContent(86,356);
hist__1->SetBinContent(87,264);
hist__1->SetBinContent(88,227);
hist__1->SetBinContent(89,176);
hist__1->SetBinContent(90,127);
hist__1->SetBinContent(91,134);
hist__1->SetBinContent(92,85);
hist__1->SetBinContent(93,60);
hist__1->SetBinContent(94,47);
hist__1->SetBinContent(95,29);
hist__1->SetBinContent(96,31);
hist__1->SetBinContent(97,18);
hist__1->SetBinContent(98,14);
hist__1->SetBinContent(99,12);
hist__1->SetBinContent(100,12);
hist__1->SetBinContent(101,63);
hist__1->SetEntries(289897);
TPaveStats *ptstats = new TPaveStats(0.78,0.775,0.98,0.935,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
TText *AText = ptstats->AddText("hist");
AText->SetTextSize(0.0368);
AText = ptstats->AddText("Entries = 289897 ");
AText = ptstats->AddText("Mean = 0.6888");
AText = ptstats->AddText("Std Dev = 33.12");
ptstats->SetOptStat(1111);
ptstats->SetOptFit(0);
ptstats->Draw();
hist__1->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(hist__1);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#000099");
hist__1->SetLineColor(ci);
hist__1->GetXaxis()->SetTitle("#DeltaZ [cm]");
hist__1->GetXaxis()->SetLabelFont(42);
hist__1->GetXaxis()->SetLabelSize(0.035);
hist__1->GetXaxis()->SetTitleSize(0.035);
hist__1->GetXaxis()->SetTitleFont(42);
hist__1->GetYaxis()->SetLabelFont(42);
hist__1->GetYaxis()->SetLabelSize(0.035);
hist__1->GetYaxis()->SetTitleSize(0.035);
hist__1->GetYaxis()->SetTitleFont(42);
hist__1->GetZaxis()->SetLabelFont(42);
hist__1->GetZaxis()->SetLabelSize(0.035);
hist__1->GetZaxis()->SetTitleSize(0.035);
hist__1->GetZaxis()->SetTitleFont(42);
hist__1->Draw("");
TPaveText *pt = new TPaveText(0.15,0.9368947,0.85,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
AText = pt->AddText("Difference in PMT-constructed position and mc-truth start position");
pt->Draw();
c->Modified();
c->cd();
c->SetSelected(c);
}
