[u/mrichter/AliRoot.git] / PHOS / macros / testSuite / TestSuite.csh

#! /bin/tcsh -f
#
# Prepare the environment
#
set TMPDIR=`mktemp -d $HOME/PHOSTestSuite.XXXXXX` ; cd $TMPDIR
set SRCDIR=$ALICE_ROOT/PHOS/macros/testSuite
set LOGFILE=$TMPDIR/TestSuite.log
set ConfigSgnl=$SRCDIR/ConfigTestSuite.C
set ConfigBckg=$SRCDIR/ConfigTestSuiteHij.C
set GRUN=$SRCDIR/grunTestSuite.C

printf "PHOS test suite is being executed in directory\n   $TMPDIR\n"

# Set mail address to send errors to
set PHOS_LIBRARIAN=yves.schutz@cern.ch

#
#begin data
set MINHITS=68
set MAXHITS=72
set MINSDIGITS=68
set MAXSDIGITS=72
set MINDIGITS=47
set MAXDIGITS=51
set MINEMCRECPOINTS=1
set MAXEMCRECPOINTS=2
set MINCPVRECPOINTS=0
set MAXCPVRECPOINTS=2
set MINTRACKSEGMENTS=1
set MAXTRACKSEGMENTS=2
set MINRECPARTICLES=1
set MAXRECPARTICLES=2
#end data

