From d66558be10478a2c8218ad0531faf2056ff14676 Mon Sep 17 00:00:00 2001 From: kleinb Date: Wed, 10 Oct 2012 19:58:57 +0000 Subject: [PATCH] adding requirement on tzero or vzero if needed --- JETAN/AliAnalysisTaskJetCluster.cxx | 28 ++++++++++++++++++++++++- JETAN/AliAnalysisTaskJetCluster.h | 4 +++- JETAN/DEV/AliAnalysisTaskJetCluster.cxx | 28 ++++++++++++++++++++++++- JETAN/DEV/AliAnalysisTaskJetCluster.h | 4 +++- 4 files changed, 60 insertions(+), 4 deletions(-) diff --git a/JETAN/AliAnalysisTaskJetCluster.cxx b/JETAN/AliAnalysisTaskJetCluster.cxx index 3fc92ac90b2..e83e0f7133d 100644 --- a/JETAN/AliAnalysisTaskJetCluster.cxx +++ b/JETAN/AliAnalysisTaskJetCluster.cxx @@ -103,6 +103,8 @@ AliAnalysisTaskJetCluster::AliAnalysisTaskJetCluster(): fUseAODMCInput(kFALSE), fUseBackgroundCalc(kFALSE), fEventSelection(kFALSE), + fRequireVZEROAC(kFALSE), + fRequireTZEROvtx(kFALSE), fFilterMask(0), fFilterMaskBestPt(0), fFilterType(0), @@ -267,7 +269,10 @@ AliAnalysisTaskJetCluster::AliAnalysisTaskJetCluster(const char* name): fUseAODTrackInput(kFALSE), fUseAODMCInput(kFALSE), fUseBackgroundCalc(kFALSE), - fEventSelection(kFALSE), fFilterMask(0), + fEventSelection(kFALSE), + fRequireVZEROAC(kFALSE), + fRequireTZEROvtx(kFALSE), + fFilterMask(0), fFilterMaskBestPt(0), fFilterType(0), fJetTypes(kJet), @@ -945,7 +950,28 @@ void AliAnalysisTaskJetCluster::UserExec(Option_t */*option*/) selectEvent = true; } } + + + Bool_t T0 = false; + Bool_t V0 = false; + const AliAODVZERO *vzero = fAOD->GetVZEROData(); + if(vzero){ + if((vzero->GetTriggerChargeA()>0)&&(vzero->GetTriggerChargeC()>0)){ + V0 = true; + } + } + const AliAODTZERO *tzero = fAOD->GetTZEROData(); + if(tzero){ + if(TMath::Abs(tzero->GetT0VertexRaw())<100){ + T0 = true; + } + } + + if(fRequireVZEROAC&&fRequireTZEROvtx)selectEvent = selectEvent&&V0&&T0; + else if(fRequireTZEROvtx)selectEvent = selectEvent&&T0; + else if(fRequireVZEROAC)selectEvent = selectEvent&&V0; + if(!selectEvent){ PostData(1, fHistList); diff --git a/JETAN/AliAnalysisTaskJetCluster.h b/JETAN/AliAnalysisTaskJetCluster.h index 7dfee8aec6f..0ed7d782c2f 100644 --- a/JETAN/AliAnalysisTaskJetCluster.h +++ b/JETAN/AliAnalysisTaskJetCluster.h @@ -148,6 +148,8 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE Bool_t fUseAODMCInput; // take MC from input AOD not from ouptu AOD Bool_t fUseBackgroundCalc; // switches on background calculations Bool_t fEventSelection; // use the event selection of this task, otherwise analyse all + Bool_t fRequireVZEROAC; // switch to require V0 AC + Bool_t fRequireTZEROvtx; // switch to require T0 vtx UInt_t fFilterMask; // filter bit for slecected tracks UInt_t fFilterMaskBestPt; // filter bit to mark jets with high quality leading tracks @@ -324,7 +326,7 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE TList *fHistList; //!leading tracks to be skipped in the randomized event Output list - ClassDef(AliAnalysisTaskJetCluster, 21) + ClassDef(AliAnalysisTaskJetCluster, 22) }; #endif diff --git a/JETAN/DEV/AliAnalysisTaskJetCluster.cxx b/JETAN/DEV/AliAnalysisTaskJetCluster.cxx index 3fc92ac90b2..e83e0f7133d 100644 --- a/JETAN/DEV/AliAnalysisTaskJetCluster.cxx +++ b/JETAN/DEV/AliAnalysisTaskJetCluster.cxx @@ -103,6 +103,8 @@ AliAnalysisTaskJetCluster::AliAnalysisTaskJetCluster(): fUseAODMCInput(kFALSE), fUseBackgroundCalc(kFALSE), fEventSelection(kFALSE), + fRequireVZEROAC(kFALSE), + fRequireTZEROvtx(kFALSE), fFilterMask(0), fFilterMaskBestPt(0), fFilterType(0), @@ -267,7 +269,10 @@ AliAnalysisTaskJetCluster::AliAnalysisTaskJetCluster(const char* name): fUseAODTrackInput(kFALSE), fUseAODMCInput(kFALSE), fUseBackgroundCalc(kFALSE), - fEventSelection(kFALSE), fFilterMask(0), + fEventSelection(kFALSE), + fRequireVZEROAC(kFALSE), + fRequireTZEROvtx(kFALSE), + fFilterMask(0), fFilterMaskBestPt(0), fFilterType(0), fJetTypes(kJet), @@ -945,7 +950,28 @@ void AliAnalysisTaskJetCluster::UserExec(Option_t */*option*/) selectEvent = true; } } + + + Bool_t T0 = false; + Bool_t V0 = false; + const AliAODVZERO *vzero = fAOD->GetVZEROData(); + if(vzero){ + if((vzero->GetTriggerChargeA()>0)&&(vzero->GetTriggerChargeC()>0)){ + V0 = true; + } + } + const AliAODTZERO *tzero = fAOD->GetTZEROData(); + if(tzero){ + if(TMath::Abs(tzero->GetT0VertexRaw())<100){ + T0 = true; + } + } + + if(fRequireVZEROAC&&fRequireTZEROvtx)selectEvent = selectEvent&&V0&&T0; + else if(fRequireTZEROvtx)selectEvent = selectEvent&&T0; + else if(fRequireVZEROAC)selectEvent = selectEvent&&V0; + if(!selectEvent){ PostData(1, fHistList); diff --git a/JETAN/DEV/AliAnalysisTaskJetCluster.h b/JETAN/DEV/AliAnalysisTaskJetCluster.h index 7dfee8aec6f..0ed7d782c2f 100644 --- a/JETAN/DEV/AliAnalysisTaskJetCluster.h +++ b/JETAN/DEV/AliAnalysisTaskJetCluster.h @@ -148,6 +148,8 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE Bool_t fUseAODMCInput; // take MC from input AOD not from ouptu AOD Bool_t fUseBackgroundCalc; // switches on background calculations Bool_t fEventSelection; // use the event selection of this task, otherwise analyse all + Bool_t fRequireVZEROAC; // switch to require V0 AC + Bool_t fRequireTZEROvtx; // switch to require T0 vtx UInt_t fFilterMask; // filter bit for slecected tracks UInt_t fFilterMaskBestPt; // filter bit to mark jets with high quality leading tracks @@ -324,7 +326,7 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE TList *fHistList; //!leading tracks to be skipped in the randomized event Output list - ClassDef(AliAnalysisTaskJetCluster, 21) + ClassDef(AliAnalysisTaskJetCluster, 22) }; #endif -- 2.39.3