hcalInf Class Reference

#include <hcalInf.h>

Inheritance diagram for hcalInf:

Collaboration diagram for hcalInf:

Public Member Functions
	hcalInf ()
Double_t	GetXPos () const
Double_t	GetYPos () const
Double_t	GetZPos () const
Double_t	GetModuleSize () const
Int_t	GetNLayers () const
Int_t	GetN1Layers () const
Double_t	GetAbsorber () const
Double_t	GetScin () const
Double_t	GetTyveec () const
Double_t	GetThicknessLayer () const
Int_t	GetXSize () const
Int_t	GetYSize () const
Double_t	GetContainerXSemiAxiss () const
Double_t	GetContainerYSemiAxiss () const
Double_t	GetElectronCut () const
Double_t	GetHadronCut () const
Double_t	GetHcalSize (Int_t num) const
char	GetType (Int_t x, Int_t y) const
Int_t	GetFastMC () const
void	DumpContainer () const
void	FreeInstance ()
Double_t	GetVariable (const char *key)
Double_t	GetVariableStrict (const char *key)
TString	GetStringVariable (const char *key)
void	AddVariable (const char *key, const char *value)

Static Public Member Functions
static hcalInf *	GetInstance (const char *filename)

Protected Member Functions
	hcalInf (const char *filename)
void	CalculateHoleSize ()
virtual	~hcalInf ()

Static Protected Attributes
static hcalInf *	fInf =NULL
static Int_t	fRefCount =0

Private Member Functions
void	InitVariables ()
Bool_t	ExcludeParameter (TString parname)
	hcalInf (const hcalInf &)
hcalInf &	operator= (const hcalInf &)
	ClassDef (hcalInf, 2)

Private Attributes
TMap *	fVariables
TObjArray	fHcalStr
Double_t	fXPos
Double_t	fYPos
Double_t	fZPos
Int_t	fNLayers
Int_t	fN1Layers
Int_t	fXSize
Int_t	fYSize
Double_t	fModuleSize
Double_t	fAbsorber
Double_t	fScin
Double_t	fTyveec
Double_t	fThicknessLayer
Double_t	fCellSize
Double_t	fHcalSize [3]
Double_t	fECut
Double_t	fHCut
Double_t	fSemiX
Double_t	fSemiY
Int_t	fFastMC
Int_t	fSuccess
TString	fFileName

Detailed Description

Definition at line 19 of file hcalInf.h.

Constructor & Destructor Documentation

hcalInf() [1/3]

hcalInf::hcalInf ( )

inline

This is ROOT requirement

Definition at line 23 of file hcalInf.h.

           : TObject(), fVariables(NULL), fHcalStr(), fXPos(0.), fYPos(0.),
    fZPos(0.), fNLayers(0), fXSize(0), fYSize(0), fModuleSize(0.), fAbsorber(0.),
    fScin(0.), fTyveec(0.), fThicknessLayer(0.), fCellSize(0.), fHcalSize(),
    fECut(0.), fHCut(0.), fSemiX(0.0), fSemiY(0.0), fFastMC(-1), 
    fSuccess(-1), fFileName("") 
    {};

hcalInf() [2/3]

hcalInf::hcalInf ( const char *  filename )

protected

If data from Parameter file differs from ours TODO: should understand a way of storring parameters in SHIP Text file constructor

Constructor reads and parses the ascii file, and fill 
the HCAL geometry container

FairRunAna* ana = FairRunAna::Instance(); if (ana!=NULL) { FairRuntimeDb* rtdb=ana->GetRuntimeDb(); rtdb->getContainer("CbmGeoHcalPar"); }

