1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
16 // LHC10x config for He3 and AntiHe3
17 // author: Eulogio Serradilla <eulogio.serradilla@cern.ch>
19 #if !defined(__CINT__) || defined(__MAKECINT__)
22 #include "AliLnDriver.h"
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 = ""
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 )
46 // lhc10b, lhc10c, lhc10d, lhc10e config for He3 and AntiHe3
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;
64 GetTriggerEfficiency(trigEff, kTrigName, period);
66 // input and output filenames
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";
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";
77 // configure the driver and run
81 driver.SetSpecies(kSpecies);
83 driver.SetInputFilenames(inputData, inputSimu, inputSimuFix, inputCorr);
84 driver.SetOutputFilenames(outputPt, outputRatio, outputSpectra);
86 driver.SetRapidityInterval(-ymax,ymax);
88 driver.SetOutputTag(outputTag);
89 driver.SetOutputCorTag(outputTag);
91 driver.SetTriggerEfficiency(trigEff);
92 driver.SetExtrapolateToINEL(inel);
93 driver.SetMCtoINEL(kMCtoINEL);
94 driver.SetPtInterval(ptmin, ptmax);
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);
105 driver.SetMakeStats(makeStats);
106 driver.SetMakeCorrections(makeCor);
107 driver.SetMakePt(makePt);
108 driver.SetMakeRatio(makeRatio);
109 driver.SetMakeSpectra(makeSpectra);
111 driver.SetDebugLevel(kDebugLevel);
117 if(!drawOutput) return 0;
119 DrawOutputCorr(kSpecies, inputCorr, driver.GetOutputCorrTag());
121 DrawPtDebug(driver.GetPtDebugFilename(), outputTag, kSpecies, 0);
122 DrawOutputRatio(outputRatio, outputTag, kSpecies);
123 DrawOutputSpectra(outputSpectra, outputTag, kSpecies);