[u/mrichter/AliRoot.git] / prod / acrcaf / qa / qa.C

void qa(Int_t runNumber) {
  TStopwatch timer;
  timer.Start();

  gEnv->SetValue("XSec.GSI.DelegProxy","2");
  // Select ROOT version
  TProof::Mgr("aliprod@alicecaf")->SetROOTVersion("v5-24-00b-caf");
  // Login to CAF
  TProof::Open("aliprod@alicecaf");

  // Enable AliRoot
  gProof->UploadPackage("/afs/cern.ch/alice/caf/sw/ALICE/PARs/v4-17-Release.rec/AF-v4-17-rec.par");
  gProof->EnablePackage("AF-v4-17-rec.par");

  // Enable analysis libs
  gSystem->Load("libANALYSIS.so");
  gSystem->Load("libANALYSISalice.so");
  gSystem->Load("libPWG1.so");
  gProof->Exec("gSystem->Load(\"libANALYSIS.so\");",kTRUE);
  gProof->Exec("gSystem->Load(\"libANALYSISalice.so\");",kTRUE);
  gProof->Exec("gSystem->Load(\"libPWG1.so\");",kTRUE);

  //____________________________________________//
  // Make the analysis manager
  AliAnalysisManager *mgr = new AliAnalysisManager("AnalysisQAManager");
  AliESDInputHandler* esdH = new AliESDInputHandler();
  mgr->SetInputEventHandler(esdH);  
  mgr->SetDebugLevel(10);

  //____________________________________________//
  // 1st Cosmic task
  AliAnalysisTaskCosmic *task1 = new AliAnalysisTaskCosmic("TaskCosmic");
  mgr->AddTask(task1);

  // Create containers for input/output
  AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("chist1",TList::Class(),AliAnalysisManager::kOutputContainer,
							    Form("run%d.root",runNumber));

  //____________________________________________//
  mgr->ConnectInput  (task1,  0, mgr->GetCommonInputContainer());
  mgr->ConnectOutput (task1,  1, coutput1);
  
  if (!mgr->InitAnalysis()) return;
  mgr->PrintStatus();
  mgr->StartAnalysis("proof",
		     Form("/ALIREC/aliprod/run%d",runNumber));

  timer.Stop();
  timer.Print();

  plot(runNumber);
}

