[u/mrichter/AliRoot.git] / PWGPP / VZERO / trending / MakeTrendingV0QA.C

#if !defined(__CINT__) || defined(__MAKECINT__)
#include <TError.h>
#include <TROOT.h>
#include <TKey.h>
#include <TH2.h>
#include <TF1.h>
#include <TH1.h>
#include <TFile.h>
#include <TCanvas.h>
#include <TPad.h>
#include <TStyle.h>
#include <TGrid.h>
#include <TGridResult.h>
#include <TEnv.h>
#include <TLegend.h>
#include <TMath.h>
#include <TSpectrum.h>
#include <TTree.h>

#include "AliCDBManager.h"
#include "AliCDBEntry.h"
#include "AliGRPObject.h"
#include "AliTriggerInput.h"
#include "AliTriggerConfiguration.h"
#endif


Int_t MakeTrendingV0QA(TString qafilename,Int_t runNumber,TString ocdbStorage = "raw://",Bool_t IsOnGrid = kFALSE,Bool_t canvasE = kFALSE)
{
	if (!qafilename) 
	{
     		Printf("Error - Invalid input file");
  		return 1;
	}
	gStyle->SetPalette(1);
	
	TString treePostFileName=Form("trending_%i.root",runNumber);
	if(IsOnGrid)
		TGrid::Connect("alien://");
	TFile*fin=TFile::Open(qafilename,"r");

	if(!fin)
	{
		Printf("ERROR: QA output not found. Exiting ...\n");
		return -1;
	}
	else
	{
		Printf("INFO: QA output file %s open. \n",fin->GetName());
	}
	AliCDBManager *man = AliCDBManager::Instance();
	man->SetDefaultStorage(ocdbStorage);
	man->SetRun(runNumber);
    	AliCDBEntry *entry2=0;
    	entry2 = man->Get("GRP/GRP/Data");
    	AliGRPObject* fGRPData=0;
    	if (entry2) 
    	{
      		printf("Found an AliGRPObject in GRP/GRP/Data, reading it\n");
      		fGRPData = dynamic_cast<AliGRPObject*>(entry2->GetObject()); 
     		entry2->SetOwner(0);
   	}
    	TString activeDetList(AliDAQ::ListOfTriggeredDetectors(fGRPData->GetDetectorMask()));
    	TString runType(fGRPData->GetRunType());
    	TString beamType(fGRPData->GetBeamType());
    	TString machineMode(fGRPData->GetMachineMode());
    	TString lhcState(fGRPData->GetLHCState());
    	printf("activeDetList %s\nrunType %s\nbeamType %s\nmachineMode %s\nlhcState %s\n",
    		activeDetList.Data(),runType.Data(),beamType.Data(),
    		machineMode.Data(),lhcState.Data());
    
    	time_t duration = fGRPData->GetTimeEnd() - fGRPData->GetTimeStart();
    
   	if(!activeDetList.Contains("VZERO"))
    	{ 
       		printf("RUN WITH VZERO NOT ACTIVE\n");
      		return 0;
    	}
    	if(!runType.Contains("PHYSICS"))
    	{ 
      		printf("RUN NO PHYSICS\n");
      		return 0;
    	}
    	if(duration<120)
    	{ 
      		printf("RUNS SHORTER THAN 2 MIN\n");
     		return 0;
    	}
    	Float_t TimesA=-9999.,TimesC=-9999., BB_BG=-9999.,BB_EE=-9999.,AdcA=-9999.;
    	Float_t AdcC=-9999.,MultA=-9999.,MultC=-9999.;
    	Int_t NumberVoieOff=0, numberBadOffset=0;
   	TTree *ttree=new TTree("trending","tree of trending variables");

   	ttree->Branch("run",&runNumber,"run/I");
   	ttree->Branch("TimesA",&TimesA,"BB Leading time;;Time (ns)/F");
	ttree->Branch("TimesC",&TimesC,"BB Leading time;;Time (ns)/F");
	ttree->Branch("BB_BG",&BB_BG,"Trigger ratio/F");
  	ttree->Branch("BB_EE",&BB_EE,"Trigger ratio/F");
  	ttree->Branch("AdcA" ,&AdcA ,"Average Charge/F");
  	ttree->Branch("AdcC" ,&AdcC ,"Average Charge/F");
  	ttree->Branch("MultA",&MultA,"Average number of Fired cell/F");
  	ttree->Branch("MultC",&MultC,"Average number of Fired cell/F");
  	ttree->Branch("NumberVoieOff",&NumberVoieOff,"Number of path off/I");
  	ttree->Branch("numberBadOffset",&numberBadOffset,"Number of bad offset /I");
  	
  	char v0QAdirName[20]="VZERO_Performance";
  	TDirectoryFile * v0QAdir=(TDirectoryFile*)fin->Get(v0QAdirName);
  	if(!v0QAdir)
  	{
  		printf("ERROR: VZERO QA directory not present in input file.\n");
  		return -1;
  	}
  	TList *list = (TList*)v0QAdir->Get("QAVZEROHists");
    	if(!list) 
    	{
      		cout << "ERROR: No list found" << endl;
		return -1;
    	}
    	TH2F *hEvents = (TH2F*)list->FindObject("hEvents");
    	TH1F *hAdcNoTimeA = (TH1F*)list->FindObject("hAdcNoTimeA");
    	TH1F *hAdcWithTimeA = (TH1F*)list->FindObject("hAdcWithTimeA");
   	TH1F *hAdcNoTimeC = (TH1F*)list->FindObject("hAdcNoTimeC");
    	TH1F *hAdcWithTimeC = (TH1F*)list->FindObject("hAdcWithTimeC");
    	TH2F *hadcpmtwithtime = (TH2F*)list->FindObject("hadcpmtwithtime");	
    	TH1F *htimepmtA = (TH1F*)list->FindObject("htimepmtA");
    	TH1F *htimepmtC = (TH1F*)list->FindObject("htimepmtC");
    	TH1F *hwidthA = (TH1F*)list->FindObject("hwidthA");
    	TH1F *hwidthC = (TH1F*)list->FindObject("hwidthC");
    	TH1F *hV0ampl = (TH1F*)list->FindObject("hV0ampl");
    	TH2F *htimepmt = (TH2F*)list->FindObject("htimepmt");	
   	TH2F *hwidthpmt = (TH2F*)list->FindObject("hwidthpmt");	
    	TH2F *hadcwidthA = (TH2F*)list->FindObject("hadcwidthA");	
    	TH2F *hadcwidthC = (TH2F*)list->FindObject("hadcwidthC");	
    	TH2F *hAdcTimeA = (TH2F*)list->FindObject("hAdcTimeA");	
    	TH2F *hAdcTimeC = (TH2F*)list->FindObject("hAdcTimeC");	
    	TH2F *htimecorr = (TH2F*)list->FindObject("htimecorr");	
    	TH1F *hV0A = (TH1F*)list->FindObject("hV0a");
    	TH1F *hV0C = (TH1F*)list->FindObject("hV0c");
    	TH1F *hV0multA = (TH1F*)list->FindObject("hV0multA");
    	TH1F *hV0multC = (TH1F*)list->FindObject("hV0multC");
    	TH2F* hVtxXYBB  =(TH2F*) list->FindObject("fhVtxXYBB");
    	TH1F* hVtxZBB   =(TH1F*) list->FindObject("fhVtxZBB");
    	TH2F* hVtxXYBGA =(TH2F*) list->FindObject("fhVtxXYBGA");
    	TH1F* hVtxZBGA  =(TH1F*) list->FindObject("fhVtxZBGA");
    	TH2F* hVtxXYBGC =(TH2F*) list->FindObject("fhVtxXYBGC");
    	TH1F* hVtxZBGC  =(TH1F*) list->FindObject("fhVtxZBGC");
    	
    	float BB = hEvents->GetBinContent(2,2);
    	float EE = hEvents->GetBinContent(1,1);
    	float BGA = hEvents->GetBinContent(3,2);
    	float BGC = hEvents->GetBinContent(2,3);
    	
    	if(hAdcWithTimeA->GetEntries()==0)
    		return 0;
    	
    	{
    		TSpectrum s;
        	float shiftA = 8.;
    		Int_t nPeaksFound = s.Search(htimepmtA);
    		Float_t *peaks = s.GetPositionY();
		Float_t *posiX = s.GetPositionX();
    		Float_t maxY = 0.;
    		Int_t index = -1;
	
    		for(int i=0;i<nPeaksFound;i++) 
    		{
    			if(peaks[i]>maxY && posiX[i]>0.) 
    			{
				maxY = peaks[i];
				index = i;
			}	
    		}
    		Float_t maxX = (index >= 0) ? s.GetPositionX()[index] : -11111;	
		
		TF1 *fgaus = new TF1("gausbbbb","gaus",maxX-1.,maxX+1.);
    		htimepmtA->Fit(fgaus,"","",maxX-1.,maxX+1.);
    		TimesA=fgaus->GetParameter(1)-shiftA;
    		delete fgaus;
    	}
    	{
   		TSpectrum s;
    		Int_t nPeaksFound = s.Search(htimepmtC);
		Float_t *peaks = s.GetPositionY();
		Float_t maxY = 0.;
		Int_t index = -1;
    		for(int i=0;i<nPeaksFound;i++) 
    		{
			if(peaks[i]>maxY) 
			{
			maxY = peaks[i];
			index = i;
			}
		}	
    		Float_t maxX = (index >= 0) ? s.GetPositionX()[index] : -11111;	
		TF1 *fgaus = new TF1("gausffff","gaus",maxX-1.,maxX+1.);
    		htimepmtC->Fit(fgaus,"","",maxX-1.,maxX+1.);
    		TimesC=fgaus->GetParameter(1);
    		delete fgaus;
    	}
    	if(BB) 
    	{
      		BB_BG=(BGA+BGC)/BB;
      		BB_EE=EE/BB;
    	}else
    	{
      		BB_BG=0;
      		BB_EE=0;
    	}
    	
    	MultA=hV0A->GetMean();
    	MultC=hV0C->GetMean();
    
    	AdcA=hAdcWithTimeA->GetMean();
    	AdcC=hAdcWithTimeC->GetMean();
    	
	double valBin=0;
	TH1D*hadcXFull=hadcpmtwithtime->ProjectionX("hadcXFull",1,hadcpmtwithtime->GetYaxis()->GetLast());
	TH1D*hadcX=hadcpmtwithtime->ProjectionX("hadcX",10,20);
	for(Int_t i=0;i<64;i++)
	{
		valBin=hadcXFull->GetBinContent(i+1);
		if(valBin==0)
			NumberVoieOff++;
		valBin=hadcX->GetBinContent(i+1);
		if(valBin==0)
			numberBadOffset++;
	}
	    	
    	TFile * trendFile = new TFile(treePostFileName.Data(),"recreate");
    	ttree->Fill();
	trendFile->cd();
	ttree->Write();
	trendFile->Close();
	
    	if(canvasE)
    	{
   	
    		TCanvas * cOut = new TCanvas("cOut",Form("Run %d",runNumber));
    		cOut->Divide(2,2);
    		cOut->cd(1); cOut->GetPad(1)->SetLogy();
    		hAdcNoTimeA->Draw("l");
    		hAdcWithTimeA->Draw("same"); hAdcWithTimeA->SetLineColor(2);
    
     		cOut->cd(2); cOut->GetPad(2)->SetLogy();
    		hAdcNoTimeC->Draw("l");
    		hAdcWithTimeC->Draw("same"); hAdcWithTimeC->SetLineColor(2);
    	
       		cOut->cd(3); cOut->GetPad(3)->SetLogz();
    		hadcpmtwithtime->Draw("colz");
    	
    		cOut->cd(4); cOut->GetPad(4)->SetLogz();
    		hEvents->Draw("colz text");
    	
    		cOut->Print(Form("QA_Run_%d.pdf(",runNumber));
	
	
    		cOut->cd(1); cOut->GetPad(1)->SetLogy();
    		htimepmtA->GetXaxis()->SetRangeUser(-25.,25.); htimepmtA->Draw();
    	
    		cOut->cd(2); cOut->GetPad(2)->SetLogy();
    		htimepmtC->GetXaxis()->SetRangeUser(-25.,25.); htimepmtC->Draw();
    	
       		cOut->cd(3); cOut->GetPad(3)->SetLogy();cOut->GetPad(3)->SetLogz(0);
    		hwidthA->GetXaxis()->SetRangeUser(0.,50.); hwidthA->Draw();
    	
    		cOut->cd(4); cOut->GetPad(4)->SetLogy();cOut->GetPad(4)->SetLogz(0);
    		hwidthC->GetXaxis()->SetRangeUser(0.,50.); hwidthC->Draw();
    	
    		cOut->Print(Form("QA_Run_%d.pdf",runNumber));
    		
    	
    		cOut->cd(1); cOut->GetPad(1)->SetLogy(0);cOut->GetPad(1)->SetLogz();
    		htimepmt->Draw("colz");
    	
    		cOut->cd(2); cOut->GetPad(2)->SetLogy(0);cOut->GetPad(2)->SetLogz();
    		hwidthpmt->GetYaxis()->SetRangeUser(0.,50.); hwidthpmt->Draw("colz");
    		
    		cOut->cd(3); cOut->GetPad(3)->SetLogy(0);cOut->GetPad(3)->SetLogz();
    		hadcwidthA->GetYaxis()->SetRangeUser(0.,50.); hadcwidthA->Draw("colz");
    	
    		cOut->cd(4); cOut->GetPad(4)->SetLogy(0);cOut->GetPad(4)->SetLogz();
    		hadcwidthC->GetYaxis()->SetRangeUser(0.,50.); hadcwidthC->Draw("colz");
    	
    		cOut->Print(Form("QA_Run_%d.pdf",runNumber));
    		
	
    		cOut->cd(1); cOut->GetPad(1)->SetLogy(0);cOut->GetPad(1)->SetLogz();
    		hAdcTimeA->Draw("colz");
	
    		cOut->cd(2); cOut->GetPad(2)->SetLogy(0);cOut->GetPad(2)->SetLogz();
    		hAdcTimeC->Draw("colz");
 	
    		cOut->cd(3); cOut->GetPad(3)->SetLogy(); cOut->GetPad(3)->SetLogz(0);
    		hV0ampl->Draw();
	
    		cOut->cd(4); cOut->GetPad(4)->SetLogy(0); cOut->GetPad(4)->SetLogz(0);
    		htimecorr->Draw("colz");
    	
    		cOut->Print(Form("QA_Run_%d.pdf",runNumber));
	
		
    		cOut->cd(1);  cOut->GetPad(1)->SetLogy(1);cOut->GetPad(1)->SetLogz(0);
    		hV0A->GetXaxis()->SetRangeUser(0.,33.);hV0A->Draw();
    	
    		cOut->cd(2); cOut->GetPad(2)->SetLogy(1);cOut->GetPad(2)->SetLogz(0);
    		hV0C->GetXaxis()->SetRangeUser(0.,33.);hV0C->Draw();
    	
    		cOut->cd(3); cOut->GetPad(3)->SetLogy(); cOut->GetPad(3)->SetLogz(0);
    		hV0multA->Draw();
		
    		cOut->cd(4); cOut->GetPad(4)->SetLogy(); cOut->GetPad(3)->SetLogz(0);
    		hV0multC->Draw();
    	
    		cOut->Print(Form("QA_Run_%d.pdf",runNumber));
	
		
    		cOut->Clear();
    		cOut->Divide(2,3);
    	
    		cOut->cd(1);  cOut->GetPad(1)->SetLogy(0);cOut->GetPad(1)->SetLogz(1);
    		hVtxXYBB->Draw("colz");
    	
    		cOut->cd(2); cOut->GetPad(2)->SetLogy(1);cOut->GetPad(2)->SetLogz(0);
    		hVtxZBB->Draw();
    	
      		cOut->cd(3); cOut->GetPad(3)->SetLogy(0); cOut->GetPad(3)->SetLogz(1);
    		hVtxXYBGA->Draw("colz");
    	
    		cOut->cd(4); cOut->GetPad(4)->SetLogy(); cOut->GetPad(3)->SetLogz(0);
    		hVtxZBGA->Draw();
    	
    		cOut->cd(5); cOut->GetPad(5)->SetLogy(0); cOut->GetPad(5)->SetLogz(1);
    		hVtxXYBGC->Draw("colz");
    	
    		cOut->cd(6); cOut->GetPad(6)->SetLogy(); cOut->GetPad(6)->SetLogz(0);
    		hVtxZBGC->Draw();
    		
    		cOut->Print(Form("QA_Run_%d.pdf)",runNumber));
    		delete cOut;
    					
 	}
	delete v0QAdir;
 	return 0;
}
Commit	Line	Data
ae186c93	1	#if !defined(__CINT__) \|\| defined(__MAKECINT__)
	2	#include <TError.h>
	3	#include <TROOT.h>
	4	#include <TKey.h>
	5	#include <TH2.h>
	6	#include <TF1.h>
	7	#include <TH1.h>
	8	#include <TFile.h>
	9	#include <TCanvas.h>
	10	#include <TPad.h>
	11	#include <TStyle.h>
	12	#include <TGrid.h>
	13	#include <TGridResult.h>
	14	#include <TEnv.h>
	15	#include <TLegend.h>
	16	#include <TMath.h>
	17	#include <TSpectrum.h>
	18	#include <TTree.h>
	19
	20	#include "AliCDBManager.h"
	21	#include "AliCDBEntry.h"
	22	#include "AliGRPObject.h"
	23	#include "AliTriggerInput.h"
	24	#include "AliTriggerConfiguration.h"
	25	#endif
	26
	27
	28	Int_t MakeTrendingV0QA(TString qafilename,Int_t runNumber,TString ocdbStorage = "raw://",Bool_t IsOnGrid = kFALSE,Bool_t canvasE = kFALSE)
	29	{
	30	if (!qafilename)
	31	{
	32	Printf("Error - Invalid input file");
	33	return 1;
	34	}
	35	gStyle->SetPalette(1);
	36
	37	TString treePostFileName=Form("trending_%i.root",runNumber);
	38	if(IsOnGrid)
	39	TGrid::Connect("alien://");
	40	TFile*fin=TFile::Open(qafilename,"r");
	41
	42	if(!fin)
	43	{
	44	Printf("ERROR: QA output not found. Exiting ...\n");
	45	return -1;
	46	}
	47	else
	48	{
	49	Printf("INFO: QA output file %s open. \n",fin->GetName());
	50	}
	51	AliCDBManager *man = AliCDBManager::Instance();
	52	man->SetDefaultStorage(ocdbStorage);
	53	man->SetRun(runNumber);
	54	AliCDBEntry *entry2=0;
	55	entry2 = man->Get("GRP/GRP/Data");
	56	AliGRPObject* fGRPData=0;
	57	if (entry2)
	58	{
	59	printf("Found an AliGRPObject in GRP/GRP/Data, reading it\n");
	60	fGRPData = dynamic_cast<AliGRPObject*>(entry2->GetObject());
	61	entry2->SetOwner(0);
	62	}
	63	TString activeDetList(AliDAQ::ListOfTriggeredDetectors(fGRPData->GetDetectorMask()));
	64	TString runType(fGRPData->GetRunType());
65	TString beamType(fGRPData->GetBeamType());
66	TString machineMode(fGRPData->GetMachineMode());
67	TString lhcState(fGRPData->GetLHCState());
68	printf("activeDetList %s\nrunType %s\nbeamType %s\nmachineMode %s\nlhcState %s\n",
69	activeDetList.Data(),runType.Data(),beamType.Data(),
70	machineMode.Data(),lhcState.Data());
71
72	time_t duration = fGRPData->GetTimeEnd() - fGRPData->GetTimeStart();
73
74	if(!activeDetList.Contains("VZERO"))
75	{
76	printf("RUN WITH VZERO NOT ACTIVE\n");
77	return 0;
78	}
79	if(!runType.Contains("PHYSICS"))
80	{
81	printf("RUN NO PHYSICS\n");
82	return 0;
83	}
84	if(duration<120)
85	{
86	printf("RUNS SHORTER THAN 2 MIN\n");
87	return 0;
88	}
89	Float_t TimesA=-9999.,TimesC=-9999., BB_BG=-9999.,BB_EE=-9999.,AdcA=-9999.;
90	Float_t AdcC=-9999.,MultA=-9999.,MultC=-9999.;
91	Int_t NumberVoieOff=0, numberBadOffset=0;
92	TTree *ttree=new TTree("trending","tree of trending variables");
93
94	ttree->Branch("run",&runNumber,"run/I");
95	ttree->Branch("TimesA",&TimesA,"BB Leading time;;Time (ns)/F");
96	ttree->Branch("TimesC",&TimesC,"BB Leading time;;Time (ns)/F");
97	ttree->Branch("BB_BG",&BB_BG,"Trigger ratio/F");
98	ttree->Branch("BB_EE",&BB_EE,"Trigger ratio/F");
99	ttree->Branch("AdcA" ,&AdcA ,"Average Charge/F");
100	ttree->Branch("AdcC" ,&AdcC ,"Average Charge/F");
101	ttree->Branch("MultA",&MultA,"Average number of Fired cell/F");
102	ttree->Branch("MultC",&MultC,"Average number of Fired cell/F");
103	ttree->Branch("NumberVoieOff",&NumberVoieOff,"Number of path off/I");
104	ttree->Branch("numberBadOffset",&numberBadOffset,"Number of bad offset /I");
105
106	char v0QAdirName[20]="VZERO_Performance";
107	TDirectoryFile * v0QAdir=(TDirectoryFile*)fin->Get(v0QAdirName);
108	if(!v0QAdir)
109	{
110	printf("ERROR: VZERO QA directory not present in input file.\n");
111	return -1;
112	}
113	TList list = (TList)v0QAdir->Get("QAVZEROHists");
114	if(!list)
115	{
116	cout << "ERROR: No list found" << endl;
117	return -1;
118	}
119	TH2F hEvents = (TH2F)list->FindObject("hEvents");
120	TH1F hAdcNoTimeA = (TH1F)list->FindObject("hAdcNoTimeA");
121	TH1F hAdcWithTimeA = (TH1F)list->FindObject("hAdcWithTimeA");
122	TH1F hAdcNoTimeC = (TH1F)list->FindObject("hAdcNoTimeC");
123	TH1F hAdcWithTimeC = (TH1F)list->FindObject("hAdcWithTimeC");
124	TH2F hadcpmtwithtime = (TH2F)list->FindObject("hadcpmtwithtime");
125	TH1F htimepmtA = (TH1F)list->FindObject("htimepmtA");
126	TH1F htimepmtC = (TH1F)list->FindObject("htimepmtC");
127	TH1F hwidthA = (TH1F)list->FindObject("hwidthA");
128	TH1F hwidthC = (TH1F)list->FindObject("hwidthC");
129	TH1F hV0ampl = (TH1F)list->FindObject("hV0ampl");
130	TH2F htimepmt = (TH2F)list->FindObject("htimepmt");
131	TH2F hwidthpmt = (TH2F)list->FindObject("hwidthpmt");
132	TH2F hadcwidthA = (TH2F)list->FindObject("hadcwidthA");
133	TH2F hadcwidthC = (TH2F)list->FindObject("hadcwidthC");
134	TH2F hAdcTimeA = (TH2F)list->FindObject("hAdcTimeA");
135	TH2F hAdcTimeC = (TH2F)list->FindObject("hAdcTimeC");
136	TH2F htimecorr = (TH2F)list->FindObject("htimecorr");
137	TH1F hV0A = (TH1F)list->FindObject("hV0a");
138	TH1F hV0C = (TH1F)list->FindObject("hV0c");
139	TH1F hV0multA = (TH1F)list->FindObject("hV0multA");
140	TH1F hV0multC = (TH1F)list->FindObject("hV0multC");
141	TH2F* hVtxXYBB =(TH2F*) list->FindObject("fhVtxXYBB");
142	TH1F* hVtxZBB =(TH1F*) list->FindObject("fhVtxZBB");
143	TH2F* hVtxXYBGA =(TH2F*) list->FindObject("fhVtxXYBGA");
144	TH1F* hVtxZBGA =(TH1F*) list->FindObject("fhVtxZBGA");
145	TH2F* hVtxXYBGC =(TH2F*) list->FindObject("fhVtxXYBGC");
146	TH1F* hVtxZBGC =(TH1F*) list->FindObject("fhVtxZBGC");
147
148	float BB = hEvents->GetBinContent(2,2);
149	float EE = hEvents->GetBinContent(1,1);
150	float BGA = hEvents->GetBinContent(3,2);
151	float BGC = hEvents->GetBinContent(2,3);
152
153	if(hAdcWithTimeA->GetEntries()==0)
154	return 0;
155
156	{
157	TSpectrum s;
158	float shiftA = 8.;
159	Int_t nPeaksFound = s.Search(htimepmtA);
160	Float_t *peaks = s.GetPositionY();
161	Float_t *posiX = s.GetPositionX();
162	Float_t maxY = 0.;
163	Int_t index = -1;
164
165	for(int i=0;i<nPeaksFound;i++)
166	{
167	if(peaks[i]>maxY && posiX[i]>0.)
168	{
169	maxY = peaks[i];
170	index = i;
171	}
172	}
173	Float_t maxX = (index >= 0) ? s.GetPositionX()[index] : -11111;
174
175	TF1 *fgaus = new TF1("gausbbbb","gaus",maxX-1.,maxX+1.);
176	htimepmtA->Fit(fgaus,"","",maxX-1.,maxX+1.);
177	TimesA=fgaus->GetParameter(1)-shiftA;
178	delete fgaus;
179	}
180	{
181	TSpectrum s;
182	Int_t nPeaksFound = s.Search(htimepmtC);
183	Float_t *peaks = s.GetPositionY();
184	Float_t maxY = 0.;
185	Int_t index = -1;
186	for(int i=0;i<nPeaksFound;i++)
187	{
188	if(peaks[i]>maxY)
189	{
190	maxY = peaks[i];
191	index = i;
192	}
193	}
194	Float_t maxX = (index >= 0) ? s.GetPositionX()[index] : -11111;
195	TF1 *fgaus = new TF1("gausffff","gaus",maxX-1.,maxX+1.);
196	htimepmtC->Fit(fgaus,"","",maxX-1.,maxX+1.);
197	TimesC=fgaus->GetParameter(1);
198	delete fgaus;
199	}
200	if(BB)
201	{
202	BB_BG=(BGA+BGC)/BB;
203	BB_EE=EE/BB;
204	}else
205	{
206	BB_BG=0;
207	BB_EE=0;
208	}
209
210	MultA=hV0A->GetMean();
211	MultC=hV0C->GetMean();
212
213	AdcA=hAdcWithTimeA->GetMean();
214	AdcC=hAdcWithTimeC->GetMean();
215
216	double valBin=0;
217	TH1D*hadcXFull=hadcpmtwithtime->ProjectionX("hadcXFull",1,hadcpmtwithtime->GetYaxis()->GetLast());
218	TH1D*hadcX=hadcpmtwithtime->ProjectionX("hadcX",10,20);
219	for(Int_t i=0;i<64;i++)
220	{
221	valBin=hadcXFull->GetBinContent(i+1);
222	if(valBin==0)
223	NumberVoieOff++;
224	valBin=hadcX->GetBinContent(i+1);
225	if(valBin==0)
226	numberBadOffset++;
227	}
228
229	TFile * trendFile = new TFile(treePostFileName.Data(),"recreate");
230	ttree->Fill();
231	trendFile->cd();
232	ttree->Write();
233	trendFile->Close();
234
235	if(canvasE)
236	{
237
238	TCanvas * cOut = new TCanvas("cOut",Form("Run %d",runNumber));
239	cOut->Divide(2,2);
240	cOut->cd(1); cOut->GetPad(1)->SetLogy();
241	hAdcNoTimeA->Draw("l");
242	hAdcWithTimeA->Draw("same"); hAdcWithTimeA->SetLineColor(2);
243
244	cOut->cd(2); cOut->GetPad(2)->SetLogy();
245	hAdcNoTimeC->Draw("l");
246	hAdcWithTimeC->Draw("same"); hAdcWithTimeC->SetLineColor(2);
247
248	cOut->cd(3); cOut->GetPad(3)->SetLogz();
249	hadcpmtwithtime->Draw("colz");
250
251	cOut->cd(4); cOut->GetPad(4)->SetLogz();
252	hEvents->Draw("colz text");
253
254	cOut->Print(Form("QA_Run_%d.pdf(",runNumber));
255
256
257	cOut->cd(1); cOut->GetPad(1)->SetLogy();
258	htimepmtA->GetXaxis()->SetRangeUser(-25.,25.); htimepmtA->Draw();
259
260	cOut->cd(2); cOut->GetPad(2)->SetLogy();
261	htimepmtC->GetXaxis()->SetRangeUser(-25.,25.); htimepmtC->Draw();
262
263	cOut->cd(3); cOut->GetPad(3)->SetLogy();cOut->GetPad(3)->SetLogz(0);
264	hwidthA->GetXaxis()->SetRangeUser(0.,50.); hwidthA->Draw();
265
266	cOut->cd(4); cOut->GetPad(4)->SetLogy();cOut->GetPad(4)->SetLogz(0);
267	hwidthC->GetXaxis()->SetRangeUser(0.,50.); hwidthC->Draw();
268
269	cOut->Print(Form("QA_Run_%d.pdf",runNumber));
270
271
272	cOut->cd(1); cOut->GetPad(1)->SetLogy(0);cOut->GetPad(1)->SetLogz();
273	htimepmt->Draw("colz");
274
275	cOut->cd(2); cOut->GetPad(2)->SetLogy(0);cOut->GetPad(2)->SetLogz();
276	hwidthpmt->GetYaxis()->SetRangeUser(0.,50.); hwidthpmt->Draw("colz");
277
278	cOut->cd(3); cOut->GetPad(3)->SetLogy(0);cOut->GetPad(3)->SetLogz();
279	hadcwidthA->GetYaxis()->SetRangeUser(0.,50.); hadcwidthA->Draw("colz");
280
281	cOut->cd(4); cOut->GetPad(4)->SetLogy(0);cOut->GetPad(4)->SetLogz();
282	hadcwidthC->GetYaxis()->SetRangeUser(0.,50.); hadcwidthC->Draw("colz");
283
284	cOut->Print(Form("QA_Run_%d.pdf",runNumber));
285
286
287	cOut->cd(1); cOut->GetPad(1)->SetLogy(0);cOut->GetPad(1)->SetLogz();
288	hAdcTimeA->Draw("colz");
289
290	cOut->cd(2); cOut->GetPad(2)->SetLogy(0);cOut->GetPad(2)->SetLogz();
291	hAdcTimeC->Draw("colz");
292
293	cOut->cd(3); cOut->GetPad(3)->SetLogy(); cOut->GetPad(3)->SetLogz(0);
294	hV0ampl->Draw();
295
296	cOut->cd(4); cOut->GetPad(4)->SetLogy(0); cOut->GetPad(4)->SetLogz(0);
297	htimecorr->Draw("colz");
298
299	cOut->Print(Form("QA_Run_%d.pdf",runNumber));
300
301
302
303	cOut->cd(1); cOut->GetPad(1)->SetLogy(1);cOut->GetPad(1)->SetLogz(0);
304	hV0A->GetXaxis()->SetRangeUser(0.,33.);hV0A->Draw();
305
306	cOut->cd(2); cOut->GetPad(2)->SetLogy(1);cOut->GetPad(2)->SetLogz(0);
307	hV0C->GetXaxis()->SetRangeUser(0.,33.);hV0C->Draw();
308
309	cOut->cd(3); cOut->GetPad(3)->SetLogy(); cOut->GetPad(3)->SetLogz(0);
310	hV0multA->Draw();
311
312	cOut->cd(4); cOut->GetPad(4)->SetLogy(); cOut->GetPad(3)->SetLogz(0);
313	hV0multC->Draw();
314
315	cOut->Print(Form("QA_Run_%d.pdf",runNumber));
316
317
318	cOut->Clear();
319	cOut->Divide(2,3);
320
321	cOut->cd(1); cOut->GetPad(1)->SetLogy(0);cOut->GetPad(1)->SetLogz(1);
322	hVtxXYBB->Draw("colz");
323
324	cOut->cd(2); cOut->GetPad(2)->SetLogy(1);cOut->GetPad(2)->SetLogz(0);
325	hVtxZBB->Draw();
326
327	cOut->cd(3); cOut->GetPad(3)->SetLogy(0); cOut->GetPad(3)->SetLogz(1);
328	hVtxXYBGA->Draw("colz");
329
330	cOut->cd(4); cOut->GetPad(4)->SetLogy(); cOut->GetPad(3)->SetLogz(0);
331	hVtxZBGA->Draw();
332
333	cOut->cd(5); cOut->GetPad(5)->SetLogy(0); cOut->GetPad(5)->SetLogz(1);
334	hVtxXYBGC->Draw("colz");
335
336	cOut->cd(6); cOut->GetPad(6)->SetLogy(); cOut->GetPad(6)->SetLogz(0);
337	hVtxZBGC->Draw();
338
339	cOut->Print(Form("QA_Run_%d.pdf)",runNumber));
340	delete cOut;
341
342	}
343	delete v0QAdir;
344	return 0;
345	}