Added functionality to read trigger selection (Ana)

diff --git a/PWG4/GammaConv/AliAnalysisTaskGammaConversion.cxx b/PWG4/GammaConv/AliAnalysisTaskGammaConversion.cxx
index 50a375f48f6dd8d81f27f7d4725c9cc66cf63765..f1d2c14a6741c32b3ff17b90536961ef5de1ea86 100644 (file)
--- a/PWG4/GammaConv/AliAnalysisTaskGammaConversion.cxx
+++ b/PWG4/GammaConv/AliAnalysisTaskGammaConversion.cxx
@@ -61,6 +61,7 @@ AliAnalysisTaskSE(),
   fOutputContainer(NULL),
   fCFManager(0x0),   // for CF
   fHistograms(NULL),
+  fTriggerCINT1B(kFALSE),
   fDoMCTruth(kFALSE),
   fDoNeutralMeson(kFALSE),
   fDoJet(kFALSE),
@@ -134,6 +135,7 @@ AliAnalysisTaskGammaConversion::AliAnalysisTaskGammaConversion(const char* name)
   fOutputContainer(0x0),
   fCFManager(0x0),   // for CF
   fHistograms(NULL),
+  fTriggerCINT1B(kFALSE),
   fDoMCTruth(kFALSE),
   fDoNeutralMeson(kFALSE),
   fDoJet(kFALSE),
@@ -310,7 +312,12 @@ void AliAnalysisTaskGammaConversion::Exec(Option_t */*option*/)
        
   //Clear the data in the v0Reader
   fV0Reader->UpdateEventByEventData();
