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