DOMMatrix*
DOMMatrix::SetMatrixValue(const nsAString& aTransformList, ErrorResult& aRv)
{
SVGTransformListParser parser(aTransformList);
if (!parser.Parse()) {
aRv.Throw(NS_ERROR_DOM_SYNTAX_ERR);
} else {
mMatrix3D = nullptr;
mMatrix2D = new gfx::Matrix();
gfxMatrix result;
const nsTArray<nsSVGTransform>& mItems = parser.GetTransformList();
for (uint32_t i = 0; i < mItems.Length(); ++i) {
result.PreMultiply(mItems[i].GetMatrix());
}
SetA(result._11);
SetB(result._12);
SetC(result._21);
SetD(result._22);
SetE(result._31);
SetF(result._32);
}
return this;
}
DOMMatrixReadOnly*
DOMMatrixReadOnly::SetMatrixValue(const nsAString& aTransformList, ErrorResult& aRv)
{
// An empty string is a no-op.
if (aTransformList.IsEmpty()) {
return this;
}
gfx::Matrix4x4 transform;
bool contains3dTransform = false;
if (!ServoCSSParser::ParseTransformIntoMatrix(aTransformList,
contains3dTransform,
transform.components)) {
aRv.Throw(NS_ERROR_DOM_SYNTAX_ERR);
return nullptr;
}
if (!contains3dTransform) {
mMatrix3D = nullptr;
mMatrix2D = new gfx::Matrix();
SetA(transform._11);
SetB(transform._12);
SetC(transform._21);
SetD(transform._22);
SetE(transform._41);
SetF(transform._42);
} else {
mMatrix3D = new gfx::Matrix4x4(transform);
mMatrix2D = nullptr;
}
return this;
}
void LLBC_MD5::MD5GroupDigest::SetInfo(uint32 a, uint32 b, uint32 c, uint32 d)
{
SetA(a);
SetB(b);
SetC(c);
SetD(d);
}
O3QuadraticEquation<TYPE>& O3QuadraticEquation<TYPE>::Set(TYPE high_x_intercept, TYPE y_intercept) {
double x_sq = -( O3recip(high_x_intercept*high_x_intercept) ) * y_intercept;
double constant = y_intercept;
SetA(x_sq);
SetB(0);
SetC(constant);
return *this;
}
LLBC_MD5::MD5GroupDigest &LLBC_MD5::MD5GroupDigest::operator +=(const LLBC_MD5::MD5GroupDigest &right)
{
SetA(GetA() + right.GetA());
SetB(GetB() + right.GetB());
SetC(GetC() + right.GetC());
SetD(GetD() + right.GetD());
return *this;
}
BOOL Pipe::Read(PipePerisist *pPerisist)
{
ms_Name.SetValue(pPerisist->Name());
ms_StartPt.SetValue1(pPerisist->StartPtX());
ms_StartPt.SetValue2(pPerisist->StartPtY());
ms_EndPt.SetValue1(pPerisist->EndPtX());
ms_EndPt.SetValue2(pPerisist->EndPtY());
ms_InJun.SetValue(pPerisist->InJunID());
ms_OutJun.SetValue(pPerisist->OutJunID());
ms_InDia.SetValue(pPerisist->InDia());
ms_InDia.SetUnit(pPerisist->InDiaUnit());
ms_Len.SetValue(pPerisist->Len());
ms_Len.SetUnit(pPerisist->LenUnit());
ms_WaveSpeed.SetValue(pPerisist->WaveSpeed());
ms_WaveSpeed.SetUnit(pPerisist->WaveSpeedUnit());
ms_Material.SetValue(pPerisist->Material());
ms_Size.SetValue(pPerisist->Size());
ms_PipeType.SetValue(pPerisist->PipeType());
ms_Thick.SetValue(pPerisist->Thick());
ms_Thick.SetUnit(pPerisist->ThickUnit());
ms_Elasticity.SetValue(pPerisist->Elasticity());
ms_Elasticity.SetUnit(pPerisist->ElasticityUnit());
ms_PossionRatio.SetValue(pPerisist->PossionRatio());
ms_InDiaReduce.SetValue(pPerisist->InDiaReduce());
ms_CalcWave.SetValue(pPerisist->CalcWave());
ms_SupportType.SetValue(pPerisist->SupportType());
ms_C1.SetValue(pPerisist->C1());
ms_LossType.SetValue(pPerisist->FrictionModel());
ms_Friction.SetValue(pPerisist->Friction());
ms_Friction.SetUnit(pPerisist->FrictionUnit());
ms_LocalK.SetValue(pPerisist->LocalK());
ms_PipeFactor.SetValue(pPerisist->PipeFactor());
ms_FitFactor.SetValue(pPerisist->FitFactor());
ms_InitFlow.SetData(pPerisist->InitFlow());
ms_SectionNum.SetValue(pPerisist->SectionNum());
ms_ParaNum.SetValue(pPerisist->ParaNum());
ms_PartFull.SetValue(pPerisist->PartFull());
SetB(pPerisist->B());
SaveAllData();
return TRUE;
}
Triangle::Triangle(Point& a, Point& b, Point& c)
{
SetA(a);
SetB(b);
SetC(c);
}
void checkdatadj(TString datafilename, TString djfilename)
{
TFile *f = new TFile(datafilename);
auto nt = (TTree *)f->Get("nt");
auto hjtpt1 = new TH1F("data_jtpt1","data_jtpt1",ptbins,ptmin,ptmax);
nt->Project(hjtpt1->GetName(),"jtpt1","weight");
auto hjtpt2 = new TH1F("data_jtpt2","data_jtpt2",ptbins,0,200);
nt->Project(hjtpt2->GetName(),"jtpt2","weight*(dijet && jtpt1>120 && jtpt2>30 && acos(cos(jtphi1-jtphi2))>2./3*3.14)");
auto hdataaj = new TH1F("data_aj","data_aj",20,0,1);
nt->Project(hdataaj->GetName(),"abs(jtpt1-jtpt2)/(jtpt1+jtpt2)",
"weight*(dijet && jtpt1>120 && jtpt2>30 && acos(cos(jtphi1-jtphi2))>2./3*3.14)");
auto hdatabtgaj = new TH1F("data_btgaj","data_btgaj",20,0,1);
nt->Project(hdatabtgaj->GetName(),"abs(jtpt1-jtpt2)/(jtpt1+jtpt2)",
"weight*(dijet && jtpt1>120 && jtpt2>30 && acos(cos(jtphi1-jtphi2))>2./3*3.14 && discr_csvSimple1>0.9 && discr_csvSimple2>0.9)");
auto hdataxj = new TH1F("data_xj","data_xj",20,0,1);
nt->Project(hdataxj->GetName(),"jtpt2/jtpt1",
"weight*(dijet && jtpt1>120 && jtpt2>30 && acos(cos(jtphi1-jtphi2))>2./3*3.14)");
auto hdatadphi = new TH1F("data_dphi","data_dphi",20,0,3.142);
nt->Project(hdatadphi->GetName(),"acos(cos(jtphi1-jtphi2))",
"weight*(dijet && jtpt1>120 && jtpt2>30)");
auto hdatadphi2 = new TH1F("data_dphi2","data_dphi2",20,2,3.142);
nt->Project(hdatadphi2->GetName(),"acos(cos(jtphi1-jtphi2))",
"weight*(dijet && jtpt1>120 && jtpt2>30)");
auto hdatabtgdphi = new TH1F("data_btgdphi","data_btgdphi",20,0,3.142);
nt->Project(hdatabtgdphi->GetName(),"acos(cos(jtphi1-jtphi2))",
"weight*(dijet && jtpt1>120 && jtpt2>30 && discr_csvSimple1>0.9 && discr_csvSimple2>0.9)");
SetData({hjtpt1,hjtpt2,hdataaj,hdatabtgaj,hdataxj,hdatadphi,hdatadphi2,hdatabtgdphi});
SetB({hdatabtgaj,hdatabtgdphi});
fout->cd();
hjtpt1->Write();
hjtpt2->Write();
hdataaj->Write();
hdataxj->Write();
hdatadphi->Write();
hdatadphi2->Write();
hdatabtgaj->Write();
hdatabtgdphi->Write();
f->Close();
f = new TFile(djfilename);
nt = (TTree *)f->Get("nt");
auto hdjjtpt1 = new TH1F("dj_jtpt1","dj_jtpt1",ptbins,ptmin,ptmax);
nt->Project(hdjjtpt1->GetName(),"jtpt1","weight");
auto hmcjtpt2 = new TH1F("mc_jtpt2","mc_jtpt2",ptbins,0,200);
nt->Project(hmcjtpt2->GetName(),"jtpt2","weight*(dijet && jtpt1>120 && jtpt2>30 && acos(cos(jtphi1-jtphi2))>2./3*3.14)");
auto hmcaj = new TH1F("mc_aj","mc_aj",20,0,1);
nt->Project(hmcaj->GetName(),"abs(jtpt1-jtpt2)/(jtpt1+jtpt2)",
"weight*(dijet && jtpt1>120 && jtpt2>30 && acos(cos(jtphi1-jtphi2))>2./3*3.14)");
auto hmcbtgaj = new TH1F("mc_btgaj","mc_btgaj",20,0,1);
nt->Project(hmcbtgaj->GetName(),"abs(jtpt1-jtpt2)/(jtpt1+jtpt2)",
"weight*(dijet && jtpt1>120 && jtpt2>30 && acos(cos(jtphi1-jtphi2))>2./3*3.14 && discr_csvSimple1>0.9 && discr_csvSimple2>0.9)");
auto hmcxj = new TH1F("mc_xj","mc_xj",20,0,1);
nt->Project(hmcxj->GetName(),"jtpt2/jtpt1",
"weight*(dijet && jtpt1>120 && jtpt2>30 && acos(cos(jtphi1-jtphi2))>2./3*3.14)");
auto hmcdphi = new TH1F("mc_dphi","mc_dphi",20,0,3.142);
nt->Project(hmcdphi->GetName(),"acos(cos(jtphi1-jtphi2))",
"weight*(dijet && jtpt1>120 && jtpt2>30)");
auto hmcdphi2 = new TH1F("mc_dphi2","mc_dphi2",20,2,3.142);
nt->Project(hmcdphi2->GetName(),"acos(cos(jtphi1-jtphi2))",
"weight*(dijet && jtpt1>120 && jtpt2>30)");
auto hmcbtgdphi = new TH1F("mc_btgdphi","mc_btgdphi",20,0,3.142);
nt->Project(hmcbtgdphi->GetName(),"acos(cos(jtphi1-jtphi2))",
"weight*(dijet && jtpt1>120 && jtpt2>30 && discr_csvSimple1>0.9 && discr_csvSimple2>0.9)");
SetMC({hdjjtpt1,hmcjtpt2,hmcaj,hmcbtgaj,hmcxj,hmcdphi,hmcdphi2,hmcbtgdphi});
SetInc({hdjjtpt1,hmcjtpt2,hmcaj,hmcxj,hmcdphi,hmcdphi2});
SetB({hmcbtgaj,hmcbtgdphi});
auto temp = new TH1F("temp","temp",10,0.5,1.5);
nt->Project("temp","dijet","weight");
int numberofdijets = temp->Integral();
cout<<"# of dijets = "<<numberofdijets<<endl;
//hmcaj->Scale(1./numberofdijets);
//hmcxj->Scale(1./numberofdijets);
//hmcdphi->Scale(1./numberofdijets);
fout->cd();
hdjjtpt1->Write();
hmcjtpt2->Write();
hmcaj->Write();
hmcxj->Write();
hmcdphi->Write();
hmcdphi2->Write();
hmcbtgaj->Write();
hmcbtgdphi->Write();
f->Close();
}
void findtruthpp()
{
seth(10,0,1);
auto hmcppxJTrue = geth("hmcppxJTrue","true;x_{J};Event fractions");
auto hmcppxJTrueTag = geth("hmcppxJTrueTag","true tagged;x_{J};Event fractions");
auto hmcppxJTrueTagCorr = geth("hmcppxJTrueTagCorr","true tagged corrected;x_{J};Event fractions");
TFile *fmcpp = new TFile(config.getFileName_djt("mcppbfa"));
Fill(fmcpp,[&] (dict &m) {
if (m["pthat"]<pthatcut) return;
if (m[pairCodeSB1]!=0) return;
// if (m["pthat"]<80) return;
// if (m["bProdCode"]!=1) return;
float w = m["weight"]*processweight((int)m["bProdCode"]);
// float w = m["weight"];
// if (m["bProdCode"]==2) return;
float corr = tageffcorrectionpp(m["jtpt1"],m["jteta1"],m[jtptSB],m[jtetaSB]);
float wb = w*corr;
if (m["jtpt1"]>pt1cut && m["refpt1"]>50 && abs(m["refparton_flavorForB1"])==5 && m[jtptSB]>pt2cut && m[dphiSB1]>PI23)
hmcppxJTrue->Fill(m[jtptSB]/m["jtpt1"],w);
if (m["jtpt1"]>pt1cut && m["refpt1"]>50 && abs(m["refparton_flavorForB1"])==5 && m[jtptSB]>pt2cut && m[dphiSB1]>PI23
&& m["discr_csvV1_1"]>0.9 && m[discr_csvV1_SB]>0.9) {
hmcppxJTrueTag->Fill(m[jtptSB]/m["jtpt1"],w);
hmcppxJTrueTagCorr->Fill(m[jtptSB]/m["jtpt1"],wb);
}
});
NormalizeAllHists();
plotputmean = true;
plotytitle = "Event fractions";
plotdivide = false;
// aktstring += "R=0.4 |#eta|<2.0";
// plotsecondline = Form("p_{T,1}>%d GeV, p_{T,2}>%d GeV", (int)pt1cut, (int)pt2cut);
// plotthirdline = "#Delta#phi>2/3#pi";
plottextposy = 0.8;
plottextposx = 0.2;
plotmeanposy = 0.43;
plotymax = 0.2;
plotymin = 0;
plotlegendpos = BottomRight;//TopLeft;
aktstring = "pp";
SetMC({hmcppxJTrue,hmcppxJTrueTag,hmcppxJTrueTagCorr});
SetData({hmcppxJTrue});
Draw({hmcppxJTrue,hmcppxJTrueTag,hmcppxJTrueTagCorr});
hmcppxJTrue->SetMinimum(0);
hmcppxJTrue->SetMaximum(0.3);
hmcppxJTrue->SetLineWidth(2);
hmcppxJTrue->SetMarkerStyle(kNone);
hmcppxJTrue->SetFillStyle(0);
hmcppxJTrueTag->SetMarkerStyle(kOpenCircle);
hmcppxJTrueTagCorr->SetMarkerStyle(kOpenSquare);
plotymax = 0.3;
SetB({hmcppxJTrue,hmcppxJTrueTag,hmcppxJTrueTagCorr});
float xjtrue = hmcppxJTrue->GetMean();
float xjtruetag = hmcppxJTrueTag->GetMean();
float xjtruetagcorr = hmcppxJTrueTagCorr->GetMean();
float exjtrue = hmcppxJTrue->GetMeanError();
float exjtruetag = hmcppxJTrueTag->GetMeanError();
float exjtruetagcorr = hmcppxJTrueTagCorr->GetMeanError();
auto c = getc();
hmcppxJTrue->Draw("hist");
hmcppxJTrueTag->Draw("E1,same");
hmcppxJTrueTagCorr->Draw("E1,same");
plotlegendpos = TopLeft;
auto l = getLegend();
l->AddEntry(hmcppxJTrue,Form("b-dijets, #LTx_{J}#GT=%.3f#pm%.3f",xjtrue,exjtrue),"L");
l->AddEntry(hmcppxJTrueTag,Form("uncorrected, #LTx_{J}#GT=%.3f#pm%.3f",xjtruetag,exjtruetag),"P");
l->AddEntry(hmcppxJTrueTagCorr,Form("corrected, #LTx_{J}#GT=%.3f#pm%.3f",xjtruetagcorr,exjtruetagcorr),"P");
l->Draw();
TLatex *Tl = new TLatex();
Tl->DrawLatexNDC(0.2, 0.8, aktstring);
SavePlot(c,"closurepp");
}
void findtruthPbPb(int binMin, int binMax)
{
TFile *fmc = new TFile(config.getFileName_djt("mcPbbfa"));
buildNamesuffix = TString::Format("_bin_%d_%d",binMin, binMax);
// buildTitlesuffix = TString::Format("%d-%d %%",binMin/2, binMax/2);
seth(10,0,1);
auto hmcPbPbxJTrue = geth("hmcPbPbxJTrue","PbPb true;x_{J};Event fractions");
auto hmcPbPbxJTrueTag = geth("hmcPbPbxJTrueTag","PbPb true tagged;x_{J};Event fractions");
auto hmcPbPbxJTrueTagCorr = geth("hmcPbPbxJTrueTagCorr","PbPb true tagged corrected;x_{J};Event fractions");
auto hmcPbPbxJTrueTagCorrPt = geth("hmcPbPbxJTrueTagCorrPt","PbPb true tagged corrected pt;x_{J};Event fractions");
auto hmcPbPbxJTrueTagCorrEta = geth("hmcPbPbxJTrueTagCorrEta","PbPb true tagged corrected eta;x_{J};Event fractions");
auto hmcPbPbxJTrueTagCorrBin = geth("hmcPbPbxJTrueTagCorrBin","PbPb true tagged corrected bin;x_{J};Event fractions");
seth(12,20,140);//10,40,100);
auto hpt2true = geth("hpt2true","true;p_{T,2} GeV");
auto hpt2truetag = geth("hpt2truetag","true tagged;p_{T,2} GeV");
auto hpt2truetagovertrue = geth("hpt2truetagovertrue","true tagged/true;p_{T,2} GeV");
auto hpt2truetagcorr = geth("hpt2truetagcorr","true tagged corrected;p_{T,2} GeV");
auto hpt2truetagcorrovertrue = geth("hpt2truetagcorrovertrue","true tagged corrected/true;p_{T,2} GeV");
seth(10,100,200);
auto hpt1true = geth("hpt1true","true;p_{T,1} GeV");
auto hpt1truetag = geth("hpt1truetag","true tagged;p_{T,1} GeV");
auto hpt1truetagovertrue = geth("hpt1truetagovertrue","true tagged/true;p_{T,1} GeV");
auto hpt1truetagcorr = geth("hpt1truetagcorr","true tagged corrected;p_{T,1} GeV");
auto hpt1truetagcorrovertrue = geth("hpt1truetagcorrovertrue","true tagged corrected/true;p_{T,1} GeV");
seth(10,0,200);
auto hbintrue = geth("hbintrue","true;bin");
auto hbintruetag = geth("hbintruetag","true tagged;bin");
auto hbintruetagovertrue = geth("hbintruetagovertrue","true tagged/true;bin");
auto hbintruetagcorr = geth("hbintruetagcorr","true tagged corrected;bin");
auto hbintruetagcorrovertrue = geth("hbintruetagcorrovertrue","true tagged corrected/true;bin");
seth(20,-2,2);
auto heta2true = geth("heta2true","true;#eta_{2}");
auto heta2truetag = geth("heta2truetag","true tagged;#eta_{2}");
auto heta2truetagovertrue = geth("heta2truetagovertrue","true tagged/true;#eta_{2}");
auto heta2truetagcorr = geth("heta2truetagcorr","true tagged corrected;#eta_{2}");
auto heta2truetagcorrovertrue = geth("heta2truetagcorrovertrue","true tagged corrected/true;#eta_{2}");
auto heta1true = geth("heta1true","true;#eta_{1}");
auto heta1truetag = geth("heta1truetag","true tagged;#eta_{1}");
auto heta1truetagovertrue = geth("heta1truetagovertrue","true tagged/true;#eta_{1}");
auto heta1truetagcorr = geth("heta1truetagcorr","true tagged corrected;#eta_{1}");
auto heta1truetagcorrovertrue = geth("heta1truetagcorrovertrue","true tagged corrected/true;#eta_{1}");
unordered_set<int> eventstodiscard = {1805770,1116573,1084397};//,
// 5755734,1599758,395810,
// 1363321,211625,3195128};
//
Fill(fmc,[&] (dict &m) {
if (m["bin"]<binMin || m["bin"]>=binMax) return;
if (m["pthat"]<pthatcut) return;
if (m[pairCodeSB1]!=0) return;
// if (m["bProdCode"]>1) return;
if (eventstodiscard.find(m["event"])!=eventstodiscard.end()) return; //kill large-weight GSP event
float w = m["weight"]*processweight((int)m["bProdCode"]); //because we have only b-dijets
float corr = tageffcorrectionPbPb(m["jtpt1"],m["jteta1"],m[jtptSB],m[jtetaSB],m["bin"]);
// float corrpt = getPbPbcorrectionPt(m["jtpt1"],m[jtptSB]);
// float correta = getPbPbcorrectionEta(m["jteta1"],m[jtetaSB]);
// float corrbin = getPbPbcorrectionBin(m["bin"]);
float wb = w*corr;
if (m["jtpt1"]>pt1cut && m["refpt1"]>50 && abs(m["refparton_flavorForB1"])==5 && m[jtptSB]>pt2cut && m[refptSB]>20 && m[dphiSB1]>PI23) {
hmcPbPbxJTrue->Fill(m[jtptSB]/m["jtpt1"],w);
hpt2true->Fill(m[jtptSB],w);
hpt1true->Fill(m["jtpt1"],w);
heta2true->Fill(m[jtetaSB],w);
heta1true->Fill(m["jteta1"],w);
hbintrue->Fill(m["bin"],w);
}
if (m["jtpt1"]>pt1cut && m["refpt1"]>50 && abs(m["refparton_flavorForB1"])==5 && m[jtptSB]>pt2cut && m[refptSB]>20 && m[dphiSB1]>PI23
&& m["discr_csvV1_1"]>0.9 && m[discr_csvV1_SB]>0.9) { //
//corrpt *= m[jtptSB] < 60 ? 1./0.7 : 1;
//wb *= m[jtptSB] < 60 ? 1./0.7 : 1;
hmcPbPbxJTrueTag->Fill(m[jtptSB]/m["jtpt1"],w);
hmcPbPbxJTrueTagCorr->Fill(m[jtptSB]/m["jtpt1"],wb);
// hmcPbPbxJTrueTagCorrPt->Fill(m[jtptSB]/m["jtpt1"],w*corrpt);
// hmcPbPbxJTrueTagCorrEta->Fill(m[jtptSB]/m["jtpt1"],w*corrpt*correta);
// hmcPbPbxJTrueTagCorrBin->Fill(m[jtptSB]/m["jtpt1"],w*corrpt*correta*corrbin);
hpt2truetag->Fill(m[jtptSB],w);
hpt1truetag->Fill(m["jtpt1"],w);
heta2truetag->Fill(m[jtetaSB],w);
heta1truetag->Fill(m["jteta1"],w);
hbintruetag->Fill(m["bin"],w);
hpt2truetagcorr->Fill(m[jtptSB],wb);
hpt1truetagcorr->Fill(m["jtpt1"],wb);
heta2truetagcorr->Fill(m[jtetaSB],wb);
heta1truetagcorr->Fill(m["jteta1"],wb);
hbintruetagcorr->Fill(m["bin"],wb);
}
});
NormalizeAllHists();
//plotymax = 9999;
aktstring = TString::Format("PbPb %d-%d %%",binMin/2, binMax/2);//TString::Format("PbPb#Delta#phi>2/3#pi %d-%d %%",binMin/2, binMax/2);
SetMC({hmcPbPbxJTrue,hmcPbPbxJTrueTag,hmcPbPbxJTrueTagCorr});
SetData({hmcPbPbxJTrue});
hmcPbPbxJTrue->SetMinimum(0);
hmcPbPbxJTrue->SetMaximum(0.3);
hmcPbPbxJTrue->SetLineWidth(2);
hmcPbPbxJTrue->SetMarkerStyle(kNone);
hmcPbPbxJTrue->SetFillStyle(0);
hmcPbPbxJTrueTag->SetMarkerStyle(kOpenCircle);
hmcPbPbxJTrueTagCorr->SetMarkerStyle(kOpenSquare);
plotymax = 0.3;
Draw({hmcPbPbxJTrue,hmcPbPbxJTrueTag,hmcPbPbxJTrueTagCorr});
SetB({hmcPbPbxJTrue,hmcPbPbxJTrueTag,hmcPbPbxJTrueTagCorr});
float xjtrue = hmcPbPbxJTrue->GetMean();
float xjtruetag = hmcPbPbxJTrueTag->GetMean();
float xjtruetagcorr = hmcPbPbxJTrueTagCorr->GetMean();
float exjtrue = hmcPbPbxJTrue->GetMeanError();
float exjtruetag = hmcPbPbxJTrueTag->GetMeanError();
float exjtruetagcorr = hmcPbPbxJTrueTagCorr->GetMeanError();
auto c = getc();
hmcPbPbxJTrue->Draw("hist");
hmcPbPbxJTrueTag->Draw("E1,same");
hmcPbPbxJTrueTagCorr->Draw("E1,same");
plotlegendpos = TopLeft;
auto l = getLegend();
l->AddEntry(hmcPbPbxJTrue,Form("b-dijets, #LTx_{J}#GT=%.3f#pm%.3f",xjtrue,exjtrue),"L");
l->AddEntry(hmcPbPbxJTrueTag,Form("uncorrected, #LTx_{J}#GT=%.3f#pm%.3f",xjtruetag,exjtruetag),"P");
l->AddEntry(hmcPbPbxJTrueTagCorr,Form("corrected, #LTx_{J}#GT=%.3f#pm%.3f",xjtruetagcorr,exjtruetagcorr),"P");
l->Draw();
TLatex *Tl = new TLatex();
Tl->DrawLatexNDC(0.2, 0.8, aktstring);
SavePlot(c,Form("closure%d%d",binMin,binMax));
// //if (binMin==0 && binMax==200) {
// Draw({hmcPbPbxJTrueTag,hmcPbPbxJTrueTagCorrPt,hmcPbPbxJTrueTagCorrEta,hmcPbPbxJTrueTagCorrBin});
// SetMC({hpt2truetag,hpt1truetag,heta2truetag,heta1truetag,hbintruetag});
// plotputmean = false;
// plotymax = 0.2;
// Draw({hpt2true,hpt2truetag,hpt2truetagcorr});
// plotymax = 0.3;
// Draw({hpt1true,hpt1truetag,hpt1truetagcorr});
// plotymax = 0.2;
// Draw({heta2true,heta2truetag,heta2truetagcorr});
// Draw({heta1true,heta1truetag,heta1truetagcorr});
// plotymax = 1;
// Draw({hbintrue,hbintruetag,hbintruetagcorr});
plotymin = 0;
plotymax = 0.2;
Draw({hpt2truetag,hpt2true});
Draw({hpt2truetagcorr,hpt2true});
hpt2truetagovertrue->Divide(hpt2truetag,hpt2true,1,1); //"B"
hpt1truetagovertrue->Divide(hpt1truetag,hpt1true,1,1); //"B"
heta2truetagovertrue->Divide(heta2truetag,heta2true,1,1); //"B"
heta1truetagovertrue->Divide(heta1truetag,heta1true,1,1); //"B"
hbintruetagovertrue->Divide(hbintruetag,hbintrue,1,1); //"B"
hpt2truetagcorrovertrue->Divide(hpt2truetagcorr,hpt2true,1,1); //"B"
hpt1truetagcorrovertrue->Divide(hpt1truetagcorr,hpt1true,1,1); //"B"
heta2truetagcorrovertrue->Divide(heta2truetagcorr,heta2true,1,1); //"B"
heta1truetagcorrovertrue->Divide(heta1truetagcorr,heta1true,1,1); //"B"
hbintruetagcorrovertrue->Divide(hbintruetagcorr,hbintrue,1,1); //"B"
NormalizeAllHists();
Draw({hpt2truetagovertrue,hpt2truetagcorrovertrue});
Draw({hpt1truetagovertrue,hpt1truetagcorrovertrue});
Draw({heta2truetagovertrue,heta2truetagcorrovertrue});
Draw({heta1truetagovertrue,heta1truetagcorrovertrue});
Draw({hbintruetagovertrue,hbintruetagcorrovertrue});
// }
}
