createPndLmdMisalignmentMatrices.C File Reference

#include <TRandom3.h>

Go to the source code of this file.

Functions
TGeoHMatrix	createRandomMatrix (double angleSigma, double shiftSigma)
int	createPndLmdMisalignmentMatrices (bool debug=false)

Variables
unsigned int	seed = 128
TRandom3 *	PRNG = new TRandom3(seed)

Function Documentation

int createPndLmdMisalignmentMatrices

(

bool

debug = false

)

Definition at line 49 of file createPndLmdMisalignmentMatrices.C.

References Cave, createRandomMatrix(), Dipole, fRun, PndLmdGeometryHelper::getAllAlignPaths(), PndLmdGeometryHelper::getInstance(), i, Magnet, Pipe, printf(), rot, PndSdsDetector::SetExclusiveSensorType(), simOutput, and TString.

                                                         {

        string misMatricesFilePath = "misalignMatrices-SensorsOnly.root";

        cout << "creating dummy geometry...\n";

        // the Geometry we want to use
        TString geometryFile = "Luminosity-Detector.root";

        // we have to setup a dummy simulation so the gGeoManager indexes the Geometry
        TString simOutput = "./dummy.root";
        FairRunSim *fRun = new FairRunSim();
        fRun->SetName("TGeant4");
        fRun->SetOutputFile(simOutput);
        fRun->SetMaterials("media_pnd.geo");
        FairModule *Cave = new PndCave("CAVE");
        Cave->SetGeometryFileName("pndcave.geo");
        fRun->AddModule(Cave);
        FairModule *Pipe = new PndPipe("PIPE");
        Pipe->SetGeometryFileName("beampipe_201309.root");
        fRun->AddModule(Pipe);
        FairModule *Magnet = new PndMagnet("MAGNET");
        Magnet->SetGeometryFileName("FullSolenoid_V842.root");
        fRun->AddModule(Magnet);
        FairModule *Dipole = new PndMagnet("MAGNET");
        Dipole->SetGeometryFileName("dipole.geo");
        fRun->AddModule(Dipole);
        PndLmdDetector *Lum = new PndLmdDetector("LUM", kTRUE);
        Lum->SetExclusiveSensorType("LumActive");  //ignore MVD
        Lum->SetGeometryFileName(geometryFile);
        fRun->AddModule(Lum);
        fRun->Init();

        // after init, the geometry can't be changed anymore, but that's okay.
        // we only want to create a matrix file that can be used in a later step.

        cout << "PndLmd: creating misalignment matrices.\n";

        PndLmdGeometryHelper &helper = PndLmdGeometryHelper::getInstance();
        vector < string > paths = helper.getAllAlignPaths(true, false, false, false, false);

        cout << "got paths.\n";

        std::map < std::string, TGeoHMatrix > matrices;

        // for testing if ONE matrix works
        if (debug) {

                misMatricesFilePath = "misalignMatrices-Sensor-0.root";

                // Misalign geometry
                std::string examplePath = "/cave_1/lmd_root_0/half_0/plane_0/module_0/sensor_0/";

                TGeoHMatrix misalignTestMat;
                double translation[3] = { 0.5, 0.5, 0.0 };
                misalignTestMat.SetTranslation(translation);
                matrices[examplePath] = misalignTestMat;

        } else {

                for (auto &i : paths) {
                        // TODO: cerate sigma parameters another way

                        double shift = 1000e-4;                         // this is 1 mm!
                        double rot = TMath::RadToDeg() * 1000e-3;       // this is... large

                        TGeoHMatrix tempMat = createRandomMatrix(rot, shift);
                        matrices[i] = tempMat;
                }
        }

        if (true) {
                // checking for good measure
                for (auto &i : matrices) {
                        cout << "\n---\n" << i.first << "\n";
                        i.second.Print();
                }
        }

        TFile *misalignmentMatrixRootfile = new TFile(misMatricesFilePath.c_str(), "NEW");
        if (misalignmentMatrixRootfile->IsOpen()) {
                printf("File opened successfully\n");

                gDirectory->WriteObject(&matrices, "PndLmdMisalignMatrices");
                misalignmentMatrixRootfile->Write();
                misalignmentMatrixRootfile->Close();

                cout << "All matrices written to file!\n";
        } else {
                cerr << "WARNING! Could not write to " << misMatricesFilePath << "\n";
        }
        return 0;
}

TGeoHMatrix createRandomMatrix	(	double	angleSigma,
		double	shiftSigma
	)

Definition at line 13 of file createPndLmdMisalignmentMatrices.C.

References mean, and PRNG.

Referenced by createPndLmdMisalignmentMatrices().

                                                                     {

        double mean = 0;

        double sigmaX, sigmaY, sigmaZ;
        double shift[3];

        // can only rotate about z
        //sigmaX = PRNG->Gaus(mean, angleSigma);
        //sigmaY = PRNG->Gaus(mean, angleSigma);
        sigmaX = 0;
        sigmaY = 0;
        sigmaZ = PRNG->Gaus(mean, angleSigma);

        // can't move in z
        shift[0] = PRNG->Gaus(mean, shiftSigma);
        shift[1] = PRNG->Gaus(mean, shiftSigma);
        //shift[2] = PRNG->Gaus(mean, shiftSigma);
        shift[2] = 0;

//      cout << "rand val: " << angles[0] << "\n";
//      cout << "rand val: " << angles[1] << "\n";
//      cout << "rand val: " << angles[2] << "\n";
//      cout << "rand val: " << shift[0] << "\n";
//      cout << "rand val: " << shift[1] << "\n";
//      cout << "rand val: " << shift[2] << "\n";

        TGeoHMatrix result;
        result.RotateZ(sigmaZ);
        result.SetTranslation(shift);

//      result.Print();

        return result;
}

Variable Documentation

TRandom3* PRNG = new TRandom3(seed)

Definition at line 10 of file createPndLmdMisalignmentMatrices.C.

Referenced by createRandomMatrix().

unsigned int seed = 128

Definition at line 5 of file createPndLmdMisalignmentMatrices.C.

Referenced by convolutionAnalysis(), emc(), emc_correction_data_production(), emc_correction_QA_data_production(), EvtRootRandomEngine::EvtRootRandomEngine(), fittest(), chigen::initialize(), poormantracks(), EvtRootRandomEngine::random(), rlxd_init(), rlxs_init(), run_sim(), run_sim_llbar_evtgen(), run_sim_stt_dpm(), run_sim_stt_evt(), run_sim_stt_evt_fix(), run_sim_sttcombi_dpm(), run_sim_sttcombi_evtgen(), run_sim_sttcombi_pgun(), run_sim_tpc_dpm(), run_sim_tpc_evt(), run_sim_tpccombi_dpm(), run_sim_tpcmvd(), runLumi0DPMbkgSte(), runLumi0SimBox(), runLumiPixel0SimBox(), runLumiPixel1bDigiNoise(), runLumiPixel1cDigiSigNoise(), runpmt(), runsim(), PndLepLepGenerator::SetSeed(), PndPiPiGenerator::SetSeed(), sim(), sim_complete_pbarA(), sim_gg(), sim_LMD(), sim_muo_dub(), sim_pid_runSimHF_GiB_DC(), sim_radgrid(), sim_sttcombi_pgun(), and simInclusive().