PWGPP/macros/plotAnalysisTaskITSTPCalignment.C

   1 void plotAnalysisTaskITSTPCalignment(const char* option = "b")
   2 {
   3   gStyle->SetCanvasBorderMode(0);
   4   gStyle->SetPadBorderMode(0);
   5   gStyle->SetFrameBorderMode(0);
   6   gStyle->SetCanvasColor(0);
   7   gStyle->SetHistFillColor(17);
   8
   9   TString optionstr(option);
  10   gROOT->LoadMacro("AliRelAlignerKalmanArray.cxx++");
  11   TFile f("outputITSTPCalignment.root");
  12
  13   //////////////////////////////////////////////////////////////////////////////////////////////////
  14   AliRelAlignerKalmanArray* fArray = dynamic_cast<AliRelAlignerKalmanArray*>(f.Get("outputArrayITSsa"));
  15   if (!fArray)
  16   {
  17     printf("fArray cannot be read!\n");
  18     return;
  19   }
  20   else
  21   {
  22     TCanvas* c1 = new TCanvas("c1","psi in time");
  23     fArray->MakeGraph(0)->Draw("A*");
  24     c1->SaveAs("graphpsi.eps");
  25
  26     TCanvas* c2 = new TCanvas("c2","theta in time");
  27     fArray->MakeGraph(1)->Draw("A*");
  28     c2->SaveAs("graphtheta.eps");
  29
  30     TCanvas* c3 = new TCanvas("c3","phi in time");
  31     fArray->MakeGraph(2)->Draw("A*");
  32     c3->SaveAs("graphphi.eps");
  33
  34     TCanvas* c4 = new TCanvas("c4","x in time");
  35     fArray->MakeGraph(3)->Draw("A*");
  36     c4->SaveAs("graphx.eps");
  37
  38     TCanvas* c5 = new TCanvas("c5","y in time");
  39     fArray->MakeGraph(4)->Draw("A*");
  40     c5->SaveAs("graphy.eps");
  41
  42     TCanvas* c6 = new TCanvas("c6","z in time");
  43     fArray->MakeGraph(5)->Draw("A*");
  44     c6->SaveAs("graphz.eps");
  45
  46     TCanvas* c7 = new TCanvas("c7","TPC vd correction in time");
  47     fArray->MakeGraph(6)->Draw("A*");
  48     c7->SaveAs("graphvd.eps");
  49
  50     TCanvas* c8 = new TCanvas("c8","TPC t0 correction in time");
  51     fArray->MakeGraph(7)->Draw("A*");
  52     c8->SaveAs("grapht0.eps");
  53
  54     TCanvas* c9 = new TCanvas("c9","TPC dv/dy in time");
  55     fArray->MakeGraph(8)->Draw("A*");
  56     c9->SaveAs("graphdvdy.eps");
  57   }
  58
  59   //////////////////////////////////////////////////////////////////////////////////////////////////
  60   TList* fList = dynamic_cast<TList*>(f.Get("outputList"));
  61   TList* pList = NULL;
  62
  63   TH1F* pMatchingEfficiency = dynamic_cast<TH1F*>(fList->At(3));
  64   TCanvas* canvasMEff = new TCanvas();
  65   pMatchingEfficiency->DrawCopy();
  66
  67
  68   //------------------------------------------------------------------------------------
  69   pList = dynamic_cast<TList*>(fList->At(0));
  70   TH2F* pZYAResidualsHistBpos = dynamic_cast<TH2F*>(pList->At(0));
  71   TH2F* pZYCResidualsHistBpos = dynamic_cast<TH2F*>(pList->At(1));
  72   TH2F* pLPAResidualsHistBpos = dynamic_cast<TH2F*>(pList->At(2));
  73   TH2F* pLPCResidualsHistBpos = dynamic_cast<TH2F*>(pList->At(3));
  74   TH2F* pPhiYAResidualsHistBpos = dynamic_cast<TH2F*>(pList->At(4));
  75   TH2F* pPhiYCResidualsHistBpos = dynamic_cast<TH2F*>(pList->At(5));
  76   TH2F* pPhiZAResidualsHistBpos = dynamic_cast<TH2F*>(pList->At(6));
  77   TH2F* pPhiZCResidualsHistBpos = dynamic_cast<TH2F*>(pList->At(7));
  78   TH2F* pPtYAResidualsHistBpos = dynamic_cast<TH2F*>(pList->At(8));
  79   TH2F* pPtYCResidualsHistBpos = dynamic_cast<TH2F*>(pList->At(9));
  80   TH2F* pPtZAResidualsHistBpos = dynamic_cast<TH2F*>(pList->At(10));
  81   TH2F* pPtZCResidualsHistBpos = dynamic_cast<TH2F*>(pList->At(11));
  82   TH2F* pLowPtYAResidualsHistBpos = dynamic_cast<TH2F*>(pList->At(12));
  83   TH2F* pLowPtYCResidualsHistBpos = dynamic_cast<TH2F*>(pList->At(13));
  84   TH2F* pLowPtZAResidualsHistBpos = dynamic_cast<TH2F*>(pList->At(14));
  85   TH2F* pLowPtZCResidualsHistBpos = dynamic_cast<TH2F*>(pList->At(15));
  86
  87   TH1D* resycBpos = pZYCResidualsHistBpos->ProjectionY();
  88   resycBpos->SetTitle("r\\phi residual distribution side C (z<0), Bpos");
  89   TH1D* reszcBpos = pZYCResidualsHistBpos->ProjectionX();
  90   reszcBpos->SetTitle("z residual distribution side C (z<0), Bpos");
  91   TH1D* respcBpos = pLPCResidualsHistBpos->ProjectionY();
  92   respcBpos->SetTitle("sin(\\phi) residual distribution side C (z<0), Bpos");
  93   TH1D* reslcBpos = pLPCResidualsHistBpos->ProjectionX();
  94   reslcBpos->SetTitle("tan(\\lambda) residual distribution side C (z<0), Bpos");
  95   TProfile* profphiycBpos = pPhiYCResidualsHistBpos->ProfileX("_pfx",1,-1,"s");
  96   profphiycBpos->SetTitle("\\phi profile of r\\phi residuals side C (z<0), Bpos");
  97   profphiycBpos->SetYTitle("\\deltar\\phi [cm]");
  98   TProfile* profphizcBpos = pPhiZCResidualsHistBpos->ProfileX("_pfx",1,-1,"s");
  99   profphizcBpos->SetTitle("\\phi profile of z residuals side C (z<0), Bpos");
 100   profphizcBpos->SetYTitle("\\deltaz [cm]");
 101   TProfile* profptycBpos = pPtYCResidualsHistBpos->ProfileX("_pfx",1,-1,"s");
 102   profptycBpos->SetTitle("pt profile of r\\phi residuals side C (z<0), Bpos");
 103   profptycBpos->SetYTitle("\\deltar\\phi [cm]");
 104   TProfile* profptzcBpos = pPtZCResidualsHistBpos->ProfileX("_pfx",1,-1,"s");
 105   profptzcBpos->SetTitle("pt profile of z residuals side C (z<0), Bpos");
 106   profptzcBpos->SetYTitle("\\deltaz [cm]");
 107   TProfile* proflowptycBpos = pLowPtYCResidualsHistBpos->ProfileX("_pfx",1,-1,"s");
 108   proflowptycBpos->SetTitle("pt profile of r\\phi residuals side C (z<0), Bpos");
 109   proflowptycBpos->SetYTitle("\\deltar\\phi [cm]");
 110   TProfile* proflowptzcBpos = pLowPtZCResidualsHistBpos->ProfileX("_pfx",1,-1,"s");
 111   proflowptzcBpos->SetTitle("pt profile of z residuals side C (z<0), Bpos");
 112   proflowptzcBpos->SetYTitle("\\deltaz [cm]");
 113
 114   TH1D* resyaBpos = pZYAResidualsHistBpos->ProjectionY();
 115   resyaBpos->SetTitle("r\\phi residual distribution side A (z>0), Bpos");
 116   TH1D* reszaBpos = pZYAResidualsHistBpos->ProjectionX();
 117   reszaBpos->SetTitle("z residual distribution side A (z>0), Bpos");
 118   TH1D* respaBpos = pLPAResidualsHistBpos->ProjectionY();
 119   respaBpos->SetTitle("sin(\\phi) residual distribution side A (z>0), Bpos");
 120   TH1D* reslaBpos = pLPAResidualsHistBpos->ProjectionX();
 121   reslaBpos->SetTitle("tan(\\lambda) residual distribution side A (z>0), Bpos");
 122   TProfile* profphiyaBpos = pPhiYAResidualsHistBpos->ProfileX("_pfx",1,-1,"s");
 123   profphiyaBpos->SetTitle("\\phi profile of r\\phi residuals side A (z>0), Bpos");
 124   profphiyaBpos->SetYTitle("\\deltar\\phi [cm]");
 125   TProfile* profphizaBpos = pPhiZAResidualsHistBpos->ProfileX("_pfx",1,-1,"s");
 126   profphizaBpos->SetTitle("\\phi profile of z residuals side A (z>0), Bpos");
 127   profphizaBpos->SetYTitle("\\deltaz [cm]");
 128   TProfile* profptyaBpos = pPtYAResidualsHistBpos->ProfileX("_pfx",1,-1,"s");
 129   profptyaBpos->SetTitle("pt profile of r\\phi residuals side A (z>0), Bpos");
 130   profptyaBpos->SetYTitle("\\deltar\\phi [cm]");
 131   TProfile* profptzaBpos = pPtZAResidualsHistBpos->ProfileX("_pfx",1,-1,"s");
 132   profptzaBpos->SetTitle("pt profile of z residuals side A (z>0), Bpos");
 133   profptzaBpos->SetYTitle("\\deltaz [cm]");
 134   TProfile* proflowptyaBpos = pLowPtYAResidualsHistBpos->ProfileX("_pfx",1,-1,"s");
 135   proflowptyaBpos->SetTitle("pt profile of r\\phi residuals side A (z>0), Bpos");
 136   proflowptyaBpos->SetYTitle("\\deltar\\phi [cm]");
 137   TProfile* proflowptzaBpos = pLowPtZAResidualsHistBpos->ProfileX("_pfx",1,-1,"s");
 138   proflowptzaBpos->SetTitle("pt profile of z residuals side A (z>0), Bpos");
 139   proflowptzaBpos->SetYTitle("\\deltaz [cm]");
 140
 141   if (pZYAResidualsHistBpos->GetEntries()>0)
 142   {
 143   TCanvas* c1Bpos = new TCanvas("c1Bpos","Residuals ZY 2D, Bpos",800,300);
 144   c1Bpos->Divide(2,1);
 145   c1Bpos->cd(1); pZYAResidualsHistBpos->DrawCopy("col");
 146   c1Bpos->cd(2); pZYCResidualsHistBpos->DrawCopy("col");
 147   c1Bpos->cd(0);
 148   c1Bpos->SaveAs("residualsZY2D-Bpos.eps");
 149
 150   TCanvas* c2Bpos = new TCanvas("c2Bpos","Residuals LP 2D, Bpos",800,300);
 151   c2Bpos->Divide(2,1);
 152   c2Bpos->cd(1); pLPAResidualsHistBpos->DrawCopy("col");
 153   c2Bpos->cd(2); pLPCResidualsHistBpos->DrawCopy("col");
 154   c2Bpos->cd(0);
 155   c2Bpos->SaveAs("residualsLP2D-Bpos.eps");
 156
 157   TCanvas* c3Bpos = new TCanvas("c3Bpos","Residuals z 1D, Bpos",800,300);
 158   c3Bpos->Divide(2,1);
 159   c3Bpos->cd(1); reszaBpos->DrawCopy();
 160   c3Bpos->cd(2); reszcBpos->DrawCopy();
 161   c3Bpos->cd(0);
 162   c3Bpos->SaveAs("residualsZ-Bpos.eps");
 163
 164   TCanvas* c4Bpos = new TCanvas("c4Bpos","Residuals y 1D, Bpos",800,300);
 165   c4Bpos->Divide(2,1);
 166   c4Bpos->cd(1); resyaBpos->DrawCopy();
 167   c4Bpos->cd(2); resycBpos->DrawCopy();
 168   c4Bpos->cd(0);
 169   c4Bpos->SaveAs("residualsY-Bpos.eps");
 170
 171   TCanvas* c5Bpos = new TCanvas("c5Bpos","Residuals phi 1D, Bpos",800,300);
 172   c5Bpos->Divide(2,1);
 173   c5Bpos->cd(1); respaBpos->DrawCopy();
 174   c5Bpos->cd(2); respcBpos->DrawCopy();
 175   c5Bpos->cd(0);
 176   c5Bpos->SaveAs("residualsPhi-Bpos.eps");
 177
 178   TCanvas* c6Bpos = new TCanvas("c6Bpos","Residuals lambda 1D, Bpos",800,300);
 179   c6Bpos->Divide(2,1);
 180   c6Bpos->cd(1); reslaBpos->DrawCopy();
 181   c6Bpos->cd(2); reslcBpos->DrawCopy();
 182   c6Bpos->cd(0);
 183   c6Bpos->SaveAs("residualsLambda-Bpos.eps");
 184
 185   TCanvas* c7Bpos = new TCanvas("c7Bpos","Profiles z in phi, Bpos",800,300);
 186   c7Bpos->Divide(2,1);
 187   c7Bpos->cd(1); profphizaBpos->DrawCopy();
 188   c7Bpos->cd(2); profphizcBpos->DrawCopy();
 189   c7Bpos->cd(0);
 190   c7Bpos->SaveAs("residualsProfilePhiZ-Bpos.eps");
 191
 192   TCanvas* c8Bpos = new TCanvas("c8Bpos","Profiles y in phi, Bpos",800,300);
 193   c8Bpos->Divide(2,1);
 194   c8Bpos->cd(1); profphiyaBpos->DrawCopy();
 195   c8Bpos->cd(2); profphiycBpos->DrawCopy();
 196   c8Bpos->cd(0);
 197   c8Bpos->SaveAs("residualsProfilePhiY-Bpos.eps");
 198
 199   TCanvas* c9Bpos = new TCanvas("c9Bpos", "Residuals Phi-Z 2D, Bpos",800,300);
 200   c9Bpos->Divide(2,1);
 201   c9Bpos->cd(1); pPhiZAResidualsHistBpos->DrawCopy("col");
 202   c9Bpos->cd(2); pPhiZCResidualsHistBpos->DrawCopy("col");
 203   c9Bpos->cd(0);
 204   c9Bpos->SaveAs("residualsPhiZ2D-Bpos.eps");
 205
 206   TCanvas* c1Bpos0 = new TCanvas("c1Bpos0", "Residuals Phi-Y 2D, Bpos",800,300);
 207   c1Bpos0->Divide(2,1);
 208   c1Bpos0->cd(1); pPhiYAResidualsHistBpos->DrawCopy("col");
 209   c1Bpos0->cd(2); pPhiYCResidualsHistBpos->DrawCopy("col");
 210   c1Bpos0->cd(0);
 211   c1Bpos0->SaveAs("residualsPhiY2D-Bpos.eps");
 212
 213   TCanvas* c11Bpos = new TCanvas("c11Bpos", "Residuals Pt-Z 2D, Bpos",800,300);
 214   c11Bpos->Divide(2,1);
 215   c11Bpos->cd(1); pPtZAResidualsHistBpos->DrawCopy("col");
 216   c11Bpos->cd(2); pPtZCResidualsHistBpos->DrawCopy("col");
 217   c11Bpos->cd(0);
 218   c11Bpos->SaveAs("residualsPtZ2D-Bpos.eps");
 219
 220   TCanvas* c12Bpos = new TCanvas("c12Bpos", "Residuals Pt-Y 2D, Bpos",800,300);
 221   c12Bpos->Divide(2,1);
 222   c12Bpos->cd(1); pPtYAResidualsHistBpos->DrawCopy("col");
 223   c12Bpos->cd(2); pPtYCResidualsHistBpos->DrawCopy("col");
 224   c12Bpos->cd(0);
 225   c12Bpos->SaveAs("residualsPtY2D-Bpos.eps");
 226
 227   TCanvas* c13Bpos = new TCanvas("c13Bpos","Profiles Pt-Z, Bpos",800,300);
 228   c13Bpos->Divide(2,1);
 229   c13Bpos->cd(1); profptzaBpos->DrawCopy();
 230   c13Bpos->cd(2); profptzcBpos->DrawCopy();
 231   c13Bpos->cd(0);
 232   c13Bpos->SaveAs("residualsProfilePtZ-Bpos.eps");
 233
 234   TCanvas* c14Bpos = new TCanvas("c14Bpos","Profiles Pt-Y, Bpos",800,300);
 235   c14Bpos->Divide(2,1);
 236   c14Bpos->cd(1); profptyaBpos->DrawCopy();
 237   c14Bpos->cd(2); profptycBpos->DrawCopy();
 238   c14Bpos->cd(0);
 239   c14Bpos->SaveAs("residualsProfilePtY-Bpos.eps");
 240
 241   TCanvas* c15Bpos = new TCanvas("c15Bpos", "Residuals low Pt-Z 2D, Bpos",800,300);
 242   c15Bpos->Divide(2,1);
 243   c15Bpos->cd(1); pLowPtZAResidualsHistBpos->DrawCopy("col");
 244   c15Bpos->cd(2); pLowPtZCResidualsHistBpos->DrawCopy("col");
 245   c15Bpos->cd(0);
 246   c15Bpos->SaveAs("residualsLowPtZ2D-Bpos.eps");
 247
 248   TCanvas* c16Bpos = new TCanvas("c16Bpos", "Residuals in low PtY 2D, Bpos",800,300);
 249   c16Bpos->Divide(2,1);
 250   c16Bpos->cd(1); pLowPtYAResidualsHistBpos->DrawCopy("col");
 251   c16Bpos->cd(2); pLowPtYCResidualsHistBpos->DrawCopy("col");
 252   c16Bpos->cd(0);
 253   c16Bpos->SaveAs("residualsLowPtY2D-Bpos.eps");
 254
 255   TCanvas* c17Bpos = new TCanvas("c17Bpos","Profiles low Pt-Z, Bpos",800,300);
 256   c17Bpos->Divide(2,1);
 257   c17Bpos->cd(1); proflowptzaBpos->DrawCopy();
 258   c17Bpos->cd(2); proflowptzcBpos->DrawCopy();
 259   c17Bpos->cd(0);
 260   c17Bpos->SaveAs("residualsProfileLowPtZ-Bpos.eps");
 261
 262   TCanvas* c18Bpos = new TCanvas("c18Bpos","Profiles low Pt-Y, Bpos",800,300);
 263   c18Bpos->Divide(2,1);
 264   c18Bpos->cd(1); proflowptyaBpos->DrawCopy();
 265   c18Bpos->cd(2); proflowptycBpos->DrawCopy();
 266   c18Bpos->cd(0);
 267   c18Bpos->SaveAs("residualsProfileLowPtY-Bpos.eps");
 268   }
 269
 270   //------------------------------------------------------------------------------------
 271   pList = dynamic_cast<TList*>(fList->At(1));
 272   TH2F* pZYAResidualsHistBneg = dynamic_cast<TH2F*>(pList->At(0));
 273   TH2F* pZYCResidualsHistBneg = dynamic_cast<TH2F*>(pList->At(1));
 274   TH2F* pLPAResidualsHistBneg = dynamic_cast<TH2F*>(pList->At(2));
 275   TH2F* pLPCResidualsHistBneg = dynamic_cast<TH2F*>(pList->At(3));
 276   TH2F* pPhiYAResidualsHistBneg = dynamic_cast<TH2F*>(pList->At(4));
 277   TH2F* pPhiYCResidualsHistBneg = dynamic_cast<TH2F*>(pList->At(5));
 278   TH2F* pPhiZAResidualsHistBneg = dynamic_cast<TH2F*>(pList->At(6));
 279   TH2F* pPhiZCResidualsHistBneg = dynamic_cast<TH2F*>(pList->At(7));
 280   TH2F* pPtYAResidualsHistBneg = dynamic_cast<TH2F*>(pList->At(8));
 281   TH2F* pPtYCResidualsHistBneg = dynamic_cast<TH2F*>(pList->At(9));
 282   TH2F* pPtZAResidualsHistBneg = dynamic_cast<TH2F*>(pList->At(10));
 283   TH2F* pPtZCResidualsHistBneg = dynamic_cast<TH2F*>(pList->At(11));
 284   TH2F* pLowPtYAResidualsHistBneg = dynamic_cast<TH2F*>(pList->At(12));
 285   TH2F* pLowPtYCResidualsHistBneg = dynamic_cast<TH2F*>(pList->At(13));
 286   TH2F* pLowPtZAResidualsHistBneg = dynamic_cast<TH2F*>(pList->At(14));
 287   TH2F* pLowPtZCResidualsHistBneg = dynamic_cast<TH2F*>(pList->At(15));
 288
 289   TH1D* resycBneg = pZYCResidualsHistBneg->ProjectionY();
 290   resycBneg->SetTitle("r\\phi residual distribution side C (z<0), Bneg");
 291   TH1D* reszcBneg = pZYCResidualsHistBneg->ProjectionX();
 292   reszcBneg->SetTitle("z residual distribution side C (z<0), Bneg");
 293   TH1D* respcBneg = pLPCResidualsHistBneg->ProjectionY();
 294   respcBneg->SetTitle("sin(\\phi) residual distribution side C (z<0), Bneg");
 295   TH1D* reslcBneg = pLPCResidualsHistBneg->ProjectionX();
 296   reslcBneg->SetTitle("tan(\\lambda) residual distribution side C (z<0), Bneg");
 297   TProfile* profphiycBneg = pPhiYCResidualsHistBneg->ProfileX("_pfx",1,-1,"s");
 298   profphiycBneg->SetTitle("\\phi profile of r\\phi residuals side C (z<0), Bneg");
 299   profphiycBneg->SetYTitle("\\deltar\\phi [cm]");
 300   TProfile* profphizcBneg = pPhiZCResidualsHistBneg->ProfileX("_pfx",1,-1,"s");
 301   profphizcBneg->SetTitle("\\phi profile of z residuals side C (z<0), Bneg");
 302   profphizcBneg->SetYTitle("\\deltaz [cm]");
 303   TProfile* profptycBneg = pPtYCResidualsHistBneg->ProfileX("_pfx",1,-1,"s");
 304   profptycBneg->SetTitle("pt profile of r\\phi residuals side C (z<0), Bneg");
 305   profptycBneg->SetYTitle("\\deltar\\phi [cm]");
 306   TProfile* profptzcBneg = pPtZCResidualsHistBneg->ProfileX("_pfx",1,-1,"s");
 307   profptzcBneg->SetTitle("pt profile of z residuals side C (z<0), Bneg");
 308   profptzcBneg->SetYTitle("\\deltaz [cm]");
 309   TProfile* proflowptycBneg = pLowPtYCResidualsHistBneg->ProfileX("_pfx",1,-1,"s");
 310   proflowptycBneg->SetTitle("pt profile of r\\phi residuals side C (z<0), Bneg");
 311   proflowptycBneg->SetYTitle("\\deltar\\phi [cm]");
 312   TProfile* proflowptzcBneg = pLowPtZCResidualsHistBneg->ProfileX("_pfx",1,-1,"s");
 313   proflowptzcBneg->SetTitle("pt profile of z residuals side C (z<0), Bneg");
 314   proflowptzcBneg->SetYTitle("\\deltaz [cm]");
 315
 316   TH1D* resyaBneg = pZYAResidualsHistBneg->ProjectionY();
 317   resyaBneg->SetTitle("r\\phi residual distribution side A (z>0), Bneg");
 318   TH1D* reszaBneg = pZYAResidualsHistBneg->ProjectionX();
 319   reszaBneg->SetTitle("z residual distribution side A (z>0), Bneg");
 320   TH1D* respaBneg = pLPAResidualsHistBneg->ProjectionY();
 321   respaBneg->SetTitle("sin(\\phi) residual distribution side A (z>0), Bneg");
 322   TH1D* reslaBneg = pLPAResidualsHistBneg->ProjectionX();
 323   reslaBneg->SetTitle("tan(\\lambda) residual distribution side A (z>0), Bneg");
 324   TProfile* profphiyaBneg = pPhiYAResidualsHistBneg->ProfileX("_pfx",1,-1,"s");
 325   profphiyaBneg->SetTitle("\\phi profile of r\\phi residuals side A (z>0), Bneg");
 326   profphiyaBneg->SetYTitle("\\deltar\\phi [cm]");
 327   TProfile* profphizaBneg = pPhiZAResidualsHistBneg->ProfileX("_pfx",1,-1,"s");
 328   profphizaBneg->SetTitle("\\phi profile of z residuals side A (z>0), Bneg");
 329   profphizaBneg->SetYTitle("\\deltaz [cm]");
 330   TProfile* profptyaBneg = pPtYAResidualsHistBneg->ProfileX("_pfx",1,-1,"s");
 331   profptyaBneg->SetTitle("pt profile of r\\phi residuals side A (z>0), Bneg");
 332   profptyaBneg->SetYTitle("\\deltar\\phi [cm]");
 333   TProfile* profptzaBneg = pPtZAResidualsHistBneg->ProfileX("_pfx",1,-1,"s");
 334   profptzaBneg->SetTitle("pt profile of z residuals side A (z>0), Bneg");
 335   profptzaBneg->SetYTitle("\\deltaz [cm]");
 336   TProfile* proflowptyaBneg = pLowPtYAResidualsHistBneg->ProfileX("_pfx",1,-1,"s");
 337   proflowptyaBneg->SetTitle("pt profile of r\\phi residuals side A (z>0), Bneg");
 338   proflowptyaBneg->SetYTitle("\\deltar\\phi [cm]");
 339   TProfile* proflowptzaBneg = pLowPtZAResidualsHistBneg->ProfileX("_pfx",1,-1,"s");
 340   proflowptzaBneg->SetTitle("pt profile of z residuals side A (z>0), Bneg");
 341   proflowptzaBneg->SetYTitle("\\deltaz [cm]");
 342
 343   if (pZYAResidualsHistBneg->GetEntries()>0)
 344   {
 345   TCanvas* c1Bneg = new TCanvas("c1Bneg","Residuals ZY 2D, Bneg",800,300);
 346   c1Bneg->Divide(2,1);
 347   c1Bneg->cd(1); pZYAResidualsHistBneg->DrawCopy("col");
 348   c1Bneg->cd(2); pZYCResidualsHistBneg->DrawCopy("col");
 349   c1Bneg->cd(0);
 350   c1Bneg->SaveAs("residualsZY2D-Bneg.eps");
 351
 352   TCanvas* c2Bneg = new TCanvas("c2Bneg","Residuals LP 2D, Bneg",800,300);
 353   c2Bneg->Divide(2,1);
 354   c2Bneg->cd(1); pLPAResidualsHistBneg->DrawCopy("col");
 355   c2Bneg->cd(2); pLPCResidualsHistBneg->DrawCopy("col");
 356   c2Bneg->cd(0);
 357   c2Bneg->SaveAs("residualsLP2D-Bneg.eps");
 358
 359   TCanvas* c3Bneg = new TCanvas("c3Bneg","Residuals z 1D, Bneg",800,300);
 360   c3Bneg->Divide(2,1);
 361   c3Bneg->cd(1); reszaBneg->DrawCopy();
 362   c3Bneg->cd(2); reszcBneg->DrawCopy();
 363   c3Bneg->cd(0);
 364   c3Bneg->SaveAs("residualsZ-Bneg.eps");
 365
 366   TCanvas* c4Bneg = new TCanvas("c4Bneg","Residuals y 1D, Bneg",800,300);
 367   c4Bneg->Divide(2,1);
 368   c4Bneg->cd(1); resyaBneg->DrawCopy();
 369   c4Bneg->cd(2); resycBneg->DrawCopy();
 370   c4Bneg->cd(0);
 371   c4Bneg->SaveAs("residualsY-Bneg.eps");
 372
 373   TCanvas* c5Bneg = new TCanvas("c5Bneg","Residuals phi 1D, Bneg",800,300);
 374   c5Bneg->Divide(2,1);
 375   c5Bneg->cd(1); respaBneg->DrawCopy();
 376   c5Bneg->cd(2); respcBneg->DrawCopy();
 377   c5Bneg->cd(0);
 378   c5Bneg->SaveAs("residualsPhi-Bneg.eps");
 379
 380   TCanvas* c6Bneg = new TCanvas("c6Bneg","Residuals lambda 1D, Bneg",800,300);
 381   c6Bneg->Divide(2,1);
 382   c6Bneg->cd(1); reslaBneg->DrawCopy();
 383   c6Bneg->cd(2); reslcBneg->DrawCopy();
 384   c6Bneg->cd(0);
 385   c6Bneg->SaveAs("residualsLambda-Bneg.eps");
 386
 387   TCanvas* c7Bneg = new TCanvas("c7Bneg","Profiles z in phi, Bneg",800,300);
 388   c7Bneg->Divide(2,1);
 389   c7Bneg->cd(1); profphizaBneg->DrawCopy();
 390   c7Bneg->cd(2); profphizcBneg->DrawCopy();
 391   c7Bneg->cd(0);
 392   c7Bneg->SaveAs("residualsProfilePhiZ-Bneg.eps");
 393
 394   TCanvas* c8Bneg = new TCanvas("c8Bneg","Profiles y in phi, Bneg",800,300);
 395   c8Bneg->Divide(2,1);
 396   c8Bneg->cd(1); profphiyaBneg->DrawCopy();
 397   c8Bneg->cd(2); profphiycBneg->DrawCopy();
 398   c8Bneg->cd(0);
 399   c8Bneg->SaveAs("residualsProfilePhiY-Bneg.eps");
 400
 401   TCanvas* c9Bneg = new TCanvas("c9Bneg", "Residuals Phi-Z 2D, Bneg",800,300);
 402   c9Bneg->Divide(2,1);
 403   c9Bneg->cd(1); pPhiZAResidualsHistBneg->DrawCopy("col");
 404   c9Bneg->cd(2); pPhiZCResidualsHistBneg->DrawCopy("col");
 405   c9Bneg->cd(0);
 406   c9Bneg->SaveAs("residualsPhiZ2D-Bneg.eps");
 407
 408   TCanvas* c1Bneg0 = new TCanvas("c1Bneg0", "Residuals Phi-Y 2D, Bneg",800,300);
 409   c1Bneg0->Divide(2,1);
 410   c1Bneg0->cd(1); pPhiYAResidualsHistBneg->DrawCopy("col");
 411   c1Bneg0->cd(2); pPhiYCResidualsHistBneg->DrawCopy("col");
 412   c1Bneg0->cd(0);
 413   c1Bneg0->SaveAs("residualsPhiY2D-Bneg.eps");
 414
 415   TCanvas* c11Bneg = new TCanvas("c11Bneg", "Residuals Pt-Z 2D, Bneg",800,300);
 416   c11Bneg->Divide(2,1);
 417   c11Bneg->cd(1); pPtZAResidualsHistBneg->DrawCopy("col");
 418   c11Bneg->cd(2); pPtZCResidualsHistBneg->DrawCopy("col");
 419   c11Bneg->cd(0);
 420   c11Bneg->SaveAs("residualsPtZ2D-Bneg.eps");
 421
 422   TCanvas* c12Bneg = new TCanvas("c12Bneg", "Residuals Pt-Y 2D, Bneg",800,300);
 423   c12Bneg->Divide(2,1);
 424   c12Bneg->cd(1); pPtYAResidualsHistBneg->DrawCopy("col");
 425   c12Bneg->cd(2); pPtYCResidualsHistBneg->DrawCopy("col");
 426   c12Bneg->cd(0);
 427   c12Bneg->SaveAs("residualsPtY2D-Bneg.eps");
 428
 429   TCanvas* c13Bneg = new TCanvas("c13Bneg","Profiles Pt-Z, Bneg",800,300);
 430   c13Bneg->Divide(2,1);
 431   c13Bneg->cd(1); profptzaBneg->DrawCopy();
 432   c13Bneg->cd(2); profptzcBneg->DrawCopy();
 433   c13Bneg->cd(0);
 434   c13Bneg->SaveAs("residualsProfilePtZ-Bneg.eps");
 435
 436   TCanvas* c14Bneg = new TCanvas("c14Bneg","Profiles Pt-Y, Bneg",800,300);
 437   c14Bneg->Divide(2,1);
 438   c14Bneg->cd(1); profptyaBneg->DrawCopy();
 439   c14Bneg->cd(2); profptycBneg->DrawCopy();
 440   c14Bneg->cd(0);
 441   c14Bneg->SaveAs("residualsProfilePtY-Bneg.eps");
 442
 443   TCanvas* c15Bneg = new TCanvas("c15Bneg", "Residuals low Pt-Z 2D, Bneg",800,300);
 444   c15Bneg->Divide(2,1);
 445   c15Bneg->cd(1); pLowPtZAResidualsHistBneg->DrawCopy("col");
 446   c15Bneg->cd(2); pLowPtZCResidualsHistBneg->DrawCopy("col");
 447   c15Bneg->cd(0);
 448   c15Bneg->SaveAs("residualsLowPtZ2D-Bneg.eps");
 449
 450   TCanvas* c16Bneg = new TCanvas("c16Bneg", "Residuals in low PtY 2D, Bneg",800,300);
 451   c16Bneg->Divide(2,1);
 452   c16Bneg->cd(1); pLowPtYAResidualsHistBneg->DrawCopy("col");
 453   c16Bneg->cd(2); pLowPtYCResidualsHistBneg->DrawCopy("col");
 454   c16Bneg->cd(0);
 455   c16Bneg->SaveAs("residualsLowPtY2D-Bneg.eps");
 456
 457   TCanvas* c17Bneg = new TCanvas("c17Bneg","Profiles low Pt-Z, Bneg",800,300);
 458   c17Bneg->Divide(2,1);
 459   c17Bneg->cd(1); proflowptzaBneg->DrawCopy();
 460   c17Bneg->cd(2); proflowptzcBneg->DrawCopy();
 461   c17Bneg->cd(0);
 462   c17Bneg->SaveAs("residualsProfileLowPtZ-Bneg.eps");
 463
 464   TCanvas* c18Bneg = new TCanvas("c18Bneg","Profiles low Pt-Y, Bneg",800,300);
 465   c18Bneg->Divide(2,1);
 466   c18Bneg->cd(1); proflowptyaBneg->DrawCopy();
 467   c18Bneg->cd(2); proflowptycBneg->DrawCopy();
 468   c18Bneg->cd(0);
 469   c18Bneg->SaveAs("residualsProfileLowPtY-Bneg.eps");
 470   }
 471
 472   //------------------------------------------------------------------------------------
 473   pList = dynamic_cast<TList*>(fList->At(2));
 474   TH2F* pZYAResidualsHistBnil = dynamic_cast<TH2F*>(pList->At(0));
 475   TH2F* pZYCResidualsHistBnil = dynamic_cast<TH2F*>(pList->At(1));
 476   TH2F* pLPAResidualsHistBnil = dynamic_cast<TH2F*>(pList->At(2));
 477   TH2F* pLPCResidualsHistBnil = dynamic_cast<TH2F*>(pList->At(3));
 478   TH2F* pPhiYAResidualsHistBnil = dynamic_cast<TH2F*>(pList->At(4));
 479   TH2F* pPhiYCResidualsHistBnil = dynamic_cast<TH2F*>(pList->At(5));
 480   TH2F* pPhiZAResidualsHistBnil = dynamic_cast<TH2F*>(pList->At(6));
 481   TH2F* pPhiZCResidualsHistBnil = dynamic_cast<TH2F*>(pList->At(7));
 482   TH2F* pPtYAResidualsHistBnil = dynamic_cast<TH2F*>(pList->At(8));
 483   TH2F* pPtYCResidualsHistBnil = dynamic_cast<TH2F*>(pList->At(9));
 484   TH2F* pPtZAResidualsHistBnil = dynamic_cast<TH2F*>(pList->At(10));
 485   TH2F* pPtZCResidualsHistBnil = dynamic_cast<TH2F*>(pList->At(11));
 486   TH2F* pLowPtYAResidualsHistBnil = dynamic_cast<TH2F*>(pList->At(12));
 487   TH2F* pLowPtYCResidualsHistBnil = dynamic_cast<TH2F*>(pList->At(13));
 488   TH2F* pLowPtZAResidualsHistBnil = dynamic_cast<TH2F*>(pList->At(14));
 489   TH2F* pLowPtZCResidualsHistBnil = dynamic_cast<TH2F*>(pList->At(15));
 490
 491   TH1D* resycBnil = pZYCResidualsHistBnil->ProjectionY();
 492   resycBnil->SetTitle("r\\phi residual distribution side C (z<0), Bnil");
 493   TH1D* reszcBnil = pZYCResidualsHistBnil->ProjectionX();
 494   reszcBnil->SetTitle("z residual distribution side C (z<0), Bnil");
 495   TH1D* respcBnil = pLPCResidualsHistBnil->ProjectionY();
 496   respcBnil->SetTitle("sin(\\phi) residual distribution side C (z<0), Bnil");
 497   TH1D* reslcBnil = pLPCResidualsHistBnil->ProjectionX();
 498   reslcBnil->SetTitle("tan(\\lambda) residual distribution side C (z<0), Bnil");
 499   TProfile* profphiycBnil = pPhiYCResidualsHistBnil->ProfileX("_pfx",1,-1,"s");
 500   profphiycBnil->SetTitle("\\phi profile of r\\phi residuals side C (z<0), Bnil");
 501   profphiycBnil->SetYTitle("\\deltar\\phi [cm]");
 502   TProfile* profphizcBnil = pPhiZCResidualsHistBnil->ProfileX("_pfx",1,-1,"s");
 503   profphizcBnil->SetTitle("\\phi profile of z residuals side C (z<0), Bnil");
 504   profphizcBnil->SetYTitle("\\deltaz [cm]");
 505   TProfile* profptycBnil = pPtYCResidualsHistBnil->ProfileX("_pfx",1,-1,"s");
 506   profptycBnil->SetTitle("pt profile of r\\phi residuals side C (z<0), Bnil");
 507   profptycBnil->SetYTitle("\\deltar\\phi [cm]");
 508   TProfile* profptzcBnil = pPtZCResidualsHistBnil->ProfileX("_pfx",1,-1,"s");
 509   profptzcBnil->SetTitle("pt profile of z residuals side C (z<0), Bnil");
 510   profptzcBnil->SetYTitle("\\deltaz [cm]");
 511   TProfile* proflowptycBnil = pLowPtYCResidualsHistBnil->ProfileX("_pfx",1,-1,"s");
 512   proflowptycBnil->SetTitle("pt profile of r\\phi residuals side C (z<0), Bnil");
 513   proflowptycBnil->SetYTitle("\\deltar\\phi [cm]");
 514   TProfile* proflowptzcBnil = pLowPtZCResidualsHistBnil->ProfileX("_pfx",1,-1,"s");
 515   proflowptzcBnil->SetTitle("pt profile of z residuals side C (z<0), Bnil");
 516   proflowptzcBnil->SetYTitle("\\deltaz [cm]");
 517
 518   TH1D* resyaBnil = pZYAResidualsHistBnil->ProjectionY();
 519   resyaBnil->SetTitle("r\\phi residual distribution side A (z>0), Bnil");
 520   TH1D* reszaBnil = pZYAResidualsHistBnil->ProjectionX();
 521   reszaBnil->SetTitle("z residual distribution side A (z>0), Bnil");
 522   TH1D* respaBnil = pLPAResidualsHistBnil->ProjectionY();
 523   respaBnil->SetTitle("sin(\\phi) residual distribution side A (z>0), Bnil");
 524   TH1D* reslaBnil = pLPAResidualsHistBnil->ProjectionX();
 525   reslaBnil->SetTitle("tan(\\lambda) residual distribution side A (z>0), Bnil");
 526   TProfile* profphiyaBnil = pPhiYAResidualsHistBnil->ProfileX("_pfx",1,-1,"s");
 527   profphiyaBnil->SetTitle("\\phi profile of r\\phi residuals side A (z>0), Bnil");
 528   profphiyaBnil->SetYTitle("\\deltar\\phi [cm]");
 529   TProfile* profphizaBnil = pPhiZAResidualsHistBnil->ProfileX("_pfx",1,-1,"s");
 530   profphizaBnil->SetTitle("\\phi profile of z residuals side A (z>0), Bnil");
 531   profphizaBnil->SetYTitle("\\deltaz [cm]");
 532   TProfile* profptyaBnil = pPtYAResidualsHistBnil->ProfileX("_pfx",1,-1,"s");
 533   profptyaBnil->SetTitle("pt profile of r\\phi residuals side A (z>0), Bnil");
 534   profptyaBnil->SetYTitle("\\deltar\\phi [cm]");
 535   TProfile* profptzaBnil = pPtZAResidualsHistBnil->ProfileX("_pfx",1,-1,"s");
 536   profptzaBnil->SetTitle("pt profile of z residuals side A (z>0), Bnil");
 537   profptzaBnil->SetYTitle("\\deltaz [cm]");
 538   TProfile* proflowptyaBnil = pLowPtYAResidualsHistBnil->ProfileX("_pfx",1,-1,"s");
 539   proflowptyaBnil->SetTitle("pt profile of r\\phi residuals side A (z>0), Bnil");
 540   proflowptyaBnil->SetYTitle("\\deltar\\phi [cm]");
 541   TProfile* proflowptzaBnil = pLowPtZAResidualsHistBnil->ProfileX("_pfx",1,-1,"s");
 542   proflowptzaBnil->SetTitle("pt profile of z residuals side A (z>0), Bnil");
 543   proflowptzaBnil->SetYTitle("\\deltaz [cm]");
 544
 545   if (pZYAResidualsHistBnil->GetEntries()>0)
 546   {
 547   TCanvas* c1Bnil = new TCanvas("c1Bnil","Residuals ZY 2D, Bnil",800,300);
 548   c1Bnil->Divide(2,1);
 549   c1Bnil->cd(1); pZYAResidualsHistBnil->DrawCopy("col");
 550   c1Bnil->cd(2); pZYCResidualsHistBnil->DrawCopy("col");
 551   c1Bnil->cd(0);
 552   c1Bnil->SaveAs("residualsZY2D-Bnil.eps");
 553
 554   TCanvas* c2Bnil = new TCanvas("c2Bnil","Residuals LP 2D, Bnil",800,300);
 555   c2Bnil->Divide(2,1);
 556   c2Bnil->cd(1); pLPAResidualsHistBnil->DrawCopy("col");
 557   c2Bnil->cd(2); pLPCResidualsHistBnil->DrawCopy("col");
 558   c2Bnil->cd(0);
 559   c2Bnil->SaveAs("residualsLP2D-Bnil.eps");
 560
 561   TCanvas* c3Bnil = new TCanvas("c3Bnil","Residuals z 1D, Bnil",800,300);
 562   c3Bnil->Divide(2,1);
 563   c3Bnil->cd(1); reszaBnil->DrawCopy();
 564   c3Bnil->cd(2); reszcBnil->DrawCopy();
 565   c3Bnil->cd(0);
 566   c3Bnil->SaveAs("residualsZ-Bnil.eps");
 567
 568   TCanvas* c4Bnil = new TCanvas("c4Bnil","Residuals y 1D, Bnil",800,300);
 569   c4Bnil->Divide(2,1);
 570   c4Bnil->cd(1); resyaBnil->DrawCopy();
 571   c4Bnil->cd(2); resycBnil->DrawCopy();
 572   c4Bnil->cd(0);
 573   c4Bnil->SaveAs("residualsY-Bnil.eps");
 574
 575   TCanvas* c5Bnil = new TCanvas("c5Bnil","Residuals phi 1D, Bnil",800,300);
 576   c5Bnil->Divide(2,1);
 577   c5Bnil->cd(1); respaBnil->DrawCopy();
 578   c5Bnil->cd(2); respcBnil->DrawCopy();
 579   c5Bnil->cd(0);
 580   c5Bnil->SaveAs("residualsPhi-Bnil.eps");
 581
 582   TCanvas* c6Bnil = new TCanvas("c6Bnil","Residuals lambda 1D, Bnil",800,300);
 583   c6Bnil->Divide(2,1);
 584   c6Bnil->cd(1); reslaBnil->DrawCopy();
 585   c6Bnil->cd(2); reslcBnil->DrawCopy();
 586   c6Bnil->cd(0);
 587   c6Bnil->SaveAs("residualsLambda-Bnil.eps");
 588
 589   TCanvas* c7Bnil = new TCanvas("c7Bnil","Profiles z in phi, Bnil",800,300);
 590   c7Bnil->Divide(2,1);
 591   c7Bnil->cd(1); profphizaBnil->DrawCopy();
 592   c7Bnil->cd(2); profphizcBnil->DrawCopy();
 593   c7Bnil->cd(0);
 594   c7Bnil->SaveAs("residualsProfilePhiZ-Bnil.eps");
 595
 596   TCanvas* c8Bnil = new TCanvas("c8Bnil","Profiles y in phi, Bnil",800,300);
 597   c8Bnil->Divide(2,1);
 598   c8Bnil->cd(1); profphiyaBnil->DrawCopy();
 599   c8Bnil->cd(2); profphiycBnil->DrawCopy();
 600   c8Bnil->cd(0);
 601   c8Bnil->SaveAs("residualsProfilePhiY-Bnil.eps");
 602
 603   TCanvas* c9Bnil = new TCanvas("c9Bnil", "Residuals Phi-Z 2D, Bnil",800,300);
 604   c9Bnil->Divide(2,1);
 605   c9Bnil->cd(1); pPhiZAResidualsHistBnil->DrawCopy("col");
 606   c9Bnil->cd(2); pPhiZCResidualsHistBnil->DrawCopy("col");
 607   c9Bnil->cd(0);
 608   c9Bnil->SaveAs("residualsPhiZ2D-Bnil.eps");
 609
 610   TCanvas* c1Bnil0 = new TCanvas("c1Bnil0", "Residuals Phi-Y 2D, Bnil",800,300);
 611   c1Bnil0->Divide(2,1);
 612   c1Bnil0->cd(1); pPhiYAResidualsHistBnil->DrawCopy("col");
 613   c1Bnil0->cd(2); pPhiYCResidualsHistBnil->DrawCopy("col");
 614   c1Bnil0->cd(0);
 615   c1Bnil0->SaveAs("residualsPhiY2D-Bnil.eps");
 616
 617   TCanvas* c11Bnil = new TCanvas("c11Bnil", "Residuals Pt-Z 2D, Bnil",800,300);
 618   c11Bnil->Divide(2,1);
 619   c11Bnil->cd(1); pPtZAResidualsHistBnil->DrawCopy("col");
 620   c11Bnil->cd(2); pPtZCResidualsHistBnil->DrawCopy("col");
 621   c11Bnil->cd(0);
 622   c11Bnil->SaveAs("residualsPtZ2D-Bnil.eps");
 623
 624   TCanvas* c12Bnil = new TCanvas("c12Bnil", "Residuals Pt-Y 2D, Bnil",800,300);
 625   c12Bnil->Divide(2,1);
 626   c12Bnil->cd(1); pPtYAResidualsHistBnil->DrawCopy("col");
 627   c12Bnil->cd(2); pPtYCResidualsHistBnil->DrawCopy("col");
 628   c12Bnil->cd(0);
 629   c12Bnil->SaveAs("residualsPtY2D-Bnil.eps");
 630
 631   TCanvas* c13Bnil = new TCanvas("c13Bnil","Profiles Pt-Z, Bnil",800,300);
 632   c13Bnil->Divide(2,1);
 633   c13Bnil->cd(1); profptzaBnil->DrawCopy();
 634   c13Bnil->cd(2); profptzcBnil->DrawCopy();
 635   c13Bnil->cd(0);
 636   c13Bnil->SaveAs("residualsProfilePtZ-Bnil.eps");
 637
 638   TCanvas* c14Bnil = new TCanvas("c14Bnil","Profiles Pt-Y, Bnil",800,300);
 639   c14Bnil->Divide(2,1);
 640   c14Bnil->cd(1); profptyaBnil->DrawCopy();
 641   c14Bnil->cd(2); profptycBnil->DrawCopy();
 642   c14Bnil->cd(0);
 643   c14Bnil->SaveAs("residualsProfilePtY-Bnil.eps");
 644
 645   TCanvas* c15Bnil = new TCanvas("c15Bnil", "Residuals low Pt-Z 2D, Bnil",800,300);
 646   c15Bnil->Divide(2,1);
 647   c15Bnil->cd(1); pLowPtZAResidualsHistBnil->DrawCopy("col");
 648   c15Bnil->cd(2); pLowPtZCResidualsHistBnil->DrawCopy("col");
 649   c15Bnil->cd(0);
 650   c15Bnil->SaveAs("residualsLowPtZ2D-Bnil.eps");
 651
 652   TCanvas* c16Bnil = new TCanvas("c16Bnil", "Residuals in low PtY 2D, Bnil",800,300);
 653   c16Bnil->Divide(2,1);
 654   c16Bnil->cd(1); pLowPtYAResidualsHistBnil->DrawCopy("col");
 655   c16Bnil->cd(2); pLowPtYCResidualsHistBnil->DrawCopy("col");
 656   c16Bnil->cd(0);
 657   c16Bnil->SaveAs("residualsLowPtY2D-Bnil.eps");
 658
 659   TCanvas* c17Bnil = new TCanvas("c17Bnil","Profiles low Pt-Z, Bnil",800,300);
 660   c17Bnil->Divide(2,1);
 661   c17Bnil->cd(1); proflowptzaBnil->DrawCopy();
 662   c17Bnil->cd(2); proflowptzcBnil->DrawCopy();
 663   c17Bnil->cd(0);
 664   c17Bnil->SaveAs("residualsProfileLowPtZ-Bnil.eps");
 665
 666   TCanvas* c18Bnil = new TCanvas("c18Bnil","Profiles low Pt-Y, Bnil",800,300);
 667   c18Bnil->Divide(2,1);
 668   c18Bnil->cd(1); proflowptyaBnil->DrawCopy();
 669   c18Bnil->cd(2); proflowptycBnil->DrawCopy();
 670   c18Bnil->cd(0);
 671   c18Bnil->SaveAs("residualsProfileLowPtY-Bnil.eps");
 672   }
 673 }
 674