///////////////////////////////////////////////////////////////////////////////
#include <TClonesArray.h>
+#include <TGeoGlobalMagField.h>
+#include <TGeoManager.h>
+#include <TStopwatch.h>
#include <TTree.h>
#include <TVirtualMC.h>
-#include <TGeoManager.h>
#include "AliACORDE.h"
#include "AliMagF.h"
#include "AliRun.h"
#include "AliACORDERawData.h"
+#include "AliACORDERawStream.h"
ClassImp(AliACORDE)
AliACORDE::AliACORDE()
: AliDetector(),
fCreateCavern(0),
- fITSGeometry(0)
+ f4CentralModulesGeometry(0)
{
//
// Default constructor
AliACORDE::AliACORDE(const char *name, const char *title)
: AliDetector(name, title),
fCreateCavern(kFALSE),
- fITSGeometry(kTRUE)
+ f4CentralModulesGeometry(kTRUE)
{
//
void AliACORDE::CreateMaterials()
{
// Magnatic field inside the pit
- Int_t isxfld = gAlice->Field()->Integ();
- Float_t sxmgmx = gAlice->Field()->Max();
+ Int_t isxfld = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Integ();
+ Float_t sxmgmx = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Max();
//Magnetic field above the Magnet.
Int_t xfield = 0; // no Magnetic field.
//
// 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
Float_t zMolasse[10] = {1., 6., 8., 11., 12., 13., 14., 19., 20., 26.};
Float_t wMolasse[10] = {0.008, 0.043, 0.485, 0.007, 0.042, 0.037, 0.215, 0.023, 0.1, 0.04};
Float_t dMolasse = 2.40;
- AliMixture(24, "Molasse$", aMolasse, zMolasse, dMolasse, 10, wMolasse);
+ AliMixture(24, "Molasse$", aMolasse, zMolasse, 2*dMolasse, 10, wMolasse); // correction to density of molasse alpine newDen=2*dmolasse
// ****************
// Defines tracking media parameters.
//
// 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
}
TBranch *branch;
char branchname[20];
- sprintf(branchname,"%s",GetName());
+ snprintf(branchname,19,"%s",GetName());
// Branch address for hit tree
TTree *treeH = fLoader->TreeH();
if (treeH ) {
}
-AliDigitizer* AliACORDE::CreateDigitizer(AliRunDigitizer* manager) const
+AliDigitizer* AliACORDE::CreateDigitizer(AliDigitizationInput* digInput) const
{
//
//
- return new AliACORDEDigitizer(manager);
+ return new AliACORDEDigitizer(digInput);
}
void AliACORDE::Digits2Raw()
// 2.3 set the array entry for each digit
AliACORDEdigit* digit = (AliACORDEdigit*) adigits->At(idig);
Int_t mod = digit->GetModule();
- Modules[mod-1]=kTRUE;
+ Modules[mod]=kTRUE;
}
// 3. Unload digits
fLoader->UnloadDigits();
// 4. Write raw data
AliACORDERawData rawdata;
- rawdata.WriteACORDERawData(Modules);
+ rawdata.WriteACORDERawData(Modules,(ndig > 1));
}
//_____________________________________________________________________________
-void AliACORDE::AddAlignableVolumes() const
+Bool_t AliACORDE::Raw2SDigits(AliRawReader* rawReader)
{
//
- // 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()));
+ // Reads the raw data stream and exracts the digits
+ //
+ // Input:
+ // rawReader : pointer to the current AliRawReader
+ // Output:
+ //
+ // Created: 31 Jan 2008 Mario Sitta
+ //
+ TStopwatch timer;
+ timer.Start();
+
+ if(!fLoader) {
+ AliError("no ACORDE loader found");
+ return kFALSE;
+ }
+
+ TTree* treeD = fLoader->TreeD();
+ if(!treeD) {
+ fLoader->MakeTree("D");
+ treeD = fLoader->TreeD();
+ }
+
+ AliACORDEdigit digit;
+ AliACORDEdigit* pdigit = &digit;
+ const Int_t kBufferSize = 4000;
+
+ treeD->Branch("ACORDE", "AliACORDEdigit", &pdigit, kBufferSize);
+
+ // rawReader->Reset();
+ AliACORDERawStream* rawStream = new AliACORDERawStream(rawReader);
+
+ if (!rawStream->Next()) return kFALSE; // No ACORDE data found
+ /*
+ for(Int_t i=0; i<64; i++) {
+ new(pdigit) AliACORDEdigit(i, (Int_t)rawStream->GetADC(i), (Int_t)rawStream->GetTime(i));
+ treeD->Fill();
}
+ */
+ fLoader->WriteDigits("OVERWRITE");
+ fLoader->UnloadDigits();
+
+ delete rawStream;
+ timer.Stop();
+ timer.Print();
+
+ return kTRUE;
}
+
+//_____________________________________________________________________________
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff