Main Page | Modules | Namespace List | Class Hierarchy | Alphabetical List | Class List | Directories | File List | Namespace Members | Class Members | File Members | Related Pages

Anp::PlotShort Class Reference

#include <PlotShort.h>

Inheritance diagram for Anp::PlotShort:

Anp::AlgSnarl Anp::Base List of all members.

Public Member Functions

 PlotShort ()
virtual ~PlotShort ()
bool Run (Record &record)
void Set (TDirectory *dir)
void Config (const Registry &reg)
bool Init (const Header &header)
void End (const DataBlock &)

Private Types

typedef std::map< std::string,
TH1 * > 		MapTH1
typedef std::map< std::string,
TH2 * > 		MapTH2
typedef std::map< std::string,
TTree * > 		MapTree
typedef std::map< int, float > BMap

Private Member Functions

void Plot (const Event &event, const Record &record)
void Plot (const Track &track, const Record &record)
void PlotStandard (const Event &event, const Record &record)
void PlotStandard (const Track &track, const Record &record)
void FillStandard (TDirectory *ptype, const Event &event, const Record &record)
void FillStandard (TDirectory *ptype, const Track &track, const Record &record)
TH1 * GetTH1 (const std::string &key, const std::string &name="")
TH2 * GetTH2 (const std::string &name="", TDirectory *dir=0, const std::string &plottype="unknown")
Particle::Particle_t Particle (const Track &track, const Record &record, int booltype) const

Private Attributes

std::string fDirName
TDirectory * fDir
TFile * fFile
TTree * fTree
std::string fPlotKey
bool fTrack
bool fEvent
bool fPlotTree
std::string fTreeName
bool fPlot
bool fPlotStandard
bool fSVMPrint
int fPassKey
bool fDebug
int fNMiss
bool fIsData
std::vector< int > fKeys
std::vector< int > fEKeys
std::vector< int > fTKeys
std::vector< int > fKNNKeys
std::map< std::string, TH1 * > fMapTH1
std::map< std::string, TH2 * > fMapTH2
std::map< std::string, TH2 * > fMapMu
std::map< std::string, TH2 * > fMapPi
std::map< std::string, TH2 * > fMapPr
std::map< std::string, TH2 * > fMapNC
std::map< std::string, TH2 * > fMapCC
std::map< std::string, TTree * > fMapTree
std::map< int, float > fTreeValues
float energy
int ptype
int ntrack
int nplane
int ntrackplane
float comp
float pur
float tpur
float tcompall
float tcompslc
int ttype
int nevent

Member Typedef Documentation

typedef std::map<int,float> Anp::PlotShort::BMap [private]
 

Definition at line 59 of file PlotShort.h.

typedef std::map<std::string, TH1 *> Anp::PlotShort::MapTH1 [private]
 

Definition at line 56 of file PlotShort.h.

typedef std::map<std::string, TH2 *> Anp::PlotShort::MapTH2 [private]
 

Definition at line 57 of file PlotShort.h.

Referenced by GetTH2().

typedef std::map<std::string, TTree *> Anp::PlotShort::MapTree [private]
 

Definition at line 58 of file PlotShort.h.

Constructor & Destructor Documentation

Anp::PlotShort::PlotShort (   ) 
 

Definition at line 31 of file PlotShort.cxx.

00032   :fDirName("2D_histos"),
00033    fDir(0),
00034    fPlotKey("knnvars"),
00035    fTrack(true),
00036    fEvent(true),
00037    fPlotTree(false),
00038    fTreeName("treefile.root"),
00039    fPlot(true),
00040    fPlotStandard(true),
00041    fSVMPrint(false),
00042    fPassKey(-1),
00043    fDebug(true),
00044    fNMiss(0)
00045 {
00046 }

Anp::PlotShort::~PlotShort (   )  [virtual]
 

Definition at line 49 of file PlotShort.cxx.

00050 {
00051 }

Member Function Documentation

void Anp::PlotShort::Config (  const Registry &  reg  )  [virtual]
 

Reimplemented from Anp::AlgSnarl.

Definition at line 144 of file PlotShort.cxx.

References fDebug, fDirName, fEKeys, fKeys, fKNNKeys, fPassKey, fPlotKey, fPlotStandard, fPlotTree, fSVMPrint, fTKeys, fTrack, fTreeName, Registry::Get(), Registry::KeyExists(), Anp::Read(), and reg.

