Removes buggy matbudget macros

All mcbm material budget map generating macros were found to have a halfing of radlen bug and a missing volume bug. To prevent these macros from being further used, they are removed.

macro/mcbm/matbudget_ana_mcbm_mvd.C

-/* Copyright (C) 2015 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
-   SPDX-License-Identifier: GPL-3.0-only
-   Authors: David Emschermann [committer] */
-
-#if !defined(__CINT__) || defined(__MAKECINT__)
-#include "FairRadLenPoint.h"
-
-#include "TCanvas.h"
-#include "TClonesArray.h"
-#include "TFile.h"
-#include "TH1.h"
-#include "TProfile2D.h"
-#include "TROOT.h"
-#include "TString.h"
-#include "TStyle.h"
-#include "TSystem.h"
-#include "TTree.h"
-#include "TVector3.h"
-
-#include <iostream>
-#include <vector>
-using std::cout;
-using std::endl;
-using std::vector;
-#endif
-
-// need to set this root include path:
-// include path: -I/opt/cbm/fairsoft_jul15p1/installation/include/root -I/opt/cbm/fairroot_v-15.07-fairsoft_jul15p1/include -I/opt/cbm/fairsoft_jul15p1/installation/include
-// .include $SIMPATH/include
-// .include $FAIRROOTPATH/include
-// .x 'matbudget_ana_mcbm_mvd.C++("mvd_v18a")'
-
-Int_t matbudget_ana_mcbm_mvd(const char* inGeo, Int_t nEvents = 10000000)
-{
-
-  // Input file (MC)
-  TString geoVersion(inGeo);
-  TString inFile = "data/matbudget." + geoVersion + ".mc.root";
-  TFile* input   = new TFile(inFile);
-  if (!input) {  // this is not checked in compiled mode
-    cout << "*** matbudget_ana: Input file " << inFile << " not found!\n"
-         << "Be sure to run matbudget_mc.C before for the respective " + geoVersion + " geometry!" << endl;
-    exit;
-  }
-
-  // Output file (material maps)
-  TString outFile = geoVersion + "_matbudget.root";
-
-  // Input tree and branch
-  TTree* tree          = (TTree*) input->Get("cbmsim");
-  TClonesArray* points = new TClonesArray("FairRadLenPoint");
-  tree->SetBranchAddress("RadLen", &points);
-
-  // Create output histograms
-  TH1D* hisRadLen = new TH1D("hisRadLen", "Radiation Length", 1000, 0, 100);
-  //  const int nStations = 8;       // number of STS stations
-  //  const int nBins     = 1000;    // number of bins in histograms (in both x and y)
-  //  const int rMax      = 55;      // maximal radius for histograms (for both x and y)
-  const int nStations = 2;    // number of STS stations
-  const int nBins     = 300;  // number of bins in histograms (in both x and y)
-  //  const int rMax      = 16;      // maximal radius for histograms (for both x and y)
-  const int rMax = 8;  // maximal radius for histograms (for both x and y)
-  TProfile2D* hStaRadLen[nStations];
-  for (int i = 0; i < nStations; ++i) {
-    TString name = "Radiation Thickness [%],";
-    name += " Station";
-    //    name += i+1;  // STS
-    name += i;  // MVD
-    hStaRadLen[i] = new TProfile2D(name, name, nBins, -rMax, rMax, nBins, -rMax, rMax);
-  }
-
-  // Auxiliary variables
-  TVector3 vecs, veco;
-  std::map<int, int> trackHitMap;
-
-
-  // Event loop
-  int firstEvent = 0;
-  for (Int_t event = firstEvent; event < (firstEvent + nEvents) && event < tree->GetEntriesFast(); event++) {
-    tree->GetEntry(event);
-    if (0 == event % 100000) cout << "*** Processing event " << event << endl;
-
-    const int nTracks = 1;
-    std::vector<double> RadLengthOnTrack(nTracks, 0.0);  //trackID, vector with points on track
-    std::vector<double> TrackLength(nTracks, 0.0);       //trackID, vector with points on track
-
-    vector<double> RadThick(nStations, 0);
-    double x, y;
-
-    // For this implementation to be correct, there should be only one MCTrack per event.
-
-    // Point loop
-    for (Int_t iPoint = 0; iPoint < points->GetEntriesFast(); iPoint++) {
-      FairRadLenPoint* point = (FairRadLenPoint*) points->At(iPoint);
-
-      // Get track position at entry and exit of material
-      TVector3 posIn, posOut, posDif;
-      posIn  = point->GetPosition();
-      posOut = point->GetPositionOut();
-      posDif = 0.5 * (posOut - posIn);
-
-      // Midpoint between in and out
-      const TVector3 middle = (posOut + posIn);
-      x                     = middle.X() / 2;
-      y                     = middle.Y() / 2;
-      const double z        = middle.Z() / 2;
-
-      // Material budget per unit length
-      const double radThick = posDif.Mag() / point->GetRadLength();
-      RadLengthOnTrack[point->GetTrackID()] += radThick;
-      TrackLength[point->GetTrackID()] += posDif.Mag();
-
-      // Determine station number
-      //      int iStation    =  posIn.Z()/10  - 3 + 0.5; // STS // suppose equidistant stations at 30-100 cm
-      //      int iStationOut =  posOut.Z()/10 - 3 + 0.5; // STS
-      int iStation    = posIn.Z() / 5 - 1 + 0.5;   // MVD
-      int iStationOut = posOut.Z() / 5 - 1 + 0.5;  // MVD
-      if (iStationOut != iStation) continue;
-      if (iStation >= nStations || iStation < 0) continue;
-      RadThick[iStation] += radThick;
-    }
-
-    // Fill material budget map for each station
-    for (int i = 0; i < nStations; ++i)
-      hStaRadLen[i]->Fill(x, y, RadThick[i] * 100);
-
-    for (int k = 0; k < RadLengthOnTrack.size(); k++)
-      if (RadLengthOnTrack[k] > 0) hisRadLen->Fill(RadLengthOnTrack[k]);
-
-  }  // event loop
-
-
-  // Plotting the results
-  TCanvas* can1 = new TCanvas();
-  //  can1->Divide(nStations/2,2);
-  can1->Divide(2, nStations / 2);
-  gStyle->SetPalette(1);
-  gStyle->SetOptStat(0);
-
-  // Open output file
-  TFile* output = new TFile(outFile, "RECREATE");
-  output->cd();
-
-  for (int iStation = 0; iStation < nStations; iStation++) {
-    can1->cd(iStation + 1);
-    hStaRadLen[iStation]->GetXaxis()->SetTitle("x [cm]");
-    hStaRadLen[iStation]->GetYaxis()->SetTitle("y [cm]");
-    //hStaRadLen[i]->GetZaxis()->SetTitle("radiation thickness [%]");
-    //    hStaRadLen[iStation]->SetAxisRange(0, 2, "Z");  // STS
-    hStaRadLen[iStation]->SetAxisRange(0, 0.5, "Z");  // MVD
-    hStaRadLen[iStation]->Draw("colz");
-    hStaRadLen[iStation]->Write();
-  }
-
-  // Plot file
-  TString plotFile = geoVersion + "_matbudget.png";
-  can1->SaveAs(plotFile);
-
-  // Close files
-  input->Close();
-  output->Close();
-  cout << "Material budget maps written to " << outFile << endl;
-
-  return 0;
-}