jdbrice

#include <assert.h>
#include <fstream>
/*********************************************************/
// class for mapping the VMM electronics
// 
// Some terminology
// in the raw data from Tonko we have (SEC, RDO, FEB ...)
// these correspond to our sTGC naming as:
//
// SEC = Plane (1, 2, 3, 4)

jdbrice

FEE[0] : 0000101110011101001001000011100100010110111010111011111011110111000010111001110100100100001110010001011011000101001100011000000000001011100111010010010000111001000101110011100011011001111001010000101110011101001001000011100100010111000100100100110001101110
FEE[1] : 0010101111000111100101011100111001001010000101111001001111010000001010111100011110010101110011100100101000111110001000010100011100101011110001111001010111001110010010100110010010101110101111100010101111000111100101011100111001001010100010110011110000110101
FEE[2] : 1100101101000111100000010111011011001111111101000011011100010100110010110100011110000001011101101101000000011010110001001000101111001011010001111000011101000101000111010010011101101001011101111100101101000111100001110100010100011101000000001101110000000000
FEE[3] : 1110101101110010101110000111001001001001010110001010000111001111111010110111001010111000011100100100100100110010000101000101100011101011011100101011100001110010010010011010010110111100101111011110101101110010101100101

jdbrice

Example:

cd /star/u/jdb/work/ssw/stgc-raw-hit-reader/StRoot/RTS/src/RTS_EXAMPLE starver dev make make ./rts_example /star/data03/pwg/jdb/FWD/daq/stgc-cosmic/st_hltcosmic_22314049_raw_0000002.daq -D stgcsctrg > /star/data03/pwg/jdb/FWD/daq/stgc-cosmic/st_hltcosmic_22314049_raw_0000002.daq.dat

If needed, copy output to wherever the above two scripts are root -b -q -l 'read_RTS.C( "/star/data03/pwg/jdb/FWD/daq/stgc-cosmic/st_hltcosmic_22314049_raw_0000002.daq.dat" )'

jdbrice

<Document  file="StarVMC/Geometry/StgmGeo/StgmGeo1.xml">

    <Replace match="GeV">*1.0E+0;</Replace>    
    <Replace match="MeV">*1.0E-3;</Replace>    
    <Replace match="keV">*1.0E-6;</Replace>

    <Module name="StgmGeo1" comment="Forward Tracking System Geometry / developed for the sTGC tracker">

        <Created date="08/19/2021" />
        <Author  name="Daniel Brandenburg" />

jdbrice

ifstream infile( starlight.c_str() );

int iEvent =0;
while (infile) {
  string sEvent, sVertex, sTrack1, sTrack2;
  if (!getline( infile, sEvent )) break;
  if (!getline( infile, sVertex )) break;
  if (!getline( infile, sTrack1 )) break;
  if (!getline( infile, sTrack2 )) break;

jdbrice

static const Int_t       nrp      = 200;       //  number of R nodes in the map
static const Int_t       nzp      = 800;       //  number of Z nodes in the map
static const Float_t     zm       = 800.0;     //  map max length
static const Float_t     rm       = 400.0;     //  map max radius
static const Int_t       nphi     =  36;
//________________________________________________________________________________
void Agufld(const Char_t *field = "FieldOn") {
//   if ( gClassTable->GetID("TGiant3") >= 0) { // root4star
//     gROOT->LoadMacro("bfc.C");
//     bfc(-1,"MagF,nodefault",0);

jdbrice

TLatex tl;
TString label = "";
TString klabel = "";

TH1* sliceFit(TH2 * h2, TString _klabel){
    TF1 * fgyRes = new TF1( "fgyRes", "gaus" );
    fgyRes->SetRange( -0.5, 0.5 );
    h2->RebinX(4);
    h2->FitSlicesY( fgyRes );

jdbrice

// Globals
TGeoManager *geom;
TGeoElementTable *table;

// Materials
TGeoMedium *mediumAir;
TGeoMedium *mediumG10;
TGeoMedium *mediumNPC02;

jdbrice

#include "TFile.h"
#include "TH1F.h"
#include "TH1D.h"
#include "TH2F.h"
#include "TH3F.h"
#include "TCanvas.h"
#include "TMath.h"
#include "TVector3.h"
#include "TLegend.h"

jdbrice

# Author: Daniel Brandenburg
#
# This program is free software: you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
# 
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the