| becfe989 |
1 | #include <Riostream.h> |
| 2 | |
| 3 | extern TRandom *gRandom; |
| 4 | extern TBenchmark *gBenchmark; |
| 5 | extern TSystem *gSystem; |
| 6 | |
| 7 | void ReadCFHeavyFlavourOutput(){ |
| 8 | |
| 9 | // example macro for reading a container coming from a CF analysis |
| 10 | // the generated and reconstructed distributions are got, |
| 11 | // the efficiencies are calculated, |
| 12 | // the reco distributions are corrected, |
| 13 | // the results plotted |
| 14 | |
| 15 | gROOT->SetStyle("Plain"); |
| 16 | gStyle->SetPalette(1); |
| 17 | gStyle->SetOptStat(1110); |
| 18 | gStyle->SetPalette(1); |
| 19 | gStyle->SetCanvasColor(0); |
| 20 | gStyle->SetFrameFillColor(0); |
| 21 | gStyle->SetOptTitle(0); |
| 22 | |
| 23 | gSystem->SetIncludePath("-I. -I$ALICE_ROOT/include -I$ROOTSYS/include"); |
| 24 | gSystem->Load("libANALYSIS.so"); |
| 25 | gSystem->Load("libANALYSISalice.so"); |
| a273622f |
26 | gSystem->Load("libCORRFW.so") ; |
| becfe989 |
27 | |
| 28 | //Setting up the container grid... |
| 29 | |
| 30 | const Int_t nstep=2; //number of selection steps (just 2 in this ex) |
| 31 | |
| 32 | const Int_t nvar=6; //number of variables on the grid: pt, y, cosThetaStar, ptPi, ptK, cT |
| 33 | |
| 34 | // Flag the sel steps. In this example, we have two, may be any nstep |
| 35 | Int_t stepGen=0; |
| 36 | Int_t stepRec=1; |
| 37 | |
| 38 | // the sensitive variables, their indeces |
| 39 | UInt_t ipt =0; |
| 40 | UInt_t iy =1; |
| 41 | UInt_t icTS = 2; // cTS stands for cosThetaStar |
| 42 | UInt_t iptPi= 3; |
| 43 | UInt_t iptK = 4; |
| 44 | UInt_t icT = 5; |
| 45 | |
| 46 | // Read the container from file |
| 47 | TFile *file = new TFile("output.root"); |
| 48 | AliCFContainer *data = (AliCFContainer*) (file->Get("container")); |
| 49 | |
| 50 | // Make some 1 & 2-D projections.. |
| 51 | // MC level |
| 52 | TCanvas *cmc1 = new TCanvas("cmc1","The MC distributions",0,800,900,1200); |
| 53 | TCanvas *cmc2 = new TCanvas("cmc2","The MC distributions",0,800,900,1200); |
| 54 | TCanvas *cmcpt = new TCanvas("cmcpt","pt distribution from MC",50,50,550,550); |
| 55 | TCanvas *cmcy = new TCanvas("cmcy","y distribution from MC",50,50,550,550); |
| 56 | TCanvas *cmccTS = new TCanvas("cmcTS","cosThetaStar distribution from MC",50,50,550,550); // cTS stands for cosThetaStar |
| 57 | TCanvas *cmcptPi = new TCanvas("cmcptPi","pt_pi distribution from MC",50,50,550,550); |
| 58 | TCanvas *cmcptK = new TCanvas("cmcptK","pt_K distribution from MC",50,50,550,550); |
| 59 | TCanvas *cmccT = new TCanvas("cmccT","cT distribution from MC",50,50,550,550); |
| 60 | |
| 61 | // Reco-aod level |
| 62 | TCanvas *cpt = new TCanvas("cpt","pt distribution from reco aod",50,50,550,550); |
| 63 | TCanvas *cy = new TCanvas("cy","y distribution from reco aod",50,50,550,550); |
| 64 | TCanvas *ccTS = new TCanvas("cTS","cosThetaStar distribution from reco aod",50,50,550,550); |
| 65 | TCanvas *cptPi = new TCanvas("cptPi","pt_pi distribution from reco aod",50,50,550,550); |
| 66 | TCanvas *cptK = new TCanvas("cptK","pt_K distribution from reco aod",50,50,550,550); |
| 67 | TCanvas *ccT = new TCanvas("ccT","cT distribution from reco aod",50,50,550,550); |
| 68 | |
| 69 | // some 2D distributions in pt and y |
| 70 | TCanvas *cGen2d = new TCanvas("cGen2d","2D distribution from MC ",50,50,550,550); |
| 71 | TCanvas *cRec2d = new TCanvas("cRec2d","2D distribution from reco aod",50,50,550,550); |
| 72 | |
| 73 | // 2D plots |
| 74 | TH2D *hMCpty2d = data->ShowProjection(ipt, iy, stepGen); |
| 75 | hMCpty2d->SetLineColor(2); |
| 76 | hMCpty2d->SetLineWidth(3); |
| 77 | hMCpty2d->SetMarkerColor(2); |
| 78 | hMCpty2d->SetMarkerStyle(20); |
| 79 | hMCpty2d->GetXaxis()->SetTitleOffset(1.2); |
| 80 | hMCpty2d->GetYaxis()->SetTitleOffset(1.5); |
| 81 | hMCpty2d->GetXaxis()->SetTitle("p_{T} (GeV/c), MC data"); |
| 82 | hMCpty2d->GetYaxis()->SetTitle("y, MC data"); |
| 83 | hMCpty2d->Draw("text"); |
| 84 | cGen2d->cd(); |
| 85 | cGen2d->SetLeftMargin(0.15); |
| 86 | cGen2d->SetRightMargin(0.05); |
| 87 | cGen2d->Update(); |
| 88 | |
| 89 | TH2D *hRECpty2d = data->ShowProjection(ipt, iy, stepRec); |
| 90 | hRECpty2d->Sumw2(); |
| 91 | hRECpty2d->SetLineColor(4); |
| 92 | hRECpty2d->SetLineWidth(3); |
| 93 | hRECpty2d->SetMarkerColor(4); |
| 94 | hRECpty2d->SetMarkerStyle(20); |
| 95 | hRECpty2d->GetXaxis()->SetTitleOffset(1.2); |
| 96 | hRECpty2d->GetYaxis()->SetTitleOffset(1.5); |
| 97 | hRECpty2d->GetXaxis()->SetTitle("p_{T} (GeV/c), AOD"); |
| 98 | hRECpty2d->GetYaxis()->SetTitle("y, AOD"); |
| 99 | hRECpty2d->Draw("text"); |
| 100 | cRec2d->cd(); |
| 101 | cRec2d->SetLeftMargin(0.15); |
| 102 | cRec2d->SetRightMargin(0.05); |
| 103 | cRec2d->Update(); |
| 104 | |
| 105 | // MC + REC 1D plots |
| 106 | // pt, y, cosThetaStar |
| 107 | cmc1->Divide(2,3); |
| 108 | |
| 109 | cmc1->cd(1); |
| 110 | TH1D *hMCpt1D = data->ShowProjection(ipt, stepGen); |
| 111 | Double_t maxpt = hMCpt1D->GetMaximum(); |
| 112 | hMCpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2); |
| 113 | hMCpt1D->Sumw2(); |
| 114 | hMCpt1D->SetMinimum(0.01); |
| 115 | hMCpt1D->SetLineColor(2); |
| 116 | hMCpt1D->SetLineWidth(3); |
| 117 | hMCpt1D->SetMarkerColor(2); |
| 118 | hMCpt1D->SetFillColor(2); |
| 119 | hMCpt1D->SetFillStyle(3005); |
| 120 | hMCpt1D->SetMarkerStyle(20); |
| 121 | hMCpt1D->GetXaxis()->SetTitleOffset(1.2); |
| 122 | hMCpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), MC data"); |
| 123 | hMCpt1D->Draw("hist"); |
| 124 | cmcpt->cd(); |
| 125 | cmcpt->SetLeftMargin(0.15); |
| 126 | cmcpt->SetRightMargin(0.05); |
| 127 | hMCpt1D->Draw("hist"); |
| 128 | hMCpt1D->Draw("err same"); |
| 129 | cmcpt->Update(); |
| 130 | |
| 131 | cmc1->cd(2); |
| 132 | TH1D *hRECpt1D = data->ShowProjection(ipt, stepRec); |
| 133 | hRECpt1D->GetYaxis()->SetRangeUser(0,maxpt*1.2); |
| 134 | hRECpt1D->SetLineColor(4); |
| 135 | hRECpt1D->SetLineWidth(3); |
| 136 | hRECpt1D->SetMarkerColor(4); |
| 137 | hRECpt1D->SetFillColor(4); |
| 138 | hRECpt1D->SetFillStyle(3004); |
| 139 | hRECpt1D->SetMarkerStyle(20); |
| 140 | hRECpt1D->GetXaxis()->SetTitleOffset(1.2); |
| 141 | hRECpt1D->Sumw2(); |
| 142 | hRECpt1D->SetMinimum(0.01); |
| 143 | hRECpt1D->GetXaxis()->SetTitle("p_{T} (GeV/c), AOD"); |
| 144 | hRECpt1D->Draw("hist"); |
| 145 | cpt->cd(); |
| 146 | cpt->SetLeftMargin(0.15); |
| 147 | cpt->SetRightMargin(0.05); |
| 148 | hRECpt1D->Draw("hist"); |
| 149 | hRECpt1D->Draw("err same"); |
| 150 | cpt->Update(); |
| 151 | |
| 152 | cmc1->cd(3); |
| 153 | TH1D *hMCy1D = data->ShowProjection(iy, stepGen); |
| 154 | Double_t maxy = hMCy1D->GetMaximum(); |
| 155 | hMCy1D->GetYaxis()->SetRangeUser(0,maxy*1.2); |
| 156 | hMCy1D->SetLineColor(2); |
| 157 | hMCy1D->SetLineWidth(3); |
| 158 | hMCy1D->SetMarkerColor(2); |
| 159 | hMCy1D->SetFillColor(2); |
| 160 | hMCy1D->SetFillStyle(3005); |
| 161 | hMCy1D->SetMarkerStyle(20); |
| 162 | hMCy1D->GetXaxis()->SetTitleOffset(1.2); |
| 163 | hMCy1D->Sumw2(); |
| 164 | hMCy1D->SetMinimum(0.01); |
| 165 | hMCy1D->GetXaxis()->SetTitle("y, MC data"); |
| 166 | hMCy1D->Draw("hist"); |
| 167 | cmcy->cd(); |
| 168 | cmcy->SetLeftMargin(0.15); |
| 169 | cmcy->SetRightMargin(0.05); |
| 170 | hMCy1D->Draw("hist"); |
| 171 | hMCy1D->Draw("err same"); |
| 172 | cmcy->Update(); |
| 173 | |
| 174 | cmc1->cd(4); |
| 175 | TH1D *hRECy1D = data->ShowProjection(iy, stepRec); |
| 176 | hRECy1D->GetYaxis()->SetRangeUser(0,maxy*1.2); |
| 177 | hRECy1D->SetLineColor(4); |
| 178 | hRECy1D->SetLineWidth(3); |
| 179 | hRECy1D->SetMarkerColor(4); |
| 180 | hRECy1D->SetFillColor(4); |
| 181 | hRECy1D->SetFillStyle(3004); |
| 182 | hRECy1D->SetMarkerStyle(20); |
| 183 | hRECy1D->GetXaxis()->SetTitleOffset(1.2); |
| 184 | hRECy1D->Sumw2(); |
| 185 | hRECy1D->SetMinimum(0.01); |
| 186 | hRECy1D->Draw("hist"); |
| 187 | cy->cd(); |
| 188 | cy->SetLeftMargin(0.15); |
| 189 | cy->SetRightMargin(0.05); |
| 190 | hRECy1D->GetXaxis()->SetTitle("y, AOD"); |
| 191 | hRECy1D->Draw("hist"); |
| 192 | hRECy1D->Draw("err same"); |
| 193 | cy->Update(); |
| 194 | |
| 195 | cmc1->cd(5); |
| 196 | TH1D *hMCcTS1D = data->ShowProjection(icTS, stepGen); |
| 197 | Double_t maxcTS = hMCcTS1D->GetMaximum(); |
| 198 | hMCcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2); |
| 199 | hMCcTS1D->SetLineColor(2); |
| 200 | hMCcTS1D->SetLineWidth(3); |
| 201 | hMCcTS1D->SetMarkerColor(2); |
| 202 | hMCcTS1D->SetFillColor(2); |
| 203 | hMCcTS1D->SetFillStyle(3005); |
| 204 | hMCcTS1D->SetMarkerStyle(20); |
| 205 | hMCcTS1D->GetXaxis()->SetTitleOffset(1.2); |
| 206 | hMCcTS1D->Sumw2(); |
| 207 | hMCcTS1D->SetMinimum(0.01); |
| 208 | hMCcTS1D->GetXaxis()->SetTitle("cosThetaStar, MC data"); |
| 209 | hMCcTS1D->Draw("hist"); |
| 210 | cmccTS->cd(); |
| 211 | cmccTS->SetLeftMargin(0.15); |
| 212 | cmccTS->SetRightMargin(0.05); |
| 213 | hMCcTS1D->Draw("hist"); |
| 214 | hMCcTS1D->Draw("err same"); |
| 215 | cmccTS->Update(); |
| 216 | |
| 217 | cmc1->cd(6); |
| 218 | TH1D *hRECcTS1D = data->ShowProjection(icTS, stepRec); |
| 219 | hRECcTS1D->GetYaxis()->SetRangeUser(0,maxcTS*1.2); |
| 220 | hRECcTS1D->SetLineColor(4); |
| 221 | hRECcTS1D->SetLineWidth(3); |
| 222 | hRECcTS1D->SetMarkerColor(4); |
| 223 | hRECcTS1D->SetFillColor(4); |
| 224 | hRECcTS1D->SetFillStyle(3004); |
| 225 | hRECcTS1D->SetMarkerStyle(20); |
| 226 | hRECcTS1D->GetXaxis()->SetTitleOffset(1.2); |
| 227 | hRECcTS1D->Sumw2(); |
| 228 | hRECcTS1D->SetMinimum(0.01); |
| 229 | hRECcTS1D->Draw("hist"); |
| 230 | ccTS->cd(); |
| 231 | ccTS->SetLeftMargin(0.15); |
| 232 | ccTS->SetRightMargin(0.05); |
| 233 | hRECcTS1D->GetXaxis()->SetTitle("cosThetaStar, AOD"); |
| 234 | hRECcTS1D->Draw("hist"); |
| 235 | hRECcTS1D->Draw("err same"); |
| 236 | ccTS->Update(); |
| 237 | |
| 238 | // ptPi, ptK, cT |
| 239 | cmc2->Divide(2,3); |
| 240 | |
| 241 | cmc2->cd(1); |
| 242 | TH1D *hMCptPi1D = data->ShowProjection(iptPi, stepGen); |
| 243 | Double_t maxptPi = hMCptPi1D->GetMaximum(); |
| 244 | hMCptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2); |
| 245 | hMCptPi1D->Sumw2(); |
| 246 | hMCptPi1D->SetMinimum(0.01); |
| 247 | hMCptPi1D->SetLineColor(2); |
| 248 | hMCptPi1D->SetLineWidth(3); |
| 249 | hMCptPi1D->SetMarkerColor(2); |
| 250 | hMCptPi1D->SetFillColor(2); |
| 251 | hMCptPi1D->SetFillStyle(3005); |
| 252 | hMCptPi1D->SetMarkerStyle(20); |
| 253 | hMCptPi1D->GetXaxis()->SetTitleOffset(1.2); |
| 254 | hMCptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), MC data"); |
| 255 | hMCptPi1D->Draw("hist"); |
| 256 | cmcptPi->cd(); |
| 257 | cmcptPi->SetLeftMargin(0.15); |
| 258 | cmcptPi->SetRightMargin(0.05); |
| 259 | hMCptPi1D->Draw("hist"); |
| 260 | hMCptPi1D->Draw("err same"); |
| 261 | cmcptPi->Update(); |
| 262 | |
| 263 | cmc2->cd(2); |
| 264 | TH1D *hRECptPi1D = data->ShowProjection(iptPi, stepRec); |
| 265 | hRECptPi1D->GetYaxis()->SetRangeUser(0,maxptPi*1.2); |
| 266 | hRECptPi1D->SetLineColor(4); |
| 267 | hRECptPi1D->SetLineWidth(3); |
| 268 | hRECptPi1D->SetMarkerColor(4); |
| 269 | hRECptPi1D->SetFillColor(4); |
| 270 | hRECptPi1D->SetFillStyle(3004); |
| 271 | hRECptPi1D->SetMarkerStyle(20); |
| 272 | hRECptPi1D->GetXaxis()->SetTitleOffset(1.2); |
| 273 | hRECptPi1D->Sumw2(); |
| 274 | hRECptPi1D->SetMinimum(0.01); |
| 275 | hRECptPi1D->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c), AOD"); |
| 276 | hRECptPi1D->Draw("hist"); |
| 277 | cptPi->cd(); |
| 278 | cptPi->SetLeftMargin(0.15); |
| 279 | cptPi->SetRightMargin(0.05); |
| 280 | hRECptPi1D->Draw("hist"); |
| 281 | hRECptPi1D->Draw("err same"); |
| 282 | cptPi->Update(); |
| 283 | |
| 284 | cmc2->cd(3); |
| 285 | TH1D *hMCptK1D = data->ShowProjection(iptK, stepGen); |
| 286 | Double_t maxptK = hMCptK1D->GetMaximum(); |
| 287 | hMCptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2); |
| 288 | hMCptK1D->SetLineColor(2); |
| 289 | hMCptK1D->SetLineWidth(3); |
| 290 | hMCptK1D->SetMarkerColor(2); |
| 291 | hMCptK1D->SetFillColor(2); |
| 292 | hMCptK1D->SetFillStyle(3005); |
| 293 | hMCptK1D->SetMarkerStyle(20); |
| 294 | hMCptK1D->GetXaxis()->SetTitleOffset(1.2); |
| 295 | hMCptK1D->Sumw2(); |
| 296 | hMCptK1D->SetMinimum(0.01); |
| 297 | hMCptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), MC data"); |
| 298 | hMCptK1D->Draw("hist"); |
| 299 | cmcptK->cd(); |
| 300 | cmcptK->SetLeftMargin(0.15); |
| 301 | cmcptK->SetRightMargin(0.05); |
| 302 | hMCptK1D->Draw("hist"); |
| 303 | hMCptK1D->Draw("err same"); |
| 304 | cmcptK->Update(); |
| 305 | |
| 306 | cmc2->cd(4); |
| 307 | TH1D *hRECptK1D = data->ShowProjection(iptK, stepRec); |
| 308 | hRECptK1D->GetYaxis()->SetRangeUser(0,maxptK*1.2); |
| 309 | hRECptK1D->SetLineColor(4); |
| 310 | hRECptK1D->SetLineWidth(3); |
| 311 | hRECptK1D->SetMarkerColor(4); |
| 312 | hRECptK1D->SetFillColor(4); |
| 313 | hRECptK1D->SetFillStyle(3004); |
| 314 | hRECptK1D->SetMarkerStyle(20); |
| 315 | hRECptK1D->GetXaxis()->SetTitleOffset(1.2); |
| 316 | hRECptK1D->Sumw2(); |
| 317 | hRECptK1D->SetMinimum(0.01); |
| 318 | hRECptK1D->Draw("hist"); |
| 319 | cptK->cd(); |
| 320 | cptK->SetLeftMargin(0.15); |
| 321 | cptK->SetRightMargin(0.05); |
| 322 | hRECptK1D->GetXaxis()->SetTitle("p_{T, K} (GeV/c), AOD"); |
| 323 | hRECptK1D->Draw("hist"); |
| 324 | hRECptK1D->Draw("err same"); |
| 325 | cptK->Update(); |
| 326 | |
| 327 | cmc2->cd(5); |
| 328 | TH1D *hMCcT1D = data->ShowProjection(icT, stepGen); |
| 329 | Double_t maxcT = hMCcT1D->GetMaximum(); |
| 330 | hMCcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2); |
| 331 | hMCcT1D->SetLineColor(2); |
| 332 | hMCcT1D->SetLineWidth(3); |
| 333 | hMCcT1D->SetMarkerColor(2); |
| 334 | hMCcT1D->SetFillColor(2); |
| 335 | hMCcT1D->SetFillStyle(3005); |
| 336 | hMCcT1D->SetMarkerStyle(20); |
| 337 | hMCcT1D->GetXaxis()->SetTitleOffset(1.2); |
| 338 | hMCcT1D->Sumw2(); |
| 339 | hMCcT1D->SetMinimum(0.01); |
| 340 | hMCcT1D->GetXaxis()->SetTitle("ct (#mum), MC data"); |
| 341 | hMCcT1D->Draw("hist"); |
| 342 | cmccT->cd(); |
| 343 | cmccT->SetLeftMargin(0.15); |
| 344 | cmccT->SetRightMargin(0.05); |
| 345 | hMCcT1D->Draw("hist"); |
| 346 | hMCcT1D->Draw("err same"); |
| 347 | cmccT->Update(); |
| 348 | |
| 349 | cmc2->cd(6); |
| 350 | TH1D *hRECcT1D = data->ShowProjection(icT, stepRec); |
| 351 | hRECcT1D->GetYaxis()->SetRangeUser(0,maxcT*1.2); |
| 352 | hRECcT1D->SetLineColor(4); |
| 353 | hRECcT1D->SetLineWidth(3); |
| 354 | hRECcT1D->SetMarkerColor(4); |
| 355 | hRECcT1D->SetFillColor(4); |
| 356 | hRECcT1D->SetFillStyle(3004); |
| 357 | hRECcT1D->SetMarkerStyle(20); |
| 358 | hRECcT1D->GetXaxis()->SetTitleOffset(1.2); |
| 359 | hRECcT1D->Sumw2(); |
| 360 | hRECcT1D->SetMinimum(0.01); |
| 361 | hRECcT1D->Draw("hist"); |
| 362 | ccT->cd(); |
| 363 | ccT->SetLeftMargin(0.15); |
| 364 | ccT->SetRightMargin(0.05); |
| 365 | hRECcT1D->GetXaxis()->SetTitle("c#t (#mum), AOD"); |
| 366 | hRECcT1D->Draw("hist"); |
| 367 | hRECcT1D->Draw("err same"); |
| 368 | ccT->Update(); |
| 369 | |
| 370 | /* |
| 371 | // printing on eps files |
| 372 | cmc1->Print("Plots/dataMC_pt_y_cTS.gif"); |
| 373 | cmc2->Print("Plots/dataMC_ptPi_ptK_cT.gif"); |
| 374 | cmcpt->Print("Plots/pt_Gen.eps"); |
| 375 | cmcy->Print("Plots/y_Gen.eps"); |
| 376 | cmccTS->Print("Plots/cTS_Gen.eps"); |
| 377 | cmcptPi->Print("Plots/ptPi_Gen.eps"); |
| 378 | cmcptK->Print("Plots/ptK_Gen.eps"); |
| 379 | cmccT->Print("Plots/cT_Gen.eps"); |
| 380 | cpt->Print("Plots/pt_Rec.eps"); |
| 381 | cy->Print("Plots/y_Rec.eps"); |
| 382 | ccTS->Print("Plots/cTS_Rec.eps"); |
| 383 | cptPi->Print("Plots/ptPi_Rec.eps"); |
| 384 | cptK->Print("Plots/ptK_Rec.eps"); |
| 385 | ccT->Print("Plots/cT_Rec.eps"); |
| 386 | cGen2d->Print("Plots/pt_y_Gen_2D.eps"); |
| 387 | cRec2d->Print("Plots/pt_y_Rec_2D.eps"); |
| 388 | |
| 389 | // printing on gif files |
| 390 | cmc1->Print("Plots/dataMC_pt_y_cTS.gif"); |
| 391 | cmc2->Print("Plots/dataMC_ptPi_ptK_cT.gif"); |
| 392 | cmcpt->Print("Plots/pt_Gen.eps"); |
| 393 | cmcy->Print("Plots/y_Gen.eps"); |
| 394 | cmccTS->Print("Plots/cTS_Gen.eps"); |
| 395 | cmcptPi->Print("Plots/ptPi_Gen.eps"); |
| 396 | cmcptK->Print("Plots/ptK_Gen.eps"); |
| 397 | cmccT->Print("Plots/cT_Gen.eps"); |
| 398 | cpt->Print("Plots/pt_Rec.eps"); |
| 399 | cy->Print("Plots/y_Rec.eps"); |
| 400 | ccTS->Print("Plots/cTS_Rec.eps"); |
| 401 | cptPi->Print("Plots/ptPi_Rec.eps"); |
| 402 | cptK->Print("Plots/ptK_Rec.eps"); |
| 403 | ccT->Print("Plots/cT_Rec.eps"); |
| 404 | cGen2d->Print("Plots/pt_y_Gen_2D.eps"); |
| 405 | cRec2d->Print("Plots/pt_y_Rec_2D.eps"); |
| 406 | */ |
| 407 | |
| 408 | //construct the efficiency grid from the data container |
| 409 | AliCFEffGrid *eff = new AliCFEffGrid("eff"," The efficiency",*data); |
| 410 | eff->CalculateEfficiency(stepRec,stepGen); //eff= step1/step0 |
| 411 | |
| 412 | //The efficiency along the variables, and some 2-D projections |
| 413 | TCanvas *ceff =new TCanvas("ceff"," Efficiency",0,0,1600,1200); |
| 414 | ceff->Divide(3,2); |
| 415 | TCanvas *ceff2D =new TCanvas("ceff2D"," Efficiency for pt and y",50,50,550,550); |
| 416 | TCanvas *ceffpt = new TCanvas("ceffpt","Efficiency vs pt",50,50,550,550); |
| 417 | TCanvas *ceffy = new TCanvas("ceffy","Efficiency vs y",50,50,550,550); |
| 418 | TCanvas *ceffcTS = new TCanvas("ceffcTS","Efficiency vs cosThetaStar",50,50,550,550); |
| 419 | TCanvas *ceffptPi = new TCanvas("ceffptPi","Efficiency vs ptPi",50,50,550,550); |
| 420 | TCanvas *ceffptK = new TCanvas("ceffptK","Efficiency vs ptK",50,50,550,550); |
| 421 | TCanvas *ceffcT = new TCanvas("ceffcT","Efficiency vs cT",50,50,550,550); |
| 422 | TCanvas *ceff2Dtext = new TCanvas("ceff2Dtext","Text plot for efficiency in pt and y",50,50,550,550); |
| 423 | |
| 424 | ceff->cd(1); |
| 425 | TH1D *hpteffCF = eff->Project(ipt); //the efficiency vs pt |
| 426 | hpteffCF->Sumw2(); |
| 427 | //hpteffCF->SetMinimum(0.01); |
| 428 | hpteffCF->SetLineColor(8); |
| 429 | hpteffCF->SetLineWidth(3); |
| 430 | hpteffCF->SetMarkerColor(8); |
| 431 | hpteffCF->SetMarkerStyle(20); |
| 432 | hpteffCF->GetXaxis()->SetTitleOffset(1.2); |
| 433 | hpteffCF->GetYaxis()->SetTitleOffset(1.5); |
| 434 | hpteffCF->Draw("hist"); |
| 435 | ceffpt->cd(); |
| 436 | ceffpt->SetLeftMargin(0.15); |
| 437 | ceffpt->SetRightMargin(0.05); |
| 438 | hpteffCF->GetXaxis()->SetTitle("p_{T} (GeV/c)"); |
| 439 | hpteffCF->GetYaxis()->SetTitle("Efficiency"); |
| 440 | hpteffCF->Draw("hist"); |
| 441 | hpteffCF->Draw("err same"); |
| 442 | |
| 443 | ceff->cd(2); |
| 444 | TH1D *hyeffCF = eff->Project(iy); //the efficiency vs y |
| 445 | hyeffCF->Sumw2(); |
| 446 | //hyeffCF->SetMinimum(0.01); |
| 447 | hyeffCF->SetLineColor(8); |
| 448 | hyeffCF->SetLineWidth(3); |
| 449 | hyeffCF->SetMarkerColor(8); |
| 450 | hyeffCF->SetMarkerStyle(20); |
| 451 | hyeffCF->GetXaxis()->SetTitleOffset(1.2); |
| 452 | hyeffCF->GetYaxis()->SetTitleOffset(1.5); |
| 453 | hyeffCF->Draw("hist"); |
| 454 | ceffy->cd(); |
| 455 | ceffy->SetLeftMargin(0.15); |
| 456 | ceffy->SetRightMargin(0.05); |
| 457 | hyeffCF->GetXaxis()->SetTitle("y"); |
| 458 | hyeffCF->GetYaxis()->SetTitle("Efficiency"); |
| 459 | hyeffCF->Draw("hist"); |
| 460 | hyeffCF->Draw("err same"); |
| 461 | |
| 462 | ceff->cd(3); |
| 463 | TH1D *hcTSeffCF = eff->Project(icTS); //the efficiency vs cosThetaStar |
| 464 | hcTSeffCF->Sumw2(); |
| 465 | //hcTSeffCF->SetMinimum(0.01); |
| 466 | hcTSeffCF->SetLineColor(8); |
| 467 | hcTSeffCF->SetLineWidth(3); |
| 468 | hcTSeffCF->SetMarkerColor(8); |
| 469 | hcTSeffCF->SetMarkerStyle(20); |
| 470 | hcTSeffCF->GetXaxis()->SetTitleOffset(1.2); |
| 471 | hcTSeffCF->GetYaxis()->SetTitleOffset(1.5); |
| 472 | hcTSeffCF->Draw("hist"); |
| 473 | ceffcTS->cd(); |
| 474 | ceffcTS->SetLeftMargin(0.15); |
| 475 | ceffcTS->SetRightMargin(0.05); |
| 476 | hcTSeffCF->GetXaxis()->SetTitle("cosThetaStar"); |
| 477 | hcTSeffCF->GetYaxis()->SetTitle("Efficiency"); |
| 478 | hcTSeffCF->Draw("hist"); |
| 479 | hcTSeffCF->Draw("err same"); |
| 480 | |
| 481 | ceff->cd(4); |
| 482 | TH1D *hptPieffCF = eff->Project(iptPi); //the efficiency vs ptPi |
| 483 | hptPieffCF->Sumw2(); |
| 484 | //hptPieffCF->SetMinimum(0.01); |
| 485 | hptPieffCF->SetLineColor(8); |
| 486 | hptPieffCF->SetLineWidth(3); |
| 487 | hptPieffCF->SetMarkerColor(8); |
| 488 | hptPieffCF->SetMarkerStyle(20); |
| 489 | hptPieffCF->GetXaxis()->SetTitleOffset(1.2); |
| 490 | hptPieffCF->GetYaxis()->SetTitleOffset(1.5); |
| 491 | hptPieffCF->Draw("hist"); |
| 492 | ceffptPi->cd(); |
| 493 | ceffptPi->SetLeftMargin(0.15); |
| 494 | ceffptPi->SetRightMargin(0.05); |
| 495 | hptPieffCF->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c)"); |
| 496 | hptPieffCF->GetYaxis()->SetTitle("Efficiency"); |
| 497 | hptPieffCF->Draw("hist"); |
| 498 | hptPieffCF->Draw("err same"); |
| 499 | |
| 500 | ceff->cd(5); |
| 501 | TH1D *hptKeffCF = eff->Project(iptK); //the efficiency vs ptK |
| 502 | hptKeffCF->Sumw2(); |
| 503 | //hptKeffCF->SetMinimum(0.01); |
| 504 | hptKeffCF->SetLineColor(8); |
| 505 | hptKeffCF->SetLineWidth(3); |
| 506 | hptKeffCF->SetMarkerColor(8); |
| 507 | hptKeffCF->SetMarkerStyle(20); |
| 508 | hptKeffCF->GetXaxis()->SetTitleOffset(1.2); |
| 509 | hptKeffCF->GetYaxis()->SetTitleOffset(1.5); |
| 510 | hptKeffCF->Draw("hist"); |
| 511 | ceffptK->cd(); |
| 512 | ceffptK->SetLeftMargin(0.15); |
| 513 | ceffptK->SetRightMargin(0.05); |
| 514 | hptKeffCF->GetXaxis()->SetTitle("p_{T, K} (GeV/c)"); |
| 515 | hptKeffCF->GetYaxis()->SetTitle("Efficiency"); |
| 516 | hptKeffCF->Draw("hist"); |
| 517 | hptKeffCF->Draw("err same"); |
| 518 | |
| 519 | ceff->cd(6); |
| 520 | TH1D *hcTeffCF = eff->Project(icT); //the efficiency vs cT |
| 521 | hcTeffCF->Sumw2(); |
| 522 | //hcTeffCF->SetMinimum(0.01); |
| 523 | hcTeffCF->SetLineColor(8); |
| 524 | hcTeffCF->SetLineWidth(3); |
| 525 | hcTeffCF->SetMarkerColor(8); |
| 526 | hcTeffCF->SetMarkerStyle(20); |
| 527 | hcTeffCF->GetXaxis()->SetTitleOffset(1.2); |
| 528 | hcTeffCF->GetYaxis()->SetTitleOffset(1.5); |
| 529 | hcTeffCF->Draw("hist"); |
| 530 | ceffcT->cd(); |
| 531 | ceffcT->SetLeftMargin(0.15); |
| 532 | ceffcT->SetRightMargin(0.05); |
| 533 | hcTeffCF->GetXaxis()->SetTitle("c#t (#mum)"); |
| 534 | hcTeffCF->GetYaxis()->SetTitle("Efficiency"); |
| 535 | hcTeffCF->Draw("hist"); |
| 536 | hcTeffCF->Draw("err same"); |
| 537 | |
| 538 | ceff2D->cd(); |
| 539 | TH2D *hptyeffCF = eff->Project(ipt,iy); //look at the numerator |
| 540 | //hptyeffCF->SetMinimum(0.01); |
| 541 | hptyeffCF->SetMarkerColor(8); |
| 542 | hptyeffCF->SetLineColor(8); |
| 543 | hptyeffCF->SetMinimum(0.01); |
| 544 | hptyeffCF->Draw("lego"); |
| 545 | ceff2Dtext->cd(); |
| 546 | hptyeffCF->GetXaxis()->SetTitle("p_{T} (GeV/c)"); |
| 547 | hptyeffCF->GetXaxis()->SetTitleOffset(1.2); |
| 548 | hptyeffCF->GetYaxis()->SetTitle("y"); |
| 549 | hptyeffCF->GetYaxis()->SetTitleOffset(1.2); |
| 550 | hptyeffCF->Draw("text"); |
| 551 | |
| 552 | /* |
| 553 | // printing eps files |
| 554 | ceff->Print("Plots/efficiencies.eps"); |
| 555 | ceffpt->Print("Plots/effpt.eps"); |
| 556 | ceffy->Print("Plots/effy.eps"); |
| 557 | ceffcTS->Print("Plots/effcTS.eps"); |
| 558 | ceffptPi->Print("Plots/effptPi.eps"); |
| 559 | ceffptK->Print("Plots/effptK.eps"); |
| 560 | ceffcT->Print("Plots/effcT.eps"); |
| 561 | ceff2D->Print("Plots/eff2D_pt_y.eps"); |
| 562 | ceff2Dtext->Print("Plots/eff2d_text_pt_y.eps"); |
| 563 | |
| 564 | // printing gif files |
| 565 | ceff->Print("Plots/efficiencies.gif"); |
| 566 | ceffpt->Print("Plots/effpt.gif"); |
| 567 | ceffy->Print("Plots/effy.gif"); |
| 568 | ceffcTS->Print("Plots/effcTS.gif"); |
| 569 | ceffptPi->Print("Plots/effptPi.gif"); |
| 570 | ceffptK->Print("Plots/effptK.gif"); |
| 571 | ceffcT->Print("Plots/effcT.gif"); |
| 572 | ceff2D->Print("Plots/eff2D_pt_y.gif"); |
| 573 | ceff2Dtext->Print("Plots/eff2d_text_pt_y.gif"); |
| 574 | */ |
| 575 | |
| 576 | // applying efficiencies - using projections, not the ApplyEffCorrection method |
| 577 | |
| 578 | TCanvas *cmultpiplypt = new TCanvas("cmultiplypt","Reco From Eff in pt distribution",50,50,550,550); |
| 579 | TCanvas *cmultpiplyy = new TCanvas("cmultiplyy","Reco From Eff in y distribution",50,50,550,550); |
| 580 | TCanvas *cmultpiplycTS = new TCanvas("cmultiplycTS","Reco From Eff in cTS distribution",50,50,550,550); |
| 581 | TCanvas *cmultpiplyptPi = new TCanvas("cmultiplyptPi","Reco From Eff in ptPi distribution",50,50,550,550); |
| 582 | TCanvas *cmultpiplyptK = new TCanvas("cmultiplyptK","Reco From Eff in ptK distribution",50,50,550,550); |
| 583 | TCanvas *cmultpiplycT = new TCanvas("cmultiplycT","Reco From Eff in cT distribution",50,50,550,550); |
| 584 | |
| 585 | TH1D *hmultiplypt = new TH1D("hmultiplypt","hmultiplypt",13,0,10); |
| 586 | cout << " bin for histo MC = " << hMCpt1D->GetNbinsX() << " while for RECO histo = " << hRECpt1D->GetNbinsX() << " while for efficiency histo = " << hpteffCF->GetNbinsX() << endl; |
| 587 | |
| 588 | const Double_t ptmin_0_4 = 0.0 ; |
| 589 | const Double_t ptmax_0_4 = 4.0 ; |
| 590 | const Double_t ptmin_4_8 = 4.0 ; |
| 591 | const Double_t ptmax_4_8 = 8.0 ; |
| 592 | const Double_t ptmin_8_10 = 8.0 ; |
| 593 | const Double_t ptmax_8_10 = 10.0 ; |
| 594 | const Int_t nbin0_0_4 = 8 ; //bins in pt |
| 595 | const Int_t nbin0_4_8 = 4 ; //bins in pt |
| 596 | const Int_t nbin0_8_10 = 1 ; //bins in pt |
| 597 | const Int_t nbins = nbin0_0_4 + nbin0_4_8 + nbin0_8_10; |
| 598 | Double_t binLim0[nbins+1]; |
| 599 | for(Int_t i=0; i<=nbin0_0_4; i++) binLim0[i]=(Double_t)ptmin_0_4 + (ptmax_0_4-ptmin_0_4)/nbin0_0_4*(Double_t)i ; |
| 600 | if (binLim0[nbin0_0_4] != ptmin_4_8) { |
| 601 | printf("Calculated bin lim for pt - 1st range - differs from expected!\n"); |
| 602 | } |
| 603 | for(Int_t i=0; i<=nbin0_4_8; i++) binLim0[i+nbin0_0_4]=(Double_t)ptmin_4_8 + (ptmax_4_8-ptmin_4_8)/nbin0_4_8*(Double_t)i ; |
| 604 | if (binLim0[nbin0_0_4+nbin0_4_8] != ptmin_8_10) { |
| 605 | printf("Calculated bin lim for pt - 2nd range - differs from expected!\n"); |
| 606 | } |
| 607 | for(Int_t i=0; i<=nbin0_8_10; i++) binLim0[i+nbin0_0_4+nbin0_4_8]=(Double_t)ptmin_8_10 + (ptmax_8_10-ptmin_8_10)/nbin0_8_10*(Double_t)i ; |
| 608 | |
| 609 | hmultiplypt->SetBins(nbins,binLim0); |
| 610 | hmultiplypt->Divide(hRECpt1D,hpteffCF,1,1,"b"); |
| 611 | hmultiplypt->SetMinimum(0); |
| 612 | hmultiplypt->SetMaximum(hMCpt1D->GetMaximum()); |
| 613 | hmultiplypt->SetLineColor(38); |
| 614 | hmultiplypt->SetLineWidth(3); |
| 615 | hmultiplypt->SetMarkerColor(38); |
| 616 | hmultiplypt->SetMarkerStyle(20); |
| 617 | hmultiplypt->GetXaxis()->SetTitleOffset(1.2); |
| 618 | hmultiplypt->GetYaxis()->SetTitleOffset(1.5); |
| 619 | hmultiplypt->GetXaxis()->SetTitle("p_{T} (GeV/c)"); |
| 620 | cmultiplypt->SetLeftMargin(0.15); |
| 621 | cmultiplypt->SetRightMargin(0.05); |
| 622 | hmultiplypt->SetFillStyle(3004); |
| 623 | hmultiplypt->SetFillColor(38); |
| 624 | cmultiplypt->cd(); |
| 625 | hmultiplypt->Draw("hist err"); |
| 626 | hMCpt1D->Draw("hist same err"); |
| 627 | |
| 628 | TH1D *hmultiplyy = new TH1D("hmultiplyy","hmultiplyy",8,-2,2); |
| 629 | hmultiplyy->Divide(hRECy1D,hyeffCF,1,1,"b"); |
| 630 | hmultiplyy->SetMinimum(0); |
| 631 | hmultiplyy->SetMaximum(hMCy1D->GetMaximum()); |
| 632 | hmultiplyy->SetLineColor(38); |
| 633 | hmultiplyy->SetLineWidth(3); |
| 634 | hmultiplyy->SetMarkerColor(38); |
| 635 | hmultiplyy->SetMarkerStyle(20); |
| 636 | hmultiplyy->GetXaxis()->SetTitleOffset(1.2); |
| 637 | hmultiplyy->GetXaxis()->SetTitle("y"); |
| 638 | hmultiplyy->GetYaxis()->SetTitleOffset(1.5); |
| 639 | cmultiplyy->SetLeftMargin(0.15); |
| 640 | cmultiplyy->SetRightMargin(0.05); |
| 641 | cmultiplyy->cd(); |
| 642 | hmultiplyy->SetFillStyle(3004); |
| 643 | hmultiplyy->SetFillColor(38); |
| 644 | hmultiplyy->Draw("hist err"); |
| 645 | hMCy1D->Draw("hist same err"); |
| 646 | |
| 647 | TH1D *hmultiplycTS = new TH1D("hmultiplycTS","hmultiplycTS",8,-1,1); |
| 648 | hmultiplycTS->Divide(hRECcTS1D,hcTSeffCF,1,1,"b"); |
| 649 | hmultiplycTS->SetMinimum(0); |
| 650 | hmultiplycTS->SetMaximum(hMCcTS1D->GetMaximum()); |
| 651 | hmultiplycTS->SetLineColor(38); |
| 652 | hmultiplycTS->SetLineWidth(3); |
| 653 | hmultiplycTS->SetMarkerColor(38); |
| 654 | hmultiplycTS->SetMarkerStyle(20); |
| 655 | hmultiplycTS->GetXaxis()->SetTitleOffset(1.2); |
| 656 | hmultiplycTS->GetXaxis()->SetTitle("cosThetaStar"); |
| 657 | hmultiplycTS->GetYaxis()->SetTitleOffset(1.5); |
| 658 | cmultiplycTS->SetLeftMargin(0.15); |
| 659 | cmultiplycTS->SetRightMargin(0.05); |
| 660 | cmultiplycTS->cd(); |
| 661 | hmultiplycTS->SetFillStyle(3004); |
| 662 | hmultiplycTS->SetFillColor(38); |
| 663 | hmultiplycTS->Draw("hist err"); |
| 664 | hMCcTS1D->Draw("hist same err"); |
| 665 | |
| 666 | TH1D *hmultiplyptPi = new TH1D("hmultiplyptPi","hmultiplyptPi",13,0,10); |
| 667 | hmultiplyptPi->SetBins(nbins,binLim0); |
| 668 | hmultiplyptPi->Divide(hRECptPi1D,hptPieffCF,1,1,"b"); |
| 669 | hmultiplyptPi->SetMinimum(0); |
| 670 | hmultiplyptPi->SetMaximum(hMCptPi1D->GetMaximum()); |
| 671 | hmultiplyptPi->SetLineColor(38); |
| 672 | hmultiplyptPi->SetLineWidth(3); |
| 673 | hmultiplyptPi->SetMarkerColor(38); |
| 674 | hmultiplyptPi->SetMarkerStyle(20); |
| 675 | hmultiplyptPi->GetXaxis()->SetTitleOffset(1.2); |
| 676 | hmultiplyptPi->GetXaxis()->SetTitle("p_{T, #pi} (GeV/c)"); |
| 677 | hmultiplyptPi->GetYaxis()->SetTitleOffset(1.5); |
| 678 | cmultiplyptPi->SetLeftMargin(0.15); |
| 679 | cmultiplyptPi->SetRightMargin(0.05); |
| 680 | cmultiplyptPi->cd(); |
| 681 | hmultiplyptPi->SetFillStyle(3004); |
| 682 | hmultiplyptPi->SetFillColor(38); |
| 683 | hmultiplyptPi->Draw("hist err"); |
| 684 | hMCptPi1D->Draw("hist same err"); |
| 685 | |
| 686 | TH1D *hmultiplyptK = new TH1D("hmultiplyptK","hmultiplyptK",13,0,10); |
| 687 | hmultiplyptK->SetBins(nbins,binLim0); |
| 688 | hmultiplyptK->Divide(hRECptK1D,hptKeffCF,1,1,"b"); |
| 689 | hmultiplyptK->SetMinimum(0); |
| 690 | hmultiplyptK->SetMaximum(hMCptK1D->GetMaximum()); |
| 691 | hmultiplyptK->SetLineColor(38); |
| 692 | hmultiplyptK->SetLineWidth(3); |
| 693 | hmultiplyptK->SetMarkerColor(38); |
| 694 | hmultiplyptK->SetMarkerStyle(20); |
| 695 | hmultiplyptK->GetXaxis()->SetTitleOffset(1.2); |
| 696 | hmultiplyptK->GetXaxis()->SetTitle("p_{T, K} (GeV/c)"); |
| 697 | hmultiplyptK->GetYaxis()->SetTitleOffset(1.5); |
| 698 | cmultiplyptK->SetLeftMargin(0.15); |
| 699 | cmultiplyptK->SetRightMargin(0.05); |
| 700 | cmultiplyptK->cd(); |
| 701 | hmultiplyptK->SetFillStyle(3004); |
| 702 | hmultiplyptK->SetFillColor(38); |
| 703 | hmultiplyptK->Draw("hist err"); |
| 704 | hMCptK1D->Draw("hist same err"); |
| 705 | |
| 706 | TH1D *hmultiplycT = new TH1D("hmultiplycT","hmultiplycT",24,0,500); |
| 707 | hmultiplycT->Divide(hRECcT1D,hcTeffCF,1,1,"b"); |
| 708 | hmultiplycT->SetMinimum(0); |
| 709 | hmultiplycT->SetMaximum(hMCcT1D->GetMaximum()); |
| 710 | hmultiplycT->SetLineColor(38); |
| 711 | hmultiplycT->SetLineWidth(3); |
| 712 | hmultiplycT->SetMarkerColor(38); |
| 713 | hmultiplycT->SetMarkerStyle(20); |
| 714 | hmultiplycT->GetXaxis()->SetTitleOffset(1.2); |
| 715 | hmultiplycT->GetXaxis()->SetTitle("c#t (#mum)"); |
| 716 | hmultiplycT->GetYaxis()->SetTitleOffset(1.5); |
| 717 | cmultiplycT->SetLeftMargin(0.15); |
| 718 | cmultiplycT->SetRightMargin(0.05); |
| 719 | cmultiplycT->cd(); |
| 720 | hmultiplycT->SetFillStyle(3004); |
| 721 | hmultiplycT->SetFillColor(38); |
| 722 | hmultiplycT->Draw("hist err"); |
| 723 | hMCcT1D->Draw("hist same err"); |
| 724 | |
| 725 | /* |
| 726 | TFile* file_histo = new TFile("fileHisto_180004.root","RECREATE"); |
| 727 | hMCpt1D->Write("hMCpt1D"); |
| 728 | hRECpt1D->Write("hRECpt1D"); |
| 729 | hMCy1D->Write("hMCy1D"); |
| 730 | hRECy1D->Write("hRECy1D"); |
| 731 | hMCcTS1D->Write("hMCcTS1D"); |
| 732 | hRECcTS1D->Write("hRECcTS1D"); |
| 733 | hMCptPi1D->Write("hMCptPi1D"); |
| 734 | hRECptPi1D->Write("hRECptPi1D"); |
| 735 | hMCptK1D->Write("hMCptK1D"); |
| 736 | hRECptK1D->Write("hRECptK1D"); |
| 737 | hMCcT1D->Write("hMCcT1D"); |
| 738 | hRECcT1D->Write("hRECcT1D"); |
| 739 | hpteffCF->Write("hpteffCF"); |
| 740 | hyeffCF->Write("hyeffCF"); |
| 741 | hcTSeffCF->Write("hcTSeffCF"); |
| 742 | hptPieffCF->Write("hptPieffCF"); |
| 743 | hptKeffCF->Write("hptKeffCF"); |
| 744 | hcTeffCF->Write("hcTeffCF"); |
| 745 | hmultiplypt->Write(""); |
| 746 | hmultiplyy->Write(""); |
| 747 | hmultiplycTS->Write(""); |
| 748 | hmultiplyptPi->Write(""); |
| 749 | hmultiplyptK->Write(""); |
| 750 | hmultiplycT->Write(""); |
| 751 | file_histo->Close(); |
| 752 | */ |
| 753 | |
| 754 | Int_t nsliceVars = 5; |
| 755 | Int_t sliceVars[5]; |
| 756 | sliceVars[0] = 1; |
| 757 | sliceVars[1] = 2; |
| 758 | sliceVars[2] = 3; |
| 759 | sliceVars[3] = 4; |
| 760 | sliceVars[4] = 5; |
| 761 | Double_t varMin[6] = {0.5, -2, -1, 0., 0., 0. }; |
| 762 | Double_t varMax[6] = {1.5, 2, 1, 10., 10., 500.}; |
| 763 | AliCFContainer* data_sliced = data->MakeSlice(nsliceVars, sliceVars, varMin, varMax); |
| 764 | cout << " the container now has " << data_sliced->GetNVar() << " dimensions " << endl; |
| 765 | |
| 766 | } |
git://git.uio.no / u / mrichter / AliRoot.git / blame