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