#include "AliESDtrackCuts.h"
#include "AliCentrality.h"
+#include "AliESDtrack.h"
+#include "AliESDpid.h"
+
#include "AliEbyEFluctuationAnalysisTask.h"
// Event by event charge fluctuation analysis
// Authors: Satyajit Jena and Panos Cristakoglou
+// (PID by R+)
//
ClassImp(AliEbyEFluctuationAnalysisTask)
fESDtrackCuts(0),
fAnalysisType("ESD"), fAnalysisMode("TPC"),
fCentralityEstimator("V0M"), fCentralityBins20(kFALSE),
- fVxMax(3.), fVyMax(3.), fVzMax(10.) {
+ fVxMax(3.), fVyMax(3.), fVzMax(10.),
+ fESDpid(NULL) {
// Constructor
+ fOutputList = new TList();
+ fOutputList->SetOwner();
+
for(Int_t iBin = 1; iBin <= nCentralityBins; iBin++) {
fHistNMult[iBin-1] = NULL;
fHistNPlusNMinus[iBin-1] = NULL;
// Define input and output slots here
// Input slot #0 works with a TChain
- DefineInput(0, TChain::Class());
+ // DefineInput(0, TChain::Class());
// Output slot #0 id reserved by the base class for AOD
// Output slot #1 writes into a TH1 container
- DefineOutput(1, TList::Class());
+ DefineOutput(2, TList::Class());
}
//________________________________________________________________________
// Create histograms
// Called once
- fOutputList = new TList();
- fOutputList->SetOwner();
+ printf("in UserCreateOutputObjects\n");
//Event stats.
TString gCutName[4] = {"Total","Offline trigger",
3000,0.5,3000.5,3000,0.5,3000.5);
fOutputList->Add(fHistNPlusNMinusMC);
}//MC analysis
+
+ PostData(2, fOutputList);
}
//________________________________________________________________________
Short_t gCharge = track->Charge();
if(gCharge > 0) nPlus += 1;
if(gCharge < 0) nMinus += 1;
+
+
}//track loop
}//TObjArray valid object
//if((nCentrality >= 1)&&(nCentrality <= 20)) {
}//MC analysis level
- PostData(1, fOutputList);
+ PostData(2, fOutputList);
}
//________________________________________________________________________
return;
}
}
+
+//___________________________________________________________
+
+Bool_t
+AliEbyEFluctuationAnalysisTask::HasTPCPID(AliESDtrack *track) const
+{
+ /*
+ * has TPC PID
+ */
+
+ /* check PID signal */
+ if (track->GetTPCsignal() <= 0. ||
+ track->GetTPCsignalN() == 0 ||
+ !track->GetInnerParam()) return kFALSE;
+ return kTRUE;
+}
+
+//___________________________________________________________
+
+Bool_t
+AliEbyEFluctuationAnalysisTask::HasTOFPID(AliESDtrack *track) const
+{
+ /*
+ * has TOF PID
+ */
+
+ /* check TOF matched track */
+ if (!(track->GetStatus() & AliESDtrack::kTOFout)||
+ !(track->GetStatus() & AliESDtrack::kTIME)) return kFALSE;
+ /* check integrated length */
+ if (track->GetIntegratedLength() < 350.) return kFALSE;
+ return kTRUE;
+}
+
+//___________________________________________________________
+
+Double_t
+AliEbyEFluctuationAnalysisTask::MakeTPCPID(AliESDtrack *track, Double_t *nSigma) const
+{
+ /*
+ * make TPC PID
+ * returns measured dEdx if PID available, otherwise -1.
+ * fills nSigma array with TPC nsigmas for e, mu, pi, K, p
+ */
+
+ /* check TPC PID */
+ if (!HasTPCPID(track)) return -1.;
+
+ /* get TPC info */
+ Double_t ptpc = track->GetInnerParam() ? track->GetInnerParam()->P() : 0.;
+ Double_t dEdx = track->GetTPCsignal();
+ Double_t dEdxN = track->GetTPCsignalN();
+
+ /* loop over particles */
+ for (Int_t ipart = 0; ipart < AliPID::kSPECIES; ipart++) {
+ Double_t bethe = fESDpid->GetTPCResponse().GetExpectedSignal(ptpc, (AliPID::EParticleType)ipart);
+ Double_t diff = dEdx - bethe;
+ Double_t sigma = fESDpid->GetTPCResponse().GetExpectedSigma(ptpc, dEdxN, (AliPID::EParticleType)ipart);
+ nSigma[ipart] = diff / sigma;
+ }
+
+ return dEdx;
+}
+
+//___________________________________________________________
+
+Double_t
+AliEbyEFluctuationAnalysisTask::MakeTOFPID(AliESDtrack *track, Double_t *nSigma) const
+{
+ /*
+ * make TOF PID
+ * returns measured beta if PID available, otherwise -1.
+ * fills nSigma array with TOF nsigmas for e, mu, pi, K, p
+ */
+
+ /* check TOF PID */
+ if (!HasTOFPID(track)) return -1.;
+
+ /* get TOF info */
+ Double_t p = track->P();
+ Double_t time = track->GetTOFsignal() - fESDpid->GetTOFResponse().GetStartTime(p);
+ Double_t length = track->GetIntegratedLength();
+ Double_t timei[5];
+ track->GetIntegratedTimes(timei);
+
+ /* loop over particles */
+ for (Int_t ipart = 0; ipart < AliPID::kSPECIES; ipart++) {
+ Double_t timez = time - timei[ipart];
+ Double_t sigma = fESDpid->GetTOFResponse().GetExpectedSigma(p, timei[ipart], AliPID::ParticleMass(ipart));
+ nSigma[ipart] = timez / sigma;
+ }
+
+ return length / (time * 2.99792457999999984e-02);
+}
+
+//___________________________________________________________
+
+void
+AliEbyEFluctuationAnalysisTask::MakePID(AliESDtrack *track, Bool_t *pidFlag) const
+{
+ /*
+ * make PID
+ * fills PID QA plots
+ * fills pidFlag array with PID flags for e, mu, pi, K, p
+ */
+
+ /* cut definitions
+ (better put them as static variables so they can be changed from outside) */
+ Double_t fgTPCPIDmomcut[AliPID::kSPECIES] = {0., 0., 0.5, 0.5, 0.7};
+ Double_t fgTPCPIDsigmacut[AliPID::kSPECIES] = {0., 0., 2., 2., 2.};
+ Double_t fgTPCPIDcompcut[AliPID::kSPECIES] = {0., 0., 3., 3., 3.};
+ Double_t fgTOFPIDmomcut[AliPID::kSPECIES] = {0., 0., 1.5, 1.5, 2.0};
+ Double_t fgTOFPIDsigmacut[AliPID::kSPECIES] = {0., 0., 2., 2., 2.};
+
+ /* make pid and check if available */
+ Double_t p = track->P();
+ Double_t pt = track->Pt();
+ Double_t ptpc = track->GetInnerParam() ? track->GetInnerParam()->P() : 0.;
+ Double_t nsigmaTPC[AliPID::kSPECIES];
+ Double_t nsigmaTOF[AliPID::kSPECIES];
+ Double_t dEdx = MakeTPCPID(track, nsigmaTPC);
+ Double_t beta = MakeTOFPID(track, nsigmaTPC);
+ Bool_t hasTPCPID = dEdx > 0.;
+ Bool_t hasTOFPID = beta > 0.;
+
+ /* fill qa histos */
+ if (hasTPCPID)
+ fHistoTPCdEdx->Fill(ptpc, dEdx);
+ if (hasTOFPID)
+ fHistoTOFbeta->Fill(p, beta);
+
+ /* loop over species */
+ for (Int_t ipart = 0; ipart < AliPID::kSPECIES; ipart++) {
+
+ /* reset PID flag */
+ pidFlag[ipart] = kFALSE;
+
+ /* fill qa histos */
+ if (hasTPCPID)
+ fHistoNSigmaTPC[ipart]->Fill(pt, nsigmaTPC[ipart]);
+ if (hasTOFPID)
+ fHistoNSigmaTOF[ipart]->Fill(pt, nsigmaTOF[ipart]);
+
+ /* combined PID cuts */
+ if (hasTPCPID && hasTOFPID) {
+ if (pt < fgTOFPIDmomcut[ipart] &&
+ TMath::Abs(nsigmaTOF[ipart]) < fgTOFPIDsigmacut[ipart] &&
+ TMath::Abs(nsigmaTPC[ipart]) < fgTPCPIDcompcut[ipart])
+ pidFlag[ipart] = kTRUE;
+
+ }
+ /* TPC-only PID cuts */
+ else if (hasTPCPID && !hasTOFPID) {
+ if (pt < fgTPCPIDmomcut[ipart] &&
+ TMath::Abs(nsigmaTPC[ipart]) < fgTPCPIDsigmacut[ipart])
+ pidFlag[ipart] = kTRUE;
+ }
+ /* TOF-only PID cuts */
+ else if (!hasTPCPID && hasTOFPID) {
+ if (pt < fgTOFPIDmomcut[ipart] &&
+ TMath::Abs(nsigmaTOF[ipart]) < fgTOFPIDsigmacut[ipart])
+ pidFlag[ipart] = kTRUE;
+ }
+
+ } /* end of loop over species */
+
+}
+
+//___________________________________________________________
+
+void
+AliEbyEFluctuationAnalysisTask::InitPID(AliESDEvent *event)
+{
+ /*
+ * init PID
+ */
+
+ /* create ESDpid object if not there yet */
+ if (!fESDpid) {
+ /* instance object */
+ Bool_t mcFlag = kFALSE; /*** WARNING: check whether is MC ***/
+ fESDpid = new AliESDpid(mcFlag);
+ /* set OADB path */
+ fESDpid->SetOADBPath("$ALICE_ROOT/OADB");
+ }
+
+ /* init ESD PID */
+ Int_t recoPass = 2; /*** WARNING: need to set the recoPass somehow better ***/
+ fESDpid->InitialiseEvent(event, recoPass); /* warning: this apparently sets TOF time
+ * resolution to some hardcoded value,
+ * therefore we have to set correct
+ * resolution value after this call */
+
+ /* set TOF resolution */
+ Double_t tofReso = 85.; /* ps */ /*** WARNING: need to set tofReso somehow better ***/
+ fESDpid->GetTOFResponse().SetTimeResolution(tofReso);
+
+}
+
+
+
+++ /dev/null
-# -*- mode: CMake -*-
-#--------------------------------------------------------------------------------#
-# Package File for PWGUDdep #
-# Author : Johny Jose (johny.jose@cern.ch) #
-# Variables Defined : #
-# #
-# SRCS - C++ source files #
-# HDRS - C++ header files #
-# DHDR - ROOT Dictionary Linkdef header file #
-# CSRCS - C source files #
-# CHDRS - C header files #
-# EINCLUDE - Include directories #
-# EDEFINE - Compiler definitions #
-# ELIBS - Extra libraries to link #
-# ELIBSDIR - Extra library directories #
-# PACKFFLAGS - Fortran compiler flags for package #
-# PACKCXXFLAGS - C++ compiler flags for package #
-# PACKCFLAGS - C compiler flags for package #
-# PACKSOFLAGS - Shared library linking flags #
-# PACKLDFLAGS - Module linker flags #
-# PACKBLIBS - Libraries to link (Executables only) #
-# EXPORT - Header files to be exported #
-# CINTHDRS - Dictionary header files #
-# CINTAUTOLINK - Set automatic dictionary generation #
-# ARLIBS - Archive Libraries and objects for linking (Executables only) #
-# SHLIBS - Shared Libraries and objects for linking (Executables only) #
-#--------------------------------------------------------------------------------#
-
-set ( SRCS
- dNdPt/AlidNdPtHelper.cxx
- dNdPt/AlidNdPtAnalysis.cxx
- dNdPt/AlidNdPtAnalysisPbPb.cxx
- dNdPt/AlidNdPtCorrection.cxx
- dNdPt/AlidNdPtAcceptanceCuts.cxx
- dNdPt/AlidNdPtEventCuts.cxx
- dNdPt/AlidNdPt.cxx
- dNdPt/AlidNdPtCutAnalysis.cxx
- dNdPt/AlidNdPtTask.cxx
- dNdPt/AlidNdPtBackgroundCuts.cxx
- dNdPt/AlidNdPtCutAnalysisPbPb.cxx
- dNdPt/AlidNdPtEfficiency.cxx
- dNdPt/AliPtResolAnalysis.cxx
- dNdPt/AliPtResolAnalysisPbPb.cxx
- dNdPt/AlidNdPtTrackDumpTask.cxx
- )
-
-string ( REPLACE ".cxx" ".h" HDRS "${SRCS}" )
-
-set ( DHDR PWGUDdepLinkDef.h)
-
-set ( EINCLUDE PWGUD/base TPC STEER/STEER STEER/ESD STEER/STEERBase)
git://git.uio.no / u / mrichter / AliRoot.git / commitdiff