#include <TClonesArray.h>
#include <TTree.h>
#include <TVirtualMC.h>
+#include <TGeoManager.h>
#include "AliACORDE.h"
-#include "AliACORDEModule.h"
#include "AliMagF.h"
#include "AliRun.h"
+#include "AliACORDERawData.h"
ClassImp(AliACORDE)
//_____________________________________________________________________________
AliACORDE::AliACORDE()
: AliDetector(),
- fModule(0)
+ fCreateCavern(0),
+ fITSGeometry(0)
{
//
// Default constructor
//_____________________________________________________________________________
AliACORDE::AliACORDE(const char *name, const char *title)
: AliDetector(name, title),
- fModule(0)
+ fCreateCavern(kFALSE),
+ fITSGeometry(kTRUE)
+
{
//
// Standard constructor
- //
- //fHits = new TClonesArray("AliACORDEhit", 400);
- //gAlice->GetMCApp()->AddHitList(fHits);
}
//_____________________________________________________________________________
//
// Default destructor
//
- if ( fModule ) { delete fModule; fModule = 0; }
}
//_____________________________________________________________________________
//
// Aluminum
AliMaterial(9, "ALUMINIUM0$", 26.98, 13., 2.7, 8.9, 37.2);
- AliMaterial(29, "ALUMINIUM1$", 26.98, 13., 2.7, 8.9, 37.2);
- AliMaterial(49, "ALUMINIUM2$", 26.98, 13., 2.7, 8.9, 37.2);
+ // AliMaterial(29, "ALUMINIUM1$", 26.98, 13., 2.7, 8.9, 37.2);
+ //AliMaterial(49, "ALUMINIUM2$", 26.98, 13., 2.7, 8.9, 37.2);
//
// Iron
- AliMaterial(10, "IRON0$ ", 55.85, 26., 7.87, 1.76, 17.1);
- AliMaterial(30, "IRON1$ ", 55.85, 26., 7.87, 1.76, 17.1);
- AliMaterial(50, "IRON2$ ", 55.85, 26., 7.87, 1.76, 17.1);
+ //AliMaterial(10, "IRON0$ ", 55.85, 26., 7.87, 1.76, 17.1);
+ //AliMaterial(30, "IRON1$ ", 55.85, 26., 7.87, 1.76, 17.1);
+ //AliMaterial(50, "IRON2$ ", 55.85, 26., 7.87, 1.76, 17.1);
//
// Air
- AliMaterial(15, "AIR0$ ", 14.61, 7.3, .001205, 30423.24, 67500.);
- AliMaterial(35, "AIR1$ ", 14.61, 7.3, .001205, 30423.24, 67500.);
- AliMaterial(55, "AIR2$ ", 14.61, 7.3, .001205, 30423.24, 67500.);
- AliMaterial(75, "AIR3$ ", 14.61, 7.3, .001205, 30423.24, 67500.);
- AliMaterial(95, "AIR4$ ", 14.61, 7.3, .001205, 30423.24, 67500.);
+ Float_t as[] = { 12.0107, 14.0067, 15.9994, 39.948 };
+ Float_t zs[] = { 6., 7., 8., 18. };
+ Float_t ws[] = { 0.000124, 0.755267, 0.231781, 0.012827 };
+ Double_t density = .00120479;
+ AliMixture(15, "AIR0$", as, zs, density, 4, ws);
+
+ //AliMaterial(15, "AIR0$ ", 14.61, 7.3, .001205, 30423.24, 67500.);
+ //AliMaterial(35, "AIR1$ ", 14.61, 7.3, .001205, 30423.24, 67500.);
+ //AliMaterial(55, "AIR2$ ", 14.61, 7.3, .001205, 30423.24, 67500.);
+ //AliMaterial(75, "AIR3$ ", 14.61, 7.3, .001205, 30423.24, 67500.);
+ //AliMaterial(95, "AIR4$ ", 14.61, 7.3, .001205, 30423.24, 67500.);
// Scintillator material polystyrene
//
// Aluminum
AliMedium(9, "ALU_C0 ", 9, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
- AliMedium(29, "ALU_C1 ", 29, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
- AliMedium(49, "ALU_C2 ", 49, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+ // AliMedium(29, "ALU_C1 ", 29, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+ //AliMedium(49, "ALU_C2 ", 49, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
//
// Iron
- AliMedium(10, "FE_C0 ", 10, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
- AliMedium(30, "FE_C1 ", 30, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
- AliMedium(50, "FE_C2 ", 50, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+ //AliMedium(10, "FE_C0 ", 10, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+ //AliMedium(30, "FE_C1 ", 30, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+ //AliMedium(50, "FE_C2 ", 50, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
//
// Air
AliMedium(15, "AIR_C0 ", 15, 0, isxfld, sxmgmx, atmaxfd, stemax, adeemax, aepsil, astmin);
- AliMedium(35, "AIR_C1 ", 35, 0, isxfld, sxmgmx, atmaxfd, stemax, adeemax, aepsil, astmin);
- AliMedium(55, "AIR_C2 ", 55, 0, isxfld, sxmgmx, atmaxfd, stemax, adeemax, aepsil, astmin);
- AliMedium(75, "AIR_C4 ", 75, 0, isxfld, sxmgmx, atmaxfd, stemax, adeemax, aepsil, astmin);
- AliMedium(95, "AIR_C5 ", 95, 0, isxfld, sxmgmx, atmaxfd, stemax, adeemax, aepsil, astmin);
+ // AliMedium(35, "AIR_C1 ", 35, 0, isxfld, sxmgmx, atmaxfd, stemax, adeemax, aepsil, astmin);
+ //AliMedium(55, "AIR_C2 ", 55, 0, isxfld, sxmgmx, atmaxfd, stemax, adeemax, aepsil, astmin);
+ //AliMedium(75, "AIR_C4 ", 75, 0, isxfld, sxmgmx, atmaxfd, stemax, adeemax, aepsil, astmin);
+ //AliMedium(95, "AIR_C5 ", 95, 0, isxfld, sxmgmx, atmaxfd, stemax, adeemax, aepsil, astmin);
// The scintillator of the CPV made of Polystyrene
// scintillator -> idtmed[1112]
- AliMedium(12 , "CPV scint.0 ", 13, 1, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);
+ //AliMedium(12 , "CPV scint.0 ", 13, 1, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);
AliMedium(13 , "CPV scint.1 ", 13, 1, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);
- AliMedium(14 , "CPV scint.2 ", 13, 1, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);
+ //AliMedium(14 , "CPV scint.2 ", 13, 1, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);
// Molasse -> idtmed[1123]
AliMedium(24 , "Molasse ", 24, 0, xfield, xfieldm, tmaxfd, stemax, deemax, xepsil, stmin);
AliMixture(17, "CONCRETE$", aconc, zconc, 2.35, 10, wconc);
// Concrete
AliMedium(17, "CC_C0 ", 17, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
- AliMedium(27, "CC_C1 ", 17, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin); // MX24
- AliMedium(37, "CC_C2 ", 17, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin); // PM25
- AliMedium(47, "CC_C3 ", 17, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin); // PGC2
+ // AliMedium(27, "CC_C1 ", 17, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin); // MX24
+ //AliMedium(37, "CC_C2 ", 17, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin); // PM25
+ //AliMedium(47, "CC_C3 ", 17, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin); // PGC2
}
}
AliDetector::MakeBranch(opt);
}
+
+AliLoader* AliACORDE::MakeLoader(const char* topfoldername)
+{
+
+ AliDebug(1,Form("Creating AliACORDELoader, Top folder is %s ",
+ topfoldername));
+ fLoader = new AliACORDELoader(GetName(),topfoldername);
+ return fLoader;
+}
+
+
+AliDigitizer* AliACORDE::CreateDigitizer(AliRunDigitizer* manager) const
+{
+ //
+ //
+ return new AliACORDEDigitizer(manager);
+}
+
+void AliACORDE::Digits2Raw()
+{
+ // Produce Raw data starting from digits
+ // 1. Get digits
+ // 2. From digits get an array with the state of the modules
+ // 3. Unload digits
+ // 4. Write raw data
+
+ // 1. Get digits
+
+ // 1.1 Get detector, load digits and set branch
+ AliACORDE* acorde = (AliACORDE*)gAlice->GetDetector("ACORDE");
+ fLoader->LoadDigits("READ");
+ TTree* treeD = fLoader->TreeD();
+ if (!treeD) {
+ Error("Digits2Raw", "no digits tree");
+ return;
+ }
+ TClonesArray *adigits = new TClonesArray ("AliACORDEdigit", 1000);
+ treeD->GetBranch("ACORDEdigit")->SetAddress(&adigits);
+ // 1.2 Get first entry (there is always only one)
+ acorde->ResetDigits();
+ treeD->GetEvent(0);
+
+ // 2. From digits get an array with the state of the modules
+ // 2.1 Define and initialize the array
+ Bool_t Modules[60];
+ for (Int_t i=0;i<60;i++) Modules[i]= kFALSE;
+ // 2.2 Loop over all digits
+ Int_t ndig = adigits->GetEntriesFast();
+ for (Int_t idig=0;idig<ndig;idig++) {
+ // 2.3 set the array entry for each digit
+ AliACORDEdigit* digit = (AliACORDEdigit*) adigits->At(idig);
+ Int_t mod = digit->GetModule();
+ Modules[mod-1]=kTRUE;
+ }
+ // 3. Unload digits
+ fLoader->UnloadDigits();
+
+ // 4. Write raw data
+ AliACORDERawData rawdata;
+ rawdata.WriteACORDERawData(Modules);
+}
+
+//_____________________________________________________________________________
+void AliACORDE::AddAlignableVolumes() const
+{
+ //
+ // Create entries for alignable volumes associating the symbolic volume
+ // name with the corresponding volume path. Needs to be syncronized with
+ // eventual changes in the geometry.
+ //
+ TString vpstr1 = "ALIC_1/ACORDE1_";
+ TString snstr1="ACORDE/Array";
+ TString volpath, symname;
+
+ for(Int_t cnt=1; cnt<=3; cnt++){
+ volpath = vpstr1;
+ volpath += cnt;
+ symname = snstr1;
+ symname += cnt;
+ if(!gGeoManager->SetAlignableEntry(symname.Data(),volpath.Data()))
+ AliFatal(Form("Alignable entry %s not created. Volume path %s not valid", symname.Data(),volpath.Data()));
+ }
+
+}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff