PWGLF/FORWARD/analysis/drawdNdeta.C

   1 void drawdNdeta(const Char_t* filename="fmdana.root",
   2                 Int_t sample = 0,
   3                 Int_t rebin =1,
   4                 Float_t v1 = -10,
   5                 Float_t v2 = 10,
   6                 Bool_t realdata=kTRUE,
   7                 Float_t energy = 900,
   8                 Int_t magfield=1){
   9
  10   gSystem->Load("libANALYSIS");
  11   gSystem->Load("libANALYSISalice");
  12   gSystem->Load("libPWG0base");
  13   gSystem->Load("libPWG0dep");
  14   gSystem->Load("libPWG2forward");
  15   gStyle->SetTextFont(132);
  16   gStyle->SetLabelFont(132,"X");
  17   gStyle->SetLabelFont(132,"Y");
  18   gStyle->SetLabelFont(132,"Z");
  19   gStyle->SetTitleFont(132,"X");
  20   gStyle->SetTitleFont(132,"Y");
  21   gStyle->SetTitleFont(132,"Z");
  22
  23   AliFMDAnaParameters* pars = AliFMDAnaParameters::Instance();
  24   if(energy == 900)
  25     pars->SetEnergy(AliFMDAnaParameters::k900);
  26   else if(energy == 10000)
  27     pars->SetEnergy(AliFMDAnaParameters::k10000);
  28   else if(energy == 7000)
  29     pars->SetEnergy(AliFMDAnaParameters::k7000);
  30   else {
  31     std::cout<<"invalid energy - quitting"<<std::endl;
  32     return; }
  33
  34   if(magfield == 0)
  35     pars->SetMagField(AliFMDAnaParameters::k0G);
  36   else if(magfield == 1)
  37     pars->SetMagField(AliFMDAnaParameters::k5G);
  38   else
  39     return;
  40   pars->SetRealData(realdata);
  41   pars->Init();
  42   pars->Print();
  43   AliFMDDndeta t;
  44   t.SetNbinsToCut(2);
  45   t.SetVtxCut1(v1);
  46   t.SetVtxCut2(v2);
  47   //  t.SetVtxCut(6);
  48   t.Init(filename);
  49   if(sample == 0) {
  50     t.GenerateMult(AliFMDDndeta::kMult);
  51     t.DrawDndeta(AliFMDDndeta::kMult,rebin,realdata);
  52   }
  53   if(sample == 1) {
  54     t.GenerateMult(AliFMDDndeta::kMultTrVtx);
  55     t.DrawDndeta(AliFMDDndeta::kMultTrVtx,rebin,realdata);
  56   }
  57   if(sample == 2) {
  58     t.GenerateMult(AliFMDDndeta::kHits);
  59     t.DrawDndeta(AliFMDDndeta::kHits,rebin);
  60   }
  61   if(sample == 3) {
  62     t.GenerateMult(AliFMDDndeta::kHitsTrVtx);
  63     t.DrawDndeta(AliFMDDndeta::kHitsTrVtx,rebin);
  64   }
  65   if(sample == 4) {
  66     t.GenerateMult(AliFMDDndeta::kMultNSD);
  67     t.DrawDndeta(AliFMDDndeta::kMultNSD,rebin,realdata);
  68   }
  69
  70
  71   if(sample == 10)
  72     t.CreateSharingEfficiency(filename,kTRUE);
  73 }