ROOT logo
/**************************************************************************
 * 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$ */

///////////////////////////////////////////////////////////////////////////////
//                                                                           //
//  T0 (T-Zero) Detector                                            //
//  This class contains the base procedures for the T0     //
//  detector                                                                 //
//                                                                           //
//Begin_Html
/*
<img src="gif/AliT0Class.gif">
</pre>
<br clear=left>
<font size=+2 color=red>
<p>The responsible person for this module is
<a href="mailto:Alla.Maevskaia@cern.ch">Alla Maevskaia</a>.
</font>
<pre>
*/
//End_Html
//                                                                           //
//                                                                           //
///////////////////////////////////////////////////////////////////////////////

#include "TClonesArray.h"
#include "TString.h"

#include "AliLoader.h"
#include "AliLog.h"
#include "AliLog.h"
#include "AliMC.h"
#include "AliRun.h"
#include "AliT0.h"
#include "AliT0Digitizer.h"
#include "AliT0RawData.h"
#include "AliT0RecPoint.h"
#include "AliT0digit.h"
#include "AliT0hit.h"
ClassImp(AliT0)

  //static  AliT0digit *digits; 

//_____________________________________________________________________________
AliT0::AliT0()
  : AliDetector(), fIdSens(0), fDigits(NULL), fRecPoints(NULL)
{
  //
  // Default constructor for class AliT0
  //
  fIshunt   = 1;
  fHits     = 0;
  fDigits   = 0;
  fRecPoints = 0;
}
 
//_____________________________________________________________________________
AliT0::AliT0(const char *name, const char *title)
  : AliDetector(name,title), fIdSens(0), fDigits(new AliT0digit()), fRecPoints(new AliT0RecPoint())
{
  //
  // Standard constructor for T0 Detector
  //

  
  //
  // Initialise Hit array
  fHits       = new TClonesArray("AliT0hit",  405);
  gAlice->GetMCApp()->AddHitList(fHits);
  //  fDigits    = new AliT0digit();
  //  fRecPoints = new AliT0RecPoint();
  fIshunt     =  1;
  //  fIdSens   =  0;
  //PH  SetMarkerColor(kRed);
}

//_____________________________________________________________________________
AliT0::~AliT0() {
  
  //destructor
  if (fHits) {
    fHits->Delete();
    delete fHits;
  }
  /*
  if (fDigits) {
    fDigits->Delete();
    delete fDigits;
    cout<<" delete fDigits; "<<endl;
  }
  if (fRecPoints) {
   fRecPoints ->Delete();
    delete fRecPoints;
    cout<<" delete fRecPoints; "<<endl;
  }
  */ 
}

//_____________________________________________________________________________
void AliT0::AddHit(Int_t track, Int_t *vol, Float_t *hits)
{
  //
  // Add a T0 hit
  //
  TClonesArray &lhits = *fHits;
  new(lhits[fNhits++]) AliT0hit(fIshunt,track,vol,hits);
}


//_____________________________________________________________________________

void AliT0::AddDigit(Int_t besttimeright, Int_t besttimeleft, Int_t meantime, 
		     Int_t timediff, Int_t sumMult, Int_t refpoint,
			TArrayI *timeCFD, TArrayI *qt0, TArrayI *timeLED, TArrayI *qt1)
{
  
  //  Add a T0 digit to the list.
 //
  
  if (!fDigits) {
    fDigits = new AliT0digit();
  }
  fDigits-> SetTimeBestA(besttimeright);
  fDigits->SetTimeBestC(besttimeleft);
  fDigits-> SetMeanTime(meantime);
  fDigits-> SetDiffTime(timediff);
  fDigits-> SetSumMult(sumMult);
  fDigits->SetTimeCFD(*timeCFD);
  fDigits->SetTimeLED(*timeLED);
  fDigits->SetQT0(*qt0);
  fDigits->SetQT1(*qt1);
  fDigits->SetRefPoint(refpoint);
}

//-------------------------------------------------------------------------
void AliT0::Init()
{
  //
  // Initialis the T0 after it has been built
  Int_t i;
  //
  if(AliLog::GetGlobalDebugLevel()>0) {
    printf("\n%s: ",ClassName());
    for(i=0;i<35;i++) printf("*");
    printf(" T0_INIT ");
    for(i=0;i<35;i++) printf("*");
    printf("\n%s: ",ClassName());
    //
    // Here the T0 initialisation code (if any!)
    for(i=0;i<80;i++) printf("*");
    printf("\n");
  }
}

//---------------------------------------------------------------------------
void AliT0::MakeBranch(Option_t* option)
{
  //
// Create Tree branches for the T0.

 // Options:
  //
  //    H          Make a branch of TClonesArray of AliT0Hit's
  //    D          Make a branch of TClonesArray of AliT0Digit's
  //
  //    R         Make a branch of  AliT0RecPoints
  //
  //  char branchname[20];
  // sprintf(branchname,"%s",GetName());
  //  strncpy(branchname, GetName(), 20);
  TString branchname = Form("%s", GetName());

  const char *cH = strstr(option,"H");
  const char *cD = strstr(option,"D");
  const char *cR = strstr(option,"R");
  const char *cS = strstr(option,"S");

    if (cH && fLoader->TreeH())
  {
     if (fHits == 0x0) fHits  = new TClonesArray("AliT0hit",  405);
     AliDetector::MakeBranch(option);
  } 
    
    
  if (cD && fLoader->TreeD())
    {
      if (fDigits == 0x0) fDigits  = new AliT0digit();
      //     MakeBranchInTree(fLoader->TreeD(), branchname,
      //		       &fDigits, 405, 0);
      fLoader->TreeD()->Branch(branchname.Data(),"AliT0digit",&fDigits);
      //   fLoader->TreeD()->Print();
    } 
  if (cR && fLoader->TreeR())
    {
      if (fRecPoints == 0x0) fRecPoints  = new AliT0RecPoint();
      MakeBranchInTree(fLoader->TreeR(), branchname, &fRecPoints);
    } 
  if (cS && fLoader->TreeS())
    {
      if (fDigits == 0x0) fDigits  = new AliT0digit();
      //     MakeBranchInTree(fLoader->TreeD(), branchname,
      //		       &fDigits, 405, 0);
      fLoader->TreeS()->Branch(branchname,"AliT0digit",&fDigits);
      //   fLoader->TreeD()->Print();
    } 
  
}    

//_____________________________________________________________________________
void AliT0::ResetHits()
{
  //
  //reset hits
  //
  AliDetector::ResetHits();
  
}
//____________________________________________________________________
void AliT0::ResetDigits()
{
  //
  // Reset number of digits and the digits array for this detector
  //
  if (fDigits) fDigits->Clear();
}

//_____________________________________________________________________________
void AliT0::SetTreeAddress()
{

  TTree    *treeH = fLoader->TreeH();
  
  if (treeH)
    {
      if (fHits == 0x0) fHits  = new TClonesArray("AliT0hit",  405);
    }
    
  AliDetector::SetTreeAddress();
  TTree *treeD = fLoader->TreeD();
  if (treeD) {
    if (fDigits == 0x0)  fDigits  = new AliT0digit();
    TBranch* branch = treeD->GetBranch ("T0");
    if (branch) branch->SetAddress(&fDigits);
  }

  TTree *treeR = fLoader->TreeR();
  if (treeR) {
    if (fRecPoints == 0x0) fRecPoints  = new  AliT0RecPoint()  ;
    TBranch* branch = treeR->GetBranch ("T0");
    if (branch) branch->SetAddress(&fRecPoints);
  }
  // SDigitizer for Federico
  TTree *treeS = fLoader->TreeS();
  if (treeS) {
    if (fDigits == 0x0)  fDigits  = new AliT0digit();
    TBranch* branch = treeS->GetBranch ("T0");
    if (branch) branch->SetAddress(&fDigits);
  }
 
}


