[u/mrichter/AliRoot.git] / PWG4 / macros / ConfigAnalysisCalorimeterQA.C

/* $Id: $ */

//------------------------------------
// Configuration macro example:
//
// Calorimeters QA
//
// Author : Gustavo Conesa Balbastre (INFN-LNF)
//------------------------------------

AliAnaPartCorrMaker*  ConfigAnalysis()
{
	//
	// Configuration goes here
	// 
	printf("======================== \n");
	printf("ConfigAnalysis() \n");
	printf("======================== \n");
	
	
	
	//-----------------------------------------------------------  
	// Reader
	//-----------------------------------------------------------
	AliCaloTrackESDReader *reader = new AliCaloTrackESDReader();
	reader->SetDebug(-1);
	
	//Switch on or off the detectors information that you want
	reader->SwitchOnEMCAL();
	reader->SwitchOffCTS();
	reader->SwitchOnPHOS();
	reader->SwitchOnEMCALCells();
	reader->SwitchOnPHOSCells();
	
	//Min particle pT
	reader->SetEMCALPtMin(0.); 
	reader->SetPHOSPtMin(0.);
	
	
	reader->Print("");
	
	
	//---------------------------------------------------------------------
	// Analysis algorithm
	//---------------------------------------------------------------------
	
	AliFidutialCut * fidCut = new AliFidutialCut();
	fidCut->DoCTSFidutialCut(kFALSE) ;
	fidCut->DoEMCALFidutialCut(kTRUE) ;
	fidCut->DoPHOSFidutialCut(kTRUE) ;
	
	
  	AliAnaCalorimeterQA *anaEMCAL = new AliAnaCalorimeterQA();
	anaEMCAL->SetDebug(-1); //10 for lots of messages
	anaEMCAL->SetCalorimeter("EMCAL");
	anaEMCAL->SwitchOnDataMC() ;//Access MC stack and fill more histograms
	//anaEMCAL->SetStyleMacro("$MACROS/style.C") ;
	anaEMCAL->AddToHistogramsName("AnaCaloQA_EMCAL_");
	anaEMCAL->SetFidutialCut(fidCut);
	anaEMCAL->SwitchOnFidutialCut();
	anaEMCAL->Print("");
	
	AliAnaCalorimeterQA *anaPHOS = new AliAnaCalorimeterQA();
	anaPHOS->SetDebug(-1); //10 for lots of messages
	anaPHOS->SetCalorimeter("PHOS");
	anaPHOS->SwitchOnDataMC() ;//Access MC stack and fill more histograms
	//anaPHOS->SetStyleMacro("$MACROS/style.C") ;
	anaPHOS->AddToHistogramsName("AnaCaloQA_PHOS_");
	anaPHOS->SetFidutialCut(fidCut);
	anaPHOS->SwitchOnFidutialCut();
	anaPHOS->Print("");
	
	
	//---------------------------------------------------------------------
	// Set  analysis algorithm and reader
	//---------------------------------------------------------------------
	maker = new AliAnaPartCorrMaker();
	maker->SetReader(reader);//pointer to reader
	maker->AddAnalysis(anaEMCAL,1);
	maker->AddAnalysis(anaPHOS,0);
	//maker->SetAnaDebug(0)  ;
	maker->SwitchOnHistogramsMaker()  ;
	maker->SwitchOffAODsMaker()  ;
	
	maker->Print("");
	//
	printf("======================== \n");
	printf("END ConfigAnalysis() \n");
	printf("======================== \n");
	return maker ;
}
Commit	Line	Data
9725fd2a	1	/* $Id: $ */
	2
	3	//------------------------------------
	4	// Configuration macro example:
	5	//
	6	// Calorimeters QA
	7	//
	8	// Author : Gustavo Conesa Balbastre (INFN-LNF)
	9	//------------------------------------
	10
	11	AliAnaPartCorrMaker* ConfigAnalysis()
	12	{
	13	//
	14	// Configuration goes here
	15	//
	16	printf("======================== \n");
	17	printf("ConfigAnalysis() \n");
	18	printf("======================== \n");
	19
	20
	21
	22	//-----------------------------------------------------------
	23	// Reader
	24	//-----------------------------------------------------------
	25	AliCaloTrackESDReader *reader = new AliCaloTrackESDReader();
	26	reader->SetDebug(-1);
	27
	28	//Switch on or off the detectors information that you want
	29	reader->SwitchOnEMCAL();
	30	reader->SwitchOffCTS();
	31	reader->SwitchOnPHOS();
	32	reader->SwitchOnEMCALCells();
	33	reader->SwitchOnPHOSCells();
	34
	35	//Min particle pT
	36	reader->SetEMCALPtMin(0.);
	37	reader->SetPHOSPtMin(0.);
	38
	39
	40	reader->Print("");
	41
	42
	43	//---------------------------------------------------------------------
	44	// Analysis algorithm
	45	//---------------------------------------------------------------------
	46
	47	AliFidutialCut * fidCut = new AliFidutialCut();
	48	fidCut->DoCTSFidutialCut(kFALSE) ;
	49	fidCut->DoEMCALFidutialCut(kTRUE) ;
	50	fidCut->DoPHOSFidutialCut(kTRUE) ;
	51
	52
	53	AliAnaCalorimeterQA *anaEMCAL = new AliAnaCalorimeterQA();
	54	anaEMCAL->SetDebug(-1); //10 for lots of messages
	55	anaEMCAL->SetCalorimeter("EMCAL");
	56	anaEMCAL->SwitchOnDataMC() ;//Access MC stack and fill more histograms
	57	//anaEMCAL->SetStyleMacro("$MACROS/style.C") ;
	58	anaEMCAL->AddToHistogramsName("AnaCaloQA_EMCAL_");
	59	anaEMCAL->SetFidutialCut(fidCut);
	60	anaEMCAL->SwitchOnFidutialCut();
	61	anaEMCAL->Print("");
	62
	63	AliAnaCalorimeterQA *anaPHOS = new AliAnaCalorimeterQA();
	64	anaPHOS->SetDebug(-1); //10 for lots of messages
65	anaPHOS->SetCalorimeter("PHOS");
66	anaPHOS->SwitchOnDataMC() ;//Access MC stack and fill more histograms
67	//anaPHOS->SetStyleMacro("$MACROS/style.C") ;
68	anaPHOS->AddToHistogramsName("AnaCaloQA_PHOS_");
69	anaPHOS->SetFidutialCut(fidCut);
70	anaPHOS->SwitchOnFidutialCut();
71	anaPHOS->Print("");
72
73
74	//---------------------------------------------------------------------
75	// Set analysis algorithm and reader
76	//---------------------------------------------------------------------
77	maker = new AliAnaPartCorrMaker();
78	maker->SetReader(reader);//pointer to reader
79	maker->AddAnalysis(anaEMCAL,1);
80	maker->AddAnalysis(anaPHOS,0);
81	//maker->SetAnaDebug(0) ;
82	maker->SwitchOnHistogramsMaker() ;
83	maker->SwitchOffAODsMaker() ;
84
85	maker->Print("");
86	//
87	printf("======================== \n");
88	printf("END ConfigAnalysis() \n");
89	printf("======================== \n");
90	return maker ;
91	}
92