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

void AddTask_GammaConvV1_pp(  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
                              TString cutnumberAODBranch = "0000000060084001001500000" // cutnumber for AOD branch
                           ) {

   // ================= 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 = "0000000002084000002200000000"; 
   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 =====
   //            find input container
   AliAnalysisTaskGammaConvV1 *task=NULL;
   task= new AliAnalysisTaskGammaConvV1(Form("GammaConvV1_%i",trainConfig));
   task->SetIsHeavyIon(0);
   task->SetIsMC(isMC);
   // Cut Numbers to use in Analysis
   Int_t numberOfCuts = 4;

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

   if(trainConfig == 1){
      cutarray[ 0] = "0000012002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , only boxes
      cutarray[ 1] = "0001012002093663003800000000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD, V0AND , only boxes
      cutarray[ 2] = "0000012002093260003800000000"; mesonCutArray[2] = "01631031009000"; //standard cut Gamma pp 2-76TeV , only boxes
      cutarray[ 3] = "0000012002093260003800000000"; mesonCutArray[3] = "01631031009000"; //standard cut Gamma pp 2-76TeV , only boxes
   } else if (trainConfig == 2) {
      cutarray[ 0] = "0000011002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , only Minbias MC
      cutarray[ 1] = "0001011002093663003800000000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD, V0AND
      cutarray[ 2] = "0000011002093260003800000000"; mesonCutArray[2] = "01631031009000"; //standard cut Gamma pp 2-76TeV
      cutarray[ 3] = "0000011002093260003800000000"; mesonCutArray[3] = "01631031009000"; //standard cut Gamma pp 2-76TeV
   } else if (trainConfig == 3) {
      cutarray[ 0] = "0002011002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , only Minbias MC
      cutarray[ 1] = "0003011002093663003800000000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD, V0AND , only Minbias MC
      cutarray[ 2] = "0002012002093663003800000000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , only Boxes MC
      cutarray[ 3] = "0003012002093663003800000000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD, V0AND, only Boxes MC
   } else if (trainConfig == 4) {
      cutarray[ 0] = "0000011002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , all photon qualities
      cutarray[ 1] = "0000011002093663003800020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 1
      cutarray[ 2] = "0000011002093663003800030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 2
      cutarray[ 3] = "0000011002093663003800040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 3
   } else if (trainConfig == 5) {
      cutarray[ 0] = "0000011007093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , all photon qualities, min R = 35 cm
      cutarray[ 1] = "0000011007093663003800020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 1, min R = 35 cm
      cutarray[ 2] = "0000011007093663003800030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 2, min R = 35 cm
	  cutarray[ 3] = "0000011007093663003800040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 3, min R = 35 cm
   } else if (trainConfig == 6) {
      cutarray[ 0] = "0000011002083663003200000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, all photon qualities
      cutarray[ 1] = "0000011002083663003200020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 1
      cutarray[ 2] = "0000011002083663003200030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 2
	  cutarray[ 3] = "0000011002083663003200040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 3   
   } else if (trainConfig == 7) {
      cutarray[ 0] = "0000011007083663003200000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, all photon qualities, min R = 35 cm
      cutarray[ 1] = "0000011007083663003200020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 1, min R = 35 cm
      cutarray[ 2] = "0000011007083663003200030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 2, min R = 35 cm
	  cutarray[ 3] = "0000011007083663003200040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 3, min R = 35 cm
   } else if (trainConfig == 8) {
      cutarray[ 0] = "0000011002083663000200000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 7TeV, all photon qualities
      cutarray[ 1] = "0000011002083663000200020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 1
      cutarray[ 2] = "0000011002083663000200030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 2
	  cutarray[ 3] = "0000011002083663000200040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 3   
   } else if (trainConfig == 9) {
      cutarray[ 0] = "0000011007083663000200000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 7TeV, all photon qualities, min R = 35 cm
      cutarray[ 1] = "0000011007083663000200020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 1, min R = 35 cm
      cutarray[ 2] = "0000011007083663000200030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 2, min R = 35 cm
	  cutarray[ 3] = "0000011007083663000200040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 3, min R = 35 cm	   
   } else if (trainConfig == 10) {
      cutarray[ 0] = "0002011002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , all photon qualities
      cutarray[ 1] = "0002011002093663003800020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 1
      cutarray[ 2] = "0002011002093663003800030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 2
      cutarray[ 3] = "0002011002093663003800040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 3
   } else if (trainConfig == 11) {
      cutarray[ 0] = "0002011007093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , all photon qualities, min R = 35 cm
      cutarray[ 1] = "0002011007093663003800020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 1, min R = 35 cm
      cutarray[ 2] = "0002011007093663003800030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 2, min R = 35 cm
      cutarray[ 3] = "0002011007093663003800040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 3, min R = 35 cm
   } else if (trainConfig == 12) {
      cutarray[ 0] = "00000010020929700282504000000"; mesonCutArray[0] = "01525065000000"; //standard cut LHC11h pp 2.76TeV 
      cutarray[ 1] = "00000010320929700282504000000"; mesonCutArray[0] = "01525065000000"; //variation eta 0.65
      cutarray[ 2] = "00000010420929700282504000000"; mesonCutArray[2] = "01525065000000"; //variation eta 0.75
      cutarray[ 3] = "00000010020929500282504000000"; mesonCutArray[3] = "01525065000000"; //variation pion p dEdx 0.3-5.
   } else if (trainConfig == 13) {
      cutarray[ 0] = "00000010024929700282504000000"; mesonCutArray[0] = "01525065000000"; //variation pt 0.075 
      cutarray[ 1] = "00000010021929700282504000000"; mesonCutArray[0] = "01525065000000"; //variation pt 0.1
      cutarray[ 2] = "00000010020629700282504000000"; mesonCutArray[2] = "01525065000000"; //variation TPC cls 0.7
      cutarray[ 3] = "00000010020829700282504000000"; mesonCutArray[3] = "01525065000000"; //variation TPC cls 0.35 
   } else if (trainConfig == 14) {
      cutarray[ 0] = "00000010020939700282504000000"; mesonCutArray[0] = "01525065000000"; //variation edEdx -4,5
      cutarray[ 1] = "00000010020969700282504000000"; mesonCutArray[0] = "01525065000000"; //variation edEdx -2.5,4
      cutarray[ 2] = "00000010020929700382504000000"; mesonCutArray[2] = "01525065000000"; //variation TOF el. PID -3,5
      cutarray[ 3] = "00000010020929700482504000000"; mesonCutArray[3] = "01525065000000"; //variation TOF el. PID -2,3
   } else if (trainConfig == 15) {
      cutarray[ 0] = "00000010020929700292504000000"; mesonCutArray[0] = "01525065000000"; //variation qt 0.03
      cutarray[ 1] = "00000010020929700222504000000"; mesonCutArray[0] = "01525065000000"; //variation qt 0.07 no2D
      cutarray[ 2] = "00000010020929700281504000000"; mesonCutArray[2] = "01525065000000"; //variation chi2 50.
      cutarray[ 3] = "00000010020929700288504000000"; mesonCutArray[3] = "01525065000000"; //variation chi2 20.
   } else if (trainConfig == 16) {
      cutarray[ 0] = "00000010020929700282604000000"; mesonCutArray[0] = "01525065000000"; //variation psi pair 0.05
      cutarray[ 1] = "00000010020929700282804000000"; mesonCutArray[0] = "01525065000000"; //variation psi pair 0.2
      cutarray[ 2] = "00000010020929700282500000000"; mesonCutArray[2] = "01525065000000"; //variation cosPA -1
      cutarray[ 3] = "00000010020929700282504000000"; mesonCutArray[3] = "01525055000000"; //variation alpha 0.75
   } 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 *Header2 = new TObjString("BOX");
   HeaderList->Add(Header2);

   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
33bbd991	1	void AddTask_GammaConvV1_pp( Int_t trainConfig = 1, //change different set of cuts
33bbd991	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
a7f0b705	5	TString fileNameInputForWeighting = "MCSpectraInput.root", // path to file for weigting input
a7f0b705	6	TString cutnumberAODBranch = "0000000060084001001500000" // cutnumber for AOD branch
33bbd991	7	) {
	8
	9	// ================= Load Librariers =================================
	10	gSystem->Load("libCore.so");
	11	gSystem->Load("libTree.so");
	12	gSystem->Load("libGeom.so");
	13	gSystem->Load("libVMC.so");
	14	gSystem->Load("libPhysics.so");
	15	gSystem->Load("libMinuit");
	16	gSystem->Load("libSTEERBase");
	17	gSystem->Load("libESD");
	18	gSystem->Load("libAOD");
	19	gSystem->Load("libANALYSIS");
	20	gSystem->Load("libANALYSISalice");
	21	gSystem->Load("libPWGGAGammaConv.so");
	22	gSystem->Load("libCDB.so");
	23	gSystem->Load("libSTEER.so");
	24	gSystem->Load("libSTEERBase.so");
	25	gSystem->Load("libTENDER.so");
	26	gSystem->Load("libTENDERSupplies.so");
	27
	28	// ================== GetAnalysisManager ===============================
	29	AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
	30	if (!mgr) {
	31	Error(Form("AddTask_GammaConvV1_%i",trainConfig), "No analysis manager found.");
	32	return ;
	33	}
	34
	35	// ================== GetInputEventHandler =============================
	36	AliVEventHandler *inputHandler=mgr->GetInputEventHandler();
	37
	38	//========= Add PID Reponse to ANALYSIS manager ====
	39	if(!(AliPIDResponse*)mgr->GetTask("PIDResponseTask")){
	40	gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
	41	AddTaskPIDResponse(isMC);
	42	}
	43
	44	//========= Set Cutnumber for V0Reader ================================
ae4f2cfb	45	TString cutnumber = "0000000002084000002200000000";
b0e0bc89	46	AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
33bbd991	47
	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	}
a7f0b705	71	if(inputHandler->IsA()==AliAODInputHandler::Class()){
	72	// AOD mode
	73	fV0ReaderV1->SetDeltaAODBranchName(Form("GammaConv_%s_gamma",cutnumberAODBranch.Data()));
	74	}
33bbd991	75	fV0ReaderV1->Init();
	76
	77	AliLog::SetGlobalLogLevel(AliLog::kInfo);
	78
33bbd991	79	//connect input V0Reader
	80	mgr->AddTask(fV0ReaderV1);
	81	mgr->ConnectInput(fV0ReaderV1,0,cinput);
	82
	83	}
	84
	85	//================================================
	86	//========= Add task to the ANALYSIS manager =====
33bbd991	87	// find input container
	88	AliAnalysisTaskGammaConvV1 *task=NULL;
	89	task= new AliAnalysisTaskGammaConvV1(Form("GammaConvV1_%i",trainConfig));
	90	task->SetIsHeavyIon(0);
	91	task->SetIsMC(isMC);
	92	// Cut Numbers to use in Analysis
	93	Int_t numberOfCuts = 4;
	94
	95	TString *cutarray = new TString[numberOfCuts];
	96	TString *mesonCutArray = new TString[numberOfCuts];
	97
b0e0bc89	98	if(trainConfig == 1){
ae4f2cfb	99	cutarray[ 0] = "0000012002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , only boxes
	100	cutarray[ 1] = "0001012002093663003800000000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD, V0AND , only boxes
	101	cutarray[ 2] = "0000012002093260003800000000"; mesonCutArray[2] = "01631031009000"; //standard cut Gamma pp 2-76TeV , only boxes
	102	cutarray[ 3] = "0000012002093260003800000000"; mesonCutArray[3] = "01631031009000"; //standard cut Gamma pp 2-76TeV , only boxes
33bbd991	103	} else if (trainConfig == 2) {
ae4f2cfb	104	cutarray[ 0] = "0000011002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , only Minbias MC
	105	cutarray[ 1] = "0001011002093663003800000000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD, V0AND
	106	cutarray[ 2] = "0000011002093260003800000000"; mesonCutArray[2] = "01631031009000"; //standard cut Gamma pp 2-76TeV
072021d3	107	cutarray[ 3] = "0000011002093260003800000000"; mesonCutArray[3] = "01631031009000"; //standard cut Gamma pp 2-76TeV
33bbd991	108	} else if (trainConfig == 3) {
ae4f2cfb	109	cutarray[ 0] = "0002011002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , only Minbias MC
	110	cutarray[ 1] = "0003011002093663003800000000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD, V0AND , only Minbias MC
	111	cutarray[ 2] = "0002012002093663003800000000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , only Boxes MC
	112	cutarray[ 3] = "0003012002093663003800000000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD, V0AND, only Boxes MC
072021d3	113	} else if (trainConfig == 4) {
	114	cutarray[ 0] = "0000011002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , all photon qualities
	115	cutarray[ 1] = "0000011002093663003800020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 1
	116	cutarray[ 2] = "0000011002093663003800030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 2
88202ecb	117	cutarray[ 3] = "0000011002093663003800040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 3
072021d3	118	} else if (trainConfig == 5) {
	119	cutarray[ 0] = "0000011007093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , all photon qualities, min R = 35 cm
	120	cutarray[ 1] = "0000011007093663003800020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 1, min R = 35 cm
	121	cutarray[ 2] = "0000011007093663003800030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 2, min R = 35 cm
	122	cutarray[ 3] = "0000011007093663003800040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , photon quality 3, min R = 35 cm
e260262f	123	} else if (trainConfig == 6) {
072021d3	124	cutarray[ 0] = "0000011002083663003200000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, all photon qualities
	125	cutarray[ 1] = "0000011002083663003200020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 1
	126	cutarray[ 2] = "0000011002083663003200030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 2
	127	cutarray[ 3] = "0000011002083663003200040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 3
e260262f	128	} else if (trainConfig == 7) {
072021d3	129	cutarray[ 0] = "0000011007083663003200000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, all photon qualities, min R = 35 cm
	130	cutarray[ 1] = "0000011007083663003200020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 1, min R = 35 cm
	131	cutarray[ 2] = "0000011007083663003200030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 2, min R = 35 cm
	132	cutarray[ 3] = "0000011007083663003200040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 7TeV, with qt cut 0.05, photon quality 3, min R = 35 cm
e260262f	133	} else if (trainConfig == 8) {
072021d3	134	cutarray[ 0] = "0000011002083663000200000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 7TeV, all photon qualities
	135	cutarray[ 1] = "0000011002083663000200020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 1
	136	cutarray[ 2] = "0000011002083663000200030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 2
	137	cutarray[ 3] = "0000011002083663000200040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 3
e260262f	138	} else if (trainConfig == 9) {
072021d3	139	cutarray[ 0] = "0000011007083663000200000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 7TeV, all photon qualities, min R = 35 cm
	140	cutarray[ 1] = "0000011007083663000200020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 1, min R = 35 cm
	141	cutarray[ 2] = "0000011007083663000200030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 2, min R = 35 cm
	142	cutarray[ 3] = "0000011007083663000200040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 7TeV, photon quality 3, min R = 35 cm
048db199	143	} else if (trainConfig == 10) {
	144	cutarray[ 0] = "0002011002093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , all photon qualities
	145	cutarray[ 1] = "0002011002093663003800020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 1
	146	cutarray[ 2] = "0002011002093663003800030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 2
88202ecb	147	cutarray[ 3] = "0002011002093663003800040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 3
048db199	148	} else if (trainConfig == 11) {
	149	cutarray[ 0] = "0002011007093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , all photon qualities, min R = 35 cm
	150	cutarray[ 1] = "0002011007093663003800020000"; mesonCutArray[1] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 1, min R = 35 cm
	151	cutarray[ 2] = "0002011007093663003800030000"; mesonCutArray[2] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 2, min R = 35 cm
88202ecb	152	cutarray[ 3] = "0002011007093663003800040000"; mesonCutArray[3] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , photon quality 3, min R = 35 cm
	153	} else if (trainConfig == 12) {
	154	cutarray[ 0] = "00000010020929700282504000000"; mesonCutArray[0] = "01525065000000"; //standard cut LHC11h pp 2.76TeV
	155	cutarray[ 1] = "00000010320929700282504000000"; mesonCutArray[0] = "01525065000000"; //variation eta 0.65
	156	cutarray[ 2] = "00000010420929700282504000000"; mesonCutArray[2] = "01525065000000"; //variation eta 0.75
	157	cutarray[ 3] = "00000010020929500282504000000"; mesonCutArray[3] = "01525065000000"; //variation pion p dEdx 0.3-5.
	158	} else if (trainConfig == 13) {
	159	cutarray[ 0] = "00000010024929700282504000000"; mesonCutArray[0] = "01525065000000"; //variation pt 0.075
	160	cutarray[ 1] = "00000010021929700282504000000"; mesonCutArray[0] = "01525065000000"; //variation pt 0.1
	161	cutarray[ 2] = "00000010020629700282504000000"; mesonCutArray[2] = "01525065000000"; //variation TPC cls 0.7
	162	cutarray[ 3] = "00000010020829700282504000000"; mesonCutArray[3] = "01525065000000"; //variation TPC cls 0.35
	163	} else if (trainConfig == 14) {
	164	cutarray[ 0] = "00000010020939700282504000000"; mesonCutArray[0] = "01525065000000"; //variation edEdx -4,5
	165	cutarray[ 1] = "00000010020969700282504000000"; mesonCutArray[0] = "01525065000000"; //variation edEdx -2.5,4
	166	cutarray[ 2] = "00000010020929700382504000000"; mesonCutArray[2] = "01525065000000"; //variation TOF el. PID -3,5
	167	cutarray[ 3] = "00000010020929700482504000000"; mesonCutArray[3] = "01525065000000"; //variation TOF el. PID -2,3
	168	} else if (trainConfig == 15) {
	169	cutarray[ 0] = "00000010020929700292504000000"; mesonCutArray[0] = "01525065000000"; //variation qt 0.03
	170	cutarray[ 1] = "00000010020929700222504000000"; mesonCutArray[0] = "01525065000000"; //variation qt 0.07 no2D
	171	cutarray[ 2] = "00000010020929700281504000000"; mesonCutArray[2] = "01525065000000"; //variation chi2 50.
	172	cutarray[ 3] = "00000010020929700288504000000"; mesonCutArray[3] = "01525065000000"; //variation chi2 20.
c323a28f	173	} else if (trainConfig == 16) {
88202ecb	174	cutarray[ 0] = "00000010020929700282604000000"; mesonCutArray[0] = "01525065000000"; //variation psi pair 0.05
	175	cutarray[ 1] = "00000010020929700282804000000"; mesonCutArray[0] = "01525065000000"; //variation psi pair 0.2
	176	cutarray[ 2] = "00000010020929700282500000000"; mesonCutArray[2] = "01525065000000"; //variation cosPA -1
	177	cutarray[ 3] = "00000010020929700282504000000"; mesonCutArray[3] = "01525055000000"; //variation alpha 0.75
	178	} else {
072021d3	179	Error(Form("GammaConvV1_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
072021d3	180	return;
33bbd991	181	}
	182
	183	TList *ConvCutList = new TList();
	184	TList *MesonCutList = new TList();
	185
	186	TList *HeaderList = new TList();
	187	TObjString *Header2 = new TObjString("BOX");
	188	HeaderList->Add(Header2);
	189
	190	ConvCutList->SetOwner(kTRUE);
	191	AliConversionCuts *analysisCuts = new AliConversionCuts[numberOfCuts];
	192	MesonCutList->SetOwner(kTRUE);
	193	AliConversionMesonCuts *analysisMesonCuts = new AliConversionMesonCuts[numberOfCuts];
	194
	195
	196	for(Int_t i = 0; i<numberOfCuts; i++){
	197	analysisCuts[i] = new AliConversionCuts();
	198	analysisCuts[i]->InitializeCutsFromCutString(cutarray[i].Data());
	199	ConvCutList->Add(analysisCuts[i]);
	200
	201	analysisCuts[i]->SetFillCutHistograms("",kFALSE);
	202	analysisMesonCuts[i] = new AliConversionMesonCuts();
	203	analysisMesonCuts[i]->InitializeCutsFromCutString(mesonCutArray[i].Data());
	204	MesonCutList->Add(analysisMesonCuts[i]);
	205	analysisMesonCuts[i]->SetFillCutHistograms("");
	206	analysisCuts[i]->SetAcceptedHeader(HeaderList);
	207	}
	208
	209	task->SetConversionCutList(numberOfCuts,ConvCutList);
	210	task->SetMesonCutList(numberOfCuts,MesonCutList);
	211	task->SetMoveParticleAccordingToVertex(kTRUE);
	212	task->SetDoMesonAnalysis(kTRUE);
ae4f2cfb	213	task->SetDoMesonQA(enableQAMesonTask); //Attention new switch for Pi0 QA
ae4f2cfb	214	task->SetDoPhotonQA(enableQAPhotonTask); //Attention new switch small for Photon QA
33bbd991	215
	216	//connect containers
	217	AliAnalysisDataContainer *coutput =
b0e0bc89	218	mgr->CreateContainer(Form("GammaConvV1_%i",trainConfig), TList::Class(),
33bbd991	219	AliAnalysisManager::kOutputContainer,Form("GammaConvV1_%i.root",trainConfig));
	220
	221	mgr->AddTask(task);
	222	mgr->ConnectInput(task,0,cinput);
	223	mgr->ConnectOutput(task,1,coutput);
	224
	225	return;
	226
	227	}