[u/mrichter/AliRoot.git] / VZERO / AliVZERO.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * 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 <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"

// --- 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 "AliDAQ.h"
#include "AliRawReader.h"
#include "AliVZERORawStream.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),
         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
  
  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->GetMCApp()->AddHitList(fHits);

//   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;
  
  
}

//_____________________________________________________________________________
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()
{
  //
  // Builds simple ROOT TNode geometry for event display
  //
}
 
//_____________________________________________________________________________
void AliVZERO::CreateGeometry()
{
  //
  // Builds simple Geant3 geometry 
  //
}
//_____________________________________________________________________________
void AliVZERO::CreateMaterials()
{
  //
  // Creates materials used for Geant3 geometry 
  //
}

//_____________________________________________________________________________
Int_t AliVZERO::DistanceToPrimitive(Int_t /*px*/, Int_t /*py*/)
{
  //
  // Calculates the distance from the mouse to the VZERO on the screen
  // Dummy routine
  //
  
  return 9999;
}
 
//_____________________________________________________________________________
void AliVZERO::Init()
{
  //
  // 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);
  
   //  Verbose level
   //  0: Silent
   //  1: cout messages
   //  2: txt files with digits 
   //  BE CAREFUL, verbose level 2 MUST be used only for debugging and
   //  it is highly suggested to use this mode only for debugging digits files
   //  reasonably small, because otherwise the size of the txt files can reach
   //  quickly several MB wasting time and disk space.
  
   ofstream ftxt;
   buffer->SetVerbose(0);
   Int_t verbose = buffer->GetVerbose();

   // 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(); 
  
   // Now retrieve the channel information: charge+time and 
   // dump it into ADC and Time arrays
   // We assume here an ordered (by PMNumber) array of
   // digits!!

   Int_t nEntries = Int_t(digits->GetEntries());
   UInt_t ADC[64];
   Int_t PMNumber[64];
   UInt_t Time[64];
  
   for (Int_t i = 0; i < nEntries; i++) {
     fVZERO->ResetDigits();
     digits->GetEvent(i);
     Int_t ndig = VZEROdigits->GetEntriesFast(); 
   
     if(ndig == 0) continue;
     if(verbose == 2) {ftxt.open("VZEROdigits.txt",ios::app);}
     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 = buffer->GetOnlineChannel(k);		
	 ADC[iChannel]       = (UInt_t) fVZERODigit->ADC();
	 PMNumber[iChannel]  = fVZERODigit->PMNumber();
	 Time[iChannel]      = (UInt_t) fVZERODigit->Time();
         if(verbose == 1) { cout <<"DDL: "<<fileName<< "\tdigit number: "<< k<<"\tPM number: "
	                    <<PMNumber[k]<<"\tADC: "<< ADC[k] << "\tTime: "<< Time[k] << endl;} 
	 if(verbose == 2) {
	      ftxt<<"DDL: "<<fileName<< "\tdigit number: "<< k<<"\tPM number: "
	                   <<PMNumber[k]<<"\tADC: "<< ADC[k] << "\tTime: "<< Time[k] << endl;}	      
//	 printf("DDL: %s, channel: %d, PM: %d, ADC: %d, Time: %d \n", 
//	            fileName,k,PMNumber[k],ADC[k],Time[k]); 
     }        
   if(verbose==2) ftxt.close();
   }

   // 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, ADC[iChannel], Time[iChannel]);       
         }
         buffer->WriteBeamFlags(); 
         buffer->WriteMBInfo(); 
         buffer->WriteMBFlags();   
         buffer->WriteBeamScalers(); 
      } 
      for(Int_t iChannel=0; iChannel < 8; iChannel++) {
          buffer->WriteTiming(iChannel, ADC[iChannel], Time[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
  
  for(Int_t i=0; i<64; i++) {
      new(pdigit) AliVZEROdigit(i, (Int_t)rawStream->GetADC(i), (Int_t)rawStream->GetTime(i)); 
      treeD->Fill();
  }
 
// Checks if everything is OK by printing results 

//   for(int i=0;i<64;i++) {
// 	printf("Channel %d : %d %d \n",i,rawStream->GetADC(i),rawStream->GetTime(i)); }
//   treeD->Print(); printf(" \n"); 
   	
  fLoader->WriteDigits("OVERWRITE");
  fLoader->UnloadDigits();	
	
  delete rawStream;

  timer.Stop();
  timer.Print();
  return kTRUE;
}
Commit	Line	Data
47890cd3	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
88cb7938	16	/* $Id$ */
47890cd3	17
	18	///////////////////////////////////////////////////////////////////////////
	19	// //
	20	// V-Zero Detector //
	21	// This class contains the base procedures for the VZERO detector //
b0d2c2d3	22	// Default geometry of November 2003 : V0R box is 4.4 cm thick //
b0d2c2d3	23	// scintillators are 2 cm thick //
47890cd3	24	// All comments should be sent to Brigitte CHEYNIS : //
61874e87	25	// b.cheynis@ipnl.in2p3.fr //
61874e87	26	// //
47890cd3	27	// //
	28	///////////////////////////////////////////////////////////////////////////
	29
88cb7938	30
b2501ea3	31	// --- Standard libraries ---
3af027ba	32	#include <Riostream.h>
b0d2c2d3	33	#include <stdlib.h>
47890cd3	34
b2501ea3	35	// --- ROOT libraries ---
b2501ea3	36	#include <TNamed.h>
b0d2c2d3	37	#include "TROOT.h"
	38	#include "TFile.h"
	39	#include "TNetFile.h"
	40	#include "TRandom.h"
	41	#include "TTree.h"
	42	#include "TBranch.h"
	43	#include "TClonesArray.h"
	44	#include "TStopwatch.h"
b2501ea3	45
	46	// --- AliRoot header files ---
	47	#include "AliRun.h"
	48	#include "AliMC.h"
47890cd3	49	#include "AliVZERO.h"
b2501ea3	50	#include "AliVZEROLoader.h"
b0d2c2d3	51	#include "AliVZERODigitizer.h"
	52	#include "AliVZEROBuffer.h"
	53	#include "AliRunDigitizer.h"
	54	#include "AliVZEROdigit.h"
362c9d61	55	#include "AliDAQ.h"
20a489c2	56	#include "AliRawReader.h"
20a489c2	57	#include "AliVZERORawStream.h"
47890cd3	58
47890cd3	59	ClassImp(AliVZERO)
0b2bea8b	60	//__________________________________________________________________
	61	AliVZERO::AliVZERO(): AliDetector(),
	62	fIdSens1(0),
	63	fThickness(0.),
	64	fThickness1(0.),
	65	fMaxStepQua(0.),
	66	fMaxStepAlu(0.),
	67	fMaxDestepQua(0.),
	68	fMaxDestepAlu(0.)
	69	{
	70	/// Default Constructor
	71
	72	AliDebug(1,Form("default (empty) ctor this=%p",this));
	73	fIshunt = 0;
	74	}
47890cd3	75	//_____________________________________________________________________________
47890cd3	76	AliVZERO::AliVZERO(const char name, const char title)
0b2bea8b	77	: AliDetector(name,title),
	78	fIdSens1(0),
	79	fThickness(4.4),
	80	fThickness1(2.0),
	81	fMaxStepQua(0.05),
	82	fMaxStepAlu(0.01),
	83	fMaxDestepQua(-1.0),
	84	fMaxDestepAlu(-1.0)
47890cd3	85	{
0b2bea8b	86
47890cd3	87	// Standard constructor for VZERO Detector
0b2bea8b	88
0b2bea8b	89	AliDebug(1,Form("ctor this=%p",this));
47890cd3	90
b0d2c2d3	91	// fIshunt = 1; // All hits are associated with primary particles
47890cd3	92
61874e87	93	fHits = new TClonesArray("AliVZEROhit", 400);
	94	fDigits = new TClonesArray("AliVZEROdigit",400);
	95
5d12ce38	96	gAlice->GetMCApp()->AddHitList(fHits);
61874e87	97
0b2bea8b	98	// fThickness = 4.4; // total thickness of the V0R box in cm
	99	// fThickness1 = 2.0; // thickness of scintillating cells in cm
	100	//
	101	// fMaxStepQua = 0.05;
	102	// fMaxStepAlu = 0.01;
	103	//
	104	// fMaxDestepQua = -1.0;
	105	// fMaxDestepAlu = -1.0;
47890cd3	106
47890cd3	107
61874e87	108	}
	109
	110	//_____________________________________________________________________________
	111	AliVZERO::~AliVZERO()
	112	{
b2501ea3	113	//
	114	// Default destructor for VZERO Detector
	115	//
	116
61874e87	117	if (fHits) {
	118	fHits->Delete();
	119	delete fHits;
b0d2c2d3	120	fHits=0; }
b2501ea3	121
	122	if (fDigits) {
	123	fDigits->Delete();
	124	delete fDigits;
b0d2c2d3	125	fDigits=0; }
61874e87	126	}
47890cd3	127
	128	//_____________________________________________________________________________
	129	void AliVZERO::BuildGeometry()
	130	{
	131	//
b2501ea3	132	// Builds simple ROOT TNode geometry for event display
47890cd3	133	//
	134	}
	135
	136	//_____________________________________________________________________________
	137	void AliVZERO::CreateGeometry()
	138	{
	139	//
b2501ea3	140	// Builds simple Geant3 geometry
47890cd3	141	//
	142	}
	143	//_____________________________________________________________________________
	144	void AliVZERO::CreateMaterials()
	145	{
	146	//
b2501ea3	147	// Creates materials used for Geant3 geometry
47890cd3	148	//
47890cd3	149	}
61874e87	150
47890cd3	151	//_____________________________________________________________________________
652e94c3	152	Int_t AliVZERO::DistanceToPrimitive(Int_t /px/, Int_t /py/)
47890cd3	153	{
47890cd3	154	//
b2501ea3	155	// Calculates the distance from the mouse to the VZERO on the screen
47890cd3	156	// Dummy routine
	157	//
	158
	159	return 9999;
	160	}
	161
b2501ea3	162	//_____________________________________________________________________________
47890cd3	163	void AliVZERO::Init()
	164	{
	165	//
b2501ea3	166	// Initialises the VZERO class after it has been built
47890cd3	167	//
	168	}
	169
61874e87	170
b2501ea3	171	//_____________________________________________________________________________
47890cd3	172	void AliVZERO::SetMaxStepQua(Float_t p1)
47890cd3	173	{
b2501ea3	174	//
	175	// Possible parametrisation of steps in active materials
	176	//
47890cd3	177	fMaxStepQua = p1;
	178	}
	179
b2501ea3	180	//_____________________________________________________________________________
47890cd3	181	void AliVZERO::SetMaxStepAlu(Float_t p1)
47890cd3	182	{
b2501ea3	183	//
	184	// Possible parametrisation of steps in Aluminum foils (not used in
	185	// version v2)
	186	//
47890cd3	187	fMaxStepAlu = p1;
	188	}
	189
b2501ea3	190	//_____________________________________________________________________________
47890cd3	191	void AliVZERO::SetMaxDestepQua(Float_t p1)
47890cd3	192	{
b2501ea3	193	//
	194	// Possible parametrisation of steps in active materials (quartz)
	195	//
47890cd3	196	fMaxDestepQua = p1;
	197	}
	198
b2501ea3	199	//_____________________________________________________________________________
47890cd3	200	void AliVZERO::SetMaxDestepAlu(Float_t p1)
47890cd3	201	{
b2501ea3	202	//
	203	// Possible parametrisation of steps in Aluminum (not used in
	204	// version v2)
	205	//
47890cd3	206	fMaxDestepAlu = p1;
47890cd3	207	}
4dbb3fd3	208
b2501ea3	209	//_____________________________________________________________________________
4dbb3fd3	210	AliLoader* AliVZERO::MakeLoader(const char* topfoldername)
4dbb3fd3	211	{
b2501ea3	212	//
b2501ea3	213	// Builds VZEROgetter (AliLoader type)
4dbb3fd3	214	// if detector wants to use customized getter, it must overload this method
b2501ea3	215	//
3e87825e	216	// Info("MakeLoader","Creating AliVZEROLoader. Top folder is %s.",topfoldername);
	217
	218	AliDebug(1,Form("Creating AliVZEROLoader, Top folder is %s ",topfoldername));
4dbb3fd3	219	fLoader = new AliVZEROLoader(GetName(),topfoldername);
	220	return fLoader;
	221	}
	222
b2501ea3	223	//_____________________________________________________________________________
	224	void AliVZERO::SetTreeAddress()
	225	{
b0d2c2d3	226	//
0e305a84	227	// Sets tree address for hits.
b2501ea3	228	//
0e305a84	229	if (fLoader->TreeH() && (fHits == 0x0))
	230	fHits = new TClonesArray("AliVZEROhit", 400);
	231
	232	AliDetector::SetTreeAddress();
	233	}
4dbb3fd3	234
b0d2c2d3	235	//_____________________________________________________________________________
	236	AliDigitizer* AliVZERO::CreateDigitizer(AliRunDigitizer* manager) const
	237	{
	238	//
	239	// Creates a digitizer for VZERO
	240	//
	241	return new AliVZERODigitizer(manager);
	242	}
	243
	244	//_____________________________________________________________________________
	245	void AliVZERO::Hits2Digits(){
	246	//
	247	// Converts hits to digits of the current event
	248	//
	249	// Inputs file name
d3d4a92f	250	const char *alifile = "galice.root";
b0d2c2d3	251
	252	// Create the run digitizer
	253	AliRunDigitizer* manager = new AliRunDigitizer(1, 1);
	254	manager->SetInputStream(0, alifile);
	255	manager->SetOutputFile("H2Dfile");
	256
	257	// Creates the VZERO digitizer
	258	AliVZERODigitizer* dig = new AliVZERODigitizer(manager);
	259
	260	// Creates the digits
	261	dig->Exec("");
	262
	263	}
	264	//_____________________________________________________________________________
	265	void AliVZERO::Digits2Raw()
	266	{
c4f89fd4	267	//
	268	// Converts digits of the current event to raw data
	269	//
b0d2c2d3	270
c4f89fd4	271	AliVZERO fVZERO = (AliVZERO)gAlice->GetDetector("VZERO");
	272	fLoader->LoadDigits();
	273	TTree* digits = fLoader->TreeD();
	274	if (!digits) {
	275	Error("Digits2Raw", "no digits tree");
	276	return;
	277	}
	278	TClonesArray * VZEROdigits = new TClonesArray("AliVZEROdigit",1000);
	279	fVZERO->SetTreeAddress();
	280	digits->GetBranch("VZERODigit")->SetAddress(&VZEROdigits);
b0d2c2d3	281
c4f89fd4	282	const char *fileName = AliDAQ::DdlFileName("VZERO",0);
c4f89fd4	283	AliVZEROBuffer* buffer = new AliVZEROBuffer(fileName);
b0d2c2d3	284
c4f89fd4	285	// Verbose level
	286	// 0: Silent
	287	// 1: cout messages
	288	// 2: txt files with digits
	289	// BE CAREFUL, verbose level 2 MUST be used only for debugging and
	290	// it is highly suggested to use this mode only for debugging digits files
	291	// reasonably small, because otherwise the size of the txt files can reach
	292	// quickly several MB wasting time and disk space.
	293
	294	ofstream ftxt;
	295	buffer->SetVerbose(0);
	296	Int_t verbose = buffer->GetVerbose();
	297
	298	// Get Trigger information first
	299	// Read trigger inputs from trigger-detector object
	300	AliDataLoader * dataLoader = fLoader->GetDigitsDataLoader();
	301	if( !dataLoader->IsFileOpen() )
	302	dataLoader->OpenFile( "READ" );
	303	AliTriggerDetector* trgdet = (AliTriggerDetector*)dataLoader->GetDirectory()->Get( "Trigger" );
	304	UInt_t triggerInfo = 0;
	305	if(trgdet) {
	306	triggerInfo = trgdet->GetMask() & 0xffff;
	307	}
	308	else {
	309	AliError(Form("There is no trigger object for %s",fLoader->GetName()));
	310	}
b0d2c2d3	311
c4f89fd4	312	buffer->WriteTriggerInfo((UInt_t)triggerInfo);
	313	buffer->WriteTriggerScalers();
	314	buffer->WriteBunchNumbers();
	315
	316	// Now retrieve the channel information: charge+time and
	317	// dump it into ADC and Time arrays
	318	// We assume here an ordered (by PMNumber) array of
	319	// digits!!
726d762c	320
c4f89fd4	321	Int_t nEntries = Int_t(digits->GetEntries());
	322	UInt_t ADC[64];
	323	Int_t PMNumber[64];
	324	UInt_t Time[64];
b0d2c2d3	325
c4f89fd4	326	for (Int_t i = 0; i < nEntries; i++) {
	327	fVZERO->ResetDigits();
	328	digits->GetEvent(i);
	329	Int_t ndig = VZEROdigits->GetEntriesFast();
b0d2c2d3	330
c4f89fd4	331	if(ndig == 0) continue;
	332	if(verbose == 2) {ftxt.open("VZEROdigits.txt",ios::app);}
	333	for(Int_t k=0; k<ndig; k++){
7c786908	334	AliVZEROdigit* fVZERODigit = (AliVZEROdigit*) VZEROdigits->At(k);
	335	// Convert aliroot channel k into FEE channel iChannel before writing data
	336	Int_t iChannel = buffer->GetOnlineChannel(k);
	337	ADC[iChannel] = (UInt_t) fVZERODigit->ADC();
	338	PMNumber[iChannel] = fVZERODigit->PMNumber();
	339	Time[iChannel] = (UInt_t) fVZERODigit->Time();
c4f89fd4	340	if(verbose == 1) { cout <<"DDL: "<<fileName<< "\tdigit number: "<< k<<"\tPM number: "
	341	<<PMNumber[k]<<"\tADC: "<< ADC[k] << "\tTime: "<< Time[k] << endl;}
	342	if(verbose == 2) {
	343	ftxt<<"DDL: "<<fileName<< "\tdigit number: "<< k<<"\tPM number: "
	344	<<PMNumber[k]<<"\tADC: "<< ADC[k] << "\tTime: "<< Time[k] << endl;}
	345	// printf("DDL: %s, channel: %d, PM: %d, ADC: %d, Time: %d \n",
	346	// fileName,k,PMNumber[k],ADC[k],Time[k]);
	347	}
	348	if(verbose==2) ftxt.close();
	349	}
b0d2c2d3	350
c4f89fd4	351	// Now fill raw data
	352
	353	for (Int_t iCIU = 0; iCIU < 8; iCIU++) {
	354
	355	// decoding of one Channel Interface Unit numbered iCIU - there are 8 channels per CIU (and 8 CIUs) :
	356
	357	for(Int_t iChannel_Offset = iCIU8; iChannel_Offset < (iCIU8)+8; iChannel_Offset=iChannel_Offset+4) {
	358	for(Int_t iChannel = iChannel_Offset; iChannel < iChannel_Offset+4; iChannel++) {
	359	buffer->WriteChannel(iChannel, ADC[iChannel], Time[iChannel]);
	360	}
	361	buffer->WriteBeamFlags();
	362	buffer->WriteMBInfo();
	363	buffer->WriteMBFlags();
	364	buffer->WriteBeamScalers();
	365	}
c61a7285	366	for(Int_t iChannel=0; iChannel < 8; iChannel++) {
c61a7285	367	buffer->WriteTiming(iChannel, ADC[iChannel], Time[iChannel]);
c4f89fd4	368	}
726d762c	369
c4f89fd4	370	// End of decoding of one CIU card
	371
	372	} // end of decoding the eight CIUs
	373
b0d2c2d3	374	delete buffer;
c4f89fd4	375	fLoader->UnloadDigits();
b0d2c2d3	376	}
b0d2c2d3	377
20a489c2	378	//_____________________________________________________________________________
	379	Bool_t AliVZERO::Raw2SDigits(AliRawReader* rawReader){
	380	// Converts the VZERO raw data into digits
	381	// The method is used for merging simulated and
	382	// real data events
	383	TStopwatch timer;
	384	timer.Start();
	385
	386	if(!fLoader) {
	387	AliError("no VZERO loader found");
	388	return kFALSE; }
	389
	390	TTree* treeD = fLoader->TreeD();
	391	if(!treeD) {
	392	fLoader->MakeTree("D");
	393	treeD = fLoader->TreeD(); }
	394
	395	AliVZEROdigit digit;
	396	AliVZEROdigit* pdigit = &digit;
	397	const Int_t kBufferSize = 4000;
	398
	399	treeD->Branch("VZERO", "AliVZEROdigit", &pdigit, kBufferSize);
	400
	401	rawReader->Reset();
	402	AliVZERORawStream* rawStream = new AliVZERORawStream(rawReader);
	403
	404	if (!rawStream->Next()) return kFALSE; // No VZERO data found
	405
	406	for(Int_t i=0; i<64; i++) {
	407	new(pdigit) AliVZEROdigit(i, (Int_t)rawStream->GetADC(i), (Int_t)rawStream->GetTime(i));
	408	treeD->Fill();
	409	}
	410
	411	// Checks if everything is OK by printing results
	412
	413	// for(int i=0;i<64;i++) {
	414	// printf("Channel %d : %d %d \n",i,rawStream->GetADC(i),rawStream->GetTime(i)); }
	415	// treeD->Print(); printf(" \n");
	416
	417	fLoader->WriteDigits("OVERWRITE");
	418	fLoader->UnloadDigits();
	419
	420	delete rawStream;
	421
	422	timer.Stop();
	423	timer.Print();
c61a7285	424	return kTRUE;
20a489c2	425	}
20a489c2	426
4dbb3fd3	427