update LMEE analyses
authorcbaumann <cbaumann@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 11 Feb 2013 11:55:10 +0000 (11:55 +0000)
committercbaumann <cbaumann@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 11 Feb 2013 11:55:10 +0000 (11:55 +0000)
PWGDQ/dielectron/macrosLMEE/AddTask_cbaumann_LMEEpp.C [new file with mode: 0644]
PWGDQ/dielectron/macrosLMEE/LMEECutLibAOD.C

diff --git a/PWGDQ/dielectron/macrosLMEE/AddTask_cbaumann_LMEEpp.C b/PWGDQ/dielectron/macrosLMEE/AddTask_cbaumann_LMEEpp.C
new file mode 100644 (file)
index 0000000..73f22f9
--- /dev/null
@@ -0,0 +1,96 @@
+AliAnalysisTask *AddTask_cbaumann_LMEEpp(){
+
+
+  //get the current analysis manager
+  AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
+  if (!mgr) {
+    Error("AddTask_cbaumann", "No analysis manager found.");
+    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: 
+  TString configBasePath("$TRAIN_ROOT/cbaumann_dielectron/");
+  TString trainRoot=gSystem->Getenv("TRAIN_ROOT");
+
+  //Base Directory for GRID / LEGO Train
+  if (trainRoot.IsNull()) configBasePath= "$ALICE_ROOT/PWGDQ/dielectron/macrosLMEE/";
+
+  TString configFile("Config_lowmasspp.C");
+  TString configFilePath(configBasePath+configFile);
+
+
+  if (!gROOT->GetListOfGlobalFunctions()->FindObject(configFile.Data()))
+         gROOT->LoadMacro(configFilePath.Data());
+
+
+  //create task and add it to the manager (MB)
+  AliAnalysisTaskMultiDielectron *taskMB = new AliAnalysisTaskMultiDielectron("MultiDieMB");
+  if (!hasMC) taskMB->UsePhysicsSelection();
+  //taskMB->SelectCollisionCandidates(AliVEvent::kMB);
+  taskMB->SelectCollisionCandidates(AliVEvent::kINT7);
+  taskMB->SetRejectPileup();
+
+          //Add event filter
+  AliDielectronEventCuts *eventCuts=new AliDielectronEventCuts("eventCuts","Vertex Track && |vtxZ|<10 && ncontrib>0");
+  eventCuts->SetRequireVertex();
+  eventCuts->SetVertexZ(-10.,10.);
+  eventCuts->SetMinVtxContributors(1);
+
+  taskMB->SetEventFilter(eventCuts);
+
+   mgr->AddTask(taskMB);
+
+
+  for (Int_t i=0; i<nDie; ++i){ //nDie defined in config file
+    //MB
+    AliDielectron *diel_lowMB = Config_lowmasspp(i);
+    if(!diel_lowMB)continue;
+    taskMB->AddDielectron(diel_lowMB);
+
+  }//loop
+  //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,
+                         "cbaumann_lowmass.root");
+
+  AliAnalysisDataContainer *cOutputHist2 =
+    mgr->CreateContainer("CF_diel_lowmass",
+                         TList::Class(),
+                         AliAnalysisManager::kOutputContainer,
+                         "cbaumann_lowmass.root");
+
+  AliAnalysisDataContainer *cOutputHist3 =
+    mgr->CreateContainer("cbauman_lowmass_EventStat",
+                         TList::Class(),
+                         AliAnalysisManager::kOutputContainer,
+                         "cbaumann_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);
+
+  return taskMB;
+
+}
index a258183..99a21a7 100644 (file)
@@ -202,14 +202,14 @@ class LMEECutLib {
          //___________________________________________
 
          AliDielectronVarCuts *pidTPCsignalWide = new AliDielectronVarCuts("pidTPCsignalWide","cut on the TPC signal");
-         pidTPCsignalWide->AddCut(AliDielectronVarManager::kTPCsignal,70.,90.);
+         pidTPCsignalWide->AddCut(AliDielectronVarManager::kTPCsignal,75.,90.);
 
 
          AliDielectronVarCuts *pTPC = new AliDielectronVarCuts("P>.4","P>.4");
-         pTPC->AddCut(AliDielectronVarManager::kPt,.4,5.);
+         pTPC->AddCut(AliDielectronVarManager::kPt,.4,2.5);
          
          AliDielectronVarCuts *pMin = new AliDielectronVarCuts("P>.2","P>.2");
-         pMin->AddCut(AliDielectronVarManager::kPt,.2,5.);
+         pMin->AddCut(AliDielectronVarManager::kPt,.2,2.5);
 
          switch (cutSet) {
                case kPbPb2011NoPID: