ROOT logo
void AddTask_GammaConvV1_pp2(  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" 
                           ) {

	// ================= 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("libCDB.so");
	gSystem->Load("libSTEER.so");
	gSystem->Load("libSTEERBase.so");
	gSystem->Load("libTENDER.so");
	gSystem->Load("libTENDERSupplies.so");
	gSystem->Load("libPWGflowBase.so");
	gSystem->Load("libPWGflowTasks.so");
	gSystem->Load("libPWGGAGammaConv.so");

	Int_t isHeavyIon = 0;	
	   
	// ================== 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 cutnumberPhoton = "002000084000002200000000";
		TString cutnumberEvent = "0000000"; 
		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;
		}

		AliConvEventCuts *fEventCuts=NULL;
		if(cutnumberEvent!=""){
			fEventCuts= new AliConvEventCuts(cutnumberEvent.Data(),cutnumberEvent.Data());
			fEventCuts->SetPreSelectionCutFlag(kTRUE);
			if(fEventCuts->InitializeCutsFromCutString(cutnumberEvent.Data())){
				fV0ReaderV1->SetEventCuts(fEventCuts);
				fEventCuts->SetFillCutHistograms("",kTRUE);
			}
		}

		// Set AnalysisCut Number
		AliConversionPhotonCuts *fCuts=NULL;
		if(cutnumberPhoton!=""){
			fCuts= new AliConversionPhotonCuts(cutnumberPhoton.Data(),cutnumberPhoton.Data());
			fCuts->SetPreSelectionCutFlag(kTRUE);
			fCuts->SetIsHeavyIon(isHeavyIon);
			if(fCuts->InitializeCutsFromCutString(cutnumberPhoton.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(isHeavyIon);
	task->SetIsMC(isMC);
	// Cut Numbers to use in Analysis
	Int_t numberOfCuts = 1;

	TString *eventCutArray = new TString[numberOfCuts];
	TString *photonCutArray = new TString[numberOfCuts];
	TString *mesonCutArray = new TString[numberOfCuts];

	if (trainConfig == 1){ 
		eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002000093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , only boxes
	} else if (trainConfig == 2) { 
		eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002000093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD, V0AND , only boxes
	} else if (trainConfig == 3) { 
		eventCutArray[ 0] = "0000012"; photonCutArray[ 0] = "002000093260003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Gamma pp 2-76TeV , only boxes
	} else if (trainConfig == 4) { 
		eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002000093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , only Minbias MC
	} else if (trainConfig == 5) { 
		eventCutArray[ 0] = "0001011"; photonCutArray[ 0] = "002000093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD, V0AND
	} else if (trainConfig == 6) { 
		eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002000093260003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Gamma pp 2-76TeV
	} else if (trainConfig == 7) {    
		eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002000093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , only Minbias MC
	} else if (trainConfig == 8) {    
		eventCutArray[ 0] = "0001311"; photonCutArray[ 0] = "002000093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD, V0AND , only Minbias MC
	} else if (trainConfig == 9) {    
		eventCutArray[ 0] = "0000312"; photonCutArray[ 0] = "002000093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , only Boxes MC
	} else if (trainConfig == 10) { 
		eventCutArray[ 0] = "0001312"; photonCutArray[ 0] = "002000093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD, V0AND, only Boxes MC
	} else if (trainConfig == 11) { 
		eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002000093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , all photon qualities
	} else if (trainConfig == 12) { 
		eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "007000093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV without SDD , all photon qualities, min R = 35 cm
	} else if (trainConfig == 13) { 
		eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002000093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , all photon qualities
	} else if (trainConfig == 14) { 
		eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "007000093663003800000000"; mesonCutArray[0] = "01631031009000"; //standard cut Pi0 pp 2.76TeV with SDD , all photon qualities, min R = 35 cm
	} else if (trainConfig == 15) { 
		eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002000092970028250400000"; mesonCutArray[0] = "01525065000000"; //standard cut LHC11h pp 2.76TeV 
	} else if (trainConfig == 16) { 
		eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002000092570028250400000"; mesonCutArray[0] = "01521065000000"; //standard cut pp 8 TeV		
	} else {
		Error(Form("GammaConvV1_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
		return;
	}

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

	TList *HeaderList = new TList();
	TObjString *Header2 = new TObjString("BOX");
	HeaderList->Add(Header2);

	EventCutList->SetOwner(kTRUE);
	AliConvEventCuts **analysisEventCuts = new AliConvEventCuts*[numberOfCuts];
	ConvCutList->SetOwner(kTRUE);
	AliConversionPhotonCuts **analysisCuts = new AliConversionPhotonCuts*[numberOfCuts];
	MesonCutList->SetOwner(kTRUE);
	AliConversionMesonCuts **analysisMesonCuts = new AliConversionMesonCuts*[numberOfCuts];


	for(Int_t i = 0; i<numberOfCuts; i++){
		analysisEventCuts[i] = new AliConvEventCuts();
		analysisEventCuts[i]->InitializeCutsFromCutString(eventCutArray[i].Data());
		EventCutList->Add(analysisEventCuts[i]);
		analysisEventCuts[i]->SetFillCutHistograms("",kFALSE);
		
		analysisCuts[i] = new AliConversionPhotonCuts();
		analysisCuts[i]->InitializeCutsFromCutString(photonCutArray[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("");
		
		analysisEventCuts[i]->SetAcceptedHeader(HeaderList);
	}

	task->SetEventCutList(numberOfCuts,EventCutList);
	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;

}
 AddTask_GammaConvV1_pp2.C:1
 AddTask_GammaConvV1_pp2.C:2
 AddTask_GammaConvV1_pp2.C:3
 AddTask_GammaConvV1_pp2.C:4
 AddTask_GammaConvV1_pp2.C:5
 AddTask_GammaConvV1_pp2.C:6
 AddTask_GammaConvV1_pp2.C:7
 AddTask_GammaConvV1_pp2.C:8
 AddTask_GammaConvV1_pp2.C:9
 AddTask_GammaConvV1_pp2.C:10
 AddTask_GammaConvV1_pp2.C:11
 AddTask_GammaConvV1_pp2.C:12
 AddTask_GammaConvV1_pp2.C:13
 AddTask_GammaConvV1_pp2.C:14
 AddTask_GammaConvV1_pp2.C:15
 AddTask_GammaConvV1_pp2.C:16
 AddTask_GammaConvV1_pp2.C:17
 AddTask_GammaConvV1_pp2.C:18
 AddTask_GammaConvV1_pp2.C:19
 AddTask_GammaConvV1_pp2.C:20
 AddTask_GammaConvV1_pp2.C:21
 AddTask_GammaConvV1_pp2.C:22
 AddTask_GammaConvV1_pp2.C:23
 AddTask_GammaConvV1_pp2.C:24
 AddTask_GammaConvV1_pp2.C:25
 AddTask_GammaConvV1_pp2.C:26
 AddTask_GammaConvV1_pp2.C:27
 AddTask_GammaConvV1_pp2.C:28
 AddTask_GammaConvV1_pp2.C:29
 AddTask_GammaConvV1_pp2.C:30
 AddTask_GammaConvV1_pp2.C:31
 AddTask_GammaConvV1_pp2.C:32
 AddTask_GammaConvV1_pp2.C:33
 AddTask_GammaConvV1_pp2.C:34
 AddTask_GammaConvV1_pp2.C:35
 AddTask_GammaConvV1_pp2.C:36
 AddTask_GammaConvV1_pp2.C:37
 AddTask_GammaConvV1_pp2.C:38
 AddTask_GammaConvV1_pp2.C:39
 AddTask_GammaConvV1_pp2.C:40
 AddTask_GammaConvV1_pp2.C:41
 AddTask_GammaConvV1_pp2.C:42
 AddTask_GammaConvV1_pp2.C:43
 AddTask_GammaConvV1_pp2.C:44
 AddTask_GammaConvV1_pp2.C:45
 AddTask_GammaConvV1_pp2.C:46
 AddTask_GammaConvV1_pp2.C:47
 AddTask_GammaConvV1_pp2.C:48
 AddTask_GammaConvV1_pp2.C:49
 AddTask_GammaConvV1_pp2.C:50
 AddTask_GammaConvV1_pp2.C:51
 AddTask_GammaConvV1_pp2.C:52
 AddTask_GammaConvV1_pp2.C:53
 AddTask_GammaConvV1_pp2.C:54
 AddTask_GammaConvV1_pp2.C:55
 AddTask_GammaConvV1_pp2.C:56
 AddTask_GammaConvV1_pp2.C:57
 AddTask_GammaConvV1_pp2.C:58
 AddTask_GammaConvV1_pp2.C:59
 AddTask_GammaConvV1_pp2.C:60
 AddTask_GammaConvV1_pp2.C:61
 AddTask_GammaConvV1_pp2.C:62
 AddTask_GammaConvV1_pp2.C:63
 AddTask_GammaConvV1_pp2.C:64
 AddTask_GammaConvV1_pp2.C:65
 AddTask_GammaConvV1_pp2.C:66
 AddTask_GammaConvV1_pp2.C:67
 AddTask_GammaConvV1_pp2.C:68
 AddTask_GammaConvV1_pp2.C:69
 AddTask_GammaConvV1_pp2.C:70
 AddTask_GammaConvV1_pp2.C:71
 AddTask_GammaConvV1_pp2.C:72
 AddTask_GammaConvV1_pp2.C:73
 AddTask_GammaConvV1_pp2.C:74
 AddTask_GammaConvV1_pp2.C:75
 AddTask_GammaConvV1_pp2.C:76
 AddTask_GammaConvV1_pp2.C:77
 AddTask_GammaConvV1_pp2.C:78
 AddTask_GammaConvV1_pp2.C:79
 AddTask_GammaConvV1_pp2.C:80
 AddTask_GammaConvV1_pp2.C:81
 AddTask_GammaConvV1_pp2.C:82
 AddTask_GammaConvV1_pp2.C:83
 AddTask_GammaConvV1_pp2.C:84
 AddTask_GammaConvV1_pp2.C:85
 AddTask_GammaConvV1_pp2.C:86
 AddTask_GammaConvV1_pp2.C:87
 AddTask_GammaConvV1_pp2.C:88
 AddTask_GammaConvV1_pp2.C:89
 AddTask_GammaConvV1_pp2.C:90
 AddTask_GammaConvV1_pp2.C:91
 AddTask_GammaConvV1_pp2.C:92
 AddTask_GammaConvV1_pp2.C:93
 AddTask_GammaConvV1_pp2.C:94
 AddTask_GammaConvV1_pp2.C:95
 AddTask_GammaConvV1_pp2.C:96
 AddTask_GammaConvV1_pp2.C:97
 AddTask_GammaConvV1_pp2.C:98
 AddTask_GammaConvV1_pp2.C:99
 AddTask_GammaConvV1_pp2.C:100
 AddTask_GammaConvV1_pp2.C:101
 AddTask_GammaConvV1_pp2.C:102
 AddTask_GammaConvV1_pp2.C:103
 AddTask_GammaConvV1_pp2.C:104
 AddTask_GammaConvV1_pp2.C:105
 AddTask_GammaConvV1_pp2.C:106
 AddTask_GammaConvV1_pp2.C:107
 AddTask_GammaConvV1_pp2.C:108
 AddTask_GammaConvV1_pp2.C:109
 AddTask_GammaConvV1_pp2.C:110
 AddTask_GammaConvV1_pp2.C:111
 AddTask_GammaConvV1_pp2.C:112
 AddTask_GammaConvV1_pp2.C:113
 AddTask_GammaConvV1_pp2.C:114
 AddTask_GammaConvV1_pp2.C:115
 AddTask_GammaConvV1_pp2.C:116
 AddTask_GammaConvV1_pp2.C:117
 AddTask_GammaConvV1_pp2.C:118
 AddTask_GammaConvV1_pp2.C:119
 AddTask_GammaConvV1_pp2.C:120
 AddTask_GammaConvV1_pp2.C:121
 AddTask_GammaConvV1_pp2.C:122
 AddTask_GammaConvV1_pp2.C:123
 AddTask_GammaConvV1_pp2.C:124
 AddTask_GammaConvV1_pp2.C:125
 AddTask_GammaConvV1_pp2.C:126
 AddTask_GammaConvV1_pp2.C:127
 AddTask_GammaConvV1_pp2.C:128
 AddTask_GammaConvV1_pp2.C:129
 AddTask_GammaConvV1_pp2.C:130
 AddTask_GammaConvV1_pp2.C:131
 AddTask_GammaConvV1_pp2.C:132
 AddTask_GammaConvV1_pp2.C:133
 AddTask_GammaConvV1_pp2.C:134
 AddTask_GammaConvV1_pp2.C:135
 AddTask_GammaConvV1_pp2.C:136
 AddTask_GammaConvV1_pp2.C:137
 AddTask_GammaConvV1_pp2.C:138
 AddTask_GammaConvV1_pp2.C:139
 AddTask_GammaConvV1_pp2.C:140
 AddTask_GammaConvV1_pp2.C:141
 AddTask_GammaConvV1_pp2.C:142
 AddTask_GammaConvV1_pp2.C:143
 AddTask_GammaConvV1_pp2.C:144
 AddTask_GammaConvV1_pp2.C:145
 AddTask_GammaConvV1_pp2.C:146
 AddTask_GammaConvV1_pp2.C:147
 AddTask_GammaConvV1_pp2.C:148
 AddTask_GammaConvV1_pp2.C:149
 AddTask_GammaConvV1_pp2.C:150
 AddTask_GammaConvV1_pp2.C:151
 AddTask_GammaConvV1_pp2.C:152
 AddTask_GammaConvV1_pp2.C:153
 AddTask_GammaConvV1_pp2.C:154
 AddTask_GammaConvV1_pp2.C:155
 AddTask_GammaConvV1_pp2.C:156
 AddTask_GammaConvV1_pp2.C:157
 AddTask_GammaConvV1_pp2.C:158
 AddTask_GammaConvV1_pp2.C:159
 AddTask_GammaConvV1_pp2.C:160
 AddTask_GammaConvV1_pp2.C:161
 AddTask_GammaConvV1_pp2.C:162
 AddTask_GammaConvV1_pp2.C:163
 AddTask_GammaConvV1_pp2.C:164
 AddTask_GammaConvV1_pp2.C:165
 AddTask_GammaConvV1_pp2.C:166
 AddTask_GammaConvV1_pp2.C:167
 AddTask_GammaConvV1_pp2.C:168
 AddTask_GammaConvV1_pp2.C:169
 AddTask_GammaConvV1_pp2.C:170
 AddTask_GammaConvV1_pp2.C:171
 AddTask_GammaConvV1_pp2.C:172
 AddTask_GammaConvV1_pp2.C:173
 AddTask_GammaConvV1_pp2.C:174
 AddTask_GammaConvV1_pp2.C:175
 AddTask_GammaConvV1_pp2.C:176
 AddTask_GammaConvV1_pp2.C:177
 AddTask_GammaConvV1_pp2.C:178
 AddTask_GammaConvV1_pp2.C:179
 AddTask_GammaConvV1_pp2.C:180
 AddTask_GammaConvV1_pp2.C:181
 AddTask_GammaConvV1_pp2.C:182
 AddTask_GammaConvV1_pp2.C:183
 AddTask_GammaConvV1_pp2.C:184
 AddTask_GammaConvV1_pp2.C:185
 AddTask_GammaConvV1_pp2.C:186
 AddTask_GammaConvV1_pp2.C:187
 AddTask_GammaConvV1_pp2.C:188
 AddTask_GammaConvV1_pp2.C:189
 AddTask_GammaConvV1_pp2.C:190
 AddTask_GammaConvV1_pp2.C:191
 AddTask_GammaConvV1_pp2.C:192
 AddTask_GammaConvV1_pp2.C:193
 AddTask_GammaConvV1_pp2.C:194
 AddTask_GammaConvV1_pp2.C:195
 AddTask_GammaConvV1_pp2.C:196
 AddTask_GammaConvV1_pp2.C:197
 AddTask_GammaConvV1_pp2.C:198
 AddTask_GammaConvV1_pp2.C:199
 AddTask_GammaConvV1_pp2.C:200
 AddTask_GammaConvV1_pp2.C:201
 AddTask_GammaConvV1_pp2.C:202
 AddTask_GammaConvV1_pp2.C:203
 AddTask_GammaConvV1_pp2.C:204
 AddTask_GammaConvV1_pp2.C:205
 AddTask_GammaConvV1_pp2.C:206
 AddTask_GammaConvV1_pp2.C:207