[u/mrichter/AliRoot.git] / PHOS / Reconstruct.C

// YS Subatech Mai 2002

//Root
#include "TString.h"

//AliRoot
#include "STEER/AliRun.h"
#include "PHOS/AliPHOSSDigitizer.h"
#include "PHOS/AliPHOSDigitizer.h"
#include "PHOS/AliPHOSClusterizerv1.h"
#include "PHOS/AliPHOSTrackSegmentMakerv1.h"
#include "PHOS/AliPHOSPIDv1.h"
#include "EMCAL/AliEMCALSDigitizer.h"
#include "EMCAL/AliEMCALDigitizer.h"
#include "EMCAL/AliEMCALClusterizerv1.h"

void Hits2SDigits( Bool_t split=kFALSE, TString fileName = "galice.root") {

  // usage : 
  // 1. write SDigits in the same file as Hits --------------- (OK)
  //root [0] .L Reconstruct.C++
  //root [1] SDigits2Digits()
  // 2. write SDigits in a separate file, one per detector, from Hits --------------- (OK)
  //root [0] .L Reconstruct.C++
  //root [1] SDigits2Digits(kTRUE) // SDigits saved in [DET}.SDigits.root (DET=PHOS, EMCAL)

  delete gAlice ; 
  gAlice = 0 ; 
  
  AliPHOSSDigitizer * sdp = new AliPHOSSDigitizer(fileName) ; 
  if (split) 
    sdp->SetSplitFile() ;
  sdp->ExecuteTask("deb") ; 

  delete sdp ;

 //  AliEMCALSDigitizer * sde = new AliEMCALSDigitizer(fileName) ; 
//   if (split) 
//     sde->SetSplitFile() ;
//   sde->ExecuteTask("deb") ; 
 
//   delete sde ; 

}

//________________________________________________________________________
void SDigits2Digits( Bool_t split=kFALSE, TString fileName = "galice.root") {
  
 // usage : 
  // 1. write SDigits in the same file as SDigits --------------- (OK)
  //root [0] .L Reconstruct.C++
  //root [1] SDigits2Digits()
  // 2. write SDigits in a separate file, one per detector, from SDigits --------------- (OK)
  //root [0] .L Reconstruct.C++
  //root [1] SDigitsDigits(kTRUE) // Digits saved in [DET}.Digits.root (DET=PHOS, EMCAL)

  delete gAlice ; 
  gAlice = 0 ; 
  
  // PHOS
  AliPHOSDigitizer * dp = 0 ; 
 
  if (split) {
    dp = new AliPHOSDigitizer("PHOS.SDigits.root") ; 
    dp->SetSplitFile() ; } 
  else 
    dp = new AliPHOSDigitizer(fileName) ; 
  
  dp->ExecuteTask("deb") ; 
  
  delete dp ;

  //EMCAL
//   AliEMCALDigitizer * de = 0 ; 

//   if (split) {
//     de = new AliEMCALDigitizer("EMCAL.SDigits.root") ;
//     de->SetSplitFile() ;
//   } else 
//     de = new AliEMCALDigitizer(fileName) ; 
  
//   de->ExecuteTask("deb") ; 
  
//   delete de ; 
}

//________________________________________________________________________
void Digits2RecPoints( Bool_t split=kFALSE, TString fileName = "galice.root") {
  
 // usage : 
  // 1. write RecPoints in the same file as Digits --------------- OK 
  //root [0] .L Reconstruct.C++
  //root [1] Digits2RecPoints()
  // 2. write RecPoints in a separate file, one per detector, from Digits --------------- OK 
  //root [0] .L Reconstruct.C++
  //root [1] Digits2RecPoints(kTRUE) // RecPoints saved in [DET}.RecPoints.root (DET=PHOS, EMCAL)

  delete gAlice ; 
  gAlice = 0 ; 
 
// PHOS
  AliPHOSClusterizer * cp = 0 ; 
 
  if (split) {
    cp = new AliPHOSClusterizerv1("PHOS.Digits.root") ; 
    cp->SetSplitFile() ; } 
  else 
    cp = new AliPHOSClusterizerv1(fileName) ; 
  
  cp->ExecuteTask("deb") ; 
  
  delete cp ;

  //EMCAL
  AliEMCALClusterizerv1 * ce = 0 ;  

  if (split) {
    ce = new AliEMCALClusterizerv1("EMCAL.Digits.root") ;
    ce->SetSplitFile() ;
  } else 
    ce = new AliEMCALClusterizerv1(fileName) ; 
  
  ce->ExecuteTask("deb") ; 
  
  delete ce ; 
}

