((TH1D*)fHistoMap->Get(particle + "_PID_Pt"))->Fill(pt);
((TH1D*)fHistoMap->Get(particle + "_PID_Phi"))->Fill(phi);
+ if(iTrack->IsOn(AliESDtrack::kTRDin))
+ {
+ ((TH1D*)fHistoMap->Get(particle + "_PID_TRDin_Pt"))->Fill(pt);
+
+ if(iTrack->IsOn(AliESDtrack::kTRDout)) ((TH1D*)fHistoMap->Get(particle + "_PID_TRDin_TRDout_Pt"))->Fill(pt);
+ if(iTrack->IsOn(AliESDtrack::kTOFout)) ((TH1D*)fHistoMap->Get(particle + "_PID_TRDin_TOFout_Pt"))->Fill(pt);
+ }
+
+ if(iTrack->IsOn(AliESDtrack::kTOFin))
+ {
+ ((TH1D*)fHistoMap->Get(particle + "_PID_TOFin_Pt"))->Fill(pt);
+
+ if(iTrack->IsOn(AliESDtrack::kTOFout)) ((TH1D*)fHistoMap->Get(particle + "_PID_TOFin_TOFout_Pt"))->Fill(pt);
+ }
+
if( this->TOFmatch(iTrack) )
{
((TH2D*)fHistoMap->Get(particle + "_PID_M2_Pt"))->Fill(pt, m2);
+ ((TH1D*)fHistoMap->Get(particle + "_PID_TOFmatch_Pt"))->Fill(pt);
}
// secondaries
// Add task B2
// author: Eulogio Serradilla <eulogio.serradilla@cern.ch>
-Double_t GetMeanNtrk(const TString& period);
-Double_t GetNSDMeanNtrk(const TString& period);
+Double_t GetMeanNtrk(const TString& period, Double_t eta);
+Double_t GetNSDMeanNtrk(const TString& period, Double_t eta);
AliAnalysisTaskB2* AddTaskB2( const TString& species
, const TString& containername
, Double_t minKNOmult = -10
, Double_t maxKNOmult = 10000
, Bool_t V0AND = kFALSE
+ , const TString& ztag = ""
, Double_t maxVz = 10
, Double_t minCentrality = 0
, Double_t maxCentrality = 20
task->SetMaxNSigmaTPC(kMaxNSigmaTPC);
task->SetMaxNSigmaTOF(kMaxNSigmaTOF);
- Double_t meanNtrk = V0AND ? GetNSDMeanNtrk(period) : GetMeanNtrk(period);
+ Double_t meanNtrk = V0AND ? GetNSDMeanNtrk(period, maxEta) : GetMeanNtrk(period, maxEta);
task->SetMeanNtrk(meanNtrk);
task->SetKNOmultInterval(minKNOmult, maxKNOmult);
gROOT->LoadMacro("$ALICE_ROOT/PWGLF/SPECTRA/Nuclei/B2/macros/PriorProbabilities.C");
Double_t prob[9];
- PriorProbabilities(prob, period, trksel);
+ PriorProbabilities(prob, period, trksel, ztag);
AliLnID* lnID = new AliLnID();
return task;
}
-Double_t GetMeanNtrk(const TString& period)
+Double_t GetMeanNtrk(const TString& period, Double_t eta)
{
//
-// average track multiplicity <Ntrk> for the given period
+// average track multiplicity <Ntrk> for the given period and eta
//
- if(period =="lhc10c900") return 3.70158; // pass3
- if(period =="lhc10b_pass2") return 9.86258;
- if(period =="lhc10c_pass2") return 9.61402;
- if(period =="lhc10b") return 5.96104; // pass3
- if(period =="lhc10c") return 5.94719; // pass3
- if(period =="lhc10d") return 5.82333; // pass2
- if(period =="lhc10e") return 5.89367; // pass2
- if(period =="lhc11a_wosdd") return 4.28597; // pass3
- if(period =="lhc11a_wsdd") return 4.69927; // pass4
-
- // MC
- if(period =="lhc10e13") return 3.13712;
- if(period =="lhc10f6a") return 4.41362;
- if(period =="lhc10e21") return 4.74991;
- if(period =="lhc11e3a_plus_wosdd") return 3.37669;
- if(period =="lhc11e3a_plus_wsdd") return 3.47885;
-
- if(period =="lhc12a5a") return 29.264;
- if(period =="lhc12a5bb") return 31.0288;
- if(period =="lhc12a5bc") return 30.6888;
- if(period =="lhc12a5bd") return 30.3528;
- if(period =="lhc12a5be") return 29.9859;
- if(period =="lhc12a5c_wsdd") return 27.5981;
-
- cerr << "Warning in GetMeanNtrk: no <Ntrk> for period " << period << endl;
+ if(TMath::Abs(eta) > 0.51) // |eta|<0.8
+ {
+ if(period =="lhc10d") return 9.47466; // pass2
+ }
+ else // |eta|<0.5
+ {
+ if(period =="lhc10c900") return 3.70158; // pass3
+ if(period =="lhc10b_pass2") return 9.86258;
+ if(period =="lhc10c_pass2") return 9.61402;
+ if(period =="lhc10b") return 5.96104; // pass3
+ if(period =="lhc10c") return 5.94719; // pass3
+ if(period =="lhc10d") return 5.82333; // pass2
+ if(period =="lhc10e") return 5.89367; // pass2
+ if(period =="lhc11a_wosdd") return 4.28597; // pass3
+ if(period =="lhc11a_wsdd") return 4.69927; // pass4
+
+ // MC
+ if(period =="lhc10e13") return 3.13712;
+ if(period =="lhc10f6a") return 4.41362;
+ if(period =="lhc10e21") return 4.74991;
+ if(period =="lhc11e3a_plus_wosdd") return 3.37669;
+ if(period =="lhc11e3a_plus_wsdd") return 3.47885;
+
+ if(period =="lhc12a5a") return 29.264;
+ if(period =="lhc12a5bb") return 31.0288;
+ if(period =="lhc12a5bc") return 30.6888;
+ if(period =="lhc12a5bd") return 30.3528;
+ if(period =="lhc12a5be") return 29.9859;
+ if(period =="lhc12a5c_wsdd") return 27.5981;
+ }
+
+ cerr << "Warning in GetMeanNtrk: no <Ntrk> for period " << period << " and |eta| < " << eta << endl;
return 1;
}
-Double_t GetNSDMeanNtrk(const TString& period)
+Double_t GetNSDMeanNtrk(const TString& period, Double_t eta)
{
//
-// average track multiplicity <Ntrk> for the given period
+// average track multiplicity <Ntrk> for the given period and eta
// (NSD events)
//
- if(period =="lhc10c900") return 3.84362; // pass3
- if(period =="lhc10b") return 6.18470; // pass3
- if(period =="lhc10c") return 6.16175; // pass3
- if(period =="lhc10d") return 6.03108; // pass2
- if(period =="lhc10e") return 6.10384; // pass2
- if(period =="lhc11a_wosdd") return 4.40312; // pass3
- if(period =="lhc11a_wsdd") return 4.87609; // pass4
-
- // MC
- if(period =="lhc10e13") return 3.33273;
- if(period =="lhc10f6a") return 4.80771;
- if(period =="lhc10e21") return 4.91967;
- if(period =="lhc11e3a_plus_wosdd") return 3.4774;
- if(period =="lhc11e3a_plus_wsdd") return 3.6467;
-
- if(period =="lhc12a5a") return 29.3414;
- if(period =="lhc12a5bb") return 31.1514;
- if(period =="lhc12a5bc") return 30.7877;
- if(period =="lhc12a5bd") return 30.4706;
- if(period =="lhc12a5be") return 30.1013;
- if(period =="lhc12a5c_wsdd") return 27.6921;
-
- cerr << "Warning in GetNSDMeanNtrk: no <Ntrk> for period " << period << endl;
+ if(TMath::Abs(eta) > 0.51) // |eta|<0.8
+ {
+ if(period =="lhc10d") return 9.77129; // pass2
+ }
+ else // |eta|<0.5
+ {
+ if(period =="lhc10c900") return 3.84362; // pass3
+ if(period =="lhc10b") return 6.18470; // pass3
+ if(period =="lhc10c") return 6.16175; // pass3
+ if(period =="lhc10d") return 6.03108; // pass2
+ if(period =="lhc10e") return 6.10384; // pass2
+ if(period =="lhc11a_wosdd") return 4.40312; // pass3
+ if(period =="lhc11a_wsdd") return 4.87609; // pass4
+
+ // MC
+ if(period =="lhc10e13") return 3.33273;
+ if(period =="lhc10f6a") return 4.80771;
+ if(period =="lhc10e21") return 4.91967;
+ if(period =="lhc11e3a_plus_wosdd") return 3.4774;
+ if(period =="lhc11e3a_plus_wsdd") return 3.6467;
+
+ if(period =="lhc12a5a") return 29.3414;
+ if(period =="lhc12a5bb") return 31.1514;
+ if(period =="lhc12a5bc") return 30.7877;
+ if(period =="lhc12a5bd") return 30.4706;
+ if(period =="lhc12a5be") return 30.1013;
+ if(period =="lhc12a5c_wsdd") return 27.6921;
+ }
+
+ cerr << "Warning in GetNSDMeanNtrk: no <Ntrk> for period " << period << " and |eta| < " << eta << endl;
return 1;
}
hMap->Add( particle[i] + "_After_Phi_Eta", 400, -2, 2, 360, 0, 2.*TMath::Pi(), "Global tracks (after track cuts)", "#eta", "#phi (rad)");
+ // TRD and TOF absorption
+
+ hMap->Add( particle[i] + "_PID_TRDin_Pt", ptBins, ptMin, ptMax, Form("%s candidates (TRDin)", particle[i].Data(),maxY), "p_{T} (GeV/c)");
+
+ hMap->Add( particle[i] + "_PID_TOFin_Pt", ptBins, ptMin, ptMax, Form("%s candidates (TOFin)", particle[i].Data(),maxY), "p_{T} (GeV/c)");
+
+ hMap->Add( particle[i] + "_PID_TRDin_TRDout_Pt", ptBins, ptMin, ptMax, Form("%s candidates (TRDin-TRDout)", particle[i].Data(),maxY), "p_{T} (GeV/c)");
+
+ hMap->Add( particle[i] + "_PID_TRDin_TOFout_Pt", ptBins, ptMin, ptMax, Form("%s candidates (TRDin-TOFout)", particle[i].Data(),maxY), "p_{T} (GeV/c)");
+
+ hMap->Add( particle[i] + "_PID_TOFin_TOFout_Pt", ptBins, ptMin, ptMax, Form("%s candidates (TOFin-TOFout)", particle[i].Data(),maxY), "p_{T} (GeV/c)");
+
+
// pid with TOF
hMap->Add( particle[i] + "_PID_M2_Pt", ptBins, ptMin, ptMax, m2Bins, m2Min, m2Max, particle[i] + " candidates", "p_{T} (GeV/c)", "m^{2} (GeV^{2}/c^{4})");
hMap->Add( particle[i] + "_PID_Pt", ptBins, ptMin, ptMax, Form("%s candidates (|y| < %0.1f, 0 < \\phi < 2\\pi)", particle[i].Data(),maxY), "p_{T} (GeV/c)");
+ hMap->Add( particle[i] + "_PID_TOFmatch_Pt", ptBins, ptMin, ptMax, Form("%s candidates (|y| < %0.1f, TOFmatch)", particle[i].Data(),maxY), "p_{T} (GeV/c)");
+
hMap->Add( Form("%s_PID_Y",particle[i].Data()), etaBins, etaMin, etaMax, Form("%s candidates (p_{T} > 0, 0 < \\phi < 2\\pi)",particle[i].Data()), "y");
hMap->Add( particle[i] + "_PID_Phi", phiBins, phiMin, phiMax, Form("%s candidates (p_{T} > 0, |y| < %0.1f)",particle[i].Data(),maxY), "\\phi (rad)");
prob[8] = c8;
}
-Bool_t PriorProbabilities(Double_t* prob, const TString& periodname, const TString& tracksel)
+Bool_t PriorProbabilities(Double_t* prob, const TString& periodname, const TString& tracksel, const TString& ztag="")
{
//
// Prior probabilities for the different samples and track selection criteria.
}
else if(trksel=="its_tpc_dca_spd")
{
- SetParameters(prob, 0.00883871, 0.00306025, 0.877528, 0.0646921, 0.0457766, 0.000102032, 2.23315e-06, 3.71013e-08, 0);
+ if(ztag=="ntrk0002")
+ {
+ SetParameters(prob, 0.00572332, 0.00235974, 0.888793, 0.060643, 0.0424185, 6.09002e-05, 1.29575e-06, 0, 0);
+ }
+ else if(ztag=="ntrk0204")
+ {
+ SetParameters(prob, 0.00605528, 0.00176134, 0.888974, 0.0614645, 0.0416846, 5.91396e-05, 1.27997e-06, 0, 0);
+ }
+ else if(ztag=="ntrk0408")
+ {
+ SetParameters(prob, 0.00679582, 0.00129216, 0.883082, 0.0641175, 0.044637, 7.44268e-05, 9.96841e-07, 4.80405e-08, 0);
+ }
+ else if(ztag=="ntrk0811")
+ {
+ SetParameters(prob, 0.00741765, 0.00105688, 0.87915, 0.0661869, 0.046099, 8.77512e-05, 1.61561e-06, 9.07646e-08, 0);
+ }
+ else if(ztag=="ntrk1120")
+ {
+ SetParameters(prob, 0.00806164, 0.000955738, 0.877048, 0.0674902, 0.0463403, 0.00010188, 1.96668e-06, 7.13305e-08, 0);
+ }
+ else if(ztag=="ntrk20xx")
+ {
+ SetParameters(prob, 0.00899995, 0.000887649, 0.876102, 0.0681763, 0.0457112, 0.00011968, 2.87651e-06, 1.16222e-07, 0);
+ }
+ else // MB
+ {
+ SetParameters(prob, 0.00883871, 0.00306025, 0.877528, 0.0646921, 0.0457766, 0.000102032, 2.23315e-06, 3.71013e-08, 0);
+ }
}
else if(trksel=="its_tpc_tof_dca" || trksel == "its_tpc_tof_dca_spd")
{
- SetParameters(prob, 0.0155752, 0.0130319, 0.813626, 0.0808387, 0.0767456, 0.000177436, 4.32771e-06, 3.93428e-07, 0);
+ if(ztag=="ntrk0002")
+ {
+ SetParameters(prob, 0.00354158, 0.0163535, 0.853117, 0.0703677, 0.0565492, 6.95016e-05, 1.68052e-06, 0, 0);
+ }
+ else if(ztag=="ntrk0204")
+ {
+ SetParameters(prob, 0.00383418, 0.0106071, 0.853307, 0.0739765, 0.0582004, 7.27388e-05, 1.58847e-06, 0, 0);
+ }
+ else if(ztag=="ntrk0408")
+ {
+ SetParameters(prob, 0.00457824, 0.00613158, 0.846351, 0.0795435, 0.0632952, 9.89517e-05, 1.19219e-06, 8.72334e-08, 0);
+ }
+ else if(ztag=="ntrk0811")
+ {
+ SetParameters(prob, 0.00517035, 0.00408579, 0.841406, 0.0833815, 0.0658357, 0.000118306, 2.28488e-06, 2.11563e-07, 0);
+ }
+ else if(ztag=="ntrk1120")
+ {
+ SetParameters(prob, 0.00577945, 0.00321126, 0.838286, 0.0861044, 0.0664752, 0.000140696, 3.25918e-06, 1.37712e-07, 0);
+ }
+ else if(ztag=="ntrk20xx")
+ {
+ SetParameters(prob, 0.00666829, 0.00269162, 0.836384, 0.0883033, 0.0657773, 0.000170996, 4.77082e-06, 1.88322e-07, 0);
+ }
+ else // MB
+ {
+ SetParameters(prob, 0.0155752, 0.0130319, 0.813626, 0.0808387, 0.0767456, 0.000177436, 4.32771e-06, 3.93428e-07, 0);
+ }
}
else // its_tpc_nsigma
{