sPHENIX discussion of software

[Abhisek Sen <sen.abhisek AT gmail.com>]

Hello,

   Thanks for the comments. I didnt had the updated geometry, thats why both inner and outer Hcal were rotated in the same direction. Below is the updated plot.

I have updated the code. Compressed geometry file is just 1.5MB which can be read from my phenix publish area (http://www.phenix.bnl.gov/WWW/publish/abhisek/sphenix/sphenix_geom.root). So only thing you need is the macro. Just download the attached macro and run. It knows where to get the geometry file. Thats it !!

 To select a subset of the geometry, click on Scenes => Geometry scene. You can change color too. To see the the scintillators in side HCal for example, click on "Hcal_envelop", then check "VizVolumeDaughters".

 Have fun!

Cheers,

Abhisek

​

On Thu, Oct 8, 2015 at 3:49 PM, Edouard Kistenev <kistenev AT bnl.gov> wrote:

No, it is just in the west, everything is fine in the east part of it.

On Oct 8, 2015, at 2:53 PM, Huang, Jin <jhuang AT bnl.gov> wrote:

Wow, that is a very useful too. Thanks Abhisek!

 

Achim: on the tile angle, is that an artifact of rendering fine structure in screen? Or the geometry source to produce the screen shot is out of date?

 

Just to double check of the official setting, attached is the scintillator energy density for the latest pre-CDR Hijing production. The sign of the tile angle between inner and outer HCal are correctly assigned. 

 

Cheers

 

Jin

 

 

______________________________

 

Jin HUANG

 

Brookhaven National Laboratory

Physics Department, Bldg 510 C

Upton, NY 11973-5000

 

Office: 631-344-5898

Cell:   757-604-9946

______________________________

 

From: sphenix-software-l-bounces AT lists.bnl.gov [mailto:sphenix-software-l-bounces AT lists.bnl.gov] On Behalf Of Achim Franz
Sent: Thursday, October 8, 2015 2:20 PM
To: Abhisek Sen <sen.abhisek AT gmail.com>
Cc: sphenix-software-l AT lists.bnl.gov
Subject: Re: [Sphenix-software-l] Visualization

 

Hi, 

 

is that just me or are the inner and outer HCal rotated in the same direction ?? 

 

Achim

 

 

On Oct 08, 2015, at 13:44 , Abhisek Sen <sen.abhisek AT gmail.com> wrote:

 

Hello all,

     After the software meeting last Tuesday, I realized it would be better to have a nice and interactive visualization of our Geant4 setup. It will give us a better understanding on whats going on. The visualization we get running the "vis.mac" is very poor. So I converted sPHENIX Geant4 geometry to a TGeo root file via GDML(http://gdml.web.cern.ch/GDML/). This also opens up the possibility of storing the geometry in GDML/TGeo object in a later stage.

 

To get the interactive visualization you only have to do, (in RCF)

root.exe /direct/phenix+sim01/sPHENIX/sunywrk/abhisek/vis/sphenix_geom.C

 

It reads the geometry file /direct/phenix+sim01/sPHENIX/sunywrk/abhisek/vis/sphenix_geom.root . This does not depend on your setup, only needs root. You can just copy the macro and root file to your laptop to get better visualization. No need of sPHENIX environment. It looks better in my mac because of better graphics card. You can see some pictures below.

 

I will show more details in next software meeting.

Cheers,

Abhisek

 

 

<2.png>
<3.png>
​

 

 

 

 

_______________________________________________
Sphenix-software-l mailing list
Sphenix-software-l AT lists.bnl.gov
https://lists.bnl.gov/mailman/listinfo/sphenix-software-l

 

<EnergyDensity.png>_______________________________________________
Sphenix-software-l mailing list
Sphenix-software-l AT lists.bnl.gov
https://lists.bnl.gov/mailman/listinfo/sphenix-software-l

void sphenix_geom()
{
   TEveManager::Create();

   TFile::SetCacheFileDir(".");
   gGeoManager = TGeoManager::Import("http://www.phenix.bnl.gov/WWW/publish/abhisek/sphenix/sphenix_geom.root";);
   gGeoManager->DefaultColors();
   TGeoVolume* top = gGeoManager->GetTopVolume();
   //Print the list of nodes
   //top->PrintNodes();

   int ngeom = 29;
   //Full list of geometries
   string subsystems[] = {"PIPE_0_25501",
			  "SVTX_0_25502",
			  "SVTXSUPPORT_0_25503",
			  "SVTX_1_25504",
			  "SVTXSUPPORT_1_25505",
			  "SVTX_2_25506",
			  "SVTXSUPPORT_2_25507",
			  "SVTX_3_25508",
			  "SVTXSUPPORT_3_25509",
			  "SVTX_4_25510",
			  "SVTXSUPPORT_4_25511",
			  "SVTX_5_25512",
			  "SVTXSUPPORT_5_25513",
			  "SVTX_6_25514",
			  "SVTXSUPPORT_6_25515",
			  "ABSORBER_PRESHOWER_1_25516",
			  "PRESHOWER_1_25517",
			  "CEMC_ELECTRONICS_0_25518",
			  "CEMC_1_25519",
			  "CEMC_SPT_0_25520",
			  "Hcal_envelope_25521",
			  "HCALIN_SPT_N1_0_25522",
			  "HCALIN_SPT_N1_1_25523",
			  "HCALIN_SPT_N1_2_25524",
			  "HCALIN_SPT_N1_3_25525",
			  "MAGNET_0_25526",
			  "MAGNET_1_25527",
			  "MAGNET_2_25528",
			  "OuterHcal_envelope_25529",
			  "BH_1_25530",
			  "BH_FORWARD_PLUS_1_25531",
			  "BH_FORWARD_NEG_1_25532" };


   for(int igeom=0; igeom<ngeom; igeom++)
   {
    cout << "Reading geometry " << subsystems[igeom] << endl;
    TGeoNode *geo_node = top->FindNode(subsystems[igeom].c_str());
    geo_node->GetVolume()->VisibleDaughters(kFALSE);
    TEveGeoTopNode* eve_node = new TEveGeoTopNode(gGeoManager,geo_node);
    eve_node->SetVisLevel(6);
    gEve->AddGlobalElement(eve_node);
   }

   gGeoManager->CloseGeometry();
   gEve->FullRedraw3D(kTRUE);

   // EClipType not exported to CINT (see TGLUtil.h):
   // 0 - no clip, 1 - clip plane, 2 - clip box
   TGLViewer *v = gEve->GetDefaultGLViewer();
   v->GetClipSet()->SetClipType(1);
   v->ColorSet().Background().SetColor(kMagenta+4);
   v->SetGuideState(TGLUtil::kAxesEdge, kTRUE, kFALSE, 0);
   v->RefreshPadEditor(v);
   v->CurrentCamera().RotateRad(-0.5, 0.5);
   v->DoDraw();
}