00145 {
00146    //
00147    // Read string configurations
00148    //
00149    Anp::Read(reg, "PlotShortDirName",  fDirName);   
00150    Anp::Read(reg, "PlotShortTreeName", fTreeName);
00151    Anp::Read(reg, "PlotShortPlotKey",  fPlotKey);
00152    Anp::Read(reg, "PlotShortTrack",    fTrack);
00153    Anp::Read(reg, "PlotShortTree",     fPlotTree);
00154    Anp::Read(reg, "PlotShortStandard", fPlotStandard);
00155    Anp::Read(reg, "PlotShortSVM",fSVMPrint);
00156    Anp::Read(reg, "PlotShortDebug",    fDebug);
00157    reg.Get("PlotShortPassKey",fPassKey);
00158    //
00159    // Get and sortkeys to be plotted
00160    //
00161 
00162    fKeys    = Anp::ReadList<int>(reg, "PlotShortList"     , " ");
00163    fEKeys   = Anp::ReadList<int>(reg, "PlotShortListEvent", " ");
00164    fTKeys   = Anp::ReadList<int>(reg, "PlotShortListTrack", " ");
00165    fKNNKeys = Anp::ReadList<int>(reg, "PlotShortListkNN"  , " ");
00166    std::sort(fEKeys.begin(), fEKeys.end());
00167    std::sort(fTKeys.begin(), fTKeys.end());
00168    std::sort(fKNNKeys.begin(), fKNNKeys.end());
00169    cout<<" Plot Short Configuration "<<endl;
00170    cout<<" kNN keys "<<endl;
00171    for(unsigned int ii=0; ii< fKNNKeys.size(); ii++) 
00172      cout<<" KNN  key :  "<<fKNNKeys[ii]<<endl; 
00173    
00174    if(reg.KeyExists("PrintConfig"))
00175    {
00176       cout << "PlotShort::Config" << endl
00177            << "   DirName = " << fDirName << endl
00178            << "   PassKey = " << fPassKey << endl;
00179    }
00180 }

void Anp::PlotShort::End (  const DataBlock &   )  [virtual]
 

Reimplemented from Anp::AlgSnarl.

Definition at line 187 of file PlotShort.cxx.

00188 {
00189 
00190 }

void Anp::PlotShort::FillStandard (  TDirectory *  ptype,
const Track &  track,
const Record &  record
)  [private]
 

Definition at line 318 of file PlotShort.cxx.

References Anp::Track::DataAt(), Anp::Track::ErrorQP(), Anp::Record::EventEnd(), Anp::Record::FindEvent(), Anp::Record::FindTruth(), fKNNKeys, Anp::Track::GetBasic(), Anp::Track::GetBegVtx(), Anp::Event::GetRecoNu(), Anp::Track::KeyExists(), Anp::Basic::MIP(), Anp::Basic::NUPlane(), Anp::Track::QP(), Anp::Vertex::R(), Anp::Track::TrackIndex(), Anp::Record::TruthEnd(), Anp::TruthIter, Anp::TruthRecoIter, Anp::Vertex::U(), Anp::Vertex::V(), Anp::RecoNu::Y(), and Anp::Vertex::Z().