//________________________________________________________________________
void RecPoints2TrackSegments( Bool_t split=kFALSE, TString fileName = "galice.root") {
  
 // usage : 
  // 1. write TrackSegments in the same file as RecPoints --------------- (OK) 
  //root [0] .L Reconstruct.C++
  //root [1] RecPoints2TrackSegments()
  // 2. write TrackSegments in a separate file, one per detector, from RecPoints --------------- (Not needed) 
  //root [0] .L Reconstruct.C++
  //root [1] RecPoints2TrackSegments(kTRUE) // TrackSegments saved in [DET}.RecData.root (DET=PHOS, EMCAL)

  delete gAlice ; 
  gAlice = 0 ; 
  
  AliPHOSTrackSegmentMaker * tmp = 0 ; 
 
  if (split) {
    tmp = new AliPHOSTrackSegmentMakerv1("PHOS.RecData.root") ; 
    tmp->SetSplitFile() ; } 
  else 
    tmp = new AliPHOSTrackSegmentMakerv1(fileName) ; 
  
  tmp->ExecuteTask("deb") ; 
  
  delete tmp ;
}

//________________________________________________________________________
void TrackSegments2RecParticles( Bool_t split=kFALSE, TString fileName = "galice.root") {
  
 // usage : 
  // 1. write RecParticles in the same file as TrackSegments ---------------  (OK)
  //root [0] .L Reconstruct.C++
  //root [1] TrackSegments2RecParticles()
  // 2. write RecParticles in a separate file, one per detector, from TrackSegments --------------- (Not needed) 
  //root [0] .L Reconstruct.C++
  //root [1] TrackSegments2RecParticles(kTRUE) // RecParticles saved in [DET}.RecData.root (DET=PHOS, EMCAL)

  delete gAlice ; 
  gAlice = 0 ; 
  
  AliPHOSPID * pp = 0 ; 
 
  if (split) 
    pp = new AliPHOSPIDv1("PHOS.RecData.root") ; 
  else 
    pp = new AliPHOSPIDv1(fileName) ; 
  
  pp->ExecuteTask("deb") ; 
  
  delete pp ;
}

//________________________________________________________________________
void Digits2RecParticles( Bool_t split=kFALSE, TString fileName = "galice.root") {
  
 // usage : 
  // 1. write RecPoints, TrackSegments and RecParticles in the same file as Digits --------------- (OK)
  //root [0] .L Reconstruct.C++
  //root [1] Digits2RecParticles()
  // 2. write RecPoints , TrackSegments and RecParticles in a separate file, one per detector, from Digits --------------- (OK)
  //root [0] .L Reconstruct.C++
  //root [1] Digits2RecParticles(kTRUE) // TrackSegments saved in [DET}.RecData.root (DET=PHOS, EMCAL)

 
  delete gAlice ; 
  gAlice = 0 ; 
 
  // PHOS
  AliPHOSClusterizer * cp = 0 ; 
 
  if (split) {
    cp = new AliPHOSClusterizerv1("PHOS.Digits.root") ; 
    cp->SetSplitFile() ; } 
  else 
    cp = new AliPHOSClusterizerv1(fileName) ; 
  
  cp->ExecuteTask("deb") ; 
  
  if (split) 
    delete cp ;
  
  AliPHOSTrackSegmentMaker * tmp = 0 ; 
  
  if (split) 
    tmp = new AliPHOSTrackSegmentMakerv1("PHOS.RecData.root") ; 
  else 
    tmp = new AliPHOSTrackSegmentMakerv1(fileName) ; 
  
  tmp->ExecuteTask("deb") ; 
  
  AliPHOSPID * pp = 0 ; 
 
  if (split) 
    pp = new AliPHOSPIDv1("PHOS.RecData.root") ; 
  else 
    pp = new AliPHOSPIDv1(fileName) ; 
  
  pp->ExecuteTask("deb") ; 
  
    delete cp ; 
    delete tmp; 
    delete pp ; 
}

