Minor changing for TB setup possibility
authorpavlinov <pavlinov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 28 Nov 2010 16:27:08 +0000 (16:27 +0000)
committerpavlinov <pavlinov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 28 Nov 2010 16:27:08 +0000 (16:27 +0000)
EMCAL/AliEMCAL.cxx
EMCAL/AliEMCALWsuCosmicRaySetUp.cxx
EMCAL/AliEMCALWsuCosmicRaySetUp.h
EMCAL/AliEMCALv0.cxx

index 93fcaeb..7c05275 100644 (file)
@@ -223,6 +223,9 @@ void AliEMCAL::CreateMaterials()
   AliMixture(5, "BondPaper$", apaper, zpaper, 0.75, 3, wpaper);
  
   // DEFINITION OF THE TRACKING MEDIA
+  // Look to the $ALICE_ROOT/data/galice.cuts for particular values
+  // of cuts.
+  // Don't forget to add a new tracking medium with non-default cuts
 
   // for EMCAL: idtmed[1599->1698] equivalent to fIdtmed[0->100]
   Int_t   isxfld = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Integ() ;
@@ -250,10 +253,10 @@ void AliEMCAL::CreateMaterials()
   AliMedium(4, "S steel$", 4, 0, 
              isxfld, sxmgmx, 10.0, 0.01, 0.1, 0.001, 0.001, 0, 0) ;
 
-  // Oct 26,2010                                                                   -> idtmed[1604]
-  Float_t stemax = 0.001;
+  // Oct 26,2010; Nov 24,2010                                                      -> idtmed[1604]
+  deemax = 0.01;
   AliMedium(5, "Paper$", 5, 0, 
-             isxfld, sxmgmx, 10.0, stemax, 0.1, 0.001, 0.001, 0, 0) ;
+             isxfld, sxmgmx, 10.0, deemax, 0.1, 0.001, 0.001, 0, 0) ;
 
 
   //set constants for Birk's Law implentation
index 5137c62..0093884 100644 (file)
@@ -77,15 +77,28 @@ void AliEMCALWsuCosmicRaySetUp::CreateGeometry()
   */
   //End_Html
 
-  Float_t dASUC[3];
+  // Master Volume
+  fMasterVolume[0] = fMasterVolume[1] = 25.0;
+  fMasterVolume[2] = 300.;
+
   Int_t *idtmed = fIdtmed->GetArray()+1;
-  int idSC = idtmed[0];
+  int idAir = idtmed[0];
+  gMC->Gsvolu(GetName(),"BOX",idAir, fMasterVolume,3); // Master volume
+  //
+  // Sc counters
   //
-  dASUC[0]=50;
-  dASUC[1]=50;
-  dASUC[2]=50;
-  //  TString tmp(GetTitle());
-  gMC->Gsvolu(GetName(),"BOX",idSC, dASUC,3); // WSUC - Wsu Cosmic Ray SetUp
+  Float_t sc[3]; // tube
+  sc[0] = 0.0;
+  sc[1] = 5.0;
+  sc[2] = 0.5;
+  Float_t zsc[3]={10.,110., 310.};
+  int idSC = idtmed[1];
+  gMC->Gsvolu("SCOU","TUBE",idSC, sc,3); // Master volume
+  Int_t idRot=0; // no rotation
+  for(Int_t i=0; i<3; i++) {
+    Float_t zpos = zsc[i] - fMasterVolume[2];
+    gMC->Gspos("SCOU", i+1, "WSUC", 0.0, 0.0, zpos, idRot, "ONLY"); 
+  }
 }
  
 //_____________________________________________________________________________
@@ -93,18 +106,26 @@ void AliEMCALWsuCosmicRaySetUp::CreateMaterials()
 {
 // Create materials and media
   Int_t   isxfld = 0;
-  Float_t sxmgmx = 0.;
-  
+  Float_t sxmgmx = 0., deemax = 0.1;  
   // AIR
   Float_t aAir[4]={12.0107,14.0067,15.9994,39.948};
   Float_t zAir[4]={6.,7.,8.,18.};
   Float_t wAir[4]={0.000124,0.755267,0.231781,0.012827};
   Float_t dAir = 1.20479E-3;
-  //  Float_t dAir1 = 1.20479E-10;
-  //
   AliMixture(1,"Air     $",aAir,zAir,dAir,4,wAir);
+
+  // --- The polysterene scintillator (CH) ---
+  Float_t aP[2] = {12.011, 1.00794} ;
+  Float_t zP[2] = {6.0, 1.0} ;
+  Float_t wP[2] = {1.0, 1.0} ;
+  Float_t dP = 1.032 ;
+  AliMixture(2, "Polystyrene$", aP, zP, dP, -2, wP) ;
+
   //
   AliMedium(1,"Air     $",1,0,isxfld,sxmgmx,10,-1,-0.1,0.1 ,-10);
+  AliMedium(2, "Scintillator$", 2, 1,
+            isxfld, sxmgmx, 10.0, 0.001, deemax, 0.001, 0.001, 0, 0) ;
+  
 }
  
 //_____________________________________________________________________________
index b3fb8ce..77728dd 100644 (file)
@@ -22,6 +22,11 @@ public:
   virtual void  CreateMaterials();
   virtual Int_t IsVersion() const {return 0;}
   void  DrawWSUC(float cxy=0.025) const; // *MENU*
+  // GetMethod
+  Float_t* GetMasterVolume() {return fMasterVolume;}
+
+  protected:
+  Float_t fMasterVolume[3]; // size of MASTER volume
 
   ClassDef(AliEMCALWsuCosmicRaySetUp,1)  // Class manager for the Wsu Cosmic Ray SetUp
 };
index 1aca666..5a831e9 100644 (file)
@@ -126,15 +126,16 @@ void AliEMCALv0::CreateGeometry()
   
   Float_t envelopA[10];
   if(gn.Contains("WSUC") ) { // TRD1 for WSUC facility
-    // 17-may-05 - just BOX
-    envelopA[0] = 26;
-    envelopA[1] = 15;
-    envelopA[2] = 30;
+    // Nov 25,2010
+    envelopA[0] = 30.;
+    envelopA[1] = 30;
+    envelopA[2] = 20;
     gMC->Gsvolu("XEN1", "BOX", fIdTmedArr[kIdSC], envelopA, 3) ;
     fEnvelop1.Set(3);
     for(int i=0; i<3; i++) fEnvelop1[i] = envelopA[i]; // 23-may-05  
-    // Position the EMCAL Mother Volume (XEN1) in WSUC  
-    gMC->Gspos("XEN1", 1, "WSUC", 0.0, 0.0, 0.0, fIdRotm, "ONLY") ;
+    // Position the EMCAL Mother Volume (XEN1) in WSUC.
+    // Look to AliEMCALWsuCosmicRaySetUp.  
+    gMC->Gspos("XEN1", 1, "WSUC", 0.0, 0.0, + 265., fIdRotm, "ONLY") ;
   } else { 
     envelopA[0] = geom->GetArm1PhiMin();                         // minimum phi angle
     envelopA[1] = geom->GetArm1PhiMax() - geom->GetArm1PhiMin(); // angular range in phi