HLT/JET/AliHLTJETReaderHeader.cxx

   1 //-*- Mode: C++ -*-
   2 // $Id: AliHLTJETReaderHeader.cxx  $
   3 /**************************************************************************
   4  * This file is property of and copyright by the ALICE HLT Project        *
   5  * ALICE Experiment at CERN, All rights reserved.                         *
   6  *                                                                        *
   7  * Primary Authors: Jochen Thaeder <thaeder@kip.uni-heidelberg.de>        *
   8  *                  for The ALICE HLT Project.                            *
   9  *                                                                        *
  10  * Permission to use, copy, modify and distribute this software and its   *
  11  * documentation strictly for non-commercial purposes is hereby granted   *
  12  * without fee, provided that the above copyright notice appears in all   *
  13  * copies and that both the copyright notice and this permission notice   *
  14  * appear in the supporting documentation. The authors make no claims     *
  15  * about the suitability of this software for any purpose. It is          *
  16  * provided "as is" without express or implied warranty.                  *
  17  **************************************************************************/
  18
  19 /** @file   AliHLTJETReaderHeader.cxx
  20     @author Jochen Thaeder
  21     @date
  22     @brief   ReaderHeader for jet finder
  23 */
  24
  25 // see header file for class documentation
  26 // or
  27 // refer to README to build package
  28 // or
  29 // visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
  30
  31 #include "AliHLTJETReaderHeader.h"
  32 #include "AliHLTJETTrackCuts.h"
  33
  34 using namespace std;
  35
  36 /** ROOT macro for the implementation of ROOT specific class methods */
  37 ClassImp(AliHLTJETReaderHeader)
  38
  39 /*
  40  * ---------------------------------------------------------------------------------
  41  *                            Constructor / Destructor
  42  * ---------------------------------------------------------------------------------
  43  */
  44
  45 // #################################################################################
  46 AliHLTJETReaderHeader::AliHLTJETReaderHeader()
  47   :
  48   AliJetReaderHeader("AliHLTJETReaderHeader"),
  49   fTrackCuts(NULL),
  50   fSeedCuts(NULL),
  51   fGridEtaBinning(0.0),
  52   fGridPhiBinning(0.0),
  53   fGridEtaRange(0.0),
  54   fGridPhiRange(0.0),
  55   fAlgorithm(AliHLTJETBase::kFFSCSquareCell),
  56   fConeRadius(0.0),
  57   fUseMC(kFALSE) {
  58   // see header file for class documentation
  59   // or
  60   // refer to README to build package
  61   // or
  62   // visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
  63
  64 }
  65
  66 // #################################################################################
  67 AliHLTJETReaderHeader::~AliHLTJETReaderHeader() {
  68   // see header file for class documentation
  69
  70 }
  71
  72 /*
  73  * ---------------------------------------------------------------------------------
  74  *                                   Initialize
  75  * ---------------------------------------------------------------------------------
  76  */
  77
  78 // #################################################################################
  79 Int_t AliHLTJETReaderHeader::Initialize() {
  80   // see header file for class documentation
  81
  82   Int_t iResult = 0;
  83
  84   // -- Set eta and phi range for grid
  85   fGridPhiRange = fFiducialPhiMin +
  86     fFiducialPhiMax + ( 2.0 * fConeRadius );
  87
  88   fGridEtaRange = TMath::Abs( fFiducialEtaMin ) + fFiducialEtaMax;
  89
  90   HLTInfo(" -= ReaderHeader =- ");
  91   HLTInfo(" Cone radius      %f", fConeRadius );
  92   HLTInfo(" Grid eta binning %f", fGridEtaBinning );
  93   HLTInfo(" Grid phi binning %f", fGridPhiBinning );
  94   HLTInfo(" Grid eta range   %f", fGridEtaRange );
  95   HLTInfo(" Grid phi range   %f", fGridPhiRange );
  96   HLTInfo(" Algorithm        %s", AliHLTJETBase::fgkJetAlgorithmType[fAlgorithm] );
  97
  98   if (fUseMC) { HLTInfo(" Use Kinematics   TRUE"); }
  99   else { HLTInfo( " Use Kinematics   FALSE"); }
 100
 101   if ( ! fTrackCuts ) {
 102     HLTError("No track cuts set in reader header");
 103     iResult = -EINPROGRESS;
 104   }
 105   else {
 106     fTrackCuts->SetEtaRange( fFiducialEtaMin, fFiducialEtaMax );
 107     fTrackCuts->SetPhiRange( fFiducialPhiMin, fFiducialPhiMax );
 108     HLTInfo(" -= TrackCuts =- " );
 109   }
 110
 111   // Check for seeds only cone based algoritms
 112   if ( fAlgorithm >= AliHLTJETBase::kFFSCSquareCell ) {
 113     if ( !fSeedCuts ) {
 114       HLTError("No seed cuts set in reader header");
 115       iResult = -EINPROGRESS;
 116     }
 117     else {
 118       fSeedCuts->SetEtaRange( fFiducialEtaMin+fConeRadius,
 119                               fFiducialEtaMax-fConeRadius );
 120       fSeedCuts->SetPhiRange( fFiducialPhiMin, fFiducialPhiMax );
 121       HLTInfo(" -= SeedCuts =- " );
 122     }
 123   }
 124
 125   return iResult;
 126 }
 127