//________________________________________________________________________
void Hits2Digits (Bool_t split=kFALSE, TString fileName = "galice.root") {
  // usage : 
  // 1. write (S)Digits in the same file as Hits --------------- (OK)
  //root [0] .L Reconstruct.C++
  //root [1] Hits2Digits()
  // 2. write (S)Digits in a separate file, one per detector, from Hits --------------- (OK)
  //root [0] .L Reconstruct.C++
  //root [1] Hits2Digits(kTRUE) // SDigits saved in [DET}.SDigits.root (DET=PHOS, EMCAL)
                                // Digits  saved in [DET}.Digits.root  (DET=PHOS, EMCAL)

  delete gAlice ; 
  gAlice = 0 ; 
  
  //PHOS
  AliPHOSSDigitizer * sdp = new AliPHOSSDigitizer(fileName) ; 
  if (split) 
    sdp->SetSplitFile() ;
  sdp->ExecuteTask("deb") ; 

  if (split) 
    delete sdp ; 

  AliPHOSDigitizer * dp = 0 ; 
 
  if (split) {
    dp = new AliPHOSDigitizer("PHOS.SDigits.root") ; 
    dp->SetSplitFile() ; } 
  else 
    dp = new AliPHOSDigitizer(fileName) ; 
  
  dp->ExecuteTask("deb") ; 

  if (split) 
    delete dp ; 

  if (!split) { 
    delete sdp ; 
    delete dp ; 
  }
  
//   //EMCAL
//   AliEMCALSDigitizer * sde = new AliEMCALSDigitizer(fileName) ; 
//   if (split) 
//     sde->SetSplitFile() ;
//   sde->ExecuteTask("deb") ; 
  
//   delete sde ; 
  
//   AliEMCALDigitizer * de = 0 ; 
//   if (split) {
//     de = new AliEMCALDigitizer("EMCAL.SDigits.root") ;
//     de->SetSplitFile() ;
//   } else 
//     de = new AliEMCALDigitizer(fileName) ; 
  
//   de->ExecuteTask("deb") ; 
  
//   delete de ; 

}
Commit	Line	Data
aca467af	1	// YS Subatech Mai 2002
	2
	3	//Root
	4	#include "TString.h"
	5
	6	//AliRoot
8d0f3f77	7	#include "STEER/AliRun.h"
aca467af	8	#include "PHOS/AliPHOSSDigitizer.h"
aca467af	9	#include "PHOS/AliPHOSDigitizer.h"
8d0f3f77	10	#include "PHOS/AliPHOSClusterizerv1.h"
	11	#include "PHOS/AliPHOSTrackSegmentMakerv1.h"
	12	#include "PHOS/AliPHOSPIDv1.h"