00318                                                                                           {
00319 
00320   std::string totname =ptype->GetName(); // Get the name of the store directory (particle type)
00321   //  std:: string totname = name+na;
00322 
00323   double purit =0;
00324   double complete =0;
00325   double true_enu =-1;
00326   if(!fIsData){
00327     const TruthIter truth = record.FindTruth(track);
00328     if(!(truth==record.TruthEnd())){
00329       true_enu = truth->ENeu();
00330       for(TruthRecoIter ireco=truth->RecoBeg(); ireco!=truth->RecoEnd(); ireco++)
00331         {
00332           if(ireco -> IsTrack() &&ireco->RecoIndex() == track.TrackIndex())
00333             {
00334               purit= ireco->Purity();
00335               complete= ireco->ComplAll();
00336             }
00337         }
00338     }
00339   }
00340   double mip = track.GetBasic().MIP();
00341   double ang =-1000;
00342   if(track.KeyExists(108)) ang = track.DataAt(108);
00343   for(unsigned int ii=0; ii< fKNNKeys.size(); ii++){
00344     stringstream knn;
00345     knn<<fKNNKeys[ii];
00346 
00347     double fknnval =-0.1;
00348     if(track.KeyExists(fKNNKeys[ii]) )fknnval= track.DataAt(fKNNKeys[ii]);
00349 
00350     TH2 *energy_v_kNN= PlotShort::GetTH2("mip_knn"+knn.str(), ptype,totname);
00351     if(energy_v_kNN)  energy_v_kNN->Fill(mip, fknnval);
00352     TH2 *qpe_v_kNN = PlotShort::GetTH2("qpe_knn"+knn.str(), ptype,totname);
00353     if(qpe_v_kNN) qpe_v_kNN->Fill(track.QP()/track.ErrorQP(), fknnval);
00354     TH2 *ang_v_kNN = PlotShort::GetTH2("ang_knn"+knn.str(), ptype,totname);
00355     if(ang_v_kNN&&track.KeyExists(108)) ang_v_kNN->Fill(ang,fknnval);
00356     
00357     if(!fIsData){
00358       TH2 *pur_v_kNN = PlotShort::GetTH2("pur_knn"+knn.str(), ptype,totname);
00359       if(pur_v_kNN)    pur_v_kNN->Fill(purit, fknnval);
00360       TH2 *comp_v_kNN = PlotShort::GetTH2("comp_knn"+knn.str(), ptype,totname);
00361       if(comp_v_kNN)      comp_v_kNN->Fill(complete, fknnval);
00362       TH2 *truE_v_kNN = PlotShort::GetTH2("true_knn"+knn.str(), ptype,totname);
00363       if(truE_v_kNN)      truE_v_kNN->Fill(true_enu, fknnval);
00364     }
00365   }
00366 
00367   TH2 *d835_v_d836 = PlotShort::GetTH2("d835_d836", ptype,totname);
00368   if(d835_v_d836 && track.KeyExists(14835) && track.KeyExists(14836) ) d835_v_d836->Fill(track.DataAt(14835), track.DataAt(14836));
00369   TH2 *d835_v_d001 = PlotShort::GetTH2("d835_d001", ptype,totname);
00370   if(d835_v_d001 && track.KeyExists(14835) && track.KeyExists(14001) ) d835_v_d001->Fill(track.DataAt(14835), track.DataAt(14001));
00371   TH2 *d835_v_d805 = PlotShort::GetTH2("d835_d805", ptype,totname);
00372   if(d835_v_d805 && track.KeyExists(14835) && track.KeyExists(14805) ) d835_v_d805->Fill(track.DataAt(14835), track.DataAt(14805));
00373   TH2 *d420_v_d805 = PlotShort::GetTH2("d420_d805", ptype,totname);
00374   if(d420_v_d805 && track.KeyExists(14420) && track.KeyExists(14805) ) d420_v_d805->Fill(track.DataAt(14420), track.DataAt(14805));
00375   TH2 *d420_v_d835 = PlotShort::GetTH2("d420_d835", ptype,totname);
00376   if(d420_v_d835 && track.KeyExists(14835) && track.KeyExists(14420) ) d420_v_d835->Fill(track.DataAt(14420), track.DataAt(14835));
00377   
00378   if(record.FindEvent(track)!=record.EventEnd())
00379   {
00380     Event evt=* ((EventIter)record.FindEvent(track));
00381     const RecoNu& rnu =( const RecoNu&) evt.GetRecoNu();
00382     double y= rnu.Y();
00383     TH2 *y_v_d835 = PlotShort::GetTH2("y_d835", ptype,totname);
00384     if(y_v_d835 && y>=0&& track.KeyExists(14420) ) y_v_d835->Fill(y, track.DataAt(14835));
00385   }
00386   
00387   double r = track.GetBegVtx().R();
00388   double z = track.GetBegVtx().Z();
00389   double u = track.GetBegVtx().U();
00390   double v = track.GetBegVtx().V();
00391   double ntrku = track.GetBasic().NUPlane();
00392 
00393   TH2 *r_v_d835 = PlotShort::GetTH2("r_d835", ptype,totname);
00394   if(r_v_d835 && r>0 && track.KeyExists(14835) ) r_v_d835->Fill(r, track.DataAt(14835));
00395   TH2 *z_v_d835 = PlotShort::GetTH2("z_d835", ptype,totname);
00396   if(z_v_d835 && z>=0 && track.KeyExists(14835) ) z_v_d835->Fill(z, track.DataAt(14835));  
00397   TH2 *u_v_d835 = PlotShort::GetTH2("u_d835", ptype,totname);
00398   if(u_v_d835 && u>=-10000 && track.KeyExists(14835) ) u_v_d835->Fill(u, track.DataAt(14835));
00399   TH2 *v_v_d835 = PlotShort::GetTH2("v_d835", ptype,totname);
00400   if(v_v_d835 && v>=-10000 && track.KeyExists(14835) ) v_v_d835->Fill(v, track.DataAt(14835));
00401   TH2 *trku_v_d835 = PlotShort::GetTH2("trku_d835", ptype,totname);
00402   if(trku_v_d835 && ntrku>0 && track.KeyExists(14835) ) trku_v_d835->Fill(ntrku, track.DataAt(14835));
00403   TH2 *d833_v_d834 = PlotShort::GetTH2("d833_d834", ptype,totname);
00404   if(d833_v_d834 && track.KeyExists(14833) &&track.KeyExists(14834)  ) d833_v_d834->Fill(track.DataAt(14833), track.DataAt(14834));
00405 
00406 }

void Anp::PlotShort::FillStandard (  TDirectory *  ptype,
const Event &  event,
const Record &  record
)  [private]
 

Definition at line 246 of file PlotShort.cxx.

References Anp::Event::DataAt(), energy, Anp::Record::FindTruth(), fKNNKeys, Anp::Event::Gev(), Anp::Event::KeyExists(), Anp::LongestTrack(), Anp::Record::TrackEnd(), Anp::TrackIter, Anp::Record::TruthEnd(), and Anp::TruthIter.

00246                                                                                           {
00247 
00248   std::string totname =ptype->GetName(); // Get the name of the store directory (particle type)    
00249   double true_enu = -1;
00250   double true_y =-1;
00251   if(!fIsData){
00252     const TruthIter truth = record.FindTruth(event);
00253     if(!(truth==record.TruthEnd())){
00254       true_enu = truth->ENeu();
00255       true_y   = truth->TrueY();
00256     }
00257   }
00258   double reco_y= 0;
00259   TH2 *a1_a2 =PlotShort::GetTH2("a1_a2",ptype,totname);
00260   TH2 *b1_b2 =PlotShort::GetTH2("b1_b2",ptype,totname);
00261   TH2 *c1_c2 =PlotShort::GetTH2("c1_c2",ptype,totname);
00262   TH2 *d1_d2 =PlotShort::GetTH2("d1_d2",ptype,totname);
00263   TH2 *e1_e2 =PlotShort::GetTH2("e1_e2",ptype,totname);
00264   TH2 *f1_f2 =PlotShort::GetTH2("f1_f2",ptype,totname);
00265 
00266   if( a1_a2 && event.KeyExists(19004) && event.KeyExists(19104))  a1_a2->Fill(event.DataAt(19004), event.DataAt(19104));
00267   if( b1_b2 && event.KeyExists(19005) && event.KeyExists(19105))  b1_b2->Fill(event.DataAt(19005), event.DataAt(19105));
00268   if( c1_c2 && event.KeyExists(19008) && event.KeyExists(19108))  c1_c2->Fill(event.DataAt(19008), event.DataAt(19108));
00269   if( d1_d2 && event.KeyExists(19208) && event.KeyExists(19204))  d1_d2->Fill(event.DataAt(19208), event.DataAt(19204));
00270   if( e1_e2 && event.KeyExists(19206) && event.KeyExists(19207))  e1_e2->Fill(event.DataAt(19206), event.DataAt(19207));
00271   if( f1_f2 && event.KeyExists(19203) && event.KeyExists(19201))  f1_f2->Fill(event.DataAt(19203), event.DataAt(19201));
00272 
00273 
00274     TH2 *recoy_truey = PlotShort::GetTH2("recoy_truey", ptype,totname);
00275     if(recoy_truey) recoy_truey->Fill(reco_y,true_y);
00276 
00277     if(!fIsData)
00278     {
00279         TH2 *truE_recoE= PlotShort::GetTH2("truE_recoE", ptype,totname);
00280         if(truE_recoE)      truE_recoE->Fill(true_enu, event.Gev());
00281     }
00282     double energy = event.Gev();
00283     TrackIter trk = Anp::LongestTrack(event,record);
00284 
00285     for(unsigned int ii=0; ii< fKNNKeys.size(); ii++){
00286       double fknnval =-0.1;
00287       if(fKNNKeys[ii]>19000 && event.KeyExists(fKNNKeys[ii])) fknnval = event.DataAt(fKNNKeys[ii]);
00288       else if(!(trk==record.TrackEnd())&&trk->KeyExists(fKNNKeys[ii]) )fknnval= trk->DataAt(fKNNKeys[ii]);
00289       
00290       stringstream knnname;
00291       knnname<<"recoE_knn"<<fKNNKeys[ii];
00292       TH2 *energy_v_kNN= PlotShort::GetTH2(knnname.str(), ptype,totname);
00293       if(energy_v_kNN)  energy_v_kNN->Fill(energy, fknnval);
00294     }
00295     if(fKNNKeys.size()>1 &&!( trk==record.TrackEnd())){
00296       double knn1 = -0.1;
00297       double knn2 = -0.1;
00298       if(fKNNKeys[0]>19000 && event.KeyExists(fKNNKeys[0])) knn1 = event.DataAt(fKNNKeys[0]);
00299       else if(trk->KeyExists(fKNNKeys[0]) )knn1= trk->DataAt(fKNNKeys[0]);
00300       if(fKNNKeys[1]>19000 && event.KeyExists(fKNNKeys[1]))  knn2 = event.DataAt(fKNNKeys[1]);
00301       else if(trk->KeyExists(fKNNKeys[1]) )knn2= trk->DataAt(fKNNKeys[1]);
00302 
00303       double knn1rad =-0.1;
00304       if(trk->KeyExists(4200)) knn1rad= trk->DataAt(4200);
00305       double knn2rad =-0.1;
00306       if(trk->KeyExists(2400)) knn2rad=trk->DataAt(2400);
00307 
00308       TH2 *kNN_v_kNN= PlotShort::GetTH2("knn_knn", ptype,totname);
00309       if(kNN_v_kNN) kNN_v_kNN->Fill(knn1,knn2);
00310       TH2 *kNN4280_v_kNN4200= PlotShort::GetTH2("knn4280_knn4200", ptype,totname);
00311       if(kNN4280_v_kNN4200) kNN4280_v_kNN4200->Fill(knn1,knn1rad);
00312       TH2 *kNN2480_v_kNN2400= PlotShort::GetTH2("knn2480_knn2400", ptype,totname);
00313       if(kNN2480_v_kNN2400) kNN2480_v_kNN2400->Fill(knn2,knn2rad);
00314 
00315     }
00316 }

TH1 * Anp::PlotShort::GetTH1 (  const std::string &  key,
const std::string &  name = ""
)  [private]
 

Definition at line 488 of file PlotShort.cxx.

References fDir, and Anp::SetDir().

00489 {
00490    TH1 *h = HistMan::Instance().CreateTH1(key, "short");
00491    if(h)
00492    {
00493       Anp::SetDir(h, fDir, name);
00494    }
00495    else
00496    {
00497       ++fNMiss;
00498    }
00499    
00500    return h;
00501 }

TH2 * Anp::PlotShort::GetTH2 (  const std::string &  name = "",
TDirectory *  dir = 0,
const std::string &  plottype = "unknown"
)  [private]
 

Definition at line 504 of file PlotShort.cxx.

References MapTH2, and Anp::SetDir().

00505 {
00506   MapTH2::iterator hit ;
00507   MapTH2::iterator end ;
00508 
00509   MapTH2* th2map ;
00510   string mapkey=name;
00511   if(fPlotStandard){
00512     if     (plottype=="mu")  th2map = &fMapMu;
00513     else if(plottype=="pion")th2map= &fMapPi;
00514     else if(plottype=="prot")th2map= &fMapPr;
00515     else if(plottype=="nc")  th2map=&fMapNC;
00516     else if(plottype=="cc")  th2map=&fMapCC;
00517     else th2map =&fMapTH2;
00518   }
00519   else th2map=&fMapTH2;
00520 
00521   hit= th2map->find(mapkey);
00522   if(hit == th2map->end())
00523   {                                                                              
00524     TH2* h = HistMan::Instance().CreateTH2(name,plottype);
00525     if(h)   Anp::SetDir(h, dir, name);    
00526 
00527     else   ++fNMiss; 
00528 
00529     hit = (th2map->insert(MapTH2::value_type(mapkey, h))).first;
00530   }
00531 
00532   return hit -> second;
00533 }

bool Anp::PlotShort::Init (  const Header &  header  )  [virtual]
 

Reimplemented from Anp::AlgSnarl.

Definition at line 182 of file PlotShort.cxx.

References fIsData, and Anp::Header::IsData().

00182                                            {
00183   fIsData = header.IsData();
00184   return true;
00185 }

Particle::Particle_t Anp::PlotShort::Particle (  const Track &  track,
const Record &  record,
int  booltype
)  const [private]
 

Definition at line 539 of file PlotShort.cxx.

References Anp::Record::FindStdHep(), Anp::Record::GetHeader(), Anp::Header::IsData(), Anp::Record::StdHepBeg(), Anp::Record::StdHepEnd(), and Anp::StdHepIter.

00540 {
00541   if(record.GetHeader().IsData())
00542     {
00543       return Particle::kUnknown;
00544     }
00545 
00546   if(record.StdHepBeg() != record.StdHepEnd())
00547     {
00548       const StdHepIter stdhep = record.FindStdHep(track);
00549       if(stdhep == record.StdHepEnd())
00550         {
00551           return Particle::kUnknown;
00552         }
00553 
00554       if(stdhep -> Particle() == Particle::kMuonPos ||
00555          stdhep -> Particle() == Particle::kMuonNeg)
00556         {
00557           return Particle::kMuon;
00558         }
00559       else if(booltype==1)
00560         {
00561           // Made to give particle type
00562           return stdhep->Particle();
00563         }
00564       else
00565         {
00566           // If you want a boolean-like answer is/not muon (default)
00567           return Particle::kNotMuon;
00568         }
00569     }
00570 
00571   return Particle::kUnknown;
00572 }

void Anp::PlotShort::Plot (  const Track &  track,
const Record &  record
)  [private]
 

Definition at line 407 of file PlotShort.cxx.

References Anp::Track::DataBeg(), Anp::Track::DataEnd(), Anp::DataIter, fDir, Anp::Record::FindTruth(), fKeys, Anp::GetDir(), ptype, Anp::Record::TruthEnd(), and Anp::TruthIter.