-       
+
+  //Take Only events with proper trigger
+  if(fTriggerCINT1B){
+    if(!fV0Reader->GetESDEvent()->IsTriggerClassFired("CINT1B-ABCE-NOPF-ALL")) return;
+  }
+
        
   // Process the MC information
   if(fDoMCTruth){
@@ -1232,6 +1239,7 @@ void AliAnalysisTaskGammaConversion::ProcessV0s(){
   }//while(fV0Reader->NextV0)
   fHistograms->FillHistogram("ESD_NumberOfSurvivingV0s", nSurvivingV0s);
   fHistograms->FillHistogram("ESD_NumberOfV0s", fV0Reader->GetNumberOfV0s());
+  fHistograms->FillHistogram("ESD_NumberOfContributorsVtx", fV0Reader->GetNumberOfContributorsVtx());
 
   fV0Reader->ResetV0IndexNumber();
 

diff --git a/PWG4/GammaConv/AliAnalysisTaskGammaConversion.h b/PWG4/GammaConv/AliAnalysisTaskGammaConversion.h
index 6832f905482e9007e285eda550838f81b39a64f2..c1dd08946d8c01b3702046773cfe0e94677a0c26 100644 (file)
--- a/PWG4/GammaConv/AliAnalysisTaskGammaConversion.h
+++ b/PWG4/GammaConv/AliAnalysisTaskGammaConversion.h
@@ -84,6 +84,7 @@ class AliAnalysisTaskGammaConversion : public AliAnalysisTaskSE
 
                
   void SetHistograms(AliGammaConversionHistograms *const histograms){fHistograms=histograms;}
+  void SetTriggerCINT1B(Bool_t flag){fTriggerCINT1B=flag;}
   void SetDoMCTruth(Bool_t flag){fDoMCTruth=flag;}
   void SetDoNeutralMeson(Bool_t flag){fDoNeutralMeson=flag;}
   void SetDoJet(Bool_t flag){fDoJet=flag;}
@@ -156,7 +157,7 @@ class AliAnalysisTaskGammaConversion : public AliAnalysisTaskSE
                
                
   AliGammaConversionHistograms *fHistograms; // Pointer to the histogram handling class
-               
+  Bool_t fTriggerCINT1B; //Flag to select trigger CINT1B
   Bool_t fDoMCTruth; // Flag to switch on/off MC truth 
   Bool_t fDoNeutralMeson; // flag
   Bool_t fDoJet; // flag

diff --git a/PWG4/GammaConv/AliV0Reader.h b/PWG4/GammaConv/AliV0Reader.h
index e7f6ccd5074600d95ea273659bf3aa7907f9cc28..34c2b5aed69368ea8074af8d9a24f5d5d7aac33a 100644 (file)
--- a/PWG4/GammaConv/AliV0Reader.h
+++ b/PWG4/GammaConv/AliV0Reader.h
@@ -96,6 +96,11 @@ class AliV0Reader : public TObject {
    */
   Int_t GetNumberOfV0s() const{return fESDEvent->GetNumberOfV0s();}
        
+  /*
+   *Returns the number of contributors to the vertex
+   */
+  Int_t GetNumberOfContributorsVtx() const{return fESDEvent->GetPrimaryVertex()->GetNContributors();}
+  
   /*
    * Check if there are any more good v0s left in the v0 stack
    * if so, fCurrent v0 is set to this v0 and can be retrieved

diff --git a/PWG4/macros/ConfigGammaConversion.C b/PWG4/macros/ConfigGammaConversion.C
index 623ac5e42c0bfc1aa09d91cb74943f626891b443..04053625be53874f47b4f12566d62442cbca0408 100644 (file)
--- a/PWG4/macros/ConfigGammaConversion.C
+++ b/PWG4/macros/ConfigGammaConversion.C
@@ -338,6 +338,7 @@ Bool_t kGCplotResolutionESDZ  = kTRUE;
 
 Bool_t kGCplotESDNumberOfV0s          = kTRUE;
 Bool_t kGCplotESDNumberOfSurvivingV0s = kTRUE;
+Bool_t kGCplotESDNumberOfContributorsVtx = kTRUE;
 
 //  debug histograms
 Bool_t kGCplotESDCutGetOnFly      = kTRUE;
@@ -685,6 +686,10 @@ Bool_t kGCrunOnTrain = kFALSE;
 Bool_t kGCdoMCTruth = kTRUE;
 /** ---------------------------- end Monte Carlo flag ---------------------------------------*/
 
+/** ------------------------------ Selecting trigger CINT1B -----------------------------------*/
+Bool_t kGCtriggerCINT1B = kFALSE;
+/** ---------------------------- end Monte Carlo flag ---------------------------------------*/
+
 /** ------------------------- Choose KFParticle OR ESDTrack  --------------------------------*/
 Bool_t kGCuseKFParticle = kTRUE;
 Bool_t kGCuseESDTrack   = kFALSE;
@@ -791,6 +796,10 @@ Bool_t scanArguments(TString arguments){
        cout<<"Switching off kGCdoMCTruth"<<endl;
        kGCdoMCTruth = kFALSE;
       }
+      else if (argument.CompareTo("-trigger-CINT1B") == 0){
+        cout<<"Selecting ONLY kGCtriggerCINT1B"<<endl;
+        kGCtriggerCINT1B = kTRUE;
+      }
       else if (argument.CompareTo("-use-own-xyz") == 0){
        cout<<"Switching on use own xyz calculation"<<endl;
        kGCdoOwnXYZCalculation = kTRUE;
@@ -1124,7 +1133,8 @@ AliAnalysisTaskGammaConversion* ConfigGammaConversion(TString arguments,AliAnaly
        
   gammaconversion->SetHistograms(histograms);
   v0Reader->SetHistograms(histograms);// also give the pointer to the v0reader, for debugging cuts
-       
+
+  gammaconversion->SetTriggerCINT1B(kGCtriggerCINT1B);
   gammaconversion->SetDoMCTruth(kGCdoMCTruth);
        
   gammaconversion->SetDoNeutralMeson(kGCrunNeutralMeson);
@@ -1609,6 +1619,7 @@ void AddHistograms(AliGammaConversionHistograms *histograms){
                
     if(kGCplotESDNumberOfV0s == kTRUE){histograms->AddHistogram("ESD_NumberOfV0s","Number of v0s",100, 0, 100,"","");}
     if(kGCplotESDNumberOfSurvivingV0s == kTRUE){histograms->AddHistogram("ESD_NumberOfSurvivingV0s","Number of surviving v0s",100, 0, 100,"","");}
+    if(kGCplotESDNumberOfContributorsVtx == kTRUE){histograms->AddHistogram("ESD_NumberOfContributorsVtx","Number of contributors to vertex",100, 0, 100,"","");}
                
     //  debug histograms
     if(kGCplotESDCutGetOnFly == kTRUE){histograms->AddHistogram("ESD_CutGetOnFly_InvMass" ,"Not GetOnFly" , kGCnXBinsGammaMass, kGCfirstXBinGammaMass, kGClastXBinGammaMass,"","");}