From a5f8985aab52f4422f3a0a24138ff4ae76209b2c Mon Sep 17 00:00:00 2001 From: cbaumann Date: Mon, 4 Feb 2013 14:23:33 +0000 Subject: [PATCH] Updated AddTask for pPb --- .../macrosLMEE/AddTask_tbroeker_lowmasspPb.C | 161 +++--------------- 1 file changed, 22 insertions(+), 139 deletions(-) diff --git a/PWGDQ/dielectron/macrosLMEE/AddTask_tbroeker_lowmasspPb.C b/PWGDQ/dielectron/macrosLMEE/AddTask_tbroeker_lowmasspPb.C index 4c92df0445a..5e22983a5e9 100644 --- a/PWGDQ/dielectron/macrosLMEE/AddTask_tbroeker_lowmasspPb.C +++ b/PWGDQ/dielectron/macrosLMEE/AddTask_tbroeker_lowmasspPb.C @@ -8,15 +8,9 @@ AliAnalysisTask *AddTask_tbroeker_lowmasspPb(Bool_t getFromAlien=kFALSE){ return 0; } - - Bool_t RunEMCtrigger = 0; - Bool_t RunHighMulttrigger = 0; - Bool_t RunMBtrigger = 1; - //Do we have an MC handler? Bool_t hasMC=(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()!=0x0); - //Get the current train configuration //Directories for GSI train: @@ -49,176 +43,65 @@ AliAnalysisTask *AddTask_tbroeker_lowmasspPb(Bool_t getFromAlien=kFALSE){ gROOT->LoadMacro(configFilePath.Data()); - 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(); -} +//taskMB->SetRejectPileup(); - //Add event filter + //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); + taskMB->SetEventFilter(eventCuts); + mgr->AddTask(taskMB); //add dielectron analysis with different cuts to the task for (Int_t i=0; iSetNoPairing(kFALSE); taskMB->AddDielectron(diel_lowMB); - } - if(RunHighMulttrigger){ - //HighMult - AliDielectron *diel_lowHighMult = Config_lowmasspPb(i); - if(!diel_lowHighMult)continue; - diel_lowHighMult->SetNoPairing(kTRUE); - taskHighMult->AddDielectron(diel_lowHighMult); -} - if(RunEMCtrigger){ - //EMC1 - AliDielectron *diel_lowEMC1 = Config_lowmasspPb(i); - if(!diel_lowEMC1)continue; - diel_lowEMC1->SetNoPairing(kTRUE); - 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, - "tbroeker_lowmass.root"); - -// AliAnalysisDataContainer *cOutputHist2 = -// mgr->CreateContainer("CF_diel_lowmass", -// TList::Class(), -// AliAnalysisManager::kOutputContainer, -// "tbroeker_lowmass.root"); -// AliAnalysisDataContainer *cOutputHist2 = -// mgr->CreateContainer("lowmass_tbroeker_EventStat", -// TList::Class(), -// AliAnalysisManager::kOutputContainer, -// "tbroeker_lowmass.root"); - - 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){ - //create output container (HighMult) - AliAnalysisDataContainer *coutputHighMult1 = - mgr->CreateContainer("tree_lowmassHighMult", + //create output container + AliAnalysisDataContainer *coutput1 = + mgr->CreateContainer("tree_lowmass", TTree::Class(), AliAnalysisManager::kExchangeContainer, "default"); - AliAnalysisDataContainer *cOutputHistHighMult1 = - mgr->CreateContainer("Histos_diel_lowmassHighMult", + AliAnalysisDataContainer *cOutputHist1 = + mgr->CreateContainer("Histos_diel_lowmass", TList::Class(), AliAnalysisManager::kOutputContainer, - "lowmass_lowmassHighMult.root"); + "tbroeker_lowmass.root"); - AliAnalysisDataContainer *cOutputHistHighMult2 = - mgr->CreateContainer("CF_diel_lowmassHighMult", + AliAnalysisDataContainer *cOutputHist2 = + mgr->CreateContainer("CF_diel_lowmass", TList::Class(), AliAnalysisManager::kOutputContainer, - "lowmass_lowmassHighMult.root"); + "tbroeker_lowmass.root"); - AliAnalysisDataContainer *cOutputHistHighMult3 = - mgr->CreateContainer("lowmass_tbroeker_EventStatHighMult", + AliAnalysisDataContainer *cOutputHist3 = + mgr->CreateContainer("tbroeker_lowmass_EventStat", TList::Class(), AliAnalysisManager::kOutputContainer, - "lowmass_lowmassHighMult.root"); - } - - if(RunEMCtrigger){ + "tbroeker_lowmass.root"); - //create output container (EMC1) - AliAnalysisDataContainer *coutputEMC11 = - mgr->CreateContainer("tree_lowmassEMC1", - TTree::Class(), - AliAnalysisManager::kExchangeContainer, - "default"); + 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); - AliAnalysisDataContainer *cOutputHistEMC11 = - mgr->CreateContainer("Histos_diel_lowmassEMC1", - TList::Class(), - AliAnalysisManager::kOutputContainer, - "lowmass_lowmassEMC1.root"); - - AliAnalysisDataContainer *cOutputHistEMC12 = - mgr->CreateContainer("CF_diel_lowmassEMC1", - TList::Class(), - AliAnalysisManager::kOutputContainer, - "lowmass_lowmassEMC1.root"); - - AliAnalysisDataContainer *cOutputHistEMC13 = - mgr->CreateContainer("lowmass_tbroeker_EventStatEMC1", - TList::Class(), - AliAnalysisManager::kOutputContainer, - "lowmass_lowmassEMC1.root"); -} - -//questionable regarding variable visibility -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); -} -if(RunMBtrigger) return taskMB; -if(RunHighMulttrigger) return taskHighMult; -if(RunEMCtrigger) return taskEMC1; + return taskMB; } -- 2.43.0