[u/mrichter/AliRoot.git] / PWG4 / macros / QA / DrawCaloQA.C

DrawCaloQA(TString calo = "EMCAL", Bool_t kine = kFALSE)
{
  gSystem->Load("libSTEERBase.so");
  gSystem->Load("libESD.so");
  gSystem->Load("libAOD.so");
  gSystem->Load("libANALYSIS.so");
  gSystem->Load("libANALYSISalice.so");
  gSystem->Load("libPHOSUtils");
  gSystem->Load("libEMCALUtils");
  gSystem->Load("libPWG4PartCorrBase.so");
  gSystem->Load("libPWG4PartCorrDep.so");
  
  TFile *f = new TFile("AnalysisResults.root","read");
  TDirectoryFile *dir = (TDirectoryFile *)f->Get("CaloQA");
  printf("Get list\n");
  TList* outputList = (TList*)dir->Get("CaloQA");
  
  AliAnaCalorimeterQA* qa = new AliAnaCalorimeterQA();
  qa->AddToHistogramsName(calo+"_");
  qa->SwitchOnPlotsMaking();
  qa->SetStyleMacro("style.C");
  if(kine) qa->SwitchOnDataMC() ;
  else     qa->SwitchOffDataMC() ;
  if(calo == "EMCAL"){
    qa->SetNumberOfModules(4);
    qa->SwitchOnCalorimetersCorrelation();

  }
  else if (calo=="PHOS"){
    qa->SetNumberOfModules(3);
    qa->SwitchOffCalorimetersCorrelation();

  }

  //Set Histrograms bins and ranges
  qa->SetHistoPtRangeAndNBins(0, 5, 20) ;
  if(calo=="EMCAL"){
    qa->SetHistoPhiRangeAndNBins(1.37, 2.23, 35) ;
    qa->SetHistoEtaRangeAndNBins(-0.8, 0.8, 10) ;
  }
  else if (calo=="PHOS"){
    qa->SetHistoPhiRangeAndNBins(255*TMath::DegToRad(), 325*TMath::DegToRad(), 200) ;
    qa->SetHistoEtaRangeAndNBins(-0.13, 0.13, 160) ;
  }
  qa->SetHistoMassRangeAndNBins(0., 0.6, 300) ;
  qa->SetHistoAsymmetryRangeAndNBins(0., 1. , 25) ;
  qa->SetHistoPOverERangeAndNBins(0,10.,100);
  qa->SetHistodEdxRangeAndNBins(0.,400.,200);
  qa->SetHistodRRangeAndNBins(0.,TMath::Pi(),300);
  qa->SetHistoTimeRangeAndNBins(0.,1000,1000);
  qa->SetHistoRatioRangeAndNBins(0.,2.,100);
  qa->SetHistoVertexDistRangeAndNBins(0.,500.,100);
  qa->SetHistoNClusterCellRangeAndNBins(0,300,300);
  qa->SetHistoXRangeAndNBins(-250,100,30);
  qa->SetHistoYRangeAndNBins(370,450,20);
  qa->SetHistoZRangeAndNBins(-350,350,80);
  qa->SetHistoRRangeAndNBins(420,460,10);

  //Make the histograms
  qa->Terminate(outputList);
  
  
}
Commit	Line	Data
255b904d	1	DrawCaloQA(TString calo = "EMCAL", Bool_t kine = kFALSE)
a0bb4dc0	2	{
255b904d	3	gSystem->Load("libSTEERBase.so");
	4	gSystem->Load("libESD.so");
	5	gSystem->Load("libAOD.so");
	6	gSystem->Load("libANALYSIS.so");
	7	gSystem->Load("libANALYSISalice.so");
	8	gSystem->Load("libPHOSUtils");
	9	gSystem->Load("libEMCALUtils");
	10	gSystem->Load("libPWG4PartCorrBase.so");
	11	gSystem->Load("libPWG4PartCorrDep.so");
	12
	13	TFile *f = new TFile("AnalysisResults.root","read");
	14	TDirectoryFile dir = (TDirectoryFile )f->Get("CaloQA");
	15	printf("Get list\n");
	16	TList* outputList = (TList*)dir->Get("CaloQA");
	17
	18	AliAnaCalorimeterQA* qa = new AliAnaCalorimeterQA();
	19	qa->AddToHistogramsName(calo+"_");
	20	qa->SwitchOnPlotsMaking();
	21	qa->SetStyleMacro("style.C");
	22	if(kine) qa->SwitchOnDataMC() ;
	23	else qa->SwitchOffDataMC() ;
	24	if(calo == "EMCAL"){
	25	qa->SetNumberOfModules(4);
	26	qa->SwitchOnCalorimetersCorrelation();
	27
	28	}
	29	else if (calo=="PHOS"){
	30	qa->SetNumberOfModules(3);
	31	qa->SwitchOffCalorimetersCorrelation();
	32
	33	}
	34
	35	//Set Histrograms bins and ranges
3a75604b	36	qa->SetHistoPtRangeAndNBins(0, 5, 20) ;
255b904d	37	if(calo=="EMCAL"){
3a75604b	38	qa->SetHistoPhiRangeAndNBins(1.37, 2.23, 35) ;
3a75604b	39	qa->SetHistoEtaRangeAndNBins(-0.8, 0.8, 10) ;
255b904d	40	}
	41	else if (calo=="PHOS"){
	42	qa->SetHistoPhiRangeAndNBins(255TMath::DegToRad(), 325TMath::DegToRad(), 200) ;
	43	qa->SetHistoEtaRangeAndNBins(-0.13, 0.13, 160) ;
	44	}
	45	qa->SetHistoMassRangeAndNBins(0., 0.6, 300) ;
	46	qa->SetHistoAsymmetryRangeAndNBins(0., 1. , 25) ;
	47	qa->SetHistoPOverERangeAndNBins(0,10.,100);
	48	qa->SetHistodEdxRangeAndNBins(0.,400.,200);
	49	qa->SetHistodRRangeAndNBins(0.,TMath::Pi(),300);
	50	qa->SetHistoTimeRangeAndNBins(0.,1000,1000);
	51	qa->SetHistoRatioRangeAndNBins(0.,2.,100);
	52	qa->SetHistoVertexDistRangeAndNBins(0.,500.,100);
	53	qa->SetHistoNClusterCellRangeAndNBins(0,300,300);
3a75604b	54	qa->SetHistoXRangeAndNBins(-250,100,30);
	55	qa->SetHistoYRangeAndNBins(370,450,20);
	56	qa->SetHistoZRangeAndNBins(-350,350,80);
	57	qa->SetHistoRRangeAndNBins(420,460,10);
255b904d	58
	59	//Make the histograms
	60	qa->Terminate(outputList);
	61
	62
a0bb4dc0	63	}