#
echo PHOS Test Suite run on `date` > $LOGFILE
echo "=======================================================================" >> $LOGFILE
echo AliROOT  version `which aliroot` >> $LOGFILE
echo Root     version `root-config --version` in `which root` >> $LOGFILE
echo "=======================================================================" >> $LOGFILE
#
# SIMULATION
#
printf "Signal simulation         ... "
echo "    **** Signal SIMULATION **** " >> $LOGFILE
mkdir -p $TMPDIR/signal
cd       $TMPDIR/signal
aliroot -b -q ${GRUN}\(100,\"${ConfigSgnl}\"\) >>& $LOGFILE 
printf "done\n"

printf "Background simulation     ... "
echo "    **** Background SIMULATION **** " >> $LOGFILE
mkdir -p $TMPDIR/background
cd       $TMPDIR/background
aliroot -b -q ${GRUN}\(2,\"${ConfigBckg}\"\) >>& $LOGFILE 
printf "done\n"

cd $TMPDIR
cp -a $TMPDIR/signal $TMPDIR/signal+background

printf "Simulation analysis       ... "
cd $TMPDIR/signal
echo '{ '  > tempo.C
echo " Int_t minhits  = $MINHITS ;"  >> tempo.C 
echo " Int_t maxhits = $MAXHITS ;"  >> tempo.C 
echo ' AliPHOSGetter * gime = AliPHOSGetter::Instance("galice.root") ;'  >> tempo.C
echo ' TH1F * hitmul = new TH1F("hitmul", "PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
echo ' Int_t max = gime->MaxEvent() ; ' >> tempo.C
echo ' Int_t evt = 0 ; ' >> tempo.C
echo ' for ( evt = 0 ; evt < max ; evt++ ) { ' >> tempo.C
echo '  gime->Event(evt,"H") ; ' >> tempo.C
echo '  hitmul->Fill(gime->Hits()->GetEntries()) ; ' >> tempo.C
echo ' } ' >> tempo.C
echo ' TF1 * gaus = new TF1("gaus", "gaus", 0., 200.) ; ' >> tempo.C 
echo ' hitmul->Fit(gaus,"", "", 40, 100) ; ' >> tempo.C 
echo ' if( gaus->GetParameter(1) > maxhits ||  gaus->GetParameter(1) < minhits ) ' >> tempo.C  
echo '  printf("ERRORSIM 1") ; ' >> tempo.C
echo ' else ' >> tempo.C
echo '  printf("ERRORSIM 0") ; ' >> tempo.C
echo '} '  >> tempo.C
aliroot -b -q tempo.C >>& $LOGFILE 
set ERRORSIM = `cat $LOGFILE | grep ERRORSIM | awk '{print $2}'`
if($ERRORSIM) then
  printf "error\n"
  uuencode  $LOGFILE  $LOGFILE | mail -s 'PHOS INSTALLATION ERROR' $PHOS_LIBRARIAN
##  exit(0)
else
  printf "done\n"
endif 
rm tempo.C 
#
# RECONSTRUCTION
#
printf "Reconstruction            ... "
echo "    **** RECONSTRUCTION **** " >> $LOGFILE
aliroot -b >>& $LOGFILE <<EOF
 AliPHOSSDigitizer sd("galice.root") ; 
 sd.ExecuteTask("deb all") ; 
 AliPHOSDigitizer d("galice.root") ; 
 d.ExecuteTask("deb all") ; 
 AliPHOSClusterizerv1 cl("galice.root") ; 
 cl.ExecuteTask("deb all") ; 
 AliPHOSTrackSegmentMakerv1 ts("galice.root") ; 
 ts.ExecuteTask("deb all") ; 
 AliPHOSPIDv1 pd("galice.root") ; 
 pd.ExecuteTask("deb all") ; 
 .q
EOF
printf "done\n"

printf "Reconstruction analysis   ... "
echo '{ '  > tempo.C
echo " Int_t minsdig  = $MINSDIGITS ;"  >> tempo.C 
echo " Int_t maxsdig  = $MAXSDIGITS ;"  >> tempo.C 
echo " Int_t mindig   = $MINDIGITS ;"  >> tempo.C 
echo " Int_t maxdig   = $MAXDIGITS ;"  >> tempo.C 
echo " Int_t minemcrp = $MINEMCRECPOINTS ;"  >> tempo.C 
echo " Int_t maxemcrp = $MAXEMCRECPOINTS ;"  >> tempo.C 
echo " Int_t mincpvrp = $MINCPVRECPOINTS ;"  >> tempo.C 
echo " Int_t maxcpvrp = $MAXCPVRECPOINTS ;"  >> tempo.C 
echo " Int_t mints    = $MINTRACKSEGMENTS ;"  >> tempo.C 
echo " Int_t maxts    = $MAXTRACKSEGMENTS ;"  >> tempo.C 
echo " Int_t minpa    = $MINRECPARTICLES ;"  >> tempo.C 
echo " Int_t maxpa    = $MAXRECPARTICLES ;"  >> tempo.C 
echo ' AliPHOSGetter * gime = AliPHOSGetter::Instance("galice.root") ;'  >> tempo.C
echo ' TH1F * sdigmul = new TH1F("sdigmul", " SDigits PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
echo ' TH1F * digmul  = new TH1F("digmul", " Digits PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
echo ' TH1F * emcrpmul= new TH1F("emcrpmul", " EMCRecPoints PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
echo ' TH1F * cpvrpmul= new TH1F("cpvrpmul", " CPVRecPoints PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
echo ' TH1F * tsmul   = new TH1F("tsmul", " TrackSegments PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
echo ' TH1F * pamul   = new TH1F("pamul", " RecParticles PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
echo ' Int_t max = gime->MaxEvent() ; ' >> tempo.C
echo ' Int_t evt = 0 ; ' >> tempo.C
echo ' for ( evt = 0 ; evt < max ; evt++ ) { ' >> tempo.C
echo '  gime->Event(evt,"SDRTP") ; ' >> tempo.C
echo '  sdigmul->Fill(gime->SDigits()->GetEntries()) ; ' >> tempo.C
echo '  digmul->Fill(gime->Digits()->GetEntries()) ; ' >> tempo.C
echo '  emcrpmul->Fill(gime->EmcRecPoints()->GetEntries()) ; ' >> tempo.C
echo '  cpvrpmul->Fill(gime->CpvRecPoints()->GetEntries()) ; ' >> tempo.C
echo '  tsmul->Fill(gime->TrackSegments()->GetEntries()) ; ' >> tempo.C
echo '  pamul->Fill(gime->RecParticles()->GetEntries()) ; ' >> tempo.C
echo ' } ' >> tempo.C
echo ' TF1 * gaus = new TF1("gaus", "gaus", 0., 200.) ; ' >> tempo.C 
echo ' sdigmul->Fit(gaus, "", "", 40., 100.) ; ' >> tempo.C 
echo ' sdigmul->Draw() ; ' >> tempo.C 
echo ' if( gaus->GetParameter(1) > maxsdig ||  gaus->GetParameter(1) < minsdig ) ' >> tempo.C  
echo '  printf("ERRORREC 1 sdigits\n") ; ' >> tempo.C
echo ' digmul->Fit(gaus, "", "", 20., 80.) ; ' >> tempo.C 
echo ' if( gaus->GetParameter(1) > maxdig ||  gaus->GetParameter(1) < mindig ) ' >> tempo.C  
echo '  printf("ERRORREC 1 digits\n") ; '  >> tempo.C
echo ' emcrpmul->Fit(gaus, "", "", 0., 4.) ; ' >> tempo.C 
echo ' if( gaus->GetParameter(1) > maxemcrp ||  gaus->GetParameter(1) < minemcrp ) ' >> tempo.C  
echo '  printf("ERRORREC 1 emc recpoints\n") ; ' >> tempo.C
echo ' cpvrpmul->Fit(gaus, "", "", 0., 4.) ; ' >> tempo.C 
echo ' if( gaus->GetParameter(1) > maxcpvrp ||  gaus->GetParameter(1) < mincpvrp ) ' >> tempo.C  
echo '  printf("ERRORREC 1 cpv recpoints\n") ; ' >> tempo.C
echo ' tsmul->Fit(gaus, "", "", 0., 4.) ; ' >> tempo.C 
echo ' if( gaus->GetParameter(1) > maxts ||  gaus->GetParameter(1) < mints ) ' >> tempo.C  
echo '  printf("ERRORREC 1 track segments\n") ; ' >> tempo.C
echo ' pamul->Fit(gaus, "", "", 0., 4. ) ; ' >> tempo.C 
echo ' if( gaus->GetParameter(1) > maxpa ||  gaus->GetParameter(1) < minpa ) ' >> tempo.C  
echo '  printf("ERRORREC 1 recparticles\n") ; ' >> tempo.C
echo '} ' >> tempo.C
aliroot -b -q tempo.C  >>& $LOGFILE 
set ERRORREC = `cat $LOGFILE | grep ERRORREC | awk '{print $2}'`
if( $?ERRORREC ) then
  printf "error\n"
  uuencode  $LOGFILE  $LOGFILE | mail -s 'PHOS INSTALLATION ERROR' $PHOS_LIBRARIAN
##  exit(0)
else
  printf "done\n"
endif 
rm tempo.C 

# Signal+background merging
printf "Signal+background merging ... "
echo "    **** Signal+background MERGING **** " >> $LOGFILE
cd $TMPDIR/background
aliroot -b >>& $LOGFILE <<EOF
 AliPHOSSDigitizer sd("galice.root") ; 
 sd.ExecuteTask("deb all") ; 
EOF

cd $TMPDIR/signal+background
cp -a $TMPDIR/signal/PHOS.SDigits*.root .
aliroot -b >>& $LOGFILE <<EOF
 AliPHOSDigitizer d("galice.root") ; 
 d.MixWith("../background/galice.root") ;
 d.ExecuteTask("deb all") ; 
 AliPHOSClusterizerv1 cl("galice.root") ; 
 cl.ExecuteTask("deb all") ; 
 AliPHOSTrackSegmentMakerv1 ts("galice.root") ; 
 ts.ExecuteTask("deb all") ; 
 AliPHOSPIDv1 pd("galice.root") ; 
 pd.ExecuteTask("deb all") ; 
 .q
EOF
printf "done\n"

#rm -fr $TMPDIR
Commit	Line	Data
be66de47	1	#! /bin/tcsh -f
	2	#
	3	# Prepare the environment
	4	#
	5	set TMPDIR=`mktemp -d $HOME/PHOSTestSuite.XXXXXX` ; cd $TMPDIR
	6	set SRCDIR=$ALICE_ROOT/PHOS/macros/testSuite
	7	set LOGFILE=$TMPDIR/TestSuite.log
	8	set ConfigSgnl=$SRCDIR/ConfigTestSuite.C
	9	set ConfigBckg=$SRCDIR/ConfigTestSuiteHij.C
	10	set GRUN=$SRCDIR/grunTestSuite.C
	11
	12	printf "PHOS test suite is being executed in directory\n $TMPDIR\n"
	13
	14	# Set mail address to send errors to
	15	set PHOS_LIBRARIAN=yves.schutz@cern.ch
	16
	17	#
	18	#begin data
	19	set MINHITS=68
	20	set MAXHITS=72
	21	set MINSDIGITS=68
	22	set MAXSDIGITS=72
	23	set MINDIGITS=47
	24	set MAXDIGITS=51
	25	set MINEMCRECPOINTS=1
	26	set MAXEMCRECPOINTS=2
	27	set MINCPVRECPOINTS=0
	28	set MAXCPVRECPOINTS=2
	29	set MINTRACKSEGMENTS=1
	30	set MAXTRACKSEGMENTS=2
	31	set MINRECPARTICLES=1
	32	set MAXRECPARTICLES=2
	33	#end data
	34
	35	#
	36	echo PHOS Test Suite run on `date` > $LOGFILE
	37	echo "=======================================================================" >> $LOGFILE
	38	echo AliROOT version `which aliroot` >> $LOGFILE
	39	echo Root version `root-config --version` in `which root` >> $LOGFILE
	40	echo "=======================================================================" >> $LOGFILE
	41	#
	42	# SIMULATION
	43	#
	44	printf "Signal simulation ... "
	45	echo " ** Signal SIMULATION ** " >> $LOGFILE
	46	mkdir -p $TMPDIR/signal
	47	cd $TMPDIR/signal
	48	aliroot -b -q ${GRUN}\(100,\"${ConfigSgnl}\"\) >>& $LOGFILE
	49	printf "done\n"
	50
	51	printf "Background simulation ... "
	52	echo " ** Background SIMULATION ** " >> $LOGFILE
	53	mkdir -p $TMPDIR/background
	54	cd $TMPDIR/background
	55	aliroot -b -q ${GRUN}\(2,\"${ConfigBckg}\"\) >>& $LOGFILE
	56	printf "done\n"
	57
	58	cd $TMPDIR
	59	cp -a $TMPDIR/signal $TMPDIR/signal+background
	60
	61	printf "Simulation analysis ... "
	62	cd $TMPDIR/signal
	63	echo '{ ' > tempo.C
	64	echo " Int_t minhits = $MINHITS ;" >> tempo.C
65	echo " Int_t maxhits = $MAXHITS ;" >> tempo.C
66	echo ' AliPHOSGetter * gime = AliPHOSGetter::Instance("galice.root") ;' >> tempo.C
67	echo ' TH1F * hitmul = new TH1F("hitmul", "PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
68	echo ' Int_t max = gime->MaxEvent() ; ' >> tempo.C
69	echo ' Int_t evt = 0 ; ' >> tempo.C
70	echo ' for ( evt = 0 ; evt < max ; evt++ ) { ' >> tempo.C
71	echo ' gime->Event(evt,"H") ; ' >> tempo.C
72	echo ' hitmul->Fill(gime->Hits()->GetEntries()) ; ' >> tempo.C
73	echo ' } ' >> tempo.C
74	echo ' TF1 * gaus = new TF1("gaus", "gaus", 0., 200.) ; ' >> tempo.C
75	echo ' hitmul->Fit(gaus,"", "", 40, 100) ; ' >> tempo.C
76	echo ' if( gaus->GetParameter(1) > maxhits \|\| gaus->GetParameter(1) < minhits ) ' >> tempo.C
77	echo ' printf("ERRORSIM 1") ; ' >> tempo.C
78	echo ' else ' >> tempo.C
79	echo ' printf("ERRORSIM 0") ; ' >> tempo.C
80	echo '} ' >> tempo.C
81	aliroot -b -q tempo.C >>& $LOGFILE
82	set ERRORSIM = `cat $LOGFILE \| grep ERRORSIM \| awk '{print $2}'`
83	if($ERRORSIM) then
84	printf "error\n"
85	uuencode $LOGFILE $LOGFILE \| mail -s 'PHOS INSTALLATION ERROR' $PHOS_LIBRARIAN
86	## exit(0)
87	else
88	printf "done\n"
89	endif
90	rm tempo.C
91	#
92	# RECONSTRUCTION
93	#
94	printf "Reconstruction ... "
95	echo " ** RECONSTRUCTION ** " >> $LOGFILE
96	aliroot -b >>& $LOGFILE <<EOF
97	AliPHOSSDigitizer sd("galice.root") ;
98	sd.ExecuteTask("deb all") ;
99	AliPHOSDigitizer d("galice.root") ;
100	d.ExecuteTask("deb all") ;
101	AliPHOSClusterizerv1 cl("galice.root") ;
102	cl.ExecuteTask("deb all") ;
103	AliPHOSTrackSegmentMakerv1 ts("galice.root") ;
104	ts.ExecuteTask("deb all") ;
105	AliPHOSPIDv1 pd("galice.root") ;
106	pd.ExecuteTask("deb all") ;
107	.q
108	EOF
109	printf "done\n"
110
111	printf "Reconstruction analysis ... "
112	echo '{ ' > tempo.C
113	echo " Int_t minsdig = $MINSDIGITS ;" >> tempo.C
114	echo " Int_t maxsdig = $MAXSDIGITS ;" >> tempo.C
115	echo " Int_t mindig = $MINDIGITS ;" >> tempo.C
116	echo " Int_t maxdig = $MAXDIGITS ;" >> tempo.C
117	echo " Int_t minemcrp = $MINEMCRECPOINTS ;" >> tempo.C
118	echo " Int_t maxemcrp = $MAXEMCRECPOINTS ;" >> tempo.C
119	echo " Int_t mincpvrp = $MINCPVRECPOINTS ;" >> tempo.C
120	echo " Int_t maxcpvrp = $MAXCPVRECPOINTS ;" >> tempo.C
121	echo " Int_t mints = $MINTRACKSEGMENTS ;" >> tempo.C
122	echo " Int_t maxts = $MAXTRACKSEGMENTS ;" >> tempo.C
123	echo " Int_t minpa = $MINRECPARTICLES ;" >> tempo.C
124	echo " Int_t maxpa = $MAXRECPARTICLES ;" >> tempo.C
125	echo ' AliPHOSGetter * gime = AliPHOSGetter::Instance("galice.root") ;' >> tempo.C
126	echo ' TH1F * sdigmul = new TH1F("sdigmul", " SDigits PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
127	echo ' TH1F * digmul = new TH1F("digmul", " Digits PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
128	echo ' TH1F * emcrpmul= new TH1F("emcrpmul", " EMCRecPoints PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
129	echo ' TH1F * cpvrpmul= new TH1F("cpvrpmul", " CPVRecPoints PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
130	echo ' TH1F * tsmul = new TH1F("tsmul", " TrackSegments PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
131	echo ' TH1F * pamul = new TH1F("pamul", " RecParticles PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
132	echo ' Int_t max = gime->MaxEvent() ; ' >> tempo.C
133	echo ' Int_t evt = 0 ; ' >> tempo.C
134	echo ' for ( evt = 0 ; evt < max ; evt++ ) { ' >> tempo.C
135	echo ' gime->Event(evt,"SDRTP") ; ' >> tempo.C
136	echo ' sdigmul->Fill(gime->SDigits()->GetEntries()) ; ' >> tempo.C
137	echo ' digmul->Fill(gime->Digits()->GetEntries()) ; ' >> tempo.C
138	echo ' emcrpmul->Fill(gime->EmcRecPoints()->GetEntries()) ; ' >> tempo.C
139	echo ' cpvrpmul->Fill(gime->CpvRecPoints()->GetEntries()) ; ' >> tempo.C
140	echo ' tsmul->Fill(gime->TrackSegments()->GetEntries()) ; ' >> tempo.C
141	echo ' pamul->Fill(gime->RecParticles()->GetEntries()) ; ' >> tempo.C
142	echo ' } ' >> tempo.C
143	echo ' TF1 * gaus = new TF1("gaus", "gaus", 0., 200.) ; ' >> tempo.C
144	echo ' sdigmul->Fit(gaus, "", "", 40., 100.) ; ' >> tempo.C
145	echo ' sdigmul->Draw() ; ' >> tempo.C
146	echo ' if( gaus->GetParameter(1) > maxsdig \|\| gaus->GetParameter(1) < minsdig ) ' >> tempo.C
147	echo ' printf("ERRORREC 1 sdigits\n") ; ' >> tempo.C
148	echo ' digmul->Fit(gaus, "", "", 20., 80.) ; ' >> tempo.C
149	echo ' if( gaus->GetParameter(1) > maxdig \|\| gaus->GetParameter(1) < mindig ) ' >> tempo.C
150	echo ' printf("ERRORREC 1 digits\n") ; ' >> tempo.C
151	echo ' emcrpmul->Fit(gaus, "", "", 0., 4.) ; ' >> tempo.C
152	echo ' if( gaus->GetParameter(1) > maxemcrp \|\| gaus->GetParameter(1) < minemcrp ) ' >> tempo.C
153	echo ' printf("ERRORREC 1 emc recpoints\n") ; ' >> tempo.C
154	echo ' cpvrpmul->Fit(gaus, "", "", 0., 4.) ; ' >> tempo.C
155	echo ' if( gaus->GetParameter(1) > maxcpvrp \|\| gaus->GetParameter(1) < mincpvrp ) ' >> tempo.C
156	echo ' printf("ERRORREC 1 cpv recpoints\n") ; ' >> tempo.C
157	echo ' tsmul->Fit(gaus, "", "", 0., 4.) ; ' >> tempo.C
158	echo ' if( gaus->GetParameter(1) > maxts \|\| gaus->GetParameter(1) < mints ) ' >> tempo.C
159	echo ' printf("ERRORREC 1 track segments\n") ; ' >> tempo.C
160	echo ' pamul->Fit(gaus, "", "", 0., 4. ) ; ' >> tempo.C
161	echo ' if( gaus->GetParameter(1) > maxpa \|\| gaus->GetParameter(1) < minpa ) ' >> tempo.C
162	echo ' printf("ERRORREC 1 recparticles\n") ; ' >> tempo.C
163	echo '} ' >> tempo.C
164	aliroot -b -q tempo.C >>& $LOGFILE
165	set ERRORREC = `cat $LOGFILE \| grep ERRORREC \| awk '{print $2}'`
166	if( $?ERRORREC ) then
167	printf "error\n"
168	uuencode $LOGFILE $LOGFILE \| mail -s 'PHOS INSTALLATION ERROR' $PHOS_LIBRARIAN
169	## exit(0)
170	else
171	printf "done\n"
172	endif
173	rm tempo.C
174
175	# Signal+background merging
176	printf "Signal+background merging ... "
177	echo " ** Signal+background MERGING ** " >> $LOGFILE
178	cd $TMPDIR/background
179	aliroot -b >>& $LOGFILE <<EOF
180	AliPHOSSDigitizer sd("galice.root") ;
181	sd.ExecuteTask("deb all") ;
182	EOF
183
184	cd $TMPDIR/signal+background
185	cp -a $TMPDIR/signal/PHOS.SDigits*.root .
186	aliroot -b >>& $LOGFILE <<EOF
187	AliPHOSDigitizer d("galice.root") ;
188	d.MixWith("../background/galice.root") ;
189	d.ExecuteTask("deb all") ;
190	AliPHOSClusterizerv1 cl("galice.root") ;
191	cl.ExecuteTask("deb all") ;
192	AliPHOSTrackSegmentMakerv1 ts("galice.root") ;
193	ts.ExecuteTask("deb all") ;
194	AliPHOSPIDv1 pd("galice.root") ;
195	pd.ExecuteTask("deb all") ;
196	.q
197	EOF
198	printf "done\n"
199
200	#rm -fr $TMPDIR