1 //--------------------------------------------------------------------------
2 // Base macro for submitting AliMUONRecoCheck analysis.
4 // The macro reads ESDs, Kinematics and TrackRefs and outputs file:
6 //--------------------------------------------------------------------------
8 void RunRecoCheck(Bool_t local = kFALSE) {
13 gSystem->Load("libANALYSIS");
14 gSystem->Load("libANALYSISalice");
16 gSystem->Load("libPWG3muondep.so");
18 TChain* chain = new TChain("esdTree");
22 printf("*** Connect to AliEn ***\n");
23 TGrid::Connect("alien://");
24 gSystem->Load("libProofPlayer.so");
26 TAlienCollection* coll = TAlienCollection::Open("wn.xml");
27 TGridResult* result = coll->GetGridResult("",0,0);
28 for(Int_t i = 0; i < result->GetEntries(); i++) {
29 printf("TURL = %s \n",result->GetKey(i,"turl"));
30 chain->Add(result->GetKey(i,"turl"));
33 chain->Add("/path_1_to/AliESDs.root");
34 chain->Add("/path_2_to/AliESDs.root");
35 chain->Add("/path_3_to/AliESDs.root");
38 // set the magnetic field for track extrapolations
39 // look in the OCDB for the value of the magnet current
40 AliCDBManager* cdb = AliCDBManager::Instance();
41 cdb->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
42 //cdb->SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-10-Release/Ideal/");
46 entry = cdb->Get("GRP/GRP/Data");
47 AliGRPObject *obj = (AliGRPObject*)entry->GetObject();
48 Float_t currentL3 = obj->GetL3Current(0);
49 printf("Loading field map for L3 current = %f A ...\n",currentL3);
51 AliLog::SetGlobalLogLevel(AliLog::kError);
53 //____________________________________________//
54 // Make the analysis manager
55 AliAnalysisManager *mgr = new AliAnalysisManager("TestManager");
56 AliVEventHandler* esdH = new AliESDInputHandler;
57 mgr->SetInputEventHandler(esdH);
58 AliMCEventHandler *mc = new AliMCEventHandler();
59 mgr->SetMCtruthEventHandler(mc);
61 //____________________________________________//
63 AliAnalysisTaskRecoCheck *task = new AliAnalysisTaskRecoCheck("TaskRecoCheck");
64 task->SetL3Current(currentL3);
67 // Create containers for input/output
70 AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
73 AliAnalysisDataContainer *coutput = mgr->CreateContainer("ctree", TTree::Class(),AliAnalysisManager::kOutputContainer,"RecoCheck.root");
75 //____________________________________________//
76 mgr->ConnectInput(task,0,cinput);
77 mgr->ConnectOutput(task,0,coutput);
79 if (mgr->InitAnalysis()) {
82 mgr->StartAnalysis("local",chain);