create_rootgeom_TOF700_support.h

Go to the documentation of this file.

//---------------------------
 
#define TOF2_MAX_UBEAMS 42
#define TOF2_UBEAMS_ORIENTATIONS 4
 
// Detector's construct parameters
Double_t Support_TOF2_long = 3.;
Double_t Support_TOF2_short = 2.;
Double_t Support_TOF2_thickness = 0.15;
 
Double_t Support_TOF2_X_center[TOF2_MAX_UBEAMS] = {0.};
Double_t Support_TOF2_Y_center[TOF2_MAX_UBEAMS] = {0.};
Double_t Support_TOF2_Z_center[TOF2_MAX_UBEAMS] = {0.};
Double_t Support_TOF2_length[TOF2_MAX_UBEAMS] = {0.};
Int_t Support_TOF2_orientation[TOF2_MAX_UBEAMS] = {0};
 
int readGeomSupport(char* geomfile);
 
void create_rootgeom_TOF700_support()
{
 
    // -------   Geometry text file name (input)   ----------------------------------
    readGeomSupport((char*)"TOF700_geometry_support.txt");
 
    TGeoTranslation* ChannelPosTOF2_trans[TOF2_MAX_UBEAMS] = {NULL};
    TGeoTranslation* AirPosTOF2_trans[4] = {NULL};
    char name[64];
    for (int i = 0; i < TOF2_MAX_UBEAMS; i++) {
        sprintf(name, "Channel_%d_PosTOF2_trans", i + 1);
        ChannelPosTOF2_trans[i] =
            new TGeoTranslation(name, Support_TOF2_X_center[i], Support_TOF2_Y_center[i], Support_TOF2_Z_center[i]);
    }
    AirPosTOF2_trans[0] = new TGeoTranslation("AirPos0", 0., +Support_TOF2_thickness, 0.);
    AirPosTOF2_trans[1] = new TGeoTranslation("AirPos1", 0., -Support_TOF2_thickness, 0.);
    AirPosTOF2_trans[2] = new TGeoTranslation("AirPos2", +Support_TOF2_thickness, 0., 0.);
    AirPosTOF2_trans[3] = new TGeoTranslation("AirPos3", -Support_TOF2_thickness, 0., 0.);
 
    char nTOF2ChannelS[TOF2_MAX_UBEAMS][32];
    char nTOF2ChannelAirS[TOF2_MAX_UBEAMS][32];
    char nTOF2ChannelV[TOF2_MAX_UBEAMS][32];
    char nTOF2ChannelAirV[TOF2_MAX_UBEAMS][32];
    TGeoBBox* TOF2ChannelS[TOF2_MAX_UBEAMS];
    TGeoBBox* TOF2ChannelAirS[TOF2_MAX_UBEAMS];
    TGeoVolume* TOF2ChannelV[TOF2_MAX_UBEAMS];
    TGeoVolume* TOF2ChannelAirV[TOF2_MAX_UBEAMS];
 
    for (int i = 0; i < TOF2_MAX_UBEAMS; i++) {
        sprintf(nTOF2ChannelS[i], "TOF2ChannelS%d", i + 1);
        sprintf(nTOF2ChannelAirS[i], "TOF2ChannelAirS%d", i + 1);
        sprintf(nTOF2ChannelV[i], "TOF2ChannelV%d", i + 1);
        sprintf(nTOF2ChannelAirV[i], "TOF2ChannelAirV%d", i + 1);
        switch (Support_TOF2_orientation[i]) {
            case 0:
                TOF2ChannelS[i] =
                    new TGeoBBox(nTOF2ChannelS[i], Support_TOF2_length[i] / 2., Support_TOF2_short, Support_TOF2_long);
                TOF2ChannelAirS[i] = new TGeoBBox(nTOF2ChannelAirS[i], Support_TOF2_length[i] / 2.,
                                                  Support_TOF2_short - Support_TOF2_thickness,
                                                  Support_TOF2_long - Support_TOF2_thickness);
 
                // Volumes
 
                TOF2ChannelV[i] = new TGeoVolume(nTOF2ChannelV[i], TOF2ChannelS[i]);
                TOF2ChannelV[i]->SetMedium(medAl);
                TOF2ChannelV[i]->SetLineColor(kYellow);
                TOF2ChannelV[i]->SetVisibility(kTRUE);
                // TOF2ChannelV[i]->SetTransparency(80);
 
                TOF2ChannelAirV[i] = new TGeoVolume(nTOF2ChannelAirV[i], TOF2ChannelAirS[i]);
                TOF2ChannelAirV[i]->SetMedium(pMedAir);
                //      TOF2ChannelAirV[i]->SetVisibility(kFALSE);
                //      TOF2ChannelAirV[i]->SetLineColor(kGreen);
                TOF2ChannelAirV[i]->SetVisibility(kTRUE);
                TOF2ChannelAirV[i]->SetLineColor(kWhite);
 
                TOF2ChannelV[i]->AddNode(TOF2ChannelAirV[i], 1, AirPosTOF2_trans[0]);
                break;
            case 1:
                TOF2ChannelS[i] =
                    new TGeoBBox(nTOF2ChannelS[i], Support_TOF2_length[i] / 2., Support_TOF2_short, Support_TOF2_long);
                TOF2ChannelAirS[i] = new TGeoBBox(nTOF2ChannelAirS[i], Support_TOF2_length[i] / 2.,
                                                  Support_TOF2_short - Support_TOF2_thickness,
                                                  Support_TOF2_long - Support_TOF2_thickness);
 
                // Volumes
 
                TOF2ChannelV[i] = new TGeoVolume(nTOF2ChannelV[i], TOF2ChannelS[i]);
                TOF2ChannelV[i]->SetMedium(medAl);
                TOF2ChannelV[i]->SetLineColor(kYellow);
                TOF2ChannelV[i]->SetVisibility(kTRUE);
                // TOF2ChannelV[i]->SetTransparency(80);
 
                TOF2ChannelAirV[i] = new TGeoVolume(nTOF2ChannelAirV[i], TOF2ChannelAirS[i]);
                TOF2ChannelAirV[i]->SetMedium(pMedAir);
                //      TOF2ChannelAirV[i]->SetVisibility(kFALSE);
                //      TOF2ChannelAirV[i]->SetLineColor(kGreen);
                TOF2ChannelAirV[i]->SetVisibility(kTRUE);
                TOF2ChannelAirV[i]->SetLineColor(kWhite);
 
                TOF2ChannelV[i]->AddNode(TOF2ChannelAirV[i], 1, AirPosTOF2_trans[1]);
                break;
            case 2:
                TOF2ChannelS[i] =
                    new TGeoBBox(nTOF2ChannelS[i], Support_TOF2_short, Support_TOF2_length[i] / 2., Support_TOF2_long);
                TOF2ChannelAirS[i] =
                    new TGeoBBox(nTOF2ChannelAirS[i], Support_TOF2_short - Support_TOF2_thickness,
                                 Support_TOF2_length[i] / 2., Support_TOF2_long - Support_TOF2_thickness);
 
                // Volumes
 
                TOF2ChannelV[i] = new TGeoVolume(nTOF2ChannelV[i], TOF2ChannelS[i]);
                TOF2ChannelV[i]->SetMedium(medAl);
                TOF2ChannelV[i]->SetLineColor(kYellow);
                TOF2ChannelV[i]->SetVisibility(kTRUE);
                // TOF2ChannelV[i]->SetTransparency(80);
 
                TOF2ChannelAirV[i] = new TGeoVolume(nTOF2ChannelAirV[i], TOF2ChannelAirS[i]);
                TOF2ChannelAirV[i]->SetMedium(pMedAir);
                //      TOF2ChannelAirV[i]->SetVisibility(kFALSE);
                //      TOF2ChannelAirV[i]->SetLineColor(kGreen);
                TOF2ChannelAirV[i]->SetVisibility(kTRUE);
                TOF2ChannelAirV[i]->SetLineColor(kWhite);
 
                TOF2ChannelV[i]->AddNode(TOF2ChannelAirV[i], 1, AirPosTOF2_trans[2]);
                break;
            case 3:
                TOF2ChannelS[i] =
                    new TGeoBBox(nTOF2ChannelS[i], Support_TOF2_short, Support_TOF2_length[i] / 2., Support_TOF2_long);
                TOF2ChannelAirS[i] =
                    new TGeoBBox(nTOF2ChannelAirS[i], Support_TOF2_short - Support_TOF2_thickness,
                                 Support_TOF2_length[i] / 2., Support_TOF2_long - Support_TOF2_thickness);
 
                // Volumes
 
                TOF2ChannelV[i] = new TGeoVolume(nTOF2ChannelV[i], TOF2ChannelS[i]);
                TOF2ChannelV[i]->SetMedium(medAl);
                TOF2ChannelV[i]->SetLineColor(kYellow);
                TOF2ChannelV[i]->SetVisibility(kTRUE);
                // TOF2ChannelV[i]->SetTransparency(80);
 
                TOF2ChannelAirV[i] = new TGeoVolume(nTOF2ChannelAirV[i], TOF2ChannelAirS[i]);
                TOF2ChannelAirV[i]->SetMedium(pMedAir);
                //      TOF2ChannelAirV[i]->SetVisibility(kFALSE);
                //      TOF2ChannelAirV[i]->SetLineColor(kGreen);
                TOF2ChannelAirV[i]->SetVisibility(kTRUE);
                TOF2ChannelAirV[i]->SetLineColor(kWhite);
 
                TOF2ChannelV[i]->AddNode(TOF2ChannelAirV[i], 1, AirPosTOF2_trans[3]);
                break;
            default:
                break;
        }
    }
 
    // Adding volumes to the TOP Volume
    for (int i = 0; i < TOF2_MAX_UBEAMS; i++) {
        TOF2Top->AddNode(TOF2ChannelV[i], i + 1, ChannelPosTOF2_trans[i]);
    }
}
 
int readGeomSupport(char* geomfile)
{
    char fname[128];
    FILE* fg = 0;
    int n = 0;
    if (strlen(geomfile) == 0) {
        printf("TOF700 geometry file name not defined!\n");
        return 0;
    }
    TString dir = getenv("VMCWORKDIR");
    sprintf(fname, "%s/geometry/%s", dir.Data(), geomfile);
    fg = fopen(fname, "r");
    if (fg == NULL) {
        printf("TOF700 geometry file %s open error!\n", fname);
        return 0;
    }
 
    int c = 0;
    while (fscanf(fg, "%d %lf %lf %lf %lf %d\n", &n, &Support_TOF2_X_center[c], &Support_TOF2_Y_center[c],
                  &Support_TOF2_Z_center[c], &Support_TOF2_length[c], &Support_TOF2_orientation[c])
           == 6)
    {
        Support_TOF2_X_center[c] += xoffs;
        Support_TOF2_Y_center[c] += yoffs;
        Support_TOF2_Z_center[c] += zoffs;
        c++;
        if (c >= TOF2_MAX_UBEAMS)
            break;
    }
    fclose(fg);
    return 1;
}