00408  {
00409 
00410   const Particle::Particle_t particle = PlotShort::Particle(track, record,1);
00411    if(!ptype || !fPlot)
00412    {
00413       return;
00414    }
00415    
00416    float it1data=-.0099;
00417    float it2data=-.0099;
00418 
00419    for(DataIter it1 = track.DataBeg(); it1 != track.DataEnd(); ++it1)
00420    {     
00421      if(it1->Key()==fKeys[0])    it1data=it1->Data();          
00422      if(it1->Key()==fKeys[1])    it2data=it1->Data();
00423    }
00424 
00425    //
00426    // Fill the 2D histograms
00427    //
00428 
00429    stringstream name;
00430    string plotname;
00431    if(!fIsData){
00432      const TruthIter truth = record.FindTruth(track);
00433      if(truth!=record.TruthEnd()){
00434      if(truth->IsCC())
00435        {
00436          plotname = "twod_cc";  
00437        }
00438      else
00439        {
00440          plotname = "twod_nc";
00441        }
00442      }
00443      else plotname="twod";
00444    }
00445    else{
00446      plotname = "twod";
00447    }
00448      name << fKeys[0] << "_" << fKeys[1];
00449 
00450    TDirectory* dir = Anp::GetDir(plotname,fDir);
00451    TH2 *h = PlotShort::GetTH2(name.str(), dir,"knnvars");
00452    
00453    if(h&&(it1data>0||it2data>0))
00454      {
00455        h -> Fill(it1data, it2data);
00456      }
00457    
00458    stringstream name2;
00459    string  plotname2;
00460 
00461    if(particle == Particle::kMuon){
00462      plotname2 = "twod_mu";
00463    }
00464    else if(particle== Particle::kProton){    
00465      plotname2 = "twod_prot";
00466    }
00467    else if(particle== Particle::kPiPlus || particle== Particle::kPiMinus){
00468      plotname2 = "twod_pion";
00469    }
00470    else{
00471      plotname2 = "twod_oth";
00472    }
00473      name2<<fKeys[0]<<"_"<<fKeys[1];
00474    
00475    
00476    TDirectory* dir2 = Anp::GetDir(plotname2,fDir);
00477    TH2 *h2 = PlotShort::GetTH2(name2.str(), dir2,"knnvars");
00478 
00479    if(h2&&(it1data>0||it2data>0)){
00480      h2 -> Fill(it1data, it2data);
00481    }
00482  
00483 }

void Anp::PlotShort::Plot (  const Event &  event,
const Record &  record
)  [private]
 

Definition at line 192 of file PlotShort.cxx.

00193 {
00194 
00195 }

void Anp::PlotShort::PlotStandard (  const Track &  track,
const Record &  record
)  [private]
 

Definition at line 223 of file PlotShort.cxx.

References fDir, and Anp::GetDir().

00224 { 
00225 
00226   if(!fDir || !fPlot)
00227     {
00228       return;
00229     }
00230   string dirname;
00231   if(!fIsData){
00232     const Particle::Particle_t particle = PlotShort::Particle(track, record,1);
00233     if(particle == Particle::kMuon)           dirname="mu";    
00234     else if(particle== Particle::kProton)     dirname="prot";
00235     else if(particle== Particle::kPiPlus 
00236             || particle== Particle::kPiMinus) dirname="pion";    
00237     else                                      dirname="other";    
00238   }
00239   else dirname="track_data";
00240 
00241   TDirectory* dir = Anp::GetDir(dirname,fDir);
00242   PlotShort::FillStandard(dir, track, record);
00243   
00244 }

void Anp::PlotShort::PlotStandard (  const Event &  event,
const Record &  record
)  [private]
 

Definition at line 197 of file PlotShort.cxx.

References fDir, Anp::Record::FindTruth(), Anp::GetDir(), Anp::Event::GetNTracks(), Anp::LongestTrack(), ntrack, Anp::Record::TrackEnd(), Anp::TrackIter, Anp::Record::TruthEnd(), and Anp::TruthIter.

00198 {
00199   
00200   if(!fDir || !fPlot)  return;    
00201   string dirname;
00202   if(!fIsData){
00203     const TruthIter truth = record.FindTruth(event); 
00204     if(truth != record.TruthEnd() &&truth->IsCC())        dirname="cc";
00205     else if(truth != record.TruthEnd() && !(truth->IsCC()))  dirname="nc";
00206     else                                                  dirname="unknown";
00207   }
00208   else dirname="event_data";
00209 
00210   TDirectory* dir = Anp::GetDir(dirname,fDir);
00211 
00212   int ntrack = event.GetNTracks();
00213   PlotShort::FillStandard( dir, event,record);
00214   if(!(ntrack==0)  ){
00215     TrackIter itrack = Anp::LongestTrack(event,record);
00216     if(itrack == record.TrackEnd()) return ;
00217     PlotShort::FillStandard( dir, (*itrack), record);    
00218   }
00219 
00220   
00221 }

bool Anp::PlotShort::Run (  Record &  record  )  [virtual]
 

Implements Anp::AlgSnarl.

Definition at line 54 of file PlotShort.cxx.

References count, Anp::Record::DataAt(), Anp::Record::EventBegIterator(), Anp::Record::EventEnd(), Anp::EventIter, Anp::Record::FindEvent(), Anp::Record::FindTruth(), fPassKey, Anp::Record::TrackBegIterator(), Anp::Record::TrackEndIterator(), Anp::TrackIterator, and Anp::TruthIter.

00055 { 
00056 
00057   if(fTrack){      
00058     for(TrackIterator itrack = record.TrackBegIterator(); itrack != record.TrackEndIterator(); ++itrack)
00059       { 
00060         if(fPassKey>0)
00061           {
00062             if(!itrack->KeyExists(fPassKey)) continue;
00063             else if(itrack->DataAt(fPassKey)<=0) continue;
00064           }
00065           if(fPlotStandard) PlotShort::PlotStandard(*itrack,record);
00066           if(fSVMPrint) 
00067             {
00068               int count =0;
00069               string cc=0;                                                             
00070               const Particle::Particle_t particle = PlotShort::Particle(*itrack, record,0);      
00071               EventIter ievent = record.FindEvent(*itrack);
00072               TruthIter truth = record.FindTruth(*ievent);
00073               if(particle==Particle::kMuon && truth->IsCC()) cc="+1";                                 
00074               else cc="-1";
00075               cout<<cc;
00076               if(itrack->KeyExists(13001)) cout<<" 1:"<< itrack->DataAt(13001);
00077               else cout<<" 1:0";
00078               if(itrack->KeyExists(14410)) cout<<" 2:"<< itrack->DataAt(14410);
00079               else cout<<" 2:0";
00080               if(itrack->KeyExists(14420)) cout<<" 3:"<< itrack->DataAt(14420);
00081               else cout<<" 3:0";
00082               if(itrack->KeyExists(14430)) cout<<" 4:"<< itrack->DataAt(14430);
00083               else cout<<" 4:0";
00084               if(itrack->KeyExists(14440)) cout<<" 5:"<< itrack->DataAt(14440);
00085               else cout<<" 5:0";
00086               if(itrack->KeyExists(14905)) cout<<" 6:"<< itrack->DataAt(14905);
00087               else cout<<" 6:0";
00088               if(itrack->KeyExists(14923)) cout<<" 7:"<< itrack->DataAt(14923);
00089               else cout<<" 7:0";
00090               if(itrack->KeyExists(14924)) cout<<" 8:"<< itrack->DataAt(14924);
00091               else cout<<" 8:0";
00092               if(itrack->KeyExists(14933)) cout<<" 9:"<< itrack->DataAt(14933);
00093               else cout<<" 9:0";
00094               if(itrack->KeyExists(14934)) cout<<" 10:"<< itrack->DataAt(14934);
00095               else cout<<" 10:0";
00096               if(itrack->KeyExists(4280))  cout<<" #"<< itrack->DataAt(4280);
00097               else cout<<" #0.0";
00098               if(itrack->KeyExists(2480))  cout<<" #"<< itrack->DataAt(2480);
00099               else cout<<" #0.0";
00100               cout<<" # SVM"<<endl;
00101             }
00102           //    else PlotShort::Plot(*itrack, record);
00103 
00104       }
00105   }
00106   if(fEvent){
00107     for(EventIter ievent = record.EventBegIterator(); ievent != record.EventEnd(); ++ievent)
00108       {
00109         if(fPassKey>0 )
00110           {
00111             if(!ievent->KeyExists(fPassKey)) continue;
00112             else if(ievent->DataAt(fPassKey)<=0) continue;
00113           }
00114           
00115         if(fPlotStandard) PlotShort::PlotStandard(*ievent,record);
00116         //        else PlotShort::Plot(*ievent, record);      
00117       }
00118   }
00119 
00120    return true;
00121 }

void Anp::PlotShort::Set (  TDirectory *  dir  )  [virtual]
 

Reimplemented from Anp::AlgSnarl.

Definition at line 124 of file PlotShort.cxx.

References fDir, fDirName, fNMiss, fPlot, and Anp::GetDir().

00125 {
00126    if(!dir)
00127    {
00128       return;
00129    }
00130    
00131    fDir = Anp::GetDir(fDirName, dir);
00132 
00133    if(fNMiss == 0)
00134    {
00135       fPlot = true;
00136    }
00137    else
00138    {
00139       cerr << "PlotShort::Set - missed " << fNMiss << " histograms" << endl;
00140    }
00141 }

Member Data Documentation

float Anp::PlotShort::comp [private]
 

Definition at line 100 of file PlotShort.h.

float Anp::PlotShort::energy [private]
 

Definition at line 95 of file PlotShort.h.

Referenced by FillStandard().

bool Anp::PlotShort::fDebug [private]
 

Definition at line 76 of file PlotShort.h.

Referenced by Config().

TDirectory* Anp::PlotShort::fDir [private]
 

Definition at line 64 of file PlotShort.h.

