ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/VHbb/python/evaluateMVA.py
(Generate patch)

Comparing UserCode/VHbb/python/evaluateMVA.py (file contents):
Revision 1.10 by nmohr, Thu Aug 9 13:52:04 2012 UTC vs.
Revision 1.14 by peller, Thu Sep 27 07:34:24 2012 UTC

# Line 31 | Line 31 | parser.add_option("-P", "--path", dest="
31                        help="path to samples")
32   parser.add_option("-S", "--samples", dest="names", default="",
33                        help="samples you want to run on")
34 < parser.add_option("-C", "--config", dest="config", default="",
34 > parser.add_option("-C", "--config", dest="config", default=[], action="append",
35                        help="configuration file")
36   (opts, args) = parser.parse_args(argv)
37   if opts.config =="":
38          opts.config = "config"
39   config = BetterConfigParser()
40 + #config.read('./config7TeV_ZZ')
41   config.read(opts.config)
42   anaTag = config.get("Analysis","tag")
43  
44   #get locations:
45   Wdir=config.get('Directories','Wdir')
46 <
46 < MVAdir=config.get('Directories','MVAdir')
46 > MVASubdir=config.get('Directories','MVAdir')
47  
48   #systematics
49   systematics=config.get('systematics','systematics')
# Line 70 | Line 70 | namelist=namelistIN.split(',')
70   doinfo=bool(int(opts.update))
71  
72   MVAlist=arglist.split(',')
73 + MVAdir=config.get('Directories','vhbbpath')
74  
75   #CONFIG
76   #factory
# Line 83 | Line 84 | factoryname=config.get('factory','factor
84   #MVAinfofiles=[]
85   MVAinfos=[]
86   for MVAname in MVAlist:
87 <    MVAinfofile = open(Wdir+'/weights/'+factoryname+'_'+MVAname+'.info','r')
87 >    MVAinfofile = open(MVAdir+'/data/'+factoryname+'_'+MVAname+'.info','r')
88      MVAinfos.append(pickle.load(MVAinfofile))
89      MVAinfofile.close()
90      
# Line 122 | Line 123 | for i in range(len( MVA_Vars['Nominal'])
123   #        reader.AddSpectator(spectators[i],MVA_spectator_buffer[i])
124   #Load raeder
125   for i in range(0,len(readers)):
126 <    readers[i].BookMVA(MVAinfos[i].MVAname,MVAinfos[i].getweightfile())
126 >    readers[i].BookMVA(MVAinfos[i].MVAname,MVAdir+'/data/'+MVAinfos[i].getweightfile())
127   #--> Now the MVA is booked
128  
129   #Apply samples
# Line 135 | Line 136 | for job in Ainfo:
136      if eval(job.active):
137          if job.name in namelist:
138              #get trees:
139 <            input = TFile.Open(job.getpath(),'read')
140 <            outfile = TFile.Open(job.path+'/'+MVAdir+job.prefix+job.identifier+'.root','recreate')
139 >            input = TFile.Open(Apath+'/'+job.getpath(),'read')
140 >            outfile = TFile.Open(Apath+'/'+MVASubdir+job.prefix+job.identifier+'.root','recreate')
141              input.cd()
142              obj = ROOT.TObject
143              for key in ROOT.gDirectory.GetListOfKeys():

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines