[u/mrichter/AliRoot.git] / START / AliSTARTDigitizer.cxx


/**************************************************************************
 * Copyright(c) 1998-2000, 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.                  *
 **************************************************************************/


#include <TTree.h> 
#include <TVector.h>
#include <TObjArray.h>
#include <TFile.h>
#include <TDirectory.h>
#include <TRandom.h>
#include <TArrayI.h>
#include <TH1.h>
#include <TError.h>


#include "AliSTARTDigitizer.h"
#include "AliSTART.h"
#include "AliSTARThit.h"
#include "AliSTARThitPhoton.h"
#include "AliSTARTdigit.h"
#include "AliRunDigitizer.h"
#include "AliHeader.h"
#include "AliGenEventHeader.h"
#include "AliRun.h"
#include "AliPDG.h"
#include "AliLoader.h"
#include "AliRunLoader.h"
#include "AliSTARTLoader.h"

#include <stdlib.h>
#include <Riostream.h>
#include <Riostream.h>

ClassImp(AliSTARTDigitizer)

//___________________________________________
  AliSTARTDigitizer::AliSTARTDigitizer()  :AliDigitizer()
{
// Default ctor - don't use it
  ;
}

//___________________________________________
AliSTARTDigitizer::AliSTARTDigitizer(AliRunDigitizer* manager) 
    :AliDigitizer(manager) 
{
	cout<<"AliSTARTDigitizer::AliSTARTDigitizer"<<endl;
// ctor which should be used
//  fDebug =0;
 // if (GetDebug()>2)
  //  cerr<<"AliSTARTDigitizer::AliSTARTDigitizer"
   //     <<"(AliRunDigitizer* manager) was processed"<<endl;

}

//------------------------------------------------------------------------
AliSTARTDigitizer::~AliSTARTDigitizer()
{
// Destructor


}

 //------------------------------------------------------------------------
Bool_t AliSTARTDigitizer::Init()
{
// Initialization
 cout<<"AliSTARTDigitizer::Init"<<endl;
 return kTRUE;
}
 

//---------------------------------------------------------------------

void AliSTARTDigitizer::Exec(Option_t* /*option*/)
{


  AliRunLoader *inRL, *outRL;//in and out Run Loaders
  AliLoader *ingime, *outgime;// in and out ITSLoaders

  outRL = AliRunLoader::GetRunLoader(fManager->GetOutputFolderName());
  outgime = outRL->GetLoader("STARTLoader");

#ifdef DEBUG
  cout<<"AliSTARTDigitizer::>SDigits2Digits start...\n";
#endif
  //
  // From hits to digits
  //
  Int_t hit, nhits;
  Int_t CountEr[13],CountEl[13];							//!!!
  Int_t volume,pmt,tr,tl,sumRight;
  Float_t timediff,timeav;
  Float_t besttimeright,besttimeleft,meanTime;
  Int_t  bestRightADC,bestLeftADC;
  Float_t besttimeleftGaus, besttimerightGaus;
  Float_t timeright[13]={13*0};
  Float_t timeleft[13]={13*0};
  Float_t channelWidth=2.5; //ps
  Int_t channelWidthADC=1; //ps
  //  Int_t thresholdAmpl=10;

  ftimeRightTDC = new TArrayI(12); 
  ftimeLeftTDC = new TArrayI(12); 
  fRightADC = new TArrayI(12); 
  fLeftADC = new TArrayI(12); 
  
  inRL = AliRunLoader::GetRunLoader(fManager->GetInputFolderName(0));
  inRL->LoadgAlice();
  
  //  fHits = new TClonesArray ("AliSTARThit", 1000);
  fPhotons = new TClonesArray ("AliSTARThitPhoton", 10000);			//!!!
  AliSTART *START  = (AliSTART*) gAlice->GetDetector("START");
  AliSTARThit  *startHit;
  //use if Cherenkov photons
  //  AliSTARThitPhoton  *startHitPhoton;						//!!!
  TBranch *brHits=0;
  TBranch *brHitPhoton=0;
  fdigits= new AliSTARTdigit();

  Int_t nFiles=fManager->GetNinputs();
  for (Int_t inputFile=0; inputFile<nFiles;  inputFile++) {

    besttimeright=9999.;
    besttimeleft=9999.;
    Int_t timeDiff=0;
    Int_t timeAv=0;
    sumRight=0;
    for (Int_t i0=0; i0<13; i0++)
      {
	timeright[i0]=0; timeleft[i0]=0;
	CountEr[i0]=0;   CountEl[i0]=0;
      }

    inRL = AliRunLoader::GetRunLoader(fManager->GetInputFolderName(inputFile));
    ingime = inRL->GetLoader("STARTLoader");
    ingime->LoadHits("READ");//probably it is necessary to load them before
    outgime->LoadDigits("UPDATE");//probably it is necessary to load them before
    //use if Cherenkov photons
    //  TClonesArray *STARThitsPhotons = START->Photons ();
    TClonesArray *fHits = START->Hits ();
    //    cout<<" Load  "<<AliSTARTLoader::LoadDigits()<<endl;

    TTree *th = ingime->TreeH();
    brHits = th->GetBranch("START");
    brHitPhoton = th->GetBranch("STARThitPhoton");
    if (brHits) {
      START->SetHitsAddressBranch(brHits,brHitPhoton);
    }else{
      cerr<<"EXEC Branch START hit not found"<<endl;
      exit(111);
    } 
    Int_t ntracks    = (Int_t) th->GetEntries();
    cout<<" ntracks "<<ntracks<<endl;
    if (ntracks<=0) return;
    // Start loop on tracks in the photon hits containers 
    // for amplitude
    /*
    if(brHitPhoton) {
      cout<<"brHitPhoton "<<endl; 
      for (Int_t track=0; track<ntracks;track++) {
	brHitPhoton -> GetEntry(track);;
	nhits = STARThitsPhotons->GetEntriesFast();
	for (hit=0;hit<nhits;hit++) {
	  startHitPhoton   = (AliSTARThitPhoton*) 
	   STARThitsPhotons ->UncheckedAt(hit);
	  pmt=startHitPhoton->fPmt;
	  volume = startHitPhoton->fArray;
	  if(RegisterPhotoE(startHitPhoton)) 
	    {
	      if (volume == 1) CountEr[pmt]++;
	      if (volume == 2) CountEl[pmt]++;
	    }
	} //hit photons
      } //track photons
    } // was photons
    */
    // Start loop on tracks in the hits containers
    for (Int_t track=0; track<ntracks;track++) {
      brHits->GetEntry(track);
      nhits = fHits->GetEntriesFast();
      //  cout<<" brHits hits "<<nhits<<endl;
      for (hit=0;hit<nhits;hit++) {
	startHit   = (AliSTARThit*) fHits->UncheckedAt(hit);
	if (!startHit) {
	  ::Error("Exec","The unchecked hit doesn't exist");
	  break;
	}
	pmt=startHit->Pmt();
	volume = startHit->Volume();
	if(volume==1){
	  timeright[pmt] = startHit->Time();
	  if(timeright[pmt]<besttimeright)
	    //&&CountEr[pmt-1]>thresholdAmpl)
	    {
	    besttimeright=timeright[pmt];
	  } //timeright
	}//time for right shoulder
	if(volume==2){            
	  timeleft[pmt] = startHit->Time();
	  if(timeleft[pmt]<besttimeleft)
	    //&&CountEl[pmt-1]>thresholdAmpl) 
	    {
	    besttimeleft=timeleft[pmt];
	    
	  } //timeleftbest
	}//time for left shoulder
      } //hit loop
    } //track loop
  
    // z position
    cout<<" right time  "<<besttimeright<<
      " right distance "<<besttimeright*30<<endl;;
    cout<<" left time  "<<besttimeleft<<
      " left distance "<<besttimeleft*30<<endl;;
  

    //folding with experimental time distribution
    
    besttimeleftGaus=gRandom->Gaus(besttimeright,0.05);
    cout<<" besttimeleftGaus "<<besttimeleftGaus<<endl;
    bestLeftADC=Int_t (besttimeleftGaus*1000/channelWidth);
    Float_t koef=69.7/350.;
    besttimeright=koef*besttimeleft;
    besttimerightGaus=gRandom->Gaus(besttimeleft,0.05);
    
    bestRightADC=Int_t (besttimerightGaus*1000/channelWidth);
    timediff=besttimerightGaus-besttimeleftGaus;
    cout<<" timediff in ns "<<timediff<<" z= "<<timediff*30<<endl;
    meanTime=(besttimerightGaus+besttimeleftGaus)/2.;
    if ( TMath::Abs(timediff)<TMath::Abs(0.3) ) 
      {
	Float_t t1=1000.*besttimeleftGaus;
	Float_t t2=1000.*besttimerightGaus;
	t1=t1/channelWidth;   //time in ps to channelWidth
	t2=t2/channelWidth;   //time in ps to channelWidth
	timeav=(t1+t2)/2.;
	
	// Time to TDC signal
	// 256 channels for timediff, range 1ns
	
	timediff=512+1000*timediff/channelWidth; // time in ps 
	
	timeAv = (Int_t)(timeav);   // time  channel numbres
	timeDiff = (Int_t)(timediff); // time  channel numbres
	//       fill digits
	fdigits->SetTimeBestLeft(bestLeftADC);
	fdigits->SetTimeBestRight(bestRightADC);
	fdigits->SetMeanTime(timeAv);
	fdigits->SetTimeDiff(timeDiff);
	for (Int_t i=0; i<12; i++)
	  {
	    //  fill TDC
	    timeright[i+1]=gRandom->Gaus(timeright[i+1],0.05);
	    timeleft[i+1]=gRandom->Gaus(timeleft[i+1],0.05);
	    tr= Int_t (timeright[i+1]*1000/channelWidth); 
	    if(tr<200) tr=0;
	    tl= Int_t (timeleft[i+1]*1000/channelWidth); 
	    if(tl<1000) tl=0;
	    
	    ftimeRightTDC->AddAt(tr,i);
	    ftimeLeftTDC->AddAt(tl,i);
	    //fill ADC
	    Int_t al=( Int_t ) CountEl[i+1]/ channelWidthADC;
	    Int_t ar=( Int_t ) CountEr[i+1]/ channelWidthADC;
	    fRightADC->AddAt(ar,i);
	    fLeftADC ->AddAt(al,i);
	    sumRight+=CountEr[i+1];
	  }
	fdigits->SetTimeRight(*ftimeRightTDC);
	fdigits->SetTimeLeft(*ftimeLeftTDC);
	fdigits->SetADCRight(*fRightADC);
	fdigits->SetADCLeft(*fLeftADC);
	// cout<<" before sum"<<endl;
	fdigits->SetSumADCRight(sumRight);
      }
    else
      {timeAv=999999; timeDiff=99999;}

// trick to find out output dir:


/*
   // trick to find out output dir:
    TTree *outTree = fManager->GetTreeD();
    if (!outTree) {
      cerr<<"something wrong with output...."<<endl;
      exit(111);
    }

    Char_t nameDigits[20];
    TDirectory *wd = gDirectory;
    outTree->GetDirectory()->cd();
    fdigits->Write(nameDigits);
    cout<<nameDigits<<endl;
    wd->cd();
*/  

     Char_t nameDigits[20];
    sprintf(nameDigits,"START_D_%d",fManager->GetOutputEventNr());
    fdigits->Write(nameDigits);

    //    outgime->WriteDigits("OVERWRITE");
  }
}


//------------------------------------------------------------------------
Bool_t AliSTARTDigitizer::RegisterPhotoE(/*AliSTARThitPhoton *hit*/)
{
    Double_t    P = 0.2;    
    Double_t    p;
    
    p = gRandom->Rndm();
    if (p > P)
      return kFALSE;
    
    return kTRUE;
}
//----------------------------------------------------------------------------
Commit	Line	Data
ede9aff7	1
	2	/**************************************************************************
	3	* Copyright(c) 1998-2000, ALICE Experiment at CERN, All rights reserved. *
	4	* *
	5	* Author: The ALICE Off-line Project. *
	6	* Contributors are mentioned in the code where appropriate. *
	7	* *
	8	* Permission to use, copy, modify and distribute this software and its *
	9	* documentation strictly for non-commercial purposes is hereby granted *
	10	* without fee, provided that the above copyright notice appears in all *
	11	* copies and that both the copyright notice and this permission notice *
	12	* appear in the supporting documentation. The authors make no claims *
	13	* about the suitability of this software for any purpose. It is *
	14	* provided "as is" without express or implied warranty. *
	15	**************************************************************************/
	16
	17
	18	#include <TTree.h>
	19	#include <TVector.h>
	20	#include <TObjArray.h>
	21	#include <TFile.h>
	22	#include <TDirectory.h>
	23	#include <TRandom.h>
e73d68f2	24	#include <TArrayI.h>
e73d68f2	25	#include <TH1.h>
2b4678e2	26	#include <TError.h>
ede9aff7	27
	28
	29	#include "AliSTARTDigitizer.h"
	30	#include "AliSTART.h"
	31	#include "AliSTARThit.h"
e73d68f2	32	#include "AliSTARThitPhoton.h"
ede9aff7	33	#include "AliSTARTdigit.h"
ede9aff7	34	#include "AliRunDigitizer.h"
e73d68f2	35	#include "AliHeader.h"
e73d68f2	36	#include "AliGenEventHeader.h"
ede9aff7	37	#include "AliRun.h"
ede9aff7	38	#include "AliPDG.h"
88cb7938	39	#include "AliLoader.h"
88cb7938	40	#include "AliRunLoader.h"
7a68e0ef	41	#include "AliSTARTLoader.h"
ede9aff7	42
ede9aff7	43	#include <stdlib.h>
ef0750c2	44	#include <Riostream.h>
ef0750c2	45	#include <Riostream.h>
ede9aff7	46
	47	ClassImp(AliSTARTDigitizer)
	48
	49	//___________________________________________
	50	AliSTARTDigitizer::AliSTARTDigitizer() :AliDigitizer()
	51	{
	52	// Default ctor - don't use it
	53	;
	54	}
	55
	56	//___________________________________________
	57	AliSTARTDigitizer::AliSTARTDigitizer(AliRunDigitizer* manager)
	58	:AliDigitizer(manager)
	59	{
	60	cout<<"AliSTARTDigitizer::AliSTARTDigitizer"<<endl;
	61	// ctor which should be used
	62	// fDebug =0;
	63	// if (GetDebug()>2)
	64	// cerr<<"AliSTARTDigitizer::AliSTARTDigitizer"
	65	// <<"(AliRunDigitizer* manager) was processed"<<endl;
	66
	67	}
	68
	69	//------------------------------------------------------------------------
	70	AliSTARTDigitizer::~AliSTARTDigitizer()
	71	{
	72	// Destructor
	73
	74
	75	}
	76
	77	//------------------------------------------------------------------------
	78	Bool_t AliSTARTDigitizer::Init()
	79	{
	80	// Initialization
	81	cout<<"AliSTARTDigitizer::Init"<<endl;
	82	return kTRUE;
	83	}
	84
	85
	86	//---------------------------------------------------------------------
	87
dc63cec5	88	void AliSTARTDigitizer::Exec(Option_t* /option/)
ede9aff7	89	{
ede9aff7	90
88cb7938	91
	92	AliRunLoader inRL, outRL;//in and out Run Loaders
	93	AliLoader ingime, outgime;// in and out ITSLoaders
	94
	95	outRL = AliRunLoader::GetRunLoader(fManager->GetOutputFolderName());
	96	outgime = outRL->GetLoader("STARTLoader");
	97
ede9aff7	98	#ifdef DEBUG
	99	cout<<"AliSTARTDigitizer::>SDigits2Digits start...\n";
	100	#endif
	101	//
	102	// From hits to digits
	103	//
e73d68f2	104	Int_t hit, nhits;
	105	Int_t CountEr[13],CountEl[13]; //!!!
	106	Int_t volume,pmt,tr,tl,sumRight;
e73d68f2	107	Float_t timediff,timeav;
ede9aff7	108	Float_t besttimeright,besttimeleft,meanTime;
e73d68f2	109	Int_t bestRightADC,bestLeftADC;
	110	Float_t besttimeleftGaus, besttimerightGaus;
	111	Float_t timeright[13]={13*0};
	112	Float_t timeleft[13]={13*0};
	113	Float_t channelWidth=2.5; //ps
	114	Int_t channelWidthADC=1; //ps
1ab4b835	115	// Int_t thresholdAmpl=10;
e73d68f2	116
	117	ftimeRightTDC = new TArrayI(12);
	118	ftimeLeftTDC = new TArrayI(12);
	119	fRightADC = new TArrayI(12);
	120	fLeftADC = new TArrayI(12);
	121
7a68e0ef	122	inRL = AliRunLoader::GetRunLoader(fManager->GetInputFolderName(0));
	123	inRL->LoadgAlice();
	124
	125	// fHits = new TClonesArray ("AliSTARThit", 1000);
e73d68f2	126	fPhotons = new TClonesArray ("AliSTARThitPhoton", 10000); //!!!
ede9aff7	127	AliSTART START = (AliSTART) gAlice->GetDetector("START");
ede9aff7	128	AliSTARThit *startHit;
dc63cec5	129	//use if Cherenkov photons
dc63cec5	130	// AliSTARThitPhoton *startHitPhoton; //!!!
ede9aff7	131	TBranch *brHits=0;
e73d68f2	132	TBranch *brHitPhoton=0;
ede9aff7	133	fdigits= new AliSTARTdigit();
	134
	135	Int_t nFiles=fManager->GetNinputs();
	136	for (Int_t inputFile=0; inputFile<nFiles; inputFile++) {
e73d68f2	137
ede9aff7	138	besttimeright=9999.;
	139	besttimeleft=9999.;
	140	Int_t timeDiff=0;
	141	Int_t timeAv=0;
e73d68f2	142	sumRight=0;
	143	for (Int_t i0=0; i0<13; i0++)
	144	{
	145	timeright[i0]=0; timeleft[i0]=0;
	146	CountEr[i0]=0; CountEl[i0]=0;
	147	}
88cb7938	148
	149	inRL = AliRunLoader::GetRunLoader(fManager->GetInputFolderName(inputFile));
	150	ingime = inRL->GetLoader("STARTLoader");
	151	ingime->LoadHits("READ");//probably it is necessary to load them before
f803559b	152	outgime->LoadDigits("UPDATE");//probably it is necessary to load them before
dc63cec5	153	//use if Cherenkov photons
dc63cec5	154	// TClonesArray *STARThitsPhotons = START->Photons ();
7a68e0ef	155	TClonesArray *fHits = START->Hits ();
7a68e0ef	156	// cout<<" Load "<<AliSTARTLoader::LoadDigits()<<endl;
ede9aff7	157
7a68e0ef	158	TTree *th = ingime->TreeH();
ede9aff7	159	brHits = th->GetBranch("START");
e73d68f2	160	brHitPhoton = th->GetBranch("STARThitPhoton");
ede9aff7	161	if (brHits) {
e73d68f2	162	START->SetHitsAddressBranch(brHits,brHitPhoton);
ede9aff7	163	}else{
00907af4	164	cerr<<"EXEC Branch START hit not found"<<endl;
00907af4	165	exit(111);
ede9aff7	166	}
ede9aff7	167	Int_t ntracks = (Int_t) th->GetEntries();
e73d68f2	168	cout<<" ntracks "<<ntracks<<endl;
ede9aff7	169	if (ntracks<=0) return;
e73d68f2	170	// Start loop on tracks in the photon hits containers
e73d68f2	171	// for amplitude
7a68e0ef	172	/*
e73d68f2	173	if(brHitPhoton) {
7a68e0ef	174	cout<<"brHitPhoton "<<endl;
e73d68f2	175	for (Int_t track=0; track<ntracks;track++) {
	176	brHitPhoton -> GetEntry(track);;
	177	nhits = STARThitsPhotons->GetEntriesFast();
	178	for (hit=0;hit<nhits;hit++) {
	179	startHitPhoton = (AliSTARThitPhoton*)
	180	STARThitsPhotons ->UncheckedAt(hit);
	181	pmt=startHitPhoton->fPmt;
	182	volume = startHitPhoton->fArray;
	183	if(RegisterPhotoE(startHitPhoton))
	184	{
	185	if (volume == 1) CountEr[pmt]++;
	186	if (volume == 2) CountEl[pmt]++;
	187	}
	188	} //hit photons
	189	} //track photons
	190	} // was photons
7a68e0ef	191	*/
ede9aff7	192	// Start loop on tracks in the hits containers
	193	for (Int_t track=0; track<ntracks;track++) {
	194	brHits->GetEntry(track);
7a68e0ef	195	nhits = fHits->GetEntriesFast();
7a68e0ef	196	// cout<<" brHits hits "<<nhits<<endl;
ede9aff7	197	for (hit=0;hit<nhits;hit++) {
7a68e0ef	198	startHit = (AliSTARThit*) fHits->UncheckedAt(hit);
2b4678e2	199	if (!startHit) {
	200	::Error("Exec","The unchecked hit doesn't exist");
	201	break;
	202	}
0b73602e	203	pmt=startHit->Pmt();
0b73602e	204	volume = startHit->Volume();
ede9aff7	205	if(volume==1){
0b73602e	206	timeright[pmt] = startHit->Time();
e73d68f2	207	if(timeright[pmt]<besttimeright)
	208	//&&CountEr[pmt-1]>thresholdAmpl)
	209	{
	210	besttimeright=timeright[pmt];
ede9aff7	211	} //timeright
	212	}//time for right shoulder
	213	if(volume==2){
0b73602e	214	timeleft[pmt] = startHit->Time();
e73d68f2	215	if(timeleft[pmt]<besttimeleft)
	216	//&&CountEl[pmt-1]>thresholdAmpl)
	217	{
	218	besttimeleft=timeleft[pmt];
	219
ede9aff7	220	} //timeleftbest
	221	}//time for left shoulder
	222	} //hit loop
	223	} //track loop
dc63cec5	224
e73d68f2	225	// z position
	226	cout<<" right time "<<besttimeright<<
	227	" right distance "<<besttimeright*30<<endl;;
	228	cout<<" left time "<<besttimeleft<<
dc63cec5	229	" left distance "<<besttimeleft*30<<endl;;
e73d68f2	230
e73d68f2	231
ede9aff7	232	//folding with experimental time distribution
e73d68f2	233
dc63cec5	234	besttimeleftGaus=gRandom->Gaus(besttimeright,0.05);
	235	cout<<" besttimeleftGaus "<<besttimeleftGaus<<endl;
	236	bestLeftADC=Int_t (besttimeleftGaus*1000/channelWidth);
7115262b	237	Float_t koef=69.7/350.;
dc63cec5	238	besttimeright=koef*besttimeleft;
dc63cec5	239	besttimerightGaus=gRandom->Gaus(besttimeleft,0.05);
e73d68f2	240
dc63cec5	241	bestRightADC=Int_t (besttimerightGaus*1000/channelWidth);
ede9aff7	242	timediff=besttimerightGaus-besttimeleftGaus;
e73d68f2	243	cout<<" timediff in ns "<<timediff<<" z= "<<timediff*30<<endl;
ede9aff7	244	meanTime=(besttimerightGaus+besttimeleftGaus)/2.;
e73d68f2	245	if ( TMath::Abs(timediff)<TMath::Abs(0.3) )
	246	{
	247	Float_t t1=1000.*besttimeleftGaus;
	248	Float_t t2=1000.*besttimerightGaus;
	249	t1=t1/channelWidth; //time in ps to channelWidth
	250	t2=t2/channelWidth; //time in ps to channelWidth
	251	timeav=(t1+t2)/2.;
	252
	253	// Time to TDC signal
	254	// 256 channels for timediff, range 1ns
	255
	256	timediff=512+1000*timediff/channelWidth; // time in ps
	257
	258	timeAv = (Int_t)(timeav); // time channel numbres
	259	timeDiff = (Int_t)(timediff); // time channel numbres
	260	// fill digits
	261	fdigits->SetTimeBestLeft(bestLeftADC);
	262	fdigits->SetTimeBestRight(bestRightADC);
	263	fdigits->SetMeanTime(timeAv);
	264	fdigits->SetTimeDiff(timeDiff);
	265	for (Int_t i=0; i<12; i++)
	266	{
	267	// fill TDC
	268	timeright[i+1]=gRandom->Gaus(timeright[i+1],0.05);
	269	timeleft[i+1]=gRandom->Gaus(timeleft[i+1],0.05);
	270	tr= Int_t (timeright[i+1]*1000/channelWidth);
	271	if(tr<200) tr=0;
	272	tl= Int_t (timeleft[i+1]*1000/channelWidth);
	273	if(tl<1000) tl=0;
	274
	275	ftimeRightTDC->AddAt(tr,i);
	276	ftimeLeftTDC->AddAt(tl,i);
	277	//fill ADC
	278	Int_t al=( Int_t ) CountEl[i+1]/ channelWidthADC;
	279	Int_t ar=( Int_t ) CountEr[i+1]/ channelWidthADC;
	280	fRightADC->AddAt(ar,i);
	281	fLeftADC ->AddAt(al,i);
	282	sumRight+=CountEr[i+1];
	283	}
	284	fdigits->SetTimeRight(*ftimeRightTDC);
	285	fdigits->SetTimeLeft(*ftimeLeftTDC);
	286	fdigits->SetADCRight(*fRightADC);
	287	fdigits->SetADCLeft(*fLeftADC);
	288	// cout<<" before sum"<<endl;
	289	fdigits->SetSumADCRight(sumRight);
	290	}
ede9aff7	291	else
ede9aff7	292	{timeAv=999999; timeDiff=99999;}
88cb7938	293
88cb7938	294	// trick to find out output dir:
7a68e0ef	295
	296
	297	/*
	298	// trick to find out output dir:
	299	TTree *outTree = fManager->GetTreeD();
	300	if (!outTree) {
	301	cerr<<"something wrong with output...."<<endl;
	302	exit(111);
	303	}
	304
	305	Char_t nameDigits[20];
	306	TDirectory *wd = gDirectory;
	307	outTree->GetDirectory()->cd();
	308	fdigits->Write(nameDigits);
	309	cout<<nameDigits<<endl;
	310	wd->cd();
	311	*/
dc63cec5	312
dc63cec5	313	Char_t nameDigits[20];
7a68e0ef	314	sprintf(nameDigits,"START_D_%d",fManager->GetOutputEventNr());
	315	fdigits->Write(nameDigits);
	316
	317	// outgime->WriteDigits("OVERWRITE");
ede9aff7	318	}
ede9aff7	319	}
	320
	321
e73d68f2	322	//------------------------------------------------------------------------
dc63cec5	323	Bool_t AliSTARTDigitizer::RegisterPhotoE(/AliSTARThitPhoton hit*/)
e73d68f2	324	{
	325	Double_t P = 0.2;
	326	Double_t p;
	327
	328	p = gRandom->Rndm();
	329	if (p > P)
	330	return kFALSE;
	331
	332	return kTRUE;
	333	}
7a68e0ef	334	//----------------------------------------------------------------------------