[u/mrichter/AliRoot.git] / PWGGA / GammaConv / macros / AddTask_GammaConvV1_PbPb2.C

void AddTask_GammaConvV1_PbPb2(  Int_t trainConfig = 1,  //change different set of cuts
                              Bool_t isMC   = kFALSE, //run MC 
                              Int_t enableQAMesonTask = 0, //enable QA in AliAnalysisTaskGammaConvV1
                              Int_t enableQAPhotonTask = 0, // enable additional QA task
                              TString fileNameInputForWeighting = "MCSpectraInput.root", // path to file for weigting input
                              Bool_t doWeighting = kFALSE,  //enable Weighting
                              TString cutnumberAODBranch = "1000000060084000001500000" 
                           ) {

   // ================= Load Librariers =================================
   gSystem->Load("libCore.so");  
   gSystem->Load("libTree.so");
   gSystem->Load("libGeom.so");
   gSystem->Load("libVMC.so");
   gSystem->Load("libPhysics.so");
   gSystem->Load("libMinuit");
   gSystem->Load("libSTEERBase");
   gSystem->Load("libESD");
   gSystem->Load("libAOD");
   gSystem->Load("libANALYSIS");
   gSystem->Load("libANALYSISalice");  
   gSystem->Load("libPWGGAGammaConv.so");
   gSystem->Load("libCDB.so");
   gSystem->Load("libSTEER.so");
   gSystem->Load("libSTEERBase.so");
   gSystem->Load("libTENDER.so");
   gSystem->Load("libTENDERSupplies.so");
      
   // ================== GetAnalysisManager ===============================
   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
   if (!mgr) {
      Error(Form("AddTask_GammaConvV1_%i",trainConfig), "No analysis manager found.");
      return ;
   }

   // ================== GetInputEventHandler =============================
   AliVEventHandler *inputHandler=mgr->GetInputEventHandler();
   
   //========= Add PID Reponse to ANALYSIS manager ====
   if(!(AliPIDResponse*)mgr->GetTask("PIDResponseTask")){
      gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
      AddTaskPIDResponse(isMC);
   }
   
   //=========  Set Cutnumber for V0Reader ================================
   TString cutnumber = "1000000000084001001500000000"; 
   AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
   //========= Add V0 Reader to  ANALYSIS manager if not yet existent =====
   if( !(AliV0ReaderV1*)mgr->GetTask("V0ReaderV1") ){
      AliV0ReaderV1 *fV0ReaderV1 = new AliV0ReaderV1("V0ReaderV1");
      
      fV0ReaderV1->SetUseOwnXYZCalculation(kTRUE);
      fV0ReaderV1->SetCreateAODs(kFALSE);// AOD Output
      fV0ReaderV1->SetUseAODConversionPhoton(kTRUE);

      if (!mgr) {
         Error("AddTask_V0ReaderV1", "No analysis manager found.");
         return;
      }

      // Set AnalysisCut Number
      AliConversionCuts *fCuts=NULL;
      if(cutnumber!=""){
         fCuts= new AliConversionCuts(cutnumber.Data(),cutnumber.Data());
         fCuts->SetPreSelectionCutFlag(kTRUE);
         if(fCuts->InitializeCutsFromCutString(cutnumber.Data())){
            fV0ReaderV1->SetConversionCuts(fCuts);
            fCuts->SetFillCutHistograms("",kTRUE);
         }
      }
      
      if(inputHandler->IsA()==AliAODInputHandler::Class()){
      // AOD mode
         fV0ReaderV1->SetDeltaAODBranchName(Form("GammaConv_%s_gamma",cutnumberAODBranch.Data()));
      }
      fV0ReaderV1->Init();

      AliLog::SetGlobalLogLevel(AliLog::kInfo);

      //connect input V0Reader
      mgr->AddTask(fV0ReaderV1);
      mgr->ConnectInput(fV0ReaderV1,0,cinput);

   }

   //================================================
   //========= Add task to the ANALYSIS manager =====
   //================================================
   AliAnalysisTaskGammaConvV1 *task=NULL;
   task= new AliAnalysisTaskGammaConvV1(Form("GammaConvV1_%i",trainConfig));
   task->SetIsHeavyIon(1);
   task->SetIsMC(isMC);
   // Cut Numbers to use in Analysis
   Int_t numberOfCuts = 1;

   TString *cutarray = new TString[numberOfCuts];
   TString *mesonCutArray = new TString[numberOfCuts];

   if (trainConfig == 1){ 
      cutarray[ 0] = "6010001042092970023220000000"; mesonCutArray[ 0] = "01522045009000"; 
   } else if (trainConfig == 2) { 
      cutarray[ 0] = "6120001042092970023220000000"; mesonCutArray[ 0] = "01522045009000"; 
   } else if (trainConfig == 3) { 
      cutarray[ 0] = "5010001042092970023220000000"; mesonCutArray[ 0] = "01522045009000"; 
   } else if (trainConfig == 4) { 
      cutarray[ 0] = "5020001042092970023220000000"; mesonCutArray[ 0] = "01522045009000";    
   } else if (trainConfig == 5) { 
      cutarray[ 0] = "5120001042092970023220000000"; mesonCutArray[ 0] = "01522045009000";    
   } else if (trainConfig == 6) { 
      cutarray[ 0] = "5240001042092970023220000000"; mesonCutArray[ 0] = "01522045009000";       
   } else if (trainConfig == 7) {    
      cutarray[ 0] = "5460001042092970023220000000"; mesonCutArray[ 0] = "01522065009000"; 
   } else if (trainConfig == 8) {    
      cutarray[ 0] = "5480001042092970023220000000"; mesonCutArray[ 0] = "01522065009000";    
   } else if (trainConfig == 9) {    
      cutarray[ 0] = "5450001042092970023220000000"; mesonCutArray[ 0] = "01522065009000"; 
   } else if (trainConfig == 10) { 
      cutarray[ 0] = "5560001042092970023220000000"; mesonCutArray[ 0] = "01522065009000";
   } else if (trainConfig == 11) { 
      cutarray[ 0] = "5680001042092970023220000000"; mesonCutArray[ 0] = "01522065009000";    
   } else if (trainConfig == 12) { 
      cutarray[ 0] = "5670001042092970023220000000"; mesonCutArray[ 0] = "01522065009000"; 
   } else if (trainConfig == 13) { 
      cutarray[ 0] = "5780001042092970023220000000"; mesonCutArray[ 0] = "01522065009000"; 
   } else if (trainConfig == 14) { 
      cutarray[ 0] = "4690001042092970023220000000"; mesonCutArray[ 0] = "01522065009000";
   } else if (trainConfig == 15) { 
      cutarray[ 0] = "5890001042092970023220000000"; mesonCutArray[ 0] = "01522065009000";    
   } else  if (trainConfig == 16){ 
      cutarray[ 0] = "6010001002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; 
   } else if (trainConfig == 17) { 
      cutarray[ 0] = "6120001002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; 
   } else if (trainConfig == 18) { 
      cutarray[ 0] = "5010001002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; 
   } else if (trainConfig == 19) { 
      cutarray[ 0] = "5020001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";    
   } else if (trainConfig == 20) { 
      cutarray[ 0] = "5120001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";    
   } else if (trainConfig == 21) { 
      cutarray[ 0] = "5240001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";       
   } else if (trainConfig == 22) {    
      cutarray[ 0] = "5460001002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; 
   } else if (trainConfig == 23) {    
      cutarray[ 0] = "5480001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";    
   } else if (trainConfig == 24) {    
      cutarray[ 0] = "5450001002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; 
   } else if (trainConfig == 25) { 
      cutarray[ 0] = "5560001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
   } else if (trainConfig == 26) { 
      cutarray[ 0] = "5680001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";    
   } else if (trainConfig == 27) { 
      cutarray[ 0] = "5670001002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; 
   } else if (trainConfig == 28) { 
      cutarray[ 0] = "5780001002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; 
   } else if (trainConfig == 29) { 
      cutarray[ 0] = "4690001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
   } else if (trainConfig == 30) { 
      cutarray[ 0] = "5890001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";    
   } else if (trainConfig == 31) { 
      cutarray[ 0] = "5080001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";    
   } else if (trainConfig == 32) { 
      cutarray[ 0] = "5250001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";       
   } else if (trainConfig == 33) { 
      cutarray[ 0] = "5350001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";       
   } else if (trainConfig == 34) { 
      cutarray[ 0] = "5450001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";       
   } else if (trainConfig == 35) { 
      cutarray[ 0] = "5340001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";       
   } else if (trainConfig == 36) { 
      cutarray[ 0] = "5230001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";       
	} else {
      Error(Form("GammaConvV1_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
      return;
   }

   TList *ConvCutList = new TList();
   TList *MesonCutList = new TList();

   TList *HeaderList = new TList();
   TObjString *Header1 = new TObjString("BOX");
   HeaderList->Add(Header1);
//    TObjString *Header3 = new TObjString("eta_2");
//    HeaderList->Add(Header3);
   
   ConvCutList->SetOwner(kTRUE);
   AliConversionCuts **analysisCuts = new AliConversionCuts*[numberOfCuts];
   MesonCutList->SetOwner(kTRUE);
   AliConversionMesonCuts **analysisMesonCuts = new AliConversionMesonCuts*[numberOfCuts];

   for(Int_t i = 0; i<numberOfCuts; i++){
      analysisCuts[i] = new AliConversionCuts();
      analysisCuts[i]->InitializeCutsFromCutString(cutarray[i].Data());
      ConvCutList->Add(analysisCuts[i]);

      analysisCuts[i]->SetFillCutHistograms("",kFALSE);
      analysisMesonCuts[i] = new AliConversionMesonCuts();
      analysisMesonCuts[i]->InitializeCutsFromCutString(mesonCutArray[i].Data());
      MesonCutList->Add(analysisMesonCuts[i]);
      analysisMesonCuts[i]->SetFillCutHistograms("");
      analysisCuts[i]->SetAcceptedHeader(HeaderList);
   }

   task->SetConversionCutList(numberOfCuts,ConvCutList);
   task->SetMesonCutList(numberOfCuts,MesonCutList);
   task->SetMoveParticleAccordingToVertex(kTRUE);
   task->SetDoMesonAnalysis(kTRUE);
   task->SetDoMesonQA(enableQAMesonTask); //Attention new switch for Pi0 QA
   task->SetDoPhotonQA(enableQAPhotonTask);  //Attention new switch small for Photon QA

   //connect containers
   AliAnalysisDataContainer *coutput =
      mgr->CreateContainer(Form("GammaConvV1_%i",trainConfig), TList::Class(),
                           AliAnalysisManager::kOutputContainer,Form("GammaConvV1_%i.root",trainConfig));

   mgr->AddTask(task);
   mgr->ConnectInput(task,0,cinput);
   mgr->ConnectOutput(task,1,coutput);

   return;

}
Commit	Line	Data
39e8c2ce	1	void AddTask_GammaConvV1_PbPb2( Int_t trainConfig = 1, //change different set of cuts
39e8c2ce	2	Bool_t isMC = kFALSE, //run MC
ae4f2cfb	3	Int_t enableQAMesonTask = 0, //enable QA in AliAnalysisTaskGammaConvV1
ae4f2cfb	4	Int_t enableQAPhotonTask = 0, // enable additional QA task
39e8c2ce	5	TString fileNameInputForWeighting = "MCSpectraInput.root", // path to file for weigting input
	6	Bool_t doWeighting = kFALSE, //enable Weighting
	7	TString cutnumberAODBranch = "1000000060084000001500000"
	8	) {
	9
	10	// ================= Load Librariers =================================
	11	gSystem->Load("libCore.so");
	12	gSystem->Load("libTree.so");
	13	gSystem->Load("libGeom.so");
	14	gSystem->Load("libVMC.so");
	15	gSystem->Load("libPhysics.so");
	16	gSystem->Load("libMinuit");
	17	gSystem->Load("libSTEERBase");
	18	gSystem->Load("libESD");
	19	gSystem->Load("libAOD");
	20	gSystem->Load("libANALYSIS");
	21	gSystem->Load("libANALYSISalice");
	22	gSystem->Load("libPWGGAGammaConv.so");
	23	gSystem->Load("libCDB.so");
	24	gSystem->Load("libSTEER.so");
	25	gSystem->Load("libSTEERBase.so");
	26	gSystem->Load("libTENDER.so");
	27	gSystem->Load("libTENDERSupplies.so");
	28
	29	// ================== GetAnalysisManager ===============================
	30	AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
	31	if (!mgr) {
	32	Error(Form("AddTask_GammaConvV1_%i",trainConfig), "No analysis manager found.");
	33	return ;
	34	}
	35
	36	// ================== GetInputEventHandler =============================
	37	AliVEventHandler *inputHandler=mgr->GetInputEventHandler();
	38
	39	//========= Add PID Reponse to ANALYSIS manager ====
	40	if(!(AliPIDResponse*)mgr->GetTask("PIDResponseTask")){
	41	gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
	42	AddTaskPIDResponse(isMC);
	43	}
	44
	45	//========= Set Cutnumber for V0Reader ================================
ae4f2cfb	46	TString cutnumber = "1000000000084001001500000000";
39e8c2ce	47	AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
	48	//========= Add V0 Reader to ANALYSIS manager if not yet existent =====
	49	if( !(AliV0ReaderV1*)mgr->GetTask("V0ReaderV1") ){
	50	AliV0ReaderV1 *fV0ReaderV1 = new AliV0ReaderV1("V0ReaderV1");
	51
	52	fV0ReaderV1->SetUseOwnXYZCalculation(kTRUE);
	53	fV0ReaderV1->SetCreateAODs(kFALSE);// AOD Output
	54	fV0ReaderV1->SetUseAODConversionPhoton(kTRUE);
	55
	56	if (!mgr) {
	57	Error("AddTask_V0ReaderV1", "No analysis manager found.");
	58	return;
	59	}
	60
	61	// Set AnalysisCut Number
	62	AliConversionCuts *fCuts=NULL;
	63	if(cutnumber!=""){
	64	fCuts= new AliConversionCuts(cutnumber.Data(),cutnumber.Data());
	65	fCuts->SetPreSelectionCutFlag(kTRUE);
	66	if(fCuts->InitializeCutsFromCutString(cutnumber.Data())){
	67	fV0ReaderV1->SetConversionCuts(fCuts);
	68	fCuts->SetFillCutHistograms("",kTRUE);
	69	}
	70	}
	71
	72	if(inputHandler->IsA()==AliAODInputHandler::Class()){
	73	// AOD mode
	74	fV0ReaderV1->SetDeltaAODBranchName(Form("GammaConv_%s_gamma",cutnumberAODBranch.Data()));
	75	}
	76	fV0ReaderV1->Init();
	77
	78	AliLog::SetGlobalLogLevel(AliLog::kInfo);
	79
	80	//connect input V0Reader
	81	mgr->AddTask(fV0ReaderV1);
	82	mgr->ConnectInput(fV0ReaderV1,0,cinput);
	83
	84	}
	85
	86	//================================================
	87	//========= Add task to the ANALYSIS manager =====
	88	//================================================
	89	AliAnalysisTaskGammaConvV1 *task=NULL;
	90	task= new AliAnalysisTaskGammaConvV1(Form("GammaConvV1_%i",trainConfig));
	91	task->SetIsHeavyIon(1);
	92	task->SetIsMC(isMC);
	93	// Cut Numbers to use in Analysis
	94	Int_t numberOfCuts = 1;
	95
	96	TString *cutarray = new TString[numberOfCuts];
	97	TString *mesonCutArray = new TString[numberOfCuts];
	98
	99	if (trainConfig == 1){
ae4f2cfb	100	cutarray[ 0] = "6010001042092970023220000000"; mesonCutArray[ 0] = "01522045009000";
39e8c2ce	101	} else if (trainConfig == 2) {
ae4f2cfb	102	cutarray[ 0] = "6120001042092970023220000000"; mesonCutArray[ 0] = "01522045009000";
39e8c2ce	103	} else if (trainConfig == 3) {
ae4f2cfb	104	cutarray[ 0] = "5010001042092970023220000000"; mesonCutArray[ 0] = "01522045009000";
39e8c2ce	105	} else if (trainConfig == 4) {
ae4f2cfb	106	cutarray[ 0] = "5020001042092970023220000000"; mesonCutArray[ 0] = "01522045009000";
39e8c2ce	107	} else if (trainConfig == 5) {
ae4f2cfb	108	cutarray[ 0] = "5120001042092970023220000000"; mesonCutArray[ 0] = "01522045009000";
39e8c2ce	109	} else if (trainConfig == 6) {
ae4f2cfb	110	cutarray[ 0] = "5240001042092970023220000000"; mesonCutArray[ 0] = "01522045009000";
39e8c2ce	111	} else if (trainConfig == 7) {
ae4f2cfb	112	cutarray[ 0] = "5460001042092970023220000000"; mesonCutArray[ 0] = "01522065009000";
39e8c2ce	113	} else if (trainConfig == 8) {
ae4f2cfb	114	cutarray[ 0] = "5480001042092970023220000000"; mesonCutArray[ 0] = "01522065009000";
39e8c2ce	115	} else if (trainConfig == 9) {
ae4f2cfb	116	cutarray[ 0] = "5450001042092970023220000000"; mesonCutArray[ 0] = "01522065009000";
39e8c2ce	117	} else if (trainConfig == 10) {
ae4f2cfb	118	cutarray[ 0] = "5560001042092970023220000000"; mesonCutArray[ 0] = "01522065009000";
39e8c2ce	119	} else if (trainConfig == 11) {
ae4f2cfb	120	cutarray[ 0] = "5680001042092970023220000000"; mesonCutArray[ 0] = "01522065009000";
39e8c2ce	121	} else if (trainConfig == 12) {
ae4f2cfb	122	cutarray[ 0] = "5670001042092970023220000000"; mesonCutArray[ 0] = "01522065009000";
39e8c2ce	123	} else if (trainConfig == 13) {
ae4f2cfb	124	cutarray[ 0] = "5780001042092970023220000000"; mesonCutArray[ 0] = "01522065009000";
39e8c2ce	125	} else if (trainConfig == 14) {
ae4f2cfb	126	cutarray[ 0] = "4690001042092970023220000000"; mesonCutArray[ 0] = "01522065009000";
39e8c2ce	127	} else if (trainConfig == 15) {
ae4f2cfb	128	cutarray[ 0] = "5890001042092970023220000000"; mesonCutArray[ 0] = "01522065009000";
c9663716	129	} else if (trainConfig == 16){
7e6c9796	130	cutarray[ 0] = "6010001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
c9663716	131	} else if (trainConfig == 17) {
7e6c9796	132	cutarray[ 0] = "6120001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
c9663716	133	} else if (trainConfig == 18) {
7e6c9796	134	cutarray[ 0] = "5010001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
c9663716	135	} else if (trainConfig == 19) {
7e6c9796	136	cutarray[ 0] = "5020001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
c9663716	137	} else if (trainConfig == 20) {
7e6c9796	138	cutarray[ 0] = "5120001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
c9663716	139	} else if (trainConfig == 21) {
7e6c9796	140	cutarray[ 0] = "5240001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
c9663716	141	} else if (trainConfig == 22) {
7e6c9796	142	cutarray[ 0] = "5460001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
c9663716	143	} else if (trainConfig == 23) {
7e6c9796	144	cutarray[ 0] = "5480001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
c9663716	145	} else if (trainConfig == 24) {
7e6c9796	146	cutarray[ 0] = "5450001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
c9663716	147	} else if (trainConfig == 25) {
7e6c9796	148	cutarray[ 0] = "5560001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
c9663716	149	} else if (trainConfig == 26) {
7e6c9796	150	cutarray[ 0] = "5680001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
c9663716	151	} else if (trainConfig == 27) {
7e6c9796	152	cutarray[ 0] = "5670001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
c9663716	153	} else if (trainConfig == 28) {
7e6c9796	154	cutarray[ 0] = "5780001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
c9663716	155	} else if (trainConfig == 29) {
7e6c9796	156	cutarray[ 0] = "4690001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
c9663716	157	} else if (trainConfig == 30) {
7e6c9796	158	cutarray[ 0] = "5890001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
955b2701	159	} else if (trainConfig == 31) {
7e6c9796	160	cutarray[ 0] = "5080001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
	161	} else if (trainConfig == 32) {
	162	cutarray[ 0] = "5250001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
	163	} else if (trainConfig == 33) {
	164	cutarray[ 0] = "5350001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
	165	} else if (trainConfig == 34) {
	166	cutarray[ 0] = "5450001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
	167	} else if (trainConfig == 35) {
	168	cutarray[ 0] = "5340001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
	169	} else if (trainConfig == 36) {
	170	cutarray[ 0] = "5230001002092970028250400000"; mesonCutArray[ 0] = "01525065000000";
	171	} else {
39e8c2ce	172	Error(Form("GammaConvV1_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
	173	return;
	174	}
	175
	176	TList *ConvCutList = new TList();
	177	TList *MesonCutList = new TList();
	178
	179	TList *HeaderList = new TList();
75f195f2	180	TObjString *Header1 = new TObjString("BOX");
39e8c2ce	181	HeaderList->Add(Header1);
	182	// TObjString *Header3 = new TObjString("eta_2");
	183	// HeaderList->Add(Header3);
	184
	185	ConvCutList->SetOwner(kTRUE);
	186	AliConversionCuts *analysisCuts = new AliConversionCuts[numberOfCuts];
	187	MesonCutList->SetOwner(kTRUE);
	188	AliConversionMesonCuts *analysisMesonCuts = new AliConversionMesonCuts[numberOfCuts];
	189
	190	for(Int_t i = 0; i<numberOfCuts; i++){
	191	analysisCuts[i] = new AliConversionCuts();
	192	analysisCuts[i]->InitializeCutsFromCutString(cutarray[i].Data());
	193	ConvCutList->Add(analysisCuts[i]);
	194
	195	analysisCuts[i]->SetFillCutHistograms("",kFALSE);
	196	analysisMesonCuts[i] = new AliConversionMesonCuts();
	197	analysisMesonCuts[i]->InitializeCutsFromCutString(mesonCutArray[i].Data());
	198	MesonCutList->Add(analysisMesonCuts[i]);
	199	analysisMesonCuts[i]->SetFillCutHistograms("");
	200	analysisCuts[i]->SetAcceptedHeader(HeaderList);
	201	}
	202
	203	task->SetConversionCutList(numberOfCuts,ConvCutList);
	204	task->SetMesonCutList(numberOfCuts,MesonCutList);
	205	task->SetMoveParticleAccordingToVertex(kTRUE);
	206	task->SetDoMesonAnalysis(kTRUE);
ae4f2cfb	207	task->SetDoMesonQA(enableQAMesonTask); //Attention new switch for Pi0 QA
ae4f2cfb	208	task->SetDoPhotonQA(enableQAPhotonTask); //Attention new switch small for Photon QA
39e8c2ce	209
	210	//connect containers
	211	AliAnalysisDataContainer *coutput =
	212	mgr->CreateContainer(Form("GammaConvV1_%i",trainConfig), TList::Class(),
	213	AliAnalysisManager::kOutputContainer,Form("GammaConvV1_%i.root",trainConfig));
	214
	215	mgr->AddTask(task);
	216	mgr->ConnectInput(task,0,cinput);
	217	mgr->ConnectOutput(task,1,coutput);
	218
	219	return;
	220
	221	}