]>
Commit | Line | Data |
---|---|---|
836ce14a | 1 | TH1D *GetHisto(float etacut = 0.12, bool eta = true, char *name, bool TPC,bool ITS, int mycase = 0, int color=1, int marker = 20, char *filename="Et.ESD.new.sim.merged.root",bool reweight = false,float kaonFactor=1.0, float lambdaFactor = 1.0, float baryonEnhancement = 1.0){ |
2 | TFile *file = new TFile(filename); | |
3 | TList *list = file->FindObject("out2"); | |
4 | char *reweightname = ""; | |
5 | if(reweight) reweightname = "Reweighted"; | |
6 | char *myname = "ITS"; | |
7 | if(TPC){ | |
8 | if(ITS) myname = "TPCITS"; | |
9 | else{ myname = "TPC";} | |
10 | } | |
bf40076b | 11 | //cout<<"Using "<<myname<<" tracks"<<endl; |
836ce14a | 12 | TH2F *signal = ((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedPiPlus",myname)))->Clone("signal"); |
13 | signal->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedPiMinus",myname))); | |
14 | signal->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedKMinus",myname)),kaonFactor); | |
15 | signal->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedKPlus",myname)),kaonFactor); | |
16 | signal->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedProton",myname)),baryonEnhancement); | |
17 | signal->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedAntiProton",myname)),baryonEnhancement); | |
18 | signal->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sUnidentifiedAssumingPion",myname))); | |
19 | ||
20 | //Et of all unidentified hadrons (plus hadrons identified as pions) calculated assuming their true mass | |
21 | TH2F *bkgd; | |
22 | switch(mycase){ | |
23 | case 0: | |
24 | bkgd = (TH2F*)((TH2F*) out2->FindObject(Form("EtReconstructed%sMisidentifiedElectrons",myname)))->Clone("bkgd"); | |
25 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sLambdaDaughters%s",myname,reweightname)),baryonEnhancement*lambdaFactor); | |
26 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sAntiLambdaDaughters%s",myname,reweightname)),baryonEnhancement*lambdaFactor); | |
27 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sK0SDaughters%s",myname,reweightname)),kaonFactor); | |
28 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sXiDaughters",myname)),baryonEnhancement); | |
29 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sAntiXiDaughters",myname)),baryonEnhancement); | |
30 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sOmegaDaughters",myname)),baryonEnhancement); | |
31 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sAntiOmegaDaughters",myname)),baryonEnhancement); | |
bf40076b | 32 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sConversionElectrons",myname)) ); |
33 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sSecondaryMuons",myname)) ); | |
34 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sSecondaryPions",myname)) ); | |
35 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sSecondaryProtons",myname)) ); | |
836ce14a | 36 | break; |
37 | case 1: | |
38 | bkgd = (TH2F*)((TH2F*) out2->FindObject(Form("EtReconstructed%sMisidentifiedElectrons",myname)))->Clone("bkgd"); | |
bf40076b | 39 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sConversionElectrons",myname)) ); |
836ce14a | 40 | break; |
41 | case 2: | |
42 | bkgd = (TH2F*)((TH2F*) out2->FindObject(Form("EtReconstructed%sLambdaDaughters%s",myname,reweightname)))->Clone("bkgd"); | |
43 | bkgd->Scale(baryonEnhancement*lambdaFactor); | |
44 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sAntiLambdaDaughters%s",myname,reweightname)),baryonEnhancement*lambdaFactor); | |
45 | break; | |
46 | case 3: | |
47 | bkgd = (TH2F*)((TH2F*) out2->FindObject(Form("EtReconstructed%sK0SDaughters%s",myname,reweightname)))->Clone("bkgd"); | |
48 | bkgd->Scale(kaonFactor); | |
49 | break; | |
50 | case 4: | |
51 | bkgd = (TH2F*)((TH2F*) out2->FindObject(Form("EtReconstructed%sLambdaDaughters%s",myname,reweightname)))->Clone("bkgd"); | |
52 | bkgd->Scale(baryonEnhancement*lambdaFactor); | |
53 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sAntiLambdaDaughters%s",myname,reweightname)),baryonEnhancement*lambdaFactor); | |
54 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sK0SDaughters%s",myname,reweightname)),kaonFactor); | |
55 | break; | |
56 | case 5: | |
57 | bkgd = (TH2F*)((TH2F*) out2->FindObject(Form("EtReconstructed%sXiDaughters",myname)))->Clone("bkgd"); | |
58 | bkgd->Scale(baryonEnhancement); | |
59 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sAntiXiDaughters",myname)),baryonEnhancement); | |
60 | break; | |
61 | case 6: | |
62 | bkgd = (TH2F*)((TH2F*) out2->FindObject(Form("EtReconstructed%sOmegaDaughters",myname)))->Clone("bkgd"); | |
63 | bkgd->Scale(baryonEnhancement); | |
64 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sAntiOmegaDaughters",myname)),baryonEnhancement); | |
65 | break; | |
66 | case 7: | |
67 | bkgd = (TH2F*)((TH2F*) out2->FindObject(Form("EtReconstructed%sXiDaughters",myname)))->Clone("bkgd"); | |
68 | bkgd->Scale(baryonEnhancement); | |
69 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sAntiXiDaughters",myname)),baryonEnhancement); | |
70 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sOmegaDaughters",myname)),baryonEnhancement); | |
71 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sAntiOmegaDaughters",myname)),baryonEnhancement); | |
72 | break; | |
b610fb37 | 73 | case 8: |
74 | bkgd = (TH2F*)((TH2F*) out2->FindObject(Form("EtReconstructed%sLambdaDaughters%s",myname,reweightname)))->Clone("bkgd"); | |
75 | bkgd->Scale(baryonEnhancement*lambdaFactor); | |
76 | break; | |
77 | case 9: | |
78 | bkgd = (TH2F*)((TH2F*) out2->FindObject(Form("EtReconstructed%sAntiLambdaDaughters%s",myname,reweightname)))->Clone("bkgd"); | |
79 | bkgd->Scale(baryonEnhancement*lambdaFactor); | |
80 | break; | |
bf40076b | 81 | case 10: |
82 | bkgd = (TH2F*) ((TH2F*) out2->FindObject(Form("EtReconstructed%sSecondaryMuons",myname)) )->Clone("bkgdOther"); | |
83 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sSecondaryPions",myname)) ); | |
84 | bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sSecondaryProtons",myname)) ); | |
85 | break; | |
836ce14a | 86 | } |
87 | TH1D *denominator; | |
88 | TH1D *numerator; | |
89 | if(eta){ | |
90 | int lowbin = bkgd->GetXaxis()->FindBin(etacut+.001);//make sure we don't accidentally get the wrong bin | |
91 | int highbin = bkgd->GetXaxis()->GetNbins(); | |
bf40076b | 92 | //cout<<"Projecting from "<<bkgd->GetXaxis()->GetBinLowEdge(lowbin)<<" to "<<bkgd->GetXaxis()->GetBinLowEdge(highbin+1)<<endl; |
836ce14a | 93 | denominator = signal->ProjectionY("name",lowbin,highbin); |
94 | numerator = bkgd->ProjectionY(name,lowbin,highbin); | |
95 | } | |
96 | else{ | |
97 | int lowbin = bkgd->GetYaxis()->FindBin(-etacut+.001);//make sure we don't accidentally get the wrong bin | |
98 | int highbin = bkgd->GetYaxis()->FindBin(etacut-.001); | |
bf40076b | 99 | //cout<<"Projecting from "<<bkgd->GetYaxis()->GetBinLowEdge(lowbin)<<" to "<<bkgd->GetYaxis()->GetBinLowEdge(highbin+1)<<endl; |
836ce14a | 100 | denominator = signal->ProjectionX("name",lowbin,highbin); |
101 | numerator = bkgd->ProjectionX(name,lowbin,highbin); | |
102 | } | |
bf40076b | 103 | denominator->Add(numerator); |
836ce14a | 104 | numerator->Divide(denominator); |
105 | numerator->SetYTitle("Ratio of E_{T}^{background}/E_{T}^{had, meas.}"); | |
106 | numerator->GetYaxis()->SetTitleOffset(1.2); | |
107 | numerator->SetLineColor(color); | |
108 | numerator->SetMarkerColor(color); | |
109 | numerator->SetMarkerStyle(marker); | |
110 | return numerator; | |
111 | ||
112 | } | |
113 | ||
bf40076b | 114 | void CorrBkgdErrors(char *prodname = "LHC11a4_bis HIJING 2.76 TeV Pb+Pb",char *shortprodname= "LHC11a4_bis", char *filename = "rootFiles/LHC11a10a_bis/Et.ESD.sim.LHC11a10a_bis.Run139465.root",bool alt = false,bool TPC = true, bool ITS=true, bool reweight = false, int scale = -1,float kaonFactor=1.0, float lambdaFactor = 1.0, float baryonEnhancement = 1.0){ |
836ce14a | 115 | gStyle->SetOptTitle(0); |
116 | gStyle->SetOptStat(0); | |
117 | gStyle->SetOptFit(0); | |
118 | TCanvas *c = new TCanvas("c","c",400,400); | |
119 | c->SetTopMargin(0.04); | |
120 | c->SetRightMargin(0.04); | |
121 | c->SetBorderSize(0); | |
122 | c->SetFillColor(0); | |
123 | c->SetFillColor(0); | |
124 | c->SetBorderMode(0); | |
125 | c->SetFrameFillColor(0); | |
126 | c->SetFrameBorderMode(0); | |
127 | c->SetLeftMargin(0.159091); | |
bf40076b | 128 | //char * filename = "Et.ESD.new.sim.LHC10d4.pp.merged.root"; |
129 | // char * filename = "rootFiles/LHC11a4_bis/Et.ESD.new.sim.LHC11a4_bis.root"; | |
130 | //char *filename = "rootFiles/LHC11a4_bis/Et.ESD.new.sim.LHC11a4_bis.root"; | |
131 | //char *filename = "rootFiles/LHC11a4_bis/Et.ESD.new.sim.LHC11a4_bis.effFromPtReco.root"; | |
132 | //char *filename = "rootFiles/LHC11a4_bis23Apr11/Et.ESD.new.sim.LHC11a4_bis.root"; | |
836ce14a | 133 | TH1D *All = GetHisto(0.1,true,"All",TPC,ITS,0,1,20,filename,reweight,kaonFactor,lambdaFactor,baryonEnhancement); |
134 | TH1D *Electrons = GetHisto(0.1,true,"Electrons",TPC,ITS,1,2,21,filename,reweight,kaonFactor,lambdaFactor,baryonEnhancement); | |
135 | TH1D *V0s = GetHisto(0.1,true,"V0s",TPC,ITS,4,4,22,filename,reweight,kaonFactor,lambdaFactor,baryonEnhancement); | |
b610fb37 | 136 | //case, color, marker |
137 | TH1D *K0S = GetHisto(0.1,true,"K0S",TPC,ITS,3,TColor::kOrange+8,33,filename,reweight,kaonFactor,lambdaFactor,baryonEnhancement); | |
138 | TH1D *Lambda = GetHisto(0.1,true,"K0S",TPC,ITS,8,TColor::kMagenta+3,29,filename,reweight,kaonFactor,lambdaFactor,baryonEnhancement); | |
139 | TH1D *AntiLambda = GetHisto(0.1,true,"K0S",TPC,ITS,9,TColor::kMagenta+3,30,filename,reweight,kaonFactor,lambdaFactor,baryonEnhancement); | |
836ce14a | 140 | TH1D *Multistrange = GetHisto(0.1,true,"Multistrange",TPC,ITS,7,TColor::kGreen+2,23,filename,reweight,kaonFactor,lambdaFactor,baryonEnhancement); |
141 | ||
142 | TH1D *Allpt = GetHisto(0.7,false,"Allpt",TPC,ITS,0,1,20,filename,reweight,kaonFactor,lambdaFactor,baryonEnhancement); | |
143 | TH1D *Electronspt = GetHisto(0.7,false,"Electronspt",TPC,ITS,1,2,21,filename,reweight,kaonFactor,lambdaFactor,baryonEnhancement); | |
144 | TH1D *V0spt = GetHisto(0.7,false,"V0spt",TPC,ITS,4,4,22,filename,reweight,kaonFactor,lambdaFactor,baryonEnhancement); | |
b610fb37 | 145 | TH1D *K0Spt = GetHisto(0.1,false,"K0S",TPC,ITS,3,TColor::kOrange+8,33,filename,reweight,kaonFactor,lambdaFactor,baryonEnhancement); |
146 | TH1D *Lambdapt = GetHisto(0.1,false,"K0S",TPC,ITS,8,TColor::kMagenta+3,29,filename,reweight,kaonFactor,lambdaFactor,baryonEnhancement); | |
147 | TH1D *AntiLambdapt = GetHisto(0.1,false,"K0S",TPC,ITS,9,TColor::kMagenta+3,30,filename,reweight,kaonFactor,lambdaFactor,baryonEnhancement); | |
836ce14a | 148 | TH1D *Multistrangept = GetHisto(0.7,false,"Multistrangept",TPC,ITS,7,TColor::kGreen+2,23,filename,reweight,kaonFactor,lambdaFactor,baryonEnhancement); |
149 | if(ITS){ | |
bf40076b | 150 | if(TPC) All->SetMaximum(0.15); |
151 | else{All->SetMaximum(0.025);} | |
836ce14a | 152 | } |
153 | else{ | |
154 | All->SetMaximum(0.2); | |
155 | } | |
156 | All->SetMinimum(0.0); | |
157 | All->SetMarkerStyle(20); | |
158 | All->GetYaxis()->SetTitleOffset(1.8); | |
159 | TF1 *func = new TF1("func","[0]",-.7,.7); | |
b610fb37 | 160 | func->SetParameter(0,0.02); |
161 | TF1 *funcLam = new TF1("funcLam","[0]",-.7,.7); | |
162 | funcLam->SetParameter(0,0.001); | |
163 | funcLam->SetLineColor(Lambda->GetMarkerColor()); | |
164 | TF1 *funcAlam = new TF1("funcAlam","[0]",-.7,.7); | |
165 | funcAlam->SetParameter(0,0.003); | |
166 | funcAlam->SetLineColor(AntiLambda->GetMarkerColor()); | |
167 | TF1 *funcK0 = new TF1("funcK0","[0]",-.7,.7); | |
168 | funcK0->SetParameter(0,0.013); | |
169 | funcK0->SetLineColor(K0S->GetMarkerColor()); | |
bf40076b | 170 | TF1 *funcElectrons = new TF1("funcElectrons","[0]",-.7,.7); |
171 | funcElectrons->SetParameter(0,0.013); | |
172 | funcElectrons->SetLineColor(Electrons->GetMarkerColor()); | |
836ce14a | 173 | All->Fit(func); |
b610fb37 | 174 | Lambda->Fit(funcLam); |
175 | AntiLambda->Fit(funcAlam); | |
176 | K0S->Fit(funcK0); | |
bf40076b | 177 | Electrons->Fit(funcElectrons); |
178 | float k0Err900 = 0.07/0.184; | |
179 | float lamErr900 = 0.005/0.048; | |
180 | float alamErr900 = 0.006/0.047; | |
181 | float k0Err7 = 0.03/0.347; | |
182 | float lamErr7 = 0.02/0.182; | |
183 | float alamErr7 = lamErr7; | |
184 | ||
185 | float lambda = funcLam->GetParameter(0); | |
186 | cout<<"Lambda contribution 900: "; | |
187 | cout<<Form("%2.2f $\\pm$ %2.2f",lambda*100.0,lambda*lamErr900*100.0)<<endl; | |
188 | cout<<"Lambda contribution 7 "; | |
189 | cout<<Form("%2.2f $\\pm$ %2.2f",lambda*100.0,lambda*lamErr7*100.0)<<endl; | |
190 | ||
191 | float antilambda = funcAlam->GetParameter(0); | |
192 | cout<<"AntiLambda contribution 900: "; | |
193 | cout<<Form("%2.2f $\\pm$ %2.2f",antilambda*100.0,antilambda*alamErr900*100.0)<<endl; | |
194 | cout<<"AntiLambda contribution 7 "; | |
195 | cout<<Form("%2.2f $\\pm$ %2.2f",antilambda*100.0,antilambda*alamErr7*100.0)<<endl; | |
196 | ||
197 | float k0 = funcK0->GetParameter(0); | |
198 | cout<<"K0 contribution 900: "; | |
199 | cout<<Form("%2.2f $\\pm$ %2.2f",k0*100.0,k0*k0Err900*100.0)<<endl; | |
200 | cout<<"K0 contribution 7 "; | |
201 | cout<<Form("%2.2f $\\pm$ %2.2f",k0*100.0,k0*k0Err7*100.0)<<endl; | |
202 | ||
203 | cout<<"Total 900: "; | |
204 | float totalerr900 = 100.0*TMath::Sqrt(lambda*lamErr900*lambda*lamErr900+antilambda*alamErr900*antilambda*alamErr900+k0*k0Err900*k0*k0Err900); | |
205 | float total900 = lambda*100.0+antilambda*100.0+k0*100.0; | |
206 | cout<<Form("%2.2f $\\pm$ %2.2f",total900,totalerr900)<<endl; | |
207 | cout<<"Total 7: "; | |
208 | float totalerr7 = 100.0*TMath::Sqrt(lambda*lamErr7*lambda*lamErr7+antilambda*alamErr7*antilambda*alamErr7+k0*k0Err7*k0*k0Err7); | |
209 | float total7 = lambda*100.0+antilambda*100.0+k0*100.0; | |
210 | cout<<Form("%2.2f $\\pm$ %2.2f",total7,totalerr7)<<endl; | |
211 | ||
212 | ||
213 | //cout<<funcLam->GetParameter(0)<<" +/- "<<funcLam->GetParError(0)<<endl; | |
214 | cout<<"AntiLambda contribution "; | |
215 | //cout<<funcAlam->GetParameter(0)<<" +/- "<<funcAlam->GetParError(0)<<endl; | |
216 | cout<<"K0 contribution "; | |
217 | //cout<<funcK0->GetParameter(0)<<" +/- "<<funcK0->GetParError(0)<<endl; | |
218 | ||
219 | ||
220 | cout<<"Electron contribution "; | |
221 | //cout<<funcElectrons->GetParameter(0)<<" +/- "<<funcElectrons->GetParError(0)<<endl; | |
836ce14a | 222 | All->Draw(); |
223 | Electrons->Draw("same"); | |
224 | V0s->Draw("same"); | |
225 | Multistrange->Draw("same"); | |
b610fb37 | 226 | K0S->Draw("same"); |
227 | Lambda->Draw("same"); | |
228 | AntiLambda->Draw("same"); | |
836ce14a | 229 | TLatex *tex = new TLatex(0.161478,1.0835,"LHC10d15: p+p, Pythia6 Perugia-0"); |
230 | tex->SetTextSize(0.0537634); | |
231 | tex->Draw(); | |
b610fb37 | 232 | TLegend *leg = new TLegend(0.636364,0.612903,0.893939,0.962366); |
836ce14a | 233 | leg->AddEntry(All,"All"); |
234 | leg->AddEntry(Electrons,"Electrons"); | |
235 | leg->AddEntry(V0s,"V0s"); | |
b610fb37 | 236 | leg->AddEntry(K0S,"K_{S}^{0}"); |
237 | leg->AddEntry(Lambda,"#Lambda"); | |
238 | leg->AddEntry(AntiLambda,"#bar{#Lambda}"); | |
836ce14a | 239 | leg->AddEntry(Multistrange,"Multistrange"); |
240 | leg->SetFillStyle(0); | |
241 | leg->SetFillColor(0); | |
242 | leg->SetBorderSize(0); | |
243 | leg->Draw(); | |
b610fb37 | 244 | TLatex *tex = new TLatex(-0.711139,0.0157696,Form("%2.5f#pm%2.5f",func->GetParameter(0),func->GetParError(0))); |
836ce14a | 245 | tex->Draw(); |
b610fb37 | 246 | TLatex *texLam = new TLatex(-0.711139,0.00201613,Form("%2.5f#pm%2.5f",funcLam->GetParameter(0),funcLam->GetParError(0))); |
247 | texLam->SetTextColor(Lambda->GetMarkerColor()); | |
248 | texLam->Draw(); | |
249 | TLatex *texAlam = new TLatex(-0.711139,0.00365716,Form("%2.5f#pm%2.5f",funcAlam->GetParameter(0),funcAlam->GetParError(0))); | |
250 | texAlam->SetTextColor(AntiLambda->GetMarkerColor()); | |
251 | texAlam->Draw(); | |
252 | TLatex *texK0 = new TLatex(-0.711139,0.008,Form("%2.5f#pm%2.5f",funcK0->GetParameter(0),funcK0->GetParError(0))); | |
253 | texK0->SetTextColor(K0S->GetMarkerColor()); | |
254 | texK0->Draw(); | |
bf40076b | 255 | |
256 | ||
836ce14a | 257 | TCanvas *c1 = new TCanvas("c1","c1",600,400); |
258 | c1->SetTopMargin(0.04); | |
259 | c1->SetRightMargin(0.04); | |
260 | c1->SetBorderSize(0); | |
261 | c1->SetFillColor(0); | |
262 | c1->SetFillColor(0); | |
263 | c1->SetBorderMode(0); | |
264 | c1->SetFrameFillColor(0); | |
265 | c1->SetFrameBorderMode(0); | |
bf40076b | 266 | if(TPC && ITS){ |
267 | Allpt->SetMaximum(0.15); | |
268 | Allpt->GetXaxis()->SetRange(Allpt->GetXaxis()->FindBin(0.15),Allpt->GetXaxis()->FindBin(1.)); | |
836ce14a | 269 | } |
270 | else{ | |
bf40076b | 271 | Allpt->GetXaxis()->SetRange(Allpt->GetXaxis()->FindBin(0.1),Allpt->GetXaxis()->FindBin(0.3)); |
836ce14a | 272 | } |
836ce14a | 273 | Allpt->SetMinimum(0.0); |
274 | Allpt->SetMarkerStyle(20); | |
275 | Allpt->Draw(); | |
276 | Electronspt->Draw("same"); | |
277 | V0spt->Draw("same"); | |
b610fb37 | 278 | K0Spt->Draw("same"); |
279 | Lambdapt->Draw("same"); | |
280 | AntiLambdapt->Draw("same"); | |
836ce14a | 281 | Multistrangept->Draw("same"); |
282 | TLatex *texpt = new TLatex(0.161478,1.0835,"LHC10d15: p+p, Pythia6 Perugia-0"); | |
283 | texpt->SetTextSize(0.0537634); | |
284 | texpt->Draw(); | |
b610fb37 | 285 | TLegend *legpt = new TLegend(0.634228,0.430108,0.892617,0.905914); |
836ce14a | 286 | legpt->AddEntry(Allpt,"All"); |
287 | legpt->AddEntry(Electronspt,"Electrons"); | |
288 | legpt->AddEntry(V0spt,"V0s"); | |
b610fb37 | 289 | legpt->AddEntry(K0Spt,"K_{S}^{0}"); |
290 | legpt->AddEntry(Lambdapt,"#Lambda"); | |
291 | legpt->AddEntry(AntiLambdapt,"#bar{#Lambda}"); | |
836ce14a | 292 | legpt->AddEntry(Multistrangept,"Multistrange"); |
293 | legpt->SetFillStyle(0); | |
294 | legpt->SetFillColor(0); | |
295 | legpt->SetBorderSize(0); | |
296 | legpt->Draw(); | |
bf40076b | 297 | float x = 0.399919; |
298 | if(reweight) x = 0.18332; | |
299 | TLatex *texprod = new TLatex(x,0.140566,prodname); | |
300 | texprod->SetTextSize(0.0537634); | |
301 | texprod->Draw(); | |
302 | if(scale!=-1){ | |
303 | TString scale900 = "Using scaling from 900 GeV p+p"; | |
304 | TString scale7 = "Using scaling from 7 TeV p+p"; | |
305 | TString scaleString; | |
306 | if(scale==2009) scaleString = scale900; | |
307 | else{scaleString = scale7;} | |
308 | TLatex *texscale = new TLatex(0.190421,0.129782,scaleString.Data()); | |
309 | texscale->SetTextSize(0.0537634); | |
310 | texscale->Draw(); | |
311 | } | |
836ce14a | 312 | |
313 | char TPCnameeps[200]; | |
314 | char TPCnamepng[200]; | |
315 | char ITSnameeps[200]; | |
316 | char ITSnamepng[200]; | |
317 | TString *None = new TString(""); | |
318 | TString *Factors = None; | |
319 | if(kaonFactor!=1.0||lambdaFactor!=1.0||baryonEnhancement!=1.0){ | |
320 | Factors = new TString(Form("Lambda%2.1fKaon%2.1fBaryon%2.1f",lambdaFactor,kaonFactor,baryonEnhancement)); | |
321 | } | |
bf40076b | 322 | TString dir = "/home/madscientistgirl/UTK/Et/aliceet/AN/pics/hadEt"; |
323 | TString detector = "TPC"; | |
324 | if(!TPC) detector = "ITS"; | |
325 | TString altName = ""; | |
326 | if(alt) altName = "Alt"; | |
327 | TString bary = ""; | |
328 | if(reweight) bary = "Enhanced"; | |
329 | gSystem->MakeDirectory(Form("%s/%s",dir.Data(),shortprodname)); | |
330 | sprintf(TPCnameeps,"%s/%s/bkgdComponentsErrors%s%s%s%s.eps",dir.Data(),shortprodname,detector.Data(),Factors->Data(),altName.Data(),bary.Data()); | |
331 | sprintf(TPCnamepng,"%s/%s/bkgdComponentsErrors%s%s%s%s.png",dir.Data(),shortprodname,detector.Data(),Factors->Data(),altName.Data(),bary.Data()); | |
332 | c->SaveAs(TPCnameeps); | |
333 | //c->SaveAs(TPCnamepng); | |
334 | sprintf(TPCnameeps,"%s/%s/bkgdComponents%s%s%s%s.eps",dir.Data(),shortprodname,detector.Data(),Factors->Data(),altName.Data(),bary.Data()); | |
335 | sprintf(TPCnamepng,"%s/%s/bkgdComponents%s%s%s%s.png",dir.Data(),shortprodname,detector.Data(),Factors->Data(),altName.Data(),bary.Data()); | |
336 | c1->SaveAs(TPCnameeps); | |
337 | //c1->SaveAs(TPCnamepng); | |
338 | // if(TPC){ | |
339 | // sprintf(TPCnameeps,"pics/bkgdComponentsErrorsTPC%s.eps",Factors->Data()); | |
340 | // sprintf(TPCnamepng,"pics/bkgdComponentsErrorsTPC%s.png",Factors->Data()); | |
341 | // c->SaveAs(TPCnameeps); | |
342 | // c->SaveAs(TPCnamepng); | |
343 | // sprintf(TPCnameeps,"pics/bkgdComponentsTPC%s.eps",Factors->Data()); | |
344 | // sprintf(TPCnamepng,"pics/bkgdComponentsTPC%s.png",Factors->Data()); | |
345 | // c1->SaveAs(TPCnameeps); | |
346 | // c1->SaveAs(TPCnamepng); | |
347 | // } | |
348 | // else{ | |
349 | // sprintf(ITSnameeps,"pics/bkgdComponentsErrorsITS%s.eps",Factors->Data()); | |
350 | // sprintf(ITSnamepng,"pics/bkgdComponentsErrorsITS%s.png",Factors->Data()); | |
351 | // c->SaveAs(ITSnameeps); | |
352 | // c->SaveAs(ITSnamepng); | |
353 | // sprintf(ITSnameeps,"pics/bkgdComponentsITS%s.eps",Factors->Data()); | |
354 | // sprintf(ITSnamepng,"pics/bkgdComponentsITS%s.png",Factors->Data()); | |
355 | // c1->SaveAs(ITSnameeps); | |
356 | // c1->SaveAs(ITSnamepng); | |
357 | // } | |
836ce14a | 358 | |
bf40076b | 359 | // TCanvas *c4 = new TCanvas("c4","c4",600,400); |
360 | // c4->SetTopMargin(0.04); | |
361 | // c4->SetRightMargin(0.04); | |
362 | // c4->SetBorderSize(0); | |
363 | // c4->SetFillColor(0); | |
364 | // c4->SetFillColor(0); | |
365 | // c4->SetBorderMode(0); | |
366 | // c4->SetFrameFillColor(0); | |
367 | // c4->SetFrameBorderMode(0); | |
368 | // Allpt->Draw(); | |
369 | // c4->SaveAs("pics/bkgd276ppV07TeVITS.C"); | |
836ce14a | 370 | } |