FairMesh.cxx

Go to the documentation of this file.

/********************************************************************************
 *    Copyright (C) 2014 GSI Helmholtzzentrum fuer Schwerionenforschung GmbH    *
 *                                                                              *
 *              This software is distributed under the terms of the             *
 *              GNU Lesser General Public Licence (LGPL) version 3,             *
 *                  copied verbatim in the file "LICENSE"                       *
 ********************************************************************************/
// -------------------------------------------------------------------------
// -----                         FairMesh source file                  -----
// -----          original author                  D.Bertini           -----
// -----          adapted april 2010               O.Hartmann          -----
// -------------------------------------------------------------------------
 
#include "FairMesh.h"
 
#include <TMath.h>     // for Sqrt
#include <TString.h>   // for TString, operator+
#include <iostream>    // for operator<<, basic_ostream, etc
 
using namespace std;
 
FairMesh::FairMesh()
    : TObject()
    , fXmin(0.0)
    , fYmin(0.)
    , fZmin(0.)
    , fXmax(0.)
    , fYmax(0.)
    , fZmax(0.)
    , NXbin(0)
    , NYbin(0)
    , NZbin(0)
    , fBinVolume(0.)
    , fDiag(0.)
    , fMeshTid(nullptr)
    , fMeshFlu(nullptr)
    , fMeshSEU(nullptr)
    , fhname("")
{}
 
FairMesh::FairMesh(const char* fname)
    : TObject()
    , fXmin(0.0)
    , fYmin(0.)
    , fZmin(0.)
    , fXmax(0.)
    , fYmax(0.)
    , fZmax(0.)
    , NXbin(0)
    , NYbin(0)
    , NZbin(0)
    , fBinVolume(0.)
    , fDiag(0.)
    , fMeshTid(nullptr)
    , fMeshFlu(nullptr)
    , fMeshSEU(nullptr)
    , fhname(fname)
{
    //  fhname = fname;
}
 
void FairMesh::calculate()
{
    // compute numbers
 
    TString tid = fhname + "TID";
    fMeshTid = new TH2D(tid.Data(), "TID", NXbin, fXmin, fXmax, NYbin, fYmin, fYmax);
    fMeshTid->Sumw2();
 
    TString flu = fhname + "FLU";
    fMeshFlu = new TH2D(flu.Data(), "Fluence", NXbin, fXmin, fXmax, NYbin, fYmin, fYmax);
    fMeshFlu->Sumw2();
 
    TString seu = fhname + "SEU";
    fMeshSEU = new TH2D(seu.Data(), "SEU", NXbin, fXmin, fXmax, NYbin, fYmin, fYmax);
    fMeshSEU->Sumw2();
 
    // reseting
    fMeshTid->Reset();
    fMeshFlu->Reset();
    fMeshSEU->Reset();
 
    fBinVolume = ((fXmax - fXmin) / NXbin) * ((fYmax - fYmin) / NYbin) * ((fZmax - fZmin) / NZbin);
 
    fDiag = TMath::Sqrt(((fXmax - fXmin) / NXbin) * ((fXmax - fXmin) / NXbin)
                        + ((fYmax - fYmin) / NYbin) * ((fYmax - fYmin) / NYbin)
                        + ((fZmax - fZmin) / NZbin) * ((fZmax - fZmin) / NZbin));
}
 
FairMesh::~FairMesh() {}
 
void FairMesh::print()
{
    cout << " Xmin " << fXmin << " Xmax " << fXmax << endl;
    cout << " Ymin " << fYmin << " Ymax " << fYmax << endl;
    cout << " Zmin " << fZmin << " Zmax " << fZmax << endl;
    cout << " NX  " << NXbin << " NY " << NXbin << endl;
}
 
ClassImp(FairMesh);