//_____________________________________________________________________________
void AliT0::MakeBranchInTreeD(TTree *treeD, const char *file)
{
    //
    // Create TreeD branches for the FMD
    //
    const Int_t kBufferSize = 4000;
    //   char branchname[20];
    // strncpy(branchname, GetName(), 20);
    TString branchname = Form("%s", GetName());
 //   sprintf(branchname,"%s",GetName());
    if(treeD)
     {
       MakeBranchInTree(treeD,  branchname.Data(),&fDigits, kBufferSize, file);
     }
}

//_____________________________________________________________________________
AliDigitizer* AliT0::CreateDigitizer(AliDigitizationInput* digInput) const
{
  return new AliT0Digitizer(digInput);
}
//____________________________________________________________________________
void AliT0::Digits2Raw()
{
//
// Starting from the T0 digits, writes the Raw Data objects
//
//  AliT0Loader* pStartLoader = (AliT0Loader*)fLoader;
  fLoader ->LoadDigits("read");
  TTree* treeD = fLoader->TreeD();
  if (!treeD) {
    AliError("no digits tree");
    return;
  }
  if (fDigits == 0x0)  fDigits  = new AliT0digit();
  
  TBranch *branch = treeD->GetBranch("T0");
  if (branch) {
    branch->SetAddress(&fDigits);
  }else{
    AliError("Branch T0 DIGIT not found");
    exit(111);
  } 
  AliT0RawData rawWriter;
  rawWriter.SetVerbose(0);
  
  AliDebug(2,Form(" Formatting raw data for T0 "));
  branch->GetEntry(0);
  //  rawWriter.RawDataT0(treeD->GetBranch("T0"));
  rawWriter.RawDataT0(fDigits);
  
  
  fLoader->UnloadDigits();
  
}

