PndRichCalDb Class Reference

#include <PndRichCalDb.h>

Inheritance diagram for PndRichCalDb:

Public Member Functions
	PndRichCalDb ()
	PndRichCalDb (const char *mapName, const char *fileType="R")
	PndRichCalDb (PndRichCalDbPar *CaldbPar)
virtual	~PndRichCalDb ()
virtual void	Init ()
virtual Double_t	GetBetaMean (dbpoint pnt)
virtual Double_t	GetBetaSig (dbpoint pnt)
virtual Double_t	GetBetaEff (dbpoint pnt)
virtual Bool_t	IsInside (Double_t p, Double_t x, Double_t y, Double_t t, Double_t f, Int_t &ip, Int_t &ix, Int_t &iy, Int_t &it, Int_t &iq, Double_t &dp, Double_t &dx, Double_t &dy, Double_t &dt, Double_t &df)
void	WriteAsciiFile (const char *fileName)
void	WriteRootFile (const char *fileName, const char *mapName)
Double_t	GetPmin () const
Double_t	GetXmin () const
Double_t	GetYmin () const
Double_t	GetTmin () const
Double_t	GetFmin () const
Double_t	GetPmax () const
Double_t	GetXmax () const
Double_t	GetYmax () const
Double_t	GetTmax () const
Double_t	GetFmax () const
Double_t	GetPstep () const
Double_t	GetXstep () const
Double_t	GetYstep () const
Double_t	GetTstep () const
Double_t	GetFstep () const
Int_t	GetNp () const
Int_t	GetNx () const
Int_t	GetNy () const
Int_t	GetNt () const
Int_t	GetNf () const
TArrayF *	GetBetaMean () const
TArrayF *	GetBetaSig () const
TArrayF *	GetBetaEff () const
const char *	GetFileName ()
virtual void	Print ()

Protected Member Functions
void	Reset ()
void	ReadAsciiFile (const char *fileName)
void	ReadRootFile (const char *fileName, const char *mapName)
void	SetCalDb (const PndRichCalDbData *data)

Protected Attributes
TString	fFileName
Double_t	fPmin
Double_t	fPmax
Double_t	fPstep
Double_t	fXmin
Double_t	fXmax
Double_t	fXstep
Double_t	fYmin
Double_t	fYmax
Double_t	fYstep
Double_t	fTmin
Double_t	fTmax
Double_t	fTstep
Double_t	fFmin
Double_t	fFmax
Double_t	fFstep
Int_t	fNp
Int_t	fNx
Int_t	fNy
Int_t	fNt
Int_t	fNf
TArrayF *	fBetaMean
TArrayF *	fBetaSig
TArrayF *	fBetaEff

Private Member Functions
	PndRichCalDb (const PndRichCalDb &L)
PndRichCalDb &	operator= (const PndRichCalDb &)

Detailed Description

Definition at line 31 of file PndRichCalDb.h.

Constructor & Destructor Documentation

PndRichCalDb::PndRichCalDb

(

)

Default constructor

Definition at line 24 of file PndRichCalDb.cxx.

References fBetaEff, fBetaMean, fBetaSig, fFmax, fFmin, fFstep, fNf, fNp, fNt, fNx, fNy, fPmax, fPmin, fPstep, fTmax, fTmin, fTstep, fXmax, fXmin, fXstep, fYmax, fYmin, fYstep, and SetName.

  : FairField(),
    fFileName(""),
    fPmin(0), fPmax(0), fPstep(0),
    fXmin(0), fXmax(0), fXstep(0),
    fYmin(0), fYmax(0), fYstep(0),
    fTmin(0), fTmax(0), fTstep(0),
    fFmin(0), fFmax(0), fFstep(0),
    fNp(0), fNx(0), fNy(0), fNt(0), fNf(0),
    fBetaMean(NULL), fBetaSig(NULL), fBetaEff(NULL)
{
  SetName("");
  fType = 0;
  fPmin=1; fPmax=11; fPstep=10;
  fXmin=2; fXmax=12; fXstep=10;
  fYmin=3; fYmax=13; fYstep=10;
  fTmin=4; fTmax=14; fTstep=10;
  fFmin=5; fFmax=15; fFstep=10;
  fNp=2; fNx=2; fNy=2; fNt=2; fNf=2;
  Int_t size = fNp*fNx*fNy*fNt*fNf;
  fBetaSig = new TArrayF(size);
  fBetaMean = new TArrayF(size);
  fBetaEff = new TArrayF(size);

}

PndRichCalDb::PndRichCalDb	(	const char *	mapName,
		const char *	fileType = "R"
	)

Standard constructor

Parameters

name	Name of field map
fileType	R = ROOT file, A = ASCII

Definition at line 54 of file PndRichCalDb.cxx.

References fFileName, SetName, and TString.

  : FairField(mapName),
    fFileName(""),
    fPmin(0), fPmax(0), fPstep(0),
    fXmin(0), fXmax(0), fXstep(0),
    fYmin(0), fYmax(0), fYstep(0),
    fTmin(0), fTmax(0), fTstep(0),
    fFmin(0), fFmax(0), fFstep(0),
    fNp(0), fNx(0), fNy(0), fNt(0), fNf(0),
    fBetaMean(NULL), fBetaSig(NULL), fBetaEff(NULL)
{
  SetName(mapName);
  TString dir = getenv("VMCWORKDIR");
  fFileName = dir + "/detectors/rich/" + mapName;
  if ( fileType[0] == 'R' ) fFileName += ".root";
  else                      fFileName += ".dat";
  fType = 1;
}

PndRichCalDb::PndRichCalDb

(

PndRichCalDbPar *

CaldbPar

)

Constructor from PndRichCalDbPar

Definition at line 77 of file PndRichCalDb.cxx.

References fFileName, PndRichCalDbPar::GetType(), PndRichCalDbPar::MapName(), SetName, and TString.

  : FairField(),
    fFileName(""),
    fPmin(0), fPmax(0), fPstep(0),
    fXmin(0), fXmax(0), fXstep(0),
    fYmin(0), fYmax(0), fYstep(0),
    fTmin(0), fTmax(0), fTstep(0),
    fFmin(0), fFmax(0), fFstep(0),
    fNp(0), fNx(0), fNy(0), fNt(0), fNf(0),
    fBetaMean(NULL), fBetaSig(NULL), fBetaEff(NULL)
{
  fType = 1;
  if ( ! calDbPar ) {
     LOG(ERROR) << "PndRichCalDb::PndRichCalDb: empty parameter container!";
     SetName("");
     fType     = -1;
  }
  else {
    TString Name=GetName();
    calDbPar->MapName(Name);
    TString dir = getenv("VMCWORKDIR");
    fFileName = dir + "/detectors/rich/" + Name + ".root";
    fFileName = Name + ".root";
    fType = calDbPar->GetType();
  }
}

PndRichCalDb::~PndRichCalDb ( )

virtual

Destructor

Definition at line 124 of file PndRichCalDb.cxx.

References fBetaEff, fBetaMean, and fBetaSig.

                            {
   if ( fBetaMean ) delete fBetaMean;
   if ( fBetaSig )  delete fBetaSig;
   if ( fBetaEff )  delete fBetaEff;
}

PndRichCalDb::PndRichCalDb ( const PndRichCalDb &  L )

private

Definition at line 106 of file PndRichCalDb.cxx.

  : FairField(),
    fFileName(L.fFileName),
    fPmin(L.fPmin), fPmax(L.fPmax), fPstep(L.fPstep),
    fXmin(L.fXmin), fXmax(L.fXmax), fXstep(L.fXstep),
    fYmin(L.fYmin), fYmax(L.fYmax), fYstep(L.fYstep),
    fTmin(L.fTmin), fTmax(L.fTmax), fTstep(L.fTstep),
    fFmin(L.fFmin), fFmax(L.fFmax), fFstep(L.fFstep),
    fNp(L.fNp), fNx(L.fNx), fNy(L.fNy), fNt(L.fNt), fNf(L.fNf),
    fBetaMean(L.fBetaMean), fBetaSig(L.fBetaSig), fBetaEff(L.fBetaEff)
{
  fType = L.fType;
}

Member Function Documentation

Double_t PndRichCalDb::GetBetaEff ( dbpoint  pnt )

virtual

Definition at line 276 of file PndRichCalDb.cxx.

References dbpoint::beta, bp, Double_t, dx, dy, f, fBetaEff, fNf, fNt, fNx, fNy, IsInside(), p, dbpoint::phi, t, dbpoint::theta, dbpoint::x, x, dbpoint::y, and y.

Referenced by PndRichResolution::Efficiency().

                                             {

  Double_t p = pnt.beta;
  Double_t x = pnt.x;
  Double_t y = pnt.y;
  Double_t t = pnt.theta*180/3.1415927;
  Double_t f = pnt.phi*180/3.1415927;
  Int_t ip    = 0;
  Int_t ix    = 0;
  Int_t iy    = 0;
  Int_t it    = 0;
  Int_t iq    = 0;
  Double_t dp = 0.;
  Double_t dx = 0.;
  Double_t dy = 0.;
  Double_t dt = 0.;
  Double_t df = 0.;

  if (x>=0&&y<0) { f = 360 - f; y = -y; }
  if (x<0&&y>=0) { f = 180 - f; x = -x; }
  if (x<0&&y<0)  { f = 180 + f; x = -x; y = -y; }
  if (f>360) f -= 360;
  if (f<0) f +=360;

  if ( IsInside(p, x, y, t, f, ip, ix, iy, it, iq, dp, dx, dy, dt, df) ) {

     Double_t pq[5][2];
     pq[0][1] = dp; pq[0][0] = 1 - pq[0][1];
     pq[1][1] = dx; pq[1][0] = 1 - pq[1][1];
     pq[2][1] = dy; pq[2][0] = 1 - pq[2][1];
     pq[3][1] = dt; pq[3][0] = 1 - pq[3][1];
     pq[4][1] = df; pq[4][0] = 1 - pq[4][1];

     Double_t value = 0;
     Int_t index;
     for(Int_t bp=0; bp<2; bp++ ) {
        for(Int_t bx=0; bx<2; bx++ ) {
           for(Int_t by=0; by<2; by++ ) {
              for(Int_t bt=0; bt<2; bt++ ) {
                 for(Int_t bf=0; bf<2; bf++ ) {
                    index =
                       (ip+bp)*fNx*fNy*fNt*fNf +
                       (ix+bx)*fNy*fNt*fNf +
                       (iy+by)*fNt*fNf +
                       (it+bt)*fNf +
                       (iq+bf);
                    value += fBetaEff->At(index)*pq[0][bp]*pq[1][bx]*pq[2][by]*pq[3][bt]*pq[4][bf];
                 }
              }
           }
        }
     }

     // Return interpolated field value
     return value;

  }

  return 0.;
}

TArrayF* PndRichCalDb::GetBetaEff ( ) const

inline

Definition at line 113 of file PndRichCalDb.h.

References fBetaEff.

{ return fBetaEff; }

Double_t PndRichCalDb::GetBetaMean ( dbpoint  pnt )

virtual

Get the field components at a certain point

Parameters

x,y,z

Point coordinates (global) [cm] Bx,By,Bz Field components [kG]

Definition at line 147 of file PndRichCalDb.cxx.

References dbpoint::beta, bp, Double_t, dx, dy, f, fBetaMean, fNf, fNt, fNx, fNy, IsInside(), p, dbpoint::phi, t, dbpoint::theta, dbpoint::x, x, dbpoint::y, and y.

Referenced by PndRichResolution::Shift().

                                              {

  Double_t p = pnt.beta;
  Double_t x = pnt.x;
  Double_t y = pnt.y;
  Double_t t = pnt.theta*180/3.1415927;
  Double_t f = pnt.phi*180/3.1415927;
  Int_t ip    = 0;
  Int_t ix    = 0;
  Int_t iy    = 0;
  Int_t it    = 0;
  Int_t iq    = 0;
  Double_t dp = 0.;
  Double_t dx = 0.;
  Double_t dy = 0.;
  Double_t dt = 0.;
  Double_t df = 0.;

  if (x>=0&&y<0) { f = 360 - f; y = -y; }
  if (x<0&&y>=0) { f = 180 - f; x = -x; }
  if (x<0&&y<0)  { f = 180 + f; x = -x; y = -y; }
  if (f>360) f -= 360;
  if (f<0) f +=360;

  if ( IsInside(p, x, y, t, f, ip, ix, iy, it, iq, dp, dx, dy, dt, df) ) {

     Double_t pq[5][2];
     pq[0][1] = dp; pq[0][0] = 1 - pq[0][1];
     pq[1][1] = dx; pq[1][0] = 1 - pq[1][1];
     pq[2][1] = dy; pq[2][0] = 1 - pq[2][1];
     pq[3][1] = dt; pq[3][0] = 1 - pq[3][1];
     pq[4][1] = df; pq[4][0] = 1 - pq[4][1];

     Double_t value = 0;
     Int_t index;
     for(Int_t bp=0; bp<2; bp++ ) {
        for(Int_t bx=0; bx<2; bx++ ) {
           for(Int_t by=0; by<2; by++ ) {
              for(Int_t bt=0; bt<2; bt++ ) {
                 for(Int_t bf=0; bf<2; bf++ ) {
                    index =
                       (ip+bp)*fNx*fNy*fNt*fNf +
                       (ix+bx)*fNy*fNt*fNf +
                       (iy+by)*fNt*fNf +
                       (it+bt)*fNf +
                       (iq+bf);
                    value += fBetaMean->At(index)*pq[0][bp]*pq[1][bx]*pq[2][by]*pq[3][bt]*pq[4][bf];
                 }
              }
           }
        }
     }

     // Return interpolated field value
     return value;

  }

  return 0.;
}

TArrayF* PndRichCalDb::GetBetaMean ( ) const

inline

Accessors to the field value arrays

Definition at line 111 of file PndRichCalDb.h.

References fBetaMean.

{ return fBetaMean; }

Double_t PndRichCalDb::GetBetaSig ( dbpoint  pnt )

virtual

Definition at line 212 of file PndRichCalDb.cxx.

References dbpoint::beta, bp, Double_t, dx, dy, f, fBetaSig, fNf, fNt, fNx, fNy, IsInside(), p, dbpoint::phi, t, dbpoint::theta, dbpoint::x, x, dbpoint::y, and y.

Referenced by PndRichResolution::Sigma().

                                             {

  Double_t p = pnt.beta;
  Double_t x = pnt.x;
  Double_t y = pnt.y;
  Double_t t = pnt.theta*180/3.1415927;
  Double_t f = pnt.phi*180/3.1415927;
  Int_t ip    = 0;
  Int_t ix    = 0;
  Int_t iy    = 0;
  Int_t it    = 0;
  Int_t iq    = 0;
  Double_t dp = 0.;
  Double_t dx = 0.;
  Double_t dy = 0.;
  Double_t dt = 0.;
  Double_t df = 0.;

  if (x>=0&&y<0) { f = 360 - f; y = -y; }
  if (x<0&&y>=0) { f = 180 - f; x = -x; }
  if (x<0&&y<0)  { f = 180 + f; x = -x; y = -y; }
  if (f>360) f -= 360;
  if (f<0) f +=360;

  if ( IsInside(p, x, y, t, f, ip, ix, iy, it, iq, dp, dx, dy, dt, df) ) {

     Double_t pq[5][2];
     pq[0][1] = dp; pq[0][0] = 1 - pq[0][1];
     pq[1][1] = dx; pq[1][0] = 1 - pq[1][1];
     pq[2][1] = dy; pq[2][0] = 1 - pq[2][1];
     pq[3][1] = dt; pq[3][0] = 1 - pq[3][1];
     pq[4][1] = df; pq[4][0] = 1 - pq[4][1];

     Double_t value = 0;
     Int_t index;
     for(Int_t bp=0; bp<2; bp++ ) {
        for(Int_t bx=0; bx<2; bx++ ) {
           for(Int_t by=0; by<2; by++ ) {
              for(Int_t bt=0; bt<2; bt++ ) {
                 for(Int_t bf=0; bf<2; bf++ ) {
                    index =
                       (ip+bp)*fNx*fNy*fNt*fNf +
                       (ix+bx)*fNy*fNt*fNf +
                       (iy+by)*fNt*fNf +
                       (it+bt)*fNf +
                       (iq+bf);
                    value += fBetaSig->At(index)*pq[0][bp]*pq[1][bx]*pq[2][by]*pq[3][bt]*pq[4][bf];
                 }
              }
           }
        }
     }

     // Return interpolated field value
     return value;

  }

  return 0.;
}

