macros/leptoquarks_fun4all/G4_Aerogel.C

0001 // $Id: G4_Aerogel.C,v 1.2 2013/10/09 01:08:17 jinhuang Exp $
0002
0003 /*!
0004  * \file G4_Aerogel.C
0005  * \brief Aerogel RICH for EIC detector
0006  * \author Jin Huang <jhuang@bnl.gov>
0007  * \version $Revision: 1.2 $
0008  * \date $Date: 2013/10/09 01:08:17 $
0009  */
0010
0011 void
0012 AerogelInit()
0013 {
0014
0015 }
0016
0017 void
0018 AerogelSetup(PHG4Reco* g4Reco, const int N_Sector = 8, //
0019     const double min_eta = 1.1 // 1.45
0020     )
0021 {
0022
0023   PHG4SectorSubsystem *ag;
0024   ag = new PHG4SectorSubsystem("Aerogel");
0025
0026   ag->get_geometry().set_normal_polar_angle(
0027       (PHG4Sector::Sector_Geometry::eta_to_polar_angle(min_eta)
0028           + PHG4Sector::Sector_Geometry::eta_to_polar_angle(2)) / 2);
0029 //  ag->get_geometry().set_normal_polar_angle(0);
0030   ag->get_geometry().set_normal_start(
0031       280 * PHG4Sector::Sector_Geometry::Unit_cm()); // 307
0032   ag->get_geometry().set_min_polar_angle(
0033       PHG4Sector::Sector_Geometry::eta_to_polar_angle(1.9));
0034   ag->get_geometry().set_max_polar_angle(
0035       PHG4Sector::Sector_Geometry::eta_to_polar_angle(min_eta));
0036   ag->get_geometry().set_min_polar_edge(
0037       PHG4Sector::Sector_Geometry::FlatEdge());
0038   ag->get_geometry().set_material("G4_AIR");
0039   ag->get_geometry().set_N_Sector(N_Sector);
0040   ag->OverlapCheck(overlapcheck);
0041
0042   // Aerogel dimensions ins cm
0043   double radiator_length = 2.;
0044   double expansion_length = 10.;
0045
0046   ag->get_geometry().AddLayers_AeroGel_ePHENIX( radiator_length * PHG4Sector::Sector_Geometry::Unit_cm(),
0047                         expansion_length * PHG4Sector::Sector_Geometry::Unit_cm() );
0048   g4Reco->registerSubsystem(ag);
0049
0050 }
0051