StCentralityMaker.cxx

/******************************************************************************
 * $Id: StCentralityMaker.cxx,v 1.4 2012/04/25 05:17:24 hmasui Exp $
 * $Log: StCentralityMaker.cxx,v $
 * Revision 1.4  2012/04/25 05:17:24  hmasui
 * Remove SetModeNBD() and corresponding data member, merged into GetNegativeBinomial() function
 *
******************************************************************************/

#include <assert.h>

#include "TError.h"

#include "StMessMgr.h"
#include "StNegativeBinomial.h"
#include "StCentrality.h"
#include "StCentralityMaker.h"

ClassImp(StCentralityMaker)

//____________________________________________________________________________________________________
// Default constructor
StCentralityMaker::StCentralityMaker()
  : mNBinomial(0), mCentrality(0)
{
  LOG_INFO << "StCentralityMaker  User needs to call Init(const Char_t* system) function to Initialize NBD (see below)" << endm;
  StCentrality::help();
}

//____________________________________________________________________________________________________
// Initialization inside the constructor
StCentralityMaker::StCentralityMaker(const Char_t* system)
  : mNBinomial(0), mCentrality(0)
{
  Init(system);
}

//____________________________________________________________________________________________________
// Default destructor
StCentralityMaker::~StCentralityMaker()
{
}

//____________________________________________________________________________________________________
const StNegativeBinomial* StCentralityMaker::GetNegativeBinomial(const UInt_t id) const
{
  if( id >= mNBinomial.size() ){
    Error("StCentralityMaker::GetNegativeBinomial", "Invalid index for negative binomial, id=%3d", id);
    assert(0);
  }

  return mNBinomial[id] ;
}

//____________________________________________________________________________________________________
StCentrality* StCentralityMaker::GetCentrality(const UInt_t id) const
{
  if( id >= mCentrality.size() ){
    Error("StCentralityMaker::GetCentrality", "Invalid index for centrality, id=%3d", id);
    assert(0);
  }

  return mCentrality[id] ;
}

//____________________________________________________________________________________________________
// Use this function to initialize centrality bin and NBD
void StCentralityMaker::Init(const Char_t* system)
{
  const TString systemName(system);

  LOG_INFO << "StCentralityMaker::Init  Initialization for " << systemName.Data() << endm;

  const TString type[] = {"default", "low", "high"};
  for(UInt_t id=0; id<3; id++){
    mCentrality.push_back( new StCentrality(systemName, type[id]) );
  }

  // Use multiplicity depdent efficiency for all systems
  for(UInt_t id=0; id<3; id++){
    Bool_t isConstEfficiency = kFALSE ;

    mNBinomial.push_back( 
        new StNegativeBinomial( mCentrality[id]->GetNpp(), mCentrality[id]->GetK(), mCentrality[id]->GetX(),
        mCentrality[id]->GetEfficiency(), mCentrality[id]->GetTriggerBias(), isConstEfficiency)
        );
  }
}