Referenced by GetTH1(), Plot(), PlotStandard(), and Set().

std::string Anp::PlotShort::fDirName [private]
 

Definition at line 62 of file PlotShort.h.

Referenced by Config(), and Set().

std::vector<int> Anp::PlotShort::fEKeys [private]
 

Definition at line 82 of file PlotShort.h.

Referenced by Config().

bool Anp::PlotShort::fEvent [private]
 

Definition at line 69 of file PlotShort.h.

TFile* Anp::PlotShort::fFile [private]
 

Definition at line 65 of file PlotShort.h.

bool Anp::PlotShort::fIsData [private]
 

Definition at line 79 of file PlotShort.h.

Referenced by Init().

std::vector<int> Anp::PlotShort::fKeys [private]
 

Definition at line 81 of file PlotShort.h.

Referenced by Config(), and Plot().

std::vector<int> Anp::PlotShort::fKNNKeys [private]
 

Definition at line 84 of file PlotShort.h.

Referenced by Config(), and FillStandard().

std::map<std::string, TH2 *> Anp::PlotShort::fMapCC [private]
 

Definition at line 91 of file PlotShort.h.

std::map<std::string, TH2 *> Anp::PlotShort::fMapMu [private]
 

Definition at line 87 of file PlotShort.h.

std::map<std::string, TH2 *> Anp::PlotShort::fMapNC [private]
 

Definition at line 90 of file PlotShort.h.

std::map<std::string, TH2 *> Anp::PlotShort::fMapPi [private]
 

Definition at line 88 of file PlotShort.h.

std::map<std::string, TH2 *> Anp::PlotShort::fMapPr [private]
 

Definition at line 89 of file PlotShort.h.

std::map<std::string, TH1 *> Anp::PlotShort::fMapTH1 [private]
 

Definition at line 85 of file PlotShort.h.

std::map<std::string, TH2 *> Anp::PlotShort::fMapTH2 [private]
 

Definition at line 86 of file PlotShort.h.

std::map<std::string, TTree *> Anp::PlotShort::fMapTree [private]
 

Definition at line 93 of file PlotShort.h.

int Anp::PlotShort::fNMiss [private]
 

Definition at line 77 of file PlotShort.h.

Referenced by Set().

int Anp::PlotShort::fPassKey [private]
 

Definition at line 75 of file PlotShort.h.

Referenced by Config(), and Run().

bool Anp::PlotShort::fPlot [private]
 

Definition at line 72 of file PlotShort.h.

Referenced by Set().

std::string Anp::PlotShort::fPlotKey [private]
 

Definition at line 67 of file PlotShort.h.

Referenced by Config().

bool Anp::PlotShort::fPlotStandard [private]
 

Definition at line 73 of file PlotShort.h.

Referenced by Config().

bool Anp::PlotShort::fPlotTree [private]
 

Definition at line 70 of file PlotShort.h.

Referenced by Config().

bool Anp::PlotShort::fSVMPrint [private]
 

Definition at line 74 of file PlotShort.h.

Referenced by Config().

std::vector<int> Anp::PlotShort::fTKeys [private]
 

Definition at line 83 of file PlotShort.h.

Referenced by Config().

bool Anp::PlotShort::fTrack [private]
 

Definition at line 68 of file PlotShort.h.

Referenced by Config().

TTree* Anp::PlotShort::fTree [private]
 

Definition at line 66 of file PlotShort.h.

std::string Anp::PlotShort::fTreeName [private]
 

Definition at line 71 of file PlotShort.h.

Referenced by Config().

std::map<int, float> Anp::PlotShort::fTreeValues [private]
 

Definition at line 94 of file PlotShort.h.

int Anp::PlotShort::nevent [private]
 

Definition at line 106 of file PlotShort.h.

int Anp::PlotShort::nplane [private]
 

Definition at line 98 of file PlotShort.h.

int Anp::PlotShort::ntrack [private]
 

Definition at line 97 of file PlotShort.h.

Referenced by PlotStandard().

int Anp::PlotShort::ntrackplane [private]
 

Definition at line 99 of file PlotShort.h.

int Anp::PlotShort::ptype [private]
 

Definition at line 96 of file PlotShort.h.

Referenced by Plot().

float Anp::PlotShort::pur [private]
 

Definition at line 101 of file PlotShort.h.

float Anp::PlotShort::tcompall [private]
 

Definition at line 103 of file PlotShort.h.

float Anp::PlotShort::tcompslc [private]
 

Definition at line 104 of file PlotShort.h.

float Anp::PlotShort::tpur [private]
 

Definition at line 102 of file PlotShort.h.

int Anp::PlotShort::ttype [private]
 

Definition at line 105 of file PlotShort.h.

The documentation for this class was generated from the following files:
Generated on Mon Feb 15 11:10:33 2010 for loon by  doxygen 1.3.9.1