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