void plot(Int_t runNumber)
{
  TFile* f = new TFile(Form("run%d.root",runNumber), "read");

  // pt, phi ....
  TCanvas* c1 = new TCanvas("c1", "pt,eta,phi", 10, 10, 1100, 800);
  c1->Divide(2,2);
  c1->cd(1);
  c1->GetPad(1)->SetLogy();
    
  fhPtP = (TH1F*) chist1->FindObject("fhPtPC");
  fhPtP->SetLineColor(2);
  fhPtP->Draw();
  fhPtN = (TH1F*) chist1->FindObject("fhPtNC");
  fhPtN->SetLineColor(4);
  fhPtN->Draw("same");

  c1->cd(2);
  fhPhiP = (TH1F*) chist1->FindObject("fhPhiPC");
  fhPhiP->SetLineColor(2);
  fhPhiP->Draw();
  fhPhiN = (TH1F*) chist1->FindObject("fhPhiNC");
  fhPhiN->SetLineColor(4);
  fhPhiN->Draw("same");

  c1->cd(3);
  fhThetaP = (TH1F*) chist1->FindObject("fhThetaPC");
  fhThetaP->SetLineColor(2);
  fhThetaP->Draw();
  fhThetaN = (TH1F*) chist1->FindObject("fhThetaNC");
  fhThetaN->SetLineColor(4);
  fhThetaN->Draw("same");

    
  // delta phi
  TCanvas* c2 = new TCanvas("c2", "#Delta#phi", 10, 10, 1100, 800);
  c2->Divide(2,2);
  // pos/neg charge
  c2->cd(1);
  c2->GetPad(1)->SetLogy();
  fhDPhiP = (TH1F*) chist1->FindObject("fhDPhiPC");
  fhDPhiN = (TH1F*) chist1->FindObject("fhDPhiNC");
  fhDPhiP->SetLineColor(2);
  fhDPhiN->SetLineColor(4);
  fhDPhiP->SetXTitle("#Delta#phi [rad]");
  fhDPhiP->Draw();
  fhDPhiN->Draw("same");

  // pos/neg z
  c2->cd(2);
  c2->GetPad(2)->SetLogy();

  fhDPhiZP = (TH1F*) chist1->FindObject("fhDPhiPZ");
  fhDPhiZN = (TH1F*) chist1->FindObject("fhDPhiNZ");
  fhDPhiBA = (TH1F*) chist1->FindObject("fhDPhi_Bad");
  fhDPhiZP->SetLineColor(2);
  fhDPhiZN->SetLineColor(4);
  fhDPhiBA->SetLineColor(6);

  fhDPhiZP->SetXTitle("#Delta#phi [rad]");
  fhDPhiZP->Draw();
  fhDPhiZN->Draw("same");
  //  fhDPhiBA->Draw("same");

    
  // delta theta
  // pos/neg charge
  c2->cd(3);
  c2->GetPad(3)->SetLogy();
  fhDThetaP = (TH1F*) chist1->FindObject("fhDThetaPC");
  fhDThetaN = (TH1F*) chist1->FindObject("fhDThetaNC");
  fhDThetaP->SetLineColor(2);
  fhDThetaN->SetLineColor(4);
 

  fhDThetaP->SetXTitle("#Delta#theta [rad]");
  fhDThetaP->Draw();
  fhDThetaN->Draw("same");

  // pos/neg z
  c2->cd(4);
  c2->GetPad(4)->SetLogy();

  fhDThetaZP = (TH1F*) chist1->FindObject("fhDThetaPZ");
  fhDThetaZN = (TH1F*) chist1->FindObject("fhDThetaNZ");
  fhDThetaBA = (TH1F*) chist1->FindObject("fhDTheta_Bad");
  fhDThetaZP->SetLineColor(2);
  fhDThetaZN->SetLineColor(4);
  fhDThetaBA->SetLineColor(6);
  fhDThetaZP->SetXTitle("#Delta#theta [rad]");
  fhDThetaZP->Draw();

  fhDThetaZN->Draw("same");
  //    fhDThetaBA->Draw("same");

    
  // delta Pt
  TCanvas* c3 = new TCanvas("c3", "#Delta p_{T}", 10, 10, 1100, 800);
  c3->Divide(2,2);
  // pos/neg charge
  c3->cd(1);
  c3->GetPad(1)->SetLogy();
  fhDPtP = (TH1F*) chist1->FindObject("fhDPtPC");
  fhDPtN = (TH1F*) chist1->FindObject("fhDPtNC");
  fhDPtP->SetLineColor(2);
  fhDPtN->SetLineColor(4);
    
  fhDPtP->Draw();
  fhDPtN->Draw("same");

  // pos/neg z
  c3->cd(2);
  c3->GetPad(2)->SetLogy();

  fhDPtZP = (TH1F*) chist1->FindObject("fhDPtPZ");
  fhDPtZN = (TH1F*) chist1->FindObject("fhDPtNZ");
  fhDPtZP->SetLineColor(2);
  fhDPtZN->SetLineColor(4);
  fhDPtZP->Draw();
  fhDPtZN->Draw("same");

  c3->cd(3);

  fpDPtP = (TH1F*) chist1->FindObject("fpDPtPC");
  fpDPtN = (TH1F*) chist1->FindObject("fpDPtNC");
  fpDPtP->SetLineColor(2);
  fpDPtN->SetLineColor(4);
  fpDPtP->Draw();
  fpDPtN->Draw("same");

  c3->cd(4);

  fpDPtZP = (TH1F*) chist1->FindObject("fpDPtPZ");
  fpDPtZN = (TH1F*) chist1->FindObject("fpDPtNZ");
  fpDPtZP->SetLineColor(2);
  fpDPtZN->SetLineColor(4);
  fpDPtZP->Draw();
  fpDPtZN->Draw("same");
    

  // dZ
  TCanvas* c4 = new TCanvas("c4", "#Delta Z", 10, 10, 1100, 800);
  c4->Divide(2,2);
  // pos/neg charge
  c4->cd(1);
  c4->GetPad(1)->SetLogy();
  fhDZP = (TH1F*) chist1->FindObject("fhDZPC");
  fhDZN = (TH1F*) chist1->FindObject("fhDZNC");
  fhDZP->SetLineColor(2);
  fhDZN->SetLineColor(4);
    
  fhDZP->Draw();
  fhDZN->Draw("same");
  
  // pos/neg z
  c4->cd(2);
  c4->GetPad(2)->SetLogy();

  fhDZZP = (TH1F*) chist1->FindObject("fhDZPZ");
  fhDZZN = (TH1F*) chist1->FindObject("fhDZNZ");
  fhDZBA = (TH1F*) chist1->FindObject("fhDZ_Bad");

  fhDZZP->SetLineColor(2);
  fhDZZN->SetLineColor(4);
  fhDZBA->SetLineColor(6);
  
  fhDZZP->Draw();
  fhDZZN->Draw("same");
  fhDZBA->Draw("same");

  // dX
  // pos/neg charge
  c4->cd(3);
  c4->GetPad(3)->SetLogy();
  fhDXP = (TH1F*) chist1->FindObject("fhDXPC");
  fhDXN = (TH1F*) chist1->FindObject("fhDXNC");
  fhDXP->SetLineColor(2);
  fhDXN->SetLineColor(4);
    
  fhDXP->Draw();
  fhDXN->Draw("same");

  // pos/neg z
  c4->cd(4);
  c4->GetPad(4)->SetLogy();

  fhDXZP = (TH1F*) chist1->FindObject("fhDXPZ");
  fhDXZN = (TH1F*) chist1->FindObject("fhDXNZ");
  fhDXBA = (TH1F*) chist1->FindObject("fhDX_Bad");

  fhDXZP->SetLineColor(2);
  fhDXZN->SetLineColor(4);
  fhDXBA->SetLineColor(6);

  fhDXZP->Draw();
  fhDXZN->Draw("same");
  fhDXBA->Draw("same");

  // dY
  TCanvas* c4a = new TCanvas("c4a", "#Delta Y", 10, 10, 1100, 800);
  c4a->Divide(2,2);
  // pos/neg charge
  c4a->cd(1);
  c4a->GetPad(1)->SetLogy();
  fhDYP = (TH1F*) chist1->FindObject("fhDYPC");
  fhDYN = (TH1F*) chist1->FindObject("fhDYNC");
  fhDYP->SetLineColor(2);
  fhDYN->SetLineColor(4);
  
  fhDYP->Draw();
  fhDYN->Draw("same");

  // pos/neg z
  c4a->cd(2);
  c4a->GetPad(2)->SetLogy();

  fhDYZP = (TH1F*) chist1->FindObject("fhDYPZ");
  fhDYZN = (TH1F*) chist1->FindObject("fhDYNZ");
  fhDYBA = (TH1F*) chist1->FindObject("fhDY_Bad");
  
  fhDYZP->SetLineColor(2);
  fhDYZN->SetLineColor(4);
  fhDYBA->SetLineColor(6);

  fhDYZP->Draw();
  fhDYZN->Draw("same");
  fhDYBA->Draw("same");

  // delta Pt
  TCanvas* c5 = new TCanvas("c5", "#Delta p_{T} (n-sigma)", 10, 10, 1100, 800);
  c5->Divide(2,2);
  c5->cd(1);

  fpDPtSP = (TH1F*) chist1->FindObject("fpDPtSPC");
  fpDPtSN = (TH1F*) chist1->FindObject("fpDPtSNC");
  fpDPtSP->SetLineColor(2);
  fpDPtSN->SetLineColor(4);
  fpDPtSP->Draw();
  fpDPtSN->Draw("same");


  c5->cd(2);

  fpDPtSPZ = (TH1F*) chist1->FindObject("fpDPtSPZ");
  fpDPtSNZ = (TH1F*) chist1->FindObject("fpDPtSNZ");
  fpDPtSPZ->SetLineColor(2);
  fpDPtSNZ->SetLineColor(4);
  fpDPtSPZ->Draw();
  fpDPtSNZ->Draw("same");

  TCanvas* c6 = new TCanvas("c6", "Dz vs z", 10, 10, 1100, 800);
  fhDZvsZ = (TH2F*) chist1->FindObject("fhDZvsZ");
  fhDZvsZ->SetXTitle("z_{in} * sign(z_{in}) * sign(z_{out}) [cm]");
  fhDZvsZ->SetYTitle("#DeltaZ [cm]");

  gStyle->SetPalette(1);
    
  fhDZvsZ->Draw("colz");
}
Commit	Line	Data
c2174181	1	void qa(Int_t runNumber) {
	2	TStopwatch timer;
	3	timer.Start();
	4
	5	gEnv->SetValue("XSec.GSI.DelegProxy","2");
	6	// Select ROOT version
	7	TProof::Mgr("aliprod@alicecaf")->SetROOTVersion("v5-24-00b-caf");
	8	// Login to CAF
	9	TProof::Open("aliprod@alicecaf");
	10
	11	// Enable AliRoot
	12	gProof->UploadPackage("/afs/cern.ch/alice/caf/sw/ALICE/PARs/v4-17-Release.rec/AF-v4-17-rec.par");
	13	gProof->EnablePackage("AF-v4-17-rec.par");
	14
a7ec3dd4	15	// Enable analysis libs
	16	gSystem->Load("libANALYSIS.so");
	17	gSystem->Load("libANALYSISalice.so");
13424c7e	18	gSystem->Load("libPWG1.so");
a7ec3dd4	19	gProof->Exec("gSystem->Load(\"libANALYSIS.so\");",kTRUE);
a7ec3dd4	20	gProof->Exec("gSystem->Load(\"libANALYSISalice.so\");",kTRUE);
13424c7e	21	gProof->Exec("gSystem->Load(\"libPWG1.so\");",kTRUE);
c2174181	22
	23	//____________________________________________//
	24	// Make the analysis manager
	25	AliAnalysisManager *mgr = new AliAnalysisManager("AnalysisQAManager");
	26	AliESDInputHandler* esdH = new AliESDInputHandler();
	27	mgr->SetInputEventHandler(esdH);
	28	mgr->SetDebugLevel(10);
	29
	30	//____________________________________________//
	31	// 1st Cosmic task
	32	AliAnalysisTaskCosmic *task1 = new AliAnalysisTaskCosmic("TaskCosmic");
	33	mgr->AddTask(task1);
	34
	35	// Create containers for input/output
	36	AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("chist1",TList::Class(),AliAnalysisManager::kOutputContainer,
	37	Form("run%d.root",runNumber));
	38
	39	//____________________________________________//
	40	mgr->ConnectInput (task1, 0, mgr->GetCommonInputContainer());
	41	mgr->ConnectOutput (task1, 1, coutput1);
	42
	43	if (!mgr->InitAnalysis()) return;
	44	mgr->PrintStatus();
	45	mgr->StartAnalysis("proof",
a7ec3dd4	46	Form("/ALIREC/aliprod/run%d",runNumber));
c2174181	47
	48	timer.Stop();
	49	timer.Print();
	50
	51	plot(runNumber);
	52	}
	53
	54	void plot(Int_t runNumber)
	55	{
	56	TFile* f = new TFile(Form("run%d.root",runNumber), "read");
	57
	58	// pt, phi ....
	59	TCanvas* c1 = new TCanvas("c1", "pt,eta,phi", 10, 10, 1100, 800);
	60	c1->Divide(2,2);
	61	c1->cd(1);
	62	c1->GetPad(1)->SetLogy();
	63
	64	fhPtP = (TH1F*) chist1->FindObject("fhPtPC");
	65	fhPtP->SetLineColor(2);
	66	fhPtP->Draw();
	67	fhPtN = (TH1F*) chist1->FindObject("fhPtNC");
	68	fhPtN->SetLineColor(4);
	69	fhPtN->Draw("same");
	70
	71	c1->cd(2);
	72	fhPhiP = (TH1F*) chist1->FindObject("fhPhiPC");
	73	fhPhiP->SetLineColor(2);
	74	fhPhiP->Draw();
	75	fhPhiN = (TH1F*) chist1->FindObject("fhPhiNC");
	76	fhPhiN->SetLineColor(4);
	77	fhPhiN->Draw("same");
	78
	79	c1->cd(3);
	80	fhThetaP = (TH1F*) chist1->FindObject("fhThetaPC");
	81	fhThetaP->SetLineColor(2);
	82	fhThetaP->Draw();
	83	fhThetaN = (TH1F*) chist1->FindObject("fhThetaNC");
	84	fhThetaN->SetLineColor(4);
	85	fhThetaN->Draw("same");
	86
	87
	88	// delta phi
	89	TCanvas* c2 = new TCanvas("c2", "#Delta#phi", 10, 10, 1100, 800);
	90	c2->Divide(2,2);
	91	// pos/neg charge
	92	c2->cd(1);
	93	c2->GetPad(1)->SetLogy();
	94	fhDPhiP = (TH1F*) chist1->FindObject("fhDPhiPC");
	95	fhDPhiN = (TH1F*) chist1->FindObject("fhDPhiNC");
	96	fhDPhiP->SetLineColor(2);
	97	fhDPhiN->SetLineColor(4);
	98	fhDPhiP->SetXTitle("#Delta#phi [rad]");
	99	fhDPhiP->Draw();
	100	fhDPhiN->Draw("same");
	101
	102	// pos/neg z
	103	c2->cd(2);
	104	c2->GetPad(2)->SetLogy();
	105
	106	fhDPhiZP = (TH1F*) chist1->FindObject("fhDPhiPZ");
	107	fhDPhiZN = (TH1F*) chist1->FindObject("fhDPhiNZ");
	108	fhDPhiBA = (TH1F*) chist1->FindObject("fhDPhi_Bad");
	109	fhDPhiZP->SetLineColor(2);
	110	fhDPhiZN->SetLineColor(4);
111	fhDPhiBA->SetLineColor(6);
112
113	fhDPhiZP->SetXTitle("#Delta#phi [rad]");
114	fhDPhiZP->Draw();
115	fhDPhiZN->Draw("same");
116	// fhDPhiBA->Draw("same");
117
118
119	// delta theta
120	// pos/neg charge
121	c2->cd(3);
122	c2->GetPad(3)->SetLogy();
123	fhDThetaP = (TH1F*) chist1->FindObject("fhDThetaPC");
124	fhDThetaN = (TH1F*) chist1->FindObject("fhDThetaNC");
125	fhDThetaP->SetLineColor(2);
126	fhDThetaN->SetLineColor(4);
127
128
129	fhDThetaP->SetXTitle("#Delta#theta [rad]");
130	fhDThetaP->Draw();
131	fhDThetaN->Draw("same");
132
133	// pos/neg z
134	c2->cd(4);
135	c2->GetPad(4)->SetLogy();
136
137	fhDThetaZP = (TH1F*) chist1->FindObject("fhDThetaPZ");
138	fhDThetaZN = (TH1F*) chist1->FindObject("fhDThetaNZ");
139	fhDThetaBA = (TH1F*) chist1->FindObject("fhDTheta_Bad");
140	fhDThetaZP->SetLineColor(2);
141	fhDThetaZN->SetLineColor(4);
142	fhDThetaBA->SetLineColor(6);
143	fhDThetaZP->SetXTitle("#Delta#theta [rad]");
144	fhDThetaZP->Draw();
145
146	fhDThetaZN->Draw("same");
147	// fhDThetaBA->Draw("same");
148
149
150	// delta Pt
151	TCanvas* c3 = new TCanvas("c3", "#Delta p_{T}", 10, 10, 1100, 800);
152	c3->Divide(2,2);
153	// pos/neg charge
154	c3->cd(1);
155	c3->GetPad(1)->SetLogy();
156	fhDPtP = (TH1F*) chist1->FindObject("fhDPtPC");
157	fhDPtN = (TH1F*) chist1->FindObject("fhDPtNC");
158	fhDPtP->SetLineColor(2);
159	fhDPtN->SetLineColor(4);
160
161	fhDPtP->Draw();
162	fhDPtN->Draw("same");
163
164	// pos/neg z
165	c3->cd(2);
166	c3->GetPad(2)->SetLogy();
167
168	fhDPtZP = (TH1F*) chist1->FindObject("fhDPtPZ");
169	fhDPtZN = (TH1F*) chist1->FindObject("fhDPtNZ");
170	fhDPtZP->SetLineColor(2);
171	fhDPtZN->SetLineColor(4);
172	fhDPtZP->Draw();
173	fhDPtZN->Draw("same");
174
175	c3->cd(3);
176
177	fpDPtP = (TH1F*) chist1->FindObject("fpDPtPC");
178	fpDPtN = (TH1F*) chist1->FindObject("fpDPtNC");
179	fpDPtP->SetLineColor(2);
180	fpDPtN->SetLineColor(4);
181	fpDPtP->Draw();
182	fpDPtN->Draw("same");
183
184	c3->cd(4);
185
186	fpDPtZP = (TH1F*) chist1->FindObject("fpDPtPZ");
187	fpDPtZN = (TH1F*) chist1->FindObject("fpDPtNZ");
188	fpDPtZP->SetLineColor(2);
189	fpDPtZN->SetLineColor(4);
190	fpDPtZP->Draw();
191	fpDPtZN->Draw("same");
192
193
194	// dZ
195	TCanvas* c4 = new TCanvas("c4", "#Delta Z", 10, 10, 1100, 800);
196	c4->Divide(2,2);
197	// pos/neg charge
198	c4->cd(1);
199	c4->GetPad(1)->SetLogy();
200	fhDZP = (TH1F*) chist1->FindObject("fhDZPC");
201	fhDZN = (TH1F*) chist1->FindObject("fhDZNC");
202	fhDZP->SetLineColor(2);
203	fhDZN->SetLineColor(4);
204
205	fhDZP->Draw();
206	fhDZN->Draw("same");
207
208	// pos/neg z
209	c4->cd(2);
210	c4->GetPad(2)->SetLogy();
211
212	fhDZZP = (TH1F*) chist1->FindObject("fhDZPZ");
213	fhDZZN = (TH1F*) chist1->FindObject("fhDZNZ");
214	fhDZBA = (TH1F*) chist1->FindObject("fhDZ_Bad");
215
216	fhDZZP->SetLineColor(2);
217	fhDZZN->SetLineColor(4);
218	fhDZBA->SetLineColor(6);
219
220	fhDZZP->Draw();
221	fhDZZN->Draw("same");
222	fhDZBA->Draw("same");
223
224	// dX
225	// pos/neg charge
226	c4->cd(3);
227	c4->GetPad(3)->SetLogy();
228	fhDXP = (TH1F*) chist1->FindObject("fhDXPC");
229	fhDXN = (TH1F*) chist1->FindObject("fhDXNC");
230	fhDXP->SetLineColor(2);
231	fhDXN->SetLineColor(4);
232
233	fhDXP->Draw();
234	fhDXN->Draw("same");
235
236	// pos/neg z
237	c4->cd(4);
238	c4->GetPad(4)->SetLogy();
239
240	fhDXZP = (TH1F*) chist1->FindObject("fhDXPZ");
241	fhDXZN = (TH1F*) chist1->FindObject("fhDXNZ");
242	fhDXBA = (TH1F*) chist1->FindObject("fhDX_Bad");
243
244	fhDXZP->SetLineColor(2);
245	fhDXZN->SetLineColor(4);
246	fhDXBA->SetLineColor(6);
247
248	fhDXZP->Draw();
249	fhDXZN->Draw("same");
250	fhDXBA->Draw("same");
251
252	// dY
253	TCanvas* c4a = new TCanvas("c4a", "#Delta Y", 10, 10, 1100, 800);
254	c4a->Divide(2,2);
255	// pos/neg charge
256	c4a->cd(1);
257	c4a->GetPad(1)->SetLogy();
258	fhDYP = (TH1F*) chist1->FindObject("fhDYPC");
259	fhDYN = (TH1F*) chist1->FindObject("fhDYNC");
260	fhDYP->SetLineColor(2);
261	fhDYN->SetLineColor(4);
262
263	fhDYP->Draw();
264	fhDYN->Draw("same");
265
266	// pos/neg z
267	c4a->cd(2);
268	c4a->GetPad(2)->SetLogy();
269
270	fhDYZP = (TH1F*) chist1->FindObject("fhDYPZ");
271	fhDYZN = (TH1F*) chist1->FindObject("fhDYNZ");
272	fhDYBA = (TH1F*) chist1->FindObject("fhDY_Bad");
273
274	fhDYZP->SetLineColor(2);
275	fhDYZN->SetLineColor(4);
276	fhDYBA->SetLineColor(6);
277
278	fhDYZP->Draw();
279	fhDYZN->Draw("same");
280	fhDYBA->Draw("same");
281
282	// delta Pt
283	TCanvas* c5 = new TCanvas("c5", "#Delta p_{T} (n-sigma)", 10, 10, 1100, 800);
284	c5->Divide(2,2);
285	c5->cd(1);
286
287	fpDPtSP = (TH1F*) chist1->FindObject("fpDPtSPC");
288	fpDPtSN = (TH1F*) chist1->FindObject("fpDPtSNC");
289	fpDPtSP->SetLineColor(2);
290	fpDPtSN->SetLineColor(4);
291	fpDPtSP->Draw();
292	fpDPtSN->Draw("same");
293
294
295	c5->cd(2);
296
297	fpDPtSPZ = (TH1F*) chist1->FindObject("fpDPtSPZ");
298	fpDPtSNZ = (TH1F*) chist1->FindObject("fpDPtSNZ");
299	fpDPtSPZ->SetLineColor(2);
300	fpDPtSNZ->SetLineColor(4);
301	fpDPtSPZ->Draw();
302	fpDPtSNZ->Draw("same");
303
304	TCanvas* c6 = new TCanvas("c6", "Dz vs z", 10, 10, 1100, 800);
305	fhDZvsZ = (TH2F*) chist1->FindObject("fhDZvsZ");
306	fhDZvsZ->SetXTitle("z_{in} * sign(z_{in}) * sign(z_{out}) [cm]");
307	fhDZvsZ->SetYTitle("#DeltaZ [cm]");
308
309	gStyle->SetPalette(1);
310
311	fhDZvsZ->Draw("colz");
312	}
313