memory consumption measures only memory used by preprocessor
[u/mrichter/AliRoot.git] / PWG3 / vertexingOld / AliBtoJPSItoEleReco_all.C
1 //--------------------------------------------------------------------------
2 // Test macro for reconstruction and analysis of B->J/psi*X
3 //                                                   |_>e+e-
4 //
5 //     Giuseppe Bruno, giuseppe.bruno@ba.infn.it
6 //      based on the for charm golden channel (D0->Kpi)
7 //--------------------------------------------------------------------------
8
9 void AliBtoJPSItoEleReco_all() {
10   
11   gSystem->Load("libANALYSIS");
12   gSystem->Load("libANALYSISalice");
13   gSystem->Load("libAOD.so");
14   gSystem->Load("libPWG3base.so");
15   gSystem->Load("libPWG3.so");
16   gSystem->Load("libPWG3vertexingOld.so");
17
18   //==============  R E C O N S T R U C T I O N ==============================
19
20   Int_t evFirst = 0;
21   Int_t evLast  = 99;
22
23   // Get field from galice.root
24   if (gAlice) {
25     delete gAlice->GetRunLoader();
26     delete gAlice; 
27     gAlice=0;
28   }  
29   AliRunLoader *rl = AliRunLoader::Open("galice.root");
30   if (rl == 0x0) {
31     cerr<<"Can not open session"<<endl;
32     return;
33   }
34   Int_t retval = rl->LoadgAlice();
35   if (retval) {
36     cerr<<"LoadgAlice returned error"<<endl;
37     delete rl;
38     return;
39   }
40   gAlice=rl->GetAliRun();
41   AliMagF *fiel = (AliMagF*)gAlice->Field();
42   // Set the conversion constant between curvature and Pt
43   AliTracker::SetFieldMap(fiel,kTRUE);
44
45   AliBtoJPSItoEleAnalysis *analysis = new AliBtoJPSItoEleAnalysis();
46   // set simulation to take info on PDG codes from kine tree
47   analysis->SetSimulation();
48   rl->LoadKinematics(); 
49   analysis->MakeTracksRefFile(gAlice,evFirst,evLast);
50   // cout << "ci arrivo 2" << endl;
51   //--- set this is you want only signal candidates in output file
52    analysis->SetOnlySignal();
53   //--- set this if you want to compute primary vertex candidate by candidate using 
54   //    other tracks in the event (for pp, broad interaction region);
55   //    it is time-consuming procedure, so it can be done after a 
56   //    preselection on invariant mass
57   //analysis->SetVertexOnTheFly();
58   //analysis->SetMassCut(0.1); // GeV
59   //--- set single-track preselections
60   analysis->SetPtCut(0.); // GeV
61   analysis->Setd0Cut(0.); // micron
62   //--- set cuts on candidates to be written to file
63   //    (see AliBtoJPSItoEleAnalysis.h for a description and for the defaults)
64   //analysis->SetBCuts(0.1,1000.,1.1,0.,0.,10000.,10000.,0.,.5);
65   analysis->SetBCuts();
66
67   //--- check the settings
68   analysis->PrintStatus();
69
70   analysis->FindCandidates(evFirst,evLast,"AliBtoJPSItoEle_all.root");
71   delete analysis;
72
73   return;
74 }
75 //==========================================================================
76 void AliBtoJPSItoEleSele() {  
77
78   gSystem->Load("libANALYSIS");
79   gSystem->Load("libANALYSISalice");
80   gSystem->Load("libAOD.so");
81   gSystem->Load("libPWG3base.so");
82   gSystem->Load("libPWG3.so");
83   gSystem->Load("libPWG3vertexingOld.so");
84
85   //========================  S E L E C T I O N ============================
86
87   AliBtoJPSItoEleAnalysis *analysis = new AliBtoJPSItoEleAnalysis();
88   analysis->SetSimulation();
89   analysis->SetOnlySignal();
90   analysis->SetBCuts(0.1,1000.,1.1,0.,0.,10000.,10000.,0.,.5);
91   analysis->ApplySelection("AliBtoJPSItoEle.root","AliBtoJPSItoEle_sele.root");
92   delete analysis;
93
94   return;
95 }
96 //==========================================================================
97  
98
99
100
101
102