1 void AliAnalysisTaskSigma1385()
6 printf("*** Connect to AliEn ***\n");
7 TGrid::Connect("alien://");
9 gSystem->Load("libTree.so");
10 gSystem->Load("libGeom.so");
11 gSystem->Load("libVMC.so");
12 gSystem->Load("libPhysics.so");
13 //____________________________________________________//
14 //_____________Setting up required packages___________//
15 //____________________________________________________//
16 gSystem->Load("libSTEERBase.so");
17 gSystem->Load("libSTEER.so");
18 gSystem->Load("libESD.so");
19 gSystem->Load("libAOD.so");
20 gSystem->Load("libANALYSIS.so");
21 gSystem->Load("libANALYSISalice.so");
22 gROOT->ProcessLine(".include $ALICE_ROOT/include");
23 gROOT->ProcessLine(".include $ALICE_ROOT/TOF");
25 gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPhysicsSelection.C");
26 gROOT->LoadMacro("AliAnalysisTaskSigma1385.cxx+g");
28 const char *collectionfile = "sigmaLHC10d1.xml";
31 TAlienCollection * myCollection = TAlienCollection::Open(collectionfile);
34 Error("AliAnalysisTaskSigma1385", Form("Cannot create an AliEn collection from %s", collectionfile)) ;
38 Info("AliAnalysisTaskSigma1385", Form("Creating the analysis chain %s", "esdTree")) ;
39 TChain* chain = new TChain("esdTree");
41 Info("AliAnalysisTaskSigma1385", "Preparing the file list") ;
42 myCollection->Reset() ;
43 while (myCollection->Next()) {
45 sprintf(esdFile, "%s", myCollection->GetTURL("")) ;
46 Info("AliAnalysisTaskSigma1385", Form("Adding %s", esdFile)) ;
51 //____________________________________________//
52 // Make the analysis manager
53 AliAnalysisManager *mgr = new AliAnalysisManager("TestManager");
54 AliVEventHandler* esdH = new AliESDInputHandler;
55 mgr->SetInputEventHandler(esdH);
58 //____________________________________________//
60 //AliPhysicsSelectionTask* physSelTask = AddTaskPhysicsSelection();
61 AliPhysicsSelectionTask* physSelTask = AddTaskPhysicsSelection(kTRUE);
62 AliAnalysisTaskSigma1385 *task1 = new AliAnalysisTaskSigma1385("TaskSigma1385");
67 task1->SetDataType("SIM");
68 task1->SelectCollisionCandidates();
70 if (!(strcmp(task1->GetDataType(), "SIM"))) {
72 AliMCEventHandler *mc = new AliMCEventHandler();
73 mc->SetReadTR(kFALSE);
74 mgr->SetMCtruthEventHandler(mc);
76 mgr->SetDebugLevel(10);
80 //------ input data ------
81 AliAnalysisDataContainer *cinput0 = mgr->GetCommonInputContainer();
83 // ----- output data -----
85 AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("chist1", TList::Class(), AliAnalysisManager::kOutputContainer, "Sigma.ESD.root");
87 //____________________________________________//
89 mgr->ConnectInput(task1, 0, cinput0);
90 mgr->ConnectOutput(task1, 1, coutput1);
91 if (!mgr->InitAnalysis()) return;
93 mgr->StartAnalysis("local", chain);