* provided "as is" without express or implied warranty. *
**************************************************************************/
-
+/* $Id$ */
///////////////////////////////////////////////////////////////////////////
// //
// V-Zero Detector //
// This class contains the base procedures for the VZERO detector //
+// Default geometry of November 2003 : V0R box is 4.4 cm thick //
+// scintillators are 2 cm thick //
// All comments should be sent to Brigitte CHEYNIS : //
// b.cheynis@ipnl.in2p3.fr //
// //
// //
///////////////////////////////////////////////////////////////////////////
+
+// --- Standard libraries ---
#include <Riostream.h>
-#include <Riostream.h>
+#include <stdlib.h>
+
+// --- ROOT libraries ---
+#include <TNamed.h>
+#include "TROOT.h"
+#include "TFile.h"
+#include "TNetFile.h"
+#include "TRandom.h"
+#include "TTree.h"
+#include "TBranch.h"
+#include "TClonesArray.h"
+#include "TStopwatch.h"
+#include "TParameter.h"
-#include "TMath.h"
-#include "TTUBE.h"
-#include "TNode.h"
-#include "TGeometry.h"
+// --- AliRoot header files ---
#include "AliRun.h"
+#include "AliMC.h"
#include "AliVZERO.h"
+#include "AliVZEROLoader.h"
+#include "AliVZERODigitizer.h"
+#include "AliVZEROBuffer.h"
+#include "AliRunDigitizer.h"
#include "AliVZEROdigit.h"
-#include "AliMC.h"
-#include "AliVZEROhit.h"
+#include "AliDAQ.h"
+#include "AliRawReader.h"
+#include "AliVZERORawStream.h"
+#include "AliVZEROCalibData.h"
ClassImp(AliVZERO)
-
-
+ //__________________________________________________________________
+AliVZERO::AliVZERO(): AliDetector(),
+ fIdSens1(0),
+ fThickness(0.),
+ fThickness1(0.),
+ fMaxStepQua(0.),
+ fMaxStepAlu(0.),
+ fMaxDestepQua(0.),
+ fMaxDestepAlu(0.)
+{
+/// Default Constructor
+
+ AliDebug(1,Form("default (empty) ctor this=%p",this));
+ fIshunt = 0;
+}
//_____________________________________________________________________________
AliVZERO::AliVZERO(const char *name, const char *title)
- : AliDetector(name,title)
+ : AliDetector(name,title),
+ fIdSens1(0),
+ fThickness(4.4),
+ fThickness1(2.0),
+ fMaxStepQua(0.05),
+ fMaxStepAlu(0.01),
+ fMaxDestepQua(-1.0),
+ fMaxDestepAlu(-1.0)
{
- //
+
// Standard constructor for VZERO Detector
- //
- fIshunt = 1; // All hits are associated with primary particles
+ AliDebug(1,Form("ctor this=%p",this));
+
+ // fIshunt = 1; // All hits are associated with primary particles
fHits = new TClonesArray("AliVZEROhit", 400);
fDigits = new TClonesArray("AliVZEROdigit",400);
- gAlice->AddHitList(fHits);
+ gAlice->GetMCApp()->AddHitList(fHits);
-// fDigits = new TClonesArray("AliVZEROdigit",400) ;
-
- fThickness = 3.1; // total thickness of the V0R box
- fThickness1 = 1.0; // thickness of the thickest cell (2.5 in version 0)
-
- fMaxStepQua = 0.05;
- fMaxStepAlu = 0.01;
+// fThickness = 4.4; // total thickness of the V0R box in cm
+// fThickness1 = 2.0; // thickness of scintillating cells in cm
+//
+// fMaxStepQua = 0.05;
+// fMaxStepAlu = 0.01;
+//
+// fMaxDestepQua = -1.0;
+// fMaxDestepAlu = -1.0;
- fMaxDestepQua = -1.0;
- fMaxDestepAlu = -1.0;
- SetMarkerColor(kRed);
}
//_____________________________________________________________________________
AliVZERO::~AliVZERO()
{
+ //
+ // Default destructor for VZERO Detector
+ //
+
if (fHits) {
fHits->Delete();
delete fHits;
- }
+ fHits=0; }
+
+ if (fDigits) {
+ fDigits->Delete();
+ delete fDigits;
+ fDigits=0; }
}
-//_____________________________________________________________________________
-void AliVZERO::BuildGeometry()
-{
- //
- // Build simple ROOT TNode geometry for event display
- //
-}
-
//_____________________________________________________________________________
void AliVZERO::CreateGeometry()
{
//
- // Build simple ROOT TNode geometry for event display
+ // Builds simple Geant3 geometry
//
}
//_____________________________________________________________________________
void AliVZERO::CreateMaterials()
{
//
- // Build simple ROOT TNode geometry for event display
+ // Creates materials used for Geant3 geometry
//
}
-
-
//_____________________________________________________________________________
-Int_t AliVZERO::DistanceToPrimitive(Int_t px, Int_t py)
-{
- //
- // Calculate the distance from the mouse to the VZERO on the screen
- // Dummy routine
- //
-
- return 9999;
-}
-
-//-------------------------------------------------------------------------
void AliVZERO::Init()
{
//
- // Initialise the VZERO after it has been built
+ // Initialises the VZERO class after it has been built
//
}
-//-------------------------------------------------------------------------
-
+//_____________________________________________________________________________
void AliVZERO::SetMaxStepQua(Float_t p1)
{
+ //
+ // Possible parametrisation of steps in active materials
+ //
fMaxStepQua = p1;
}
-//___________________________________________
+//_____________________________________________________________________________
void AliVZERO::SetMaxStepAlu(Float_t p1)
{
+ //
+ // Possible parametrisation of steps in Aluminum foils (not used in
+ // version v2)
+ //
fMaxStepAlu = p1;
}
-//___________________________________________
+//_____________________________________________________________________________
void AliVZERO::SetMaxDestepQua(Float_t p1)
{
+ //
+ // Possible parametrisation of steps in active materials (quartz)
+ //
fMaxDestepQua = p1;
}
-//___________________________________________
+//_____________________________________________________________________________
void AliVZERO::SetMaxDestepAlu(Float_t p1)
{
+ //
+ // Possible parametrisation of steps in Aluminum (not used in
+ // version v2)
+ //
fMaxDestepAlu = p1;
}
+
+//_____________________________________________________________________________
+AliLoader* AliVZERO::MakeLoader(const char* topfoldername)
+{
+ //
+ // Builds VZEROgetter (AliLoader type)
+ // if detector wants to use customized getter, it must overload this method
+ //
+// Info("MakeLoader","Creating AliVZEROLoader. Top folder is %s.",topfoldername);
+
+ AliDebug(1,Form("Creating AliVZEROLoader, Top folder is %s ",topfoldername));
+ fLoader = new AliVZEROLoader(GetName(),topfoldername);
+ return fLoader;
+}
+
+//_____________________________________________________________________________
+void AliVZERO::SetTreeAddress()
+{
+ //
+ // Sets tree address for hits.
+ //
+ if (fLoader->TreeH() && (fHits == 0x0))
+ fHits = new TClonesArray("AliVZEROhit", 400);
+
+ AliDetector::SetTreeAddress();
+}
+
+//_____________________________________________________________________________
+AliDigitizer* AliVZERO::CreateDigitizer(AliRunDigitizer* manager) const
+{
+ //
+ // Creates a digitizer for VZERO
+ //
+ return new AliVZERODigitizer(manager);
+}
+
+//_____________________________________________________________________________
+void AliVZERO::Hits2Digits(){
+ //
+ // Converts hits to digits of the current event
+ //
+ // Inputs file name
+ const char *alifile = "galice.root";
+
+ // Create the run digitizer
+ AliRunDigitizer* manager = new AliRunDigitizer(1, 1);
+ manager->SetInputStream(0, alifile);
+ manager->SetOutputFile("H2Dfile");
+
+ // Creates the VZERO digitizer
+ AliVZERODigitizer* dig = new AliVZERODigitizer(manager);
+
+ // Creates the digits
+ dig->Exec("");
+
+}
+//_____________________________________________________________________________
+void AliVZERO::Digits2Raw()
+{
+ //
+ // Converts digits of the current event to raw data
+ //
+
+ AliVZERO *fVZERO = (AliVZERO*)gAlice->GetDetector("VZERO");
+ fLoader->LoadDigits();
+ TTree* digits = fLoader->TreeD();
+ if (!digits) {
+ Error("Digits2Raw", "no digits tree");
+ return;
+ }
+ TClonesArray * VZEROdigits = new TClonesArray("AliVZEROdigit",1000);
+ fVZERO->SetTreeAddress();
+ digits->GetBranch("VZERODigit")->SetAddress(&VZEROdigits);
+
+ const char *fileName = AliDAQ::DdlFileName("VZERO",0);
+ AliVZEROBuffer* buffer = new AliVZEROBuffer(fileName);
+
+ // Get Trigger information first
+ // Read trigger inputs from trigger-detector object
+ AliDataLoader * dataLoader = fLoader->GetDigitsDataLoader();
+ if( !dataLoader->IsFileOpen() )
+ dataLoader->OpenFile( "READ" );
+ AliTriggerDetector* trgdet = (AliTriggerDetector*)dataLoader->GetDirectory()->Get( "Trigger" );
+ UInt_t triggerInfo = 0;
+ if(trgdet) {
+ triggerInfo = trgdet->GetMask() & 0xffff;
+ }
+ else {
+ AliError(Form("There is no trigger object for %s",fLoader->GetName()));
+ }
+
+ buffer->WriteTriggerInfo((UInt_t)triggerInfo);
+ buffer->WriteTriggerScalers();
+ buffer->WriteBunchNumbers();
+
+ Int_t aBBflagsV0A = 0;
+ Int_t aBBflagsV0C = 0;
+ Int_t aBGflagsV0A = 0;
+ Int_t aBGflagsV0C = 0;
+
+ if (digits->GetUserInfo()->FindObject("BBflagsV0A")) {
+ aBBflagsV0A = ((TParameter<int>*)digits->GetUserInfo()->FindObject("BBflagsV0A"))->GetVal();
+ }
+ else
+ AliWarning("V0A beam-beam flags not found in digits tree UserInfo! The flags will not be written to the raw-data stream!");
+
+ if (digits->GetUserInfo()->FindObject("BBflagsV0C")) {
+ aBBflagsV0C = ((TParameter<int>*)digits->GetUserInfo()->FindObject("BBflagsV0C"))->GetVal();
+ }
+ else
+ AliWarning("V0C beam-beam flags not found in digits tree UserInfo! The flags will not be written to the raw-data stream!");
+
+ if (digits->GetUserInfo()->FindObject("BGflagsV0A")) {
+ aBGflagsV0A = ((TParameter<int>*)digits->GetUserInfo()->FindObject("BGflagsV0A"))->GetVal();
+ }
+ else
+ AliWarning("V0A beam-gas flags not found in digits tree UserInfo! The flags will not be written to the raw-data stream!");
+
+ if (digits->GetUserInfo()->FindObject("BGflagsV0C")) {
+ aBGflagsV0C = ((TParameter<int>*)digits->GetUserInfo()->FindObject("BGflagsV0C"))->GetVal();
+ }
+ else
+ AliWarning("V0C beam-gas flags not found in digits tree UserInfo! The flags will not be written to the raw-data stream!");
+
+ // Now retrieve the channel information: charge smaples + time and
+ // dump it into ADC and Time arrays
+ Int_t nEntries = Int_t(digits->GetEntries());
+ Short_t aADC[64][AliVZEROdigit::kNClocks];
+ Float_t aTime[64];
+ Float_t aWidth[64];
+ Bool_t aIntegrator[64];
+ Bool_t aBBflag[64];
+ Bool_t aBGflag[64];
+
+ for (Int_t i = 0; i < nEntries; i++) {
+ fVZERO->ResetDigits();
+ digits->GetEvent(i);
+ Int_t ndig = VZEROdigits->GetEntriesFast();
+
+ if(ndig == 0) continue;
+ for(Int_t k=0; k<ndig; k++){
+ AliVZEROdigit* fVZERODigit = (AliVZEROdigit*) VZEROdigits->At(k);
+ // Convert aliroot channel k into FEE channel iChannel before writing data
+ Int_t iChannel = AliVZEROCalibData::GetBoardNumber(fVZERODigit->PMNumber()) * 8 +
+ AliVZEROCalibData::GetFEEChannelNumber(fVZERODigit->PMNumber());
+ for(Int_t iClock = 0; iClock < AliVZEROdigit::kNClocks; ++iClock) aADC[iChannel][iClock] = fVZERODigit->ChargeADC(iClock);
+ aTime[iChannel] = fVZERODigit->Time();
+ aWidth[iChannel] = fVZERODigit->Width();
+ aIntegrator[iChannel]= fVZERODigit->Integrator();
+ if(fVZERODigit->PMNumber() < 32) {
+ aBBflag[iChannel] = (aBBflagsV0C >> fVZERODigit->PMNumber()) & 0x1;
+ aBGflag[iChannel] = (aBGflagsV0C >> fVZERODigit->PMNumber()) & 0x1;
+ }
+ else {
+ aBBflag[iChannel] = (aBBflagsV0A >> (fVZERODigit->PMNumber()-32)) & 0x1;
+ aBGflag[iChannel] = (aBGflagsV0A >> (fVZERODigit->PMNumber()-32)) & 0x1;
+ }
+ AliDebug(1,Form("DDL: %s\tdigit number: %d\tPM number: %d\tADC: %f\tTime: %f",
+ fileName,k,fVZERODigit->PMNumber(),aADC[k],aTime[k]));
+ }
+ }
+
+ // Now fill raw data
+
+ for (Int_t iCIU = 0; iCIU < 8; iCIU++) {
+
+ // decoding of one Channel Interface Unit numbered iCIU - there are 8 channels per CIU (and 8 CIUs) :
+
+ for(Int_t iChannel_Offset = iCIU*8; iChannel_Offset < (iCIU*8)+8; iChannel_Offset=iChannel_Offset+4) {
+ for(Int_t iChannel = iChannel_Offset; iChannel < iChannel_Offset+4; iChannel++) {
+ buffer->WriteChannel(iChannel, aADC[iChannel], aIntegrator[iChannel]);
+ }
+ buffer->WriteBeamFlags(&aBBflag[iChannel_Offset],&aBGflag[iChannel_Offset]);
+ buffer->WriteMBInfo();
+ buffer->WriteMBFlags();
+ buffer->WriteBeamScalers();
+ }
+
+ for(Int_t iChannel = iCIU*8 + 7; iChannel >= iCIU*8; iChannel--) {
+ buffer->WriteTiming(aTime[iChannel], aWidth[iChannel]);
+ }
+
+ // End of decoding of one CIU card
+
+ } // end of decoding the eight CIUs
+
+ delete buffer;
+ fLoader->UnloadDigits();
+}
+
+//_____________________________________________________________________________
+Bool_t AliVZERO::Raw2SDigits(AliRawReader* rawReader){
+ // Converts the VZERO raw data into digits
+ // The method is used for merging simulated and
+ // real data events
+ TStopwatch timer;
+ timer.Start();
+
+ if(!fLoader) {
+ AliError("no VZERO loader found");
+ return kFALSE; }
+
+ TTree* treeD = fLoader->TreeD();
+ if(!treeD) {
+ fLoader->MakeTree("D");
+ treeD = fLoader->TreeD(); }
+
+ AliVZEROdigit digit;
+ AliVZEROdigit* pdigit = &digit;
+ const Int_t kBufferSize = 4000;
+
+ treeD->Branch("VZERO", "AliVZEROdigit", &pdigit, kBufferSize);
+
+ rawReader->Reset();
+ AliVZERORawStream* rawStream = new AliVZERORawStream(rawReader);
+
+ if (!rawStream->Next()) return kFALSE; // No VZERO data found
+
+ fLoader->WriteDigits("OVERWRITE");
+ fLoader->UnloadDigits();
+
+ delete rawStream;
+
+ timer.Stop();
+ timer.Print();
+ return kTRUE;
+}
+
+