TArrayF* PndRichCalDb::GetBetaSig ( ) const

inline

Definition at line 112 of file PndRichCalDb.h.

References fBetaSig.

{ return fBetaSig; }

const char* PndRichCalDb::GetFileName ( )

inline

Accessor to field map file

Definition at line 117 of file PndRichCalDb.h.

References fFileName.

{ return fFileName.Data(); }

Double_t PndRichCalDb::GetFmax ( ) const

inline

Definition at line 97 of file PndRichCalDb.h.

References fFmax.

{ return fFmax; }

Double_t PndRichCalDb::GetFmin ( ) const

inline

Definition at line 92 of file PndRichCalDb.h.

References fFmin.

{ return fFmin; }

Double_t PndRichCalDb::GetFstep ( ) const

inline

Definition at line 102 of file PndRichCalDb.h.

References fFstep.

{ return fFstep; }

Int_t PndRichCalDb::GetNf ( ) const

inline

Definition at line 107 of file PndRichCalDb.h.

References fNf.

{ return fNf; }

Int_t PndRichCalDb::GetNp ( ) const

inline

Definition at line 103 of file PndRichCalDb.h.

References fNp.

{ return fNp; }

Int_t PndRichCalDb::GetNt ( ) const

inline

Definition at line 106 of file PndRichCalDb.h.

References fNt.

{ return fNt; }

Int_t PndRichCalDb::GetNx ( ) const

inline

Definition at line 104 of file PndRichCalDb.h.

References fNx.

{ return fNx; }

Int_t PndRichCalDb::GetNy ( ) const

inline

Definition at line 105 of file PndRichCalDb.h.

References fNy.

{ return fNy; }

Double_t PndRichCalDb::GetPmax ( ) const

inline

Definition at line 93 of file PndRichCalDb.h.

References fPmax.

{ return fPmax; } 

Double_t PndRichCalDb::GetPmin ( ) const

inline

Accessors to field parameters in local coordinate system

Definition at line 88 of file PndRichCalDb.h.

References fPmin.

{ return fPmin; } 

Double_t PndRichCalDb::GetPstep ( ) const

inline

Definition at line 98 of file PndRichCalDb.h.

References fPstep.

{ return fPstep; } 

Double_t PndRichCalDb::GetTmax ( ) const

inline

Definition at line 96 of file PndRichCalDb.h.

References fTmax.

{ return fTmax; } 

Double_t PndRichCalDb::GetTmin ( ) const

inline

Definition at line 91 of file PndRichCalDb.h.

References fTmin.

{ return fTmin; } 

Double_t PndRichCalDb::GetTstep ( ) const

inline

Definition at line 101 of file PndRichCalDb.h.

References fTstep.

{ return fTstep; } 

Double_t PndRichCalDb::GetXmax ( ) const

inline

Definition at line 94 of file PndRichCalDb.h.

References fXmax.

{ return fXmax; } 

Double_t PndRichCalDb::GetXmin ( ) const

inline

Definition at line 89 of file PndRichCalDb.h.

References fXmin.

{ return fXmin; } 

Double_t PndRichCalDb::GetXstep ( ) const

inline

Definition at line 99 of file PndRichCalDb.h.

References fXstep.

{ return fXstep; } 

Double_t PndRichCalDb::GetYmax ( ) const

inline

Definition at line 95 of file PndRichCalDb.h.

References fYmax.

{ return fYmax; } 

Double_t PndRichCalDb::GetYmin ( ) const

inline

Definition at line 90 of file PndRichCalDb.h.

References fYmin.

{ return fYmin; } 

Double_t PndRichCalDb::GetYstep ( ) const

inline

Definition at line 100 of file PndRichCalDb.h.

References fYstep.

{ return fYstep; } 

void PndRichCalDb::Init ( )

virtual

Initialisation (read map from file)

Definition at line 134 of file PndRichCalDb.cxx.

References fFileName, ReadAsciiFile(), and ReadRootFile().

Referenced by PndRichResolution::PndRichResolution().

                        {
  if      (fFileName.EndsWith(".root")) ReadRootFile(fFileName, fName);
  else if (fFileName.EndsWith(".dat"))  ReadAsciiFile(fFileName);
  else {
    LOG(ERROR) << "-E- PndRichCalDb::Init: No proper file name defined! ("<<fFileName.Data()<<")";
    Fatal("Init", "No proper file name");
  }
}

Bool_t PndRichCalDb::IsInside ( Double_t  p, Double_t  x, Double_t  y, Double_t  t, Double_t  f, Int_t &  ip, Int_t &  ix, Int_t &  iy, Int_t &  it, Int_t &  iq, Double_t &  dp, Double_t &  dx, Double_t &  dy, Double_t &  dt, Double_t &  df  )

virtual

Determine whether a point is inside the field map

Parameters

x,y,z	Point coordinates (global) [cm]
ix,iy,iz	(return) Grid cell
dx,dy,dz	(return) Distance from grid point [cm] if inside kTRUE if inside map, else kFALSE

Definition at line 342 of file PndRichCalDb.cxx.

