add qa check
authormverweij <marta.verweij@cern.ch>
Fri, 28 Mar 2014 13:12:35 +0000 (14:12 +0100)
committermverweij <marta.verweij@cern.ch>
Fri, 28 Mar 2014 13:13:48 +0000 (14:13 +0100)
PWG/EMCAL/AliAnalysisTaskEMCALClusterizeFast.cxx
PWG/EMCAL/AliEmcalClusterMaker.cxx
PWG/EMCAL/AliEmcalClusterMaker.h

index 7dfd455..1e30462 100644 (file)
@@ -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;
 }
index 7f4d5a1..94c1964 100644 (file)
@@ -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);
index 477fc95..0f03ec6 100644 (file)
@@ -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