add qa check

diff --git a/PWG/EMCAL/AliAnalysisTaskEMCALClusterizeFast.cxx b/PWG/EMCAL/AliAnalysisTaskEMCALClusterizeFast.cxx
index 7dfd45508f0ef0790ce6e180b1179b54265bb702..1e304629a47b8e20c413100bd6585d49373efe5e 100644 (file)
--- a/PWG/EMCAL/AliAnalysisTaskEMCALClusterizeFast.cxx
+++ b/PWG/EMCAL/AliAnalysisTaskEMCALClusterizeFast.cxx
@@ -512,11 +512,12 @@ Bool_t AliAnalysisTaskEMCALClusterizeFast::AcceptCell(Int_t cellNumber) {
   Bool_t accept = kTRUE;
   if(fRejectExoticCells) {
     //Remove exotic cells before making digits
+    Bool_t exRemoval = fRecoUtils->IsRejectExoticCell();
     fRecoUtils->SwitchOnRejectExoticCell();//switch on and off
     Int_t bunchCrossNo = InputEvent()->GetBunchCrossNumber();
     Bool_t isEx = fRecoUtils->IsExoticCell(cellNumber, fCaloCells, bunchCrossNo);
     if(isEx) accept = kFALSE;
-    fRecoUtils->SwitchOffRejectExoticCell();//switch on and off
+    if(!exRemoval) fRecoUtils->SwitchOffRejectExoticCell();//switch on and off
   }
   return accept;
 }

diff --git a/PWG/EMCAL/AliEmcalClusterMaker.cxx b/PWG/EMCAL/AliEmcalClusterMaker.cxx
index 7f4d5a1a2a2f83ca270a730479411dfea4129296..94c1964e8b079e7ae146bef5639658e0f59c7bfe 100644 (file)
--- a/PWG/EMCAL/AliEmcalClusterMaker.cxx
+++ b/PWG/EMCAL/AliEmcalClusterMaker.cxx
@@ -28,7 +28,8 @@ AliEmcalClusterMaker::AliEmcalClusterMaker() :
   fEnergyTimeHistBefore(0),
   fEnergyDistAfter(0),
   fEtaPhiDistAfter(0),
-  fEnergyTimeHistAfter(0)
+  fEnergyTimeHistAfter(0),
+  fEnergyExoticClusters(0)
 {
   // Default constructor.
 }
@@ -45,7 +46,8 @@ AliEmcalClusterMaker::AliEmcalClusterMaker(const char *name, Bool_t histo) :
   fEnergyTimeHistBefore(0),
   fEnergyDistAfter(0),
   fEtaPhiDistAfter(0),
-  fEnergyTimeHistAfter(0)
+  fEnergyTimeHistAfter(0),
+  fEnergyExoticClusters(0)
 {
   // Standard constructor.
   
@@ -86,6 +88,8 @@ void AliEmcalClusterMaker::UserCreateOutputObjects()
   fOutput->Add(fEtaPhiDistAfter);
   fEnergyTimeHistAfter = new TH2F("hEnergyTimeDistAfter","hEnergyTimeDistAfter;E_{clus} (GeV);time",60,0,30,500,0,1e-6);
   fOutput->Add(fEnergyTimeHistAfter);
+  fEnergyExoticClusters = new TH1F("fEnergyExoticClusters","fEnergyExoticClusters;E_{ex clus} (GeV)",60,0,30);
+  fOutput->Add(fEnergyExoticClusters);
   PostData(1, fOutput);
 }
 
@@ -156,8 +160,14 @@ Bool_t AliEmcalClusterMaker::Run()
     }
     if (fRecoUtils) {
       if (fRecoUtils->IsRejectExoticCluster()) {
-       if (fRecoUtils->IsExoticCluster(oc,fCaloCells))
-       continue;
+       Bool_t exRemoval = fRecoUtils->IsRejectExoticCell();
+       fRecoUtils->SwitchOnRejectExoticCell();                  //switch on temporarily
+       Bool_t exResult = fRecoUtils->IsExoticCluster(oc,fCaloCells);
+       if (!exRemoval) fRecoUtils->SwitchOffRejectExoticCell(); //switch back off
+       if(exResult) {
+         fEnergyExoticClusters->Fill(oc->E());
+         continue;
+       }
       }
       if (fRecoUtils->GetNonLinearityFunction()!=AliEMCALRecoUtils::kNoCorrection) {
        Double_t energy = fRecoUtils->CorrectClusterEnergyLinearity(oc);

diff --git a/PWG/EMCAL/AliEmcalClusterMaker.h b/PWG/EMCAL/AliEmcalClusterMaker.h
index 477fc953c53531b031d9f4ae1747847593c1400d..0f03ec66b512b82309bd34b3ad856a0f4f111195 100644 (file)
--- a/PWG/EMCAL/AliEmcalClusterMaker.h
+++ b/PWG/EMCAL/AliEmcalClusterMaker.h
@@ -35,11 +35,12 @@ class AliEmcalClusterMaker : public AliAnalysisTaskEmcal {
   TH1F                  *fEnergyDistAfter;           //!energy distribution after
   TH2F                  *fEtaPhiDistAfter;           //!eta/phi distribution after
   TH2F                  *fEnergyTimeHistAfter;       //!energy/time distribution after
+  TH1F                  *fEnergyExoticClusters;      //!energy of exotic clusters
 
  private:
   AliEmcalClusterMaker(const AliEmcalClusterMaker&);            // not implemented
   AliEmcalClusterMaker &operator=(const AliEmcalClusterMaker&); // not implemented
 
-  ClassDef(AliEmcalClusterMaker, 1) // Emcal cluster maker
+  ClassDef(AliEmcalClusterMaker, 2) // Emcal cluster maker
 };
 #endif