Definition at line 139 of file hcalInf.cxx.

  : TObject(), 
    fVariables(new TMap(200)), 
    fHcalStr(), 
    fXPos(0.), 
    fYPos(0.),
    fZPos(0.), 
    fNLayers(0), 
    fN1Layers(0), 
    fXSize(0), 
    fYSize(0), 
    fModuleSize(0.), 
    fAbsorber(0.), 
    fScin(0.), 
    fTyveec(0.), 
    fThicknessLayer(0.),
    fCellSize(0.), 
    fHcalSize(), 
    fECut(0.), 
    fHCut(0.), 
    fSemiX(0.),
    fSemiY(0.),
    fFastMC(-1),
    fSuccess(1), 
    fFileName(filename) 
{
  std::ifstream file(filename);
  Int_t linenum;
  Double_t val;
  string buffer;
  string message;
  string variable;
  string value;
  TObjString* str=NULL;
  char** err=NULL;
  char winend[2]={13, 0};
  int ssize=-1;
 
  if (!file) {
    cerr << "hcalInf: Can't open information file " << filename << "!" << endl;
    cerr << "hcalInf: Hcal will not be constructed correctly." << endl;
    fSuccess=0;
    return;
  }
 
  linenum=0;
  while(getline(file,buffer)) {
    linenum++;
    message=buffer.substr(buffer.find_first_not_of("    "));    //Skiping initial spaces
    message=message.substr(0,message.find("#"));    //Removing comments
    // Threat windows end of strings correctly
    message=message.substr(0,message.find(winend));
    if (message.empty()) continue;      //This was just a comment
    variable=message.substr(0,message.find("="));
    if (variable=="structure") {
      while(getline(file,buffer)) {
    linenum++;
    if (buffer.empty()) break;
    message=buffer.substr(buffer.find_first_not_of("    "));    //Skiping initial spaces
    message=message.substr(0,message.find("#"));    //Removing comments
    message=message.substr(0,message.find_last_not_of("     ")+1);  //Skiping ending spaces
 
        // Threat windows end of strings correctly
        message=message.substr(0,message.find(winend));
 
    if (!message.empty()) {
      if (-1==ssize)
        ssize=message.size();
      else
        if (ssize!=(Int_t)message.size()) {
          cerr << "Error in HCAL structure at line " << linenum;
          cerr << "." << endl;
          cerr << "Line length differs from previous one" << endl;
              fSuccess=0;
          file.close();
          return;
          
        }
      
      str=new TObjString(message.c_str());
      fHcalStr.Add(str);
    }
      }
      break;
    }
    if (variable==message) {
      cerr << "Syntax error: File " << filename << ".Line " << linenum << "." << endl;
      fSuccess=0;
      file.close();
      return;
    }
    variable=variable.substr(0,variable.find_first_of("     "));
    value=message.substr(message.find("=")+1);
    value=value.substr(value.find_first_not_of("    "));    //Skiping initial spaces
    value=value.substr(0,value.find_first_of("  "));
/*
    value=value.substr(0,value.find_first_not_of("1234567890-+e."));
    val=strtod(value.c_str(),err);
    if (err) {
      cerr << "Syntax error after =: File " << filename << ".Line " << linenum << "." << endl;
      fSuccess=0;
      file.close();
      return;
    }
*/
    AddVariable(variable.c_str(), value.c_str());
  }
  file.close();
  InitVariables();
}

~hcalInf()

hcalInf::~hcalInf ( )

protectedvirtual

Definition at line 28 of file hcalInf.cxx.

{
  for(Int_t i=0;i<fHcalStr.GetSize();i++)
    delete (TObjString*)(fHcalStr.At(i));
  fHcalStr.Clear();
}

hcalInf() [3/3]

hcalInf::hcalInf ( const hcalInf &  )

private

Member Function Documentation

AddVariable()

void hcalInf::AddVariable	(	const char *	key,
		const char *	value
	)

Definition at line 128 of file hcalInf.cxx.

{
  TObjString* skey=(TObjString*)fVariables->FindObject(key);
  //Value for this key already exists!!!
  if (skey!=NULL) return;
  skey=new TObjString(key);
  skey->String().ToLower();
  TObjString* svalue=new TObjString(value);
  fVariables->Add(skey, svalue);
}

CalculateHoleSize()

void hcalInf::CalculateHoleSize ( )

protected

ClassDef()

hcalInf::ClassDef ( hcalInf  , 2    )

private

DumpContainer()

void hcalInf::DumpContainer

(

)

const

Definition at line 363 of file hcalInf.cxx.