aca467af	13	#include "EMCAL/AliEMCALSDigitizer.h"
aca467af	14	#include "EMCAL/AliEMCALDigitizer.h"
8d0f3f77	15	#include "EMCAL/AliEMCALClusterizerv1.h"
aca467af	16
	17	void Hits2SDigits( Bool_t split=kFALSE, TString fileName = "galice.root") {
	18
	19	// usage :
	20	// 1. write SDigits in the same file as Hits --------------- (OK)
	21	//root [0] .L Reconstruct.C++
	22	//root [1] SDigits2Digits()
	23	// 2. write SDigits in a separate file, one per detector, from Hits --------------- (OK)
	24	//root [0] .L Reconstruct.C++
	25	//root [1] SDigits2Digits(kTRUE) // SDigits saved in [DET}.SDigits.root (DET=PHOS, EMCAL)
	26
8d0f3f77	27	delete gAlice ;
	28	gAlice = 0 ;
	29
aca467af	30	AliPHOSSDigitizer * sdp = new AliPHOSSDigitizer(fileName) ;
	31	if (split)
	32	sdp->SetSplitFile() ;
	33	sdp->ExecuteTask("deb") ;
	34
e9fa1121	35	delete sdp ;
e9fa1121	36
8d0f3f77	37	// AliEMCALSDigitizer * sde = new AliEMCALSDigitizer(fileName) ;
	38	// if (split)
	39	// sde->SetSplitFile() ;
	40	// sde->ExecuteTask("deb") ;
aca467af	41
8d0f3f77	42	// delete sde ;
aca467af	43
	44	}
	45
	46	//________________________________________________________________________
	47	void SDigits2Digits( Bool_t split=kFALSE, TString fileName = "galice.root") {
	48
	49	// usage :
8d0f3f77	50	// 1. write SDigits in the same file as SDigits --------------- (OK)
aca467af	51	//root [0] .L Reconstruct.C++
8d0f3f77	52	//root [1] SDigits2Digits()
8d0f3f77	53	// 2. write SDigits in a separate file, one per detector, from SDigits --------------- (OK)
aca467af	54	//root [0] .L Reconstruct.C++
8d0f3f77	55	//root [1] SDigitsDigits(kTRUE) // Digits saved in [DET}.Digits.root (DET=PHOS, EMCAL)
aca467af	56
8d0f3f77	57	delete gAlice ;
	58	gAlice = 0 ;
	59
e9fa1121	60	// PHOS
8d0f3f77	61	AliPHOSDigitizer * dp = 0 ;
e9fa1121	62
aca467af	63	if (split) {
	64	dp = new AliPHOSDigitizer("PHOS.SDigits.root") ;
	65	dp->SetSplitFile() ; }
	66	else
	67	dp = new AliPHOSDigitizer(fileName) ;
	68
	69	dp->ExecuteTask("deb") ;
	70
e9fa1121	71	delete dp ;
	72
	73	//EMCAL
8d0f3f77	74	// AliEMCALDigitizer * de = 0 ;
	75
	76	// if (split) {
	77	// de = new AliEMCALDigitizer("EMCAL.SDigits.root") ;
	78	// de->SetSplitFile() ;
	79	// } else
	80	// de = new AliEMCALDigitizer(fileName) ;
	81
	82	// de->ExecuteTask("deb") ;
	83
	84	// delete de ;
	85	}
	86
	87	//________________________________________________________________________
	88	void Digits2RecPoints( Bool_t split=kFALSE, TString fileName = "galice.root") {
	89
	90	// usage :
	91	// 1. write RecPoints in the same file as Digits --------------- OK
	92	//root [0] .L Reconstruct.C++
	93	//root [1] Digits2RecPoints()
	94	// 2. write RecPoints in a separate file, one per detector, from Digits --------------- OK
	95	//root [0] .L Reconstruct.C++
	96	//root [1] Digits2RecPoints(kTRUE) // RecPoints saved in [DET}.RecPoints.root (DET=PHOS, EMCAL)
e9fa1121	97
8d0f3f77	98	delete gAlice ;
	99	gAlice = 0 ;
	100
	101	// PHOS
	102	AliPHOSClusterizer * cp = 0 ;
	103
aca467af	104	if (split) {
8d0f3f77	105	cp = new AliPHOSClusterizerv1("PHOS.Digits.root") ;
	106	cp->SetSplitFile() ; }
	107	else
	108	cp = new AliPHOSClusterizerv1(fileName) ;
	109
	110	cp->ExecuteTask("deb") ;
	111
	112	delete cp ;
	113
	114	//EMCAL
	115	AliEMCALClusterizerv1 * ce = 0 ;
	116
	117	if (split) {
	118	ce = new AliEMCALClusterizerv1("EMCAL.Digits.root") ;
	119	ce->SetSplitFile() ;
aca467af	120	} else
8d0f3f77	121	ce = new AliEMCALClusterizerv1(fileName) ;
	122
	123	ce->ExecuteTask("deb") ;
	124
	125	delete ce ;
	126	}
	127
	128	//________________________________________________________________________
	129	void RecPoints2TrackSegments( Bool_t split=kFALSE, TString fileName = "galice.root") {
	130
	131	// usage :
	132	// 1. write TrackSegments in the same file as RecPoints --------------- (OK)
	133	//root [0] .L Reconstruct.C++
	134	//root [1] RecPoints2TrackSegments()
	135	// 2. write TrackSegments in a separate file, one per detector, from RecPoints --------------- (Not needed)
	136	//root [0] .L Reconstruct.C++
	137	//root [1] RecPoints2TrackSegments(kTRUE) // TrackSegments saved in [DET}.RecData.root (DET=PHOS, EMCAL)
	138
	139	delete gAlice ;
	140	gAlice = 0 ;
	141
	142	AliPHOSTrackSegmentMaker * tmp = 0 ;
	143
	144	if (split) {
	145	tmp = new AliPHOSTrackSegmentMakerv1("PHOS.RecData.root") ;
	146	tmp->SetSplitFile() ; }
	147	else
	148	tmp = new AliPHOSTrackSegmentMakerv1(fileName) ;
	149
	150	tmp->ExecuteTask("deb") ;
	151
	152	delete tmp ;
	153	}
	154
	155	//________________________________________________________________________
	156	void TrackSegments2RecParticles( Bool_t split=kFALSE, TString fileName = "galice.root") {
	157
	158	// usage :
	159	// 1. write RecParticles in the same file as TrackSegments --------------- (OK)
	160	//root [0] .L Reconstruct.C++
	161	//root [1] TrackSegments2RecParticles()
	162	// 2. write RecParticles in a separate file, one per detector, from TrackSegments --------------- (Not needed)
	163	//root [0] .L Reconstruct.C++
	164	//root [1] TrackSegments2RecParticles(kTRUE) // RecParticles saved in [DET}.RecData.root (DET=PHOS, EMCAL)
	165
	166	delete gAlice ;
	167	gAlice = 0 ;
	168
	169	AliPHOSPID * pp = 0 ;
	170
	171	if (split)
	172	pp = new AliPHOSPIDv1("PHOS.RecData.root") ;
	173	else
	174	pp = new AliPHOSPIDv1(fileName) ;
e9fa1121	175
8d0f3f77	176	pp->ExecuteTask("deb") ;
e9fa1121	177
8d0f3f77	178	delete pp ;
	179	}
	180
	181	//________________________________________________________________________
	182	void Digits2RecParticles( Bool_t split=kFALSE, TString fileName = "galice.root") {
	183
	184	// usage :
	185	// 1. write RecPoints, TrackSegments and RecParticles in the same file as Digits --------------- (OK)
	186	//root [0] .L Reconstruct.C++
	187	//root [1] Digits2RecParticles()
	188	// 2. write RecPoints , TrackSegments and RecParticles in a separate file, one per detector, from Digits --------------- (OK)
	189	//root [0] .L Reconstruct.C++
	190	//root [1] Digits2RecParticles(kTRUE) // TrackSegments saved in [DET}.RecData.root (DET=PHOS, EMCAL)
	191
	192
	193	delete gAlice ;
	194	gAlice = 0 ;
	195
	196	// PHOS
	197	AliPHOSClusterizer * cp = 0 ;
	198
	199	if (split) {
	200	cp = new AliPHOSClusterizerv1("PHOS.Digits.root") ;
	201	cp->SetSplitFile() ; }
	202	else
	203	cp = new AliPHOSClusterizerv1(fileName) ;
	204
	205	cp->ExecuteTask("deb") ;
	206
	207	if (split)
	208	delete cp ;
	209
	210	AliPHOSTrackSegmentMaker * tmp = 0 ;
	211
	212	if (split)
	213	tmp = new AliPHOSTrackSegmentMakerv1("PHOS.RecData.root") ;
	214	else
	215	tmp = new AliPHOSTrackSegmentMakerv1(fileName) ;
	216
	217	tmp->ExecuteTask("deb") ;
	218
	219	AliPHOSPID * pp = 0 ;
	220
	221	if (split)
	222	pp = new AliPHOSPIDv1("PHOS.RecData.root") ;
	223	else
	224	pp = new AliPHOSPIDv1(fileName) ;
	225
	226	pp->ExecuteTask("deb") ;
	227
	228	delete cp ;
	229	delete tmp;
	230	delete pp ;
aca467af	231	}
	232
	233	//________________________________________________________________________
	234	void Hits2Digits (Bool_t split=kFALSE, TString fileName = "galice.root") {
	235	// usage :
	236	// 1. write (S)Digits in the same file as Hits --------------- (OK)
	237	//root [0] .L Reconstruct.C++
	238	//root [1] Hits2Digits()
	239	// 2. write (S)Digits in a separate file, one per detector, from Hits --------------- (OK)
	240	//root [0] .L Reconstruct.C++
	241	//root [1] Hits2Digits(kTRUE) // SDigits saved in [DET}.SDigits.root (DET=PHOS, EMCAL)
	242	// Digits saved in [DET}.Digits.root (DET=PHOS, EMCAL)
	243
8d0f3f77	244	delete gAlice ;
	245	gAlice = 0 ;
	246
e9fa1121	247	//PHOS
	248	AliPHOSSDigitizer * sdp = new AliPHOSSDigitizer(fileName) ;
	249	if (split)
	250	sdp->SetSplitFile() ;
	251	sdp->ExecuteTask("deb") ;
	252
8d0f3f77	253	if (split)
8d0f3f77	254	delete sdp ;
e9fa1121	255
8d0f3f77	256	AliPHOSDigitizer * dp = 0 ;
e9fa1121	257
	258	if (split) {
	259	dp = new AliPHOSDigitizer("PHOS.SDigits.root") ;
	260	dp->SetSplitFile() ; }
	261	else
	262	dp = new AliPHOSDigitizer(fileName) ;
	263
	264	dp->ExecuteTask("deb") ;
e9fa1121	265
e9fa1121	266	if (split)
8d0f3f77	267	delete dp ;
	268
	269	if (!split) {
	270	delete sdp ;
	271	delete dp ;
	272	}
e9fa1121	273
8d0f3f77	274	// //EMCAL
	275	// AliEMCALSDigitizer * sde = new AliEMCALSDigitizer(fileName) ;
	276	// if (split)
	277	// sde->SetSplitFile() ;
	278	// sde->ExecuteTask("deb") ;
e9fa1121	279
8d0f3f77	280	// delete sde ;
	281
	282	// AliEMCALDigitizer * de = 0 ;
	283	// if (split) {
	284	// de = new AliEMCALDigitizer("EMCAL.SDigits.root") ;
	285	// de->SetSplitFile() ;
	286	// } else
	287	// de = new AliEMCALDigitizer(fileName) ;
e9fa1121	288
8d0f3f77	289	// de->ExecuteTask("deb") ;
e9fa1121	290
8d0f3f77	291	// delete de ;
8d0f3f77	292
aca467af	293	}
e9fa1121	294
8d0f3f77	295
e9fa1121	296