ROOT logo


AliAnalysisTask *AddTask_mkoehler_lowmass(){


  //get the current analysis manager
  AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
  if (!mgr) {
    Error("AddTask_lowmass", "No analysis manager found.");
    return 0;
  }


  Bool_t RunEMCtrigger = 0;
  Bool_t RunHighMulttrigger = 0;
  Bool_t RunMBtrigger = 1;

  Bool_t hasMC=(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()!=0x0);

  //Get the current train configuration
  TString trainConfig=gSystem->Getenv("CONFIG_FILE");

  //set config file name
  TString configFile("$TRAIN_ROOT/mkoehler_lowmassee/Config_lowmass.C");
  TString trainRoot=gSystem->Getenv("TRAIN_ROOT");

  if (trainRoot.IsNull()) configFile="$ALICE_ROOT/PWGDQ/dielectron/macrosLMEE/Config_lowmass.C";

  if(RunMBtrigger){
  //create task and add it to the manager (MB)
  AliAnalysisTaskMultiDielectron *taskMB = new AliAnalysisTaskMultiDielectron("MultiDieMB");
  if (!hasMC) taskMB->UsePhysicsSelection();
  taskMB->SelectCollisionCandidates(AliVEvent::kMB);
  taskMB->SetRejectPileup();
	}

  if(RunHighMulttrigger){
  //create task and add it to the manager (HighMult)
  AliAnalysisTaskMultiDielectron *taskHighMult = new AliAnalysisTaskMultiDielectron("MultiDieHighMult");
  if (!hasMC) taskHighMult->UsePhysicsSelection();
  taskHighMult->SelectCollisionCandidates(AliVEvent::kHighMult);
  taskHighMult->SetRejectPileup();
	}

  if(RunEMCtrigger){
  //create task and add it to the manager (EMC1)
  AliAnalysisTaskMultiDielectron *taskEMC1 = new AliAnalysisTaskMultiDielectron("MultiDieEMC1");
  if (!hasMC) taskEMC1->UsePhysicsSelection();
  taskEMC1->SelectCollisionCandidates(AliVEvent::kEMC1);
  taskEMC1->SetRejectPileup();
}

          //Add event filter
  AliDielectronEventCuts *eventCuts=new AliDielectronEventCuts("eventCuts","Vertex Track && |vtxZ|<10 && ncontrib>0");
  eventCuts->SetRequireVertex();
  eventCuts->SetVertexZ(-10.,10.);
  eventCuts->SetMinVtxContributors(1);

   if(RunMBtrigger)taskMB->SetEventFilter(eventCuts);
   if(RunHighMulttrigger) taskHighMult->SetEventFilter(eventCuts);
   if(RunEMCtrigger) taskEMC1->SetEventFilter(eventCuts);

    if(RunMBtrigger) mgr->AddTask(taskMB);
    if(RunHighMulttrigger)mgr->AddTask(taskHighMult);
    if(RunEMCtrigger)mgr->AddTask(taskEMC1);

  //load dielectron configuration file
  gROOT->LoadMacro(configFile.Data());

  //add dielectron analysis with different cuts to the task
  for (Int_t i=0; i<nDie; ++i){ //nDie defined in config file
  if(RunMBtrigger){
    //MB
    AliDielectron *diel_lowMB = Config_lowmass(i);
    if(!diel_lowMB)continue;
    taskMB->AddDielectron(diel_lowMB);
	}
  if(RunHighMulttrigger){
    //HighMult
    AliDielectron *diel_lowHighMult = Config_lowmass(i);
    if(!diel_lowHighMult)continue;
    taskHighMult->AddDielectron(diel_lowHighMult);
}
	  if(RunEMCtrigger){
    //EMC1 
    AliDielectron *diel_lowEMC1 = Config_lowmass(i);
    if(!diel_lowEMC1)continue;
    taskEMC1->AddDielectron(diel_lowEMC1);
	}

  }//loop

  if(RunMBtrigger){
  //create output container
  AliAnalysisDataContainer *coutput1 =
    mgr->CreateContainer("tree_lowmass",
                         TTree::Class(),
                         AliAnalysisManager::kExchangeContainer,
                         "default");
  
  AliAnalysisDataContainer *cOutputHist1 =
    mgr->CreateContainer("Histos_diel_lowmass",
                         TList::Class(),
                         AliAnalysisManager::kOutputContainer,
                         "mkoehler_lowmass.root");

  AliAnalysisDataContainer *cOutputHist2 =
    mgr->CreateContainer("CF_diel_lowmass",
                         TList::Class(),
                         AliAnalysisManager::kOutputContainer,
                         "mkoehler_lowmass.root");

  AliAnalysisDataContainer *cOutputHist3 =
    mgr->CreateContainer("mkoehler_lowmass_EventStat",
                         TList::Class(),
                         AliAnalysisManager::kOutputContainer,
                         "mkoehler_lowmass.root");
	}
  if(RunHighMulttrigger){
  //create output container (HighMult)
  AliAnalysisDataContainer *coutputHighMult1 =
    mgr->CreateContainer("tree_lowmassHighMult",
                         TTree::Class(),
                         AliAnalysisManager::kExchangeContainer,
                         "default");
  
  AliAnalysisDataContainer *cOutputHistHighMult1 =
    mgr->CreateContainer("Histos_diel_lowmassHighMult",
                         TList::Class(),
                         AliAnalysisManager::kOutputContainer,
                         "mkoehler_lowmassHighMult.root");

  AliAnalysisDataContainer *cOutputHistHighMult2 =
    mgr->CreateContainer("CF_diel_lowmassHighMult",
                         TList::Class(),
                         AliAnalysisManager::kOutputContainer,
                         "mkoehler_lowmassHighMult.root");

  AliAnalysisDataContainer *cOutputHistHighMult3 =
    mgr->CreateContainer("mkoehler_lowmass_EventStatHighMult",
                         TList::Class(),
                         AliAnalysisManager::kOutputContainer,
                         "mkoehler_lowmassHighMult.root");
	}

	  if(RunEMCtrigger){

  //create output container (EMC1)
  AliAnalysisDataContainer *coutputEMC11 =
    mgr->CreateContainer("tree_lowmassEMC1",
                         TTree::Class(),
                         AliAnalysisManager::kExchangeContainer,
                         "default");
  
  AliAnalysisDataContainer *cOutputHistEMC11 =
    mgr->CreateContainer("Histos_diel_lowmassEMC1",
                         TList::Class(),
                         AliAnalysisManager::kOutputContainer,
                         "mkoehler_lowmassEMC1.root");

  AliAnalysisDataContainer *cOutputHistEMC12 =
    mgr->CreateContainer("CF_diel_lowmassEMC1",
                         TList::Class(),
                         AliAnalysisManager::kOutputContainer,
                         "mkoehler_lowmassEMC1.root");

  AliAnalysisDataContainer *cOutputHistEMC13 =
    mgr->CreateContainer("mkoehler_lowmass_EventStatEMC1",
                         TList::Class(),
                         AliAnalysisManager::kOutputContainer,
                         "mkoehler_lowmassEMC1.root");
}

  if(RunMBtrigger){
  mgr->ConnectInput(taskMB,  0, mgr->GetCommonInputContainer());
  mgr->ConnectOutput(taskMB, 0, coutput1 );
  mgr->ConnectOutput(taskMB, 1, cOutputHist1);
  mgr->ConnectOutput(taskMB, 2, cOutputHist2);
  mgr->ConnectOutput(taskMB, 3, cOutputHist3);
}
  if(RunHighMulttrigger){
  mgr->ConnectInput(taskHighMult,  0, mgr->GetCommonInputContainer());
  mgr->ConnectOutput(taskHighMult, 0, coutputHighMult1 );
  mgr->ConnectOutput(taskHighMult, 1, cOutputHistHighMult1);
  mgr->ConnectOutput(taskHighMult, 2, cOutputHistHighMult2);
  mgr->ConnectOutput(taskHighMult, 3, cOutputHistHighMult3);
}
  if(RunEMCtrigger){
  mgr->ConnectInput(taskEMC1,  0, mgr->GetCommonInputContainer());
  mgr->ConnectOutput(taskEMC1, 0, coutputEMC11 );
  mgr->ConnectOutput(taskEMC1, 1, cOutputHistEMC11);
  mgr->ConnectOutput(taskEMC1, 2, cOutputHistEMC12);
  mgr->ConnectOutput(taskEMC1, 3, cOutputHistEMC13);
}


  return taskMB;

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