//____________________________________________________________________________
void AliT0::Raw2Digits(AliRawReader *rawReader,TTree* digitsTree)
{

 //T0 raw data-> digits conversion
 // reconstruct time information from raw data
 // cout<<"  AliT0::Raw2Digits(AliRawReader *rawReader,TTree* digitsTree) "<<
  // rawReader<<" "<<digitsTree<<endl;

 
  //  AliT0RawReader myrawreader(rawReader,digitsTree);
   AliT0RawReader myrawreader(rawReader);
   if (!myrawreader.Next())
     AliDebug(1,Form(" no raw data found!! %i", myrawreader.Next()));
   Int_t allData[110][5];
   for (Int_t i=0; i<110; i++) {
     allData[i][0]=myrawreader.GetData(i,0);
   }


   fDigits = new AliT0digit();
   digitsTree->Branch("T0","AliT0digit",&fDigits);
   
   
   TArrayI *timeLED = new TArrayI(24);
   TArrayI * timeCFD = new TArrayI(24);
   TArrayI *chargeQT0 = new TArrayI(24);
   TArrayI *chargeQT1 = new TArrayI(24);
   
   for (Int_t in=0; in<24; in++)
     {
       timeLED->AddAt(allData[in+1][0],in);
       timeCFD->AddAt(allData[in+25][0],in);
       chargeQT0->AddAt(allData[in+55][0],in);
       chargeQT1->AddAt(allData[in+79][0],in);
       AliDebug(2, Form(" readed Raw %i %i %i %i %i", in, timeLED->At(in),timeCFD->At(in),chargeQT0->At(in),chargeQT1->At(in)));
     }
  
   fDigits->SetTimeCFD(*timeCFD);
   fDigits->SetQT0(*chargeQT1);

   fDigits->SetTimeLED(*timeLED);
   fDigits->SetQT1(*chargeQT1);

   fDigits->SetMeanTime(allData[49][0]);
   fDigits->SetDiffTime(allData[50][0]);
   fDigits->SetTimeBestA(allData[51][0]);
   fDigits->SetTimeBestC(allData[52][0]);
   digitsTree->Fill();
   fDigits->Write();
 
   delete timeCFD ;
   delete chargeQT0;
   delete timeLED ;
   delete chargeQT1;


}
 AliT0.cxx:1
 AliT0.cxx:2
 AliT0.cxx:3
 AliT0.cxx:4
 AliT0.cxx:5
 AliT0.cxx:6
 AliT0.cxx:7
 AliT0.cxx:8
 AliT0.cxx:9
 AliT0.cxx:10
 AliT0.cxx:11
 AliT0.cxx:12
 AliT0.cxx:13
 AliT0.cxx:14
 AliT0.cxx:15
 AliT0.cxx:16
 AliT0.cxx:17
 AliT0.cxx:18
 AliT0.cxx:19
 AliT0.cxx:20
 AliT0.cxx:21
 AliT0.cxx:22
 AliT0.cxx:23
 AliT0.cxx:24
 AliT0.cxx:25
 AliT0.cxx:26
 AliT0.cxx:27
 AliT0.cxx:28
 AliT0.cxx:29
 AliT0.cxx:30
 AliT0.cxx:31
 AliT0.cxx:32
 AliT0.cxx:33
 AliT0.cxx:34
 AliT0.cxx:35
 AliT0.cxx:36
 AliT0.cxx:37
 AliT0.cxx:38
 AliT0.cxx:39
 AliT0.cxx:40
 AliT0.cxx:41
 AliT0.cxx:42
 AliT0.cxx:43
 AliT0.cxx:44
 AliT0.cxx:45
 AliT0.cxx:46
 AliT0.cxx:47
 AliT0.cxx:48
 AliT0.cxx:49
 AliT0.cxx:50
 AliT0.cxx:51
 AliT0.cxx:52
 AliT0.cxx:53
 AliT0.cxx:54
 AliT0.cxx:55
 AliT0.cxx:56
 AliT0.cxx:57
 AliT0.cxx:58
 AliT0.cxx:59
 AliT0.cxx:60
 AliT0.cxx:61
 AliT0.cxx:62
 AliT0.cxx:63
 AliT0.cxx:64
 AliT0.cxx:65
 AliT0.cxx:66
 AliT0.cxx:67
 AliT0.cxx:68
 AliT0.cxx:69
 AliT0.cxx:70
 AliT0.cxx:71
 AliT0.cxx:72
 AliT0.cxx:73
 AliT0.cxx:74
 AliT0.cxx:75
 AliT0.cxx:76
 AliT0.cxx:77
 AliT0.cxx:78
 AliT0.cxx:79
 AliT0.cxx:80
 AliT0.cxx:81
 AliT0.cxx:82
 AliT0.cxx:83
 AliT0.cxx:84
 AliT0.cxx:85
 AliT0.cxx:86
 AliT0.cxx:87
 AliT0.cxx:88
 AliT0.cxx:89
 AliT0.cxx:90
 AliT0.cxx:91
 AliT0.cxx:92
 AliT0.cxx:93
 AliT0.cxx:94
 AliT0.cxx:95
 AliT0.cxx:96
 AliT0.cxx:97
 AliT0.cxx:98
 AliT0.cxx:99
 AliT0.cxx:100
 AliT0.cxx:101
 AliT0.cxx:102
 AliT0.cxx:103
 AliT0.cxx:104
 AliT0.cxx:105
 AliT0.cxx:106
 AliT0.cxx:107
 AliT0.cxx:108
 AliT0.cxx:109
 AliT0.cxx:110
 AliT0.cxx:111
 AliT0.cxx:112
 AliT0.cxx:113
 AliT0.cxx:114
 AliT0.cxx:115
 AliT0.cxx:116
 AliT0.cxx:117
 AliT0.cxx:118
 AliT0.cxx:119
 AliT0.cxx:120
 AliT0.cxx:121
 AliT0.cxx:122
 AliT0.cxx:123
 AliT0.cxx:124
 AliT0.cxx:125
 AliT0.cxx:126
 AliT0.cxx:127
 AliT0.cxx:128
 AliT0.cxx:129
 AliT0.cxx:130
 AliT0.cxx:131
 AliT0.cxx:132
 AliT0.cxx:133
 AliT0.cxx:134
 AliT0.cxx:135
 AliT0.cxx:136
 AliT0.cxx:137
 AliT0.cxx:138
 AliT0.cxx:139
 AliT0.cxx:140
 AliT0.cxx:141
 AliT0.cxx:142
 AliT0.cxx:143
 AliT0.cxx:144
 AliT0.cxx:145
 AliT0.cxx:146
 AliT0.cxx:147
 AliT0.cxx:148
 AliT0.cxx:149
 AliT0.cxx:150
 AliT0.cxx:151
 AliT0.cxx:152
 AliT0.cxx:153
 AliT0.cxx:154
 AliT0.cxx:155
 AliT0.cxx:156
 AliT0.cxx:157
 AliT0.cxx:158
 AliT0.cxx:159
 AliT0.cxx:160
 AliT0.cxx:161
 AliT0.cxx:162
 AliT0.cxx:163
 AliT0.cxx:164
 AliT0.cxx:165
 AliT0.cxx:166
 AliT0.cxx:167
 AliT0.cxx:168
 AliT0.cxx:169
 AliT0.cxx:170
 AliT0.cxx:171
 AliT0.cxx:172
 AliT0.cxx:173
 AliT0.cxx:174
 AliT0.cxx:175
 AliT0.cxx:176
 AliT0.cxx:177
 AliT0.cxx:178
 AliT0.cxx:179
 AliT0.cxx:180
 AliT0.cxx:181
 AliT0.cxx:182
 AliT0.cxx:183
 AliT0.cxx:184
 AliT0.cxx:185
 AliT0.cxx:186
 AliT0.cxx:187
 AliT0.cxx:188
 AliT0.cxx:189
 AliT0.cxx:190
 AliT0.cxx:191
 AliT0.cxx:192
 AliT0.cxx:193
 AliT0.cxx:194
 AliT0.cxx:195
 AliT0.cxx:196
 AliT0.cxx:197
 AliT0.cxx:198
 AliT0.cxx:199
 AliT0.cxx:200
 AliT0.cxx:201
 AliT0.cxx:202
 AliT0.cxx:203
 AliT0.cxx:204
 AliT0.cxx:205
 AliT0.cxx:206
 AliT0.cxx:207
 AliT0.cxx:208
 AliT0.cxx:209
 AliT0.cxx:210
 AliT0.cxx:211
 AliT0.cxx:212
 AliT0.cxx:213
 AliT0.cxx:214
 AliT0.cxx:215
 AliT0.cxx:216
 AliT0.cxx:217
 AliT0.cxx:218
 AliT0.cxx:219
 AliT0.cxx:220
 AliT0.cxx:221
 AliT0.cxx:222
 AliT0.cxx:223
 AliT0.cxx:224
 AliT0.cxx:225
 AliT0.cxx:226
 AliT0.cxx:227
 AliT0.cxx:228
 AliT0.cxx:229
 AliT0.cxx:230
 AliT0.cxx:231
 AliT0.cxx:232
 AliT0.cxx:233
 AliT0.cxx:234
 AliT0.cxx:235
 AliT0.cxx:236
 AliT0.cxx:237
 AliT0.cxx:238
 AliT0.cxx:239
 AliT0.cxx:240
 AliT0.cxx:241
 AliT0.cxx:242
 AliT0.cxx:243
 AliT0.cxx:244
 AliT0.cxx:245
 AliT0.cxx:246
 AliT0.cxx:247
 AliT0.cxx:248
 AliT0.cxx:249
 AliT0.cxx:250
 AliT0.cxx:251
 AliT0.cxx:252
 AliT0.cxx:253
 AliT0.cxx:254
 AliT0.cxx:255
 AliT0.cxx:256
 AliT0.cxx:257
 AliT0.cxx:258
 AliT0.cxx:259
 AliT0.cxx:260
 AliT0.cxx:261
 AliT0.cxx:262
 AliT0.cxx:263
 AliT0.cxx:264
 AliT0.cxx:265
 AliT0.cxx:266
 AliT0.cxx:267
 AliT0.cxx:268
 AliT0.cxx:269
 AliT0.cxx:270
 AliT0.cxx:271
 AliT0.cxx:272
 AliT0.cxx:273
 AliT0.cxx:274
 AliT0.cxx:275
 AliT0.cxx:276
 AliT0.cxx:277
 AliT0.cxx:278
 AliT0.cxx:279
 AliT0.cxx:280
 AliT0.cxx:281
 AliT0.cxx:282
 AliT0.cxx:283
 AliT0.cxx:284
 AliT0.cxx:285
 AliT0.cxx:286
 AliT0.cxx:287
 AliT0.cxx:288
 AliT0.cxx:289
 AliT0.cxx:290
 AliT0.cxx:291
 AliT0.cxx:292
 AliT0.cxx:293
 AliT0.cxx:294
 AliT0.cxx:295
 AliT0.cxx:296
 AliT0.cxx:297
 AliT0.cxx:298
 AliT0.cxx:299
 AliT0.cxx:300
 AliT0.cxx:301
 AliT0.cxx:302
 AliT0.cxx:303
 AliT0.cxx:304
 AliT0.cxx:305
 AliT0.cxx:306
 AliT0.cxx:307
 AliT0.cxx:308
 AliT0.cxx:309
 AliT0.cxx:310
 AliT0.cxx:311
 AliT0.cxx:312
 AliT0.cxx:313
 AliT0.cxx:314
 AliT0.cxx:315
 AliT0.cxx:316
 AliT0.cxx:317
 AliT0.cxx:318
 AliT0.cxx:319
 AliT0.cxx:320
 AliT0.cxx:321
 AliT0.cxx:322
 AliT0.cxx:323
 AliT0.cxx:324
 AliT0.cxx:325
 AliT0.cxx:326
 AliT0.cxx:327
 AliT0.cxx:328
 AliT0.cxx:329
 AliT0.cxx:330
 AliT0.cxx:331
 AliT0.cxx:332
 AliT0.cxx:333
 AliT0.cxx:334
 AliT0.cxx:335
 AliT0.cxx:336
 AliT0.cxx:337
 AliT0.cxx:338
 AliT0.cxx:339
 AliT0.cxx:340
 AliT0.cxx:341
 AliT0.cxx:342
 AliT0.cxx:343
 AliT0.cxx:344
 AliT0.cxx:345
 AliT0.cxx:346
 AliT0.cxx:347
 AliT0.cxx:348
 AliT0.cxx:349
 AliT0.cxx:350
 AliT0.cxx:351
 AliT0.cxx:352
 AliT0.cxx:353
 AliT0.cxx:354
 AliT0.cxx:355
 AliT0.cxx:356
 AliT0.cxx:357
 AliT0.cxx:358
 AliT0.cxx:359
 AliT0.cxx:360
 AliT0.cxx:361
 AliT0.cxx:362
 AliT0.cxx:363
 AliT0.cxx:364
 AliT0.cxx:365
 AliT0.cxx:366
 AliT0.cxx:367
 AliT0.cxx:368
 AliT0.cxx:369
 AliT0.cxx:370
 AliT0.cxx:371
 AliT0.cxx:372
 AliT0.cxx:373
 AliT0.cxx:374
 AliT0.cxx:375
 AliT0.cxx:376
 AliT0.cxx:377
 AliT0.cxx:378
 AliT0.cxx:379
 AliT0.cxx:380
 AliT0.cxx:381
 AliT0.cxx:382
 AliT0.cxx:383
 AliT0.cxx:384
 AliT0.cxx:385
 AliT0.cxx:386
 AliT0.cxx:387
 AliT0.cxx:388
 AliT0.cxx:389
 AliT0.cxx:390
 AliT0.cxx:391