[u/mrichter/AliRoot.git] / PWGLF / SPECTRA / Nuclei / B2 / macros / Config_Deuteron_TOF_LHC10x.C

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

// LHC10x config for deuterons and antideuterons
// author: Eulogio Serradilla <eulogio.serradilla@cern.ch>

#if !defined(__CINT__) || defined(__MAKECINT__)
#include <TROOT.h>
#include <TString.h>
#include "AliLnDriver.h"
#endif

#include "Config.h"

Int_t Config_Deuteron_TOF_LHC10x(  const TString& inputDir    = "~/alice/input"
                                 , const TString& outputDir   = "~/alice/output"
                                 , const TString& period      = "lhc10bcde"
                                 , const TString& outputTag   = "lhc10bcde"
                                 , const TString& trkselTag   = "-tpc3-nsd-moc"
                                 , const TString& multTag     = ""
                                 , const TString& multCorTag  = ""
                                 , Double_t       ymax        = 0.5
                                 , Bool_t         inel        = 0  // for mult
                                 , Bool_t         drawOutput  = 1  // for batch
                                 , Double_t       ptmin       = 0.7
                                 , Double_t       ptmax       = 3.0
                                 , Double_t       ptpid       = 0.8
                                 , Bool_t         makeStats   = 1
                                 , Bool_t         makeCor     = 1
                                 , Bool_t         makePt      = 1
                                 , Bool_t         makeRatio   = 1
                                 , Bool_t         makeSpectra = 1)
{
//
// lhc10b, lhc10c, lhc10d, lhc10e config for deuterons and antideuterons
//
	const TString  kSpecies         = "Deuteron";
	const TString  kTrkSel          = "its_tpc_tof_dca";
	const TString  kTrigName        = "mband";
	const Bool_t   kMCtoINEL        = 1;
	const Bool_t   kPid             = 1;
	const Int_t    kPidProc         = 0; // 0 m2, 1 dm2, 2 time
	const Double_t kPidEff          = 1.;
	const Bool_t   kSecondaries     = 1;
	const Int_t    kSecProc         = 0; // 0 tff, 1 mc
	const Int_t    kMatDCAxyMod     = 1; // 0 geant, 1 flat
	const Bool_t   kAntiNucTemplate = 0;
	const Int_t    kNbin            = 5;
	const Double_t kDCAxy[2]        = {-0.2,0.2};
	const Bool_t   kEfficiency      = 1;
	const Bool_t   kFitFrac         = 1;
	const Int_t    kDebugLevel      = 1;
	
	Double_t bkgLimitToF[2]    = {-2, 2};
	Double_t pidLimitToF[2]    = {-2.,6.};
	
	Double_t bkgLimitM2[2]     = {1.8,6.0};
	Double_t pidLimitM2[2]     = {1.8,6.};
	
	if(kPidProc==1)
	{
		for(Int_t i=0; i<2; ++i)
		{
			bkgLimitM2[i] -= 3.51792;
			pidLimitM2[i] -= 3.51792;
		}
	}
	
	Double_t trigEff[3];
	GetTriggerEfficiency(trigEff, kTrigName, period);
	
	// input and output filenames
	
	TString inputData     = inputDir + "/" + period + "/" + MakeInputName(kSpecies, period, kTrkSel+trkselTag+multTag) + ".root";
	TString inputSimu     = inputDir + "/" + period + "/" + MakeSimuName(kSpecies, period, kTrkSel+trkselTag+multCorTag) + ".root";
	TString inputSimuFix  = inputDir + "/" + period + "/" + MakeSimuFixName(kSpecies, period, kTrkSel+trkselTag+multCorTag) + ".root";
	TString inputCorr     = inputDir + "/" + period + "/" + MakeInputName(kSpecies, period, kTrkSel+trkselTag+multTag) + "-corr.root";
	
	TString outputPt      = outputDir + "/" + MakeOutputName(kSpecies, outputTag) + "-Pt.root";
	TString outputRatio   = outputDir + "/" + MakeOutputName(kSpecies, outputTag) + "-Ratio.root";
	TString outputSpectra = outputDir + "/" + MakeOutputName(kSpecies, outputTag) + "-Spectra.root";
	
	// configure the driver and run
	
	AliLnDriver driver;
	
	driver.SetSpecies(kSpecies);
	
	driver.SetInputFilenames(inputData, inputSimu, inputSimuFix, inputCorr);
	driver.SetOutputFilenames(outputPt, outputRatio, outputSpectra);
	
	driver.SetRapidityInterval(-ymax,ymax);
	
	driver.SetOutputTag(outputTag);
	driver.SetOutputCorTag(outputTag);
	
	driver.SetTriggerEfficiency(trigEff);
	driver.SetExtrapolateToINEL(inel);
	driver.SetMCtoINEL(kMCtoINEL);
	driver.SetPtInterval(ptmin, ptmax);
	driver.SetPid(kPid);
	driver.SetPidProcedure(kPidProc);
	driver.SetPidEfficiency(kPidEff);
	driver.SetPidPt(ptpid);
	
	if(kPidProc==2)
	{
		driver.SetBkgInterval(bkgLimitToF[0], bkgLimitToF[1]);
		driver.SetPidInterval(pidLimitToF[0], pidLimitToF[1]);
	}
	else
	{
		driver.SetBkgInterval(bkgLimitM2[0], bkgLimitM2[1]);
		driver.SetPidInterval(pidLimitM2[0], pidLimitM2[1]);
	}
	
	driver.SetSecondaries(kSecondaries);
	driver.SetSecProcedure(kSecProc);
	driver.SetMatDCAxyModel(kMatDCAxyMod);
	driver.SetAntiNucleusAsTemplate(kAntiNucTemplate);
	driver.SetNBin(kNbin);
	driver.SetDCAxyInterval(kDCAxy[0], kDCAxy[1]);
	driver.SetEfficiency(kEfficiency,0);
	driver.SetFitFractionCorr(kFitFrac);
	
	driver.SetDebugLevel(kDebugLevel);
	
	driver.SetMakeStats(makeStats);
	driver.SetMakeCorrections(makeCor);
	driver.SetMakePt(makePt);
	driver.SetMakeRatio(makeRatio);
	driver.SetMakeSpectra(makeSpectra);
	
	driver.Run();
	
	// draw output
	
	if(!drawOutput) return 0;
	
	DrawOutputCorr(kSpecies, inputCorr, driver.GetOutputCorrTag());
	
	if(kSecProc == 0) gROOT->ProcessLine(Form(".x DrawSec.C+g(\"%s\",\"%s\",\"Deuteron\", %f, %f, %f, %f)", driver.GetPtCorrDebugFilename().Data(), driver.GetOutputCorrTag().Data(), ptmin, ptmax, kDCAxy[0], kDCAxy[1]));
	
	DrawPtDebug(driver.GetPtDebugFilename(), outputTag, kSpecies, kPid, ptmax, ptpid);
	DrawOutputRatio(outputRatio, outputTag, kSpecies);
	DrawOutputSpectra(outputSpectra, outputTag, kSpecies);
	
	return 0;
}
Commit	Line	Data
c683985a	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 deuterons and antideuterons
	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_Deuteron_TOF_LHC10x( const TString& inputDir = "~/alice/input"
	28	, const TString& outputDir = "~/alice/output"
	29	, const TString& period = "lhc10bcde"
	30	, const TString& outputTag = "lhc10bcde"
	31	, const TString& trkselTag = "-tpc3-nsd-moc"
	32	, const TString& multTag = ""
	33	, const TString& multCorTag = ""
	34	, Double_t ymax = 0.5
	35	, Bool_t inel = 0 // for mult
	36	, Bool_t drawOutput = 1 // for batch
	37	, Double_t ptmin = 0.7
	38	, Double_t ptmax = 3.0
	39	, Double_t ptpid = 0.8
	40	, Bool_t makeStats = 1
	41	, Bool_t makeCor = 1
	42	, Bool_t makePt = 1
	43	, Bool_t makeRatio = 1
	44	, Bool_t makeSpectra = 1)
	45	{
	46	//
	47	// lhc10b, lhc10c, lhc10d, lhc10e config for deuterons and antideuterons
	48	//
	49	const TString kSpecies = "Deuteron";
	50	const TString kTrkSel = "its_tpc_tof_dca";
	51	const TString kTrigName = "mband";
	52	const Bool_t kMCtoINEL = 1;
	53	const Bool_t kPid = 1;
	54	const Int_t kPidProc = 0; // 0 m2, 1 dm2, 2 time
	55	const Double_t kPidEff = 1.;
	56	const Bool_t kSecondaries = 1;
	57	const Int_t kSecProc = 0; // 0 tff, 1 mc
	58	const Int_t kMatDCAxyMod = 1; // 0 geant, 1 flat
	59	const Bool_t kAntiNucTemplate = 0;
	60	const Int_t kNbin = 5;
	61	const Double_t kDCAxy[2] = {-0.2,0.2};
	62	const Bool_t kEfficiency = 1;
	63	const Bool_t kFitFrac = 1;
	64	const Int_t kDebugLevel = 1;
65
66	Double_t bkgLimitToF[2] = {-2, 2};
67	Double_t pidLimitToF[2] = {-2.,6.};
68
69	Double_t bkgLimitM2[2] = {1.8,6.0};
70	Double_t pidLimitM2[2] = {1.8,6.};
71
72	if(kPidProc==1)
73	{
74	for(Int_t i=0; i<2; ++i)
75	{
76	bkgLimitM2[i] -= 3.51792;
77	pidLimitM2[i] -= 3.51792;
78	}
79	}
80
81	Double_t trigEff[3];
82	GetTriggerEfficiency(trigEff, kTrigName, period);
83
84	// input and output filenames
85
86	TString inputData = inputDir + "/" + period + "/" + MakeInputName(kSpecies, period, kTrkSel+trkselTag+multTag) + ".root";
87	TString inputSimu = inputDir + "/" + period + "/" + MakeSimuName(kSpecies, period, kTrkSel+trkselTag+multCorTag) + ".root";
88	TString inputSimuFix = inputDir + "/" + period + "/" + MakeSimuFixName(kSpecies, period, kTrkSel+trkselTag+multCorTag) + ".root";
89	TString inputCorr = inputDir + "/" + period + "/" + MakeInputName(kSpecies, period, kTrkSel+trkselTag+multTag) + "-corr.root";
90
91	TString outputPt = outputDir + "/" + MakeOutputName(kSpecies, outputTag) + "-Pt.root";
92	TString outputRatio = outputDir + "/" + MakeOutputName(kSpecies, outputTag) + "-Ratio.root";
93	TString outputSpectra = outputDir + "/" + MakeOutputName(kSpecies, outputTag) + "-Spectra.root";
94
95	// configure the driver and run
96
97	AliLnDriver driver;
98
99	driver.SetSpecies(kSpecies);
100
101	driver.SetInputFilenames(inputData, inputSimu, inputSimuFix, inputCorr);
102	driver.SetOutputFilenames(outputPt, outputRatio, outputSpectra);
103
104	driver.SetRapidityInterval(-ymax,ymax);
105
106	driver.SetOutputTag(outputTag);
107	driver.SetOutputCorTag(outputTag);
108
109	driver.SetTriggerEfficiency(trigEff);
110	driver.SetExtrapolateToINEL(inel);
111	driver.SetMCtoINEL(kMCtoINEL);
112	driver.SetPtInterval(ptmin, ptmax);
113	driver.SetPid(kPid);
114	driver.SetPidProcedure(kPidProc);
115	driver.SetPidEfficiency(kPidEff);
116	driver.SetPidPt(ptpid);
117
118	if(kPidProc==2)
119	{
120	driver.SetBkgInterval(bkgLimitToF[0], bkgLimitToF[1]);
121	driver.SetPidInterval(pidLimitToF[0], pidLimitToF[1]);
122	}
123	else
124	{
125	driver.SetBkgInterval(bkgLimitM2[0], bkgLimitM2[1]);
126	driver.SetPidInterval(pidLimitM2[0], pidLimitM2[1]);
127	}
128
129	driver.SetSecondaries(kSecondaries);
130	driver.SetSecProcedure(kSecProc);
131	driver.SetMatDCAxyModel(kMatDCAxyMod);
132	driver.SetAntiNucleusAsTemplate(kAntiNucTemplate);
133	driver.SetNBin(kNbin);
134	driver.SetDCAxyInterval(kDCAxy[0], kDCAxy[1]);
135	driver.SetEfficiency(kEfficiency,0);
136	driver.SetFitFractionCorr(kFitFrac);
137
138	driver.SetDebugLevel(kDebugLevel);
139
140	driver.SetMakeStats(makeStats);
141	driver.SetMakeCorrections(makeCor);
142	driver.SetMakePt(makePt);
143	driver.SetMakeRatio(makeRatio);
144	driver.SetMakeSpectra(makeSpectra);
145
146	driver.Run();
147
148	// draw output
149
150	if(!drawOutput) return 0;
151
152	DrawOutputCorr(kSpecies, inputCorr, driver.GetOutputCorrTag());
153
154	if(kSecProc == 0) gROOT->ProcessLine(Form(".x DrawSec.C+g(\"%s\",\"%s\",\"Deuteron\", %f, %f, %f, %f)", driver.GetPtCorrDebugFilename().Data(), driver.GetOutputCorrTag().Data(), ptmin, ptmax, kDCAxy[0], kDCAxy[1]));
155
156	DrawPtDebug(driver.GetPtDebugFilename(), outputTag, kSpecies, kPid, ptmax, ptpid);
157	DrawOutputRatio(outputRatio, outputTag, kSpecies);
158	DrawOutputSpectra(outputSpectra, outputTag, kSpecies);
159
160	return 0;
161	}