Fix for Bug: 58726 including a new QA histogram for ACORDE ESDs (Multiplicity of...

[u/mrichter/AliRoot.git] / ACORDE / AliACORDEReconstructor.cxx

/**************************************************************************\r
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *\r
 *                                                                        *\r
 * Author: The ALICE Off-line Project.                                    *\r
 * Contributors are mentioned in the code where appropriate.              *\r
 *                                                                        *\r
 * Permission to use, copy, modify and distribute this software and its   *\r
 * documentation strictly for non-commercial purposes is hereby granted   *\r
 * without fee, provided that the above copyright notice appears in all   *\r
 * copies and that both the copyright notice and this permission notice   *\r
 * appear in the supporting documentation. The authors make no claims     *\r
 * about the suitability of this software for any purpose. It is          *\r
 * provided "as is" without express or implied warranty.                  *\r
 **************************************************************************/\r
\r
/* $Id: AliACORDEReconstructor.cxx 20956 2007-09-26 14:22:18Z mrodrigu $ */\r
//////////////////////////////////////////////////////////////////////////////\r
//                                                                          //\r
//  Class for ACORDE reconstruction                                         //\r
//////////////////////////////////////////////////////////////////////////////\r
\r
#include "AliRawReader.h"\r
\r
#include "AliACORDEReconstructor.h"\r
#include "AliACORDERawStream.h"\r
#include "AliESDEvent.h"\r
#include "AliACORDEdigit.h"\r
#include "AliACORDERecoParam.h"\r
\r
ClassImp(AliACORDEReconstructor)\r
\r
AliACORDEReconstructor:: AliACORDEReconstructor():\r
  AliReconstructor(),\r
  fESDACORDE(0x0),\r
  fAcordeRecoParam(0x0),\r
  fCalibData(0x0)\r
{\r
  // Default constructor  \r
  // Get calibration data\r
\r
  fCalibData = GetCalibData();\r
  fAcordeRecoParam = GetRecoParam();\r
}\r
\r
//_______________________________________________________________________\r
AliACORDECalibData *AliACORDEReconstructor::GetCalibData() const\r
{\r
  return 0x0;\r
}\r
//____________________________________________________________________________\r
AliACORDERecoParam *AliACORDEReconstructor::GetRecoParam() const\r
{\r
  return 0x0;\r
}\r
//_____________________________________________________________________________\r
AliACORDEReconstructor& AliACORDEReconstructor::operator = \r
  (const AliACORDEReconstructor& /*reconstructor*/)\r
{\r
// assignment operator\r
\r
  Fatal("operator =", "assignment operator not implemented");\r
  return *this;\r
}\r
\r
//_____________________________________________________________________________\r
AliACORDEReconstructor::~AliACORDEReconstructor()\r
{\r
// destructor\r
  delete fESDACORDE; \r
}\r
\r
//_____________________________________________________________________________\r
void AliACORDEReconstructor::Init()\r
{\r
// initializer\r
    fESDACORDE  = new AliESDACORDE;\r
}\r
\r
void AliACORDEReconstructor::ConvertDigits(AliRawReader* rawReader, TTree* digitsTree) const\r
{\r
\r
  if (!digitsTree) {\r
    AliError("No digits tree!");\r
    return;\r
  }\r
\r
  TClonesArray* digitsArray = new TClonesArray("AliACORDEdigit");\r
  digitsTree->Branch("ACORDEdigit", &digitsArray);\r
\r
  rawReader->Reset();\r
  AliACORDERawStream rawStream(rawReader);\r
  if (rawStream.Next()) {\r
    for(Int_t iChannel = 0; iChannel < 60; iChannel++) {\r
      Int_t  index = iChannel / 30;\r
      Int_t  bit   = iChannel % 30;\r
      if (rawStream.GetWord(index) & (1 << bit))\r
        new ((*digitsArray)[digitsArray->GetEntriesFast()]) AliACORDEdigit(iChannel+1,0);\r
    }\r
  }\r
\r
  digitsTree->Fill();\r
			\r
}\r
\r
void AliACORDEReconstructor::FillESD(TTree* digitsTree, TTree* /*clustersTree*/,AliESDEvent* esd) const\r
{\r
\r
  // fills ESD with ACORDE Digits\r
\r
  if (!digitsTree)\r
    {\r
      AliError("No digits tree!");\r
      return;\r
    }\r
\r
  TClonesArray* digitsArray = NULL;\r
  TBranch* digitBranch = digitsTree->GetBranch("ACORDEdigit");\r
  if (!digitBranch) {\r
    AliError("No ACORDE digits branch found!");\r
    return;\r
  }\r
  digitBranch->SetAddress(&digitsArray);\r
\r
  digitsTree->GetEvent(0);\r
\r
  Bool_t AcoHitSingle[60],AcoHitMulti[60];\r
  for(Int_t i = 0; i < 60; i++) { AcoHitSingle[i] = AcoHitMulti[i] = kFALSE; }\r
\r
  Int_t nDigits = digitsArray->GetEntriesFast();\r
    \r
  for (Int_t d=0; d<nDigits; d++) {    \r
    AliACORDEdigit* digit = (AliACORDEdigit*)digitsArray->At(d);\r
    Int_t module = digit->GetModule();\r
\r
    AcoHitSingle[module-1] = kTRUE;\r
    AcoHitMulti[module-1] = kTRUE;\r
  }  \r
  TString ActiveTriggerDetector = esd->GetFiredTriggerClasses();\r
  if (ActiveTriggerDetector.Contains("ASL")) fESDACORDE->SetACORDEBitPattern(AcoHitSingle);\r
  else if (ActiveTriggerDetector.Contains("AMU")) fESDACORDE->SetACORDEBitPattern(AcoHitMulti);\r
	else fESDACORDE->SetACORDEBitPattern(AcoHitSingle);\r
\r
  if (esd)\r
    {\r
      AliDebug(1, Form("Writing ACORDE data to ESD Tree"));\r
      esd->SetACORDEData(fESDACORDE);\r
    }	\r
}\r
\r
\r