Changeset 493aa9b in git for examples

Timestamp:

Oct 10, 2014, 10:10:54 PM (10 years ago)

Author:

Christophe Delaere <christophe.delaere@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

c4d879d

Parents:

e0c1ff0

Message:

First working version

Many features are still missing, but that version provides the minimal
features in EventDisplay and in addition: event control and transluscent
geometry.

File:

• examples/geometry.C (modified) (8 diffs)

examples/geometry.C

@@ -157 +157 @@
    TGeoMaterial *matVacuum = new TGeoMaterial("Vacuum", 0,0,0);
    TGeoMaterial *matAl = new TGeoMaterial("Al", 26.98,13,2.7); // placeholder
+   matVacuum->SetTransparency(85); //TODO: tune
+   matAl->SetTransparency(85); //TODO: tune
 
    //--- define some media
@@ -188 +190 @@
    confReader->ReadFile(configFile);
 
-   tk_radius_ = confReader->GetDouble(Form("%s::Radius",ParticlePropagator), 1.0)*100;          // tk_radius
-   tk_length_ = confReader->GetDouble(Form("%s::HalfLength",ParticlePropagator), 3.0)*100;      // tk_length
+   tk_radius_ = confReader->GetDouble(Form("%s::Radius",ParticlePropagator), 1.0)*100.;         // tk_radius
+   tk_length_ = confReader->GetDouble(Form("%s::HalfLength",ParticlePropagator), 3.0)*100.;     // tk_length
    tk_Bz_     = confReader->GetDouble("ParticlePropagator::Bz", 0.0);                           // tk_Bz
 
@@ -491 +493 @@
    gEtaAxis = det3D.getCaloAxes().first;
    gPhiAxis = det3D.getCaloAxes().second;
+
+
+/*
+   // make the top container volume -> designed to contain a "big" detector (ATLAS)
+   TGeoVolume *top = geom->MakeBox("TOP", 0, 1500, 1500, 2300);
+   geom->SetTopVolume(top);
+
+   // build the detector
+   top->AddNode(det3D.getDetector(true),1);
+   geom->CloseGeometry();
+
+   gGeoManager->DefaultColors();
+
+   TGeoVolume* top = gGeoManager->GetTopVolume()->FindNode("Delphes3DGeometry_1")->GetVolume();
+
+   TEveGeoTopNode* trk = new TEveGeoTopNode(gGeoManager, top->FindNode("tracker_1"));
+   trk->SetVisLevel(6);
+   gEve->AddGlobalElement(trk);
+
+   TEveGeoTopNode* calo = new TEveGeoTopNode(gGeoManager, top->FindNode("Calorimeter_barrel_1"));
+   calo->SetVisLevel(3);
+   gEve->AddGlobalElement(calo);
+   calo = new TEveGeoTopNode(gGeoManager, top->FindNode("Calorimeter_endcap_1"));
+   calo->SetVisLevel(3);
+   calo->UseNodeTrans();
+   gEve->AddGlobalElement(calo);
+   calo = new TEveGeoTopNode(gGeoManager, top->FindNode("Calorimeter_endcap_2"));
+   calo->SetVisLevel(3);
+   calo->UseNodeTrans();
+   gEve->AddGlobalElement(calo);
+
+   TEveGeoTopNode* muon = new TEveGeoTopNode(gGeoManager, top->FindNode("muons_barrel_1"));
+   muon->SetVisLevel(4);
+   gEve->AddGlobalElement(muon);
+   muon = new TEveGeoTopNode(gGeoManager, top->FindNode("muons_endcap_1"));
+   muon->SetVisLevel(4);
+   muon->UseNodeTrans();
+   gEve->AddGlobalElement(muon);
+   muon = new TEveGeoTopNode(gGeoManager, top->FindNode("muons_endcap_2"));
+   muon->SetVisLevel(4);
+   muon->UseNodeTrans();
+   gEve->AddGlobalElement(muon);
+
+*/
+
+
    //TGeoVolume* top = gGeoManager->GetTopVolume()->FindNode("Delphes3DGeometry_1")->GetVolume();
-   TGeoVolume* top = det3D.getDetector(true);
+   //TGeoVolume* top = det3D.getDetector(true);
+   TGeoVolume* top = det3D.getDetector(false); //TODO: can this be set in the GUI?
+   geom->SetTopVolume(top);
    TEveElementList *geometry = new TEveElementList("Geometry");
    TEveGeoTopNode* trk = new TEveGeoTopNode(gGeoManager, top->FindNode("tracker_1"));
@@ -519 +569 @@
    muon->UseNodeTrans();
    geometry->AddElement(muon);
-   gGeoManager->DefaultColors();
+   //gGeoManager->DefaultColors();
 
    // Create chain of root trees
@@ -576 +626 @@
    TEveTrackPropagator *trkProp = gTrackList->GetPropagator();
    trkProp->SetMagField(0.0, 0.0, -gBz);
-   trkProp->SetMaxR(gRadius*100.0);
-   trkProp->SetMaxZ(gHalfLength*100.0);
+   trkProp->SetMaxR(gRadius);
+   trkProp->SetMaxZ(gHalfLength);
 
    // viewers and scenes
 
    TEveCalo3D *calo3d = new TEveCalo3D(gCaloData);
-   calo3d->SetBarrelRadius(gRadius*100.0);
-   calo3d->SetEndCapPos(gHalfLength*100.0);
+   calo3d->SetBarrelRadius(gRadius);
+   calo3d->SetEndCapPos(gHalfLength);
 
    //gStyle->SetPalette(1, 0);
@@ -699 +749 @@
     eveJetCone->SetMainTransparency(60);
     eveJetCone->SetLineColor(kYellow);
-    eveJetCone->SetCylinder(gRadius*100.0 - 10, gHalfLength*100.0 - 10);
+    eveJetCone->SetCylinder(gRadius - 10, gHalfLength - 10);
     eveJetCone->SetPickable(kTRUE);
     eveJetCone->AddEllipticCone(jet->Eta, jet->Phi, jet->DeltaEta, jet->DeltaPhi);
@@ -793 +843 @@
    gSystem->Load("libGeom");
    gSystem->Load("../libDelphes");
+
    delphes_event_display("delphes_card_CMS.tcl", "../delphes_output.root"); //TODO propagate parameters
-   return;
+   make_gui();
+   load_event();
+   gEve->Redraw3D(kTRUE); // Reset camera after the first event has been shown.
+
+   // EClipType not exported to CINT (see TGLUtil.h):
+   // 0 - no clip, 1 - clip plane, 2 - clip box
+   TGLViewer *v = gEve->GetDefaultGLViewer();
+   Double_t plane[4] = { 0., 1., 0., 0. };
+   //v->GetClipSet()->SetClipState(1,plane);
+   //v->GetClipSet()->SetClipType(1);
+   //v->ColorSet().Background().SetColor(kMagenta+4);
+   //v->SetGuideState(TGLUtil::kAxesEdge, kTRUE, kFALSE, 0);
+   v->RefreshPadEditor(v);
+   v->CurrentCamera().RotateRad(-1.2, 0.5);
+   v->DoDraw();
+
 /*
    TGeoManager *geom = new TGeoManager("delphes", "Delphes geometry");
@@ -863 +929 @@
    v->CurrentCamera().RotateRad(-1.2, 0.5);
    v->DoDraw();
-
-   make_gui();
-   // load_event();
-   // gEve->Redraw3D(kTRUE); // Reset camera after the first event has been shown.
-   */
-}
-
+*/
+}
+