PWGLF/SPECTRA/Nuclei/B2/macros/Config_He3_TPC_LHC10x.C

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16 // LHC10x config for He3 and AntiHe3
  17 // author: Eulogio Serradilla <eulogio.serradilla@cern.ch>
  18
  19 #if !defined(__CINT__) || defined(__MAKECINT__)
  20 #include <TROOT.h>
  21 #include <TString.h>
  22 #include "AliLnDriver.h"
  23 #endif
  24
  25 #include "Config.h"
  26
  27 Int_t Config_He3_TPC_LHC10x(  const TString& inputDir   = "~/alice/input"
  28                             , const TString& outputDir  = "~/alice/output"
  29                             , const TString& period     = "lhc10bcde"
  30                             , const TString& outputTag  = "pp-mband-7TeV"
  31                             , const TString& trkselTag  = "-tpc3-nsd-moc-vbin"
  32                             , const TString& multTag    = ""
  33                             , const TString& multCorTag = ""
  34                             , Double_t ymax             = 0.5
  35                             , Bool_t   inel             = kFALSE
  36                             , Bool_t   drawOutput       = kTRUE
  37                             , Double_t ptmin            = 0.4  // GeV/c
  38                             , Double_t ptmax            = 10.  // GeV/c
  39                             , Bool_t   makeStats        = kTRUE
  40                             , Bool_t   makeCor          = kTRUE
  41                             , Bool_t   makePt           = kTRUE
  42                             , Bool_t   makeRatio        = kTRUE
  43                             , Bool_t   makeSpectra      = kTRUE )
  44 {
  45 //
  46 // lhc10b, lhc10c, lhc10d, lhc10e config for He3 and AntiHe3
  47 // (TPC)
  48 //
  49         const TString  kSpecies         = "He3";
  50         const TString  kTrkSel          = "its_tpc_dca";
  51         const TString  kTrigName        = "mband";
  52         const Bool_t   kMCtoINEL        = kTRUE;
  53         const Double_t kPidEff          = 1.;
  54         const Bool_t   kSecondaries     = kFALSE;
  55         const Int_t    kSecProc         = 1; // 0 tff, 1 mc
  56         const Int_t    kMatDCAxyMod     = 1; // 0 geant, 1 flat
  57         const Bool_t   kAntiNucTemplate = kFALSE;
  58         const Int_t    kNbin            = 5;
  59         const Double_t kDCAxy[2]        = {-0.2, 0.2};
  60         const Bool_t   kEfficiency      = kTRUE;
  61         const Int_t    kDebugLevel      = 1;
  62
  63         Double_t trigEff[3];
  64         GetTriggerEfficiency(trigEff, kTrigName, period);
  65
  66         // input and output filenames
  67
  68         TString inputData     = inputDir + "/" + period + "/" + MakeInputName(kSpecies, period, kTrkSel+trkselTag+multTag) + ".root";
  69         TString inputSimu     = inputDir + "/" + period + "/" + MakeSimuName(kSpecies, period, kTrkSel+trkselTag+multCorTag) + ".root";
  70         TString inputSimuFix  = inputDir + "/" + period + "/" + MakeSimuFixName(kSpecies, period, kTrkSel+trkselTag+multCorTag) + ".root";
  71         TString inputCorr     = inputDir + "/" + period + "/" + MakeInputName(kSpecies, period, kTrkSel+trkselTag+multTag) + "-corr.root";
  72
  73         TString outputPt      = outputDir + "/" + MakeOutputName(kSpecies, outputTag) + "-Pt.root";
  74         TString outputRatio   = outputDir + "/" + MakeOutputName(kSpecies, outputTag) + "-Ratio.root";
  75         TString outputSpectra = outputDir + "/" + MakeOutputName(kSpecies, outputTag) + "-Spectra.root";
  76
  77         // configure the driver and run
  78
  79         AliLnDriver driver;
  80
  81         driver.SetSpecies(kSpecies);
  82
  83         driver.SetInputFilenames(inputData, inputSimu, inputSimuFix, inputCorr);
  84         driver.SetOutputFilenames(outputPt, outputRatio, outputSpectra);
  85
  86         driver.SetRapidityInterval(-ymax,ymax);
  87
  88         driver.SetOutputTag(outputTag);
  89         driver.SetOutputCorTag(outputTag);
  90
  91         driver.SetTriggerEfficiency(trigEff);
  92         driver.SetExtrapolateToINEL(inel);
  93         driver.SetMCtoINEL(kMCtoINEL);
  94         driver.SetPtInterval(ptmin, ptmax);
  95         driver.SetPid(0);
  96         driver.SetPidEfficiency(kPidEff);
  97         driver.SetSecondaries(kSecondaries);
  98         driver.SetSecProcedure(kSecProc);
  99         driver.SetMatDCAxyModel(kMatDCAxyMod);
 100         driver.SetAntiNucleusAsTemplate(kAntiNucTemplate);
 101         driver.SetNBin(kNbin);
 102         driver.SetDCAxyInterval(kDCAxy[0], kDCAxy[1]);
 103         driver.SetEfficiency(kEfficiency,0);
 104
 105         driver.SetMakeStats(makeStats);
 106         driver.SetMakeCorrections(makeCor);
 107         driver.SetMakePt(makePt);
 108         driver.SetMakeRatio(makeRatio);
 109         driver.SetMakeSpectra(makeSpectra);
 110
 111         driver.SetDebugLevel(kDebugLevel);
 112
 113         driver.Run();
 114
 115         // draw output
 116
 117         if(!drawOutput) return 0;
 118
 119         DrawOutputCorr(kSpecies, inputCorr, driver.GetOutputCorrTag());
 120
 121         DrawPtDebug(driver.GetPtDebugFilename(), outputTag, kSpecies, 0);
 122         DrawOutputRatio(outputRatio, outputTag, kSpecies);
 123         DrawOutputSpectra(outputSpectra, outputTag, kSpecies);
 124
 125         return 0;
 126 }