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

#! /bin/tcsh -f
#
# Prepare the environment
#
mkdir $HOME/PHOSTestSuite ; cd $HOME/PHOSTestSuite
cp $ALICE_ROOT/PHOS/grunTestSuite.C . 
cp $ALICE_ROOT/PHOS/ConfigTestSuite.C . 
#
#begin data
set MINHITS=68
set MAXHITS=72
set MINSDIGITS=68
set MAXSDIGITS=72
set MINDIGITS=40
set MAXDIGITS=48
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` > TestSuite.log
echo "=======================================================================" >> TestSuite.log
echo AliROOT  version `which aliroot` >> TestSuite.log
echo Root     version `which root` >> TestSuite.log
echo "=======================================================================" >> TestSuite.log
#
# SIMULATION
#
echo "    **** SIMULATION **** " >> TestSuite.log
aliroot -b -q grunTestSuite.C\(100\) >>& TestSuite.log 
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 >>& TestSuite.log 
set ERRORSIM = `cat TestSuite.log | grep ERRORSIM | awk '{print $2}'`
if($ERRORSIM) then
  echo ERROR
  uuencode  TestSuite.log  TestSuite.log | mail -s 'PHOS INSTALLATION ERROR' yves.schutz@cern.ch
  exit(0)
endif 
rm tempo.C 
#
# RECONSTRUCTION
#
echo "    **** RECONSTRUCTION **** " >> TestSuite.log
aliroot -b >>& TestSuite.log <<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
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  >>& TestSuite.log 
set ERRORREC = `cat TestSuite.log | grep ERRORREC | awk '{print $2}'`
if($ERRORREC) then
  echo ERROR
  uuencode  TestSuite.log  TestSuite.log | mail -s 'PHOS INSTALLATION ERROR' yves.schutz@cern.ch
  exit(0)
endif 
#rm tempo.C 
#/cd $HOME
#rm -fr $HOME/PHOSTestSuite
Commit	Line	Data
88cb7938	1	#! /bin/tcsh -f
	2	#
	3	# Prepare the environment
	4	#
	5	mkdir $HOME/PHOSTestSuite ; cd $HOME/PHOSTestSuite
	6	cp $ALICE_ROOT/PHOS/grunTestSuite.C .
	7	cp $ALICE_ROOT/PHOS/ConfigTestSuite.C .
	8	#
	9	#begin data
	10	set MINHITS=68
	11	set MAXHITS=72
	12	set MINSDIGITS=68
	13	set MAXSDIGITS=72
b4cb3c78	14	set MINDIGITS=40
b4cb3c78	15	set MAXDIGITS=48
88cb7938	16	set MINEMCRECPOINTS=1
	17	set MAXEMCRECPOINTS=2
	18	set MINCPVRECPOINTS=0
	19	set MAXCPVRECPOINTS=2
	20	set MINTRACKSEGMENTS=1
	21	set MAXTRACKSEGMENTS=2
	22	set MINRECPARTICLES=1
	23	set MAXRECPARTICLES=2
	24	#end data
	25	#
	26	echo PHOS Test Suite run on `date` > TestSuite.log
	27	echo "=======================================================================" >> TestSuite.log
	28	echo AliROOT version `which aliroot` >> TestSuite.log
	29	echo Root version `which root` >> TestSuite.log
	30	echo "=======================================================================" >> TestSuite.log
	31	#
	32	# SIMULATION
	33	#
	34	echo " ** SIMULATION ** " >> TestSuite.log
	35	aliroot -b -q grunTestSuite.C\(100\) >>& TestSuite.log
	36	echo '{ ' > tempo.C
	37	echo " Int_t minhits = $MINHITS ;" >> tempo.C
	38	echo " Int_t maxhits = $MAXHITS ;" >> tempo.C
	39	echo ' AliPHOSGetter * gime = AliPHOSGetter::Instance("galice.root") ;' >> tempo.C
	40	echo ' TH1F * hitmul = new TH1F("hitmul", "PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
	41	echo ' Int_t max = gime->MaxEvent() ; ' >> tempo.C
	42	echo ' Int_t evt = 0 ; ' >> tempo.C
	43	echo ' for ( evt = 0 ; evt < max ; evt++ ) { ' >> tempo.C
	44	echo ' gime->Event(evt,"H") ; ' >> tempo.C
	45	echo ' hitmul->Fill(gime->Hits()->GetEntries()) ; ' >> tempo.C
	46	echo ' } ' >> tempo.C
	47	echo ' TF1 * gaus = new TF1("gaus", "gaus", 0., 200.) ; ' >> tempo.C
	48	echo ' hitmul->Fit(gaus,"", "", 40, 100) ; ' >> tempo.C
	49	echo ' if( gaus->GetParameter(1) > maxhits \|\| gaus->GetParameter(1) < minhits ) ' >> tempo.C
	50	echo ' printf("ERRORSIM 1") ; ' >> tempo.C
	51	echo ' else ' >> tempo.C
	52	echo ' printf("ERRORSIM 0") ; ' >> tempo.C
	53	echo '} ' >> tempo.C
	54	aliroot -b -q tempo.C >>& TestSuite.log
	55	set ERRORSIM = `cat TestSuite.log \| grep ERRORSIM \| awk '{print $2}'`
	56	if($ERRORSIM) then
	57	echo ERROR
	58	uuencode TestSuite.log TestSuite.log \| mail -s 'PHOS INSTALLATION ERROR' yves.schutz@cern.ch
	59	exit(0)
	60	endif
	61	rm tempo.C
	62	#
	63	# RECONSTRUCTION
	64	#
	65	echo " ** RECONSTRUCTION ** " >> TestSuite.log
	66	aliroot -b >>& TestSuite.log <<EOF
	67	AliPHOSSDigitizer sd("galice.root") ;
	68	sd.ExecuteTask("deb all") ;
	69	AliPHOSDigitizer d("galice.root") ;
	70	d.ExecuteTask("deb all") ;
	71	AliPHOSClusterizerv1 cl("galice.root") ;
	72	cl.ExecuteTask("deb all") ;
	73	AliPHOSTrackSegmentMakerv1 ts("galice.root") ;
	74	ts.ExecuteTask("deb all") ;
	75	AliPHOSPIDv1 pd("galice.root") ;
	76	pd.ExecuteTask("deb all") ;
	77	.q
	78	EOF
	79	echo '{ ' > tempo.C
80	echo " Int_t minsdig = $MINSDIGITS ;" >> tempo.C
81	echo " Int_t maxsdig = $MAXSDIGITS ;" >> tempo.C
82	echo " Int_t mindig = $MINDIGITS ;" >> tempo.C
83	echo " Int_t maxdig = $MAXDIGITS ;" >> tempo.C
84	echo " Int_t minemcrp = $MINEMCRECPOINTS ;" >> tempo.C
85	echo " Int_t maxemcrp = $MAXEMCRECPOINTS ;" >> tempo.C
86	echo " Int_t mincpvrp = $MINCPVRECPOINTS ;" >> tempo.C
87	echo " Int_t maxcpvrp = $MAXCPVRECPOINTS ;" >> tempo.C
88	echo " Int_t mints = $MINTRACKSEGMENTS ;" >> tempo.C
89	echo " Int_t maxts = $MAXTRACKSEGMENTS ;" >> tempo.C
90	echo " Int_t minpa = $MINRECPARTICLES ;" >> tempo.C
91	echo " Int_t maxpa = $MAXRECPARTICLES ;" >> tempo.C
92	echo ' AliPHOSGetter * gime = AliPHOSGetter::Instance("galice.root") ;' >> tempo.C
93	echo ' TH1F * sdigmul = new TH1F("sdigmul", " SDigits PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
94	echo ' TH1F * digmul = new TH1F("digmul", " Digits PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
95	echo ' TH1F * emcrpmul= new TH1F("emcrpmul", " EMCRecPoints PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
96	echo ' TH1F * cpvrpmul= new TH1F("cpvrpmul", " CPVRecPoints PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
97	echo ' TH1F * tsmul = new TH1F("tsmul", " TrackSegments PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
98	echo ' TH1F * pamul = new TH1F("pamul", " RecParticles PHOS Test Suite", 100, 0., 200.) ; ' >> tempo.C
99	echo ' Int_t max = gime->MaxEvent() ; ' >> tempo.C
100	echo ' Int_t evt = 0 ; ' >> tempo.C
101	echo ' for ( evt = 0 ; evt < max ; evt++ ) { ' >> tempo.C
102	echo ' gime->Event(evt,"SDRTP") ; ' >> tempo.C
103	echo ' sdigmul->Fill(gime->SDigits()->GetEntries()) ; ' >> tempo.C
104	echo ' digmul->Fill(gime->Digits()->GetEntries()) ; ' >> tempo.C
105	echo ' emcrpmul->Fill(gime->EmcRecPoints()->GetEntries()) ; ' >> tempo.C
106	echo ' cpvrpmul->Fill(gime->CpvRecPoints()->GetEntries()) ; ' >> tempo.C
107	echo ' tsmul->Fill(gime->TrackSegments()->GetEntries()) ; ' >> tempo.C
108	echo ' pamul->Fill(gime->RecParticles()->GetEntries()) ; ' >> tempo.C
109	echo ' } ' >> tempo.C
110	echo ' TF1 * gaus = new TF1("gaus", "gaus", 0., 200.) ; ' >> tempo.C
111	echo ' sdigmul->Fit(gaus, "", "", 40., 100.) ; ' >> tempo.C
112	echo ' sdigmul->Draw() ; ' >> tempo.C
113	echo ' if( gaus->GetParameter(1) > maxsdig \|\| gaus->GetParameter(1) < minsdig ) ' >> tempo.C
114	echo ' printf("ERRORREC 1 sdigits\n") ; ' >> tempo.C
115	echo ' digmul->Fit(gaus, "", "", 20., 80.) ; ' >> tempo.C
116	echo ' if( gaus->GetParameter(1) > maxdig \|\| gaus->GetParameter(1) < mindig ) ' >> tempo.C
117	echo ' printf("ERRORREC 1 digits\n") ; ' >> tempo.C
118	echo ' emcrpmul->Fit(gaus, "", "", 0., 4.) ; ' >> tempo.C
119	echo ' if( gaus->GetParameter(1) > maxemcrp \|\| gaus->GetParameter(1) < minemcrp ) ' >> tempo.C
120	echo ' printf("ERRORREC 1 emc recpoints\n") ; ' >> tempo.C
121	echo ' cpvrpmul->Fit(gaus, "", "", 0., 4.) ; ' >> tempo.C
122	echo ' if( gaus->GetParameter(1) > maxcpvrp \|\| gaus->GetParameter(1) < mincpvrp ) ' >> tempo.C
123	echo ' printf("ERRORREC 1 cpv recpoints\n") ; ' >> tempo.C
124	echo ' tsmul->Fit(gaus, "", "", 0., 4.) ; ' >> tempo.C
125	echo ' if( gaus->GetParameter(1) > maxts \|\| gaus->GetParameter(1) < mints ) ' >> tempo.C
126	echo ' printf("ERRORREC 1 track segments\n") ; ' >> tempo.C
127	echo ' pamul->Fit(gaus, "", "", 0., 4. ) ; ' >> tempo.C
128	echo ' if( gaus->GetParameter(1) > maxpa \|\| gaus->GetParameter(1) < minpa ) ' >> tempo.C
129	echo ' printf("ERRORREC 1 recparticles\n") ; ' >> tempo.C
130	echo '} ' >> tempo.C
131	aliroot -b -q tempo.C >>& TestSuite.log
132	set ERRORREC = `cat TestSuite.log \| grep ERRORREC \| awk '{print $2}'`
133	if($ERRORREC) then
134	echo ERROR
135	uuencode TestSuite.log TestSuite.log \| mail -s 'PHOS INSTALLATION ERROR' yves.schutz@cern.ch
136	exit(0)
137	endif
138	#rm tempo.C
139	#/cd $HOME
140	#rm -fr $HOME/PHOSTestSuite