{
  // Print out the HCAL geometry container
 
  if (fVariables)
  {
    TObjString* key;
    TIterator* iter=fVariables->MakeIterator();
    while((key=(TObjString*)iter->Next())!=NULL)
    {
      TObjString* str=(TObjString*)fVariables->GetValue(key);
      cout << key->String() << "=" << str->String() << endl;
    }
  }
  TObjString* key;
  TIterator* iter=fHcalStr.MakeIterator();
 
  Int_t modules=0;
  Int_t channels=0;
  Int_t i;
  Int_t j;
  Int_t m[10];
  char stri[2]={0, 0};
  TString st;
  for(i=0;i<10;i++) m[i]=0;
 
  while((key=(TObjString*)iter->Next())!=NULL)
  {
    st=key->String();
    cout << key->String() << endl;
    for(i=0;i<st.Length();i++)
    {
      stri[0]=st[i];
      j=atoi(stri);
      m[j]++;
      if (j) modules++;
      channels+=j*j;
    }
  }
  cout << "Total modules:  " << modules << endl;
  cout << "Total channels: " << channels << endl;
  for(i=1;i<10;i++)
  {
    if (m[i]==0) continue;
    cout << "   Type " << i << " : modules=" << m[i] << ", channels=" << m[i]*i*i << endl;
  }
 
}

ExcludeParameter()

Bool_t hcalInf::ExcludeParameter ( TString  parname )

private

Ignore a parameter during comparision

Definition at line 261 of file hcalInf.cxx.

{
  if (parname.CompareTo("hcalversion")==0) return kTRUE;
  return kFALSE;
}

FreeInstance()

void hcalInf::FreeInstance ( )

inline

Definition at line 168 of file hcalInf.h.

{
  fRefCount--;
  if (fRefCount==0)
  {
    delete this;
    fInf=NULL;
  }
}

GetAbsorber()

Double_t hcalInf::GetAbsorber ( ) const

inline

Definition at line 40 of file hcalInf.h.

{return fAbsorber;}

GetContainerXSemiAxiss()

Double_t hcalInf::GetContainerXSemiAxiss ( ) const

inline

Definition at line 48 of file hcalInf.h.

{return fSemiX;}

GetContainerYSemiAxiss()

Double_t hcalInf::GetContainerYSemiAxiss ( ) const

inline

Definition at line 49 of file hcalInf.h.

{return fSemiY;}

GetElectronCut()

Double_t hcalInf::GetElectronCut ( ) const

inline

Geant cuts information

Definition at line 52 of file hcalInf.h.

{return fECut;}

GetFastMC()

Int_t hcalInf::GetFastMC ( ) const

inline

Definition at line 60 of file hcalInf.h.

{return fFastMC;}

GetHadronCut()

Double_t hcalInf::GetHadronCut ( ) const

inline

Definition at line 53 of file hcalInf.h.

{return fHCut;}

GetHcalSize()

Double_t hcalInf::GetHcalSize ( Int_t  num ) const

inline

Definition at line 55 of file hcalInf.h.

                                               {
    if (num>-1&&num<3) return fHcalSize[num];
    return -1;
  }

GetInstance()

hcalInf * hcalInf::GetInstance ( const char *  filename )

static

Definition at line 35 of file hcalInf.cxx.

{
    if (filename==NULL)
    {
        if (fInf!=NULL)
            fRefCount++;
        return fInf;
    }
    TString newname=gSystem->Getenv("VMCWORKDIR");
    newname+="/geometry/";
    newname+=filename;
    if (fInf!=NULL) {
        if (fInf->fFileName==newname) {
            fRefCount++;
            return fInf;
        } else {
          cerr << "hcalInf: Trying create ";
          cerr << "instance of hcalInf with";
          cerr << " name " << filename;
          cerr << ", which is different from ";
          cerr << fInf->fFileName << "." << endl;
          return NULL;
        }
    }
    fInf=new hcalInf(newname);
    //Is something wrong?
    if (fInf->fSuccess==0)
    {
        delete fInf;
        return NULL;
    }
    fRefCount++;
    return fInf;
}

GetModuleSize()

Double_t hcalInf::GetModuleSize ( ) const

inline

Definition at line 36 of file hcalInf.h.

{return fModuleSize;}

GetN1Layers()

Int_t hcalInf::GetN1Layers ( ) const

inline

Definition at line 39 of file hcalInf.h.

{return fN1Layers;}

GetNLayers()

Int_t hcalInf::GetNLayers ( ) const

inline

Definition at line 38 of file hcalInf.h.

{return fNLayers;}

GetScin()

Double_t hcalInf::GetScin ( ) const

inline

Definition at line 41 of file hcalInf.h.

{return fScin;}

GetStringVariable()

TString hcalInf::GetStringVariable

(

const char *

key

)

Definition at line 103 of file hcalInf.cxx.

{
  TObjString* value=(TObjString*)fVariables->GetValue(key);
  if (value==NULL)
  {
    Fatal("GetStringVariable","Can't find variable named %s.", key);
  }
  return value->GetString();
}

GetThicknessLayer()

Double_t hcalInf::GetThicknessLayer ( ) const

inline

Definition at line 43 of file hcalInf.h.

{return fThicknessLayer;}

GetType()

char hcalInf::GetType ( Int_t  x, Int_t  y  ) const

inline

Returns the type of the module with position (x,y)

Definition at line 158 of file hcalInf.h.

{
  if (x<0||y<0||y>fHcalStr.GetLast()) return 0;
  TObjString* str=(TObjString*)fHcalStr.At(y);
  if (str->GetString().Length()<x) return 0;
  char stri[2]={str->GetString()[x],0};
  return atoi(stri);
}

GetTyveec()

Double_t hcalInf::GetTyveec ( ) const

inline

Definition at line 42 of file hcalInf.h.

{return fTyveec;}

GetVariable()

Double_t hcalInf::GetVariable

(

const char *

key

)

key must be lower case. For example, if have in geo file AaaA=90, then you should call GetVariableStrict("aaaa"). If variable not found, will return -1111

Definition at line 115 of file hcalInf.cxx.

{
  TObjString* value=(TObjString*)fVariables->GetValue(key);
  if (value==NULL)
    return -1111;
  Double_t val;
  char* err=NULL;
  val=strtod(value->GetString(),&err);
  if (err[0]!='\0')
    return -1111;
  return val;
}

GetVariableStrict()

Double_t hcalInf::GetVariableStrict

(

const char *

key

)

key must be lower case. For example, if have in geo file AaaA=90, then you should call GetVariableStrict("aaaa"). If variable not found, will generate Fatal

Definition at line 82 of file hcalInf.cxx.

{
  TObjString* value=(TObjString*)fVariables->GetValue(key);
  if (value==NULL)
  {
    cerr << "Can't find variable named \"" << key << "\"";
    Fatal("GetVariableStrict","Exiting...");
  }
  Double_t val;
  char* err=NULL;
  val=strtod(value->GetString(),&err);
  if (err[0]!='\0')
  {
    cerr << "Can't convert variable named \"" << key ;
    cerr << "\" to floating point. Value is \"";
    cerr << value->GetString() << "\"." << endl;
    Fatal("GetVariableStrict","Exiting...");
  }
  return val;
}

GetXPos()

Double_t hcalInf::GetXPos ( ) const

inline

Getters

Definition at line 32 of file hcalInf.h.

{return fXPos;} 

GetXSize()

Int_t hcalInf::GetXSize ( ) const

inline

Size of HCAL in super modules

Definition at line 46 of file hcalInf.h.

{return fXSize;}

GetYPos()

Double_t hcalInf::GetYPos ( ) const

inline

Definition at line 33 of file hcalInf.h.

{return fYPos;} 

GetYSize()

Int_t hcalInf::GetYSize ( ) const

inline

Definition at line 47 of file hcalInf.h.

{return fYSize;}

GetZPos()

Double_t hcalInf::GetZPos ( ) const

inline

Definition at line 34 of file hcalInf.h.

{return fZPos;}

InitVariables()

void hcalInf::InitVariables ( )

private

Init all other variables from fVariables and fHcalStr

Definition at line 325 of file hcalInf.cxx.

{
  TString stri;
  TObjString* str=(TObjString*)fHcalStr.At(0);
 
  fXPos=GetVariableStrict("xpos");
  fYPos=GetVariableStrict("ypos");
  fZPos=GetVariableStrict("zpos");
  fNLayers=(Int_t)GetVariableStrict("nlayers");
  fN1Layers=(Int_t)GetVariableStrict("n1layers");
  fModuleSize=GetVariableStrict("modulesize");
  fAbsorber=GetVariableStrict("absorber");
  fScin=GetVariableStrict("scin");
  fTyveec=GetVariableStrict("tyveec");
  fHcalSize[2]=GetVariableStrict("hcalzsize");
  fECut=GetVariableStrict("ecut");
  fHCut=GetVariableStrict("hcut");
  fFastMC=(Int_t)GetVariableStrict("fastmc");
  fSemiX=GetVariableStrict("xsemiaxis");
  fSemiY=GetVariableStrict("ysemiaxis");
 
  stri=""; stri+=str->GetString().Length();
  AddVariable("xsize", stri);
  stri=""; stri+=fHcalStr.GetLast()+1;
  AddVariable("ysize", stri);
 
  fXSize=(Int_t)GetVariableStrict("xsize");
  fYSize=(Int_t)GetVariableStrict("ysize");
  fThicknessLayer = 2.0*GetTyveec()+GetScin()+GetAbsorber();
  fHcalSize[0]    = GetXSize()*GetModuleSize();
  fHcalSize[1]    = GetYSize()*GetModuleSize();
  stri=""; stri+=fHcalSize[0];
  AddVariable("xhcalsize",stri);
  stri=""; stri+=fHcalSize[1];
  AddVariable("yhcalsize",stri);
}