References Double_t, fFmax, fFmin, fFstep, fPmax, fPmin, fPstep, fTmax, fTmin, fTstep, fXmin, fXstep, fYmax, fYmin, and fYstep.

Referenced by GetBetaEff(), GetBetaMean(), and GetBetaSig().

                                                                                                    {

  // ---  Check for being outside the map range
  p = p < fPmax ? p : fPmax-0.00001;
  t = t < fTmax ? t : fTmax-0.01;
  if ( ! ( p >= fPmin && p < fPmax &&
           x >= fXmin && x < fXmax && y >= fYmin && y < fYmax &&
           t >= fTmin && t < fTmax && f >= fFmin && f < fFmax) ) {
    ip = ix = iy = it = iq = 0;
    dp = dx = dy = dt = df = 0.;
    return kFALSE;
  }

  dp = ( p - fPmin ) / fPstep;
  dx = ( x - fXmin ) / fXstep;
  dy = ( y - fYmin ) / fYstep;
  dt = ( t - fTmin ) / fTstep;
  df = ( f - fFmin ) / fFstep;

  // --- Determine grid cell
  ip = Int_t( dp );
  ix = Int_t( dx );
  iy = Int_t( dy );
  it = Int_t( dt );
  iq = Int_t( df );


  // Relative distance from grid point (in units of cell size)
  dp -= Double_t(ip);
  dx -= Double_t(ix);
  dy -= Double_t(iy);
  dt -= Double_t(it);
  df -= Double_t(iq);

  return kTRUE;

}

PndRichCalDb& PndRichCalDb::operator= ( const PndRichCalDb &  )

inlineprivate

Definition at line 125 of file PndRichCalDb.h.

{return *this;};

void PndRichCalDb::Print ( )

virtual

Screen output

Definition at line 464 of file PndRichCalDb.cxx.

References fFmax, fFmin, fFstep, fNf, fNp, fNt, fNx, fNy, fPmax, fPmin, fPstep, fTmax, fTmin, fTstep, fXmax, fXmin, fXstep, fYmax, fYmin, fYstep, and TString.

Referenced by PndRichResolution::PndRichResolution().

                         {
  TString type = "Map";
  if ( fType == 2 ) type = "Soleniod Map ";
  if ( fType == 3 ) type = "Dipole Map ";
  if ( fType == 4 ) type = "Trans Map ";
  cout << "======================================================" << endl;
  cout.precision(4);
  cout << showpoint;
  cout << "----  " << fTitle << " : " << fName << endl;
  cout << "----" << endl;
  cout << "----  Field type     : " << type << endl;
  cout << "----" << endl;
  cout << "----  Field map grid : " << endl;
  cout << "----  p = " << setw(4) << fPmin << " to " << setw(4) << fPmax
       << "   , " << fNp << " grid points, dp = " << fPstep << "   " << endl;
  cout << "----  x = " << setw(4) << fXmin << " to " << setw(4) << fXmax
       << " cm, " << fNx << " grid points, dx = " << fXstep << " cm" << endl;
  cout << "----  y = " << setw(4) << fYmin << " to " << setw(4) << fYmax
       << " cm, " << fNy << " grid points, dy = " << fYstep << " cm" << endl;
  cout << "----  t = " << setw(4) << fTmin << " to " << setw(4) << fTmax
       << "   , " << fNt << " grid points, dt = " << fTstep << "   " << endl;
  cout << "----  f = " << setw(4) << fFmin << " to " << setw(4) << fFmax
       << "   , " << fNf << " grid points, df = " << fFstep << "   " << endl;
  cout << endl;
/*  cout << "----  Field scaling factor: " << fScale << endl;
  Double_t bx = GetBx(0.,0.,0.);
  Double_t by = GetBy(0.,0.,0.);
  Double_t bz = GetBz(0.,0.,0.);
  cout << "----" << endl;
  cout << "----  Field at origin is ( " << setw(6) << bx << ", " << setw(6)
       << by << ", " << setw(6) << bz << ") kG" << endl;*/
 cout << "======================================================" << endl;
}

void PndRichCalDb::ReadAsciiFile ( const char *  fileName )

protected

Read the field map from an ASCII file

Definition at line 516 of file PndRichCalDb.cxx.

References Double_t, fBetaEff, fBetaMean, fBetaSig, fFmax, fFmin, fFstep, fNf, fNp, fNt, fNx, fNy, fPmax, fPmin, fPstep, fTmax, fTmin, fTstep, fXmax, fXmin, fXstep, fYmax, fYmin, fYstep, m, and s.

Referenced by Init().

                                                     {

  // Open file
  cout << "-I- PndRichCalDb: Reading field map from ASCII file "
       << fileName << endl;
  ifstream mapFile(fileName);
  if ( ! mapFile.is_open() ) {
    cerr << "-E- PndRichCalDb:ReadAsciiFile: Could not open file! " << endl;
    Fatal("ReadAsciiFile","Could not open file");
  }

  // Read grid parameters

  mapFile >>fPmin >> fPmax >> fNp;
  mapFile >>fXmin >> fXmax >> fNx;
  mapFile >>fYmin >> fYmax >> fNy;
  mapFile >>fTmin >> fTmax >> fNt;
  mapFile >>fFmin >> fFmax >> fNf;
  fPstep = ( fPmax - fPmin ) / Double_t( fNp - 1 );
  fXstep = ( fXmax - fXmin ) / Double_t( fNx - 1 );
  fYstep = ( fYmax - fYmin ) / Double_t( fNy - 1 );
  fTstep = ( fTmax - fTmin ) / Double_t( fNt - 1 );
  fFstep = ( fFmax - fFmin ) / Double_t( fNf - 1 );

  // Create field arrays
  fBetaMean = new TArrayF(fNp * fNx * fNy * fNt * fNf);
  fBetaSig  = new TArrayF(fNp * fNx * fNy * fNt * fNf);
  fBetaEff  = new TArrayF(fNp * fNx * fNy * fNt * fNf);

  // Read the field values
  cout << right;
  Int_t nTot = fNp * fNx * fNy * fNt * fNf;
  cout << "-I- PndRichCalDb: " << nTot << " entries to read... "
       << setw(3) << 0 << " % ";
  Int_t index = 0;
//  div_t modul;
//  Int_t iDiv = TMath::Nint(nTot/100.);
  Double_t m, s, e;
  for (Int_t ip=0; ip<fNp; ip++) {
     for (Int_t ix=0; ix<fNx; ix++) {
        for (Int_t iy = 0; iy<fNy; iy++) {
           for (Int_t it = 0; it<fNt; it++) {
              for (Int_t iq = 0; iq<fNf; iq++) {
                 if (! mapFile.good()) cerr << "-E- PndRichCalDb::ReadAsciiFile: "
                    << "I/O Error at " << ip << " " << ix << " "
                    << iy << " " << it << " " << iq << endl;
                 index =
                    ip*fNx*fNy*fNt*fNf +
                    ix*fNy*fNt*fNf +
                    iy*fNt*fNf +
                    it*fNf +
                    iq;
/*                 if (iDiv!=0)
                 {
                    modul = div(index,iDiv);
                    if ( modul.rem == 0 ) {
                       Double_t perc = TMath::Nint(100.*index/nTot);
                       cout << "\b\b\b\b\b\b" << setw(3) << perc << " % " << flush;
                    }
                 }*/
                 mapFile >>  m >> s >> e;
                 if (s<0||e<0||e>1)
                    cout << index << " " << m << " " << s << " " << e << endl;
                 fBetaMean->AddAt(m, index);
                 fBetaSig->AddAt(s, index);
                 fBetaEff->AddAt(e, index);
                 if ( mapFile.eof() ) {
                    cerr << endl << "-E- PndRichCalDb::ReadAsciiFile: EOF"
                       << " reached at "  << ip << " " << ix << " " << iy << " " << it << " " << iq << endl;
                    mapFile.close();
                    break;
                 }
              }  // f-Loop
           }  // t-Loop
        }  // y-Loop
     }  // x-Loop
  }  // p-Loop

  cout << "   " << index+1 << " read" << endl;

  mapFile.close();

}

void PndRichCalDb::ReadRootFile ( const char *  fileName, const char *  mapName  )

protected

Read field map from a ROOT file

Definition at line 604 of file PndRichCalDb.cxx.

References exit(), file, and SetCalDb().

Referenced by Init().

                                                     {

  // Store gFile pointer
  TFile* oldFile = gFile;

  // Open root file
  LOG(INFO) << "PndRichCalDb: Reading field map from ROOT file  "<<fileName;
  TFile* file = new TFile(fileName, "READ");
  if (file->IsZombie()) {
    LOG(ERROR) << "-E- PndRichCalDb::ReadRootfile: Cannot read from file! (" << fileName <<")";
    Fatal("ReadRootFile","Cannot read from file");
  }

  // Get the field data object
  PndRichCalDbData* data = NULL;
  file->GetObject(mapName, data);
  if ( ! data ) {
     LOG(ERROR) << "PndRichCalDb::ReadRootFile: data object " << fileName << " not found in file! ";
     exit(-1);
  }

  // Get the field parameters
  SetCalDb(data);

  // Close the root file and delete the data object
  file->Close();
  delete data;
  delete file;
  if ( oldFile ) oldFile->cd();

}

void PndRichCalDb::Reset ( )

protected

Reset the field parameters and data

Definition at line 502 of file PndRichCalDb.cxx.

References fBetaEff, fBetaMean, fBetaSig, fFmax, fFmin, fFstep, fNf, fNp, fNt, fNx, fNy, fPmax, fPmin, fPstep, fTmax, fTmin, fTstep, fXmax, fXmin, fXstep, fYmax, fYmin, and fYstep.

                         {
   fPmin = fXmin = fYmin = fTmin = fFmin = 0;
   fPmax = fXmax = fYmax = fTmax = fFmax = 0;
   fPstep = fXstep = fYstep = fTstep = fFstep = 0;
   fNp = fNx = fNy = fNt = fNf = 0;
   if ( fBetaMean ) { delete fBetaMean; fBetaMean = NULL; }
   if ( fBetaSig ) { delete fBetaSig; fBetaSig = NULL; }
   if ( fBetaEff ) { delete fBetaEff; fBetaEff = NULL; }
}

void PndRichCalDb::SetCalDb ( const PndRichCalDbData *  data )

protected

Set field parameters and data

Definition at line 641 of file PndRichCalDb.cxx.

References Double_t, fBetaEff, fBetaMean, fBetaSig, fFmax, fFmin, fFstep, fNf, fNp, fNt, fNx, fNy, fPmax, fPmin, fPstep, fTmax, fTmin, fTstep, fXmax, fXmin, fXstep, fYmax, fYmin, fYstep, PndRichCalDbData::GetBetaEff(), PndRichCalDbData::GetBetaMean(), PndRichCalDbData::GetBetaSig(), PndRichCalDbData::GetFmax(), PndRichCalDbData::GetFmin(), PndRichCalDbData::GetNf(), PndRichCalDbData::GetNp(), PndRichCalDbData::GetNt(), PndRichCalDbData::GetNx(), PndRichCalDbData::GetNy(), PndRichCalDbData::GetPmax(), PndRichCalDbData::GetPmin(), PndRichCalDbData::GetTmax(), PndRichCalDbData::GetTmin(), PndRichCalDbData::GetType(), PndRichCalDbData::GetXmax(), PndRichCalDbData::GetXmin(), PndRichCalDbData::GetYmax(), and PndRichCalDbData::GetYmin().

Referenced by ReadRootFile().

                                                        {

  // Check compatibility
   std::cout << "fType = " << fType << " " << data->GetType() << std::endl;
   fType = data->GetType();
  if ( data->GetType() != fType ) {
    LOG(ERROR) << "PndRichCalDb::SetField: Incompatible map types Field map is of type "<<fType<<" \n but map on file is of type " << data->GetType();
    Fatal("SetField","Incompatible map types");
  }


  fPmin = data->GetPmin();
  fXmin = data->GetXmin();
  fYmin = data->GetYmin();
  fTmin = data->GetTmin();
  fFmin = data->GetFmin();
  fPmax = data->GetPmax();
  fXmax = data->GetXmax();
  fYmax = data->GetYmax();
  fTmax = data->GetTmax();
  fFmax = data->GetFmax();
  fNp = data->GetNp();
  fNx = data->GetNx();
  fNy = data->GetNy();
  fNt = data->GetNt();
  fNf = data->GetNf();
  fPstep = ( fPmax - fPmin ) / Double_t( fNp - 1 );
  fXstep = ( fXmax - fXmin ) / Double_t( fNx - 1 );
  fYstep = ( fYmax - fYmin ) / Double_t( fNy - 1 );
  fTstep = ( fTmax - fTmin ) / Double_t( fNt - 1 );
  fFstep = ( fFmax - fFmin ) / Double_t( fNf - 1 );
  if ( fBetaMean ) delete fBetaMean;
  if ( fBetaSig ) delete fBetaSig;
  if ( fBetaEff ) delete fBetaEff;
  fBetaMean = new TArrayF(*(data->GetBetaMean()));
  fBetaSig = new TArrayF(*(data->GetBetaSig()));
  fBetaEff = new TArrayF(*(data->GetBetaEff()));

}

void PndRichCalDb::WriteAsciiFile

(

const char *

fileName

)

Write the field map to an ASCII file

Definition at line 386 of file PndRichCalDb.cxx.

References fBetaEff, fBetaMean, fBetaSig, fFmax, fFmin, fNf, fNp, fNt, fNx, fNy, fPmax, fPmin, fTmax, fTmin, fXmax, fXmin, fYmax, and fYmin.

                                                      {

  // Open file
  LOG(INFO) << "PndRichCalDb: Writing field map to ASCII file "<<fileName;
  ofstream mapFile(fileName);
  if ( ! mapFile.is_open() ) {
    LOG(ERROR) << "PndRichCalDb:ReadAsciiFile: Could not open file! ";
    return;
  }

  // Write field map grid parameters
  mapFile.precision(4);
  mapFile << showpoint;

  mapFile << fPmin << " " << fPmax << " " << fNp << endl;
  mapFile << fXmin << " " << fXmax << " " << fNx << endl;
  mapFile << fYmin << " " << fYmax << " " << fNy << endl;
  mapFile << fTmin << " " << fTmax << " " << fNt << endl;
  mapFile << fFmin << " " << fFmax << " " << fNf << endl;

  // Write field values
  cout << right;
  //Int_t nTot = fNp * fNx * fNy * fNt * fNf; //[R.K. 01/2017] unused variable?
  cout << "-I- PndRichCalDb: " << fNp*fNx*fNy*fNt*fNf << " entries to write... "
       << setw(3) << 0 << " % ";
  Int_t index=0;
//  div_t modul;
//  Int_t iDiv = TMath::Nint(nTot/100.);
  for(Int_t ip=0; ip<fNp; ip++) {
     for(Int_t ix=0; ix<fNx; ix++) {
        for(Int_t iy=0; iy<fNy; iy++) {
           for(Int_t it=0; it<fNt; it++) {
              for(Int_t iq=0; iq<fNf; iq++) {
                 index =
                    ip*fNx*fNy*fNt*fNf +
                    ix*fNy*fNt*fNf +
                    iy*fNt*fNf +
                    it*fNf +
                    iq;
/*                 if (iDiv!=0)
                 {
                    modul = div(index,iDiv);
                    if ( modul.rem == 0 ) {
                       Double_t perc = TMath::Nint(100.*index/nTot);
                       cout << "\b\b\b\b\b\b" << setw(3) << perc << " % " << flush;
                    }
                 }*/
                 mapFile << fBetaMean->At(index) << " " << fBetaSig->At(index)
                       << " " << fBetaEff->At(index) << endl;
              }  // f-Loop
           }  // t-Loop
        }  // y-Loop
     }  // x-Loop
  }  // p-Loop
   mapFile.close();

}

void PndRichCalDb::WriteRootFile	(	const char *	fileName,
		const char *	mapName
	)

Write field map data to a ROOT file

Definition at line 448 of file PndRichCalDb.cxx.

References file.

                                                      {

  PndRichCalDbData* data = new PndRichCalDbData(mapName, *this);
  TFile* oldFile = gFile;
  TFile* file = new TFile(fileName, "RECREATE");
  data->Write();
  file->Close();
  if(oldFile) oldFile->cd();

}

Member Data Documentation

TArrayF* PndRichCalDb::fBetaEff

protected

Definition at line 166 of file PndRichCalDb.h.

Referenced by GetBetaEff(), PndRichCalDb(), ReadAsciiFile(), Reset(), SetCalDb(), WriteAsciiFile(), and ~PndRichCalDb().

TArrayF* PndRichCalDb::fBetaMean

protected

Arrays with the field values

Definition at line 164 of file PndRichCalDb.h.

Referenced by GetBetaMean(), PndRichCalDb(), ReadAsciiFile(), Reset(), SetCalDb(), WriteAsciiFile(), and ~PndRichCalDb().

TArrayF* PndRichCalDb::fBetaSig

protected

Definition at line 165 of file PndRichCalDb.h.

Referenced by GetBetaSig(), PndRichCalDb(), ReadAsciiFile(), Reset(), SetCalDb(), WriteAsciiFile(), and ~PndRichCalDb().

TString PndRichCalDb::fFileName

protected

Map file name

Definition at line 148 of file PndRichCalDb.h.

Referenced by GetFileName(), Init(), and PndRichCalDb().

Double_t PndRichCalDb::fFmax

protected

Definition at line 156 of file PndRichCalDb.h.

Referenced by GetFmax(), IsInside(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), SetCalDb(), and WriteAsciiFile().

Double_t PndRichCalDb::fFmin

protected

Definition at line 156 of file PndRichCalDb.h.

Referenced by GetFmin(), IsInside(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), SetCalDb(), and WriteAsciiFile().

Double_t PndRichCalDb::fFstep

protected

Definition at line 156 of file PndRichCalDb.h.

Referenced by GetFstep(), IsInside(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), and SetCalDb().

Int_t PndRichCalDb::fNf

protected

Definition at line 160 of file PndRichCalDb.h.

Referenced by GetBetaEff(), GetBetaMean(), GetBetaSig(), GetNf(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), SetCalDb(), and WriteAsciiFile().

Int_t PndRichCalDb::fNp

protected

Number of grid points

Definition at line 160 of file PndRichCalDb.h.

Referenced by GetNp(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), SetCalDb(), and WriteAsciiFile().

Int_t PndRichCalDb::fNt

protected

Definition at line 160 of file PndRichCalDb.h.

Referenced by GetBetaEff(), GetBetaMean(), GetBetaSig(), GetNt(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), SetCalDb(), and WriteAsciiFile().

Int_t PndRichCalDb::fNx

protected

Definition at line 160 of file PndRichCalDb.h.

Referenced by GetBetaEff(), GetBetaMean(), GetBetaSig(), GetNx(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), SetCalDb(), and WriteAsciiFile().

Int_t PndRichCalDb::fNy

protected

Definition at line 160 of file PndRichCalDb.h.

Referenced by GetBetaEff(), GetBetaMean(), GetBetaSig(), GetNy(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), SetCalDb(), and WriteAsciiFile().

Double_t PndRichCalDb::fPmax

protected

Definition at line 152 of file PndRichCalDb.h.

Referenced by GetPmax(), IsInside(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), SetCalDb(), and WriteAsciiFile().

Double_t PndRichCalDb::fPmin

protected

Field limits in local coordinate system

Definition at line 152 of file PndRichCalDb.h.

Referenced by GetPmin(), IsInside(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), SetCalDb(), and WriteAsciiFile().

Double_t PndRichCalDb::fPstep

protected

Definition at line 152 of file PndRichCalDb.h.

Referenced by GetPstep(), IsInside(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), and SetCalDb().

Double_t PndRichCalDb::fTmax

protected

Definition at line 155 of file PndRichCalDb.h.

Referenced by GetTmax(), IsInside(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), SetCalDb(), and WriteAsciiFile().

Double_t PndRichCalDb::fTmin

protected

Definition at line 155 of file PndRichCalDb.h.

Referenced by GetTmin(), IsInside(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), SetCalDb(), and WriteAsciiFile().

Double_t PndRichCalDb::fTstep

protected

Definition at line 155 of file PndRichCalDb.h.

Referenced by GetTstep(), IsInside(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), and SetCalDb().

Double_t PndRichCalDb::fXmax

protected

Definition at line 153 of file PndRichCalDb.h.

Referenced by GetXmax(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), SetCalDb(), and WriteAsciiFile().

Double_t PndRichCalDb::fXmin

protected

Definition at line 153 of file PndRichCalDb.h.

Referenced by GetXmin(), IsInside(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), SetCalDb(), and WriteAsciiFile().

Double_t PndRichCalDb::fXstep

protected

Definition at line 153 of file PndRichCalDb.h.

Referenced by GetXstep(), IsInside(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), and SetCalDb().

Double_t PndRichCalDb::fYmax

protected

Definition at line 154 of file PndRichCalDb.h.

Referenced by GetYmax(), IsInside(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), SetCalDb(), and WriteAsciiFile().

Double_t PndRichCalDb::fYmin

protected

Definition at line 154 of file PndRichCalDb.h.

Referenced by GetYmin(), IsInside(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), SetCalDb(), and WriteAsciiFile().

Double_t PndRichCalDb::fYstep

protected

Definition at line 154 of file PndRichCalDb.h.

Referenced by GetYstep(), IsInside(), PndRichCalDb(), Print(), ReadAsciiFile(), Reset(), and SetCalDb().

---

The documentation for this class was generated from the following files:

• PndRichCalDb.h
• PndRichCalDb.cxx