.so cleanup: removed from gSystem->Load()

[u/mrichter/AliRoot.git] / JETAN / JetAnalysisManagerLoc.C

diff --git a/JETAN/JetAnalysisManagerLoc.C b/JETAN/JetAnalysisManagerLoc.C
index fc966a9838784bf9be621c847e65cd22aca495c9..87f4efa3d246e6ad4e1140b4f779ab588a9d0537 100644 (file)
--- a/JETAN/JetAnalysisManagerLoc.C
+++ b/JETAN/JetAnalysisManagerLoc.C
@@ -1,21 +1,23 @@
 void JetAnalysisManagerLoc()
 {
-      gSystem->Load("libTree.so");
-      gSystem->Load("libGeom.so");
-      gSystem->Load("libVMC.so");
-      gSystem->Load("libANALYSIS.so");
-      gSystem->Load("libSTEERBase.so");
-      gSystem->Load("libAOD.so");
-      gSystem->Load("libESD.so");
-      gSystem->Load("libJETAN.so");
-      gSystem->Load("libPWG0base.so");
+      gSystem->Load("libTree");
+      gSystem->Load("libPhysics");
+      gSystem->Load("libGeom");
+      gSystem->Load("libVMC");
+      gSystem->Load("libANALYSIS");
+      gSystem->Load("libSTEERBase");
+      gSystem->Load("libAOD");
+      gSystem->Load("libESD");
+      gSystem->Load("libANALYSISalice");
+      gSystem->Load("libJETAN");
+
      //
     if (gApplication) gApplication->InitializeGraphics();
     // Create the chain
     //
     gROOT->LoadMacro("CreateESDChain.C");
     TChain* chain = new TChain("esdTree");
-    chain->Add("/home/morsch/AliRoot/data/data_jets102/AliESDs.root");
+    chain->Add("~/alice/data/highpt/kPythia6Jets125_150/030/AliESDs.root");
 
     /////////////////////////////////////////////////////////////////////////////////// 
     // Create the analysis manager
@@ -33,22 +35,25 @@ void JetAnalysisManagerLoc()
     mgr->SetInputEventHandler  (inpHandler);
     mgr->SetOutputEventHandler (aodHandler);
     mgr->SetMCtruthEventHandler(mcHandler);
-    mgr-> SetDebugLevel(10);
+    mgr->SetDebugLevel(10);
     /////////////////////////////////////////////////////////////////////////////////// 
     
     
     //
     // Set of cuts
     // 
-    // standard
-    AliESDtrackCuts* esdTrackCutsL = new AliESDtrackCuts("AliESDtrackCuts", "Loose");
+    // standard   
+
+    AliESDtrackCuts* esdTrackCutsL = new AliESDtrackCuts("Standard Track Cuts", "ESD Track Cuts");
     esdTrackCutsL->SetMinNClustersTPC(50);
     esdTrackCutsL->SetMaxChi2PerClusterTPC(3.5);
-    esdTrackCutsL->SetMaxCovDiagonalElements(2,2,0.5,0.5,2);
     esdTrackCutsL->SetRequireTPCRefit(kTRUE);
-    esdTrackCutsL->SetMinNsigmaToVertex(3);
-    esdTrackCutsL->SetRequireSigmaToVertex(kTRUE);
-    esdTrackCutsL->SetAcceptKingDaughters(kFALSE);
+    esdTrackCutsL->SetMaxDCAToVertexXY(2.4);
+    esdTrackCutsL->SetMaxDCAToVertexZ(3.2);
+    esdTrackCutsL->SetDCAToVertex2D(kTRUE);
+    esdTrackCutsL->SetRequireSigmaToVertex(kFALSE);
+    esdTrackCutsL->SetAcceptKinkDaughters(kFALSE);
+    
     //
     // hard
     AliESDtrackCuts* esdTrackCutsH = new AliESDtrackCuts("AliESDtrackCuts", "Hard");
@@ -56,9 +61,9 @@ void JetAnalysisManagerLoc()
     esdTrackCutsH->SetMaxChi2PerClusterTPC(2.0);
     esdTrackCutsH->SetMaxCovDiagonalElements(2,2,0.5,0.5,2);
     esdTrackCutsH->SetRequireTPCRefit(kTRUE);
-    esdTrackCutsH->SetMinNsigmaToVertex(2);
+    esdTrackCutsH->SetMaxNsigmaToVertex(2);
     esdTrackCutsH->SetRequireSigmaToVertex(kTRUE);
-    esdTrackCutsH->SetAcceptKingDaughters(kFALSE);
+    esdTrackCutsH->SetAcceptKinkDaughters(kFALSE);
     //
     AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter");
     trackFilter->AddCuts(esdTrackCutsL);
@@ -71,18 +76,26 @@ void JetAnalysisManagerLoc()
 
     AliAnalysisTaskJets *jetana = new AliAnalysisTaskJets("JetAnalysis");
     jetana->SetDebugLevel(10);
+
+
+
+    AliAnalysisTaskJets *jetanaMC = new AliAnalysisTaskJets("JetAnalysisMC");
+    jetanaMC->SetDebugLevel(10);
+    jetanaMC->SetConfigFile("ConfigJetAnalysisMC.C");
+    jetanaMC->SetNonStdBranch("jetsMC");
+    mgr->AddTask(jetanaMC);
     mgr->AddTask(jetana);
 
     //
     // Create containers for input/output
-    AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("cchain",TChain::Class(), 
-                                                            AliAnalysisManager::kInputContainer);
-
-    AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("tree", TTree::Class(),
-                                                             AliAnalysisManager::kOutputContainer, "default");
+    AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
+    AliAnalysisDataContainer *coutput1 = mgr->GetCommonOutputContainer();
     AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TList::Class(),
                                                              AliAnalysisManager::kOutputContainer, "histos.root");
 
+    AliAnalysisDataContainer *coutputMC2 = mgr->CreateContainer("histosMC", TList::Class(),
+                                                             AliAnalysisManager::kOutputContainer, "histosMC.root");
+
     mgr->ConnectInput  (esdfilter,  0, cinput1  );
     mgr->ConnectOutput (esdfilter,  0, coutput1 );
 
@@ -90,6 +103,10 @@ void JetAnalysisManagerLoc()
     mgr->ConnectOutput (jetana,     0, coutput1 );
     mgr->ConnectOutput (jetana,     1, coutput2 );
 
+    mgr->ConnectInput  (jetanaMC,     0, cinput1  );
+    mgr->ConnectOutput (jetanaMC,     0, coutput1 );
+    mgr->ConnectOutput (jetanaMC,     1, coutputMC2 );
+
 
     //
     // Run the analysis