operator=()

hcalInf & hcalInf::operator= ( const hcalInf &  )

private

Member Data Documentation

fAbsorber

Double_t hcalInf::fAbsorber

private

Thickness of absorber in one layer [cm]

Definition at line 125 of file hcalInf.h.

fCellSize

Double_t hcalInf::fCellSize

private

transverse size of HCAL cell for simulation [cm]

Definition at line 133 of file hcalInf.h.

fECut

Double_t hcalInf::fECut

private

Electron cut for HCAL

Definition at line 137 of file hcalInf.h.

fFastMC

Int_t hcalInf::fFastMC

private

Flag to run Fast (1) or Full (0) MC code

Definition at line 145 of file hcalInf.h.

fFileName

TString hcalInf::fFileName

private

Definition at line 150 of file hcalInf.h.

fHcalSize

Double_t hcalInf::fHcalSize[3]

private

Size of HCAL container [cm]

Definition at line 135 of file hcalInf.h.

fHcalStr

TObjArray hcalInf::fHcalStr

private

Structure of HCAL as array of strings This variable should be saved in parameter file

Definition at line 103 of file hcalInf.h.

fHCut

Double_t hcalInf::fHCut

private

Hadron cut for HCAL

Definition at line 139 of file hcalInf.h.

fInf

hcalInf * hcalInf::fInf =NULL

staticprotected

Definition at line 88 of file hcalInf.h.

fModuleSize

Double_t hcalInf::fModuleSize

private

transverse supermodule size in cm

Definition at line 122 of file hcalInf.h.

fN1Layers

Int_t hcalInf::fN1Layers

private

Number of layers in first hcal section

Definition at line 114 of file hcalInf.h.

fNLayers

Int_t hcalInf::fNLayers

private

Number of HCAL layers

Definition at line 112 of file hcalInf.h.

fRefCount

Int_t hcalInf::fRefCount =0

staticprotected

Definition at line 89 of file hcalInf.h.

fScin

Double_t hcalInf::fScin

private

Thickness of scintillator in one layer [cm]

Definition at line 127 of file hcalInf.h.

fSemiX

Double_t hcalInf::fSemiX

private

Semiaxises of keeping volume

Definition at line 141 of file hcalInf.h.

fSemiY

Double_t hcalInf::fSemiY

private

Definition at line 142 of file hcalInf.h.

fSuccess

Int_t hcalInf::fSuccess

private

1 if evething Ok

Definition at line 148 of file hcalInf.h.

fThicknessLayer

Double_t hcalInf::fThicknessLayer

private

Total thickness of layer [cm]

Definition at line 131 of file hcalInf.h.

fTyveec

Double_t hcalInf::fTyveec

private

Thickness of tyvec in one layer [cm]

Definition at line 129 of file hcalInf.h.

fVariables

TMap* hcalInf::fVariables

private

A map containing all variables This variable should be saved in parameter file

Definition at line 100 of file hcalInf.h.

fXPos

Double_t hcalInf::fXPos

private

x-position of HCAL center [cm]

Definition at line 105 of file hcalInf.h.

fXSize

Int_t hcalInf::fXSize

private

x-size of hcal in supermodules

Definition at line 117 of file hcalInf.h.

fYPos

Double_t hcalInf::fYPos

private

y-position of HCAL center [cm]

Definition at line 107 of file hcalInf.h.

fYSize

Int_t hcalInf::fYSize

private

y-size of hcal in supermodules

Definition at line 119 of file hcalInf.h.

fZPos

Double_t hcalInf::fZPos

private

z-position of HCAL front [cm]

Definition at line 109 of file hcalInf.h.

---

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

• hcal/hcalInf.h
• hcal/hcalInf.cxx