AliEMCalTriggerPatchAnalysisComponent.cxx

Go to the documentation of this file.

/**************************************************************************
 * Copyright(c) 1998-2014, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/
#include <TClonesArray.h>

#include "AliEmcalTriggerPatchInfo.h"
#include "AliEMCalTriggerBinningComponent.h"
#include "AliEMCalTriggerEventData.h"
#include "AliEMCalTriggerPatchAnalysisComponent.h"

ClassImp(EMCalTriggerPtAnalysis::AliEMCalTriggerPatchAnalysisComponent)

namespace EMCalTriggerPtAnalysis {

AliEMCalTriggerPatchAnalysisComponent::AliEMCalTriggerPatchAnalysisComponent() :
  AliEMCalTriggerTracksAnalysisComponent(),
  fSwapOnlineThresholds(kFALSE),
  fSwapOfflineThresholds(kFALSE)
{
}

AliEMCalTriggerPatchAnalysisComponent::AliEMCalTriggerPatchAnalysisComponent(const char *name) :
  AliEMCalTriggerTracksAnalysisComponent(name),
  fSwapOnlineThresholds(kFALSE),
  fSwapOfflineThresholds(kFALSE)
{
}

void AliEMCalTriggerPatchAnalysisComponent::CreateHistos() {
  AliEMCalTriggerTracksAnalysisComponent::CreateHistos();

  AliEMCalTriggerBinningDimension *etabinning = fBinning->GetBinning("eta"),
      *phibinning = fBinning->GetBinning("phi");
  const TAxis *patchaxes[6] = {
    DefineAxis("energy", 100, 0., 100),
    DefineAxis("energyRough", 100, 0., 100),
    DefineAxis("amplitude", 5000, 0., 5000.),       // limit for the moment
    DefineAxis("eta", etabinning),
    DefineAxis("phi", phibinning),
    DefineAxis("isMain", 2, -0.5, 1.5)
  };

  const TAxis *ampaxes[5] = {
     DefineAxis("amplitudeOnline", 5000, 0., 5000.),
     DefineAxis("amplitudeOffline", 5000, 0., 5000.),
     DefineAxis("eta", etabinning),
     DefineAxis("phi", phibinning),
     DefineAxis("isMain", 2, -0.5, 1.5)
  };

  std::string patchnames[] = {"Level0", "JetHigh", "JetLow", "GammaHigh", "GammaLow"};
  std::string triggermodes[] = {"Online", "Offline"};
  for(std::string * triggerpatch = patchnames; triggerpatch < patchnames + sizeof(patchnames)/sizeof(std::string); ++triggerpatch){
    for(std::string *triggermode = triggermodes; triggermode < triggermodes + sizeof(triggermodes)/sizeof(std::string); ++triggermode){
      if((!strcmp(triggermode->c_str(), "Offline")) && (!strcmp(triggerpatch->c_str(), "Level0"))) continue; // Don't process L0 in case of offline
      printf("Adding patch for trigger %s in case of %s\n", triggerpatch->c_str(), triggermode->c_str());
      fHistos->CreateTHnSparse(Form("PatchInfo%s%s", triggerpatch->c_str(), triggermode->c_str()), Form("Patch energy for %s %s trigger patches", triggerpatch->c_str(), triggermode->c_str()), 6, patchaxes, "s");
      if(strcmp(triggerpatch->c_str(), "Level0")){
        // Add histogram for online-offline of amplitudes
        fHistos->CreateTHnSparse(Form("PatchAmplitude%s%s", triggerpatch->c_str(), triggermode->c_str()), Form("Patch amplitudes for %s %s trigger patches", triggerpatch->c_str(), triggermode->c_str()), 5, ampaxes, "s");
      }
    }
  }
}

void AliEMCalTriggerPatchAnalysisComponent::Process(const AliEMCalTriggerEventData* const data) {
  AliEmcalTriggerPatchInfo *triggerpatch(NULL);
  TIter patchIter(data->GetTriggerPatchContainer());
  while((triggerpatch = dynamic_cast<AliEmcalTriggerPatchInfo *>(patchIter()))){
  bool isMain = triggerpatch->IsOfflineSimple() ? triggerpatch->IsMainTriggerSimple() : triggerpatch->IsMainTrigger();
    double triggerpatchinfo[6] = {triggerpatch->GetPatchE(),triggerpatch->GetADCAmpGeVRough(),
        static_cast<double>(triggerpatch->IsOfflineSimple() ? triggerpatch->GetADCOfflineAmp() : triggerpatch->GetADCAmp()), triggerpatch->GetEtaGeo(),
        triggerpatch->GetPhiGeo(), isMain ? 1. : 0.};
    double amplitudeinfo[5] = {triggerpatch->GetADCAmp(), triggerpatch->GetADCOfflineAmp(), triggerpatch->GetEtaGeo(), triggerpatch->GetPhiGeo(), isMain ? 1. : 0.};
    if(triggerpatch->IsOfflineSimple()){
      if((!fSwapOfflineThresholds && triggerpatch->IsJetHighSimple()) || (fSwapOfflineThresholds && triggerpatch->IsJetLowSimple())){
        fHistos->FillTHnSparse("PatchInfoJetHighOffline", triggerpatchinfo);
        fHistos->FillTHnSparse("PatchAmplitudeJetHighOffline", amplitudeinfo);
      }
      if((!fSwapOfflineThresholds && triggerpatch->IsJetLowSimple()) || (fSwapOfflineThresholds && triggerpatch->IsJetHighSimple())){
        fHistos->FillTHnSparse("PatchInfoJetLowOffline", triggerpatchinfo);
        fHistos->FillTHnSparse("PatchAmplitudeJetLowOffline", amplitudeinfo);
      }
      if((!fSwapOfflineThresholds && triggerpatch->IsGammaHighSimple()) || (fSwapOfflineThresholds && triggerpatch->IsGammaLowSimple())){
        fHistos->FillTHnSparse("PatchInfoGammaHighOffline", triggerpatchinfo);
        fHistos->FillTHnSparse("PatchAmplitudeGammaHighOffline", amplitudeinfo);
      }
      if((!fSwapOfflineThresholds && triggerpatch->IsGammaLowSimple()) || (fSwapOfflineThresholds && triggerpatch->IsGammaHighSimple())){
        fHistos->FillTHnSparse("PatchInfoGammaLowOffline", triggerpatchinfo);
        fHistos->FillTHnSparse("PatchAmplitudeGammaLowOffline", amplitudeinfo);
      }
    } else{
      if((!fSwapOnlineThresholds && triggerpatch->IsJetHigh()) || (fSwapOnlineThresholds && triggerpatch->IsJetLow())){
        fHistos->FillTHnSparse("PatchInfoJetHighOnline", triggerpatchinfo);
        fHistos->FillTHnSparse("PatchAmplitudeJetHighOnline", amplitudeinfo);
      }
      if((!fSwapOnlineThresholds && triggerpatch->IsJetLow()) || (fSwapOnlineThresholds && triggerpatch->IsJetHigh())){
        fHistos->FillTHnSparse("PatchInfoJetLowOnline", triggerpatchinfo);
        fHistos->FillTHnSparse("PatchAmplitudeJetLowOnline", amplitudeinfo);
      }
      if((!fSwapOnlineThresholds && triggerpatch->IsGammaHigh()) || (fSwapOnlineThresholds && triggerpatch->IsGammaLow())){
        fHistos->FillTHnSparse("PatchInfoGammaHighOnline", triggerpatchinfo);
        fHistos->FillTHnSparse("PatchAmplitudeGammaHighOnline", amplitudeinfo);
      }
      if((!fSwapOnlineThresholds && triggerpatch->IsGammaLow()) || (fSwapOnlineThresholds && triggerpatch->IsGammaHigh())){
        fHistos->FillTHnSparse("PatchInfoGammaLowOnline", triggerpatchinfo);
        fHistos->FillTHnSparse("PatchAmplitudeGammaLowOnline", amplitudeinfo);
      }
      if(triggerpatch->IsLevel0()){
        fHistos->FillTHnSparse("PatchInfoLevel0Online", triggerpatchinfo);
      }
    }
  }
}

} /* namespace EMCalTriggerPtAnalysis */