adding requirement on tzero or vzero if needed

diff --git a/JETAN/AliAnalysisTaskJetCluster.cxx b/JETAN/AliAnalysisTaskJetCluster.cxx
index 3fc92ac90b289d0aedb8b1056108e41ad6b3e046..e83e0f7133df850432fa33cacfe22308a20bb448 100644 (file)
--- 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 7dfee8aec6fcf83abf5132b0201ce24810e24684..0ed7d782c2f1853a06fbcf333b4342ffbf3a98c6 100644 (file)
--- 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 3fc92ac90b289d0aedb8b1056108e41ad6b3e046..e83e0f7133df850432fa33cacfe22308a20bb448 100644 (file)
--- 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 7dfee8aec6fcf83abf5132b0201ce24810e24684..0ed7d782c2f1853a06fbcf333b4342ffbf3a98c6 100644 (file)
--- 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